Gene discovery and gene expression in the rice blast fungus, Magnaporthe grisea: analysis of expressed sequence tags

Alternative splicing diversifies the transcriptome and proteome of the rice blast fungus during host infection

Alternative splicing (AS) contributes to diversifying and regulating cellular responses to environmental conditions and developmental cues by differentially producing multiple mRNA and protein isoforms from a single gene. Previous studies on AS in pathogenic fungi focused on profiling AS isoforms under a limited number of conditions. We analysed AS profiles in the rice blast fungus Magnaporthe oryzae, a global threat to rice production, using high-quality transcriptome data representing its vegetative growth (mycelia) and multiple host infection stages. We identified 4,270 AS isoforms derived from 2,413 genes, including 499 genes presumably regulated by infection-specific AS. AS appears to increase during infection, with 32.7% of the AS isoforms being produced during infection but absent in mycelia. Analysis of the isoforms observed at each infection stage showed that 636 AS isoforms were more abundant than corresponding annotated mRNAs, especially after initial hyphal penetration into host cell. Many such dominant isoforms were predicted to encode regulatory proteins such as transcription factors and phospho-transferases. We also identified the genes encoding distinct proteins via AS and confirmed the translation of some isoforms via a proteomic analysis, suggesting potential AS-mediated neo-functionalization of some genes during infection. Comprehensive profiling of the pattern of genome-wide AS during multiple stages of rice-M. oryzae interaction established a foundational resource that will help investigate the role and regulation of AS during rice infection.

Transcriptional Landscapes of Long Non-coding RNAs and Alternative Splicing in Pyricularia oryzae Revealed by RNA-Seq

Pyricularia oryzae causes the rice blast, which is one of the most devastating crop diseases worldwide, and is a model fungal pathogen widely used for dissecting the molecular mechanisms underlying fungal virulence/pathogenicity. Although the whole genome sequence of P. oryzae is publicly available, its current transcriptomes remain incomplete, lacking the information on non-protein coding genes and alternative splicing. Here, we performed and analyzed RNA-Seq of conidia and hyphae, resulting in the identification of 3,374 novel genes. Interestingly, the vast majority of these novel genes likely transcribed long non-coding RNAs (lncRNAs), and most of them were localized in the intergenic regions. Notably, their expressions were concomitant with the transcription of neighboring genes thereof in conidia and hyphae. In addition, 2,358 genes were found to undergo alternative splicing events. Furthermore, we exemplified that a lncRNA was important for hyphal growth likely by regulating the neighboring protein-coding gene and that alternative splicing of the transcription factor gene CON7 was required for appressorium formation. In summary, results from this study indicate that lncRNA transcripts and alternative splicing events are two important mechanisms for regulating the expression of genes important for conidiation, hyphal growth, and pathogenesis, and provide new insights into transcriptomes and gene regulation in the rice blast fungus.

Stress-Induced Changes in Alternative Splicing Landscape in Rice: Functional Significance of Splice Isoforms in Stress Tolerance

Improvements in yield and quality of rice are crucial for global food security. However, global rice production is substantially hindered by various biotic and abiotic stresses. Making further improvements in rice yield is a major challenge to the rice research community, which can be accomplished through developing abiotic stress-resilient rice varieties and engineering durable agrochemical-independent pathogen resistance in high-yielding elite rice varieties. This, in turn, needs increased understanding of the mechanisms by which stresses affect rice growth and development. Alternative splicing (AS), a post-transcriptional gene regulatory mechanism, allows rapid changes in the transcriptome and can generate novel regulatory mechanisms to confer plasticity to plant growth and development. Mounting evidence indicates that AS has a prominent role in regulating rice growth and development under stress conditions. Several regulatory and structural genes and splicing factors of rice undergo different types of stress-induced AS events, and the functional significance of some of them in stress tolerance has been defined. Both rice and its pathogens use this complex regulatory mechanism to devise strategies against each other. This review covers the current understanding and evidence for the involvement of AS in biotic and abiotic stress-responsive genes, and its relevance to rice growth and development. Furthermore, we discuss implications of AS for the virulence of different rice pathogens and highlight the areas of further research and potential future avenues to develop climate-smart and disease-resistant rice varieties.

Transcriptome analysis reveals the host selection fitness mechanisms of the Rhizoctonia solani AG1IA pathogen

Rhizoctonia solani AG1IA is a major generalist pathogen that causes sheath blight. Its genome, which was the first to be sequenced from the Rhizoctonia genus, may serve as a model for studying pathogenic mechanisms. To explore the pathogen-host fitness mechanism of sheath-blight fungus, a comprehensive comparative transcriptome ecotype analysis of R. solani AG1IA isolated from rice, soybean and corn during infection was performed. Special characteristics in gene expression, gene ontology terms and expression of pathogenesis-associated genes, including genes encoding secreted proteins, candidate effectors, hydrolases, and proteins involved in secondary metabolite production and the MAPK pathway, were revealed. Furthermore, as an important means of pathogenic modulation, diverse alternative splicing of key pathogenic genes in Rhizoctonia solani AG1IA during infections of the abovementioned hosts was uncovered for the first time. These important findings of key factors in the pathogenicity of R. solani AG1IA ecotypes during infection of various hosts explain host preference and provide novel insights into the pathogenic mechanisms and host-pathogen selection. Furthermore, they provide information on the fitness of Rhizoctonia, a severe pathogen with a wide host range.

Transcriptome analysis reveals the complexity of alternative splicing regulation in the fungus Verticillium dahliae

Background

Alternative splicing (AS) regulation is extensive and shapes the functional complexity of higher organisms. However, the contribution of alternative splicing to fungal biology is not well studied.

Results

This study provides sequences of the transcriptomes of the plant wilt pathogen Verticillium dahliae, using two different strains and multiple methods for cDNA library preparations. We identified alternatively spliced mRNA isoforms in over a half of the multi-exonic fungal genes. Over one-thousand isoforms involve TopHat novel splice junction; multiple types of combinatory alternative splicing patterns were identified. We showed that one Verticillium gene could use four different 5′ splice sites and two different 3′ donor sites to produce up to five mature mRNAs, representing one of the most sophisticated alternative splicing model in eukaryotes other than animals. Hundreds of novel intron types involving a pair of new splice sites were identified in the V. dahliae genome. All the types of AS events were validated by using RT-PCR. Functional enrichment analysis showed that AS genes are involved in most known biological functions and enriched in ATP biosynthesis, sexual/asexual reproduction, morphogenesis, signal transduction etc., predicting that the AS regulation modulates mRNA isoform output and shapes the V. dahliae proteome plasticity of the pathogen in response to the environmental and developmental changes.

Conclusions

These findings demonstrate the comprehensive alternative splicing mechanisms in a fungal plant pathogen, which argues the importance of this fungus in developing complicate genome regulation strategies in eukaryotes.

Electronic supplementary material

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

Transcriptome sequencing and analysis of the entomopathogenic fungus Hirsutella sinensis isolated from Ophiocordyceps sinensis

BACKGROUND

Ophiocordyceps sinensis, a worm and fungus combined mixture which Hirsutella sinensis is parasitic on the caterpillar body, has been used as a traditional medicine or healthy food in China for thousands of years. H. sinensis is reported as the only correct anamorph of O. sinensis and its main active ingredients are similar to the natural O. sinensis.

RESULTS

