Skip to main content

Rapid transcriptome characterization and parsing of sequences in a non-model host-pathogen interaction; pea-Sclerotinia sclerotiorum



White mold, caused by Sclerotinia sclerotiorum, is one of the most important diseases of pea (Pisum sativum L.), however, little is known about the genetics and biochemistry of this interaction. Identification of genes underlying resistance in the host or pathogenicity and virulence factors in the pathogen will increase our knowledge of the pea-S. sclerotiorum interaction and facilitate the introgression of new resistance genes into commercial pea varieties. Although the S. sclerotiorum genome sequence is available, no pea genome is available, due in part to its large genome size (~3500 Mb) and extensive repeated motifs. Here we present an EST data set specific to the interaction between S. sclerotiorum and pea, and a method to distinguish pathogen and host sequences without a species-specific reference genome.


10,158 contigs were obtained by de novo assembly of 128,720 high-quality reads generated by 454 pyrosequencing of the pea-S. sclerotiorum interactome. A method based on the tBLASTx program was modified to distinguish pea and S. sclerotiorum ESTs. To test this strategy, a mixture of known ESTs (18,490 pea and 17,198 S. sclerotiorum ESTs) from public databases were pooled and parsed; the tBLASTx method successfully separated 90.1% of the artificial EST mix with 99.9% accuracy. The tBLASTx method successfully parsed 89.4% of the 454-derived EST contigs, as validated by PCR, into pea (6,299 contigs) and S. sclerotiorum (2,780 contigs) categories. Two thousand eight hundred and forty pea ESTs and 996 S. sclerotiorum ESTs were predicted to be expressed specifically during the pea-S. sclerotiorum interaction as determined by homology search against 81,449 pea ESTs (from flowers, leaves, cotyledons, epi- and hypocotyl, and etiolated and light treated etiolated seedlings) and 57,751 S. sclerotiorum ESTs (from mycelia at neutral pH, developing apothecia and developing sclerotia). Among those ESTs specifically expressed, 277 (9.8%) pea ESTs were predicted to be involved in plant defense and response to biotic or abiotic stress, and 93 (9.3%) S. sclerotiorum ESTs were predicted to be involved in pathogenicity/virulence. Additionally, 142 S. sclerotiorum ESTs were identified as secretory/signal peptides of which only 21 were previously reported.


We present and characterize an EST resource specific to the pea-S. sclerotiorum interaction. Additionally, the tBLASTx method used to parse S. sclerotiorum and pea ESTs was demonstrated to be a reliable and accurate method to distinguish ESTs without a reference genome.


White mold, caused by Sclerotinia sclerotiorum (Lib.) de Bary, is a devastating disease of over 400 reported dicotyledonous hosts [1]. The disease causes economically significant losses of many crop plants including pea (Pisum sativum L.) under the appropriate environmental conditions [2]. Currently, little is known about the genetic control of pathogenicity in the fungus and mechanisms of resistance in pea. Although hundreds of pea cultivars have been screened for white mold resistance in replicated greenhouse and laboratory tests [3], only partial resistance has been identified to date.

The identification of genes underlying S. sclerotiorum pathogenicity and resistance in pea would increase our knowledge of the pea-S. sclerotiorum interaction and facilitate the introgression of resistance into pea varieties. However, progress in these areas has been hampered by the lack of sequence information regarding the pea genome. Although other legume genomes, including the models Medicago truncatula, Lotus japonicus and economically important Glycine max (soybean) are available [4], Pisum sativum is still genome resource-poor in part due to the large genome size and large fraction of highly repetitive DNA [5].

The performance of Next-Generation sequencing (NGS) technologies continue to rise while costs continue to fall which enables researchers to conduct whole transcriptome sequencing (RNAseq) studies of interactions between plants and pathogenic fungi [6]. The application of NGS in plant-fungal interaction research promises to shorten the overall time of development of molecular genetic information necessary for functional and translational studies. However, RNAseq has rarely been used to study plant-pathogen interactions, particularly in non-model systems. One reason for this is the difficulty in distinguishing plant and fungal ESTs and even virus or viroid contamination, particularly when reference genomes are not available. Here we report novel transcriptome sequence information from the pea-S. sclerotiorum interaction obtained by 454 pyrosequencing and propose a method of rapid and efficient transcriptome characterization in a non-model species with little prior molecular information. We also report on the development and validation of a strategy to distinguish plant and fungal ESTs using the tBLASTx program and “proxy-reference” genomes in the absence of true reference genomes.


Contiguous EST assembly

10,158 contigs were obtained by de novo assembly of 128,720 high-quality reads produced on a Roche 454 GS FLX sequencer (see Additional file 1). Minimum contig length was 50 bp, maximum length was 1,015 bp and average length was 200 bp. Average read coverage of contigs was 4.5X, and the maximum read coverage was 2,303X (Figure 1).

Figure 1
figure 1

Depth of 454 pyrosequence read coverage (X) for the 10,158 pea- S. sclerotiorum transcriptome EST contigs.

Filter for virus or viroid contamination

The tBLASTx program identified 51 contigs with a BLAST hit (alignment identity) to virus or viroid DNA with an e-value cutoff threshold of less than 1e-3. Further evaluation of these 51 EST contigs with tBLASTx against 3 legume and 7 fungal genome databases (which acted as proxy-reference genomes) revealed that 46 contigs showed significant alignment with the proxy plant genome database, 43 showed significant alignment with the fungal genome database and 40 showed significant alignment with both databases, 2 contigs showed significant alignment only with virus genomes. By comparing the e-value ratio (virus/fungi or virus/plant) of all proxy-reference genome alignments, 10 contigs were assigned to pea and 9 contigs were assigned to S. sclerotiorum based on an e-value ratio >1e20, 30 contigs were difficult to distinguish, with e-value ratios between 1e-20 and 1e20. BLASTn analysis of the 32 unassigned contigs against pea and S. sclerotiorum ESTs from known sources revealed that 20 contigs, including the 2 that only aligned with virus genomes, had high identity matches to pea with 95% accuracy and 95% query coverage, and 1 to S. sclerotiorum. The 11 contigs which were difficult to assign to a genome database by either tBLASTx or BLASTn methods were far more similar to plant or fungi than to virus genomes.

Development and testing a method to distinguish pea and S. sclerotiorum ESTs using an artificially mixed pool

Pea and S. sclerotiorum ESTs were downloaded from GenBank to test the tBLASTx sorting method. Three hundred twenty-one ESTs with vector contamination and 71 ESTs highly similar to virus or viroids were removed from the total pool of 36,080 known pea and S. sclerotiorum ESTs. Using an e-value threshold of 1e-3, 35,688 mixed ESTs from pea and S. sclerotiorum were compared to legume and fungal proxy-reference genome databases and parsed using the tBLASTx program (Figure 2). 11,191 ESTs only aligned with the legume proxy-reference genomes, 11,259 ESTs only aligned with the fungal proxy-reference genomes, 11,266 ESTs similar to both plant and fungal proxy-reference genomes and 1,972 ESTs did not match either proxy-reference genome database. The ESTs with tBLASTx results to both plant and fungal genomes were analyzed further by comparing the e-value ratio (fungi/plant) of fungal and plant proxy-reference genome alignments. 4,098 ESTs were assigned to pea based on an e-value ratio >1e20, 5,649 ESTs were assigned to S. sclerotiorum based on an e-value ratio <1e-20, while 1,519 ESTs were difficult to distinguish due to high e-value alignments to both proxy-reference genome databases, with e-value ratios between 1e-20 and 1e20. This method successfully separated 90.1% of the known ESTs into pea or S. sclerotiorum categories, with only a 0.1% misallocation rate. Only 5.5% of ESTs had zero similarity to either of the proxy-reference genomes, and 4.3% of ESTs had high similarity to both the plant and fungal proxy-reference genome databases (Table 1).

