Identification and analysis of in planta expressed genes of Magnaporthe oryzae
© Kim et al. 2010
Received: 4 September 2009
Accepted: 10 February 2010
Published: 10 February 2010
Skip to main content
© Kim et al. 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.
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.
Classification of 2,315 uniESTs according to origin of the sequences
Contigs (No. ESTs)
Infection library specific
Subtraction library specific
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.
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.
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.
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.
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 .
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.
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.
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.
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.