H. sinensis L0106, asexual strain of O. sinensis, was isolated and identified in this study. Three transcriptomes of H. sinensis at different cultivation periods (growth period 3d, pre-stable period 6d and stable period 9d) were sequenced for the first time by RNA-Seq method, and 25,511 unigenes (3d), 25,214 unigenes (6d) and 16,245 unigenes (9d) were assembled and obtained, respectively. These unigenes of the three samples were further assembled into 20,822 unigenes (All), and 62.3 percent of unigenes (All) could be annotated based on protein databases. Subsequently, the genes and enzymes involved in the biosynthesis of the active ingredients according to the sequencing and annotation results were predicted. Based on the predictions, we further investigated the interaction of different pathway networks and the corresponding enzymes. Furthermore, the differentially expressed genes (DEGs) of H. sinensis grown during different developmental stages (3d-VS-6d, 3d-VS-9d and 6d-VS-9d) were globally detected and analyzed based on the data from RNA-Seq, and 764 DEGs between 3d and 6d, 1,869 DEGs between 3d and 9d, and 770 DEGs between 6d and 9d were found, respectively.

CONCLUSIONS

This work presented here would aid in understanding and carrying out future studies on the genetic basis of H. sinensis and contribute to the further artificial production and application of this organism. This study provided a substantial contribution and basis to further characterize the gene expression profiles of H. sinensis in the metabolic pathways of active ingredients.

Deep RNA sequencing reveals a high frequency of alternative splicing events in the fungus Trichoderma longibrachiatum

BACKGROUND

Alternative splicing is crucial for proteome diversity and functional complexity in higher organisms. However, the alternative splicing landscape in fungi is still elusive.

RESULTS

The transcriptome of the filamentous fungus Trichoderma longibrachiatum was deep sequenced using Illumina Solexa technology. A total of 14305 splice junctions were discovered. Analyses of alternative splicing events revealed that the number of all alternative splicing events (10034), intron retentions (IR, 9369), alternative 5' splice sites (A5SS, 167), and alternative 3' splice sites (A3SS, 302) is 7.3, 7.4, 5.1, and 5.9-fold higher, respectively, than those observed in the fungus Aspergillus oryzae using Illumina Solexa technology. This unexpectedly high ratio of alternative splicing suggests that alternative splicing is important to the transcriptome diversity of T. longibrachiatum. Alternatively spliced introns had longer lengths, higher GC contents, and lower splice site scores than constitutive introns. Further analysis demonstrated that the isoform relative frequencies were correlated with the splice site scores of the isoforms. Moreover, comparative transcriptomics determined that most enzymes related to glycolysis and the citrate cycle and glyoxylate cycle as well as a few carbohydrate-active enzymes are transcriptionally regulated.

CONCLUSIONS

This study, consisting of a comprehensive analysis of the alternative splicing landscape in the filamentous fungus T. longibrachiatum, revealed an unexpectedly high ratio of alternative splicing events and provided new insights into transcriptome diversity in fungi.

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.

New gene models and alternative splicing in the maize pathogen Colletotrichum graminicola revealed by RNA-Seq analysis

BACKGROUND

An annotated genomic sequence of the corn anthracnose fungus Colletotrichum graminicola has been published previously, but correct identification of gene models by means of automated gene annotation remains a challenge. RNA-Seq offers the potential for substantially improved gene annotations and for the identification of posttranscriptional RNA modifications, such as alternative splicing and RNA editing.

RESULTS

Based on the nucleotide sequence information of transcripts, we identified 819 novel transcriptionally active regions (nTARs) and revised 906 incorrectly predicted gene models, including revisions of exon-intron structure, gene orientation and sequencing errors. Among the nTARs, 146 share significant similarity with proteins that have been identified in other species suggesting that they are hitherto unidentified genes in C. graminicola. Moreover, 5'- and 3'-UTR sequences of 4378 genes have been retrieved and alternatively spliced variants of 69 genes have been identified. Comparative analysis of RNA-Seq data and the genome sequence did not provide evidence for RNA editing in C. graminicola.

CONCLUSIONS

We successfully employed deep sequencing RNA-Seq data in combination with an elaborate bioinformatics strategy in order to identify novel genes, incorrect gene models and mechanisms of transcript processing in the corn anthracnose fungus C. graminicola. Sequence data of the revised genome annotation including several hundreds of novel transcripts, improved gene models and candidate genes for alternative splicing have been made accessible in a comprehensive database. Our results significantly contribute to both routine laboratory experiments and large-scale genomics or transcriptomic studies in C. graminicola.

Fungal alternative splicing is associated with multicellular complexity and virulence: a genome-wide multi-species study

Alternative splicing (AS) is a cellular process that increases a cell's coding capacity from a limited set of genes. Although AS is common in higher plants and animals, its prevalence in other eukaryotes is mostly unknown. In fungi the involvement of AS in gene expression and its effect on multi-cellularity and virulence is of great medical and economic interest. We present a genome-wide comparative study of AS in 23 informative fungi of different taxa, based on alignments of public transcript sequences. Random sampling of expressed sequence tags allows for robust and comparable estimations of AS rates. We find that a greater fraction of fungal genes than previously expected is associated with AS. We estimate that on average, 6.4% of the annotated genes are affected by AS, with Cryptococcus neoformans showing an extraordinary rate of 18%. The investigated Basidiomycota show higher average AS rates (8.6%) than the Ascomycota (6.0%), although not significant. We find that multi-cellular complexity and younger evolutionary age associate with higher AS rates. Furthermore, AS affects genes involved in pathogenic lifestyle, particularly in functions of stress response and dimorphic switching. Together, our analysis strongly supports the view that AS is a rather common phenomenon in fungi and associates with higher multi-cellular complexity.

Identification and characterization of in planta-expressed secreted effector proteins from Magnaporthe oryzae that induce cell death in rice

Interactions between rice and Magnaporthe oryzae involve the recognition of cellular components and the exchange of complex molecular signals from both partners. How these interactions occur in rice cells is still elusive. We employed robust-long serial analysis of gene expression, massively parallel signature sequencing, and sequencing by synthesis to examine transcriptome profiles of infected rice leaves. A total of 6,413 in planta-expressed fungal genes, including 851 genes encoding predicted effector proteins, were identified. We used a protoplast transient expression system to assess 42 of the predicted effector proteins for the ability to induce plant cell death. Ectopic expression assays identified five novel effectors that induced host cell death only when they contained the signal peptide for secretion to the extracellular space. Four of them induced cell death in Nicotiana benthamiana. Although the five effectors are highly diverse in their sequences, the physiological basis of cell death induced by each was similar. This study demonstrates that our integrative genomic approach is effective for the identification of in planta-expressed cell death-inducing effectors from M. oryzae that may play an important role facilitating colonization and fungal growth during infection.

Transcriptional profiling of ESTs from the biocontrol fungus Chaetomium cupreum

Comparative analysis was applied to two cDNA/ESTs libraries (C1 and C2) from Chaetomium cupreum. A total of 5538 ESTs were sequenced and assembled into 2162 unigenes including 585 contigs and 1577 singletons. BlastX analysis enabled the identification of 1211 unigenes with similarities to sequences in the public databases. MFS monosaccharide transporter was found as the gene expressed at the highest level in library C2, but no expression in C1. The majority of unigenes were library specific. Comparative analysis of the ESTs further revealed the difference of C. cupreum in gene expression and metabolic pathways between libraries. Two different sequences similar to the 48-KDa endochitinase and 46-KDa endochitinase were identified in libraries C1 and C2, respectively.

MGOS: development of a community annotation database for Magnaporthe oryzae

Magnaporthe oryzae causes rice blast disease, which is the most serious disease of cultivated rice worldwide. We previously developed the Magnaporthe grisea-Orzya sativa (MGOS) database as a repository for the M. oryzae and rice genome sequences together with a comprehensive set of functional interaction data generated by a major consortium of U.S. researchers. The MGOS database has now undergone a major redesign to include data from the international blast research community, accessible with a new intuitive, easy-to-use interface. Registered database users can manually annotate gene sequences and features as well as add mutant data and literature on individual gene pages. Over 900 genes have been manually curated based on various biological databases and the scientific literature. Gene names and descriptions, gene ontology annotations, published and unpublished information on mutants and their phenotypes, responses in diverse microarray analyses, and related literature have been incorporated. Thus far, 362 M. oryzae genes have associated information on mutants. MGOS is now poised to become a one-stop repository for all structural and functional data available on all genes of this critically important rice pathogen.

