FTFD: an informatics pipeline supporting phylogenomic analysis of fungal transcription factors

Molecular Analysis of UV-C Induced Resveratrol Accumulation in Polygonum cuspidatum Leaves

Resveratrol is one of the most studied plant secondary metabolites owing to its numerous health benefits. It is accumulated in some plants following biotic and abiotic stress pressures, including UV-C irradiation. Polygonum cuspidatum represents the major natural source of concentrated resveratrol but the underlying mechanisms as well as the effects of UV-C irradiation on resveratrol content have not yet been documented. Herein, we found that UV-C irradiation significantly increased by 2.6-fold and 1.6-fold the resveratrol content in irradiated leaf samples followed by a dark incubation for 6 h and 12 h, respectively, compared to the untreated samples. De novo transcriptome sequencing and assembly resulted into 165,013 unigenes with 98 unigenes mapped to the resveratrol biosynthetic pathway. Differential expression analysis showed that P.cuspidatum strongly induced the genes directly involved in the resveratrol synthesis, including phenylalanine ammonia-lyase, cinnamic acid 4-hydroxylase, 4-coumarate-CoA ligase and stilbene synthase (STS) genes, while strongly decreased the chalcone synthase (CHS) genes after exposure to UV-C. Since CHS and STS share the same substrate, P. cuspidatum tends to preferentially divert the substrate to the resveratrol synthesis pathway under UV-C treatment. We identified several members of the MYB, bHLH and ERF families as potential regulators of the resveratrol biosynthesis genes.

rmtA-Dependent Transcriptome and Its Role in Secondary Metabolism, Environmental Stress, and Virulence in Aspergillus flavus

Aspergillus flavus colonizes numerous oil seed crops such as maize, peanuts, treenuts and cottonseed worldwide, contaminating them with aflatoxins and other harmful toxins. Previously our lab characterized the gene rmtA, which encodes an arginine methyltransferase in A. flavus, and demonstrated its role governing the expression of regulators in the aflatoxin gene cluster and subsequent synthesis of toxin. Furthermore, our studies revealed that rmtA also controls conidial and sclerotial development implicating it as an epigenetic regulator in A. flavus. To confirm this, we performed a RNA sequencing analysis to ascertain the extent of rmtA’s influence on the transcriptome of A. flavus. In this analysis we identified over 2000 genes that were rmtA-dependent, including over 200 transcription factor genes, as well as an uncharacterized secondary metabolite gene cluster possibly responsible for the synthesis of an epidithiodiketopiperazine-like compound. Our results also revealed rmtA-dependent genes involved in multiple types of abiotic stress response in A. flavus. Importantly, hundreds of genes active during maize infection were also regulated by rmtA. In addition, in the animal infection model, rmtA was dispensable for virulence, however forced overexpression of rmtA increased mortality with respect to the wild type.

Identification and characterization of circular RNAs in Ganoderma lucidum

Circular RNAs (circRNAs) play important roles in animals, plants, and fungi. However, no circRNAs have been reported in Ganoderma lucidum. Here, we carried out a genome-wide identification of the circRNAs in G.lucidum using RNA-Seq data, and analyzed their features. In total, 250 and 2193 circRNAs were identified from strand-specific RNA-seq data generated from the polyA(−) and polyA(−)/RNase R-treated libraries, respectively. Six of 131 (4.58%) predicted circRNAs were experimentally confirmed. Across three developmental stages, 731 exonic circRNAs (back spliced read counts ≥ 5) and their parent genes were further analyzed. CircRNAs were preferred originating from exons with flanking introns, and the lengths of the flanking intron were longer than those of the control introns. A total of 200 circRNAs were differentially expressed across the three developmental stages of G. lucidum. The expression profiles of 119 (16.3%) exonic circRNAs and their parent genes showed significant positive correlations (r ≥ 0.9, q < 0.01), whereas 226 (30.9%) exonic circRNAs and their parent genes exhibited significant negative correlations (r ≤ −0.9, q < 0.01), in which 53 parent genes are potentially involved in the transcriptional regulation, polysaccharide biosynthesis etc. Our results indicated that circRNAs are present in G. lucidum, with potentially important regulatory roles.

A MYST family histone acetyltransferase, MoSAS3, is required for development and pathogenicity in the rice blast fungus

Histone acetylation has been established as a principal epigenetic regulatory mechanism in eukaryotes. Sas3, a histone acetyltransferase belonging to the largest family of acetyltransferase, MYST, is the catalytic subunit of a conserved histone acetyltransferase complex. To date, the functions of Sas3 and its orthologues have been extensively studied in yeast, humans and flies in relation to global acetylation and transcriptional regulation. However, its precise impact on development and pathogenicity in fungal plant pathogens has yet to be elucidated. Considering the importance of Sas3 in H3K14 acetylation, here we investigate the roles of its orthologue in the rice blast fungus, Magnaporthe oryzae (Pyricularia oryzae). Unlike a previously reported Sas3 deletion in yeast, which led to no remarkable phenotypic changes, we found that MoSAS3 deletion alone had a profound effect on fungal growth and development, including asexual reproduction, germination and appressorium formation in M. oryzae. Such defects in pre-penetration development resulted in complete loss of pathogenicity in the deletion mutant. Furthermore, genetic analysis of MoSAS3 and MoGCN5 encoding a Gcn5-related N-acetyltransferase family histone acetyltransferase suggested that two conserved components of histone acetylation are integrated differently into epigenetic regulatory mechanisms in the yeast and a filamentous fungus. RNA-seq analysis of ΔMosas3 showed two general trends: many DNA repair and DNA damage response genes are up-regulated, while carbon and nitrogen metabolism genes are down-regulated in ΔMosas3. Our work demonstrates the importance of MYST family histone acetyltransferase as a developmental regulator and illuminates a degree of functional variation in conserved catalytic subunits among different fungal species.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata

BACKGROUND

Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution.

RESULTS

A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting.

CONCLUSIONS

The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

Draft Genome Sequence of Amphirosellinia nigrospora JS-1675, an Endophytic Fungus from Pteris cretica

The fungus Amphirosellinia nigrospora strain JS-1675 has been reported to exert antimicrobial effects against various plant-pathogenic bacteria and fungi. Here, we report the draft genome sequence of A. nigrospora for the first time.

The fungus Amphirosellinia nigrospora strain JS-1675 has been reported to exert antimicrobial effects against various plant-pathogenic bacteria and fungi. Here, we report the draft genome sequence of A. nigrospora for the first time. The assembly comprises 48,177,783 bp with 18 scaffolds.

Regulation of conidiation, dimorphic transition, and microsclerotia formation by MrSwi6 transcription factor in dimorphic fungus Metarhizium rileyi

Microsclerotia (MS) produced in the liquid culture of the dimorphic insect pathogen Metarhizium rileyi can be used as a mycoinsecticide. Bioinformatics analysis demonstrated that the cell cycle signaling pathway was involved in regulating MS formation. To investigate the mechanisms by which the signaling pathway is regulated, a cell cycle box binding transcription factor MrSwi6 of M. rileyi was characterized. MrSwi6 was highly expressed during periods of yeast-hypha transition and conidia and MS formation. When compared with wild-type and complemented strains, disruption of MrSwi6 significantly reduced conidia (15-36%) and MS formation (96.2%), and exhibited decreased virulence levels. Digital expression profiling revealed that genes involved in antioxidation, pigment biosynthesis, and ion transport and storage were regulated by MrSwi6 during conidia and MS development. These results confirmed the significance of MrSwi6 in dimorphic transition, conidia and MS formation, and virulence in M. rileyi.

