A bioinformatic tool for analysis of EST transcript abundance during infection-related development by Magnaporthe grisea

Untangling the transcriptome from fungus-infected plant tissues

The development of sequencing technology allows low-cost generation of sequence data. The huge amount of raw sequence data now available has introduced many challenges associated with analysis of these large-scale data banks. For example, it is very important to distinguish materials of plant and fungal origin in fungus-infected plant tissue. The origin of transcripts that were sequenced from Library 895-M6 (poplar tissue infected by Marssonina brunnea) on Illumina/Solexa GA IIx was determined by combining three methods: (1) based on the taxonomic information of homologous sequences; (2) based on the reference genome sequence; (3) based on the transcriptome sequence of the host and its pathogen obtained from Library 895 (poplar) and Library M6 (M. brunnea) as well as Library 895-M6 (mixture of poplar and M. brunnea). We idenified accurately the origin of 80,978 (99.5%) contigs in the mixed poplar and M. brunnea sample (Library 895-M6) by integrating the results from the three methods. The results of this study demonstrate that a combination of these three approaches described here is an effective strategy for determining the origin of sequences in a mixed pool, and provides a basis for further transcriptome analysis of the mixed sample.

Peroxisomal alanine: glyoxylate aminotransferase AGT1 is indispensable for appressorium function of the rice blast pathogen, Magnaporthe oryzae

The role of β-oxidation and the glyoxylate cycle in fungal pathogenesis is well documented. However, an ambiguity still remains over their interaction in peroxisomes to facilitate fungal pathogenicity and virulence. In this report, we characterize a gene encoding an alanine, glyoxylate aminotransferase 1 (AGT1) in Magnaporthe oryzae, the causative agent of rice blast disease, and demonstrate that AGT1 is required for pathogenicity of M. oryzae. Targeted deletion of AGT1 resulted in the failure of penetration via appressoria; therefore, mutants lacking the gene were unable to induce blast symptoms on the hosts rice and barley. This penetration failure may be associated with a disruption in lipid mobilization during conidial germination as turgor generation in the appressorium requires mobilization of lipid reserves from the conidium. Analysis of enhanced green fluorescent protein expression using the transcriptional and translational fusion with the AGT1 promoter and open reading frame, respectively, revealed that AGT1 expressed constitutively in all in vitro grown cell types and during in planta colonization, and localized in peroxisomes. Peroxisomal localization was further confirmed by colocalization with red fluorescent protein fused with the peroxisomal targeting signal 1. Surprisingly, conidia produced by the Δagt1 mutant were unable to form appressoria on artificial inductive surfaces, even after prolonged incubation. When supplemented with nicotinamide adenine dinucleotide (NAD⁺)+pyruvate, appressorium formation was restored on an artificial inductive surface. Taken together, our data indicate that AGT1-dependent pyruvate formation by transferring an amino group of alanine to glyoxylate, an intermediate of the glyoxylate cycle is required for lipid mobilization and utilization. This pyruvate can be converted to non-fermentable carbon sources, which may require reoxidation of NADH generated by the β-oxidation of fatty acids to NAD⁺ in peroxisomes. Therefore, it may provide a means to maintain redox homeostasis in appressoria.

A novel protein Com1 is required for normal conidium morphology and full virulence in Magnaporthe oryzae

In Magnaporthe oryzae, pyriform conidia are the primary inoculum and the main source for disease dissemination in the field. In this study, we identified and characterized the COM1 gene that was disrupted in three insertional mutants producing slender conidia. COM1 encodes a putative transcription regulator unique to filamentous ascomycetes. The com1 disruption and deletion mutants had similar defects in conidium morphology and were significantly reduced in virulence on rice and barley seedlings. Microscopic examination revealed that the Deltacom1 mutants were defective in appressorium turgor generation, penetration, and infectious growth. COM1 was expressed constitutively in M. oryzae. The Com1 protein had putative helix-loop-helix structures and three predicted nuclear localization signal sequences. In transformants expressing COM1(335-613)-enhanced green fluorescent protein fusion constructs, fluorescence signals were observed in the nucleus. Our data indicated that the COM1 gene may encode a novel transcription regulator that regulates conidial development and invasive growth in M. oryzae.

Computational techniques for elucidating plant-pathogen interactions from large-scale experiments on fungi and oomycetes