Figure 2
figure 2

Flow diagram of host and pathogen EST parsing method for mixed transcriptomes from plant and fungi based on a modified tBLASTx method.

Table 1 Preliminary testing of the tBLASTx method with pooled known pea and S . sclerotiorum ESTs

Parsing 454 pyrosequence pea and S. sclerotiorum ESTs with tBLASTx and BLASTn

Initial tBLASTx parsing of 10,158 contigs with an e-value threshold of 1e-3 resulted in identification of 4,523 pea ESTs, 2,304 S. sclerotiorum ESTs, 1,974 ESTs that matched both pea and S. sclerotiorum, and 1,357 ESTs that did not matched either proxy-reference genome database. The 1,974 ESTs that matched both proxy-reference genomes at the e-value threshold of 1e-3 were further subdivided using the e-value ratio method into 544 pea ESTs (fungi/plant e-value ratio >1e20), 355 S. sclerotiorum ESTs (fungi/plant e-value ratio <1e-20) and 1,075 that were ambiguous (fungi/plant e-value ratio <1e20 and >1e-20). This brought the number of classified ESTs for each category to 5,067 for pea, 2,659 for S. sclerotiorum, 1,075 as ambiguous with high matches to both proxy-reference genomes, and 1,357 with no significant alignment. The remaining 2,432 EST contigs that were ambiguous or showed no significant alignment were further parsed with BLASTn analysis against known pea and S. sclerotiorum ESTs if identity and query coverage were both equal to or greater than 95%. 1,232 ESTs of this pool were assigned with BLASTn to pea and 121 ESTs were assigned to S. sclerotiorum, leaving 310 ambiguous and 769 EST contigs with no significant alignment. In total with tBLASTx and BLASTn, 10,158 contigs were separated into 6,299 pea ESTs, 2,780 S. sclerotiorum ESTs, 310 ambiguous ESTs and 769 unassigned ESTs (Figure 3).

Figure 3
figure 3

Distribution and assignment of 10,158 EST contigs derived from pyrosequencing of the pea- Sclerotinia sclerotiorum transcriptome using the tBLASTx method. EST contigs with an e-value score <1e-3 to only plant or fungal EST reference databases were automatically assigned. Those EST contigs with hits (<1e-3) to both databases were assigned based on a cutoff ratio of e-values (fungi/plant), those ESTs with a ratio > =1e20 were classified as pea, those < =1e-20 as Sclerotinia and those between >1e-20 and <1e20 as ambiguous. ESTs without any hits (>1e-3) were classified as unassigned ESTs.

Validation of tBLASTx and BLASTn EST parsing results by PCR

Validation of the tBLASTx and BLASTn assignment was performed for 50 S. sclerotiorum and 50 pea EST contigs randomly sampled from the two assigned categories. All 50 primer sets designed to the pea EST contigs amplified the expected amplicon size in both the pea-S. sclerotiorum and non-inoculated pea cDNA indicating correct parsing assignment of the pea ESTs (Figure 4). Of the 50 PCR primers designed to the S. sclerotiorum ESTs, 47 amplified a PCR product from both the pea-S. sclerotiorum and S. sclerotiorum only cDNA samples and most of them amplified the same size amplicon in both cDNA samples. Two of the 50 S. sclerotiorum PCR primer pairs amplified the expected PCR products from the pea-S. sclerotiorum cDNA sample but not the S. sclerotiorum only cDNA, perhaps indicating that this transcript is only expressed during the interaction with pea. One S. sclerotiorum primer set failed to amplify any PCR product from either template.

Figure 4
figure 4

PCR results for validation of 50 pea and 50S. sclerotiorum pyrosequence-classified EST contigs against a mix of host and pathogen interaction cDNA, host cDNA only and pathogen cDNA only. All 50 PCR pea primer sets amplified the expected amplicons from the pea-S. sclerotiorum and pea cDNA samples. Forty-seven of the 50 PCR S. sclerotiorum primer sets amplified expected amplicons from the pea-S. sclerotiorum and S. sclerotiorum cDNA samples, one primer pair failed to amplify from either cDNA and the remaining two only amplified product from the pea-S. sclerotiorum cDNA possibly indicating interaction induced transcripts.

Unique ESTs expressed in the pea-S. sclerotiorum interaction

To detect unique genes expressed in our pea-S. sclerotiorum interaction, the 6,299 classified pea ESTs in our data set were compared with BLASTn against 81,449 recently published pea ESTs from flowers, leaves, cotyledons, epi- and hypocotyl, and etiolated and light treated etiolated seedlings [7]. Of these 6,299 ESTs, 3,459 ESTs had significant alignments with an e-value cutoff of 1e-10, in which 1,668 contigs had a percentage identity threshold of 95% for 95% or more of the query sequence, leaving 2,840 potentially unique pea ESTs to the pea-S. sclerotiorum interaction. It was possible to annotate 1,631 of these ESTs of which 67 contigs encode transcription factors (Table 2), 69 were involved in signaling pathways (Table 3) and 82 contigs were involved in encoding defense-associated proteins (Table 4).

Table 2 Description of genes encoded by 67 unique pea ESTs encoding transcription factors
Table 3 Description of genes encoded by 69 unique pea ESTs involved in the signaling pathways
Table 4 Description of 82 unique pea ESTs encoding defense-associated proteins

The 2,780 S. sclerotiorum EST contigs were also assessed with BLASTn against 57,751 S. sclerotiorum ESTs (from mycelia at neutral pH, developing apothecia and developing sclerotia). Of these, 1,784 ESTs matched with an e-value cutoff of 1e-10, in which 294 ESTs matched with 95% identity for 95% of more of the query length to the S. sclerotiorum EST growth libraries. Of the remaining 996 unique ESTs, it was possible to annotate 438 ESTs of which 95 ESTs were described as being related to pathogen virulence or pathogenicity (Table 5).

Table 5 Description of 95 unique S . sclerotiorum EST contigs with putative involvement in virulence or pathogenicity

Prediction of secretory/signal peptides for the S. sclerotiorum contigs

A total of 2,754 coding regions were predicted with OrfPredictor from the set of 2,780 S. sclerotiorum ESTs. The peptide sequences were then used as a query for SignalP 3.0, which predicts the presence and location of signal peptide cleavage sites in amino acid sequences and identifies them as secretory proteins. The neural network (NN) method predicted 244 secretory signals, and the Hidden Markov Model (HMM) predicted 216. A total of 142 ESTs were identified by both NN and HMM and can be considered putative secretory peptides with high confidence (see Additional file 2). Of these 142 predicted secretory proteins, 21 were reported to be involved in pathogen virulence or pathogenicity (Table 6).