A High-Quality Draft Genome Sequence of Colletotrichum gloeosporioides sensu stricto SMCG1#C, a Causal Agent of Anthracnose on Cunninghamia lanceolata in China

Colletotrichum has a broad host range and causes major yield losses of crops. The fungus Colletotrichum gloeosporioides is associated with anthracnose on Chinese fir. In this study, we present a high-quality draft genome sequence of C. gloeosporioides sensu stricto SMCG1#C, providing a reference genomic data for further research on anthracnose of Chinese fir and other hosts.

The Transcriptional Regulator Hbx1 Affects the Expression of Thousands of Genes in the Aflatoxin-Producing Fungus Aspergillus flavus

In filamentous fungi, homeobox proteins are conserved transcriptional regulators described to control conidiogenesis and fruiting body formation. Eight homeobox (hbx) genes are found in the genome of the aflatoxin-producing ascomycete, Aspergillus flavus While loss-of-function of seven of the eight genes had little to no effect on fungal growth and development, disruption of hbx1, resulted in aconidial colonies and lack of sclerotial production. Furthermore, the hbx1 mutant was unable to produce aflatoxins B1 and B2, cyclopiazonic acid and aflatrem. In the present study, hbx1 transcriptome analysis revealed that hbx1 has a broad effect on A. flavus gene expression, and the effect of hbx1 increases overtime, impacting more than five thousand protein-coding genes. Among the affected genes, those in the category of secondary metabolism (SM), followed by that of cellular transport, were the most affected. Specifically, regarding the effect of hbx1 on SM, we found that genes in 44 SM gene clusters where upregulated while 49 were downregulated in the absence of hbx1, including genes in the SM clusters responsible for the synthesis of asparasone, piperazine and aflavarin, all known to be associated with sclerotia. In addition, our study revealed that hbx1 affects the expression of other transcription factor genes involved in development, including the conidiation central regulatory pathway and flb genes.

The methionine biosynthesis regulator AaMetR contributes to oxidative stress tolerance and virulence in Alternaria alternata

The tangerine pathotype of A. alternata, which produces a unique host-selective ACT toxin causes brown spots on citrus leaves and fruits. In this study, we report a methionine biosynthesis regulator (MetR), which belong to bZIP transcription factor, is required for methionine metabolism, oxidative stress tolerance and pathogenicity. We generated two ΔAaMetR mutants in the tangerine pathotype of Alternaria alternata and investigated the resulting mutant phenotypes. The ΔAaMetR disruption mutant grew poorly in the absence of methionine and unable to produce conidia. Furthermore, pathogenicity tests have shown that ΔAaMetR mutant on their tangerine host can neither penetrate nor cause disease. These ΔAaMetR mutants exhibit an increased sensitivity to exogenous H₂O₂ and many ROS generating oxidants. To elucidate the transcription network of AaMetR, we performed RNA-Seq experiments on wild-type and ΔAaMetR mutant and identified genes that were differentially expressed between the two genotypes. Transcriptome data demonstrated that AaMetR contributes in many other biological processes including ROS detoxification, sulfur transfer, and amino acid metabolism. Comparative transcriptome analysis indicated that the ΔAaMetR mutant up-regulated several genes involved in cysteine and methionine metabolism. In conclusion, our results highlight the global regulatory role of AaMetR in cysteine and methionine metabolism and provide new insights into the crucial role of ROS detoxification, sporulation and pathogenicity in the tangerine pathotype of A. alternata.

A transcription factor, MrMsn2, in the dimorphic fungus Metarhizium rileyi is essential for dimorphism transition, aggravated pigmentation, conidiation and microsclerotia formation

Microsclerotia (MS) are pseudoparenchymatous aggregations of hyphae of fungi that can be induced in liquid culture for biocontrol applications. Previously, we determined that the high-osmolarity glycerol (HOG) signalling pathway was involved in regulating MS development in the dimorphic insect pathogen Metarhizium rileyi. To further investigate the mechanisms by which the signalling pathway is regulated, we characterized the transcriptional factor MrMsn2, a homologue of the yeast C2 H2 transcriptional factor Msn2, which is predicted to function downstream of the HOG pathway in M. rileyi. Compared with wild-type and complemented strains, disruption of MrMsn2 increased the yeast-to-hypha transition rate, enhanced conidiation capacity and aggravated pigmentation in M. rileyi. The ▵MrMsn2 mutants were sensitive to stress, produced morphologically abnormal clones and had significantly reduced MS formation and decreased virulence levels. Digital expression profiling revealed that genes involved in antioxidation, pigment biosynthesis and ion transport and storage were regulated by MrMsn2 during conidia and MS development. Taken together, our findings confirm that MrMsn2 controlled the yeast-to-hypha transition, conidia and MS formation, and virulence.

Comparative Transcriptome Analysis Between a Spontaneous Albino Mutant and Its Sibling Strain of Cordyceps militaris in Response to Light Stress

Albinism has been used for new variety screening in some edible mushrooms and the underlying mechanisms are fascinating. Albino fruiting body of Cordyceps militaris, a well-known edible fungus and model organism for Cordyceps, has the potential to be a nutraceutical or functional food due to its high content of metabolites and antioxidant activities. In this study, a spontaneous albino mutant strain (505) of C. militaris was obtained. In comparison to its normal sibling strain (498), the albino strain stably remained white in response to light and had significantly decreased conidia and carotenoid production but accumulated more cordycepin. Transcriptome analysis of both strains revealed that all the seven photoreceptors were expressed similarly in response to light. However, many more genes in the albino strain were differentially expressed in response to light than its sibling strain. The significantly enriched pathways in 498L vs. 505L were mainly associated with replication and repair. Some secondary metabolite backbone genes including those encoding DMAT, two NRPS-like proteins, three NPRS, and lanosterol synthase were differentially expressed in the albino when compared with that of the normal strains. Transcriptome and real-time quantitative PCR analyses indicated that some cytochrome P450s and methyltransferases might be related to the phenotypic differences observed between the two strains. This study compared the genome-wide transcriptional responses to light irradiation in a spontaneous albino mutant and its normal sibling strain of an edible fungus, and these findings potentially pave the way for further investigation of the pigment biosynthetic pathway.

Enhanced hypocrellin production via coexpression of alpha-amylase and hemoglobin genes in Shiraia bambusicola

Shiraia bambusicola is an important and valuable macrofungus and hypocrellins are its main secondary metabolites which have been widely applied in many medical fields. However, during SSF process of this filamentous fungus, use ratio of corn substrate and dissolved oxygen supply are two main limiting factors, which influence production cost, yield and product quality. To solve these problems, overexpressions of amy365-1 and vgb in S. bambusicola were investigated and three overexpression transformants were constructed. Results demonstrated that expressions and coexpression of AMY365-1 and VHb not only increased the productions of biomass, amylase, hypocrellin, but also up-regulated relative expression levels of four central carbon metabolism genes (pdc, ald, acs, acc) and seven hypocrellin biosynthesis genes (fad, mono, zftf, omef, msf, pks, mco). Furthermore, expression of VHb decreased SSF period. When amy365-1 and vgb were coexpressed, relative expression levels of zftf and pks reached their highest levels at 72 h under liquid fermentation, hypocrellin production reached the highest level 75.85 mg/gds which was 2.99-fold compared with wild type strain within 11 days under SSF, and residual starch of solid substrates was decreased from 35.47 to 14.57%.