A Pmk1-interacting gene is involved in appressorium differentiation and plant infection in Magnaporthe oryzae

In the rice blast fungus Magnaporthe oryzae, the PMK1 mitogen-activated protein (MAP) kinase gene regulates appressorium formation and infectious growth. Its homologs in many other fungi also play critical roles in fungal development and pathogenicity. However, the targets of this important MAP kinase and its interacting genes are not well characterized. In this study, we constructed two yeast two-hybrid libraries of M. oryzae and screened for Pmk1-interacting proteins. Among the nine Pmk1-interacting clones (PICs) identified, two of them, PIC1 and PIC5, were selected for further characterization. Pic1 has one putative nuclear localization signal and one putative MAP kinase phosphorylation site. Pic5 contains one transmembrane domain and two functionally unknown CTNS (cystinosin/ERS1p repeat) motifs. The interaction of Pmk1 with Pic1 or Pic5 was confirmed by coimmunoprecipitation assays. Targeted gene deletion of PIC1 had no apparent effects on vegetative growth and pathogenicity but resulted in a significant reduction in conidiation and abnormal germ tube differentiation on onion epidermal cells. Deletion of PIC5 led to a reduction in conidiation and hyphal growth. Autolysis of aerial hyphae became visible in cultures older than 4 days. The pic5 mutant was defective in germ tube growth and appressorium differentiation. It was reduced in appressorial penetration and virulence on the plant. Both PIC1 and PIC5 are conserved in filamentous ascomycetes, but none of their orthologs have been functionally characterized. Our data indicate that PIC5 is a novel virulence factor involved in appressorium differentiation and pathogenesis in M. oryzae.

Application of high-throughput DNA sequencing in phytopathology

The new sequencing technologies are already making a big impact in academic research on medically important microbes and may soon revolutionize diagnostics, epidemiology, and infection control. Plant pathology also stands to gain from exploiting these opportunities. This manuscript reviews some applications of these high-throughput sequencing methods that are relevant to phytopathology, with emphasis on the associated computational and bioinformatics challenges and their solutions. Second-generation sequencing technologies have recently been exploited in genomics of both prokaryotic and eukaryotic plant pathogens. They are also proving to be useful in diagnostics, especially with respect to viruses.

Genome-wide characterization of methylguanosine-capped and polyadenylated small RNAs in the rice blast fungus Magnaporthe oryzae

Small RNAs are well described in higher eukaryotes such as mammals and plants; however, knowledge in simple eukaryotes such as filamentous fungi is limited. In this study, we discovered and characterized methylguanosine-capped and polyadenylated small RNAs (CPA-sRNAs) by using differential RNA selection, full-length cDNA cloning and 454 transcriptome sequencing of the rice blast fungus Magnaporthe oryzae. This fungus causes blast, a devastating disease on rice, the principle food staple for over half the world’s population. CPA-sRNAs mapped primarily to the transcription initiation and termination sites of protein-coding genes and were positively correlated with gene expression, particularly for highly expressed genes including those encoding ribosomal proteins. Numerous CPA-sRNAs also mapped to rRNAs, tRNAs, snRNAs, transposable elements and intergenic regions. Many other 454 sequence reads could not be mapped to the genome; however, inspection revealed evidence for non-template additions and chimeric sequences. CPA-sRNAs were independently confirmed using a high affinity variant of eIF-4E to capture 5′-methylguanosine-capped RNA followed by 3′-RACE sequencing. These results expand the repertoire of small RNAs in filamentous fungi.

Detection and analysis of alternative splicing in Yarrowia lipolytica reveal structural constraints facilitating nonsense-mediated decay of intron-retaining transcripts

Background

Hemiascomycetous yeasts have intron-poor genomes with very few cases of alternative splicing. Most of the reported examples result from intron retention in Saccharomyces cerevisiae and some have been shown to be functionally significant. Here we used transcriptome-wide approaches to evaluate the mechanisms underlying the generation of alternative transcripts in Yarrowia lipolytica, a yeast highly divergent from S. cerevisiae.

Results

Experimental investigation of Y. lipolytica gene models identified several cases of alternative splicing, mostly generated by intron retention, principally affecting the first intron of the gene. The retention of introns almost invariably creates a premature termination codon, as a direct consequence of the structure of intron boundaries. An analysis of Y. lipolytica introns revealed that introns of multiples of three nucleotides in length, particularly those without stop codons, were underrepresented. In other organisms, premature termination codon-containing transcripts are targeted for degradation by the nonsense-mediated mRNA decay (NMD) machinery. In Y. lipolytica, homologs of S. cerevisiae UPF1 and UPF2 genes were identified, but not UPF3. The inactivation of Y. lipolytica UPF1 and UPF2 resulted in the accumulation of unspliced transcripts of a test set of genes.

Conclusions

Y. lipolytica is the hemiascomycete with the most intron-rich genome sequenced to date, and it has several unusual genes with large introns or alternative transcription start sites, or introns in the 5' UTR. Our results suggest Y. lipolytica intron structure is subject to significant constraints, leading to the under-representation of stop-free introns. Consequently, intron-containing transcripts are degraded by a functional NMD pathway.

Survey of the transcriptome of Aspergillus oryzae via massively parallel mRNA sequencing

Aspergillus oryzae, an important filamentous fungus used in food fermentation and the enzyme industry, has been shown through genome sequencing and various other tools to have prominent features in its genomic composition. However, the functional complexity of the A. oryzae transcriptome has not yet been fully elucidated. Here, we applied direct high-throughput paired-end RNA-sequencing (RNA-Seq) to the transcriptome of A. oryzae under four different culture conditions. With the high resolution and sensitivity afforded by RNA-Seq, we were able to identify a substantial number of novel transcripts, new exons, untranslated regions, alternative upstream initiation codons and upstream open reading frames, which provide remarkable insight into the A. oryzae transcriptome. We were also able to assess the alternative mRNA isoforms in A. oryzae and found a large number of genes undergoing alternative splicing. Many genes and pathways that might be involved in higher levels of protein production in solid-state culture than in liquid culture were identified by comparing gene expression levels between different cultures. Our analysis indicated that the transcriptome of A. oryzae is much more complex than previously anticipated, and these results may provide a blueprint for further study of the A. oryzae transcriptome.

Identification and analysis of in planta expressed genes of Magnaporthe oryzae

Background

Infection of plants by pathogens and the subsequent disease development involves substantial changes in the biochemistry and physiology of both partners. Analysis of genes that are expressed during these interactions represents a powerful strategy to obtain insights into the molecular events underlying these changes. We have employed expressed sequence tag (EST) analysis to identify rice genes involved in defense responses against infection by the blast fungus Magnaporthe oryzae and fungal genes involved in infectious growth within the host during a compatible interaction.

Results

A cDNA library was constructed with RNA from rice leaves (Oryza sativa cv. Hwacheong) infected with M. oryzae strain KJ201. To enrich for fungal genes, subtraction library using PCR-based suppression subtractive hybridization was constructed with RNA from infected rice leaves as a tester and that from uninfected rice leaves as the driver. A total of 4,148 clones from two libraries were sequenced to generate 2,302 non-redundant ESTs. Of these, 712 and 1,562 ESTs could be identified to encode fungal and rice genes, respectively. To predict gene function, Gene Ontology (GO) analysis was applied, with 31% and 32% of rice and fungal ESTs being assigned to GO terms, respectively. One hundred uniESTs were found to be specific to fungal infection EST. More than 80 full-length fungal cDNA sequences were used to validate ab initio annotated gene model of M. oryzae genome sequence.

