RNA-Seq Analysis of the Sclerotinia homoeocarpa – Creeping Bentgrass Pathosystem

Sclerotinia homoeocarpa causes dollar spot disease, the predominate disease on highly-maintained turfgrass. Currently, there are major gaps in our understanding of the molecular interactions between S. homoeocarpa and creeping bentgrass. In this study, 454 sequencing technology was used in the de novo assembly of S. homoeocarpa and creeping bentgrass transcriptomes. Transcript sequence data obtained using Illumina's first generation sequencing-by-synthesis (SBS) were mapped to the transcriptome assemblies to estimate transcript representation in different SBS libraries. SBS libraries included a S. homoeocarpa culture control, a creeping bentgrass uninoculated control, and a library for creeping bentgrass inoculated with S. homoeocarpa and incubated for 96 h. A Fisher's exact test was performed to determine transcripts that were significantly different during creeping bentgrass infection with S. homoeocarpa. Fungal transcripts of interest included glycosyl hydrolases, proteases, and ABC transporters. Of particular interest were the large number of glycosyl hydrolase transcripts that target a wide range of plant cell wall compounds, corroborating the suggested wide host range and saprophytic abilities of S. homoeocarpa. Several of the multidrug resistance ABC transporters may be important for resistance to both fungicides and plant defense compounds. Creeping bentgrass transcripts of interest included germins, ubiquitin transcripts involved in proteasome degradation, and cinnamoyl reductase, which is involved in lignin production. This analysis provides an extensive overview of the S. homoeocarpa-turfgrass pathosystem and provides a starting point for the characterization of potential virulence factors and host defense responses. In particular, determination of important host defense responses may assist in the development of highly resistant creeping bentgrass varieties.


Introduction
Sclerotinia homoeocarpa (Bennett) is an ascomycete fungus that causes dollar spot disease on turfgrass world-wide [1]. S. homoeocarpa can affect all species of turfgrass as well as some dicot plants [1]. Creeping bentgrass (Agrostis stolonifera L.) is a cool-season turfgrass that is common on golf course greens in the United States and Canada [2]. Many commonly used creeping bentgrass cultivars are highly susceptible to dollar spot disease and are subject to frequent, low mowing practices that promote disease outbreaks [1,3].
Symptoms of dollar spot include straw-colored, hourglassshaped lesions with characteristic reddish-brown borders. Diseased areas grow to about 2.5 cm wide [1]. Under high disease pressure and favorable conditions the diseased areas will merge to form larger patches of diseased turf. Environmental conditions favoring disease development have been well documented [4,5,6]; however, prediction models designed to reduce fungicide inputs have been ineffective due to a lack of understanding of the S. homoeocarpa lifestyle, epidemiology, and disease etiology [1,4].
S. homoeocarpa is thought to overwinter in the thatch of turf swards as stroma. In the spring, S. homoeocarpa becomes active and infects the newly emerging leaf tissue through wounds, stomates, and directly with appressorium formation [7,8,9]. Only infertile apothecia have been recorded for North American isolates [10]; however, population studies suggest that genetic recombination in this fungus is possible [11]. Therefore, the possibility of ascospore production and the role of sexual or asexual spores as initial inoculum for dispersal of S. homoeocarpa cannot be discounted. Early studies noted root browning and cell death through the production of diffusible toxins [12]. Recent studies have identified several tetranorditerpenoid compounds that could be responsible for the root-browning and exhibit extremely phytotoxic properties; however, a correlation between the production of these compounds and disease symptoms was not established [13].
Cultural management strategies for dollar spot include maintaining adequate nitrogen balance, promotion of air flow to assist in dew removal, and using moderately resistant cultivars or species of turf [14,15]. On highly maintained areas such as golf courses, cultural practices are not sufficient for management of dollar spot, and fungicides are often applied biweekly to weekly. High amounts of fungicide use have resulted in resistance to several chemical classes commonly used on turf [1]. Other nonfungicide products that have been marketed for dollar spot control include plant defense activators that work by activating two different plant defense pathways: induced systemic resistance (ISR) and systemic acquired resistance (SAR). However, it is unclear which pathway would be most beneficial for preventing dollar spot epidemics.
A better understanding of the molecular interactions between S. homoeocarpa and creeping bentgrass will be essential for the development of more sustainable and practical management strategies, including the use of plant defense activators and the development of cultivars with increased resistance to S. homoeocarpa.
The introduction of next generation sequencing, also termed massively parallel sequencing (MPS), has enabled researchers to sequence the genomes and transcriptomes of organisms at a relatively low cost in return for a vast amount of data with quantitative properties [16,17,18,19]. In this paper, two MPS technologies were used to generate sequence data for RNA-Sequence (RNA-Seq) analysis: Illumina's sequencing-by-synthesis (SBS) and Roche's 454-pyrosequencing. The 454 reads were used for the de novo assembly of S. homoeocarpa and creeping bentgrass transcriptome libraries. SBS reads were mapped to the 454 assemblies to calculate transcript levels from S. homoeocarpa and creeping bentgrass during dollar spot disease development. The objective of this study was to identify transcripts that may be important for fungal virulence and creeping bentgrass defense. The results of the analysis will be used to form testable hypotheses for future studies on dollar spot etiology and turfgrass defense mechanisms.