Verticillium dahliae transcription factor VdFTF1 regulates the expression of multiple secreted virulence factors and is required for full virulence in cotton

Fungal transcription factors (TFs) implicated in the regulation of virulence gene expression have been identified in a number of plant pathogens. In Verticillium dahliae, despite its agricultural importance, few regulators of transcription have been characterized. In this study, a T-DNA insertion mutant with significantly reduced virulence towards cotton was identified. The T-DNA was traced to VdFTF1, a gene encoding a TF containing a Fungal_trans domain. Transient expression in onion epidermal cells indicated that VdFTF1 is localized to the nucleus. The VdFTF1-deletion strains displayed normal vegetative growth, mycelial pigmentation and conidial morphology, but exhibited significantly reduced virulence on cotton, suggesting that VdFTF1 is required exclusively for pathogenesis. Comparisons of global transcription patterns of wild-type and VdFTF1-deletion strains indicated that VdFTF1 affected the expression of 802 genes, 233 of which were associated with catalytic processes. These genes encoded 69 potentially secreted proteins, 43 of which contained a carbohydrate enzyme domain known to participate in pathogenesis during infection of cotton. Targeted gene deletion of one VdFTF1-regulated gene resulted in significantly impaired vascular colonization, as measured by quantitative polymerase chain reaction, as well as aggressiveness and symptom severity in cotton. In conclusion, VdFTF1, which encodes a TF containing a Fungal_trans domain, regulates the gene expression of plant cell wall degradation enzymes in V. dahliae, which are required for full virulence on cotton.

Collection and Curation of Transcriptional Regulatory Interactions in Aspergillus nidulans and Neurospora crassa Reveal Structural and Evolutionary Features of the Regulatory Networks

Transcriptional regulation has important roles in various biological processes (e.g., development and metabolism) in filamentous fungi. However, regulatory interactions between transcription factors (TFs) and their target genes in these species have only been described in different forms by primary scientific literature, which limits the integrated analysis of these data. Here, we extensively curated the reported transcriptional regulatory interactions in Aspergillus nidulans and Neurospora crassa. For each interaction, the identifiers of involved proteins or genes were unified, and the types of supporting experiments were recorded. Then, transcriptional regulatory networks were reconstructed from the interactions supported by classical low-throughput experiments. Analysis of the networks revealed the presence of hub targets regulated by multiple TFs and network motifs of other structures (e.g., regulatory loops). Comparison of the regulatory interactions between the two species identified 33 conserved interactions supported by classical experiments in both species, most of which are involved in the regulation of metabolic genes. We anticipate the curated data would serve as a catalog for the studies of transcriptional regulation in filamentous fungi.

Genome-wide identification and analysis of the basic leucine zipper (bZIP) transcription factor gene family in Ustilaginoidea virens

The basic leucine zipper (bZIP) transcription factor (TF) family is one of the largest and most diverse TF families widely distributed across the eukaryotes. The bZIP TF family plays an important role in growth, development, and response to abiotic or biotic stresses, which have been well characterized in plants, but not in plant pathogenic fungi. In this study, we performed genome-wide and systematic bioinformatics analysis of bZIP genes in the fungus Ustilaginoidea virens, the causal agent of rice false smut disease. We identified 28 bZIP family members in the U. virens genome by searching for the bZIP domain in predicted genes. The gene structures, motifs, and phylogenetic relationships were analyzed for bZIP genes in U. virens (UvbZIP). Together with bZIP proteins from two other fungi, the bZIP genes can be divided into eight groups according to their phylogenetic relationships. Based on RNA-Seq data, the expression profiles of UvbZIP genes at different infection stages were evaluated. Results showed that 17 UvbZIP genes were up-regulated during the infection period. Furthermore, 11 infection-related UvbZIP genes were investigated under H2O2 stress and the expression level of eight genes were changed, which confirmed their role in stress tolerance and pathogenicity. In summary, our genome-wide systematic characterization and expression analysis of UvbZIP genes provided insight into the molecular function of these genes in U. virens and provides a reference for other pathogens.

Genome Sequence of an Endophytic Fungus, Fusarium solani JS-169, Which Has Antifungal Activity

An endophytic fungus, Fusarium solani strain JS-169, isolated from a mulberry twig, showed considerable antifungal activity. Here, we report the draft genome sequence of this strain. The assembly comprises 17 scaffolds, with an N₅₀ value of 4.93 Mb. The assembled genome was 45,813,297 bp in length, with a G+C content of 49.91%.

Profiling of secondary metabolite gene clusters regulated by LaeA in Aspergillus niger FGSC A1279 based on genome sequencing and transcriptome analysis

The global regulator LaeA controls the production of many fungal secondary metabolites, possibly via chromatin remodeling. Here we aimed to survey the secondary metabolite profile regulated by LaeA in Aspergillus niger FGSC A1279 by genome sequencing and comparative transcriptomics between the laeA deletion (ΔlaeA) and overexpressing (OE-laeA) mutants. Genome sequencing revealed four putative polyketide synthase genes specific to FGSC A1279, suggesting that the corresponding polyketide compounds might be unique to FGSC A1279. RNA-seq data revealed 281 putative secondary metabolite genes upregulated in the OE-laeA mutants, including 22 secondary metabolite backbone genes. LC-MS chemical profiling illustrated that many secondary metabolites were produced in OE-laeA mutants compared to wild type and ΔlaeA mutants, providing potential resources for drug discovery. KEGG analysis annotated 16 secondary metabolite clusters putatively linked to metabolic pathways. Furthermore, 34 of 61 Zn₂Cys₆ transcription factors located in secondary metabolite clusters were differentially expressed between ΔlaeA and OE-laeA mutants. Three secondary metabolite clusters (cluster 18, 30 and 33) containing Zn₂Cys₆ transcription factors that were upregulated in OE-laeA mutants were putatively linked to KEGG pathways, suggesting that Zn₂Cys₆ transcription factors might play an important role in synthesizing secondary metabolites regulated by LaeA. Taken together, LaeA dramatically influences the secondary metabolite profile in FGSC A1279.

Draft Genome Sequence of Aspergillus persii NIBRFGC000004109, Which Has Antibacterial Activity against Plant-Pathogenic Bacteria

The fungus Aspergillus persii strain NIBRFGC000004109 is capable of producing penicillic acid and showed antibacterial activity against various plant-pathogenic bacteria, including Xanthomonas arboricola pv. pruni. Here, we report the first draft whole-genome sequence of A. persii. The assembly comprises 38,414,373 bp, with 12 scaffolds.

Gapless genome assembly of Colletotrichum higginsianum reveals chromosome structure and association of transposable elements with secondary metabolite gene clusters

BACKGROUND

The ascomycete fungus Colletotrichum higginsianum causes anthracnose disease of brassica crops and the model plant Arabidopsis thaliana. Previous versions of the genome sequence were highly fragmented, causing errors in the prediction of protein-coding genes and preventing the analysis of repetitive sequences and genome architecture.

RESULTS

Here, we re-sequenced the genome using single-molecule real-time (SMRT) sequencing technology and, in combination with optical map data, this provided a gapless assembly of all twelve chromosomes except for the ribosomal DNA repeat cluster on chromosome 7. The more accurate gene annotation made possible by this new assembly revealed a large repertoire of secondary metabolism (SM) key genes (89) and putative biosynthetic pathways (77 SM gene clusters). The two mini-chromosomes differed from the ten core chromosomes in being repeat- and AT-rich and gene-poor but were significantly enriched with genes encoding putative secreted effector proteins. Transposable elements (TEs) were found to occupy 7% of the genome by length. Certain TE families showed a statistically significant association with effector genes and SM cluster genes and were transcriptionally active at particular stages of fungal development. All 24 subtelomeres were found to contain one of three highly-conserved repeat elements which, by providing sites for homologous recombination, were probably instrumental in four segmental duplications.

