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.
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.
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.
One of the most daunting challenges in the post-genomic era is to identify all transcribed regions in a sequenced genome. Although computational programs have played an important role in genome annotation, highly accurate prediction of the coding regions of many genes is still challenging . Therefore, the experimental approaches such as ESTs (expressed sequenced tags) , full-length cDNA sequencing , SAGE (serial analysis of gene expression) [4, 5], MPSS (massively parallel signature sequencing) [6, 7], RATE (robust analysis of transcript ends  and microarrays [9, 10] are essential tools to validate the annotated putative transcriptional units (TUs) in the sequenced genomes.
Although EST and full-length cDNA sequencing are effective approaches for gene discovery and expression profiling, neither is cost-effective nor comprehensive enough to isolate the majority of rare transcripts. Since the fabrication of most types of microarrays relies on sequences that are derived from ESTs, full-length cDNAs or computer program-predicted coding regions, many rare transcripts may be absent from these arrays. SAGE and MPSS are two powerful methods for genome-wide expression profiling and novel gene identification [4, 5, 7]. Both methods determine the sequence of short tags derived from a defined position at the 3' regions of expressed mRNAs in a cell. SAGE uses traditional sequencing technologies, while MPSS uses a novel bead-based hybridization procedure . The main advantage of MPSS over SAGE is that over one million signatures can be obtained from a library in comparison to 100,000 to 200,000 tags in a SAGE library. However, SAGE method is technically simple and can be performed in most molecular laboratories, whereas MPSS involves many complex cloning steps and library construction and sequencing is only performed at Solexa, Inc. The newly developed 5'-RATE method simplifies tag cloning and sequencing procedures and should have a broad application for transcriptome analysis of complex genomes .
Magnaporthe grisea, the causal agent of rice blast disease, causes severe yield losses to cereals, including the economically important crops, such as rice, barley, wheat, and millets . This fungus completes its life cycle by conidia landing on the leaf surface, conidia germination, appressorium formation, penetration, establishment and new conidia formation within seven to ten days after infection [11, 12]. In addition to its economic importance, the recent release of the complete genome draft sequence  and characterization of several important avirulence/pathogenecity-related genes , makes M. grisea a model plant pathogen for studying host-pathogen interactions. Although 12,841 genes (Annotation release 5, ) are predicted in the genome of M. grisea, many of these genes have no experimental support. Identification and characterization of all expressed genes from the infectious (appressoria) and non-infectious (vegetative mycelia) structures are indispensable for understanding the mechanism of fungal pathogenesis and designing novel strategies to combat the disease.
The main objective of this study was to identify all expressed transcripts encoded in the M. grisea genome using MPSS, RL-SAGE and oligoarrays, and to compare the qualitative and quantitative measurements produced by these three methods. Using these three technologies, the M. grisea transcriptome in the vegetative mycelial stage and in the infectious appressorial stage was analyzed in depth, in which many novel sense and antisense transcripts and alternative splicing variants were discovered. We also compared the correlation coefficients of the gene expression patterns revealed by MPSS, RL-SAGE and oligoarray. Our results provide the most comprehensive analysis of the sense, antisense and alternative transcripts of a fungal genome published to date and represent a useful genomic resource for further detailed functional analysis of the expressed transcripts involved in M. grisea growth and pathogenesis.
Characterization of mycelium and appressorium transcriptomes of M. grisea by MPSS and RL-SAGE methods
Using MPSS, a total of 1,326,720 transcript tags from mycelia and 1,459,146 transcript tags from appressoria were obtained. As per Meyers et al. , we filtered the total MPSS tags into two groups: significant (= 4 copies in either library) and non-significant (= 3 copies in both libraries) tags. After performing a clustering analysis with the total tags, 12,531 (61.7%) and 12,927(76.7%) significant tags and 7,768 (38.3%) and 3,500 (21.3%) non-significant tags were identified from mycelial and appressorial tissues, respectively (Table 1). Clustering analysis of the significant tags from both mycelium and appressorium MPSS libraries resulted in 20,062 distinct tags corresponding to 8,669 annotated genes. Using both significant and non-significant tags, a total of 9,467 unique genes were identified. Approximately, 37.5% (7,531) of the appressorial significant tags and 35.6% (7,135) of the mycelial significant tags were specific to each stage. When the tags matched to the annotated genes, 1,784 and 2,088 genes specific to appressorium and mycelium tissues, respectively, were identified [see Additional File 1]. These results show that the transcriptomes of appressorium and mycelium tissues are quite different.
An RL-SAGE library was constructed using the same RNA sample that was used for the construction of the mycelium MPSS library (Table 1). The average insert size of the RL-SAGE clones was approximately 1.1 kb. A total of 245,873 individual tags and 51,925 distinct tags were obtained from 7,292 sequence reads. Out of the total distinct tags, 16,580 (31.9%) were significant (= 2 copies) and 35,345 (68.1%) were non-significant (single copy) transcript tags.
Matching analysis of the MPSS and RL-SAGE tags to M. grisea genomic and expressed sequences
We compared the MPSS and RL-SAGE data to both genomic and expressed sequences of M. grisea to identify the genes from which these tags were obtained. The whole genome sequence from the Broad Institute (Release 5.0 on January 2006 at ) and the ESTs from TIGR (Release 5.0 on April 28, 2004 ) and the COGEME database ) were used to match MPSS and RL-SAGE tags. The matching results are summarized below.
1) MPSS tags
The matching rate of the significant tags to the genomic sequence from both libraries was about 85% (Table 2). In contrast, the matching rate to M. grisea ESTs for the significant tags from the mycelial (61.8%) and appressorial (51.2%) libraries was much lower than that to the genomic sequence. Similarly, 40 to 50% of the significant tags matched to the predicted coding sequences (CDS) . These results suggest that a large proportion of significant MPSS tags (at least 50% of tags in both the libraries) were not present in the public EST collections and could be considered as novel transcripts. Given that 5–10% of the ESTs may have no matches in the genome sequence, we conservatively estimate that at least 40% of the significant MPSS tags represent previously unidentified transcripts. We also observed that ~ 25 to 32% of the significant MPSS tags from both the libraries matched to regions within 500 bp downstream of the 3' region of annotated genes and about 7% of the significant tags from both the libraries matched to regions within 500 bp upstream of the 5' region of annotated genes. About 10 to 15% of the significant tags located in the intergenic regions matched to the genomic sequence but not to the CDS, 500 bp downstream or 500 bp upstream sequences (Table 2). Combining the CDS, 500 bp upstream and 500 bp downstream of CDS, we identified 7,686 genes in appressoria and 6,735 genes in mycelia by MPSS. A total of 3,362 genes were commonly expressed in both appressoria and mycelia. While 1,784 and 2,088 genes were specifically expressed at the appressorial and mycelial stages, respectively.
2) RL-SAGE tags
In comparison with the MPSS libraries, more significant tags were identified in the mycelial RL-SAGE library. About 65% of significant RL-SAGE tags matched the M. grisea genomic sequences (Table 2). Unlike the MPSS tags, only 31% (5,067) of the significant tags from the RL-SAGE library matched to the M. grisea EST sequences (TIGR) and a small portion (16%) of the tags matched to the CDS of the annotated genes (Table 2). We also observed that about 31% of the mycelial RL-SAGE significant tags matched to the 500 bp downstream of the 3' region of the annotated genes and about 6% of the significant tags matched within 500 bp upstream of the 5' region of the annotated genes (Table 2). The reason for the high percentage of the RL-SAGE tags located within the 500 bp downstream region might be due to more NlaIII sites (RL-SAGE) than DpnII sites (MPSS) in the 3' UTR region. Combining the tags hit to CDS, 500 bp upstream and 500 bp downstream regions, we identified 6,028 genes in mycelia by the RL-SAGE method.
3) Identification of novel transcript tags in the MPSS and RL-SAGE libraries
A detailed analysis of the novel transcripts that did not match to any sequences in the ESTs and annotated genes was performed (Table 2). The percentage of the novel tags in the appressorium MPSS, mycelium MPSS, and mycelium RL-SAGE libraries was about 35%, 26%, and 35%, respectively [see Additional File 2]. A total of 3,339 genes in the RL-SAGE library identified by CDS and their 500 bp upstream and down stream regions did not match the ESTs either in TIGR or COGEME database. Similarly, 3,186 and 2,298 genes identified in the MPSS appressorial and mycelial library, respectively, also did not show any matches in TIGR or COGEME M. grisea ESTs. In the two MPSS libraries, about 67% and 76% of the non-significant appressorial and mycelial tags matched to the genome sequence, respectively (Table 2). Even among the tags that did not match the genome sequence, many of them were significant tags, which might be derived from the un-sequenced or intron-exon junctions in the genome. Some of them may be true transcripts encoded in the genome but could not be matched due to sequencing errors in the genome or in the EST, MPSS or RL-SAGE tag sequencing. To validate the MPSS and RL-SAGE results, two genes (MGG_04847.5 and MGG_0490.5) without any EST matches in the public databases were amplified and cloned (data not shown), demonstrating that the majority of the identified novel tags might be true transcripts.
Antisense transcripts for the annotated genes of M. grisea
1) Antisense transcript tags from the MPSS libraries
Numerous antisense tags were identified in this study that matched M. grisea annotated genes in an antisense orientation. A total of 3,747 significant antisense tags (1,825 from the coding region, 1,452 from the 500 bp downstream regions, and 470 from the 500 bp upstream regions) corresponding to 2,958 genes from the mycelial library and 2,849 significant antisense tags (653 from the coding region, 1,879 from the 500 bp downstream regions, and 317 from the 500 bp upstream regions) corresponding to 2,629 genes from the appressorial library were identified. Of which, 1,586 and 1,257 genes were mycelium and appressorium-specific, respectively. Among them, 1,372 antisense genes were commonly expressed in both appressoria and mycelia. Total antisense tags from significant and non-significant tags are shown in Figure 1. Only fewer significant antisense tags (653 tags) were identified in the appressorium library as compared to the mycelial library (1,825 tags). Interestingly, we observed antisense transcripts for some of the well-known genes that are involved in appressorium formation and pathogenesis such as hydrophobin (MPG1), and calmodulin [see Additional File 3]. When the tags only matching the antisense sequences and having a single hit to the annotated genes were chosen, 232 genes in appressoria and 274 genes in mycelia were identified as antisense transcripts. The antisense tags and their frequency are listed in Additional File 4.
2) Antisense transcript tags from the RL-SAGE library
Although the total number of significant antisense transcript tags (3,558) corresponding to 3,100 genes identified from the mycelial RL-SAGE library was similar to that of the mycelial MPSS library, the tag distribution in the three regions of the annotated genes was different (Figure 1). The proportion of the RL-SAGE antisense tags located in the 500 bp downstream regions was twice (71.2%) that of the MPSS library (38.7%) at the mycelial stage. Similar with the sense tags, this difference might be due to the use of different anchoring enzymes in MPSS (Dpn II) and RL-SAGE (Nla III) library construction. Among the identified antisense tags, 364 were present in both RL-SAGE and MPSS libraries, 1,730 were specific to the MPSS libraries, and 431 were specific to the RL-SAGE library. The antisense tags and their frequency are listed in Additional File 4.
3) Annotated CDS regions with sense and antisense MPSS tags
Nearly 19% of the annotated genes in the M. grisea genome had sense and antisense tags matching to both strands of the annotated genes. Among these genes, 2,075 (16% of the total annotated genes) were from the mycelium MPSS library and 774 (6.0%) were from the appressorium MPSS library. There were 1,627 genes with sense and antisense pairs specifically present in the mycelium library and 326 genes with sense and antisense pairs specifically present in the appressorium library. Four hundred and forty eight genes (3.5%) were present in the two libraries [see Additional File 5]. To support that the MPSS sense and antisense tags are reliable, we found that six annotated genes have both sense and antisense transcripts from the RL-SAGE library and public EST collections (Table 3).
4) Annotated CDS regions with sense and antisense RL-SAGE tags
We identified that nearly 10% of the annotated genes in the M. grisea genome had bothsense and antisense tags from the mycelial RL-SAGE library(data not shown). The significant lower number of tag pairs at the protein coding region in the RL-SAGE library is likely because that the majority of the RL-SAGE tags were located in the 500 bp downstream region of the annotated genes (Figure 1). Some of the sense and antisense RL-SAGE tag pairs are also present in the MPSS libraries and public EST collections (Table 3).
Identification of alternative transcript tags in mycelia and appressoria
Since MPSS and RL-SAGE tags are derived from the 3' end of each transcript, the presence of more than one tag in the CDS of a gene suggests an alternative termination of the gene. These termination differences could be derived from either alternative polyadenylation or alternative splicing at the 3' end. All the transcripts from both cases are called alternative transcripts in this study. To determine the extent of transcriptional diversity in M. grisea, we assessed the proportion of the genes with alternative terminations in both libraries.
1) Alternative sense and antisense MPSS tags in mycelia and appressoria
When the sense tags were matched to the annotated genes, ~ 20–35% of the genes were found to have at least two alternative transcripts (Table 4). In general, more genes expressed in mycelia had more alternative transcript tags in comparison to that in appressoria. The hypothetical protein similar to linoleate diol synthase precursor (MGG_13239.5) in the mycelium MPSS library matched to 29 alternative tags and another hypothetical gene (MGG_01625.5) in the appressorium library had 12 alternative tags (data not shown). A total of 2,542 mycelial and 1,517 appressorial genes were found to encode alternative transcripts. Among them, 911 genes were common between the two tissues, 1,631 were mycelia-specific and 606 were appressoria-specific. Further, we used M. grisea TIGR EST database from which alternative splice forms from M. grisea were clustered. The genes undergoing alternative termination or splicing have more than one MPSS tags. A total of 47 appressorial and 55 mycelial unique clusters have two or more MPSS tags [see Additional File 6].
Similarly, many annotated genes had more than one antisense tags in the MPSS libraries (Table 4). Twenty percent of the annotated genes in appressoria and 32% in mycelia had more than one antisense tag. As with the sense transcript tags in mycelia, the same gene (MGG_13239.5) had thirteen alternative antisense tags (data not shown).
2) Alternative sense and antisense RL-SAGE tags in mycelia
About one-fourth (27.5%) of the annotated genes in mycelia were found to produce at least two alternative sense tags in the RL-SAGE library, which is less than that observed in the MPSS mycelial library (Table 4). In the RL-SAGE mycelial library, a hypothetical protein similar to reverse transcriptase (MGG_13890.5) was found to encode sixteen alternative transcript tags (data not shown). Many genes with known functions were found to encode alternative transcript tags [see Additional File 7]. A total of 10,629 alternative tags were commonly present in both RL-SAGE and MPSS [see Additional File 8].
Among the antisense RL-SAGE tags in the mycelium library, nearly a quarter (24%) had at least two alternative antisense transcript tags per gene (Table 4). A hypothetical gene (MGG_00329.5) was found to generate seven alternative antisense transcript tags (data not shown). Several genes with multiple sense and antisense alternative tags were also identified. For example, four sense and one antisense tags were obtained for the HSP70 gene (data not shown).
Characterization of the appressorium and mycelium transcriptomes by oligoarray hybridizations
To compare the transcriptional profiles generated from MPSS and RL-SAGE with that from oligoarray analysis, the same RNA samples used in MPSS and RL-SAGE library construction were hybridized to the M. grisea/O. sativa oligoarray . Using a stringent cut off at false discovery rate (FDR) = 0.05 that corresponds to a p- value of 0.001, 9,138 genes (43.9%) were identified to be statistically significant expressed in mycelium and appressorium tissues. Among them, 8,569 probes are from M. grisea genes and 569 probes from rice genes. The hybridizations with the rice genes were likely due to sequence similarity between housekeeping genes in both organisms. Of the 8,569 M. grisea genes, 4,652 (54%) and 3,917 (46%) were differentially (2 fold) up-regulated in appressoria, and down-regulated in mycelia, respectively. We identified 846 M. grisea genes that were ≥ 3.0 fold significantly up-regulated in appressoria, and 792 genes that were ≤ 3.0 fold significantly down-regulated in mycelia. The top 20 highly and specifically up- and down-regulated genes in appressoria and mycelia are shown in [see Additional File 9].
To gain more insight on the molecular mechanisms involved in M. grisea pathogenesis, we tried to functionally categorize 8,569 significant expressed genes, that were induced or repressed in the appressorium (4,652 genes) and also in the mycelium (3,917 genes), respectively, into different functional classes using KOGs analysis based on putative function of proteins . Functional classification and percentage of genes represented in appressorial and mycelial tissues are shown in the Additional File 10 and Additional File 11. The results indicate that a significant proportion of the genes (58% in appressoria and 64% in mycelia) were unclassified. The relative categories of genes expressed at the mycelial and appressorial stages are shown in Figure 2. The abundance of the gene category of cell cycle control, cell division, chromosome partitioning, cytoskeleton, lipid transport, and metabolism was significantly high in appressoria as compared to mycelia (Figure 2). On the contrary, the abundance of the genes involved in translation, ribosomal structure and biogenesis, RNA processing and modification, nuclear structure, coenzyme transport and metabolism, was significantly expressed in the mycelia when compared to appressoria.
Qualitative analysis of the M. grisea transcriptome using MPSS, RL-SAGE, and oligoarray platforms
To compare the expression pattern of the transcripts identified from the three different methods, a comparative analysis was performed. Because of the differences in the three platforms, a direct comparison of gene expression was not feasible. Therefore, an indirect comparison was made in which the oligoarray sequences (13,666) were used as the common targets to match the mycelial MPSS and RL-SAGE tags. The appressorium tags were not included in the analysis since no appressorial RL-SAGE library was made. Using the mycelial RL-SAGE (51,925) and MPSS (20,299) tags, 5,720 MPSS and 3,824 RL-SAGE tags matched the gene sequences on the microarray (13,660). When the expression profiles of the mycelial and appressorial tissues were compared, 3,793 genes were significantly expressed in the mycelial tissue. Clustering analysis of the genes from the three groups resulted in the identification of 7,741 genes. A total of 1,521 genes were commonly expressed in all three platforms (Figure 3). Oligoarray and RL-SAGE shared 32.58%, oligoarray and MPSS shared 34.30% and MPSS and RL-SAGE shared about 41.83% of the 7,741 genes. Oligoarray and MPSS together identified 91.49%, oligoarray and RL-SAGE together identified 74.21% and RL-SAGE and MPSS together identified 86.92% of the 7,741 unique genes. These results demonstrated that each method could identify specific groups of expressed genes and a combination of either of the two methods can identify the majority of the transcripts.
Quantitative analysis of the mycelium and appressorium transcriptomes using MPSS, RL-SAGE, and oligoarray platforms
To quantitatively assess the transcripts identified in the three different platforms, a Pearson correlation coefficient analysis was performed. In a comparative analysis between MPSS or RL-SAGE mycelial tags and the oligoarray mycelial mean signal intensity, a poor correlation was observed when unfiltered data was used (Table 5). In the MPSS and oligoarray data comparison, regression analysis was performed after removing the ten genes that were within the regression standard residual cutoff of = -1.75 to = 5.3. The removal of this small number of outliers had increased the correlation co-efficiency from 0.18 to 0.51 (Table 5). Similarly, after removal of fourteen outlying genes (based on a residual cutoff of = -1.75 to = 5.3), the correlation coefficiency between RL-SAGE tag frequency and the oligoarray mean signal intensity was improved significantly (Table 5). For example, the correlation coefficiency was increased from 0.29 to 0.64 when RL-SAGE tags with copy number = 10 were used. In general, the correlation coefficiency was higher in the comparisons for the genes with high copy numbers.
Conversely, a moderate correlation coefficiency (0.59) was obtained between MPSS and the oligoarray appressorial data without data filtering (Table 5). After the removal of the four outlying genes, the correlation coefficiency was increased from 0.65 to 0.74 for the MPSS tags with copy number = 10 (Table 5). In mycelia, the correlation was only about 0.5 for the genes with 25 or more copies, even after the removal of ten outlying genes. Comparing mycelial MPSS and RL-SAGE tags data, a low correlation was observed without data filtering. After the removal of the 18 outliers from the dataset, the correlation coefficiency was increased from 0.068 to 0.4 between mycelial MPSS and RL-SAGE tags (Table 5).
In summary, we found a low to moderate correlation among the expression data from the three platforms, especially those data between MPSS and oligoarray. In general, a better expression correlation was observed for high copy number tags in the MPSS and RL-SAGE libraries with their corresponding genes on the oligoarray.
Recent technological innovations have advanced genomics in an unprecedented way. Several complex genomes have been sequenced in recent years providing an excellent starting point to fully understand the genetic blueprint of an organism. However, identification of all the expressed portions of a sequenced genome is a challenging task, yet critical to the understanding of gene regulation and metabolic networks. The public availability of the whole genome sequence of M. grisea has established a solid foundation to further understand the pathogenicity mechanism of this notorious fungal plant pathogen which causes severe yield losses in rice growing countries . Elucidating the transcriptome of M. grisea may ultimately lead to the development of novel approaches for combating rice blast disease. In the last few years, many researchers have adopted various gene expression profiling techniques to characterize the M. grisea transcriptome under various conditions or in different cell types including EST sequencing , microarrays , and SAGE . However, these methods have only provided partial information about the M. grisea transcriptome due to technical limitations and the depth of the surveys performed in these studies. In this study, we employed three global and quantitative expression tools, namely MPSS, RL-SAGE and oligoarrays, to profile the M. grisea transcriptome at two developmental stages. In these experiments, the same RNA samples isolated from mycelia and appressoria were used so that the results from three different platforms could be readily compared. A total of 12,531 and 16,580 significant tags in mycelia have been identified by MPSS and RL-SAGE, respectively. In appressoria, 12,927 significant MPSS tags were identified. Many identified transcripts were not present in the existing EST or cDNA collections of M. grisea and many of them matched unannotated regions of the genome.
Both RL-SAGE and MPSS are tag-based approaches for transcriptome analysis and genome annotation. They are different from the conventional approach that focuses on the large clone collections following the principle of collecting a "representative clone" for each gene. Although conventional approaches are useful to catalog the expressed genes in certain tissues, especially moderately or highly expressed genes, many weakly expressed genes might have been missed in these collections. More importantly, it is not possible for the conventional approaches to address the questions regarding the dynamics of transcriptional regulation and regulatory principles like alternative promoter usage and splicing . In contrast, RL-SAGE and MPSS methods isolate and sequence short tags (17–21 bp) from the 3' regions of most transcripts. At least 100,000 RL-SAGE tags or a million MPSS tags can be easily obtained from these libraries.
In this study, we obtained approximately 1.3 and 1.4 million tags from the mycelial and appressorial MPSS libraries, respectively. The matching rate of the significant tags from each library to the M. grisea draft genome sequence was about 85%, suggesting that the MPSS data have a very deep coverage of the transcriptome. In contrast, only about 50 to 60% of the significant tags matched to the existing EST collections in the public databases. A similar result was also obtained from the RL-SAGE library, suggesting that current M. grisea EST collections are incomplete. One possible explanation is that most M. grisea ESTs were sequenced from the 5' region of the transcripts , whereas most of the MPSS and RL-SAGE tags were derived from the 3' region of expressed genes. Due to sequencing cost limitations, we only sequenced only 7,000 clones and obtained a quarter million tags from the mycelium RL-SAGE library. The transcripts recovered from the MPSS and RL-SAGE methods were overlapping but not identical due to the use of different anchoring enzymes in the library construction. Interestingly, we found that the genome matching rate of RL-SAGE tags is lower than that of the MPSS tags. These are two possible reasons. First, sequencing errors might generate unmatched tags, especially for singleton tags. Second, the significant MPSS tags used for matching have = 4 copies whereas the significant RL-SAGE tags have = 2 copies, suggesting the MPSS tags selected for matching may be more reliable. The last reason is that most of the RL-SAGE tags matched putative 3' UTR region, which may frequently targeted for RNA variation as reported in mammalian system . Nevertheless, our results demonstrated that MPSS and RL-SAGE methods are powerful techniques for deep transcriptome analysis and novel gene discovery. The two methods are complementary and different types of transcripts could be identified from each of these methods.
One of the advantages of tag-based techniques is the detection of alternatively terminated transcripts in the RNA population. From both MPSS and RL-SAGE libraries, we found many annotated genes have alternative transcript tags. Some of them have corresponding EST transcripts. The percentage of genes with evidence of alternative termination ranged from 27% in the RL-SAGE mycelium library to 35% in the MPSS mycelium library. The higher rate of alternative transcript tags in the MPSS library may be due to the fact that more MPSS tags (66%) matched to the coding regions of the annotated genes than that of the RL-SAGE tags (37%). It has been previously shown that a high rate of alternative transcripts was found in the protein coding regions (74%), and a low rate (4%) of alternative transcripts was found in the 3' UTR . The reason for the lower rate of alternative transcripts in appressoria than in mycelia is not clear. In addition, cloning and sequencing confirmation of some alternative transcripts without ETS support is required. Regulation of the alternative transcripts and functions of these sense and antisense alternative transcripts in M. grisea warrant more detailed analyses.
Naturally occurring antisense transcripts were first observed in prokaryotes and viruses and later found in eukaryotes. There is evidence for the involvement of antisense transcripts in alternative transcription [26, 27], RNA editing [28, 29], DNA methylation [30, 31], genomic imprinting [32, 33], and X-chromosome inactivation . In this study, many antisense tags were identified in mycelium or appressorium libraries, corresponding to more than 25% of the annotated genes. It is not clear at this point that why M. grisea genome encode so many antisense transcripts. Like the sense tags, almost two-thirds of antisense transcripts detected in the RL-SAGE library were located in the 500 bp downstream regions, whereas the majority of antisense tags from the MPSS library were located within the protein coding regions. Surprisingly, we found that almost one-fifth of the annotated genes encode both sense and antisense transcripts. In these genes, antisense transcripts could form double stranded RNA (dsRNA) with their sense transcripts. If a dsRNA is formed, it could be degraded to form small interfering RNAs that could decrease sense RNA abundance . Alternatively, interference by RNA polymerase II transcription activity on the antisense strand could restrict sense-strand transcription .
Recently, few studies demonstrated the function of antisense transcripts in fungal growth and development. For example, the circadian clock gene in the fungus Neurospora, a close relative of M. grisea, is regulated by the presence of natural antisense transcripts [37, 38]. Casas-Flores et al  expressed an antisense version of the pkr-1 gene of Trichoderma atroviride, encoding the regulatory subunit of protein kinase A (PKA), resulted in a non-sporulating phenotype.
Through data mining for EST, MPSS and RL-SAGE antisense tags, we have identified longer antisense transcripts in M. grisea for transcription factors MST12 (MGG_12958.5) and DEAD box-containing protein (MGG_12894.5), ribosomal protein, S9 (MGG_12892.5), ribosomal protein, L34 (MGG_05296.5) and ATP synthase alpha chain (MGG_07752.5). MST12 (a yeast homologs of STE12) is essential for host penetration and invasive growth, but was not required for appressorium formation [40, 41]. MST12 is regulated by MAP kinase at the downstream of signal transduction cascade during pathogenesis [40, 41]. It will be interesting to know if the MST12-mediated signal transduction cascade is regulated by antisense mechanism or not. A detailed characterization of MST12 antisense transcript may reveal its novel role in pathogenesis.
MPSS, SAGE, and oligoarrays are three widely used methods for transcriptome profiling. We performed qualitative and quantitative comparative analysis of the mycelial and appressorial transcriptomes revealed by the three methods. More than 40% of the annotated genes were detected by both MPSS and RL-SAGE methods. There was a good correlation in gene expression levels between the appressorium MPSS expression data and appressorium oligoarray data (0.67) and a moderate correlation between the mycelium MPSS and the mycelium oligoarray data (0.49) after removing several outlying genes in the datasets. However, the correlation between RL-SAGE and MPSS or oligoarray data was not significant. The low correlation between RL-SAGE and oligoarray might be because that oligoarray probes are designed from the protein coding regions of the annotated genes and a large number of RL-SAGE tags are located in the 500 bp downstream regions (putative 3' UTRs). We speculate that a low correlation between MPSS and RL-SAGE is due to the following two reasons. First, the use of different anchoring enzymes would change the location of the MPSS and RL-SAGE tags within a given transcript. We found that the majority of the MPSS tags matched the protein coding regions and the majority of RL-SAGE tags matched within 500 bp downstream (putative 3' UTRs). In some cases, because the M. grisea annotation is incomplete, the MPSS and RL-SAGE tags derived from the same transcript may be mapped to two different predicted genes. Second, MPSS and RL-SAGE use two completely different library construction and sequencing procedures. Any bias in the PCR amplification of synthesized cDNAs could lead to generation of different tag populations. Nonetheless, these data provide the first detailed analysis of transcriptional activity in an important fungal pathogen of plants, and constitutes a starting point for large-scale functional analysis of many novel fungal genes identified in the study.
We sequenced one RL-SAGE library of mycelia and two MPSS libraries of appressoria and mycelia of M. grisea. Using the same RNA samples of appressoria and mycelia, oligoarray hybridization was performed to check if these three approaches can detect similar sets of expressed genes in M. grisea. The distinct transcripts detected by MPSS and RL-SAGE in appressoria and mycelia ranged from 12,000 to 16,000, which correspond to about 9,000 genes, representing 80% of the predicted genes in M. grisea . A low to moderate correlation among the expression data from the three platforms was observed. MPSS and RL-SAGE methods identified many novel sense and antisense transcripts, which are differentially expressed at the two important growth stages of M. grisea. The identified novel transcripts, especially those specifically expressed in appressoria, are valuable genomic resource for a better understanding of the molecular basis of M. grisea pathogenicity. The established MPSS and RL-SAGE websites provide useful genomics resource for the public. The microarray (GSM126989) and RL-SAGE (GSM127012) data were deposited at the NCBI-GEO website and also the MPSS data can be downloaded at the M. grisea MPSS website .
Fungal strains, growth conditions and RNA isolation
The M. grisea strain 70–15 was chosen for the transcriptome profiling because of the availability of its whole genome sequence . The mycelia of 70–15 was cultured on a liquid medium [0.2% (w/v) yeast extract and 1% (w/v) sucrose] for 72 h (28°C at 200 rpm). The harvested mycelia were filtered and grinded for RNA isolation using the TRIzol method (Invitrogen, CA). For isolation of total RNA from the germinating appressoria, the mycelia of 70–15 were grown for two weeks on oatmeal agar plates and then the conidia were induced under white fluorescence light for five days. About 2 ml of conidia suspension (5 × 105 spores ml-1) was sprayed on Falcon Petri plates (150 mm × 15 mm) (Falcon, NJ) and the lids were covered with moist filter papers. These plates were incubated at 28°C and appressorium formation was monitored under a microscope at 6 h intervals. Fungal tissue was harvested 24 h after incubation since over 90% of the conidia extended a germ tube from the basal and/or the apical cell. Using a sterile blade, appressoria were scrapped and transferred quickly to the TRIzol solution, and then the suspension was centrifuged at 12,000 × g at 5°C for 5 min. About 2 g of appressorium pellet was collected and subjected for total RNA isolation. The poly (A+) mRNA was isolated from the total RNA using the Oligotex mRNA midi kit (Qiagen, CA).
MPSS library construction and sequence analysis
MPSS library construction was carried out at Solexa, Inc. (Hayward, CA) as described by Brenner et al  and Meyers et al . About 500 μg of total RNA isolated from mycelium and appressorium tissues (described above) were used in the MPSS library construction. The entire data set is available at . All tags were normalized to tags per million (TPM) as described by Meyers et al .
RL-SAGE library construction and sequence analysis
About 50 ng of mRNA isolated from mycelium tissue was used for RL-SAGE library construction as described by Gowda et al . A total of 7,292 sequence reads from the RL-SAGE library were sequenced at Arizona Genomics Institute. The ditags and distinct tags were extracted from these sequences using SAGEspy program developed at the Ohio Supercomputer Center . All RL-SAGE tags from mycelium library are available from the MGOS database .
Annotation of MPSS and RL-SAGE tags
The distinct MPSS and RL-SAGE tag sequences were matched to the M. grisea reference sequences including the genomic DNA, annotated genes (CDS), and 500 bp upstream (putative 5'UTR) and downstream (putative 3'UTR) regions that are available from the Broad Institute (version 5.0, release in January, 2006 ). The EST dataset of M. grisea from the TIGR database was used for matching MPSS and RL-SAGE tags release 5.0 on April 28, 2004 ). We also used a number of tools developed at the Ohio Supercomputer Center , the Magnaporthe grisea Oryza sativa (MGOS) database  and the University of Delaware MPSS database for data analyses. We identified the antisense transcripts from MPSS and RL-SAGE tags by converting all of the tags to antisense orientation using a reverse-complementation procedure before matching to the various sequences of M. grisea.
The M. grisea oligoarray and the hybridization procedures
The M. grisea oliogoarray chip containing 60-mer oligos representing a total of 21,885 probes was obtained from Agilent (G4137A; Wilmington, Delaware). Of the 21,885 probes, 13,666 are from the annotated genes of M. grisea and 7,144 probes are from the rice ESTs . The remaining 1,075 probes include quality controls, positive controls and negative controls. The total RNA was isolated from mycelia and appressoria using the TRIzol reagent (Invitrogen) following the manufacturers suggested protocol. Prior to hybridizations, quality and quantity of the total RNA sample was confirmed by running an agarose gel electrophoresis and by using a spectrophotometer. In the oligoarray hybridization experiments, we included six technical replicates of one RNA sample from mycelial and appressorial tissue, of which three were dye-reversal. About 500 ng of total RNA was used as template for cRNA production, and Cyanine dyes were incorporated using the Agilent low RNA input linear amp kit (5198-3523; Agilent). Normal yields from 500 ng total RNA input using an in vitro transcription were 15 μg cRNA (15 pmole cyanine dye incorporated/ug cRNA). One μg of labeled cRNA (Cy3 and Cy5 labeled sample) was diluted to 175 μl and defragmented at 60°C for 30 min following the Agilent hybridization protocols (5184-3568; Agilent). Defragmented samples were diluted to 500 μl (30% formamide final concentration) and hybridized for 20 h at 40°C. Arrays were washed, dried and scanned on an Agilent G2565BA microarray scanner described by . The raw TIFF images were analyzed using the Agilent Feature Extraction software v 8.1 using the recommended default settings.
Microarray data and KOG analysis
To minimize the variation in probe labeling and detection, intensities of Cy3- and Cy5-labeled probes were normalized using subgrid LOWESS normalization. Spots with lower signal intensity than the negative controls or with intensities less than twice the average background for the channel were manually blocked (flagged) from further analysis. This corrected and normalized dye bias data eliminated larger component of the variance, thus giving a greater confidence to the evaluated (treatment) and reference (control) data for testing the statistical significance. The genes with valid signals in all six replicates were exported to Partek Pro v.6.0 software (Partek Inc., Missouri). The normalized values were used to calculate the ratio of channel intensities (Cy5/Cy3), which were then log2 transformed. The transformed ratio was plotted in a histogram with ± two standard deviations away from the mean. A ± 1.7-fold increase or decrease in signal intensity or ± 0.77 on the Log2 scale from the histogram was considered to indicate genes that are differentially expressed. The normalized data were then subjected to ANOVA model using Boenferoni method [45, 46]. FDR (≤ 0.05) [47, 48] was calculated based on the p-value (≤ 0.001) from ANOVA. This stringent criterion limits the ability of the oligoarray experiment to detect small but biologically important changes between the appressorium and mycelium at approximately a 95% confidence interval. The genes that are significantly and differentially expressed in the appressorial and mycelial oligoarray data were used to compare the appressorial MPSS data with mycelial MPSS and RL SAGE data. Appressorium and mycelium-differentially expressed genes from the microarray analysis were functionally categorized using the euKaryotic Orthologous Groups (KOGs) database . The gene sequences were blasted against the KOG database with E-value of 1.0e-20.
Chen J, Sun M, Lee S, Zhou G, Rowley JD, Wang SM: Identifying novel transcripts and novel genes in the human genome by using novel SAGE tags. Proc Natl Acad Sci USA. 2002, 99: 12257-12262. 10.1073/pnas.192436499.
Adams MD, Kerlavage AR, Fleischmann RD, Fuldner RA, Bult CJ, Lee NH, Kirkness EF, Weinstock KG, Gocayne JD, White O, et al: Initial assessment of human gene diversity and expression patterns based upon 83 million nucleotides of cDNA sequence. Nature. 1995, 377: 3-174.
Kikuchi S, Satoh K, Nagata T, Kawagashira N, Doi K, Kishimoto N, Yazaki J, Ishikawa M, Yamada H, Ooka H, Hotta I, Kojima K, Namiki T, Ohneda E, Yahagi W, Suzuki K, Li CJ, Ohtsuki K, Shishiki T, Otomo Y, Murakami K, Iida Y, Sugano S, Fujimura T, Suzuki Y, Tsunoda Y, Kurosaki T, Kodama T, Masuda H, Kobayashi M, Xie Q, Lu M, Narikawa R, Sugiyama A, Mizuno K, Yokomizo S, Niikura J, Ikeda R, Ishibiki J, Kawamata M, Yoshimura A, Miura J, Kusumegi T, Oka M, Ryu R, Ueda M, Matsubara K, Kawai J, Carninci P, Adachi J, Aizawa K, Arakawa T, Fukuda S, Hara A, Hashizume W, Hayatsu N, Imotani K, Ishii Y, Itoh M, Kagawa I, Kondo S, Konno H, Miyazaki A, Osato N, Ota Y, Saito R, Sasaki D, Sato K, Shibata K, Shinagawa A, Shiraki T, Yoshino M, Hayashizaki Y, Yasunishi A: Collection, mapping, and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301: 376-379. 10.1126/science.1081288.
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: 890-897. 10.1104/pp.103.034496.
Brenner S, Johnson M, Bridgham J, Golda G, Lloyd DH, Johnson D, Luo S, McCurdy S, Foy M, Ewan M, Roth R, George D, Eletr S, Albrecht G, Vermaas E, Williams SR, Moon K, Burcham T, Pallas M, DuBridge RB, Kirchner J, Fearon K, Mao J, Corcoran K: Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays. Nat Biotechnol. 2000, 18: 630-634. 10.1038/76469.
Meyers BC, Tej SS, Vu TH, Haudenschild CD, Agrawal V, Edberg SB, Ghazal H, Decola S: The use of MPSS for whole-genome transcriptional analysis in Arabidopsis. Genome Res. 2004, 14: 1641-1653. 10.1101/gr.2275604.
Gowda M, Li H, Alessi J, Chen F, Pratt R, Wang GL: Robust analysis of 5'-transcript ends (5'-RATE): a novel technique for transcriptome analysis and genome annotation. Nucleic AcidsRes. 2006, 34: e126-10.1093/nar/gkl522.
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: 1337-1347.
Soanes DM, Talbot NJ: A bioinformatics tool for analysis of EST abundance during infection-related development by Magnaporthe grisea. Mol Plant Pathol. 2005, 6: 503-512. 10.1111/j.1364-3703.2005.00298.x.
Kim DD, Kim TT, Walsh T, Kobayashi Y, Matise TC, Buyske S, Gabriel A: A. Widespread RNA editing of embedded Alu elements in the human transcriptome. Genome Res. 2004, 14: 1719-1725. 10.1101/gr.2855504.
Tufarelli C, Stanley JA, Garrick D, Sharpe JA, Ayyub H, Wood WG, Higgs DR: Transcription of antisense RNA leading to gene silencing and methylation as a novel cause of human genetic disease. Nat Genet. 2003, 34: 157-165. 10.1038/ng1157.
Lewis A, Mitsuya K, Umlauf D, Smith P, Dean W, Walter J, Higgins M, Feil R, Reik W: Imprinting on distal chromosome 7 in the placenta involves repressive histone methylation independent of DNAmethylation. Nat Genet. 2004, 36: 1291-1295. 10.1038/ng1468.
Moore T, Constancia M, Zubair M, Bailleul B, Feil R, Sasaki H, Reik W: Multiple imprinted sense and antisense transcripts, differential methylation and tandem repeats in a putative imprinting control region upstream of mouse Igf2. Proc Natl Acad Sci USA. 1997, 94: 12509-12514. 10.1073/pnas.94.23.12509.
Eszterhas SK, Bouhassira EE, Martin DI, Fiering S: Transcriptional interference by independently regulated genes occurs in any relative arrangement of the genes and is influenced by chromosomal integration position. Mol Cell Biol. 2002, 22: 469-479. 10.1128/MCB.22.2.469-479.2002.
This work is supported by the Plant Genome Research Program of the National Science Foundation #0321437 and #0115642). We thank Dr. Jin-Rong Xu for his advice on the harvesting germinating appressoria and RNA isolation and Tae-Ho Lee for functional classification analysis using KOG programs.
Authors and Affiliations
Department of Plant Pathology, Ohio State University, Columbus, OH, 43210, USA
Malali Gowda, RC Venu, Mohan B Raghupathy, Huameng Li & Guo-Liang Wang
Department of Plant and Soil Sciences, University of Delaware, DE, 19711, USA
Kan Nobuta & Blake C Meyers
Arizona Genomics Institute, Department of Plant Sciences, University of Arizona, Tucson, AZ, 85721, USA
Ohio Supercomputer Center, The Ohio State University, Columbus, OH, 43210, USA
Little Falls Site, Agilent Technologies Inc, 2850 Centerville Road, Delaware, 19711, USA
Solexa, Inc., 25861 Industrial Blvd, Hayward, CA, 94545, USA
Christian D Haudenschild
Fungal Genomics Laboratory, Department of Plant Pathology, North Carolina State University, Raleigh, NC, 27606, USA
Department of Biological Science, Myongji University, Kyonggido, 449728, Korea
MG designed the study, performed the experiments, constructed RL-SAGE library, conducted data mining, and wrote the article. RCV performed MPSS, RL-SAGE, and oligoarray data analysis and assisted in writing the article, MBR performed oligoarray data analysis and revised the article, KN constructed the M. grisea MPSS database, HL and ES designed MPSS and RL-SAGE tag analysis software, HRK and RW sequenced the RL-SAGE libraries, SC performed oligoarray experiment, CDH performed MPSS library construction, RBD assisted in RL-SAGE clone sequencing, BHN performed KOG analysis, BM and GLW contributed to conception and design of the study and revising and writing the article. All authors read and approved the manuscript.
Malali Gowda, RC Venu contributed equally to this work.
Additional file 10: The up-regulated genes (4,649) in appressoria identified by oligoarray analysis were subjected for KOGs analysis based on the putative function in the KOGs protein database. Functional classification and percentage of genes represented in appressorial tissue is shown. (PPT 40 KB)
Additional file 11: The up-regulated genes (3,784) in mycelia identified by oligoarray analysis were subjected for KOGs analysis based on the putative function in the KOGs protein database. Functional classification and percentage of genes represented in mycelia tissue is shown. (PPT 41 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License (
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Gowda, M., Venu, R., Raghupathy, M.B. et al. Deep and comparative analysis of the mycelium and appressorium transcriptomes of Magnaporthe grisea using MPSS, RL-SAGE, and oligoarray methods.
BMC Genomics7, 310 (2006). https://doi.org/10.1186/1471-2164-7-310