Eukaryotic plant pathogens are responsible for the destruction of billions of dollars worth of crops each year. With large-scale genomics of both pathogens and hosts and the corresponding computational analysis, biologists are now able to gain knowledge about many pathogenic and defense genes concurrently. To study the interactions between these two organism groups, it is necessary to design experiments to elucidate the genes being expressed during the invasion of the pathogen into the host. For the most part, this does not require new software development, though it does require the use of existing software in novel ways. This article provides a broad overview of several key and illustrative experiments and the corresponding computational analyses, outlining the knowledge gained in each. It goes on to describe databases for plant-pathogen data and important initiatives such as Plant-Associated Microbe Gene Ontology. It discusses how various emerging approaches will increase the power of computers in host-pathogen interaction studies.

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.

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/.

Functional characterization of CgCTR2, a putative vacuole copper transporter that is involved in germination and pathogenicity in Colletotrichum gloeosporioides

Copper is a cofactor and transition metal involved in redox reactions that are essential in all eukaryotes. Here, we report that a vacuolar copper transporter that is highly expressed in resting spores is involved in germination and pathogenicity in the plant pathogen Colletotrichum gloeosporioides. A screen of C. gloeosporioides transformants obtained by means of a promoterless green fluorescent protein (GFP) construct led to the identification of transformant N159 in which GFP signal was observed in spores. The transforming vector was inserted 70 bp upstream of a putative gene with homology to the Saccharomyces cerevisiae vacuolar copper transporter gene CTR2. The C. gloeosporioides CTR2 (CgCTR2) gene fully complemented growth defects of yeast ctr2Delta mutants, and a CgCTR2-cyan fluorescent protein (CFP) fusion protein accumulated in vacuole membranes, confirming the function of the protein as a vacuolar copper transporter. Expression analysis indicated that CgCTR2 transcript is abundant in resting conidia and during germination in rich medium and downregulated during "pathogenic" germination and the early stages of plant infection. CgCTR2 overexpression and silencing mutants were generated and characterized. The Cgctr2 mutants had markedly reduced Cu superoxide dismutase (SOD) activity, suggesting that CgCTR2 is important in providing copper to copper-dependent cytosolic activities. The Cgctr2-silenced mutants had increased sensitivity to H2O2 and reduced germination rates. The mutants were also less virulent to plants, but they did not display any defects in appressorium formation and penetration efficiency. An external copper supply compensated for the hypersensitivity to H2O2 but not for the germination and pathogenicity defects of the mutants. Similarly, overexpression of CgCTR2 enhanced resistance to H2O2 but had no effect on germination or pathogenicity. Our results show that copper is necessary for optimal germination and pathogenicity and that CgCTR2 is involved in regulating cellular copper balance during these processes.

Cloning, sequencing and expression analysis of the NAR promoter activated during hyphal stage of Magnaporthe grisea

The promoter of NAR gene in Magnaporthe grisea was isolated and sequenced. The promoter sequences contained the "TATA" box, the "CAAT" box, and binding sites for fungal regulatory proteins. Programs that predict promoter sequences indicated that promoter sequence lies between locations 430 and 857 of the NAR promoter fragment. GFP expression under the NAR promoter and NAR transcript analysis revealed that this promoter is activated primarily at the mycelial stage in the rice blast fungus and could be used to express native or extrinsic genes in the mycelia of the rice blast fungus.

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.

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.

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.

The dawn of fungal pathogen genomics

Recent advances in sequencing technologies have led to a remarkable increase in the number of sequenced fungal genomes. Several important plant pathogenic fungi are among those that have been sequenced or are being sequenced. Additional fungal pathogens are likely to be sequenced in the near future. Analysis of the available genomes has provided useful information about genes that may be important for plant infection and colonization. Genome features, such as repetitive sequences, telomeres, conserved syntenic blocks, and expansion of pathogenicity-related genes, are discussed in detail with Magnaporthe oryzae (M. grisea) and Fusarium graminearum as examples. Functional and comparative genomic studies in plant pathogenic fungi, although still in the early stages and limited to a few pathogens, have enormous potential to improve our understanding of the molecular mechanisms involved in host-pathogen interactions. Development of advanced genomics tools and infrastructure is critical for efficient utilization of the vast wealth of available genome sequence information and will form a solid foundation for systems biology studies of plant pathogenic fungi.