CONCLUSION

The gapless genome of C. higginsianum provides access to repeat-rich regions that were previously poorly assembled, notably the mini-chromosomes and subtelomeres, and allowed prediction of the complete SM gene repertoire. It also provides insights into the potential role of TEs in gene and genome evolution and host adaptation in this asexual pathogen.

Draft Genome Sequence of the Fungus Associated with Oak Wilt Mortality in South Korea, Raffaelea quercus-mongolicae KACC44405

The fungus Raffaelea quercus-mongolicae is the causal agent of Korean oak wilt, a disease associated with mass mortality of oak trees (e.g., Quercus spp.). The fungus is vectored and dispersed by the ambrosia beetle, Platypus koryoensis Here, we present the 27.0-Mb draft genome sequence of R. quercus-mongolicae strain KACC44405.

Draft Genome Sequence of an Endophytic Fungus, Gaeumannomyces sp. Strain JS-464, Isolated from a Reed Plant, Phragmites communis

An endophytic fungus, Gaeumannomyces sp. strain JS-464, is capable of producing a number of secondary metabolites which showed significant nitric oxide reduction activity. The draft genome assembly has a size of 53,151,282 bp, with a G+C content of 53.11% consisting of 80 scaffolds with an N₅₀ of 7.46 Mbp.

Transcription Factor Repertoire of Necrotrophic Fungal Phytopathogen Ascochyta rabiei: Predominance of MYB Transcription Factors As Potential Regulators of Secretome

Transcription factors (TFs) are the key players in gene expression and their study is highly significant for shedding light on the molecular mechanisms and evolutionary history of organisms. During host-pathogen interaction, extensive reprogramming of gene expression facilitated by TFs is likely to occur in both host and pathogen. To date, the knowledge about TF repertoire in filamentous fungi is in infancy. The necrotrophic fungus Ascochyta rabiei, that causes destructive Ascochyta blight (AB) disease of chickpea (Cicer arietinum), demands more comprehensive study for better understanding of Ascochyta-legume pathosystem. In the present study, we performed the genome-wide identification and analysis of TFs in A. rabiei. Taking advantage of A. rabiei genome sequence, we used a bioinformatic approach to predict the TF repertoire of A. rabiei. For identification and classification of A. rabiei TFs, we designed a comprehensive pipeline using a combination of BLAST and InterProScan software. A total of 381 A. rabiei TFs were predicted and divided into 32 fungal specific families of TFs. The gene structure, domain organization and phylogenetic analysis of abundant families of A. rabiei TFs were also carried out. Comparative study of A. rabiei TFs with that of other necrotrophic, biotrophic, hemibiotrophic, symbiotic, and saprotrophic fungi was performed. It suggested presence of both conserved as well as unique features among them. Moreover, cis-acting elements on promoter sequences of earlier predicted A. rabiei secretome were also identified. With the help of published A. rabiei transcriptome data, the differential expression of TF and secretory protein coding genes was analyzed. Furthermore, comprehensive expression analysis of few selected A. rabiei TFs using quantitative real-time polymerase chain reaction revealed variety of expression patterns during host colonization. These genes were expressed in at least one of the time points tested post infection. Overall, this study illustrates the first genome-wide identification and analysis of TF repertoire of A. rabiei. This work would provide the basis for further studies to dissect role of TFs in the molecular mechanisms during A. rabiei-chickpea interactions.

Lentinula edodes Genome Survey and Postharvest Transcriptome Analysis

Lentinula edodes is a popular, cultivated edible and medicinal mushroom. Lentinula edodes is susceptible to postharvest problems, such as gill browning, fruiting body softening, and lentinan degradation. We constructed a de novo assembly draft genome sequence and performed gene prediction for Lentinula edodesDe novo assembly was carried out using short reads from paired-end and mate-paired libraries and by using long reads by PacBio, resulting in a contig number of 1,951 and an N₅₀ of 1 Mb. Furthermore, we predicted genes by Augustus using transcriptome sequencing (RNA-seq) data from the whole life cycle of Lentinula edodes, resulting in 12,959 predicted genes. This analysis revealed that Lentinula edodes lacks lignin peroxidase. To reveal genes involved in the loss of quality of Lentinula edodes postharvest fruiting bodies, transcriptome analysis was carried out using serial analysis of gene expression (SuperSAGE). This analysis revealed that many cell wall-related enzymes are upregulated after harvest, such as β-1,3-1,6-glucan-degrading enzymes in glycoside hydrolase (GH) families GH5, GH16, GH30, GH55, and GH128, and thaumatin-like proteins. In addition, we found that several chitin-related genes are upregulated, such as putative chitinases in GH family 18, exochitinases in GH20, and a putative chitosanase in GH family 75. The results suggest that cell wall-degrading enzymes synergistically cooperate for rapid fruiting body autolysis. Many putative transcription factor genes were upregulated postharvest, such as genes containing high-mobility-group (HMG) domains and zinc finger domains. Several cell death-related proteins were also upregulated postharvest.IMPORTANCE Our data collectively suggest that there is a rapid fruiting body autolysis system in Lentinula edodes The genes for the loss of postharvest quality newly found in this research will be targets for the future breeding of strains that keep fresh longer than present strains. De novoLentinula edodes genome assembly data will be used for the construction of a complete Lentinula edodes chromosome map for future breeding.

Transcriptome Analyses Shed New Insights into Primary Metabolism and Regulation of Blumeria graminis f. sp. tritici during Conidiation

Conidia of the obligate biotrophic fungal pathogen Blumeria graminis f. sp. tritici (Bgt) play a vital role in its survival and rapid dispersal. However, little is known about the genetic basis for its asexual reproduction. To uncover the primary metabolic and regulatory events during conidiation, we sequenced the transcriptome of Bgt epiphytic structures at 3 (vegetative hyphae growth), 4 (foot cells initiation), and 5 (conidiophore erection) days post-inoculation (dpi). RNA-seq analyses identified 556 and 404 (combined 685) differentially expressed genes (DEGs) at 4 and 5 dpi compared with their expression levels at 3 dpi, respectively. We found that several genes involved in the conversion from a variety of sugars to glucose, glycolysis, the tricarboxylic acid cycle (TAC), the electron transport chain (ETC), and unsaturated fatty acid oxidation were activated during conidiation, suggesting that more energy supply is required during this process. Moreover, we found that glucose was converted into glycogen, which was accumulated in developing conidiophores, indicating that it could be the primary energy storage molecule in Bgt conidia. Clustering for the expression profiles of 91 regulatory genes showed that calcium (Ca²⁺), H₂O₂, and phosphoinositide (PIP) signaling were involved in Bgt conidiation. Furthermore, a strong accumulation of H₂O₂ in developing conidiophores was detected. Application of EGTA, a Ca²⁺ chelator, and trifluoperazine dihydrochloride (TFP), a calmodulin (CaM) antagonist, markedly suppressed the generation of H₂O₂, affected foot cell and conidiophore development and reduced conidia production significantly. These results suggest that Ca²⁺ and H₂O₂ signaling play important roles in conidiogenesis and a crosslink between them is present. In addition to some conidiation-related orthologs known in other fungi, such as the velvet complex components, we identified several other novel B. graminis-specific genes that have not been previously found to be implicated in fungal conidiation, reflecting a unique molecular mechanism underlying asexual development of cereal powdery mildews.