Transcriptomic Analysis
The S. homoeocarpa (SH) 454 sequencing data contained 600,760 reads total, and the A. stolonifera (AS) 454 sequence data contained 205,403 reads total ( Table 1). The 454 read lengths ranged from 50 bp to .500 bp, with a majority of the reads between 400 and 500 bp in length. The transcriptome coverage for each of the 454 assemblies was 3.36 coverage for AS and 17.26 coverage for the SH assembly. These sequencing coverages were calculated by dividing the total number of sequence reads by the size of the respective, assembled transcriptome libraries.
The SBS data, which was used for calculating significant differences in transcript levels between libraries, resulted in 4.3-7.2 million reads (Table 2); however, these reads were only 16 bp long. This is due to the use of first generation Illumina sequencing protocols. The SBS reads provided ample coverage of the entire transcript assemblies: 146 coverage for SH and 126 for AS. The SBS library construction resulted in 9, 319 SH transcripts with lengths ranging from 400 bp to 3, 500 bp and a distribution peak at 500 bp. Construction of the SBS transcript library for AS resulted in 20, 293 transcript sequences with lengths ranging from 200 bp to 8,500 bp and distribution peak at 480 bp. The full length of transcripts was not estimated since very few annotated coding sequences from either SH or AS were available at Genbank. The mapping percentage of SBS transcripts was quite low. Only 67.4% of SH PDA reads, 63.8% of SH PDB reads, and 33.3% of Interaction reads were mapped to the SH transcripts assembly, respectively (Table 2). Only 16.8% of AS only reads, and 5.7% of Interaction reads mapped to the AS transcripts library (Table 2). This is likely due to the short length of the SBS reads and stringent mapping parameters used in this experiment, which eliminated reads mapping to more than one location and did not allow for any base mismatches.
A Venn diagram was constructed to identify the distribution of the unique and common reads in both SBS and 454 SH, AS, and interaction libraries (Figure 1). There are discrepancies in the read distribution that appear between the SBS and the 454 libraries. The most obvious example is the increase in the number of reads that are unique to the SH-only library from 6.2% in the SBS to 58.9% in the 454 data. This is because there were more total SH reads in the 454 dataset, where four of the six libraries (69.8% of the reads) were constructed using RNA from S. homoeocarpa grown in pure culture. In the SBS dataset, only two out of four (31.7% of the reads) of the libraries were constructed from S. homoeocarpa grown in pure culture. The number of reads common to the interaction and control libraries was consistently low in both SBS and 454 datasets ( Figure 1). The common reads are likely a result of similar sequences in both organisms, including those of transposable elements found in both organisms.