Conclusion

This study shows the power of ESTs to refine genome annotation and functional characterization. Results of this work have advanced our understanding of the molecular mechanisms underpinning fungal-plant interactions and formed the basis for new hypothesis.

Generation and analysis of expression sequence tags from haustoria of the wheat stripe rust fungus Puccinia striiformis f. sp. Tritici

BACKGROUND

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most destructive diseases of wheat (Triticum aestivum L.) worldwide. In spite of its agricultural importance, the genomics and genetics of the pathogen are poorly characterized. Pst transcripts from urediniospores and germinated urediniospores have been examined previously, but little is known about genes expressed during host infection. Some genes involved in virulence in other rust fungi have been found to be specifically expressed in haustoria. Therefore, the objective of this study was to generate a cDNA library to characterize genes expressed in haustoria of Pst.

RESULTS

A total of 5,126 EST sequences of high quality were generated from haustoria of Pst, from which 287 contigs and 847 singletons were derived. Approximately 10% and 26% of the 1,134 unique sequences were homologous to proteins with known functions and hypothetical proteins, respectively. The remaining 64% of the unique sequences had no significant similarities in GenBank. Fifteen genes were predicted to be proteins secreted from Pst haustoria. Analysis of ten genes, including six secreted protein genes, using quantitative RT-PCR revealed changes in transcript levels in different developmental and infection stages of the pathogen.

CONCLUSIONS

The haustorial cDNA library was useful in identifying genes of the stripe rust fungus expressed during the infection process. From the library, we identified 15 genes encoding putative secreted proteins and six genes induced during the infection process. These genes are candidates for further studies to determine their functions in wheat-Pst interactions.

Candidate effector gene identification in the ascomycete fungal phytopathogen Venturia inaequalis by expressed sequence tag analysis

The hemi-biotrophic fungus Venturia inaequalis infects members of the Maloideae, causing the economically important apple disease, scab. The plant-pathogen interaction of Malus and V. inaequalis follows the gene-for-gene model. cDNA libraries were constructed, and bioinformatic analysis of the resulting expressed sequence tags (ESTs) was used to characterize potential effector genes. Effectors are small proteins, secreted in planta, that are assumed to facilitate infection. Therefore, a cDNA library was constructed from a compatible interaction. To distinguish pathogen from plant sequences, the library was probed with genomic DNA from V. inaequalis to enrich for pathogen genes, and cDNA libraries were constructed from in vitro-grown material. A suppression subtractive hybridization library enriched for cellophane-induced genes was included, as growth on cellophane may mimic that in planta, with the differentiation of structures resembling those formed during plant colonization. Clustering of ESTs from the in planta and in vitro libraries indicated a fungal origin of the resulting non-redundant sequence. A total of 937 ESTs was classified as putatively fungal, which could be assembled into 633 non-redundant sequences. Sixteen new candidate effector genes were identified from V. inaequalis based on features common to characterized effector genes from filamentous fungi, i.e. they encode a small, novel, cysteine-rich protein, with a putative signal peptide. Three of the 16 candidates, in particular, conformed to most of the protein structural characteristics expected of fungal effectors and showed significant levels of transcriptional up-regulation during in planta growth. In addition to candidate effector genes, this collection of ESTs represents a valuable genomic resource for V. inaequalis.

Interaction transcriptome analysis identifies Magnaporthe oryzae BAS1-4 as Biotrophy-associated secreted proteins in rice blast disease

Biotrophic invasive hyphae (IH) of the blast fungus Magnaporthe oryzae secrete effectors to alter host defenses and cellular processes as they successively invade living rice (Oryza sativa) cells. However, few blast effectors have been identified. Indeed, understanding fungal and rice genes contributing to biotrophic invasion has been difficult because so few plant cells have encountered IH at the earliest infection stages. We developed a robust procedure for isolating infected-rice sheath RNAs in which approximately 20% of the RNA originated from IH in first-invaded cells. We analyzed these IH RNAs relative to control mycelial RNAs using M. oryzae oligoarrays. With a 10-fold differential expression threshold, we identified known effector PWL2 and 58 candidate effectors. Four of these candidates were confirmed to be fungal biotrophy-associated secreted (BAS) proteins. Fluorescently labeled BAS proteins were secreted into rice cells in distinct patterns in compatible, but not in incompatible, interactions. BAS1 and BAS2 proteins preferentially accumulated in biotrophic interfacial complexes along with known avirulence effectors, BAS3 showed additional localization near cell wall crossing points, and BAS4 uniformly outlined growing IH. Analysis of the same infected-tissue RNAs with rice oligoarrays identified putative effector-induced rice susceptibility genes, which are highly enriched for sensor-transduction components rather than typically identified defense response genes.

A Cdc42 ortholog is required for penetration and virulence of Magnaporthe grisea

Cdc42, a member of the Rho-family small GTP-binding proteins, is a pivotal signaling switch that cycles between active GTP-bound and inactive GDP-bound forms, controlling actin cytoskeleton organization and cell polarity. In this report, we show that MgCdc42, a Cdc42 ortholog in Magnaporthe grisea, is required for its plant penetration. Consequently, the deletion mutants show dramatically decreased virulence to rice due to the arrest of penetration and infectious growth, which may be attributed to the defect of turgor and superoxide generation during the appressorial development in Mgcdc42 deletion mutants. In addition, the mutants also exhibit pleotropic defects including gherkin-shaped conidia, delayed germination as well as decreased sporulation. Furthermore, dominant negative mutation leads to a similar phenotype to that of the deletion mutants, lending further support to the conclusion that MgCdc42 is required for the penetration and virulence of M. grisea.

Genome-wide validation of Magnaporthe grisea gene structures based on transcription evidence

Accurate cDNA data is useful to validate gene structures in a genome. We sequenced 35189 expressed sequence tags (ESTs) obtained from the highly destructive rice blast fungus, Magnaporthe grisea. Our custom-made computational programs mapped these ESTs on the M. grisea genome sequence, and reconstructed gene structures as well as protein-coding regions. As a result, we predicted 4480 protein-coding sequences, which were more accurate than ab initio predictions. Moreover, cross-species comparisons suggested that our predicted proteins were nearly complete. The cDNA clones obtained in this study will be important for further experimental studies. Our genome annotation is available at http://www.mg.dna.affrc.go.jp/.

Comparative EST analysis of a Zoophthora radicans isolate derived from Pieris brassicae and an isogenic strain adapted to Plutella xylostella

Zoophthora radicansis an entomopathogenic fungus with the potential to be used as an insect biological control agent. To better understand the mechanisms used byZ. radicansto infect different hosts, we generated expressed sequence tag (EST) datasets from aZ. radicansstrain originally isolated fromPieris brassicae, and an isogenic strain passaged throughPlutella xylostella. In total, 1839 ESTs were generated which clustered into 466 contigs and 433 singletons to provide a set of 899 unique sequences. Approximately 85 % of the ESTs were significantly similar (E≤e−03) to other fungal genes, of which 69.6 % encoded proteins with a reported function. Proteins involved in protein synthesis and metabolism were encoded by 38.3 % of the ESTs, while 26.3 % encoded proteins involved in cell-cycle regulation, DNA synthesis, protein fate, transport, cell defence, transcription and RNA synthesis, and 4.9 % encoded proteins associated with cellular transport, signal transduction, control of cellular organization and cell-wall degradation. Several proteinases, including aspartic proteinases, trypsins, trypsin-like serine proteases and metalloproteases, with the potential to degrade insect cuticle were expressed by the two isolates.