CYPome of the conifer pathogen Heterobasidion irregulare: Inventory, phylogeny, and transcriptional analysis of the response to biocontrol

The molecular mechanisms underlying the interaction of the pathogen, Heterobasidion annosum s.l., the conifer tree and the biocontrol fungus, Phlebiopsis gigantea have not been fully elucidated. Members of the cytochrome P450 (CYP) protein family may contribute to the detoxification of components of chemical defence of conifer trees by H. annosum during infection. Additionally, they may also be involved in the interaction between H. annosum and P. gigantea. A genome-wide analysis of CYPs in Heterobasidion irregulare was carried out alongside gene expression studies. According to the Standardized CYP Nomenclature criteria, the H. irregulare genome has 121 CYP genes and 17 CYP pseudogenes classified into 11 clans, 35 families, and 64 subfamilies. Tandem CYP arrays originating from gene duplications and belonging to the same family and subfamily were found. Phylogenetic analysis showed that all the families of H. irregulare CYPs were monophyletic groups except for the family CYP5144. Microarray analysis revealed the transcriptional pattern for 130 transcripts of CYP-encoding genes during growth on culture filtrate produced by P. gigantea. The high level of P450 gene diversity identified in this study could result from extensive gene duplications presumably caused by the high metabolic demands of H. irregulare in its ecological niches.

Transcriptome-Based Discovery of Fusarium graminearum Stress Responses to FgHV1 Infection

Fusarium graminearum hypovirus 1 (FgHV1), which is phylogenetically related to Cryphonectria hypovirus 1 (CHV1), is a virus in the family Hypoviridae that infects the plant pathogenic fungus F. graminearum. Although hypovirus FgHV1 infection does not attenuate the virulence of the host (hypovirulence), it results in defects in mycelial growth and spore production. We now report that the vertical transmission rate of FgHV1 through asexual spores reached 100%. Using RNA deep sequencing, we performed genome-wide expression analysis to reveal phenotype-related genes with expression changes in response to FgHV1 infection. A total of 378 genes were differentially expressed, suggesting that hypovirus infection causes a significant alteration of fungal gene expression. Nearly two times as many genes were up-regulated as were down-regulated. A differentially expressed gene enrichment analysis identified a number of important pathways. Metabolic processes, the ubiquitination system, and especially cellular redox regulation were the most affected categories in F. graminearum challenged with FgHV1. The p20, encoded by FgHV1 could induce H₂O₂ accumulation and hypersensitive response in Nicotiana benthamiana leaves. Moreover, hypovirus FgHV1 may regulate transcription factors and trigger the RNA silencing pathway in F. graminearum.

Schizophyllum commune has an extensive and functional alternative splicing repertoire

Recent genome-wide studies have demonstrated that fungi possess the machinery to alternatively splice pre-mRNA. However, there has not been a systematic categorization of the functional impact of alternative splicing in a fungus. We investigate alternative splicing and its functional consequences in the model mushroom forming fungus Schizophyllum commune. Alternative splicing was demonstrated for 2,285 out of 12,988 expressed genes, resulting in 20% additional transcripts. Intron retentions were the most common alternative splicing events, accounting for 33% of all splicing events, and 43% of the events in coding regions. On the other hand, exon skipping events were rare in coding regions (1%) but enriched in UTRs where they accounted for 57% of the events. Specific functional groups, including transcription factors, contained alternatively spliced genes. Alternatively spliced transcripts were regulated differently throughout development in 19% of the 2,285 alternatively spliced genes. Notably, 69% of alternatively spliced genes have predicted alternative functionality by loss or gain of functional domains, or by acquiring alternative subcellular locations. S. commune exhibits more alternative splicing than any other studied fungus. Taken together, alternative splicing increases the complexity of the S. commune proteome considerably and provides it with a rich repertoire of alternative functionality that is exploited dynamically.

Genome Sequence of the Edible Cultivated Mushroom Lentinula edodes (Shiitake) Reveals Insights into Lignocellulose Degradation

Lentinula edodes, one of the most popular, edible mushroom species with a high content of proteins and polysaccharides as well as unique aroma, is widely cultivated in many Asian countries, especially in China, Japan and Korea. As a white rot fungus with lignocellulose degradation ability, L. edodes has the potential for application in the utilization of agriculture straw resources. Here, we report its 41.8-Mb genome, encoding 14,889 predicted genes. Through a phylogenetic analysis with model species of fungi, the evolutionary divergence time of L. edodes and Gymnopus luxurians was estimated to be 39 MYA. The carbohydrate-active enzyme genes in L. edodes were compared with those of the other 25 fungal species, and 101 lignocellulolytic enzymes were identified in L. edodes, similar to other white rot fungi. Transcriptome analysis showed that the expression of genes encoding two cellulases and 16 transcription factor was up-regulated when mycelia were cultivated for 120 minutes in cellulose medium versus glucose medium. Our results will foster a better understanding of the molecular mechanism of lignocellulose degradation and provide the basis for partial replacement of wood sawdust with agricultural wastes in L. edodes cultivation.

Characterization of 47 Cys2 -His2 zinc finger proteins required for the development and pathogenicity of the rice blast fungus Magnaporthe oryzae

The Cys2 -His2 (C2H2) zinc finger protein family is the second-largest family of transcription factors (TFs) in Magnaporthe oryzae, the causal fungus responsible for the destructive rice blast disease. However, little is known about the roles of most C2H2 TFs in the development and pathogenicity of M. oryzae. The roles of 47 C2H2 genes in development and pathogenicity were investigated by gene deletion in M. oryzae. The TF-dependent genes in mycelia or appressoria were analyzed with RNA sequencing and quantitative PCR (qPCR). Forty-four C2H2 genes are involved in growth (20 genes), conidiation (28 genes), appressorium formation (four genes) and pathogenicity (22 genes) in M. oryzae. Of these, MGG_14931, named as VRF1, is required for pathogenicity, specifically controlling appressorium maturation by affecting the expression of genes related to appressorial structure and function, including melanin biosynthesis, chitin catabolism, lipid metabolism, proteolysis, transmembrane transport, and response to oxidative stress; MGG_01776, named as VRF2, is required for plant penetration and invasive growth; conidiation-related gene CON7 is required for conidial differentiation; and MoCREA, encoding a carbon catabolite repression protein, is a novel repressor of lipid catabolism when glucose obtainable in M. oryzae. This study provides many insights into the regulation of growth, asexual development, appressorium formation, and pathogenicity by C2H2 TFs in M. oryzae.

Compartmentalized gene regulatory network of the pathogenic fungus Fusarium graminearum

Head blight caused by Fusarium graminearum threatens world-wide wheat production, resulting in both yield loss and mycotoxin contamination. We reconstructed the global F. graminearum gene regulatory network (GRN) from a large collection of transcriptomic data using Bayesian network inference, a machine-learning algorithm. This GRN reveals connectivity between key regulators and their target genes. Focusing on key regulators, this network contains eight distinct but interwoven modules. Enriched for unique functions, such as cell cycle, DNA replication, transcription, translation and stress responses, each module exhibits distinct expression profiles. Evolutionarily, the F. graminearum genome can be divided into core regions shared with closely related species and variable regions harboring genes that are unique to F. graminearum and perform species-specific functions. Interestingly, the inferred top regulators regulate genes that are significantly enriched from the same genomic regions (P < 0.05), revealing a compartmentalized network structure that may reflect network rewiring related to specific adaptation of this plant pathogen. This first-ever reconstructed filamentous fungal GRN primes our understanding of pathogenicity at the systems biology level and provides enticing prospects for novel disease control strategies involving the targeting of master regulators in pathogens. The program can be used to construct GRNs of other plant pathogens.