Comparative analysis
The transcripts were blasted against the NCBI protein database for functional annotation using Blast2GO [20] blastx function. Fifty seven percent (5,257) of the 9,319 SH transcripts returned at least one annotation hit. Likewise, 57% (11, 551) of the 20, 293 AS transcripts returned at least one hit. From the blastn taxonomy report, top hits of the S. homoeocarpa sequences were primarily from Sclerotinia sclerotiorum (Lib.) de Bary (44. 1%) and Botryotinia fuckeliana (de Bary) Whetzel (34.9%) (Figure 2 a). This is not surprising considering their taxonomic similarity to S. homoeocarpa, and the abundance of S. sclerotiorum and B. fuckeliana sequence data available on GenBank. The 26 fungal species represented in the taxonomy report belong to various classes of the subphylum Pezizomycotina, with seven of the 26 species being pathogens of monocot plants. There were a small number of sequences related to monocots including rice (Oryza sativa Linnaeus, 1.6%), bamboo (Phyllostachys edulis (Carr.) Houz., 0.3%), and barley (Hordeum vulgare Linnaeus, 0.4%). We attributed this either to contamination of the fungal sequences by A. stolonifera sequences in the interaction library, to the presence of very highly conserved sequences, or to the presence of sequences occurring in both organisms, such as those of transposable elements. The top blast hit species for the A. stolonifera sequences included predominantly monocot grass species (Figure 2 b). Fungal sequences that were found included those of Pyrenophora teres Smed.-Pet. and Phaeosphaeria nodorum (E. Müll) Hedjar, pathogens of barley and wheat, respectively.

Comparison of RPKM ratios to Relative Transcription ratios
Reads Per Kilobase of exon model per Million mapped reads (RPKM) value was calculated of all SH and AS transcripts in different libraries and the Log 2 ratio was calculated to reflect the fold changes of interaction library compared to control libraries (SH PDA and AS only). Relative transcription data from real-time RT-PCR assays was obtained for selected transcripts with significantly different transcription levels during infection. The relative transcription data verified that the RPKM calculating and statistical analysis of transcriptome data was accurate ( Figure 3). The real time RT-PCR data confirms the upregulation of two fungal xylanases (SH_5411 and SH_5726), a fungal laccase (SH_7961), and a creeping bentgrass germin protein (AS_608) ( Figure 2). Downregulated transcripts SH_6925 (scyatalone dehydrogenase) and SH_8369 (b-mannosidase) were also verified ( Figure 2).

Transcripts Expression Analysis
To identify transcripts that were up-regulated in the interaction condition, transcripts with statistically significant changes of RPKM values during infection were examined to determine if there were patterns in the types of genes with increased or decreased transcription level. Table 3 summarizes some of the types of SH transcripts of interest during infection including a long list of glycosyl hydrolase transcripts and various proteases. Table 4 summarizes transcripts of interest in the AS library including a long list of retrotransposons, defense-related proteins including 13 different germin-like proteins, and various enzymes involved in plant hormone synthesis and regulation. Interestingly, all jasmonate-induced protein transcripts were significantly down regulated at 96 hpi. A full list of SH and AS transcripts that are significantly over or under expressed during infection are available as supplementary data (Tables S1 and S2, respectively).

Gene Ontology Analysis
The list of significantly upregulated SH and AS transcripts were sorted into gene ontology (GO) term categories for molecular function, biological processes, and cellular component using the Blast2GO [20] program (Figures 4 and 5). The molecular function of GO terms associated with upregulated SH transcripts include a large percentage of plant cell-wall degrading enzymes categories (Figure 4 a) including polygalacturonase activity (4.0%), cellulase activity (5.1%), endopeptidase activity (4.0%), a-arabinofuranosidase activity (2.9%), and serine-type peptidase activity (5.1%). Other upregulated molecular function of GO terms included those for xenobiotic-transporting activity (5.1%) and oxidoreductase

