Identification and analysis of in planta expressed genes of Magnaporthe oryzae
© Kim et al; licensee BioMed Central Ltd. 2010
Received: 4 September 2009
Accepted: 10 February 2010
Published: 10 February 2010
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.
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.
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.
Rice blast, caused by Magnaporthe oryzae, is one of the most devastating diseases in rice growing regions worldwide, causing 11-15% yield loss annually . Genetic tractability as well as economic importance makes this disease a model pathosystem to understand plant-microbe interactions. The genome sequences of both organisms are available [2–4], and both forward and reverse genomic studies to understand molecular mechanisms for pathogenesis on a genome scale have been undertaken . Understanding the precise molecular mechanisms of infection will facilitate design of novel control strategies.
The process of M. oryzae infection starts when a conidium lands on the rice leaf surface. After germination by hydration, an appressorium develops at the tip of the germ tube, from which a penetration peg emerges to penetrate the cuticle layer into the rice cell using mechanical force. Within the plant cell, the fungus faces two different fates. In an incompatible interaction, resistance gene products recognize corresponding avirulence gene products from the invading pathogen and invoke a series of defense responses to restrict pathogen growth . In a compatible interaction, however, the host plant mobilizes defense responses much later, resulting in visible coalescing lesions. From these expanding lesions, a conidiophore emerges, releasing tens of thousands of new conidia able to start a second round of infection. In this respect, the gene products of the pathogen expressed during infectious growth may play roles as pathogenicity factors required for evasion of the host's defense responses and successful colonization. This theory is supported by our previous pilot study in which analysis of in planta expressed genes identified by EST sequencing of an M.oryzae infected rice cDNA library . By sequencing 511 randomly selected cDNA clones, 72 of 293 ESTs could be assigned as fungal genes with sequence similarity to NCBI entries. Among them, MHP1, encoding a Class II hydrophobin, was preferentially expressed during infectious growth and was involved in full virulence of rice blast fungus, acting in the late stages of infection . We also identified novel pathogenesis-related protein genes that are upregulated in rice after pathogen infection [9, 10]. Jantasuriyarat et al.  performed large-scale analysis of ESTs during M.oryzae/rice interactions at the early stage of infection. Plant defense genes were well represented in this study. However, fungal genes involved in interactions with the host were scarce, primarily because they harvested infected tissue at 6 and 24 h after inoculation, at which time fungal spores had just started to penetrate. Genome-wide transcriptome analyses have also been conducted on this fungus within the last few years through the use of EST, SAGE, and RL-SAGE [7, 11–16]. Most of these studies focused on pre-penetration stages, such as conidiation and appressorium formation, and on growth in a variety of in vitro conditions, including complete medium, minimal medium, nitrogen starvation, and rice cell wall medium. From these global efforts, 28,682 EST sequences comprising 8,821 uniESTs of M. oryzae are registered in the dbEST of NCBI, which covers about 80% of the genes electronically annotated for this fungus . However, our knowledge of the molecular basis of the establishment of disease and proliferation within the plant cells remains limited.
During recent years, studies aimed at cataloging pathogen genes expressed during interactions with their host have been conducted in many pathosystems. Large-scale EST analysis with cDNA libraries from infected plant tissues has been used for this purpose . More than 1,000 fungal genes could be identified out of 1,869 EST sequences from senescent leaves (21-25 days after inoculation) with lesions containing visible pycnidia. However, few fungal genes could be recovered from Fusarium graminearum infected wheat leaves (86 of 3,546 uniESTs) , from Brassica napus stem segments with expanding lesions caused by Sclerotinia sclerotiorum (52 of 767 uniESTs) , from wheat roots infected by take-all fungus (9 of 114 ESTs) , or M.oryzae infected rice leaves (4 of 13,570 uniESTs) . Low-level representation of fungal genes may be because they harvested infected tissues at early time points. Suppression subtractive hybridization was also applied to identify genes expressed during the pathogen's interaction with its host [20, 21]. Broeker et al.  identified 81 Puccinia graminis f. sp. tritici genes by sequencing 454 random clones from rust-infected wheat plants subtracted from those of healthy leaves. In this study, we used two complementary approaches: EST analysis and suppression subtractive hybridization and subsequent sequencing to catalog genes expressed during rice-M.oryzae interactions at the genome level. 712 fungal genes were identified from 2,315 uniESTs and their putative functions were assessed with GO. We were able to identify 100 novel fungal genes specific to infection. cDNAs with full-length ORFs were used to correct ab initio gene models annotated from genome sequences, and to analyze fungal splice site context.
Generation of infection ESTs
Summary of ESTs generated from infection and subtractrion libraries
Number of sequences analyzed
Number of contigs
Number of singletons
Number of cDNAs contained in contigs
The most frequently represented sequence was contig203, which showed sequence similarity to elongation factor 1-α of the dimorphic fungus Ajellomyces capsulatus . It occurred 76 times from SL and five times from IL. Contig226 was sequenced 63 times only from SL, and Contig200, showing sequence similarity to senescence-associated protein of Pisum sativum , appeared 29 times in SL and once in IL. The repertoire of highly redundant clones differed between the two libraries (Additional file 1, Table S1). In IL, genes involved in photosynthesis, such as RuBisCo activase (Contig 526) and RuBisCo small subunit C (Contig 549), and plant defense-related genes, such as metallothionein (Contig 496) and pathogenesis-related protein 1 (Contig 597), were included in the 10 most redundant contigs. Only two contigs (Contig 310 and Contig 370) of the 10 most redundant contigs were considered to be of fungal origin. However, fungal genes made up 50% of the most redundant 10 contigs of SL. Additionally, after subtracting, genes involved in photosynthesis, defense-related genes including senescence-associated protein (Contig 200), probenazole-induced protein (Contig 173), and metallothionein I (Contig 153) were the most abundant sequences.
Annotation of the origin of genes
Classification of 2,315 uniESTs according to origin of the sequences
Contigs (No. ESTs)
Infection library specific
Subtraction library specific
GO-based functional classification
Proteins assigned as extracellular were poorly represented in this analysis, although secretion proteins are thought to play important roles in host interactions [26, 27] and virulence [8, 28, 29]. Thus, we screened 2,302 uniESTs for signal peptides using SignalP3.0 . Because most EST sequences did not bear a start codon, genomic resources were used. For fungal genes, ab initio annotated protein sequences matched to infection ESTs were subjected to SignalP analysis. In the case of rice genes, where 32,127 full-length (fl-) cDNA sequences are available, fl cDNA sequences matched to infection ESTs were translated and in-frame amino acid sequences were screened with SignalP. As a result, 63 fungal genes and 246 rice genes of 2,302 uniESTs were believed to have a signal peptide. Among these were three hydrophobins (F204, F233 encoding MHP1, and F238 encoding MPG1), two of which have been identified as secreted and associated with the cell wall [8, 29]. ESTs encoding spore coat proteins (F207) and cell wall degrading enzymes, such as chitinase3 (F037) and 1,4-beta-D-glucan cellobiohydrolase (F298), were also identified.
Search for novel genes uniquely or preferentially expressed during rice-M. oryzae interactions
Expression profiles of infection-enriched genes
Use of ESTs in validation and improvement of genome annotation
Splicing site context
Fungal genes expressed during infectious growth within their host and plant genes expressed after attack by their pathogen may have important role(s) in virulence and defense responses. In this study, we identified 2,302 uniESTs expressed during rice-M. oryzae interactions using two complementary approaches. The first approach was EST analysis (IL) of blast fungus-infected rice leaves with rapidly expanding typical blast lesions. We harvested leaves with typical lesions and collected these daily from 410 days after inoculation. 1,539 uniESTs were identified from the EST analysis. The second strategy was using SSH (SL) to enrich fungal genes and rice genes, whose expression was induced during interactions, resulting in the generation of 963 uniESTs.
Fungal genes expressed during infectious growth were identified
cDNA libraries were constructed with RNA from rice leaves infected with rice blast fungus. Thus, pathogen genes, as well as host plant genes expressed during their interaction, were contained in these libraries. Assigning the origin of ESTs can be complex. The G+C content was successfully used to identify oomycete genes from soybean genes infected with Phytophthora sojae ; however, this approach is not applicable for pathosystems involving ascomycete fungi. Sequence similarity to previously identified genes is widely used as a major criterion to distinguish their origin [7, 11, 17, 19, 20]. However, ESTs that did not show significant sequence similarity to database entries could not be annotated. This ranged from 30% to more than 40% of ESTs, depending on the organism. We could not identify the origin of 42.7% of ESTs in our previous study . Li et al.  adapted codon usage pattern as a subsidiary tool to differentiate Sclerotinia sclerotiorum genes from Brassica napus genes. However, rice and M. oryzae genes do not show significant differences in their codon usage pattern when we analyzed full-length cDNA sequences (data not shown). In this study, we took advantage of genome sequences and many EST sequences of both organisms to distinguish the origin of the ESTs. With the BLASTX results against NCBI nr, BLASTN data against genomic and EST sequences of both organisms, 99% of uniESTs could be annotated for their origin (Additional file 1 Table S1). Among these, 712 uniESTs were identified as fungal genes, comprising 31% of the 2,302 uniESTs. 402 uniESTs (533 ESTs) were from the infected library, comprising 24% of the uniESTs (27% of total ESTs), and 310 uniESTs were added from the subtraction library.
A relatively high proportion of ESTs were of fungal origin. This result is consistent with our previous study , where about 24.6% of uniESTs were believed to be fungal genes. Interaction transcriptome  approaches have been used for understanding molecular mechanisms of defense and/or pathogenicity in a variety of pathosystems [7, 11, 17–20, 37]. However, only a limited number (or proportion) of fungal genes were identified. This is primarily because previous authors used plant material at early stages of infection, or partially resistant host cultivars. Jantasuriyarat et al.  analyzed 68,920 ESTs from eight cDNA libraries including six constructed from rice leaves infected with rice blast fungus, and found only four fungal genes. They collected rice leaves at 6 and 24 h after blast fungus inoculation, which are the very early stages of infection when fungal conidia have just germinated or just started to penetrate into rice epidermal cells. Talbot et al.  estimated fungal biomass to be 10% at 72 h after inoculation, at which time lesions start to appear. Recovering a high portion of fungal genes in this study might be attributed to the leaf material from which total RNA was isolated. We collected diseased leaves at late stages of infection daily up to 10 days after inoculation, when the whole leaf blade was covered with expanding lesions and had started to wither. Thus, fungal genes found comprised 27% of uniESTs in the infection library; this result was not unexpected. These fungal genes will provide a good resource to understand fungal pathogenicity, especially at later stages of infection.
Subtraction efficiency is high
43.5% of fungal genes (310 uniESTs of 712 fungal genes) were obtained from a subtraction library. 50% (1,128 of 2,259) of the ESTs from the subtraction library were identified as fungal genes. We used uninfected rice RNA as a driver to enrich fungal genes as well as rice genes that were up-regulated on fungal infection. Gilleroux and Osbourn  used two driver cDNA populations (i.e., one derived from mock-inoculated roots and the other from a G. graminis culture grown in complete medium) to subtract cDNAs from G.graminis-infected wheat roots. However, they did not identify any fungal gene using this strategy, possibly because fungal genes comprise a low proportion of infected root tissue and/or the fungal genes expressed during plant infection are also expressed during growth in vitro, to be subtracted out by the fungal driver. Using uninfected plant RNA as a driver turned out to be a useful strategy to identify candidate fungal genes to be functionally characterized. In addition to cataloging fungal genes expressed during infection, we could identify suites of rice genes encoding ABC transporter families specifically obtained from the subtraction library. Forty-two ESTs (18 uniESTs) showed sequence similarity to 11 different ABC transporter families, of which 38 ESTs were from the subtraction library. Seventeen ESTs encoded PDR9, expression of which is induced by environmental stresses, such as heavy metals, hypoxic stress, and redox perturbations . Direct evidence for a role of ABC transporters in plant defense was obtained from the study of NpPDR1 (formerly named NpABC1), a plasma membrane PDR type ABC transporter of Nicotiana plumbaginifolia, which was induced by the elicitor analog sclareolide, sclareol, a virulent strain of Pseudomonas syringe pv. tabaci and by non-pathogenic pseudomonads, P. fluorescens and P.marginalis pv marginalis , and is believed to be involved in the secretion of the antimicrobial diterpene, sclareol . Transgenic plants in which NpPDR1 expression was prevented by RNA interference showed reduced resistance to the fungal pathogen Botrytis cinerea . Thus, ABC transporters identified in this study may play a role in rice defense against the invading fungus. These results lead us to suggest that various materials toxic to plant cells might be excreted from invading fungal cells into the plant cytosol, and that plant cells struggle to pump out and/or detoxify these deleterious materials. This was further supported by the fact that a high proportion (10.4%, 74 genes) of in planta expressed fungal genes encoded proteins with a signal peptide.
There was little overlap in the fungal genes identified from the two libraries in terms of uniEST. This may be attributed to the low depth of coverage in this study. Another explanation may be drawn if we consider the number of total ESTs, as most overlapped sequences are highly redundant so that 587 ESTs comprised 75 uniESTs.
Novel genes identified
Although EST collections form the foundation for various genome-scale experiments within as yet unsequenced genomes , they provide valuable information for gene discovery, genome annotation, and expression profiles, even in organisms where the genome sequence is available. They have also been successfully applied to identify genes expressed during the interaction of the plant pathogen with its host [7, 35]. Although 12,465 uniESTs of M. oryzae from nine different libraries  were already available, we found 100 previously unidentified unigenes through analysis with strict criteria. This exactly identified 14% of the 712 fungal genes. Sixty seven unigenes were sequenced from the infection library, while 37 genes were added from the subtraction library. The nine libraries that Ebbole et al.  used were constructed with RNA covering the fungal life cycle, including conidia, appressoria forming germlings, and mixed culture undergoing mating, which are stress conditions thought to mimic the environment the fungus might encounter in nature, such as rice cell walls, complete medium, minimal medium, nitrogen starvation, and mixed stress conditions, and mycelia from mutants that could not elaborate the infection structure, the appressorium. Nutrient starvation, especially nitrogen-deprived conditions, is regarded as an environmental cue for the fungus to infect the plant [41, 42]. M. oryzae secrets proteins that cause senescence of rice leaves, reminiscent of the symptoms caused by the fungus itself when it is under nitrogen starvation . NPR1 and NPR2, nitrogen-regulatory genes non-allelic to NUT1, were reported to positively regulate both nitrogen metabolism and pathogenicity . However, many genes were newly identified only from in planta ESTs, indicating that the repertories of genes expressed during infectious growth and under nutrient stress might differ, and that nutrient stress might not fully represent the conditions that the fungus would meet within the plant cell, as suggested in our previous work .
Resources for improvement of genome annotation
ESTs can be used in gene discovery and genome annotation. Of 712 ESTs from this study, 125 (17.5% of fungal ESTs) did not match any electronically annotated gene. These sequences will help in further correction of genome annotations. Validation of genome annotation could further be confirmed by comparing full-length cDNA sequences with genome sequences. Full-length or near full-length cDNA encompassing the complete ORF can serve as a valuable resource for accurate genome annotation and functional analysis. Large-scale full-length cDNA sequencing projects with the aim of identifying the complete transcriptome have been conducted and used for improving genome annotation in many model organisms, including Arabidopsis thaliana [44–47] and rice , as well as humans  and mice [49, 50]. For example, 32% of the first version of the Arabidopsis gene model was found to be inaccurately annotated, of 10,507 genes having fl-cDNA sequences, and was improved by incorporation of fl-cDNA sequences . The rice genome community has recognized the value of fl-cDNA sequences and deep EST resources, and acquired more than 28,000 fl-cDNA clones prior to completion of the genome sequence . Few efforts have been made to collect and use full-length cDNA sequences in filamentous fungi, in which targeted gene deletion mutants could readily be obtained by homologous recombination due to the haploid nuclear status. However, full-length cDNA sequences from filamentous fungi still have the same significance for accurate genome annotation, functional analysis, and for proteomics as in other eukaryotes, because fungal genomes are also interspersed with non-coding DNA sequences. In this respect, full-length cDNA clones acquired in this study could be a valuable resource for the community. Of 83 full-length cDNAs, 24 (28.6%) showed different gene structures from the ab initio annotated gene model. Electronic annotation started at different positions in 11 cases and erroneously stopped in eight cases, resulting in longer or shorter ORFs. In the remaining cases, electronic annotation started and stopped at the exact sites, but depicted erroneous internal exonintron boundaries. In most cases, erroneous splicing was attributed to the wrong gene model. The intron cis elements (i.e., 5'-GU...AG-3' donor-acceptor splicing site pairs, 5'- and 3'-ss consensus sequence, internal branch site, and polypyrimidine tracts) were also conserved in the M. oryzae genome, and were consistent with results for other filamentous fungi, with more degenerate 5'-splice sites, as suggested by Kupfer et al. . These results can be used to further improve genome annotation.
In summary, genome-wide identification of genes expressed during the interaction between rice and the invading fungal pathogen M. oryzae was carried out and resulted in cataloging of 2,302 uniESTs. More than 700 fungal genes were identified, of which 100 genes were newly identified in this study, and more than 80 genes with full-length ORF sequences were used to validate the genome annotation. Further characterization of the genes reported here will help to unravel the mechanisms of pathogenicity, especially at later stages of infection, as well as the defense responses of the host plant.
Infection of plants by pathogens and the subsequent disease development are accompanied by substantial changes in the physiology of both partners. Analysis of genes expressed during these interactions represents a powerful strategy to obtain insights into the molecular events underlying these changes. Here, we used two complementary strategies, EST and SSH, to catalog genes during a compatible rice and blast fungus interaction at late stages of infection. Fungal genes constituting as many as 30% of uniESTs, as well as plant defense genes, were identified. We identified 100 previously unreported fungal genes despite the wealth of genome resources on this fungus. Additionally, whole sequences from clones having full-length cDNAs were determined by primer walking, which was used to verify the ab initio annotated gene model of M. oryzae genome sequences and splicing context analysis. Taken together, the data obtained in this study will serve as a valuable resource to expand our understanding of genome organization and molecular aspects during rice M. oryzae interactions.
Construction of cDNA libraries
Rice plants (Oryza sativa L. cv. Hwacheong) at the 3rd to 4th leaf stages were spray-inoculated with 2 × 105 spores/ml. After keeping in a dew chamber for 24 h at 25°C, the inoculated plants were transferred to the greenhouse to let the lesions develop. Rice leaves on which typical susceptible type lesions spread through the entire leaves were cut and immediately frozen in liquid nitrogen daily up to 10 days after inoculation (in planta sample). Total RNA was isolated from the frozen plant tissues through phenol/chloroform extraction, followed by lithium chloride precipitation . Poly A(+) RNA was purified using the PolyATrack mRNA isolation system (Promega) and used for construction of two libraries. cDNA was synthesized using the lamda ZAP cDNA synthesis kit and cDNAs larger than 0.5 kb after gel filtration were used to construct the cDNA library (infection library: IL) using the ZAP Express cDNA Gigapack III Gold Cloning kit (Stratagene, La Jolla, CA, USA). Mass in vivo excision was conducted to rescue phagemid clones. Individual colonies were grown in LB medium amended with 100 μg/ml ampicillin in 96 well plates and used for plasmid preparation. A subtracted cDNA library was also constructed by suppression subtractive hybridization using PCR-select™ cDNA subtraction kit (Clontech, Palo Alto, CA, USA) following the manufacturer's instructions. cDNA obtained from M. oryzae infected rice leaves was used as a 'tester' and cDNA from uninfected healthy leaves was used as a 'driver'. Subtracted PCR products were cloned into pCR4 blunt-TOPO vector (Invitrogen) and transformed into E. coli strain TOP10 (Invitrogen) with blue-white selection. White colonies were grown in LB plates amended with ampicillin and used for plasmid preparation.
EST sequencing and analysis
Plasmid DNAs were purified from 1 ml each of the overnight grown bacterial cultures and subjected to automated sequencing using ABI PRISM BigDye Terminator on ABI Prism 3700 sequencer (Applied Biosystems, Foster City, CA, USA) or using ET terminator dye (Amersham, Uppsala, Sweden) on the multicapillary sequencer of RISA384 system (Shimadzu, Tokyo, Japan). Sequence analysis was conducted with an automated pipeline embedded in the web-based bioinformatic systems, Comparative Fungal Genomics Platform (http://cfgp.snu.ac.kr)  and Fungal EST Database (FED; http://fedb.snu.ac.kr, Park et al., in preparation). Sequence data in the chromatogram file were processed using Phred  to call bases with cutoff value >20 and Crossmatch to mask vector sequences. The output sequences were further trimmed using a computer script to get rid of the poly A or poly T tracks. The sequences longer than 100 bp were collected and clustered using CAP3 . The resulting uniESTs were searched against GenBank nr database using a BLASTX algorithm, and against rice [3, 4] and M. oryzae  genome sequences with the BLASTN algorithm. Functional categorization was conducted with GO term assignment. Ab initio annotated ORFs from M. oryzae genome and the translated rice full-length cDNA sequences, which have corresponding infectionESTs were used to search the InterPro database (version12), and the InterPro terms were mapped to GO terms as described in Ebbole et al. . For full-length cDNA sequences, representing cDNA clones were retrieved and their sequences read by primer walking. The intron and exon databases were created by using FELINES  with the parameters used by Kupfer et al. . Splicing site consensus sequences for 5' exon-intron junctions, branch point, and 3' intron-exon junctions were calculated using WebLogo server at (http://weblogo.threeplusone.com/). EST sequences were deposited into dbEST of NCBI under the serial accession number from GT930449 to GT934681, GT968229, and GT968230. The full length cDNA sequences were deposited into NCBI under the serial accession numbers from GU395207 to GU395289 and listed in Table S3.
Real-time RT-PCR was conducted following the protocols described previously  with slight modification. Total RNA for cDNA library was used to assess the expression level during infectious growth. Fungal tissue for expression dynamics during development and nutrient starvation was prepared according to the procedures as described previously . Conidia bearing mycelia grown on oatmeal agar for 10 days at 25°C were harvested by scraping the surface with razor blade. A total of 2 × 107 conidia were incubated in liquid complete media for 8 hours at 25°C by shaking to get germinated conidia. Mycelia for nutrient starvation conditions were prepared as described in Talbot et al. . Samples were frozen in liquid nitrogen immediately after harvest, ground with pre-chilled mortar and pestle, and stored at -80°C before use. Total RNA was isolated from frozen mycelial powder using an Easy-Spin RNA extraction kit (iNtRON Biotechnology, Seoul, Korea). Five micrograms of total RNA was reverse transcribed into first-strand cDNA with oligo (dT) primer using SuperScript™ First-Strand Synthesis System (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instruction, and diluted to the concentration of input RNA to be 12.5 ng/μl with nuclease-free water. Reactions were done in a 10 μl volume containing 100 nM of each primer, 2 μl cDNA (25 ng of input RNA) and 5 μl of 2× iQ™ SYBR® Green Supermix (Bio-rad, Hercules, CA, USA). Real-time PCR was run on an iCycler iQ Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). After 3 min denaturation at 95°C, samples were run for 40 cycles of 15 s at 95°C, 30 s at 60°C, and 30 s at 72°C. After each run, amplification specificity was checked with a dissociation curve acquired by heating the samples from 60 to 95°C. Relative abundance of transcripts in complete media was calculated with 2-ΔCt, where ΔCt = (Ct, gene of interest- Ct, cyclophilin)CM. Fold changes during infectious growth and growth under nutrient starvation compared to growing in liquid complete medium were calculated with 2-ΔΔCt, where ΔΔCt = (Ct, gene of interest- Ct, cyclophilin)test condition - (Ct, gene of interest - Ct, cyclophilin)control . Real-time PCR was conducted with 3 replicates. The primer pairs for transcripts amplification were denoted in Additional file 4 (Table S4).
This work was supported by the National Research Foundation of Korea grants (2009-0063340 and 2009-0080161) and the grants from the Biogreen21 project (20080401-034-044-009-01-00), the TDPAF (309015-04-SB020), and the Crop Functional Genomics Center (2009K001198). J. Park is grateful for the graduate fellowship through the Brain Korea 21 Program.
- Baker B, Zambryski P, Staskawicz B, Dinesh-Kumar SP: Signaling in plant-microbe interactions. Science. 1997, 276: 726-733. 10.1126/science.276.5313.726.PubMedView Article
- Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu JR, Pan H: The genome sequence of the rice blast fungus Magnaporthe grisea. Nature. 2005, 434 (7036): 980-986. 10.1038/nature03449.PubMedView Article
- Goff SA, Ricke D, Lan TH, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H: A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science. 2002, 296 (5565): 92-100. 10.1126/science.1068275.PubMedView Article
- Yu J, Hu S, Wang J, Wong GK, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X: A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002, 296 (5565): 79-92. 10.1126/science.1068037.PubMedView Article
- Jeon J, Park SY, Chi MH, Choi J, Park J, Rho HS, Kim S, Goh J, Yoo S, Park JY: Genome-wide functional analysis of pathogenicity genes in the rice blast fungus. Nat Genet. 2007, 39 (4): 561-565. 10.1038/ng2002.PubMedView Article
- Ahn IP, Kim S, Kang S, Suh SC, Lee YH: Rice Defense Mechanisms Against Cochliobolus miyabeanus and Magnaporthe grisea Are Distinct. Phytopathology. 2005, 95 (11): 1248-1255. 10.1094/PHYTO-95-1248.PubMedView Article
- Kim S, Ahn IP, Lee YH: Analysis of genes expressed during rice-Magnaporthe grisea interactions. Mol Plant-Microbe Interact. 2001, 14 (11): 1340-1346. 10.1094/MPMI.2001.14.11.1340.PubMedView Article
- Kim S, Ahn IP, Rho HS, Lee YH: MHP1, a Magnaporthe grisea hydrophobin gene, is required for fungal development and plant colonization. Mol Microbiol. 2005, 57 (5): 1224-1237. 10.1111/j.1365-2958.2005.04750.x.PubMedView Article
- Kim S, Ahn IP, Park CH, Park SG, Park SY, Jwa NS, Lee YH: Molecular characterization of the cDNA encoding an acidic isoform of PR-1 protein in rice. Mol Cells. 2001, 11 (1): 115-121.PubMed
- Park CH, Kim S, Park JY, Ahn IP, Jwa NS, Im KH, Lee YH: Molecular characterization of a pathogenesis-related protein 8 gene encoding a class III chitinase in rice. Mol Cells. 2004, 17 (1): 144-150.PubMed
- Jantasuriyarat C, Gowda M, Haller K, Hatfield J, Lu G, Stahlberg E, Zhou B, Li H, Kim H, Yu Y: Large-scale identification of expressed sequence tags involved in rice and rice blast fungus interaction. Plant Physiol. 2005, 138 (1): 105-115. 10.1104/pp.104.055624.PubMed CentralPubMedView Article
- Ebbole DJ, Jin Y, Thon M, Pan H, Bhattarai E, Thomas T, Dean R: Gene discovery and gene expression in the rice blast fungus, Magnaporthe grisea: analysis of expressed sequence tags. Mol Plant-Microbe Interact. 2004, 17 (12): 1337-1347. 10.1094/MPMI.2004.17.12.1337.PubMedView Article
- Gowda M, Jantasuriyarat C, Dean RA, Wang GL: Robust-LongSAGE (RL-SAGE): a substantially improved LongSAGE method for gene discovery and transcriptome analysis. Plant Physiol. 2004, 134 (3): 890-897. 10.1104/pp.103.034496.PubMed CentralPubMedView Article
- Irie T, Matsumura H, Terauchi R, Saitoh H: Serial Analysis of Gene Expression (SAGE) of Magnaporthe grisea: genes involved in appressorium formation. Mol Genet Genomics. 2003, 270 (2): 181-9. 10.1007/s00438-003-0911-6.PubMedView Article
- Lu JP, Liu TB, Lin FC: Identification of mature appressorium-enriched transcripts in Magnaporthe grisea, the rice blast fungus, using suppression subtractive hybridization. FEMS Microbiol Lett. 2005, 245 (1): 131-137. 10.1016/j.femsle.2005.02.032.PubMedView Article
- Rauyaree P, Choi WB, Fang EG, Blackmon B, Dean RA: Genes expressed during early stages of rice infection with the rice blast fungus Magnaporthe grisea. Mol Plant Pathol. 2001, 2 (6): 347-354. 10.1046/j.1464-6722.2001.00085.x.PubMedView Article
- Keon J, Antoniw J, Rudd J, Skinner W, Hargreaves J, Hammond-Kosack K: Analysis of expressed sequence tags from the wheat leaf blotch pathogen Mycosphaerella graminicola (anamorph Septoria tritici). Fungal Genet Biol. 2005, 42 (5): 376-389. 10.1016/j.fgb.2004.12.005.PubMedView Article
- Kruger WM, Pritsch C, Chao S, Muehlbauer GJ: Functional and comparative bioinformatic analysis of expressed genes from wheat spikes infected with Fusarium graminearum. Mol Plant-Microbe Interact. 2002, 15 (5): 445-455. 10.1094/MPMI.2002.15.5.445.PubMedView Article
- Li R, Rimmer R, Buchwaldt L, Sharpe AG, Seguin-Swartz G, Coutu C, Hegedus DD: Interaction of Sclerotinia sclerotiorum with a resistant Brassica napus cultivar: expressed sequence tag analysis identifies genes associated with fungal pathogenesis. Fungal Genet Biol. 2004, 41 (8): 735-753. 10.1016/j.fgb.2004.03.001.PubMedView Article
- Guilleroux M, Osbourn A: Gene expression during infection of wheat roots by the 'take-all' fungus Gaeumannomyces graminis. Mol Plant Pathol. 2004, 5 (3): 203-216. 10.1111/j.1364-3703.2004.00219.x.PubMedView Article
- Broeker K, Bernard F, Moerschbacher BM: An EST library from Puccinia graminis f. sp. tritici reveals genes potentially involved in fungal differentiation. FEMS Microbiol Lett. 2006, 256 (2): 273-281. 10.1111/j.1574-6968.2006.00127.x.PubMedView Article
- Diatchenko L, Lau YF, Campbell AP, Chenchik A, Moqadam F, Huang B, Lukyanov S, Lukyanov K, Gurskaya N, Sverdlov ED: Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc Natl Acad Sci USA. 1996, 93 (12): 6025-6030. 10.1073/pnas.93.12.6025.PubMed CentralPubMedView Article
- Shearer G: Cloning and analysis of cDNA encoding an elongation factor 1 alpha from the dimorphic fungus Histoplasma capsulatum. Gene. 1995, 161 (1): 119-123. 10.1016/0378-1119(95)00269-C.PubMedView Article
- Pariasca JAT, Sunaga A, Miyazaki T, Hisaka H, Sonoda M, Nakagawa H, Sato T: Cloning of cDNAs encoding senescence-associated genes, ACC synthase and ACC oxidase from stored snow pea pods (Pisum sativum L. var saccharatum) and their expression during pod storage. Postharvest Biol and Technol. 2001, 22: 239-247. 10.1016/S0925-5214(01)00094-1.View Article
- Kikuchi S, Satoh K, Nagata T, Kawagashira N, Doi K, Kishimoto N, Yazaki J, Ishikawa M, Yamada H, Ooka H: Collection, mapping, and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301 (5631): 376-379. 10.1126/science.1081288.PubMedView Article
- Hugot K, Riviere MP, Moreilhon C, Dayem MA, Cozzitorto J, Arbiol G, Barbry P, Weiss C, Galiana E: Coordinated regulation of genes for secretion in tobacco at late developmental stages: association with resistance against oomycetes. Plant Physiol. 2004, 134 (2): 858-870. 10.1104/pp.103.034173.PubMed CentralPubMedView Article
- Oh IS, Park AR, Bae MS, Kwon SJ, Kim YS, Lee JE, Kang NY, Lee S, Cheong H, Park OK: Secretome analysis reveals an Arabidopsis lipase involved in defense against Alternaria brassicicola. Plant Cell. 2005, 17 (10): 2832-2847. 10.1105/tpc.105.034819.PubMed CentralPubMedView Article
- Ahn N, Kim S, Choi W, Im KH, Lee YH: Extracellular matrix protein gene, EMP1, is required for appressorium formation and pathogenicity of the rice blast fungus, Magnaporthe grisea. Mol Cells. 2004, 17 (1): 166-173.PubMed
- Talbot NJ, Ebbole DJ, Hamer JE: Identification and characterization of MPG1, a gene involved in pathogenicity from the rice blast fungus Magnaporthe grisea. Plant Cell. 1993, 5 (11): 1575-1590. 10.1105/tpc.5.11.1575.PubMed CentralPubMedView Article
- Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340 (4): 783-795. 10.1016/j.jmb.2004.05.028.PubMedView Article
- Katagiri T, Shinozaki K: Disruption of a gene encoding phosphatidic acid phosphatase causes abnormal phenotypes in cell growth and abnormal cytokinesis in Saccharomyces cerevisiae. Biochem Biophys Res Commun. 1998, 248 (1): 87-92. 10.1006/bbrc.1998.8914.PubMedView Article
- van Schooten B, Testerink C, Munnik T: Signalling diacylglycerol pyrophosphate, a new phosphatidic acid metabolite. Biochim Biophys Acta. 2006, 1761 (2): 151-159.PubMedView Article
- Rudd S: Expressed sequence tags: alternative or complement to whole genome sequences?. Trends Plant Sci. 2003, 8 (7): 321-329. 10.1016/S1360-1385(03)00131-6.PubMedView Article
- Kupfer DM, Drabenstot SD, Buchanan KL, Lai H, Zhu H, Dyer DW, Roe BA, Murphy JW: Introns and splicing elements of five diverse fungi. Eukaryot Cell. 2004, 3 (5): 1088-1100. 10.1128/EC.3.5.1088-1100.2004.PubMed CentralPubMedView Article
- Qutob D, Hraber PT, Sobral BWS, Gijzen M: Comparative analysis of expressed sequences in Phytophthora sojae. Plant Physiol. 2000, 123: 243-253. 10.1104/pp.123.1.243.PubMed CentralPubMedView Article
- Birch PR, Kamoun S: Studying interaction transcriptomes: coordinated analyses of gene expression during plant-microorganism interactions. New Technologies for Life Sciences: A Trends Guide. Edited by: Wood R. 2000, New York: Elsevier Science, 77-82.
- Randall TA, Dwyer RA, Huitema E, Beyer K, Cvitanich C, Kelkar H, Fong AM, Gates K, Roberts S, Yatzkan E: Large-scale gene discovery in the oomycete Phytophthora infestans reveals likely components of phytopathogenicity shared with true fungi. Mol Plant-Microbe Interact. 2005, 18 (3): 229-243. 10.1094/MPMI-18-0229.PubMedView Article
- Moons A: Ospdr9, which encodes a PDR-type ABC transporter, is induced by heavy metals, hypoxic stress and redox perturbations in rice roots. FEBS Lett. 2003, 553 (3): 370-376. 10.1016/S0014-5793(03)01060-3.PubMedView Article
- Stukkens Y, Bultreys A, Grec S, Trombik T, Vanham D, Boutry M: NpPDR1, a pleiotropic drug resistance-type ATP-binding cassette transporter from Nicotiana plumbaginifolia, plays a major role in plant pathogen defense. Plant Physiol. 2005, 139 (1): 341-352. 10.1104/pp.105.062372.PubMed CentralPubMedView Article
- Jasinski M, Stukkens Y, Degand H, Purnelle B, Marchand-Brynaert J, Boutry M: A plant plasma membrane ATP binding cassette-type transporter is involved in antifungal terpenoid secretion. Plant Cell. 2001, 13 (5): 1095-1107. 10.1105/tpc.13.5.1095.PubMed CentralPubMedView Article
- Snoeijers SS, P'erez-Garc'ýa A, Joosten MHAJ, Wit PJGMD: The effect of nitrogen on disease development and gene expression in bacterial and fungal plant pathogens. Eur J Plant Pathol. 2000, 106: 493-506. 10.1023/A:1008720704105.View Article
- Talbot NJ, McCafferty HRK, Ma M, Moore K, Hamer JE: Nitrogen starvation of the rice blast fungus Magnaporthe grisea may act as an environmental cue for disease symptom expression. Physiol Mol Plant Pathol. 1997, 50 (3): 179-195. 10.1006/pmpp.1997.0081.View Article
- Lau G, Hamer JE: Regulatory genes controlling MPG1 expression and pathogenicity in the rice blast fungus Magnaporthe grisea. Plant Cell. 1996, 8 (5): 771-781. 10.1105/tpc.8.5.771.PubMed CentralPubMedView Article
- Castelli V, Aury JM, Jaillon O, Wincker P, Clepet C, Menard M, Cruaud C, Quetier F, Scarpelli C, Schachter V: Whole genome sequence comparisons and "full-length" cDNA sequences: a combined approach to evaluate and improve Arabidopsis genome annotation. Genome Res. 2004, 14 (3): 406-413. 10.1101/gr.1515604.PubMed CentralPubMedView Article
- Haas BJ, Volfovsky N, Town CD, Troukhan M, Alexandrov N, Feldmann KA, Flavell RB, White O, Salzberg SL: Full-length messenger RNA sequences greatly improve genome annotation. Genome Biol. 2002, 3 (6): RESEARCH0029-10.1186/gb-2002-3-6-research0029.PubMed CentralPubMedView Article
- Seki M, Narusaka M, Kamiya A, Ishida J, Satou M, Sakurai T, Nakajima M, Enju A, Akiyama K, Oono Y: Functional annotation of a full-length Arabidopsis cDNA collection. Science. 2002, 296 (5565): 141-145. 10.1126/science.1071006.PubMedView Article
- Yamada K, Lim J, Dale JM, Chen H, Shinn P, Palm CJ, Southwick AM, Wu HC, Kim C, Nguyen M: Empirical analysis of transcriptional activity in the Arabidopsis genome. Science. 2003, 302 (5646): 842-846. 10.1126/science.1088305.PubMedView Article
- Imanishi T, Itoh T, Suzuki Y, O'Donovan C, Fukuchi S, Koyanagi KO, Barrero RA, Tamura T, Yamaguchi-Kabata Y, Tanino M: Integrative annotation of 21,037 human genes validated by full-length cDNA clones. PLoS Biol. 2004, 2 (6): e162-10.1371/journal.pbio.0020162.PubMed CentralPubMedView Article
- Hayashizaki Y: RIKEN mouse genome encyclopedia. Mech Ageing Dev. 2003, 124 (1): 93-102. 10.1016/S0047-6374(02)00173-2.PubMedView Article
- Okazaki Y, Furuno M, Kasukawa T, Adachi J, Bono H, Kondo S, Nikaido I, Osato N, Saito R, Suzuki H: Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs. Nature. 2002, 420 (6915): 563-573. 10.1038/nature01266.PubMedView Article
- Davis KR, Ausubel FM: Characterization of elicitor-induced defense responses in suspension-cultured cells of Arabidopsis. Mol Plant-Microbe Interact. 1989, 2: 363-368.View Article
- Park J, Park B, Jung K, Jang S, Yu K, Choi J, Kong S, Kim S, Kim H, Kim JF: CFGP: a web-based, comparative fungal genomics platform. Nucleic Acids Res. 2008, 36 (36 Database): 562-571.
- Ewing B, Green P: Base-calling of automated sequencer traces using Phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.PubMedView Article
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.PubMed CentralPubMedView Article
- Drabenstot SD, Kupfer DM, White JD, Dyer DW, Roe BA, Buchanan KL, Murphy JW: FELINES: a utility for extracting and examining EST-defined introns and exons. Nucleic Acids Res. 2003, 31 (22): e141-10.1093/nar/gng141.PubMed CentralPubMedView Article
- Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14 (6): 1188-1190. 10.1101/gr.849004.PubMed CentralPubMedView Article
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta DeltaC(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.