Large-scale gene discovery in the septoria tritici blotch fungus Mycosphaerella graminicola with a focus on in planta expression

The foliar disease septoria tritici blotch, caused by the fungus Mycosphaerella graminicola, is currently the most important wheat disease in Europe. Gene expression was examined under highly different conditions, using 10 expressed sequence tag libraries generated from M. graminicola isolate IPO323 using seven in vitro and three in planta growth conditions. To identify fungal clones in the interaction libraries, we developed a selection method based on hybridization with the entire genomic DNA of M. graminicola, to selectively enrich these libraries for fungal genes. Assembly of the 27,007 expressed sequence tags resulted in 9,190 unigenes, representing 5.2 Mb of the estimated 39-Mb genome size of M. graminicola. All libraries contributed significantly to the number of unigenes, especially the in planta libraries representing different stages of pathogenesis, which covered 15% of the library-specific unigenes. Even under presymptomatic conditions (5 days postinoculation), when fungal biomass is less than 5%, this method enabled us to efficiently capture fungal genes expressed during pathogenesis. Many of these genes were uniquely expressed in planta, indicating that in planta gene expression significantly differed from in vitro expression. Examples of gene discovery included a number of cell wall-degrading enzymes, a broad set of genes involved in signal transduction (n=11) and a range of ATP-binding cassette (n=20) and major facilitator superfamily transporter genes (n=12) potentially involved in protection against antifungal compounds or the secretion of pathogenicity factors. In addition, evidence is provided for a mycovirus in M. graminicola that is highly expressed under various stress conditions, in particular, under nitrogen starvation. Our analyses provide a unique window on in vitro and in planta gene expression of M. graminicola.

Stage-specific gene expression during urediniospore germination in Puccinia striiformis f. sp tritici

Background

Puccinia striiformis f. sp. tritici is an obligate biotrophic pathogen that causes leaf stripe rust on wheat. Although it is critical to understand molecular mechanisms of pathogenesis in the wheat stripe rust fungus for developing novel disease management strategies, little is known about its genome and gene functions due to difficulties in molecular studies with this important pathogen. To identify genes expressed during early infection stages, in this study we constructed a cDNA library with RNA isolated from urediniospores of P. striiformis f. sp. tritici germinated for 10 h.

Results

A total of 4798 ESTs were sequenced from the germinated urediniospore library and assembled into 315 contigs and 803 singletons. About 23.9% and 13.3% of the resulting 1118 unisequences were homologous to functionally characterized proteins and hypothetical proteins, respectively. The rest 62.8% unisequences had no significant homologs in GenBank. Several of these ESTs shared significant homology with known fungal pathogenicity or virulence factors, such as HESP767 of the flax rust and PMK1, GAS1, and GAS2 of the rice blast fungus. We selected six ESTs (Ps28, Ps85, Ps87, Ps259, Ps261, and Ps159) for assaying their expression patterns during urediniospore germination and wheat infection by quantitative real-time PCR. All of them had the highest transcript level in germinated urediniospores and a much less transcript level in un-germinated urediniospores and infected wheat tissues (1–7 dpi). The transcript level of Ps159 increased at later infection stages (6–7 dpi). Our data indicated that these genes were highly expressed in germinated urediniospores and may play important roles in fungal-plant interactions during early infection stages in the wheat stripe rust fungus.

Conclusion

Genes expressed in germinated urediniospores of P. striiformis f. sp. tritici were identified by EST analysis. Six of them were confirmed by quantitative real-time PCR assays to be highly expressed in germinated urediniospores.

Alternative splicing: current perspectives

Alternative splicing is a well-characterized mechanism by which multiple transcripts are generated from a single mRNA precursor. By allowing production of several protein isoforms from one pre-mRNA, alternative splicing contributes to proteomic diversity. But what do we know about the origin of this mechanism? Do the same evolutionary forces apply to alternatively and constitutively splice exons? Do similar forces act on all types of alternative splicing? Are the products generated by alternative splicing functional? Why is "improper" recognition of exons and introns allowed by the splicing machinery? In this review, we summarize the current knowledge regarding these issues from an evolutionary perspective.

MADS-box transcription factor mig1 is required for infectious growth in Magnaporthe grisea

Magnaporthe grisea is a model fungus for studying fungus-plant interactions. Two mitogen-activated protein (MAP) kinase genes, PMK1 and MPS1, have been implicated in regulating plant infection processes in M. grisea. However, transcription factors activated by these MAP kinases are not well studied. In this study we functionally characterized the MIG1 gene that encodes a MADS-box transcription factor homologous to Saccharomyces cerevisiae Rlm1. In yeast two-hybrid assays, MIG1 interacts with MPS1, suggesting that MIG1 may function downstream from the MPS1 pathway. The mig1 deletion mutant had a normal growth rate and formed melanized appressoria, but it was nonpathogenic and failed to infect rice leaves through wounds. Appressoria formed by the mig1 mutant developed penetration pegs and primary infectious hyphae, but further differentiation of the secondary infectious hyphae inside live plant cells was blocked. However, the mig1 mutant formed infectious hypha-like structures in heat-killed plant cells or cellophane membranes. In transformants expressing the MIG1-GFP fusion, green fluorescent protein (GFP) signals were not detectable in vegetative hyphae and conidiophores. Mig1-GFP was localized to nuclei in conidia, appressoria, and infectious hyphae. Deletion of the MADS box had no effect on the expression and localization of the MIG1-GFP fusion but eliminated its ability to complement the mig1 mutant. These results suggest that MIG1 may be required for overcoming plant defense responses and the differentiation of secondary infectious hyphae in live plant cells. The MADS-box domain is essential for the function of MIG1 but dispensable for its nuclear localization, which may be associated with the activation of MIG1 by MPS1 during conidiation and plant infection.

Gene discovery and transcript analyses in the corn smut pathogen Ustilago maydis: expressed sequence tag and genome sequence comparison

Background

Ustilago maydis is the basidiomycete fungus responsible for common smut of corn and is a model organism for the study of fungal phytopathogenesis. To aid in the annotation of the genome sequence of this organism, several expressed sequence tag (EST) libraries were generated from a variety of U. maydis cell types. In addition to utility in the context of gene identification and structure annotation, the ESTs were analyzed to identify differentially abundant transcripts and to detect evidence of alternative splicing and anti-sense transcription.

Results

Four cDNA libraries were constructed using RNA isolated from U. maydis diploid teliospores (U. maydis strains 518 × 521) and haploid cells of strain 521 grown under nutrient rich, carbon starved, and nitrogen starved conditions. Using the genome sequence as a scaffold, the 15,901 ESTs were assembled into 6,101 contiguous expressed sequences (contigs); among these, 5,482 corresponded to predicted genes in the MUMDB (MIPS Ustilago maydis database), while 619 aligned to regions of the genome not yet designated as genes in MUMDB. A comparison of EST abundance identified numerous genes that may be regulated in a cell type or starvation-specific manner. The transcriptional response to nitrogen starvation was assessed using RT-qPCR. The results of this suggest that there may be cross-talk between the nitrogen and carbon signalling pathways in U. maydis. Bioinformatic analysis identified numerous examples of alternative splicing and anti-sense transcription. While intron retention was the predominant form of alternative splicing in U. maydis, other varieties were also evident (e.g. exon skipping). Selected instances of both alternative splicing and anti-sense transcription were independently confirmed using RT-PCR.

Conclusion

Through this work: 1) substantial sequence information has been provided for U. maydis genome annotation; 2) new genes were identified through the discovery of 619 contigs that had previously escaped annotation; 3) evidence is provided that suggests the regulation of nitrogen metabolism in U. maydis differs from that of other model fungi, and 4) Alternative splicing and anti-sense transcription were identified in U. maydis and, amid similar observations in other basidiomycetes, this suggests these phenomena may be widespread in this group of fungi. These advances emphasize the importance of EST analysis in genome annotation.