Conserved Domain Analysis
Protein domains were identified in the S. homoeocarpa and A. stolonifera transcripts libraries using HMMER software (v3.0) [21]. Also among SH transcripts, 1588 were predicted to contain signal peptide by signalP (v4.0) [22], and thus could be candidates for secreted proteins. Three proteases with secretion signals included two subtilisin/sedolisin proteases with log-fold increases of 6.4 and 6.2 as well as a cuticle degrading serine protease with a log-fold increase of 4.0. Conserved domains significantly enriched in the interaction library (P#0.01) are listed in Table 5 (SH conserved domains) and Table 6 (AS conserved domains). Domains enriched in the SH interaction library include a variety of glycosyl hydrolases, proteases domains, and transporters (Table 5). Domains of interest in the enriched AS interaction library include cytochrome P450, various ABC transporters, and cupin domains, which include germin-class enzymes ( Table 6). A complete list of SH and AS conserved domains can be found as supplementary data in tables S3 and S4, respectively.

Discussion
Previous analyses using RNA-Seq have evaluated the reproducibility and accuracy of SBS [18] and 454 data for RNA transcription studies [23]. Both technologies have demonstrated higher-quality results than microarray and EST data [18,19]. However, SBS technologies were found to be superior due to the greater coverage and depth of sequences for a much lower cost. In this study, the mapping percentage of the SBS data to the S. homoeocarpa assembly ranged from 33.3-67.4% and A. stolonifera mapping percentages ranged from 5.7-16.8%. The mapping percentages were low due to the very stringent parameters that eliminated all reads with even a single mismatch and reads that mapped to more than a single location, as well as possible transcriptome misassembly. The use of these parameters and their resulting influence on the statistical analysis of the RPKM data was verified with real-time relative expression data. Mapping uncertainty in RNA-Seq analyses is one of the pitfalls of the SBS sequencing method. Short sequence lengths, as was the case in this study, make it difficult to account for paralogous genes, alternatively spliced isoforms, and low complexity sequences, all of which result in multireads that are eliminated in the mapping process [24]. Low mapping percentage may also be due to the lower coverage of 454 transcriptome assemblies compared to the high SBS coverage.
RNA-Seq analysis in this study has made it possible to identify potential pathogenicity factors transcribed by S. homoeocarpa during infection and colonization of A. stolonifera. The most striking discovery was the number of enriched S. homoeocarpa transcripts within the interaction library that encoded glycosyl hydrolase enzymes. In fact, 52 of the upregulated SH transcripts encoded glycosyl hydrolase genes and 22.3% of molecular function GO terms were associated with glycosyl hydrolase activity. Within the total SH transcript library, a keyword search using search terms such as glycosyl hydrolase, cellulase and pectinase revealed over 100 transcripts annotated as glycosyl hydrolases in the SH library, representing 85 different families of the total 125 glycosyl hydrolase families [25,26]. The numerous transcripts representing such a varied group of glycosyl hydrolase families corroborate the ability of this fungus to infect and cause disease on a wide range of plant hosts. In this study, S. homoeocarpa glycosyl hydrolase genes showing upregulation were predominantly xylanases and arabi- nases. This is not surprising since grass cell walls consist of approximately 40% xylans as well as glucuronoarabinoxylans (GAX) that make up a majority of monocot hemicelluloses [27]. Enzymes such as pectinases, xylanases, and cellulases are often virulence and pathogenicity factors for plant pathogenic fungi and bacteria. For example, Botrytis cinerea Pers. and Septoria nodorum (Berk) Berk., demonstrate reduced virulence when a single gene encoding a cell wall degrading enzyme is knocked-out [28,29,30]. However, the large number of upregulated and functionally redundant cell wall degrading enzymes expressed by S. homoeocarpa interacting with creeping bentgrass will make verifying the role of any one enzyme extremely difficult. Numerous proteinases, mostly serine proteases, were also enriched in the fungal interaction library. Secretion signals were detected on three of the proteases of interest, including two subtilisin/sedolisin proteases and one identified as a cuticle degrading serine protease. These three secreted proteases were upregulated in the interactoin library by 4.0-6.4 log-fold. Of the molecular function GO terms, 5.1% were associated specifically with serine-protease activity. Serine proteases are pathogenicity determinants for a variety of pathogenic fungi including nematode-parasitizing fungi, entomopathogens, and plant pathogenic fungi [31]. For example, protease deficient mutants of the hemibiotroph Pyrenopeziza brassicae B. Sutton & Rawl. were either avirulent or had dramatically reduced virulence [32]. Proteases in other plant pathogenic fungi have been characterized, but their role in pathogenicity has not necessarily been confirmed in each case [33,34].
In general, the presence of a variety of hydrolytic enzymes including glycosyl hydrolases and proteinases supports the broad host range reported for S. homoeocarpa [1]. The diversity and number of enzymes also support saprophytic ability of S. homoeocarpa [35]. For example, fungi with superior saprotrophic abilities, such as Aspergillus spp. were found to produce a more diverse collection of hydrolytic enzymes under a variety of conditions, while phytopathogens with inferior saprophytic abilities such as Verticillium albo-atrum Reinke & Berthold and V. lecanii (Zimm.) Viégas, produced a more targeted array of enzymes suited to their particular host range [35]. The large number of hydrolytic enzymes expressed by S. homoeocarpa supports the hypothesis of superior saprophytic abilities of the dollar spot pathogen, which may play a role in the overwintering of the fungus within the thatch during the winter months.
Multi-drug resistance (MDR) ATP binding cassette (ABC) transporter proteins were also upregulated in the S. homoeocarpa interaction library, with seven of the MDR-ABC transporters being upregulated while only one MDR-ABC transporter transcript was downregulated. Members of the MDR-ABC transporter family have important roles in fungicide resistance, resistance to plant defense compounds and efflux of endogeneous toxins [36,37]. In B. cinerea, both ABC and MFS transporters have been implicated in resistance of the pathogen to a variety of fungicides including phenylpyrroles, azoles, anilopyrimidines, dicarboximides, and strobilurins [38,39,40]. The ABC transporter BcatrB of B. cinerea is strongly induced in the presence of fludioxonil, cyprodinil, cyproconazole, tebuconazole, and trifloxystrobin suggesting a role for this single transporter in resistance to multiple fungicide types [41]. Enrichment of a variety of transport proteins with similar function in the S. homoeocarpa interaction library may provide an explanation for the development of isolates that are resistant to more than one fungicide class. For example, isolates are known that are resistant to demethylation inhibitors as well as benzimidazoles and dicarboximides [1,42,43]. Aside from fungicide resistance, ABC transporters are also implicated in resistance to plant defense compounds and secretion of fungal toxins. For example, BcatrB not only provides resistance to fungicides but also to the plant defense compounds resveratrol and camalexin [39,40]. Characterization of these transporter proteins could lead to a better understanding of fungicide resistance mechanisms and S. homoeocarpa-creeping bentgrass interactions, resulting in more effective fungicide application programs and fungicide resistance management strategies.
Integrated fungicide resistance management recommendations often include using cultivars of grass that have increased pathogen resistance. Several cultivars of creeping bentgrass have decreased susceptibility to dollar spot disease [44,45]; however, the molecular defense mechanisms responsible for the decrease in susceptibility are unknown. In this study, some of the most prominent, upregulated defense-related transcripts were the germin and germin-like proteins. Germins belong to the cupin superfamily and include the important defense enzyme, oxalate oxidase [46].
Germin-type oxalate oxidases are enzymes that are unique to Gramineae plants [46]. In creeping bentgrass cultivars, oxalate oxidase has been shown to be more active in the moderately resistant creeping bentgrass cultivar, L-93, compared to Crenshaw, a highly susceptible cultivar [47]. Oxalic acid is a pathoge- nicity factor for Sclerotinia sclerotiorum (Lib.) de Bary [48], and oxalic acid is produced in large quantities by S. homoeocarpa as well [49]. Cloning oxalate oxidase into dicot hosts such as sunflower, soybean, canola, and tomato renders them resistant to S. sclerotiorum [50,51,52]. It is possible that increased production of oxalate oxidases by some cultivars of creeping bentgrass is responsible for their decreased susceptibility to S. homoeocarpa. Furthermore, the product of oxalate oxidase activity is H 2 O 2 , a plant defense signal that is capable of passing through cell membranes, is directly toxic to invading microorganisms, cross- links glycoproteins and phenolic molecules to form defense structures, and can induce systemic acquired resistance [53,54]. Clearly, more research on germin-like proteins and oxalate oxidase activity in creeping bentgrass is warranted. In addition to the plant defense-related transcripts, other notable creeping bentgrass transcripts were identified as transposons and retrotransposons including pong, em, mutator, gypsy, and copia subclasses. This is not surprising since it has been documented that transposons make up approximately 80% of the DNA in cereal crops [55]. In a recent EST analysis of creeping bentgrass and colonial bentgrass, 1.4% and 0.18% of the total ESTs were retrotransposon related, and retrotransposon transcripts were eight times higher in creeping bentgrass under disease stress than in healthy colonial bentgrass plants that were not under disease stress [56]. In the RNA-Seq data analysed here, 878 of 20,493 (4%) of the total AS transcripts were identified as transposons and 731 (3.5%) were identified as retrotransposons. Seventy-two of the 1017 transcripts upregulated in the interaction library were identified as retrotransposons: that is 7% of all upregulated AS transcripts. This is in clear support of the previous study concluding that retrotransposons in creeping bentgrass are likely activated by disease stress [56]. In the shared library (reads present in both AS and SH), only 0.1% of the reads mapped to transcripts identified as transposons. This means that it is not likely that identical transposon sequences are shared between SH and AS. However, several of the classes of retrotransposons identified in the study, such as the Ty and copia classes, are known to exist in both fungi and plants [57,58]. The upregulation of 72 retrotransposons in this study provides direct support for activation of transposons in creeping bentgrass in response to disease stress as demonstrated previously [56]. It has been reported that retrotransposons are transcriptionally activated in plants that are undergoing stress such as during protoplast formation, wounding, and in the presence of salicylic acid and jasmonate [59]. The influence of retrotransposons, if any, on plant defense pathways has not yet been determined; however, their increased transcription during infection should be investigated. For example, it would be interesting to determine if any of these retroelements are located within or near genes involved in plant defense responses, where they could affect the expression of these genes.
The RNA-Seq analysis presented in this study demonstrates the immense value of MPS technology for developing a comprehensive understanding of disease pathosystems. Advances in SBS technology since the start of this study would now provide adequate read length and depth for the de novo transcriptome assemblies, making it even more powerful and cost efficient. The resulting data analysis will be used to support future research efforts to characterize S. homoeocarpa virulence factors and creeping bentgrass defense-related transcripts.

Preparation of RNA for SBS Libraries
Sclerotinia homoeocarpa isolate MB01 was isolated from creeping bentgrass at The Ohio State University Turfgrass Research and Education Facility in Columbus, Ohio. The fungus was cultured on potato dextrose agar (PDA, Difco, Franklin Lakes, NJ) and 5 mm plugs of actively growing mycelium were transferred to potato dextrose broth (PDB, Difco) medium and incubated at 26uC with shaking at 160 rpm for 96 h. Cultures grown on PDA overlayed with cellophane membranes were also used for SBS library construction. At 96 h, mycelia were filtered from PDB or scraped from the cellophane membrane on PDA, and total RNA was extracted using TriReagent (Sigma, St. Louis, MO) according to manufacturer's instructions. For interaction and creeping bentgrass libraries, creeping bentgrass cv. Crenshaw was grown in a growth chamber at 26uC day and 22uC night temperatures with a 12 h photoperiod and at 70% relative humidity. Threeweek old seedlings were challenged with millet seeds colonized by S. homoeocarpa MB01. These plants were incubated in clear bags to create humidity for 48 h after which the plants were uncovered during the day and bagged again for the night. After 96 h, leaves from inoculated and noninoculated plants were harvested, homogenized in liquid nitrogen, and RNA was extracted using TriReagent (Sigma) according to the manufacturer's instructions.
The SBS library templates were prepared using the Illumina Duplex-Specific thermostable nuclease (DSN) normalization kit and were analyzed using the Illumina Genome Analyzer I (GAI) at the Molecular

Sequencing, Assembly, and Library Construction
The 454 raw sequence reads were assembled using GS Data Analysis Software (v2.5, Roche, Indianapolis, IN) after removal of adaptor sequences. S. homoeocarpa (SH) reads from the PDA, PDB, and 96 hpi Interaction libraries, and the A. stolonifera (AS) reads from the 96 h control and 96 h Interaction libraries were assembled using GS de novo assembler (v2.5, Roche) with parameter setting as 90% identity and a minimum 40 bp overlap. To eliminate SH reads in AS assembly library, reads from the S. homoeocarpa (SH) Interaction libraries that matched creeping bentgrass (Agrostis stolonifera) GenBank EST data were removed from the SH transcript assembly. Then reads from the A. stolonifera (AS) Interaction library matching SH assembly were also removed. Both S. homoeocarpa and A. stolonifera transcriptome assembly as used in this study was submitted to GenBank TSA database (Accession: PRJNA84359).
Full length SBS reads were directly mapped to the SH and AS transcriptome libraries using the Maq (v 0.7) [60] alignment. Mapping parameters allowed for only unique mapped reads with no mismatches. The number of mapped SBS reads for each transcript was counted by ''soap.coverage'' in SOAPAligner package.

Transcript Expression Analysis
RPKM values were calculated for each transcript by dividing number of mapped reads by length of the transcripts and number of total sequenced reads in this library. Fisher's Exact Test, a test used for analyzing unreplicated transcript data [61] was used to determine statistical significance of transcript RPKM value change. The test was done using a pairwise comparison of Interaction RPKM values versus RPKM from fungus grown on PDA, or for Interaction RPKM values versus uninoculated control grass RPKM values. The fungal data contained 9,319 unique annotated transcripts and the grass data contained 20,293 unique annotated transcripts. Therefore, the statistical significance of the Fisher's Test was evaluated against a Bonferroni corrected P-value of 1.86610 25 and 2.43610 26 for fungal and grass transcripts, respectively. Furthermore, only transcripts with at least 2-log fold change in transcript abundance were selected as upregulated. These conservative criteria were applied to avoid false positives.
The log-fold RPKM values of selected transcripts were validated using real time reverse transcription polymerase chain (RT-PCR) reaction on an iQ5 real time thermocycler (Biorad, Hercules, CA). RNA was extracted from fungus grown on PDA, noninoculated grass, and infected grass according to the previously described methods. RNA was extracted using Trizol (Invitrogen, Carlsbad, CA) and RNA extracts were treated with RQ1 DNAse (Promega, St. Louis Obispo, CA) according to the manufacturer's directions, except that RNase Inhibitor (Promega) was added to each 100 ml of RNA extract. Complementary DNA was created from 1 mg of total RNA using iScript (BioRad, Hercules, CA). PCR was conducted using the SYBR supermix (Biorad) as per manufacturer's directions. The cycle conditions were 95uC for 3 minutes followed by 40 cycles of 95uC for 30 s, 58.5uC for 30 s, and 72uC for 30 s. The melting temperature profiles and gel electrophoresis was used to evaluate the specificity of the reactions and the absence of primer dimers. Real time RT-PCR was conducted with three biological replicates for each of grass, interaction, and fungal cDNA samples. Three technical replicates of each cDNA sample were used in each experiment.
Reaction efficiency and relative expression data were analyzed using the relative expression software tool (REST) program (Qiagen; [62]). Log base two values of relative expression ratios calculated in the REST program and the corresponding log base two RPKM expression ratios were compared graphically.

Library annotation
Taxonomic and functional annotation of the SH and AS transcripts were conducted by using Blast2GO [20] software to run blastx and blastn algorithms against non-redundant nucleotide/protein database from the National Center for Biotechnology Information (National Institutes of Health). Within the SH library, any transcripts that resulted in top blastn hit species of Hordeum, Oryza, Sorghum, Agrostis, or Vitis were removed. Similarly, any transcripts in the AS library that had top blastn hit species of Sclerotinia, Botryotinia, Glomerella, Ajellomyces, or Nectria were removed. Combined graphs for GO terms associated with statistically upregulated transcripts are presented in the results section.

Conserved Domain Analysis
Transcript sequences were translated to proteins in all 6 possible frames, and conserved protein domains were identified using HMMER (v3.0) [21] to screen the Pfam-A database. All the predicted domains with E-value ,0.01 were reported. There were a total of 14,161 domains identified in the SH transcripts assembly and 18,591 domains identified in the AS transcripts assembly. The enrichment ratio and significant value was calculated based on a Fisher's Exact Test. Descriptions for enriched domains were found on the Pfam website.