Dry and wet approaches for genome-wide functional annotation of conventional and unconventional transcriptional activators

Transcription factors (TFs) are master gene products that regulate gene expression in response to a variety of stimuli. They interact with DNA in a sequence-specific manner using a variety of DNA-binding domain (DBD) modules. This allows to properly position their second domain, called "effector domain", to directly or indirectly recruit positively or negatively acting co-regulators including chromatin modifiers, thus modulating preinitiation complex formation as well as transcription elongation. At variance with the DBDs, which are comprised of well-defined and easily recognizable DNA binding motifs, effector domains are usually much less conserved and thus considerably more difficult to predict. Also not so easy to identify are the DNA-binding sites of TFs, especially on a genome-wide basis and in the case of overlapping binding regions. Another emerging issue, with many potential regulatory implications, is that of so-called "moonlighting" transcription factors, i.e., proteins with an annotated function unrelated to transcription and lacking any recognizable DBD or effector domain, that play a role in gene regulation as their second job. Starting from bioinformatic and experimental high-throughput tools for an unbiased, genome-wide identification and functional characterization of TFs (especially transcriptional activators), we describe both established (and usually well affordable) as well as newly developed platforms for DNA-binding site identification. Selected combinations of these search tools, some of which rely on next-generation sequencing approaches, allow delineating the entire repertoire of TFs and unconventional regulators encoded by the any sequenced genome.

De Novo Sequencing and Transcriptome Analysis of Pleurotus eryngii subsp. tuoliensis (Bailinggu) Mycelia in Response to Cold Stimulation

Cold stimulation of Bailinggu's mycelia is the main factor that triggers primordia initiation for successful production of fruiting bodies under commercial cultivation. Yet, the molecular-level mechanisms involved in mycelia response to cold stimulation are still unclear. Here, we performed comparative transcriptomic analysis using RNA-Seq technology to better understand the gene expression regulation during different temporal stages of cold stimulation in Bailinggu. A total of 21,558 Bailinggu mycelia unigenes were de novo assembled and annotated from four libraries (control at 25 °C, plus cold stimulation treatments at -3 °C for a duration of 1-2 days, 5-6 days, and 9-10 days). GO and KEGG pathway analysis indicated that functional groups of differentially expressed unigenes associated with cell wall and membrane stabilization, calcium signaling and mitogen-activated protein kinases (MAPK) pathways, and soluble sugars and protein biosynthesis and metabolism pathways play a vital role in Bailinggu's response to cold stimulation. Six hundred and seven potential EST-based SSRs loci were identified in these unigenes, and 100 EST-SSR primers were randomly selected for validation. The overall polymorphism rate was 92% by using 10 wild strains of Bailinggu. Therefore, these results can serve as a valuable resource for a better understanding of the molecular mechanisms associated with Bailinggu's response to cold stimulation.

The Genomes of Three Uneven Siblings: Footprints of the Lifestyles of Three Trichoderma Species

The genus Trichoderma contains fungi with high relevance for humans, with applications in enzyme production for plant cell wall degradation and use in biocontrol. Here, we provide a broad, comprehensive overview of the genomic content of these species for "hot topic" research aspects, including CAZymes, transport, transcription factors, and development, along with a detailed analysis and annotation of less-studied topics, such as signal transduction, genome integrity, chromatin, photobiology, or lipid, sulfur, and nitrogen metabolism in T. reesei, T. atroviride, and T. virens, and we open up new perspectives to those topics discussed previously. In total, we covered more than 2,000 of the predicted 9,000 to 11,000 genes of each Trichoderma species discussed, which is >20% of the respective gene content. Additionally, we considered available transcriptome data for the annotated genes. Highlights of our analyses include overall carbohydrate cleavage preferences due to the different genomic contents and regulation of the respective genes. We found light regulation of many sulfur metabolic genes. Additionally, a new Golgi 1,2-mannosidase likely involved in N-linked glycosylation was detected, as were indications for the ability of Trichoderma spp. to generate hybrid galactose-containing N-linked glycans. The genomic inventory of effector proteins revealed numerous compounds unique to Trichoderma, and these warrant further investigation. We found interesting expansions in the Trichoderma genus in several signaling pathways, such as G-protein-coupled receptors, RAS GTPases, and casein kinases. A particularly interesting feature absolutely unique to T. atroviride is the duplication of the alternative sulfur amino acid synthesis pathway.

A Survey of the Gene Repertoire of Gigaspora rosea Unravels Conserved Features among Glomeromycota for Obligate Biotrophy

Arbuscular mycorrhizal (AM) fungi are a diverse group of soil fungi (Glomeromycota) that form the most ancient mutualistic association termed AM symbiosis with a majority of land plants, improving their nutrition uptake and resistance to stresses. In contrast to their great ecological implications, the knowledge of the molecular biological mechanisms involved is still scant, partly due to the limited genomic resources available. Here, we describe the gene repertoire of a new AM fungus Gigaspora rosea (Diversisporales). Among the 86332 non-redundant virtual transcripts assembled, 15346 presented similarities with proteins in the Refseq database and 10175 were assigned with GO terms. KOG and Interpro domain annotations clearly showed an enrichment of genes involved in signal transduction in G. rosea. KEGG pathway analysis indicates that most primary metabolic processes are active in G. rosea. However, as for Rhizophagus irregularis, several metabolic genes were not found, including the fatty acid synthase (FAS) gene. This finding supports the hypothesis that AM fungi depend on the lipids produced by their hosts. Furthermore, the presence of a large number of transporters and 100s of secreted proteins, together with the reduced number of plant cell wall degrading enzymes could be interpreted as an evolutionary adaptation to its mutualistic obligate biotrophy. The detection of meiosis-related genes suggests that G. rosea might use a cryptic sexual process. Lastly, a phylogeny of basal fungi clearly shows Glomeromycota as a sister clade to Mucoromycotina, not only to the Mucorales or Mortierellales. The characterization of the gene repertoire from an AM fungal species belonging to the order of Diversisporales and its comparison with the gene sets of R. irregularis (Glomerales) and Gigaspora margarita (Diversisporales), reveal that AM fungi share several features linked to mutualistic obligate biotrophy. This work contributes to lay the foundation for forthcoming studies into the genomics of Diversisporales, and also illuminates the utility of comparing gene repertoires of species from Diversisporales and other clades of Glomeromycota to gain more insights into the genetics and evolution of this fungal group.

Comparative Transcriptomic Analysis of Virulence Factors in Leptosphaeria maculans during Compatible and Incompatible Interactions with Canola