Table 6 S . sclerotiorum EST contigs encoding potential secretory/signal peptides involved in virulence or pathogenicity


Significance of study and summary of the main findings

Despite Pisum sativum being used by Gregor Mendel to propose a model of particulate inheritance and being a highly nutritious food source for populations worldwide, few genomic resources exist for pea. One of the pathogens of pea, S. sclerotiorum is not only capable of causing devastating disease of pea but is able to infect over 400 plant species [1]. By sequencing a normalized cDNA pool of the pea-S. sclerotiorum interaction with next generation sequencing we have catalogued a number of novel genes putatively involved in pathogenicity and resistance. To our knowledge this is the first study to examine the pea-S. sclerotiorum “interactome”. Sequencing the transcriptome (RNA-seq) is the method of choice in non-model systems for transcript discovery and genome annotation [8]. However, it has rarely been used to study plant-fungal interactions; one reason for this is the difficulty in distinguishing plant and fungal ESTs, particularly when reference genomes are not available. Using genomes of closely related species and tBLASTx to parse pea and S. sclerotiorum ESTs we demonstrated that Roche 454-pyrosequencing is a useful technique to characterize the host-pathogen interactome when genome resources are limited.

tBLASTx parsing method

Two different strategies have been utilized previously to identify transcript origins in mixed plant and fungal EST datasets. One is a predictive method based on triplet nucleotide usage frequencies [9] and the other is a homology method using the BLASTp algorithm [10]. One shortcoming of the BLASTp method is that it could not be applied to novel genes or sequences from the non-coding regions of genes. Although the triplet nucleotide frequency method extends the application of the algorithm to both coding and non-coding sequences, the classification accuracy is approximately 90%, and required the use of a training set of ESTs to develop the nucleotide frequency for separation. A combined method was also used by Fernandez et al. [11], although this method distinguished 91% of the ESTs from the Coffea arabica Hemileia vastatrix interaction no validation of the method was presented [11].

Classification of genes from a pool of mixed cDNA by traditional sequence similarity analysis (BLAST) is of interest to many investigations into plant-pathogen interactions. DNA sequencing is becoming more affordable and whole genome sequences of many organisms are becoming available and will aid in plant-pathogen interaction studies. However, in pea these resources are not available, therefore, we used a standalone BLAST approach against proxy-reference genome databases with high genetic similarity to pea or S. sclerotiorum to distinguish mixed transcripts. Using an artificial mixture of known pea and Sclerotinia ESTs, we found the error rate using the BLAST method was significantly lower than the triplet nucleotide frequencies method (Table 7). We also demonstrated that the tBLASTx algorithm provided improved sorting of contigs relative to the BLASTn algorithm, and results in fewer ambiguous reads (see Additional file 3). In addition, although one individual genome of S. sclerotiorum (strain 1980) has been sequenced [12], there are still 1.6 Mb of predicted gaps in the 39.6 Mb assembly. To avoid ignoring unique genes between two different strains of the same species, a multi-fungal genome approach was adopted in this study. It was demonstrated that the assignment error rate based on 7 closely related fungal genomes was slightly decreased relative to assignment based on the single S. sclerotiorum genome (see Additional file 4). The e-value and e-value ratio utilized in our study to differentiate pea and S. sclerotiorum reads chosen selected after comparing several e-values, to maximize discrimination while reducing the error rate (see Additional file 5). Additionally, we determined error rates for this method using the artificial EST mix and validated the technique using our EST data set. We found that the percentage of unassigned ESTs (23.9%) in the 454 data set was higher than in the test EST data set (9%). One hypothesis to explain this difference is the average sequence length in the 454 data (200 bp) was shorter than known pea (496 bp) or Sclerotinia ESTs (674 bp) used in test data, which may result in no significant alignment against the proxy-reference plant or fungal genome databases, particularly when non-coding mRNA is considered. The remaining unassigned EST contigs (21%) were parsed using BLASTn searches against known pea and S. sclerotiorum ESTs, which assigned 44.4% of the unassigned contigs. Using this combined tBLASTx and BLASTn approach 89.4% of the 10,158 contigs were identified as pea (6,299) or S. sclerotiorum EST (2,780). Additionally, the EST parsing method was validated by PCR demonstrating that the parsing method was able to correctly assign ESTs from the 454 data set with a low error rate.

Table 7 Comparison of triplet nucleotide frequencies method and tBlastx method to assign artificial EST mixture of pea- S . sclerotiorum (including 18,490 pea ESTs and 17,198S . sclerotiorum ESTs)

Pea ESTs unique to the pea-S. sclerotiorum interaction

In response to pathogen attack, plants have evolved complex signaling and defense pathways. Putatively unique ESTs in our pea-S. sclerotiorum interactome were defined and identified by comparing EST contigs in our library against those of non-interaction EST libraries of pea and S. sclerotiorum. Although we identified a total of 2,840 (45.1%) putatively unique pea ESTs it was only possible to annotate 1,631 of these and only 451 had annotations suggesting roles in defense or response to biotic and abiotic stress. Most of the annotated genes are consistent with previous expression profiling analyses in Brassica napus infected with Sclerotinia sclerotiorum[13]. Following infection, many genes, including those encoding defense-associated proteins, enzymes involved in signaling pathways, and genes encoding transcription factors were induced.

Transcriptional control of the expression of stress-responsive genes is a crucial part of plant response to a range of abiotic and biotic stresses [14]. We demonstrated that 67 putative transcription factors were detected. These genes were classified into the MYB family, the Apetala2/Ethylene responsive element binding protein (AP2/EREBP) family, WRKY family and others (Table 2). Seven MYB family transcription factors were detected in our data and they play a key role in hormone signal transduction and disease resistance [15]. Eight AP2/EREBP transcription factors, including 3 ethylene insensitive transcription factors (contig 1311, 7671 and 8261) and 3 AP2/ERF genes (7472, 8894 and 9486), are key regulatory elements for ethylene signaling and response for biotic or abiotic stresses [16, 17].WRKY40 act as negative regulators of defense signaling and have been associated with negatively regulating resistance to P. syringae in Arabidopsis [18].

Plant defenses are regulated through a complex network of transduction pathways [19]. Sixty-nine unique pea ESTs involved in signaling pathways were detected in this study. The signaling pathways were mediated by different signaling molecules, like abscisic acid (ABA), auxin, brassinosteroid, calcium ion, ethylene (ET), gibberellic acid (GA), jasmonic acid (JA), salicylic acid (SA) and small GTPase (Table 3). Those results were consistent with previous studies of signaling pathways involved in plant resistance to Sclerotinia sclerotiorum[20, 21].