Alternate intron processing of family 5 endoglucanase transcripts from the genus Phytophthora

Twenty-one homologs of family 5 endo-(1-4)-beta-glucanase genes (EGLs) were identified and characterized in the oomycete plant pathogens Phytophthora infestans, P. sojae, and P. ramorum, providing the first comprehensive analysis of this family in Phytophthora. Phylogenetic analysis revealed that these genes constitute a unique eukaryotic group, with closest similarity to bacterial endoglucanases. Many of the identified EGL copies were clustered in a few genomic regions, and contained from zero to three introns. Using reverse transcription PCR to study in vitro and in planta gene expression levels of P. sojae, we detected partially processed RNA transcripts retaining one or more of their introns. In some cases, the positions of intron/exon splicing sites were also found to be variable. The relative proportions of these transcripts remain apparently unchanged under various growing conditions, but differ among orthologous copies of the three Phytophthora species. The alternate processing of introns in this group of EGLs generates both coding and non-coding RNA isoforms. This is the first report on Phytophthora family 5 endoglucanases, and the first record for alternative intron processing of oomycete transcripts.

Expressed sequence tags from phytophthora sojae reveal genes specific to development and infection

Six unique expressed sequence tag (EST) libraries were generated from four developmental stages of Phytophthora sojae P6497. RNA was extracted from mycelia, swimming zoospores, germinating cysts, and soybean (Glycine max (L.) Merr.) cv. Harosoy tissues heavily infected with P. sojae. Three libraries were created from mycelia growing on defined medium, complex medium, and nutrient-limited medium. The 26,943 high-quality sequences obtained clustered into 7,863 unigenes composed of 2,845 contigs and 5,018 singletons. The total number of P. sojae unigenes matching sequences in the genome assembly was 7,412 (94%). Of these unigenes, 7,088 (90%) matched gene models predicted from the P. sojae sequence assembly, but only 2,047 (26%) matched P. ramorum gene models. Analysis of EST frequency from different growth conditions and morphological stages revealed genes that were specific to or highly represented in particular growth conditions and life stages. Additionally, our results indicate that, during infection, the pathogen derives most of its carbon and energy via glycolysis of sugars in the plant. Sequences identified with putative roles in pathogenesis included avirulence homologs possessing the RxLR motif, elicitins, and hydrolytic enzymes. This large collection of P. sojae ESTs will serve as a valuable public genomic resource.

Generation of a wheat leaf rust, Puccinia triticina, EST database from stage-specific cDNA libraries

Thirteen cDNA libraries constructed from small amounts of leaf rust mRNA using optimized methods served as the source for the generation of 25 558 high-quality DNA sequence reads. Five life-cycle stages were sampled: resting urediniospores, urediniospores germinated over water or plant extract, compatible, interactive stages during appressorium or haustorium formation just before sporulation, and an incompatible interaction. mRNA populations were subjected to treatments such as full-length cDNA production, subtractive and normalizing hybridizations, and size selection methods combined with PCR amplification. Pathogen and host sequences from interactive libraries were differentiated in silico using cereal and fungal sequences, codon usage analyses, and by means of a partial prototype cDNA microarray hybridized with genomic DNAs. This yielded a non-redundant unigene set of 9760 putative fungal sequences consisting of 6616 singlets and 3144 contigs, representing 4.7 Mbp. At an E-value 10(-5), 3670 unigenes (38%) matched sequences in various databases and collections but only 694 unigenes (7%) were similar to genes with known functions. In total, 296 unigenes were identified as most probably wheat and ten as rRNA sequences. Annotation rates were low for germinated urediniospores (4%) and appressoria (2%). Gene sets obtained from the various life-cycle stages appear to be remarkably different, suggesting drastic reprogramming of the transcriptome during these major differentiation processes. Redundancy within contigs yielded information about possible expression levels of certain genes among stages. Many sequences were similar to genes from other rusts such as Uromyces and Melampsora species; some of these genes have been implicated in pathogenicity and virulence.

Mirl is highly upregulated and localized to nuclei during infectious hyphal growth in the rice blast fungus

Rice blast, caused by Magnaporthe grisea, is a devastating disease of rice throughout the world. Many recent molecular studies have focused on the early infection stages, but our knowledge about molecular events at the infectious hyphae stage is limited. In this study, 750 hygromycin-resistant transformants were isolated by transforming M. grisea Guyll with a promoterless enhanced green fluorescent protein (EGFP) construct. In one of the transformants, L1320, EGFP signals were observed in the nuclei of infectious hyphae. The transforming vector was inserted in a predicted gene named MIR1 and resulted in a Mir1 1-107-EGFP fusion. Mir1 is a low-complexity protein with no known protein domain and has no homolog in GenBank or other sequenced fungal genomes. Quantitative real-time reverse-transcriptase polymerase chain reaction analysis and expression assays of MIR1-EGFP fusion constructs indicated that the expression of MIR1 was highly induced during plant infection. Deletion analyses identified a 458-bp region that was sufficient for the MIR1 promoter activity. Further characterization revealed that a 96-bp sequence was essential for the enhanced in planta expression. MIR1 is an M. grisea-specific gene that is highly conserved among the field isolates belonging to the M. grisea species complex. The mir1 mutants had no obvious defects in appressorial penetration and rice infection. When overexpressed with the RP27 promoter, nuclear localization of the Mir1-EGFP fusion was observed in conidia and vegetative hyphae. These data suggest that the expression but not the nuclear localization of MIR1 is specific to infectious hyphae and that reporter genes based on MIR1 may be suitable for monitoring infectious growth in M. grisea.

Generation, annotation, and analysis of ESTs from four different Trichoderma strains grown under conditions related to biocontrol

The functional genomics project "TrichoEST" was developed focused on different taxonomic groups of Trichoderma with biocontrol potential. Four cDNA libraries were constructed, using similar growth conditions, from four different Trichoderma strains: Trichoderma longibrachiatum T52, Trichoderma asperellum T53, Trichoderma virens T59, and Trichoderma sp. T78. In this study, we present the analysis of the 8,160 expressed sequence tags (ESTs) generated. Each EST library was independently assembled and 1,000-1,300 unique sequences were identified in each strain. First, we queried our collection of ESTs against the NCBI nonredundant database using the BLASTX algorithm. Moreover, using the Gene Ontology hierarchy, we performed the annotation of 40.9% of the unique sequences. Later, based on the EST abundance, we examined the highly expressed genes in the four strains. A hydrophobin was found as the gene expressed at the highest level in two of the strains, but we also found that other unique sequences similar to the HEX1, QID3, and NMT1 proteins were highly represented in at least two of the Trichoderma strains.

Expression of Magnaporthe grisea avirulence gene ACE1 is connected to the initiation of appressorium-mediated penetration

Magnaporthe grisea is responsible for a devastating fungal disease of rice called blast. Current control of this disease relies on resistant rice cultivars that recognize M. grisea signals corresponding to specific secreted proteins encoded by avirulence genes. The M. grisea ACE1 avirulence gene differs from others, since it controls the biosynthesis of a secondary metabolite likely recognized by rice cultivars carrying the Pi33 resistance gene. Using a transcriptional fusion between ACE1 promoter and eGFP, we showed that ACE1 is only expressed in appressoria during fungal penetration into rice and barley leaves, onion skin, and cellophane membranes. ACE1 is almost not expressed in appressoria differentiated on Teflon and Mylar artificial membranes. ACE1 expression is not induced by cellophane and plant cell wall components, demonstrating that it does not require typical host plant compounds. Cyclic AMP (cAMP) signaling mutants delta cpkA and delta mac1 sum1-99 and tetraspanin mutant delta pls1::hph differentiate melanized appressoria with normal turgor but are unable to penetrate host plant leaves. ACE1 is normally expressed in these mutants, suggesting that it does not require cAMP signaling or a successful penetration event. ACE1 is not expressed in appressoria of the buf1::hph mutant defective for melanin biosynthesis and appressorial turgor. The addition of hyperosmotic solutes to buf1::hph appressoria restores appressorial development and ACE1 expression. Treatments of young wild-type appressoria with actin and tubulin inhibitors reduce both fungal penetration and ACE1 expression. These experiments suggest that ACE1 appressorium-specific expression does not depend on host plant signals but is connected to the onset of appressorium-mediated penetration.