Leptosphaeria maculans is a hemibiotrophic fungus that causes blackleg of canola (Brassica napus), one of the most devastating diseases of this crop. In the present study, transcriptome profiling of L. maculans was performed in an effort to understand and define the pathogenicity genes that govern both the biotrophic and the necrotrophic phase of the fungus, as well as those that separate a compatible from an incompatible interaction. For this purpose, comparative RNA-seq analyses were performed on L. maculans isolate D5 at four different time points following inoculation on susceptible cultivar Topas-DH16516 or resistant introgression line Topas-Rlm2. Analysis of 1.6 billion Illumina reads readily identified differentially expressed genes that were over represented by candidate secretory effector proteins, CAZymes, and other pathogenicity genes. Comparisons between the compatible and incompatible interactions led to the identification of 28 effector proteins whose chronology and level of expression suggested a role in the establishment and maintenance of biotrophy with the plant. These included all known Avr genes of isolate D5 along with eight newly characterized effectors. In addition, another 15 effector proteins were found to be exclusively expressed during the necrotrophic phase of the fungus, which supports the concept that L. maculans has a separate and distinct arsenal contributing to each phase. As for CAZymes, they were often highly expressed at 3 dpi but with no difference in expression between the compatible and incompatible interactions, indicating that other factors were necessary to determine the outcome of the interaction. However, their significantly higher expression at 11 dpi in the compatible interaction confirmed that they contributed to the necrotrophic phase of the fungus. A notable exception was LysM genes whose high expression was singularly observed on the susceptible host at 7 dpi. In the case of TFs, their higher expression at 7 and 11 dpi on susceptible Topas support an important role in regulating the genes involved in the different pathogenic phases of L. maculans. In conclusion, comparison of the transcriptome of L. maculans during compatible and incompatible interactions has led to the identification of key pathogenicity genes that regulate not only the fate of the interaction but also lifestyle transitions of the fungus.

Synergistic and Dose-Controlled Regulation of Cellulase Gene Expression in Penicillium oxalicum

Filamentous fungus Penicillium oxalicum produces diverse lignocellulolytic enzymes, which are regulated by the combinations of many transcription factors. Here, a single-gene disruptant library for 470 transcription factors was constructed and systematically screened for cellulase production. Twenty transcription factors (including ClrB, CreA, XlnR, Ace1, AmyR, and 15 unknown proteins) were identified to play putative roles in the activation or repression of cellulase synthesis. Most of these regulators have not been characterized in any fungi before. We identified the ClrB, CreA, XlnR, and AmyR transcription factors as critical dose-dependent regulators of cellulase expression, the core regulons of which were identified by analyzing several transcriptomes and/or secretomes. Synergistic and additive modes of combinatorial control of each cellulase gene by these regulatory factors were achieved, and cellulase expression was fine-tuned in a proper and controlled manner. With one of these targets, the expression of the major intracellular β-glucosidase Bgl2 was found to be dependent on ClrB. The Bgl2-deficient background resulted in a substantial gene activation by ClrB and proved to be closely correlated with the relief of repression mediated by CreA and AmyR during cellulase induction. Our results also signify that probing the synergistic and dose-controlled regulation mechanisms of cellulolytic regulators and using it for reconstruction of expression regulation network (RERN) may be a promising strategy for cellulolytic fungi to develop enzyme hyper-producers. Based on our data, ClrB was identified as focal point for the synergistic activation regulation of cellulase expression by integrating cellulolytic regulators and their target genes, which refined our understanding of transcriptional-regulatory network as a “seesaw model” in which the coordinated regulation of cellulolytic genes is established by counteracting activators and repressors.

Cellulolytic fungi have evolved into sophisticated lignocellulolytic systems to adapt to their natural habitat. This trait is important for filamentous fungi, which are the main source of cellulases utilized to degrade lignocellulose to fermentable sugars. Penicillium oxalicum, which produces lignocellulolytic enzymes with more diverse components than Trichoderma reesei, has the capacity to secrete large amounts of cellulases. Meanwhile, cellulase expression is regulated by a complex network involved in many transcription factors in this organism. To better understand how cellulase genes are systematically regulated in P. oxalicum, we employed molecular genetics to uncover the cellulolytic transcription factors on a genome-wide scale. We discovered the synergistic and tunable regulation of cellulase expression by integrating cellulolytic regulators and their target genes, which refined our understanding of transcriptional-regulatory network as a “seesaw model” in which the coordinated regulation of cellulolytic genes is established by counteracting activators and repressors.

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.

Systematic functional profiling of transcription factor networks in Cryptococcus neoformans