Expression of downstream proteins, including defense-associated proteins, was induced through signal transduction and transcription factor regulation after pathogen infection. In this study, 82 unique pea ESTs encoding defense related proteins were detected (Table 4). Four contigs (5422, 6766, 7235 and 9781), encoding putative pathogenesis-related (PR) proteins involved in the response to pathogen attack were prominent. Numbers of cell-wall-related genes were also detected; those contigs involved in the biosynthesis of plant cell wall structures and the disassembly of fungal cell walls. Chitinase (encoding by contig 491, 589, 1622, 2311 and 5742), beta-1, 3-glucanase (encoding by 196, 1243 and 4388) and other glycoside hydrolases are known to possess anti-fungal activity by degrading fungal cell walls [22].

S. sclerotiorum ESTs unique to the pea-S. sclerotiorum interaction

Pathogens have evolved a number of strategies to gain entry into the host cell and to overcome the plant defense system. In this study, we identified 996 S. sclerotiorum contigs as specifically expressed during pea-S. sclerotiorum interaction through comparison of EST contigs against S. sclerotiorum ESTs from growth libraries. Ninety-five of 438 annotated contigs were described as being involved in pathogen virulence or pathogenicity (Table 5).

Fungi produce enzymes that degrade the cell wall and wall-associated polymers to penetrate plant cells. There were 39 specifically expressed contigs involved in the penetration of the plant cuticle and cell wall. Contig 6412 encodes an exoglucanase 2 precursor, which has cellulolytic activity [23] and is involved in cellulose degradation; enzymes encoded by 11 contigs (473, 555, 717, 970, 3115, 3513, 5562, 5918, 6759, 6383 and 9655) are involved in hemicellulose degradation; enzymes encoded by 11 contigs (1533, 3147, 3861, 5679, 6431, 6824, 7085, 7661, 7834, 8501 and 10069) are involved with pectin degradation. In addition, carbohydrate esterase encoded by contig 1562 was also involved in plant polysaccharide degradation. Integrity of the fungal cell wall is also very important for pathogenesis and some reports showed the deletion of biosynthetic cell wall enzymes resulted in dramatically reduced virulence [24]. In our data, 18 contigs were identified as affecting biosynthesis and integrity of fungal cell walls. Enzymes encoded by contigs 6654, 6797, 8474 and 8681 were involved in chitin synthesis; contigs 2424, 3362, 5737, 6956, 7397 and 9176 were involved in glucan synthesis; and 8 contigs (239, 330, 2930, 3825, 4843, 5773, 8785 and 9570) were involved in mannan synthesis.

Sclerotinia sclerotiorum differentiates appressoria into infection cushions prior to invasion and we found 12 genes involved in the formation of infection structures. Eight contigs were involved in response to the host immune system, of which 3 efflux transporters encoded by contigs 1220, 4783 and 6180 are responsible not only for export of compounds involved in pathogenesis such as secondary metabolites, but also export of host-derived antimicrobial compounds [2527]. Contig 1769 had similarity to the guanine nucleotide-binding protein (G protein) alpha subunit which is an important signal transducing molecule in cells, essential for growth, asexual and sexual development, and virulence in both animal and plant pathogenic filamentous fungal species [28]. Importin beta-2 encoded by contig 3623 belongs to the importin β family which mediates transport between the nucleus and cytoplasm of macromolecules that contain nuclear import or export signals. All importin β members have the ability to recognize and bind specific cargo involved in the recognition of the host and signaling [29, 30].

Secreted/signaling proteins

Proteins secreted by fungi play a key role in the development of plant disease and the evolution of pathogenicity [31]. Some secreted proteins can degrade polymers encountered, such as cellulose, lipid, protein, and lignin, and transport the resulting simple sugars, amino acids, and fatty acids into the growing cell for use [32]. Using the SignalP3.0 program with stringent criteria, 142 contigs encoding putative secreted proteins were identified in the 2,780 S. sclerotiorum contigs. Twenty-one of the 66 annotated contigs were described as involved in pathogen virulence/pathogenicity in previous research (Table 6). Contig 355 encodes an enolase which is usually present on the cell surface or even secreted and is a potential virulence factor. In bacterial systems enolase has been demonstrated to contribute to pathogenicity by binding plasminogen in the infected host, potentially allowing the bacteria to acquire surface-associated proteolytic activity [3335]. The basic leucine zipper transcription factor, encoded by contig 395, is a member of the bZIP family, one bZIP family member (Moatf1) from the rice fungus Magnaporthe oryzae mediates oxidative stress responses and is necessary for full virulence [36]. Contig 1352 encoding fkbp-type peptidylprolyl isomerase, with high homology to the Mip (macro-phage infectivity potentiator) protein, has been shown to be an essential virulence factor in Legionella pneumophila[3739]. Chitin synthase 1 (contig 1434) plays a major role in cell wall biogenesis. Disruption of Botrytis cinerea class I chitin synthase gene Bcchs1 results in cell wall weakening and reduced virulence [40, 41]. Autophagy is necessary for turnover of organic matter during the formation of conidia and appressoria and for normal development and pathogenicity in Magnaporthe grisea. Autophagy is required for the virulence of some eukaryotic pathogens [4244]. Contig 6759 encodes endo-β-1,4 xylanase which plays a significant role in the virulence of Magnaporthe oryzae, affects both penetration and expansion of M. oryzae in infected plants [45]. Pectin methylesterase (PME) produced by phytopathogenic bacteria and fungi catalyses the demethoxylation of pectin, a major plant cell wall polysaccharide [46]. The possible role of secreted adenylate kinase (AK), encoded by contig 9219, as a virulence factor is in producing and keeping an intact pool of toxic mixtures of AMP, ADP, and ATP, which allows Pseudomonas aeruginosa to exert its full virulence [47]. Glutathione reductase is important to nitric oxide and macrophage resistance and is essential for virulence [48] and in Candida albicans GRX2, a putative glutaredoxin, is required for virulence in a murine model [49].


Here we present an EST resource that is specific for the pea-S. sclerotiorum interaction. We demonstrate and validate a method to reliably parse host and pathogen ESTs without the need for reference genomes. The ESTs were compared to non-interaction EST libraries to identify candidate resistance and pathogenicity genes. We also catalogued 145 proteins putatively secreted by S. sclerotiorum. The EST dataset will be a useful reference for further plant-fungus interaction studies, particularly for the Sclerotinia and legume research communities. Additionally, the S. sclerotiorum ESTs will be a valuable resource for the annotation of the S. sclerotiorum genome. Although the depth of our sequencing was not sufficient to obtain a global view of transcripts expressed during the pea-S. sclerotiorum interaction, the results are still very useful for the identification of plant resistance, fungal pathogenicity and virulence genes. This study sets the ground work and will be a resource for our current pea-S. sclerotiorum RNAseq expression profiling studies.


Plant, fungal growth and inoculation