Application of cDNA array for studying the gene expression profile of mature appressoria of Magnaporthe grisea

Appressorium is an infection structure of the phytopathogenic fungus Magnaporthe grisea. Analysis of gene expression profiles of appressorium development provides insight into the molecular basis of pathogenicity and control of this fungal plant disease. A cDNA array representing 2927 unique genes based on a large EST (expressed sequence tag) database of M. grisea strain Y34 was constructed and used to profile the gene expression patterns at mycelium and appressorium maturation stages. Compared with mycelia, 55 up-regulated and 22 down-regulated genes were identified in mature appressoria. Among 77 genes, 16 genes showed no similarity to the genome sequences of M. grisea. A novel homologue of peptidyl-prolyl cis-trans isomerase was found to be expressed at low-level in mature appressoria of M. grisea. The results indicated that the genes such as pyruvate carboxylase, phospholipid metabolism-related protein and glyceraldehyde 3-phosphate dehydrogenase involved in gluconeogenesis, lipid metabolism and glycolysis, showed differential expression in mature appressoria. Furthermore, genes such as PTH11, beta subunit of G protein and SGT1 involved in cell signalling, were expressed differentially in mature appressoria. Northern blot analysis was used to confirm the cDNA array results.

From genes to genomes: a new paradigm for studying fungal pathogenesis in Magnaporthe oryzae

Magnaporthe oryzae is the most destructive fungal pathogen of rice worldwide and because of its amenability to classical and molecular genetic manipulation, availability of a genome sequence, and other resources it has emerged as a leading model system to study host-pathogen interactions. This chapter reviews recent progress toward elucidation of the molecular basis of infection-related morphogenesis, host penetration, invasive growth, and host-pathogen interactions. Related information on genome analysis and genomic studies of plant infection processes is summarized under specific topics where appropriate. Particular emphasis is placed on the role of MAP kinase and cAMP signal transduction pathways and unique features in the genome such as repetitive sequences and expanded gene families. Emerging developments in functional genome analysis through large-scale insertional mutagenesis and gene expression profiling are detailed. The chapter concludes with new prospects in the area of systems biology, such as protein expression profiling, and highlighting remaining crucial information needed to fully appreciate host-pathogen interactions.

Magnaporthe as a model for understanding host-pathogen interactions

The rice blast pathosystem has been the subject of intense interest in part because of the importance of the disease to world agriculture, but also because both Magnaporthe oryzae and its host are amenable to advanced experimental approaches. The goal of this review is to provide an overview of the system and to point out recent significant studies that update our understanding of the biology of M. oryzae. The genome sequence of M. oryzae has provided insight into how genome structure and pathogen population genetic variability has been shaped by transposable elements. The sequence allows systematic approaches to long-standing areas of investigation, including pathogen development and the molecular basis of compatible and incompatible interactions with its host. Rice blast provides an integrated system to illustrate most of the important concepts governing fungal/plant interactions and serves as an excellent starting point for gaining a broad perspective of issues in plant pathology.

Deep and comparative analysis of the mycelium and appressorium transcriptomes of Magnaporthe grisea using MPSS, RL-SAGE, and oligoarray methods

Background

Rice blast, caused by the fungal pathogen Magnaporthe grisea, is a devastating disease causing tremendous yield loss in rice production. The public availability of the complete genome sequence of M. grisea provides ample opportunities to understand the molecular mechanism of its pathogenesis on rice plants at the transcriptome level. To identify all the expressed genes encoded in the fungal genome, we have analyzed the mycelium and appressorium transcriptomes using massively parallel signature sequencing (MPSS), robust-long serial analysis of gene expression (RL-SAGE) and oligoarray methods.

Results

The MPSS analyses identified 12,531 and 12,927 distinct significant tags from mycelia and appressoria, respectively, while the RL-SAGE analysis identified 16,580 distinct significant tags from the mycelial library. When matching these 12,531 mycelial and 12,927 appressorial significant tags to the annotated CDS, 500 bp upstream and 500 bp downstream of CDS, 6,735 unique genes in mycelia and 7,686 unique genes in appressoria were identified. A total of 7,135 mycelium-specific and 7,531 appressorium-specific significant MPSS tags were identified, which correspond to 2,088 and 1,784 annotated genes, respectively, when matching to the same set of reference sequences. Nearly 85% of the significant MPSS tags from mycelia and appressoria and 65% of the significant tags from the RL-SAGE mycelium library matched to the M. grisea genome. MPSS and RL-SAGE methods supported the expression of more than 9,000 genes, representing over 80% of the predicted genes in M. grisea. About 40% of the MPSS tags and 55% of the RL-SAGE tags represent novel transcripts since they had no matches in the existing M. grisea EST collections. Over 19% of the annotated genes were found to produce both sense and antisense tags in the protein-coding region. The oligoarray analysis identified the expression of 3,793 mycelium-specific and 4,652 appressorium-specific genes. A total of 2,430 mycelial genes and 1,886 appressorial genes were identified by both MPSS and oligoarray.

Conclusion

The comprehensive and deep transcriptome analysis by MPSS and RL-SAGE methods identified many novel sense and antisense transcripts in the M. grisea genome at two important growth stages. The differentially expressed transcripts that were identified, especially those specifically expressed in appressoria, represent a genomic resource useful for gaining a better understanding of the molecular basis of M. grisea pathogenicity. Further analysis of the novel antisense transcripts will provide new insights into the regulation and function of these genes in fungal growth, development and pathogenesis in the host plants.

MGOS: A resource for studying Magnaporthe grisea and Oryza sativa interactions

The MGOS (Magnaporthe grisea Oryza sativa) web-based database contains data from Oryza sativa and Magnaporthe grisea interaction experiments in which M. grisea is the fungal pathogen that causes the rice blast disease. In order to study the interactions, a consortium of fungal and rice geneticists was formed to construct a comprehensive set of experiments that would elucidate information about the gene expression of both rice and M. grisea during the infection cycle. These experiments included constructing and sequencing cDNA and robust long-serial analysis gene expression libraries from both host and pathogen during different stages of infection in both resistant and susceptible interactions, generating >50,000 M. grisea mutants and applying them to susceptible rice strains to test for pathogenicity, and constructing a dual O. sativa-M. grisea microarray. MGOS was developed as a central web-based repository for all the experimental data along with the rice and M. grisea genomic sequence. Community-based annotation is available for the M. grisea genes to aid in the study of the interactions.

Review of fungal chitinases

Chitin is the second most abundant organic and renewable source in nature, after cellulose. Chitinases are chitin-degrading enzymes. Chitinases have important biophysiological functions and immense potential applications. In recent years, researches on fungal chitinases have made fast progress, especially in molecular levels. Therefore, the present review will focus on recent advances of fungal chitinases, containing their nomenclature and assays, purification and characterization, molecular cloning and expression, family and structure, regulation, and function and application.

Generation, annotation and analysis of ESTs from Trichoderma harzianum CECT 2413

BACKGROUND