Cryptococcus neoformans causes life-threatening meningoencephalitis in humans, but its overall biological and pathogenic regulatory circuits remain elusive, particularly due to the presence of an evolutionarily divergent set of transcription factors (TFs). Here, we report the construction of a high-quality library of 322 signature-tagged gene-deletion strains for 155 putative TF genes previously predicted using the DNA-binding domain TF database, and examine their in vitro and in vivo phenotypic traits under 32 distinct growth conditions. At least one phenotypic trait is exhibited by 145 out of 155 TF mutants (93%) and ∼85% of them (132/155) are functionally characterized for the first time in this study. The genotypic and phenotypic data for each TF are available in the C. neoformans TF phenome database (http://tf.cryptococcus.org). In conclusion, our phenome-based functional analysis of the C. neoformans TF mutant library provides key insights into transcriptional networks of basidiomycetous fungi and human fungal pathogens.

The Cell Factory Aspergillus Enters the Big Data Era: Opportunities and Challenges for Optimising Product Formation

Living with limits. Getting more from less. Producing commodities and high-value products from renewable resources including waste. What is the driving force and quintessence of bioeconomy outlines the lifestyle and product portfolio of Aspergillus, a saprophytic genus, to which some of the top-performing microbial cell factories belong: Aspergillus niger, Aspergillus oryzae and Aspergillus terreus. What makes them so interesting for exploitation in biotechnology and how can they help us to address key challenges of the twenty-first century? How can these strains become trimmed for better growth on second-generation feedstocks and how can we enlarge their product portfolio by genetic and metabolic engineering to get more from less? On the other hand, what makes it so challenging to deduce biological meaning from the wealth of Aspergillus -omics data? And which hurdles hinder us to model and engineer industrial strains for higher productivity and better rheological performance under industrial cultivation conditions? In this review, we will address these issues by highlighting most recent findings from the Aspergillus research with a focus on fungal growth, physiology, morphology and product formation. Indeed, the last years brought us many surprising insights into model and industrial strains. They clearly told us that similar is not the same: there are different ways to make a hypha, there are more protein secretion routes than anticipated and there are different molecular and physical mechanisms which control polar growth and the development of hyphal networks. We will discuss new conceptual frameworks derived from these insights and the future scientific advances necessary to create value from Aspergillus Big Data.

Systematic characterization of the bZIP transcription factor gene family in the rice blast fungus, Magnaporthe oryzae

Regulatory roles of the basic leucine zipper (bZIP) transcription factors (TFs) in fungi have been identified in diverse cellular processes such as development, nutrient utilization and various stress responses. In this study, the 22 Magnaporthe oryzae genes encoding bZIP TFs were systematically characterized. Phylogenetic analysis of fungal bZIP TFs revealed that seven MobZIPs are Magnaporthe-specific, while others belongs to 15 clades of orthologous Ascomycota genes. Expression patterns of MobZIPs under various conditions showed that they are highly stress responsive. We generated deletion mutants for 13 MobZIPs: nine with orthologues in other fungal species and four Magnaporthe-specific ones. Seven of them exhibited defects in mycelial growth, development and/or pathogenicity. Consistent with the conserved functions of the orthologues, MobZIP22 and MobZIP13 played a role in sulfur metabolism and iron homeostasis respectively. Along with MobZIP22 and MobZIP13, one Magnaporthe-specific gene, MobZIP11 is essential for pathogenicity in a reactive oxygen species-dependent manner. Taken together, our results will contribute to understanding the regulatory mechanisms of the bZIP TF gene family in fungal development, adaptation to environmental stresses and pathogenicity in the rice blast fungus.

The genome of the generalist plant pathogen Fusarium avenaceum is enriched with genes involved in redox, signaling and secondary metabolism

Fusarium avenaceum is a fungus commonly isolated from soil and associated with a wide range of host plants. We present here three genome sequences of F. avenaceum, one isolated from barley in Finland and two from spring and winter wheat in Canada. The sizes of the three genomes range from 41.6–43.1 MB, with 13217–13445 predicted protein-coding genes. Whole-genome analysis showed that the three genomes are highly syntenic, and share>95% gene orthologs. Comparative analysis to other sequenced Fusaria shows that F. avenaceum has a very large potential for producing secondary metabolites, with between 75 and 80 key enzymes belonging to the polyketide, non-ribosomal peptide, terpene, alkaloid and indole-diterpene synthase classes. In addition to known metabolites from F. avenaceum, fuscofusarin and JM-47 were detected for the first time in this species. Many protein families are expanded in F. avenaceum, such as transcription factors, and proteins involved in redox reactions and signal transduction, suggesting evolutionary adaptation to a diverse and cosmopolitan ecology. We found that 20% of all predicted proteins were considered to be secreted, supporting a life in the extracellular space during interaction with plant hosts.

Systematic analysis of Zn2Cys6 transcription factors required for development and pathogenicity by high-throughput gene knockout in the rice blast fungus

Because of great challenges and workload in deleting genes on a large scale, the functions of most genes in pathogenic fungi are still unclear. In this study, we developed a high-throughput gene knockout system using a novel yeast-Escherichia-Agrobacterium shuttle vector, pKO1B, in the rice blast fungus Magnaporthe oryzae. Using this method, we deleted 104 fungal-specific Zn(2)Cys(6) transcription factor (TF) genes in M. oryzae. We then analyzed the phenotypes of these mutants with regard to growth, asexual and infection-related development, pathogenesis, and 9 abiotic stresses. The resulting data provide new insights into how this rice pathogen of global significance regulates important traits in the infection cycle through Zn(2)Cys(6)TF genes. A large variation in biological functions of Zn(2)Cys(6)TF genes was observed under the conditions tested. Sixty-one of 104 Zn(2)Cys(6) TF genes were found to be required for fungal development. In-depth analysis of TF genes revealed that TF genes involved in pathogenicity frequently tend to function in multiple development stages, and disclosed many highly conserved but unidentified functional TF genes of importance in the fungal kingdom. We further found that the virulence-required TF genes GPF1 and CNF2 have similar regulation mechanisms in the gene expression involved in pathogenicity. These experimental validations clearly demonstrated the value of a high-throughput gene knockout system in understanding the biological functions of genes on a genome scale in fungi, and provided a solid foundation for elucidating the gene expression network that regulates the development and pathogenicity of M. oryzae.

In silico analysis of the core signaling proteome from the barley powdery mildew pathogen (Blumeria graminis f.sp. hordei)

BACKGROUND

Compared to other ascomycetes, the barley powdery mildew pathogen Blumeria graminis f.sp. hordei (Bgh) has a large genome (ca. 120 Mbp) that harbors a relatively small number of protein-coding genes (ca. 6500). This genomic assemblage is thought to be the result of numerous gene losses, which likely represent an evolutionary adaptation to a parasitic lifestyle in close association with its host plant, barley (Hordeum vulgare). Approximately 8% of the Bgh genes are predicted to encode virulence effectors that are secreted into host tissue and/or cells to promote pathogenesis; the remaining proteome is largely uncharacterized at present.

RESULTS

We provide a comparative analysis of the conceptual Bgh proteome, with an emphasis on proteins with known roles in fungal development and pathogenicity, for example heterotrimeric G proteins and G protein coupled receptors; components of calcium and cAMP signaling; small monomeric GTPases; mitogen-activated protein cascades and transcription factors. The predicted Bgh proteome lacks a number of proteins that are otherwise conserved in filamentous fungi, including two proteins that are required for the formation of anastomoses (somatic hyphal connections). By contrast, apart from minor modifications, all major canonical signaling pathways are retained in Bgh. A family of kinases that preferentially occur in pathogenic species of the fungal clade Leotiomyceta is unusually expanded in Bgh and its close relative, Blumeria graminis f.sp. tritici.

CONCLUSIONS

Our analysis reveals characteristic features of the proteome of a fungal phytopathogen that occupies an extreme habitat: the living plant cell.

Genome-wide identification and characterization of long intergenic non-coding RNAs in Ganoderma lucidum

Ganoderma lucidum is a white-rot fungus best-known for its medicinal activities. We have previously sequenced its genome and annotated the protein coding genes. However, long non-coding RNAs in G. lucidum genome have not been analyzed. In this study, we have identified and characterized long intergenic non-coding RNAs (lincRNA) in G. lucidum systematically. We developed a computational pipeline, which was used to analyze RNA-Seq data derived from G. lucidum samples collected from three developmental stages. A total of 402 lincRNA candidates were identified, with an average length of 609 bp. Analysis of their adjacent protein-coding genes (apcGenes) revealed that 46 apcGenes belong to the pathways of triterpenoid biosynthesis and lignin degradation, or families of cytochrome P450, mating type B genes, and carbohydrate-active enzymes. To determine if lincRNAs and these apcGenes have any interactions, the corresponding pairs of lincRNAs and apcGenes were analyzed in detail. We developed a modified 3' RACE method to analyze the transcriptional direction of a transcript. Among the 46 lincRNAs, 37 were found unidirectionally transcribed, and 9 were found bidirectionally transcribed. The expression profiles of 16 of these 37 lincRNAs were found to be highly correlated with those of the apcGenes across the three developmental stages. Among them, 11 are positively correlated (r>0.8) and 5 are negatively correlated (r<-0.8). The co-localization and co-expression of lincRNAs and those apcGenes playing important functions is consistent with the notion that lincRNAs might be important regulators for cellular processes. In summary, this represents the very first study to identify and characterize lincRNAs in the genomes of basidiomycetes. The results obtained here have laid the foundation for study of potential lincRNA-mediated expression regulation of genes in G. lucidum.

Roles of Forkhead-box Transcription Factors in Controlling Development, Pathogenicity, and Stress Response in Magnaporthe oryzae

Although multiple transcription factors (TFs) have been characterized via mutagenesis to understand their roles in controlling pathogenicity and infection-related development in Magnaporthe oryzae, the causal agent of rice blast, if and how forkhead-box (FOX) TFs contribute to these processes remain to be characterized. Four putative FOX TF genes were identified in the genome of M. oryzae, and phylogenetic analysis suggested that two of them (MoFKH1 and MoHCM1) correspond to Ascomycota-specific members of the FOX TF family while the others (MoFOX1 and MoFOX2) are Pezizomycotina-specific members. Deletion of MoFKH1 (ΔMofkh1) resulted in reduced mycelial growth and conidial germination, abnormal septation and stress response, and reduced virulence. Similarly, ΔMohcm1 exhibited reduced mycelial growth and conidial germination. Conidia of ΔMofkh1 and ΔMohcm1 were more sensitive to one or both of the cell cycle inhibitors hydroxyurea and benomyl, suggesting their role in cell cycle control. On the other hand, loss of MoFOX1 (ΔMofox1) did not show any noticeable changes in development, pathogenicity, and stress response. Deletion of MoFOX2 was not successful even after repeated attempts. Taken together, these results suggested that MoFKH1 and Mo-HCM1 are important in fungal development and that MoFKH1 is further implicated in pathogenicity and stress response in M. oryzae.