Three plants of pea cultivar ‘Lifter’ (PI628276) were established per 1 gallon plastic pot in Sunshine LA 4 potting mix (Sun Gro Horticulture, Bellevue, WA). The plants were maintained in a greenhouse for 4 weeks with supplemental lighting extending the day length to approximately 14 h (October). Day and night temperatures were 22 ± 2°C and 16 ± 2°C, respectively. S. sclerotiorum isolate WMA-1 was isolated from a diseased pea plant in 2003 from a pea field (Washington, USA) with white mold disease symptoms and stored as air dry sclerotia at room temperature. Isolate WMA-1 (=ATCC MYA-4521) was demonstrated to be genetically representative of eight S. sclerotiorum strains sampled from legume hosts from various geographic locations using randomly amplified polymorphic DNA (RAPD) analysis (Kawabe and Peever, unpublished). Plants were inoculated with a 5 mm plug collected from the leading edge of an actively growing colony on a potato dextrose agar (PDA). The plug was placed fungal side down on the stem between the 4th and 5th detectable nodes and held in place by wrapping with Parafilm. Plants were transferred to a growth chamber with a 12 h photoperiod, an approximate 60% relative humidity, temperature of 20 ± 1°C and a 12 h photoperiod, for 72 hours to allow disease lesion development prior to RNA extraction.

Total RNA extraction and purification of mRNA from total RNA

A 1 cm stem section was collected from each of 18 infected plants by cutting above and below the lesion front advancing toward the base of the plant. The stem section included both necrotic and green tissue with the advancing lesion front located in the center of the section. Stem sections were snap-frozen in liquid nitrogen and ground to a fine powder with a mortar and pestle. A total of 3 ml of TRIzol (Invitrogen, Carlsbad, CA, USA) was added to the ground tissue and the sample was split in half for column purification with the TRIzol Plus RNA purification kit (Invitrogen, Carlsbad, CA, USA). The additional step of on-column DNA digestion was performed with DNase I (Invitrogen, Carlsbad, CA, USA) to remove contaminating DNA. RNA was eluted in 250 μl of water per spin column. Poly-A RNA was isolated from total RNA with the Oligotex kit using the mRNA spin-column protocol (Qiagen, Valencia, CA, USA). Purified mRNA was eluted in a total of 100 μl of 5 mM Tris (pH 7.5). RNA and mRNA quantity was determined with a spectrophotometer (NanoDrop Technologies Inc., Wilmington, DE, USA). Total RNA and mRNA quality was assessed with an RNA Nano LabChip on an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).

cDNA synthesis, normalization and 454 pyrosequencing

Purified mRNA was used to construct a full length normalized cDNA pool through the services of Evrogen [50]. Briefly, the service utilized the SMART cDNA cloning methodology to generate a full length cDNA pool [51], which was normalized using a duplex-specific nuclease [52]. The double stranded normalized cDNA pool was sheared by nebulization and prepared for and sequenced as per manufacturer’s instructions on a Roche 454 GS FLX sequencer using an entire plate at Washington State University.

Data filtering and de novo assembly