The filamentous fungus Trichoderma harzianum is used as biological control agent of several plant-pathogenic fungi. In order to study the genome of this fungus, a functional genomics project called "TrichoEST" was developed to give insights into genes involved in biological control activities using an approach based on the generation of expressed sequence tags (ESTs).

RESULTS

Eight different cDNA libraries from T. harzianum strain CECT 2413 were constructed. Different growth conditions involving mainly different nutrient conditions and/or stresses were used. We here present the analysis of the 8,710 ESTs generated. A total of 3,478 unique sequences were identified of which 81.4% had sequence similarity with GenBank entries, using the BLASTX algorithm. Using the Gene Ontology hierarchy, we performed the annotation of 51.1% of the unique sequences and compared its distribution among the gene libraries. Additionally, the InterProScan algorithm was used in order to further characterize the sequences. The identification of the putatively secreted proteins was also carried out. Later, based on the EST abundance, we examined the highly expressed genes and a hydrophobin was identified as the gene expressed at the highest level. We compared our collection of ESTs with the previous collections obtained from Trichoderma species and we also compared our sequence set with different complete eukaryotic genomes from several animals, plants and fungi. Accordingly, the presence of similar sequences in different kingdoms was also studied.

CONCLUSION

This EST collection and its annotation provide a significant resource for basic and applied research on T. harzianum, a fungus with a high biotechnological interest.

Comparative EST analysis provides insights into the basal aquatic fungus Blastocladiella emersonii

Background

Blastocladiella emersonii is an aquatic fungus of the Chytridiomycete class, which is at the base of the fungal phylogenetic tree. In this sense, some ancestral characteristics of fungi and animals or fungi and plants could have been retained in this aquatic fungus and lost in members of late-diverging fungal species. To identify in B. emersonii sequences associated with these ancestral characteristics two approaches were followed: (1) a large-scale comparative analysis between putative unigene sequences (uniseqs) from B. emersonii and three databases constructed ad hoc with fungal proteins, animal proteins and plant unigenes deposited in Genbank, and (2) a pairwise comparison between B. emersonii full-length cDNA sequences and their putative orthologues in the ascomycete Neurospora crassa and the basidiomycete Ustilago maydis.

Results

Comparative analyses of B. emersonii uniseqs with fungi, animal and plant databases through the two approaches mentioned above produced 166 B. emersonii sequences, which were identified as putatively absent from other fungi or not previously described. Through these approaches we found: (1) possible orthologues of genes previously identified as specific to animals and/or plants, and (2) genes conserved in fungi, but with a large difference in divergence rate in B. emersonii. Among these sequences, we observed cDNAs encoding enzymes from coenzyme B₁₂-dependent propionyl-CoA pathway, a metabolic route not previously described in fungi, and validated their expression in Northern blots.

Conclusion

Using two different approaches involving comparative sequence analyses, we could identify sequences from the early-diverging fungus B. emersonii previously considered specific to animals or plants, and highly divergent sequences from the same fungus relative to other fungi.

High level production of the Magnaporthe grisea fructose 1,6-bisphosphate aldolase enzyme in Escherichia coli using a small volume bench-top fermentor

The Class II fructose 1,6-bisphosphate aldolase from the Rice Blast causative agent Magnaporthe grisea was subcloned in the Escherichia coli vector pT7-7. The enzyme was overexpressed using fed-batch fermentation in a small bench-top reactor. A total of 275 g of cells and 1.3 g of highly purified enzyme with a specific activity of 70 U/mg were obtained from a 1.5L culture. The purified enzyme is a homodimer of 39.6 kDa subunits with a zinc ion at the active site. Kinetic characterization indicates that the enzyme has a K(m) of 51 microM, a k(cat) of 46 s(-1), and a pH optimum of 7.8 for fructose 1,6-bisphosphate cleavage. The fermentation system procedure reported exemplifies the potential of using a lab-scale bioreactor for the large scale production of recombinant enzymes.

Identification of an endo-beta-1,4-D-xylanase from Magnaporthe grisea by gene knockout analysis, purification, and heterologous expression

Magnaporthe grisea, a destructive ascomycetous pathogen of rice, secretes cell wall-degrading enzymes into a culture medium containing purified rice cell walls as the sole carbon source. From M. grisea grown under the culture conditions described here, we have identified an expressed sequenced tag, XYL-6, a gene that is also expressed in M. grisea-infected rice leaves 24 h postinoculation with conidia. This gene encodes a protein about 65% similar to endo-beta-1,4-D-glycanases within glycoside hydrolase family GH10. A M. grisea knockout mutant for XYL-6 was created, and it was shown to be as virulent as the parent strain in infecting the rice host. The proteins secreted by the parent strain and by the xyl-6Delta mutant were each fractionated by liquid chromatography, and the collected fractions were assayed for endo-beta-1,4-D-glucanase or endo-beta-1,4-D-xylanase activities. Two protein-containing peaks with endo-beta-1,4-D-xylanase activity secreted by the parent strain are not detectable in the column eluant of the proteins secreted by the mutant. The two endoxylanases (XYL-6alpha and XYL-6beta) from the parent were each purified to homogeneity. N-terminal amino acid sequencing indicated that XYL-6alpha is a fragment of XYL-6beta and that XYL-6beta is identical to the deduced protein sequence encoded by the XYL-6 gene. Finally, XYL-6 was introduced into Pichia pastoris for heterologous expression, which resulted in the purification of a fusion protein, XYL-6H, from the Pichia pastoris culture filtrate. XYL-6H is active in cleaving arabinoxylan. These experiments unequivocally established that the XYL-6 gene encodes a secreted endo-beta-1,4-D-xylanase.

The presence of GC-AG introns in Neurospora crassa and other euascomycetes determined from analyses of complete genomes: implications for automated gene prediction

A combination of experimental and computational approaches was employed to identify introns with noncanonical GC-AG splice sites (GC-AG introns) within euascomycete genomes. Evaluation of 2335 cDNA-confirmed introns from Neurospora crassa revealed 27 such introns (1.2%). A similar frequency (1.0%) of GC-AG introns was identified in Fusarium graminearum, in which 3 of 292 cDNA-confirmed introns contained GC-AG splice sites. Computational analyses of the N. crassa genome using a GC-AG intron consensus sequence identified an additional 20 probable GC-AG introns in this fungus. For 8 of the 47 GC-AG introns identified in N. crassa a GC donor site is also present in a homolog from Magnaporthe grisea, F. graminearum, or Aspergillus nidulans. In most cases, however, homologs in these fungi contain a GT-AG intron or no intron at the corresponding position. These findings have important implications for fungal genome annotation, as the automated annotations of euascomycete genomes incorrectly identified intron boundaries for all of the confirmed and probable GC-AG introns reported here.

Comparative genomic analysis of phytopathogenic fungi using expressed sequence tag (EST) collections

SUMMARY We describe the analysis of 57 727 unique expressed sequence tags (ESTs) from 15 species of phytopathogenic and three species of saprophytic fungi. This resource is held within the COGEME phytopathogen EST database (http://cogeme.ex.ac.uk/). Comparative analysis was performed to investigate the differences between pathogenic and free-living fungi based on a substantial collection of expressed gene sequences and available, completed fungal genome sequences. We report that the expressed gene inventories of pathogenic fungi were not significantly more similar to each other than to those of free-living filamentous fungi. As expected, however, filamentous fungi as a group share more sequences in common than with the free-living yeast species Saccharomyces cerevisiae. Interestingly, ESTs of the obligate biotrophic fungus Blumeria graminis f. sp. hordei were more dissimilar to those of all other fungal species assessed, having a lower number of sequences in common with filamentous ascomycetes studied to date and also possessing a larger proportion of unisequences of unknown function. Our analysis of ESTs in the COGEME database enabled identification of a set of functional groups of genes that are more highly represented in the genomes of pathogenic fungi than non-pathogenic species.