35 Mb of sequence data representing 162,729 reads were generated by 454 sequencing. Quality trimming, adaptor sequence removal and size selection of reads was performed with Galaxy software ( [53]. After trimming adaptors, 128,720 reads with quality scores over 20 and sequence length longer than 50 bp were assembled with Abyss [54]. Parameters were adjusted for optimal assembly as measured by N50 statistic (a weighted median statistic such that 50% of the entire assembly is contained in contigs or scaffolds equal to or larger than this value).

Virus or viroid contamination detection

To determine whether any viruses or viroids were present in the fungi-infected plant cDNA sample, viroid and virus databases [55], including 41 complete viroid genomes and 2628 virus genomes, were downloaded from NCBI (released in April 2011). All EST contigs were analyzed with tBLASTx against viroid and virus databases. The e-value cutoff threshold was set at 1e-3. Contigs with a BLAST hit to viroid and virus databases were further analyzed by tBLASTx program against 3 legume genomes database and 7 fungi genomes database individually using the same cutoff threshold (see next section “development of a S. sclerotiorum and P. sativum parsing method”).

Development of a S. sclerotiorum and P. sativum parsing method

To separate S. sclerotiorum and pea ESTs from the mixed pool, a procedure based on that proposed by Hsiang et al. in 2003 [56] was employed with modifications (Figure 2). Briefly, the mixed ESTs were compared with tBLASTx (NCBI-BLAST-2.2.24+) to fungal and plant “proxy-reference” genome databases (Table 8). These proxy reference databases were established as the pea genome is not available and the inclusion of additional ascomycetes genomes to S. sclerotiorum (strain 1980) improved the assignment rate. The proxy-fungal genome database was a mixture of Sclerotinia sclerotiorum (strain 1980) and 6 closely related Ascomycete fungi (Botrytis cinerea, Chaetomium globosum, Fusarium graminearum, Magnaporthe grisea, Neurospora crassa and Verticillium dahlia) and a plant genome database including 3 sequenced legume genomes (Glycine max, Lotus japonicus and Medicago truncatula). ESTs that only matched to fungal or plant genome database with an e-value of 1e-03 or better were automatically classified into S. sclerotiorum or pea ESTs, respectively. ESTs, which matched (e-value <1e-3) to both fungi and plant databases, were further analyzed by comparing the e-value of best-hit from fungi and plant genome results. An e-value ratio was determined by dividing the best-hit e-value to fungi and plant genomes from the tBLASTx searches. A cutoff ratio were set at > =1e20 for pea ESTs, <=1e-20 for S. sclerotiorum ESTs and those that fell between 1e-20 and <1e20 were considered to be ambiguous. To acquire a final sort of results, those ESTs without a BLAST hit or those found to be ambiguous were assigned with BLASTn against known S. sclerotiorum or pea ESTs if their identity was above 95% in similarity across 95% of the sequence length. 81,449 pea ESTs (from flowers, leaves, cotyledons, epi- and hypocotyl, and etiolated and light treated etiolated seedlings) [7] and 57,751 S. sclerotiorum ESTs (from mycelia growing at neutral pH, developing apothecia and developing sclerotia---downloaded from BROAD database) were used to assist in the classification and annotation of contigs.

Table 8 Source of fungal plant genome databases used for tBLASTx EST assignment

To verify the feasibility of the EST parsing method, 17,533 S. sclerotiorum ESTs derived from developing S. sclerotiorum libraries were downloaded from BROAD institute and 18,547 P. sativum ESTs were obtained from the GenBank EST database by search keyword ‘Pisum sativum’. Vector contamination was removed from the downloaded ESTs by BLAST search with UniVec database (GenBank) in P. sativum and S. sclerotiorum ESTs were trimmed. After vector trimming, tBLASTx analysis of the downloaded ESTs was performed separately against the proxy-reference fungal and plant databases (Table 8). The following relevant data from tBLASTx output were extracted to an Excel file: query sequence name, query sequence length, fungi database target name, fungi database e-value for top match, total query sequence length for all match to fungi database, plant database target name, plant database top match e-value, total query sequence length for all match to plant database.

PCR to confirm validity of classified contigs

Fifty contigs from S. sclerotiorum and 50 contigs from pea were randomly sampled to check the validity of EST contig classification. Primers were designed for each contig using the program Primer3 [57]. cDNA from pea inoculated with S. sclerotiorum, cDNA from non-inoculated pea, cDNA from S. sclerotiorum growing on PDA medium, and genomic DNA extracted from pea and S. sclerotiorum using DNeasy plant mini kit (Qiagen, Valencia, CA, USA ) were used as template in PCR with primer pairs for each contig. PCR contained 4 μl of 5 × GoTaq PCR Buffer (Promega, Madison, WI, USA), 200 μM each dNTP, 2.5 μM each primer, 0.4 U of GoTaq polymerase, and approximately 50 ng of DNA template in a final volume of 20 μl. PCR were held at 94°C for 2 min; followed by 40 cycles of 94°C for 30 s, 60°C for 30 s, and 72°C for 1 min; with a final extension at 72°C for 10 min. PCR products from each contig were separated on a 1% agarose gel and visualized with ethidium bromide.

Gene annotation and analysis

The biological function of EST contigs was predicted with gene ontology (GO) terms based on BLASTx analysis using the program BLAST2GO [58, 59]. Default BLASTx parameters with an e-value threshold of 1e-3 and a high-scoring segment pairs (hsp) filter of 33 were retained so as to assign function to as many contigs as possible while ensuring short matching sequences less than 100 nucleotides were excluded. An annotation configuration with e-value-hit-filter 1.0E-6, annotation cut off “55” and GO weight “10” was selected.

Prediction of secretory/signal peptides for the S. sclerotiorum ESTs

The secretory/signal peptides for each S. sclerotiorum EST were analyzed using prediction algorithms. Firstly, OrfPredictor [60] was used to predict protein coding regions for the assembled ESTs. The output for OrfPredictor was a file of predicted coding regions from the ESTs in FASTA format, where the definition line contains the query identifier, the frame, the beginning and the end position of the predicted coding region, and the predicted protein peptide sequences. The peptide sequences were then used as a query for SignalP 3.0 ( with default settings, which predicts the presence and location of signal peptide cleavage sites [13]. Both a neural network (NN) and Hidden Markov Model (HMM) approach were used. EST contigs identified by both NN and HMM were considered to be secretory/signal peptides with high confidence.

Authors’ information

XF is a Research Associate at Michigan State University with a background and interest in expression profiling. KM is an Associate Professor and legume breeder in the Department of Plant Sciences at North Dakota State University. TC is Trait Production Manager at Dow AgroSciences LLC. TP is an Associate Professor of Plant Pathology at Washington State University, with particular interests in the mechanisms of fungal speciation and the genetics of host specificity. MC is an Assistant Professor of Plant Pathology at Michigan State University.


  1. Bolton MD, Thomma BPHJ, Nelson BD: Sclerotinia sclerotiorum (Lib.) de Bary: biology and molecular traits of a cosmopolitan pathogen. Mol Plant Pathol. 2006, 7 (1): 1-16. 10.1111/j.1364-3703.2005.00316.x.

    Article  CAS  PubMed  Google Scholar 

  2. Boland GJ, Hall R: Index of plant hosts of Sclerotinia sclerotiorum. Can J Plant Pathol. 1994, 16 (2): 93-108. 10.1080/07060669409500766.

    Article  Google Scholar 

  3. Porter LD, Hoheisel G, Coffman VA: Resistance of peas to Sclerotinia sclerotiorum in the Pisum core collection. Plant Pathol. 2009, 58 (1): 52-60. 10.1111/j.1365-3059.2008.01937.x.

    Article  Google Scholar 

  4. Cannon SB, Sterck L, Rombauts S, Sato S, Cheung F, Gouzy J, Wang X, Mudge J, Vasdewani J, Schiex T, et al: Legume genome evolution viewed through the Medicago truncatula and Lotus japonicus genomes. Proc Natl Acad Sci. 2006, 103 (40): 14959-14964. 10.1073/pnas.0603228103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Neumann P, Nouzova M, Macas J: Molecular and cytogenetic analysis of repetitive DNA in pea (Pisum sativum L.). Genome. 2001, 44 (4): 716-728.

    Article  CAS  PubMed  Google Scholar 

  6. Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63. 10.1038/nrg2484.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Franssen S, Shrestha R, Brautigam A, Bornberg-Bauer E, Weber A: Comprehensive transcriptome analysis of the highly complex Pisum sativum genome using next generation sequencing. BMC Genomics. 2011, 12 (1): 227-10.1186/1471-2164-12-227.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Baginsky S, Hennig L, Zimmermann P, Gruissem W: Gene expression analysis, proteomics, and network discovery. Plant Physiol. 2010, 152 (2): 402-410. 10.1104/pp.109.150433.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Emmersen J, Rudd S, Mewes HW, Tetko IV: Separation of sequences from host-pathogen interface using triplet nucleotide frequencies. Fungal Genet Biol. 2007, 44 (4): 231-241. 10.1016/j.fgb.2006.11.010.

    Article  CAS  PubMed  Google Scholar 

  10. Bowen JK, Mesarich CH, Rees-George J, Cui W, Fitzgerald A, Win J, Plummer KM, Templeton MD: Candidate effector gene identification in the ascomycete fungal phytopathogen Venturia inaequalis by expressed sequence tag analysis. Mol Plant Pathol. 2009, 10 (3): 431-448. 10.1111/j.1364-3703.2009.00543.x.

    Article  CAS  PubMed  Google Scholar 

  11. Fernandez D, Tisserant E, Talhinhas P, Azinheira H, Vieira A, Petitot AS, Loureiro A, Poulain J, DAS C, Silva MD, Duplessis S: 454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemileia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant-rust interaction. Mol Plant Pathol. 2012, 13 (1): 17-37. 10.1111/j.1364-3703.2011.00723.x.

    Article  CAS  PubMed  Google Scholar 

  12. Amselem J, Cuomo CA, van Kan JA, Viaud M, Benito EP, Couloux A, Coutinho PM, de Vries RP, Dyer PS, Fillinger S, et al: Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea. PLoS Genet. 2011, 7 (8): e1002230-10.1371/journal.pgen.1002230.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Zhao J, Buchwaldt L, Rimmer SR, Sharpe A, McGregor L, Bekkaoui D, Hegedus D: Patterns of differential gene expression in Brassica napus cultivars infected with Sclerotinia sclerotiorum. Mol Plant Pathol. 2009, 10 (5): 635-649. 10.1111/j.1364-3703.2009.00558.x.

    Article  CAS  PubMed  Google Scholar 

  14. Singh K, Foley RC, Onate-Sanchez L: Transcription factors in plant defense and stress responses. Curr Opin Plant Biol. 2002, 5 (5): 430-436. 10.1016/S1369-5266(02)00289-3.

    Article  CAS  PubMed  Google Scholar 

  15. Du H, Zhang L, Liu L, Tang XF, Yang WJ, Wu YM, Huang YB, Tang YX: Biochemical and molecular characterization of plant MYB transcription factor family. Biochemistry (Mosc). 2009, 74 (1): 1-11. 10.1134/S0006297909010015.

    Article  CAS  Google Scholar 

  16. Guo H, Ecker JR: The ethylene signaling pathway: new insights. Curr Opin Plant Biol. 2004, 7 (1): 40-49. 10.1016/j.pbi.2003.11.011.

    Article  CAS  PubMed  Google Scholar 

  17. Stepanova AN, Alonso JM: Ethylene signaling and response: where different regulatory modules meet. Curr Opin Plant Biol. 2009, 12 (5): 548-555. 10.1016/j.pbi.2009.07.009.

    Article  CAS  PubMed  Google Scholar 

  18. Xu X, Chen C, Fan B, Chen Z: Physical and functional interactions between pathogen-induced arabidopsis WRKY18, WRKY40, and WRKY60 transcription factors. Plant Cell Online. 2006, 18 (5): 1310-1326. 10.1105/tpc.105.037523.

    Article  CAS  Google Scholar 

  19. Kunkel BN, Brooks DM: Cross talk between signaling pathways in pathogen defense. Curr Opin Plant Biol. 2002, 5 (4): 325-331. 10.1016/S1369-5266(02)00275-3.

    Article  CAS  PubMed  Google Scholar 

  20. Perchepied L, Balagué C, Riou C, Claudel-Renard C, Rivière N, Grezes-Besset B, Roby D: Nitric oxide participates in the complex interplay of defense-related signaling pathways controlling disease resistance to Sclerotinia sclerotiorum in Arabidopsis thaliana. Mol Plant Microbe Interact. 2010, 23 (7): 846-860. 10.1094/MPMI-23-7-0846.

    Article  CAS  PubMed  Google Scholar 

  21. Wang Z, Tan X, Zhang Z, Gu S, Li G, Shi H: Defense to Sclerotinia sclerotiorum in oilseed rape is associated with the sequential activations of salicylic acid signaling and jasmonic acid signaling. Plant Sci. 2012, 184: 75-82.

    Article  CAS  PubMed  Google Scholar 

  22. Gruber S, Seidl-Seiboth V: Self versus non-self: fungal cell wall degradation in Trichoderma. Microbiology. 2012, 158 (Pt 1): 26-34.

    Article  CAS  PubMed  Google Scholar 

  23. Przybył K, Dahm H, Ciesielska A, Moliński K: Cellulolytic activity and virulence of Ophiostoma ulmi and O. novo-ulmi isolates. For Pathol. 2006, 36 (1): 58-67. 10.1111/j.1439-0329.2006.00429.x.

    Article  Google Scholar 

  24. Navarro-García F, Sánchez M, Nombela C, Pla J: Virulence genes in the pathogenic yeast Candida albicans. FEMS Microbiol Rev. 2001, 25 (2): 245-268.

    Article  PubMed  Google Scholar 

  25. Davidson AL, Chen J: ATP-Binding cassette transporters in bacteria. Annu Rev Biochem. 2004, 73 (1): 241-268. 10.1146/annurev.biochem.73.011303.073626.

    Article  CAS  PubMed  Google Scholar 

  26. Coleman JJ, Mylonakis E: Efflux in fungi: la Pièce de résistance. PLoS Pathog. 2009, 5 (6): e1000486-10.1371/journal.ppat.1000486.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Coleman JJ, White GJ, Rodriguez-Carres M, VanEtten HD: An ABC transporter and a cytochrome P450 of Nectria haematococca MPVI are virulence factors on Pea and are the major tolerance mechanisms to the Phytoalexin Pisatin. Mol Plant Microbe Interact. 2010, 24 (3): 368-376.

    Article  Google Scholar 

  28. Li L, Wright SJ, Krystofova S, Park G, Borkovich KA: Heterotrimeric G protein signaling in filamentous fungi*. Annu Rev Microbiol. 2007, 61 (1): 423-452. 10.1146/annurev.micro.61.080706.093432.

    Article  CAS  PubMed  Google Scholar 

  29. Fulcher AJ, Jans DA: Regulation of nucleocytoplasmic trafficking of viral proteins: an integral role in pathogenesis?. Biochim Biophys Acta. 2011, 1813 (12): 2176-2190. 10.1016/j.bbamcr.2011.03.019.

    Article  CAS  PubMed  Google Scholar 

  30. Lott K, Cingolani G: The importin β binding domain as a master regulator of nucleocytoplasmic transport. Biochim et Biophysica Acta (BBA) - Mol Cell Res. 2011, 1813 (9): 1578-1592. 10.1016/j.bbamcr.2010.10.012.

    Article  CAS  Google Scholar 

  31. Rep M: Small proteins of plant-pathogenic fungi secreted during host colonization. FEMS Microbiol Lett. 2005, 253 (1): 19-27. 10.1016/j.femsle.2005.09.014.

    Article  CAS  PubMed  Google Scholar 

  32. Soanes DM, Richards TA, Talbot NJ: Insights from sequencing fungal and oomycete genomes: what can we learn about plant disease and the evolution of pathogenicity?. Plant Cell. 2007, 19 (11): 3318-3326. 10.1105/tpc.107.056663.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Pancholi V, Fischetti VA: α-enolase, a novel strong plasmin(ogen) binding protein on the surface of pathogenic streptococci. J Biol Chem. 1998, 273 (23): 14503-14515. 10.1074/jbc.273.23.14503.

    Article  CAS  PubMed  Google Scholar 

  34. Bergmann S, Rohde M, Chhatwal GS, Hammerschmidt S: α-Enolase of Streptococcus pneumoniae is a plasmin(ogen)-binding protein displayed on the bacterial cell surface. Mol Microbiol. 2001, 40 (6): 1273-1287. 10.1046/j.1365-2958.2001.02448.x.

    Article  CAS  PubMed  Google Scholar 

  35. Antúnez K, Anido M, Arredondo D, Evans JD, Zunino P: Paenibacillus larvae enolase as a virulence factor in honeybee larvae infection. Vet Microbiol. 2011, 147 (1–2): 83-89.

    Article  PubMed  Google Scholar 

  36. Guo M, Guo W, Chen Y, Dong S, Zhang X, Zhang H, Song W, Wang W, Wang Q, Lv R, et al: The basic leucine zipper transcription factor Moatf1 mediates oxidative stress responses and is necessary for full virulence of the rice blast fungus Magnaporthe oryzae. Mol Plant Microbe Interact. 2010, 23 (8): 1053-1068. 10.1094/MPMI-23-8-1053.

    Article  CAS  PubMed  Google Scholar 

  37. Fischer G, Bang H, Ludwig B, Mann K, Hacker J: Mip protein of Legionella pneumophila exhibits peptidyl-prolyl-cis/trans isomerase (PPIase) activity. Mol Microbiol. 1992, 6 (10): 1375-1383. 10.1111/j.1365-2958.1992.tb00858.x.

    Article  CAS  PubMed  Google Scholar 

  38. Hacker J, Fischer G: Immunophilins: structure—function relationship and possible role in microbial pathogenicity. Mol Microbiol. 1993, 10 (3): 445-456. 10.1111/j.1365-2958.1993.tb00917.x.

    Article  CAS  PubMed  Google Scholar 

  39. Schmidt B, Tradler T, Rahfeld J-U, Ludwig B, Jain B, Mann K, Rücknagel KP, Janowski B, Schierhorn A, Küllertz G, et al: A cyclophilin-like peptidyl-prolyl cis/trans isomerase from Legionella pneumophila – characterization, molecular cloning and overexpression. Mol Microbiol. 1996, 21 (6): 1147-1160. 10.1046/j.1365-2958.1996.00061.x.

    Article  CAS  PubMed  Google Scholar 

  40. Yarden O, Yanofsky C: Chitin synthase 1 plays a major role in cell wall biogenesis in Neurospora crassa. Genes Dev. 1991, 5 (12B): 2420-2430. 10.1101/gad.5.12b.2420.

    Article  CAS  PubMed  Google Scholar 

  41. Soulie MC, Piffeteau A, Choquer M, Boccara M, Vidal-Cros A: Disruption of Botrytis cinerea class I chitin synthase gene Bcchs1 results in cell wall weakening and reduced virulence. Fungal Genet Biol. 2003, 40 (1): 38-46. 10.1016/S1087-1845(03)00065-3.

    Article  CAS  PubMed  Google Scholar 

  42. Besteiro S, Williams RA, Morrison LS, Coombs GH, Mottram JC: Endosome sorting and autophagy are essential for differentiation and virulence of Leishmania major. J Biol Chem. 2006, 281 (16): 11384-11396. 10.1074/jbc.M512307200.

    Article  CAS  PubMed  Google Scholar 

  43. Liu XH, Lu JP, Zhang L, Dong B, Min H, Lin FC: Involvement of a Magnaporthe grisea serine/threonine kinase gene, MgATG1, in appressorium turgor and pathogenesis. Eukaryot Cell. 2007, 6 (6): 997-1005. 10.1128/EC.00011-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Hu G, Hacham M, Waterman SR, Panepinto J, Shin S, Liu X, Gibbons J, Valyi-Nagy T, Obara K, Jaffe HA, et al: PI3K signaling of autophagy is required for starvation tolerance and virulence of Cryptococcus neoformans. J Clin Invest. 2008, 118 (3): 1186-1197. 10.1172/JCI32053.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Nguyen QB, Itoh K, Van Vu B, Tosa Y, Nakayashiki H: Simultaneous silencing of endo-β-1,4 xylanase genes reveals their roles in the virulence of Magnaporthe oryzae. Mol Microbiol. 2011, 81 (4): 1008-1019. 10.1111/j.1365-2958.2011.07746.x.

    Article  CAS  PubMed  Google Scholar 

  46. Jolie RP, Duvetter T, Van Loey AM, Hendrickx ME: Pectin methylesterase and its proteinaceous inhibitor: a review. Carbohydr Res. 2010, 345 (18): 2583-2595. 10.1016/j.carres.2010.10.002.

    Article  CAS  PubMed  Google Scholar 

  47. Markaryan A, Zaborina O, Punj V, Chakrabarty AM: Adenylate kinase as a virulence factor of Pseudomonas aeruginosa. J Bacteriol. 2001, 183 (11): 3345-3352. 10.1128/JB.183.11.3345-3352.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Missall TA, Pusateri ME, Donlin MJ, Chambers KT, Corbett JA, Lodge JK: Posttranslational, translational, and transcriptional responses to nitric oxide stress in Cryptococcus neoformans: implications for virulence. Eukaryot Cell. 2006, 5 (3): 518-529. 10.1128/EC.5.3.518-529.2006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Chaves GM, Bates S, Maccallum DM, Odds FC: Candida albicans GRX2, encoding a putative glutaredoxin, is required for virulence in a murine model. Genet Mol Res. 2007, 6 (4): 1051-1063.

    CAS  PubMed  Google Scholar 

  50. Evrogen:,

  51. Zhu YY, Machleder EM, Chenchik A, Li R, Siebert PD: Reverse transcriptase template switching: a SMART (TM) approach for full-length cDNA library construction. Biotechniques. 2001, 30 (4): 892-897.

    CAS  PubMed  Google Scholar 

  52. Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, Kozhemyako VB, Matz MV, Meleshkevitch E, Moroz LL, Lukyanov SA, Shagin DA: Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucleic Acids Res. 2004, 32 (3): E37-10.1093/nar/gnh031.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Goecks J, Nekrutenko A, Taylor J: Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010, 11 (8): R86-10.1186/gb-2010-11-8-r86.

    Article  PubMed Central  PubMed  Google Scholar 

  54. Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJ, Birol I: ABySS: a parallel assembler for short read sequence data. Genome Res. 2009, 19 (6): 1117-1123. 10.1101/gr.089532.108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. NCBI VGi:,

  56. Hsiang T, Goodwin PH: Distinguishing plant and fungal sequences in ESTs from infected plant tissues. J Microbiol Methods. 2003, 54 (3): 339-351. 10.1016/S0167-7012(03)00067-8.

    Article  CAS  PubMed  Google Scholar 

  57. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    CAS  PubMed  Google Scholar 

  58. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.

    Article  CAS  PubMed  Google Scholar 

  59. Conesa A, Gotz S: Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008, 2008: 619832-

    Article  PubMed Central  PubMed  Google Scholar 

  60. Min XJ, Butler G, Storms R, Tsang A: OrfPredictor: predicting protein-coding regions in EST-derived sequences. Nucleic Acids Res. 2005, 33 (Web Server issue): W677-W680.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  61. Gastebois A, Clavaud C, Aimanianda V, Latgé J-P: Aspergillus fumigatus: cell wall polysaccharides, their biosynthesis and organization. Future Microbiol. 2009, 4 (5): 583-595. 10.2217/fmb.09.29.

    Article  CAS  PubMed  Google Scholar 

Download references


This project was supported by USDA-ARS Specific Cooperative Agreement #58-5442-9-239 (National Sclerotinia Initiative). We thank Ahmit Dhingra, Department of Horticulture, Washington State University for 454 pyrosequencing.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Martin I Chilvers.

Additional information

Competing interests

The authors declare that they have no competing interest.

Authors’ contributions

XF performed bioinformatics analysis, validated the EST parsing and drafted the manuscript. KM provided pea lines for analysis and contributed to direction of the study. TC participated in design of the study and conducted signal peptide analysis. TP initiated the project and directed the study. MC directed the project, performed inoculations and initial experiments to develop the normalized cDNA and drafted the manuscript. All authors read, edited and approved the final manuscript.

Electronic supplementary material

Additional file 1: Fasta files for 10158 contigs (pea and Sclerotinia ESTs), parsed into ambiguous contigs (1a), no blast-hit contigs (1b), pea contigs (1c) and S. sclerotiorum contigs (1d). (ZIP 628 KB)

Additional file 2: S. sclerotiorum EST contigs encoding potential secretory/signal peptides. (DOCX 40 KB)


Additional file 3: Comparison of BLASTn and tBLASTx method to assign artificial EST mixture of pea- S. sclerotiorum (including 18,490 pea ESTs and 17,198  S. sclerotiorum ESTs). (DOCX 17 KB)


Additional file 4: The assignment results of an artificial EST mixture using the tBlastx method against 7 fungal genome databases and the S. sclerotiorum genome only. (DOCX 17 KB)


Additional file 5: Comparison of different e-value ratios (fungi/plant) to distinguish species ESTs from the artificial EST mixture. (DOCX 18 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access 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.

Reprints and Permissions

About this article

Cite this article

Zhuang, X., McPhee, K.E., Coram, T.E. et al. Rapid transcriptome characterization and parsing of sequences in a non-model host-pathogen interaction; pea-Sclerotinia sclerotiorum. BMC Genomics 13, 668 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: