Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Transcriptional Reprogramming of Wheat and the Hemibiotrophic Pathogen Septoria tritici during Two Phases of the Compatible Interaction

  • Fen Yang ,

    Contributed equally to this work with: Fen Yang, Wanshun Li

    yangf@plen.ku.dk

    Affiliation Department of Plant and Environmental Sciences, Faculty of Science, University of Copenhagen, Frederiksberg, Denmark

  • Wanshun Li ,

    Contributed equally to this work with: Fen Yang, Wanshun Li

    Affiliation BGI-Tech, BGI, Shenzhen, China

  • Hans J. L. Jørgensen

    Affiliation Department of Plant and Environmental Sciences, Faculty of Science, University of Copenhagen, Frederiksberg, Denmark

Abstract

The disease septoria leaf blotch of wheat, caused by fungal pathogen Septoria tritici, is of worldwide concern. The fungus exhibits a hemibiotrophic lifestyle, with a long symptomless, biotrophic phase followed by a sudden transition to necrotrophy associated with host necrosis. Little is known about the systematic interaction between fungal pathogenicity and host responses at specific growth stages and the factors triggering the transition. In order to gain some insights into global transcriptome alterations in both host and pathogen during the two phases of the compatible interaction, disease transition was monitored using pathogenesis-related gene markers and H2O2 signature prior to RNA-Seq. Transcriptome analysis revealed that the slow symptomless growth was accompanied by minor metabolic responses and slightly suppressed defences in the host, whereas necrotrophic growth was associated with enhanced host responses involving energy metabolism, transport, signalling, defence and oxidative stress as well as a decrease in photosynthesis. The fungus expresses distinct classes of stage-specific genes encoding potential effectors, probably first suppressing plant defence responses/facilitating the symptomless growth and later triggering life style transition and inducing host necrosis/facilitating the necrotrophic growth. Transport, signalling, anti-oxidative stress mechanisms and apoplastic nutrient acquisition play important roles in the entire infection process of S. tritici. Our findings uncover systematic S. tritici-induced expression profiles of wheat related to specific fungal infection strategies and provide a transcriptome resource for studying both hosts and pathogens in plant-Dothideomycete interactions.

Introduction

Septoria tritici (teleomorph Mycosphaerella graminicola) is the causal agent of septoria leaf blotch, a foliar disease of wheat (Triticum aestivum) that poses a significant threat to global food production [1]. S. tritici penetrates host leaves through the stomata and grows slowly as filamentous hyphae in the intercellular spaces between the wheat mesophyll cells, typically up to 9–11 days [2,3]. Subsequently, the fungus exhibits a rapid switch to necrotrophy immediately prior to symptom expression at 12–20 days after penetration [35]. This transition is associated with induction of host defence responses that have been suggested to share characteristics with host programmed cell death (PCD) and differential regulation of signalling pathways [4,6]. The necrotrophic process involves cell-wall degradation and accumulation of H2O2, resulting in massive collapse of mesophyll tissue, leakage of nutrients from dying plant cells into the apoplastic spaces, rapid increase of fungal biomass and sporulation in characteristic necrotic foliar blotches [3,5,7,8]. The key features of this fungus, distinguishing it from most other current models, are a very long period of symptomless growth prior to a rapid transition to necrotrophy and that it remains extracellular with respect to host cells during the entire infection process without forming any specialised feeding structures [2,3,6]. However, the molecular basis underlying this transition remains unclear. This basis includes what the pathogen relies on to facilitate the long initial symptomless growth and subsequently trigger the transition of its lifestyle associated with host cell death and how the fungus copes with the activation of plant responses during the different infection phases. To address these issues, several investigations have attempted to understand S. tritici gene function, particularly the role of small secreted protein effectors. Likewise, fungal transcriptome profiles have been studied in adaption to environmental changes as well as host responses in cultivars with different susceptibility to infection. The research has been facilitated by the publication of the full genome sequence of S. tritici, revealing 10933 genes [9].

Effectors are believed to enable colonisation by the pathogen by interfering with a variety of plant-encoded virulence targets and to effectively suppress or activate defence responses. They can be recognised by disease resistance proteins, which subsequently trigger plant immunity [8]. In S. tritici, it has been suggested that the effectors may trigger disease transition, appearance of disease lesions and host cell death [10]. The only effector described in S. tritici is a LysM homologue shown to contribute to virulence in wheat [8]. Genes encoding secreted proteins with intragenic coding repeats within S. tritici populations were suggested to play potential roles as effectors [11]. In silico prediction of the secretome based on the full genome sequence identified 492 candidate virulence effectors and revealed a protein family encoding secreted (chloro)peroxidases, which is expanded within all family members of Mycosphaerellaceae [12]. Transcription profiling of S. tritici in planta using a microarray containing 2563 genes or expressed sequence tag (EST) libraries identified genes encoding cell-wall-degrading enzymes (CWDEs) and genes involved in signal transduction and transport [5] and revealed fungal physiological adaption with respect to membrane transport, chemical and oxidative stress mechanisms and metabolism during the rapid growth transition [4,13]. However, these transcriptome studies provide limited or ambiguous information on fungal pathogenicity based on a low number of identifications (150 genes in the biotrophic phase and max. 450 genes in the necrotrophic phase).

Studies on host responses in the wheat-S. tritici interaction revealed that H2O2 is important in the defence of wheat against the pathogen. During a compatible interaction, H2O2 levels increased dramatically and peaked at the late necrotrophic stage, implying that this was likely a stress-response and not involved in defence [3,7]. Other responses reported include expression of a wheat protein disulfide isomerase gene and pathogenesis-related (PR) genes [14,15] as well as structural defence responses [15]. Biochemical investigations furthermore showed DNA laddering, translocation of cytochrome c from mitochondria to the cytosol, loss of host cell membrane integrity, degradation of host total RNA and differential regulation of host mitogen-activated protein kinase (MAPK) pathways during symptom development in a compatible interaction [4,6]. Recently, large-scale shotgun proteomics and phosphoproteomics have been conducted on S. tritici-infected wheat leaves from the symptomless stage, suggesting that resistance is likely related to rapidly and intensively triggered signal transduction cascades resulting in a multiple-level activation of transcription and translation processes of defence responses [16].

Transcriptomics is an important tool to investigate the regulatory mechanisms between the host and the pathogen as it can shed light on the shifts of metabolic and cellular processes due to the interaction. Among the various technologies of transcriptome profiling, RNA-Seq has become a revolutionary tool and has recently been employed in the study of plant-pathogen interactions. It can distinguish between paralogous genes, detect and quantify transcripts of low or high abundance and identify transcript sequence polymorphisms and novel trans-splicing and splice isoforms [17]. Moreover, there is no strict requirement for a reference genome sequence, which is required in hybridization-based approaches such as microarray. Although transcriptome studies have been carried out in S. tritici in vitro and in planta, limited information about fungal gene expression has been obtained, in particular during the initial symptomless stage, due to low fungal biomass hardly detected by previous technologies. In addition, lack of large-scale systematic studies of host responses to S. tritici at RNA level strongly encouraged us to investigate the host transcriptome responses to infection.

Therefore, we used RNA-Seq in the present study to follow the transcriptional reprogramming with a key focus on both host and fungus simultaneously at distinct stages of the compatible interaction. These stages and disease transition were clearly defined by characterisation of fungal biomass, two essential plant PR genes/proteins (β-1,4-glucanase and chitinase) and ROS prior to transcriptome analysis. For the first time, we were not only able to study transcriptome changes (approx. 1800 transcripts) across the different growth phases in S. tritici, but also to obtain broad-spectrum insights into the host molecular responses corresponding to the distinct lifestyle phases of this pathogen. To our best knowledge, this is the first RNA-Seq study performed in a member of the Dothideomycetes. Our study has considerably expanded the knowledge obtained from previous gene expression investigations in S. tritici.

Materials and Methods

Plant growth and inoculation

Growth of the susceptible wheat cv. Sevin, preparation of inoculum of S. tritici isolate IPO323 and inoculation were performed as described by Shetty et al. [3]. Control plants were mock-inoculated with water. Approximately 20 leaves were collected from two separate pots, serving as one biological replicate. Two biological replicates for inoculated and control samples were harvested every day from 3 to 14 days after inoculation (dai) and immediately frozen in liquid nitrogen. The leaf samples were ground in liquid nitrogen and stored at -80 °C until use.

Histochemical staining for H2O2 and visualization of fungal structures

In vivo detection of H2O2 was carried out using 3,3’-diaminobenzidine (DAB, Sigma) as described by Shetty et al. [3]. The leaves were then cleared, stained by 0.1% Evans blue in lactoglycerol to visualise the fungal surface structures and studied by light microscopy [3].

Quantification of fungal biomass

Total DNA was extracted from infected and control samples using the DNeasy Plant Mini Kit (Qiagen, Venlo, The Netherlands). Fungal DNA was determined by qPCR using primers for S. tritici mating type gene 1-1 [18] with total DNA as template as previously described [19].

Enzyme assay and western blotting

Water-soluble protein was extracted in 50 mM sodium acetate (pH 5.2) at 4 °C. Protein concentration in the extracts was determined by the Bio-Rad Protein Assay (Bio-Rad) with bovine serum albumin as standard. β-1,3-glucanase activity was assayed as described by Shetty et al. [15]. Briefly, protein samples were incubated with 0.1% (w/v) laminarin in 0.05 M sodium acetate buffer (pH 5.2) for 15 min at 37 °C. The reaction was stopped by adding DNS reagent containing 0.5% (w/v) 3,5-dinitrosalicylic acid (Sigma) and 15% (w/v) potassium sodium tartrate tetrahydrate (Sigma) followed by boiling for 10 min. The absorbance was measured in an ELISA-reader at 540 nm. A standard curve relating the amount of glucose equivalents to the absorbance at 540 nm was used to determine the activity. Chitinase activity was measured by the Chitinase Assay Kit (Sigma) following the manufacturer’s protocol. Five µg protein were separated on CriterionTM XT Precast Gels (12% Bis-Tris, Bio-Rad) followed by blotting to nitrocellulose membranes (Whatman). Western blotting using rabbit antibodies against barley chitinase (kindly provided by David B. Collinge, University of Copenhagen) was performed as described by Yang et al. [20].

RNA extraction and qPCR

RNA extraction, gDNA removal, cDNA synthesis and qPCR were performed as previously described [20,21]. Primers for genes encoding wheat β-1,3-glucanase and chitinase were used according to Shetty et al. [15]. The wheat 25S rRNA gene (Forward primer: 5’-AAGGCCGAAGAGGAGAAAGGT-3’; Reversed primer: 5’-CGTCCCTTAGGATCGGCTTAC-3’) served as reference gene. Relative expression of genes was determined using the formula: relative expression = 2-(ΔCt target gene - ΔCt reference gene) where Ct refers to the threshold cycle in PCR [20]. All relative expression values of genes are reported as means ±SD. Analysis to reveal significant differences in gene expression in infected Sevin versus the respective controls at each time point was performed by Student’s t-test.

Construction of cDNA library, data filtering and transcript assembly

Total RNA (minimal 10 µg, 400 ng µl-1) was extracted from pooled leaf powder from two biological samples for each treatment using the SpectrumTM Plant Total RNA Kit (Sigma). RNA samples were subjected to DNase (Sigma) digestion to remove remaining DNA. mRNA was purified and enriched from total RNA samples by using poly-T oligo-coated magnetic beads (Invitrogen). Following purification and fragmentation, first-strand cDNA was generated using Superscript II reverse transcriptase (Invitrogen) and random hexamer primers. The cDNA was further converted into double stranded DNA using a strand-specific method according to Wang et al. [22]. After quality control of the cDNA libraries, pair-end sequencing analysis was carried out via Illumina HiSeq™ 2000 at the Beijing Genomics Institute (Shenzhen, China) according to the Illumina manufacturer’s protocol.

To ensure the accuracy of subsequent analysis, raw sequences were cleaned by removal of adaptors and sequencing errors. The reads containing the sequencing adaptor, more than 5% unknown nucleotides and more than 30% bases of quality value less than 10, were eliminated. This output was termed ‘clean reads’, which was used for alignment. All the clean reads were deposited in the National Center for Biotechnology Information (NCBI) and can be accessed in the Short Read Archive (SRA) linking to BioProject accession number: 196595. Prior to de novo assembly, the low-quality nucleotides at 3-ends were trimmed using a customer Perl script (CONDETRI: http://code.google.com/p/condetri) with the following parameters: hq=20, lq=10, frac=0.8, lfrac=0.1, minlen=40, mh=4, ml=1 and sc=64 [23]. To reduce the data complexity, two groups were generated from the clean reads from control (C4, C10 and C13) and fungal infected (I4, I10 and I13) samples for separate de novo assembly. The publicly available program Trinity (trinityrnaseq_r2012-05-18; http://trinityrnaseq.sourceforge.net/) was used for de novo transcriptome assembly of either group of trimmed clean reads to generate a set of contigs/transcripts with the following parameters: min_glue=2, SS_lib_type=RF, V=10, edge-thr=0.05, min_kmer_cov=3, path_reinforcement_distance=80 and group_pairs_distance=250 [24]. After the contigs from all samples were combined, the redundancy of contigs was removed and the contigs with partly overlap were elongated by the TGICL and Phrap assemblers [25]. The following parameters were used to ensure quality of assembly: a minimum of 95% identity between contigs, a minimum of 35 overlapping bases, a minimum of 35 scores and a maximum of 20 unmatched overhanging bases at sequence ends. Finally, based on the sequence similarity, the transcripts were divided into two classes: clusters (prefixed with ‘CL’) and singletons (prefixed with ‘unigene’). In each ‘CL’, the sequence similarity between transcripts was more than 70% and the transcripts were splice isoforms from one gene or paralogous gene.

Identification of fungal transcripts, sequence annotation and bioinformatics analysis

To identify the fungal transcripts, we BLAST searched all transcripts in the S. tritici genome (http://genome.jgi-psf.org/Mycgr3/Mycgr3.home.html). More than 80% in sequence similarity with regard to the alignment length was required to be considered as a fungal transcript. BLAST search against the fungal genome using a threshold of 50% similarity was performed as well. Fungal proteins were assessed for signal peptides using SignalP (http://www.cbs.dtu.dk/services/SignalP). For annotation, wheat and fungal transcripts, respectively, were aligned to four public databases [non-redundant protein database (Nr) in NCBI, non-redundant nucleotide database (Nt) in NCBI, Swiss-Prot and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database] by BLAST (E-value <10-5). Gene ontology (GO) classification was analysed by the Blast2GO software (v2.5.0) based on Nr annotation.

After annotation, the clean reads of each sample were mapped to all fungal and wheat transcripts using the SOAPaligner/soap2 software (BGI, China), allowing mismatches of no more than three bases, to obtain the transcripts from each sample. To minimize the influence of sequencing depth between samples, the total number of mapped reads was normalised by multiplying with a normalisation factor [26]. Expression level of the transcript was based on the number of unique match reads calculated and normalised to fragments per kilobase per million mapped fragments (FPKM) [27]. Analysis of differentially expressed wheat transcripts was performed as described by Chen et al. [28], which was based on the Poisson distribution [29] and normalisation for differences in sequencing depth between samples and gene length. A threshold of FPKM >2 in at least one sample, false discovery rate (FDR) <0.05 and at least two-fold change in FPKM of the transcript from infected versus control samples at minimum one time point was used to define differentially expressed wheat transcripts (performed by the in-house Perl script).

Principal component analysis (PCA) of sequencing data from all samples was performed by R script using log2FPKM values of transcripts in the samples. The expression profiles of the differentially expressed wheat transcripts [log2Ratio(infected/control)] and fungal transcripts (log2FPKM; FPKM>2 at minimum one time point) over the time course were determined by clustering analysis based on the k-means method using Euclidean distance. Heat map analysis of the differentially expressed wheat transcripts [log2Ratio(infected/control)] over the time course was performed using R script.

Results

Development of a bioassay to determine the disease transition

Given that the period of symptomless growth of S. tritici in wheat leaves can vary between experiments ranging from one week to several weeks and disease symptoms (Figure S1) are not uniform between leaves, we used the characteristic features of the biotrophic and necrotrophic phases in a bioassay to enable an interaction transcriptome analysis in the two distinct phases. The bioassay was based on characterisation of fungal biomass, H2O2 accumulation and plant PR-genes/proteins (β-1,4-glucanase and chitinase). The onset of a measurable increase in S. tritici growth has been shown to occur at day 9 and 10 [4,7]. Shetty et al. [3,15] additionally reported dramatically enhanced expression of β-1,3-glucanase and chitinase genes from 9 dai followed by extensive H2O2 accumulation, cell collapse and appearance of macroscopic symptoms from 12 dai during the compatible interaction. Therefore, we performed these assays over time in the compatible interaction to define the disease transition state. In the present study, fungal hyphae penetrated the stomata in the epidermis directly causing relatively little H2O2 accumulation up to day 8 (Figure 1). Increased accumulation of H2O2 in the substomatal cavities correlated with arrest of the progressing hypha was observed by day 9 and 10 (Figure 1), when strongly enhanced PR-gene expression and activity occurred in response to the fungus (Figure 2B–D). From 12 dai, the massive H2O2 accumulation in the mesophyll and hyphal development (Figure 1) were consistent with a significant increase in the levels of fungal biomass (Figure 2A). These results collectively show that following a long period of slow symptomless growth, the fungus switched to necrotrophy by day 9–10, advancing the increase in the biomass and massive H2O2 production by day 12.

thumbnail
Figure 1. Detection of hydrogen peroxide in S. tritici-infected wheat leaves from 3 to 14 dai.

DAB was used for staining H2O2 (red-brown staining). Fungal surface structures were stained with Evans blue. Arrows indicate attempted fungal penetration through the stomata.

https://doi.org/10.1371/journal.pone.0081606.g001

thumbnail
Figure 2. The bioassay to follow disease transition from 3 to 14 dai.

(A) Fungal biomass expressed as content of S. tritici DNA. Significant increases (P < 0.05) in levels of fungal DNA are indicated by asterisks. (B) Wheat β-1,3-glucanase and chitinase gene expression analysed by qPCR. (C) β-1,3-glucanase and chitinase protein activities expressed as released glucose (μM) and OD405 measurement based on the assays, respectively. Significant differences (P < 0.05) between infected and control samples are indicated by asterisks. (D) Chitinase protein expression in protein extracts from control (upper panel) and S. tritici-infected leaves (lower panel) analysed by western blotting. Bands indicate chitinase. Representatives of two biological replicates at each time point are shown for the control and infected plants.

https://doi.org/10.1371/journal.pone.0081606.g002

A snapshot of transcript profiles during the wheat-S. tritici interaction

It is believed that the changes of lifestyle in a fungal pathogen and process of disease development are accompanied by stage-specific molecular events [30]. To gain comprehensive information on the molecular and physiological bases during particular stages of the interaction between wheat and S. tritici, we explored the complexity of transcriptional reprogramming of both host and pathogen at three time points (day 4, 10 and 13) by RNA-Seq analysis. The time points were selected based on the bioassay of marker genes/proteins and H2O2 accumulation and represented the initial symptomless phase, the intermediate transition stage and the necrotrophic phase, respectively.

Table 1 summarises the sequencing data from S. tritici-infected wheat leaves and their respective controls. The mRNA sequencing workflow is shown in Figure S2. After the pre-processing of the reads, we obtained more than 12 million clean pair-end sequencing reads and 2G high-quality bases (Q20 bases approximately 97%) from both control samples and the mixed transcriptome samples (Table 1). In order to conduct dual sequencing analysis without the host reference genome, we used the de novo assembly method for control and infected samples separately. In total, 88770 non-redundant transcripts were obtained from all the samples including 12695 clusters (38680 transcripts) and 50090 singletons, the total length and N50 length of which were 63937118 bp and 1109 bp, respectively. Subsequently, we filtered the fungal transcripts by mapping all assembled transcripts to the fungal genome, resulting in identification of 1829 and 1811 fungal transcripts with a threshold of 50% and 80% sequence similarity, respectively (Table S1). In total, 313 fungal transcripts encoded proteins that were predicted to contain signal peptides. Functional categories were assigned to all annotated transcripts using GO classification and KEGG pathway analysis, resulting in 67610 and 1626 annotated wheat and fungal transcripts, respectively. Fungal (Table S1) and wheat (Tables S2, S3) transcripts with known function were grouped into 14 and 19 GO functional categories based on GO biological processes as well as 90 and 127 pathways, respectively.

C4C10C13I4I10I13
Number of reads123304231237547812583655127577791242676413020331
Total base pairs (Mbp)222022282265229722372344
Q20 (%)95.6597.8495.6197.4396.0698.12
Total number of contigs8931559389
Length of contigs (bp)7055183531332619
N50 length of contigs (bp)1233692
Number of mapped fungal transcripts91132157816921822
Number of wheat transcripts 819767940779642790327792775461
Highly expressed fungal transcripts 000107142761
Highly expressed wheat transcripts 1015811908956393721271513690

Table 1. Overview of sequencing transcriptome data from wheat leaves.

The transcripts with FPKM>10 are considered highly expressed. C4, control at 4 dai; C10, control at 10 dai; C13, control at 13 dai; I4, S. tritici-infected wheat at 4 dai; I10, S. tritici-infected wheat at 10 dai; I13, S. tritici-infected wheat at 13 dai.
CSV
Download CSV

Mapping the reads of each sample to wheat and fungal transcripts resulted in the identification of 1578, 1692 and 1822 fungal transcripts from infected samples at 4, 10 and 13 dai, respectively, as well as approximately 80000 non-redundant plant transcripts from each sample (Table 1). Nine, eleven and thirty-two low-expressed transcripts encoding tubulins, ATPases, ribosomal proteins and cytochrome b, respectively, were mapped as fungal transcripts from control samples, which is very likely due to the homology of these conserved genes between plant and the fungus. To further examine an overall relationship of the transcriptome data from different samples and effects of both fungal infection and interaction phases, PCA was conducted and clearly revealed that the greatest changes in host transcriptome were caused by the fungal infection following the interaction process (Figure 3). In particular, the 13 day-infected sample was separated from all other samples by principal component 1, accounting for 73.9% of the variance, strongly suggesting the accelerated host-pathogen interplay during the necrotrophic stage.

thumbnail
Figure 3. Principal component analysis of the transcriptome data.

Principal components (PCs) 1, 2 and 3 account for 73.9, 11.6 and 7% of the variance, respectively. C4, control at 4 dai; C10, control at 10 dai; C13, control at 13 dai; I4, S. tritici-infected wheat at 4 dai; I10, S. tritici-infected wheat at 10 dai; I13, S. tritici-infected wheat at 13 dai.

https://doi.org/10.1371/journal.pone.0081606.g003

Fungal transcriptome during the compatible interaction

Figure 4 presents ten clusters of expression profiles of fungal transcripts as well as functional category distribution frequencies of the transcripts in each cluster. Not surprisingly, the transcripts with known function annotated in metabolic and cellular processes were dominating in each cluster. The clusters fell into three groups based on the major profiling over the time course: group A consisted of cluster 1, 5, 9 and 10, the majority of 660 transcripts showing constant abundance; group B consisted of cluster 2, 3, 6, 7 and 8, the majority of 963 transcripts showing increased abundance; group C consisted of cluster 4, the majority of 206 transcripts showing decreased abundance at 10 dai and highest abundance at 13 dai. A substantially higher number of transcripts annotated in signalling, responses to stimuli, reproduction, localisation, metabolic and cellular processes and cellular component organisation or biogenesis were represented in group B compared to other two groups (Figure 4). In addition, several transcripts involved in primary and energy metabolism and a considerable number of highly abundant transcripts annotated as ribosomal proteins were present in group B (Table S1). The FPKM values of housekeeping transcripts such as ribosomal and tubulin transcripts increased approximately 10-fold at 13 dai compared to the symptomless stage. These data suggest a highly active metabolic growth and an increase in fungal biomass at the destructive necrotrophic phase.

thumbnail
Figure 4. Functional category distribution in the ten expression clusters of fungal transcripts.

The fungal transcripts were identified from S. tritici-infected wheat at 4, 10 and 13 dai. Expression level of the transcripts was calculated by log2FPKM. Histogram representation of the category distribution is expressed as percentage of the amount of transcripts belonging to the cluster. Transcripts coding for unknown products were included in the analysis. The clusters fall into three groups (A, B and C) based on the major profiling over the time course.

https://doi.org/10.1371/journal.pone.0081606.g004

Since the transcripts with the highest expression levels at specific growth stages likely play important roles in the control of the disease transition, we selected several transcripts with the highest expression levels at 4 dai (biotrophy) or 13 dai (necrotrophy) and listed them in Table 2. Besides the transcripts annotated as hypothetical proteins, most of which contained signal peptides and more than six cysteines, the transcripts were involved in pathogenesis and virulence, metabolism, signalling, transport, stress and defence, cell adhesion and cell wall degradation and remodelling (Table 2). Notably, these highly expressed transcripts at 4 and 13 dai precisely belonged to group A and B, respectively, indicating that the transition of S. tritici lifestyle was likely accompanied and regulated by the expression of stage-specific sets of genes. More interestingly, although seven transcripts involved in host cell wall degradation including glycoside hydrolase families 5, 16 and 76 in group A and glycoside hydrolase families 3 and clan CH-D in group C were identified, none of them were of high abundance at either stage (FPKM <8; Table S1).

FPKM
IDAnnotation4 d10 d13 dBiological processSC
Transcripts highly expressed (FPKM ≥ 8) at 4 d
96543Hydrophobin 29.923.115.6Pathogenesis and virulenceYes5
70312Peptidase S8 and S5310.212.423.3MetabolismYes5
74298Glucoamylase15.916.019.0MetabolismYes5
77672Esterase/lipase/thioesterase19.515.88.8MetabolismNo5
81692β-1,6-N-acetylglucosaminyltransferase13.25.76.9MetabolismYes5
104538β-1,6-N-acetylglucosaminyltransferase52.044.547.8MetabolismNo1
86867NADPH-dependent FMN reductase33.736.0102.0MetabolismYes1
100165Fumarate lyase39.311.832.3MetabolismNo1
108219Autophagy-related protein 811.99.928.6MetabolismNo5
109517Carbohydrate kinase, PfkB17.06.410.1MetabolismNo5
102849Homeobox transcription factor prospero14.211.237.2MetabolismYes1
44262Thiamine biosynthesis Thi4 protein11.06.913.0MetabolismNo5
104937Pyridoxamine 5'-phosphate oxidase17.85.314.2MetabolismNo5
104409Similar to bacterial rhodopsin29.340.370.2SignallingNo1
109991 93710Serine/threonine protein kinase62.457.3247.0SignallingYes1
99917Serine/threonine protein kinase19.912.916.9SignallingYes5
97868Tyrosine protein kinase16.625.472.0SignallingYes1
110667PAS-domain protein63.870.2222.0SignallingYes1
39898Aquaporin24.910.811.8TransportNo5
42411GPR1/FUN34/yaaH membrane protein20.920.112.4TransportNo5
58200ABC transporter11.16.810.2TransportNo5
76071Phosphate transporter19.913.612.3TransportYes5
105313Sugar transporter superfamily20.812.912.4TransportNo5
64952Cyclophilin10.511.940.6Stress and defenceNo1
72265Cyclophilin12.24.511.0Stress and defenceNo5
106153Carbohydrate-binding WSC18.219.125.2Stress and defenceNo1
59219Haem peroxidase9.75.72.7Stress and defenceNo10
94368Chloroperoxidase10.914.56.7Stress and defenceYes5
101235Chloroperoxidase12.116.514.8Stress and defenceYes5
103393Allergen V5/Tpx-1 related23.124.863.8Stress and defenceYes1
103593Cu2+/Zn2+ superoxide dismutase 112.820.157.3Stress and defenceNo1
105791Aldehyde dehydrogenase28.619.420.4Stress and defenceNo1
101210Zinc-containing alcohol dehydrogenase9.62.34.3Stress and defenceNo10
56432Fasciclin 9.11.85.2Cell adhesionYes10
104730von Willebrand factor 24.421.631.4Cell adhesionYes1
100647Chitinase9.61.49.4Cell wall degradation and remodellingYes10
106219Glycolipid anchored surface protein GAS112.43.321.3Cell wall degradation and remodellingYes5
90001Yeast PIR protein repeat-like38.436.773.6Cell wall degradation and remodellingYes1
69789Hypothetical protein21.69.919--Yes5
94290Hypothetical protein11288.154.4--Yes1
106329Hypothetical protein25.422.59.5--Yes5
103950Hypothetical protein14.77.626.1--Yes5
111606Hypothetical protein8.210.014.3--Yes5
Transcripts highly expressed (FPKM ≥ 10) at 13 d
39947Cerato-platanin1.22.944.3Host immune systemYes3
105487Peptidoglycan-binding LysM1.212.7119.8Pathogenesis and virulenceYes3
111221Peptidoglycan-binding LysM3.012.1105.3Pathogenesis and virulenceYes7
88451Necrosis-inducing protein NPP10.73.510.5Pathogenesis and virulence
104853Ubiquinol-cytochrome C reductase4.04.523.5MetabolismYes7
107092β-1,6-N-acetylglucosaminyltransferase2.62.420.9MetabolismYes8
108724Cysteine peptidase3.33.929.1MetabolismNo7
95071Nucleoside diphosphate kinase5.13.433.3MetabolismNo7
57734G-protein beta WD-40 repeat3.11.323.9SignallingNo8
86648G-protein beta WD-40 repeat2.94.922.1SignallingNo7
69026Translationally controlled tumor protein5.77.331.9SignallingNo7
78859EGF-like0.92.534.4SignallingNo3
103900EGF-like0.221.2120.3SignallingYes3
106452Serine/threonine protein kinase3.511.871.2SignallingYes7
100205Phox-like5.55.224.2SignallingNo7
65946Amino acid permease0.65.329.7TransportNo3
73144Ammonium transporter0.17.124.4TransportNo3
77435Porin, eukaryotic type6.03.126.8TransportNo7
65963Cytochrome b561 0.10.726.6Transport Yes2
90089Ankyrin0.17.042.7Transport Yes3
82936Plastocyanin-like4.96.332.8TransportYes7
98411Cytochrome P4500.43.420.6Stress and defenceNo3
55916Glutathione S-transferase3.42.520.6Stress and defenceNo8
71387Alkyl hydroperoxide reductase3.63.629.4Stress and defenceNo7
101587Thioredoxin2.25.041.4Stress and defenceYes7
104975Peroxiredoxin3.44.626.3Stress and defenceNo7
72449Heat shock protein Hsp702.76.223.3Stress and defenceNo7
83835Heat shock protein Hsp701.63.324.3Stress and defenceYes8
103427von Willebrand factor4.512.8141.1Cell adhesionYes7
106780von Willebrand factor2.63.336.7Cell adhesionYes7
102341Glycolipid anchored surface protein GAS17.34.728.3Cell wall degradation and remodellingYes7
102481Yeast PIR protein repeat-like0.72.343.1Cell wall degradation and remodellingYes3
104794Chitinase0.42.036.0Cell wall degradation and remodellingYes3
44596Kinesin, motor region4.121.540.4LocomotionNo7
83064Hypothetical protein0.63.456.2--Yes3
79161Hypothetical protein1.43.836.8--Yes3
89705Hypothetical protein0023.0--Yes7
92097Hypothetical protein0.76.724.3--Yes3
97031Hypothetical protein0.63.622.4--Yes3
110220Hypothetical protein2.49.868.2--Yes7
102792Hypothetical protein5.463.7277.0--Yes1
104404Hypothetical protein0.412.260.0--Yes3
104444Hypothetical protein1.024.9102.7--Yes3
105825Hypothetical protein0.814.974.9--Yes3
105826Hypothetical protein1.59.735.7--Yes3
108482Hypothetical protein0.10.320.7--Yes2
107904Hypothetical protein0.33.743.8--Yes3
107286Hypothetical protein0.93.329.8--Yes3

Table 2. Expression patterns of selected S. tritici transcripts during infection.

FPKM indicates the expression level of the transcript. Most of the hypothetical proteins are small cysteine-rich proteins.
S, SignalP; C, cluster number in Figure 4; d, days after inoculation.
CSV
Download CSV

Furthermore, since fungal secreted proteins may have potential roles as effectors, 313 transcripts encoding proteins containing signal peptides, the majority of which were hypothetical proteins, peptidases, glycoside hydrolases, antioxidants, kinases, heat shock proteins and oxidoreductases, were clustered to obtain the expression patterns (Figure S3A). Their expression levels increased (clusters 1, 3, 4, 7–10) or remained constant (clusters 2, 5 and 6) over the time course. Likewise, of 492 candidate virulence effectors previously predicted in silico [12], 246 fungal transcripts were identified in the present study with increased or constant expression levels during infection (Figure S3B; Table S1).

Wheat leaf transcriptome in response to infection

In total, 33314 wheat transcripts, fulfilling the selection criteria, were identified in response to the pathogen (Table S2). Cluster analysis and category distribution were performed to assess the overall patterns of host transcriptome changes over the time course of infection (Figure S4). Clusters 1 and 3 contained 3266 transcripts mainly regulated by the fungus at 4 dai. Clusters 5, 6, 7 and 9 contained 16543 transcripts mainly regulated at 13 dai. Clusters 2, 4 and 8 contained 11324 transcripts mainly regulated at 10 dai. Additionally, cluster 10 contained 2181 transcripts down-regulated and up-regulated at 4 and 13 dai, respectively. The major transcripts in each cluster were annotated in signalling, response to stimuli, reproduction, regulation of biological process, metabolic and cellular processes, locomotion, localisation, development, cellular component organisation or biogenesis and biological adhesion. Considerably higher numbers of transcripts involved in metabolism, transport, defence and stress, signalling and cellular component organisation changed in abundance at 10 and 13 dai (Figure S4; Table S2), suggesting a strong induction of metabolic and defence-activities from 10 dai, which was associated with increased H2O2 production (Figure 1).

To gain further insights into the responses of specific functions of host genes in fungal infection at different stages, we carried out a heat map analysis of the transcripts. This included genes specifically annotated in primary metabolism (e.g., metabolism of carbohydrates, nitrogen and amino acids as well as photosynthesis), transport (e.g., ATPases, ABC transporters, aquaporins and sugar, amino acid, peptide, nucleic acid and lipid transporters) and signalling [e.g., membrane-bond receptors, receptor kinases, G proteins, phospholipids, calcium-dependent protein kinases, MAPK pathways, WRKYs and other transcription factors, hormone pathways (salicylic acid, jasmonate, ethylene, auxin, gibberellic acid and brassinosteroid)] as well as genes involved in stress and defence (e.g., oxidative stress, detoxification, PR-genes, heat shock proteins, the phenylpropanoid pathway, LRR domain proteins and other stress-related genes) (Figure 5). The analysis revealed two distinct gene expression patterns of the host in response to the fungus during the biotrophic and necrotrophic stages. At the early symptomless stage, the fungus caused minor alterations in the expression of host transcripts involved in primary metabolism, transport and signalling, although some transcripts significantly changing in abundance were observed. By contrast, a considerable number of transcripts involved in defence and defence-related signalling, encoding proteins such as antioxidants, PR-proteins, heat shock proteins, disease resistance proteins, WRKY transcription factors and NAC domain proteins were slightly down-regulated at this stage, even though some, such as receptor-like kinases and phytohormone-inducing proteins, could be up-regulated (Figure 5; Table S2). At the necrotrophic stage, great changes were observed in expression of host transcripts involved in host metabolism, transport, signalling and defence with up-regulation of a substantial number of transcripts (Figure 5). Furthermore, it is worth to note that levels of the major transcripts annotated in carbohydrate metabolism were up-regulated at this stage, whereas major transcripts annotated in photosynthesis were down-regulated. Additionally, several transcripts involved in anti-oxidative stress were strongly down-regulated at the late necrotrophic stage (Table S2), which was correlated to massive accumulation of H2O2 (Figure 1).

thumbnail
Figure 5. Heat map visualization of differentially expressed wheat transcripts.

The plant transcripts are mainly involved in primary metabolism, signalling, transport and stress and defence in response to S. tritici. Differential expression patterns are based on the log2 fold changes of transcript abundance in S. tritici-infected wheat versus controls at 4, 10 and 13 dai.

https://doi.org/10.1371/journal.pone.0081606.g005

Discussion

Understanding the molecular mechanisms underlying fungal infection and host plant defence is a prerequisite for understanding the host-pathogen interaction and can contribute to development of new strategies of crop protection against pathogen. A global study of host and fungal transcriptional reprogramming during the interaction allows us to gain insights into such mechanisms. Here, we report the first use of RNA-Seq technique in the compatible interaction between wheat and S. tritici, with particular focus on the metabolic and physiological alterations of both host and pathogen during the shift from biotrophic to necrotrophic growth stages.

It has been suggested that the transition from slow biotrophic growth to rapid necrotrophic growth is associated with distinct expression patterns of genes regulated by nutrient status, revealed by comparison of transcriptomes between in planta transition/necrotrophic stage and in vitro nutrition-rich/limited conditions [4,13]. Our interaction transcriptome analysis directly demonstrates two different phases of the compatible interaction: i) slow symptomless growth associated with minor host molecular responses and ii) necrotrophic feeding/high metabolic growth associated with enhanced metabolic and defence responses in plant. The symptomless growth of S. tritici was achieved by suppressing or attenuating host defence responses at transcription level (Figure 5), despite up-regulation of some defence-related genes such as receptor kinases, phytohormone-inducing proteins and PR-2 and 3 as well as accumulation of stress-related compounds such as ROS. This symptomless growth is analogous with most biotroph lifestyles of both repression of plant defence responses and induction of host specific genes. However, accumulation of several defence-related proteins and down-regulation of several proteins involved in photosynthesis have been observed at the biotrophic stage as revealed by proteomics [16]. Down-regulation of several defence-related transcripts and no alteration in the levels of several photosynthesis-related transcripts appear not to be correlated to the regulation behaviour of their protein products. The poor correlation between transcript and protein levels has been well demonstrated due to the regulation of transcription, RNA processing, translation and protein modifications and turnover and furthermore suggests that integrated ‘omics’ studies are required to comprehensively understand the molecular processes in the biological system. On the other hand, the necrotrophic growth of S. tritici triggers the plant cells to create a highly defensive and metabolic active environment, generating ROS and strongly inducing defence, signalling and metabolism-related genes. The crosstalk between induced signalling events, including receptor-mediated signal perception, protein phosphorylation, ion fluxes, production of ROS, generation and regulation of secondary signalling molecules (e.g., phytohormones and WRKYs), leading to the activation of defence genes, has been well established as the mechanisms of host plant defence to different pathogens [31]. More interestingly, the highly oxidative cell environment in wheat during the necrotrophic stage was accompanied by down-regulation of large number of plant transcripts encoding ROS-scavenging enzymes, thus allowing the ROS to accumulate. For some necrotrophic pathogens, it has been suggested that ROS accumulation and depletion of host antioxidants can result in increased susceptibility as well as host necrosis [32]. Interestingly, this is not the case for S. tritici that is able to tolerate the high amounts of H2O2 accumulating during the necrotrophic growth phase, but does not thrive in this environment since scavenging of H2O2 boosts fungal growth [3,7]. Furthermore, levels of transcripts involved in photosynthesis were suppressed during the necrotrophic phase, whereas carbohydrate and energy metabolism was up-regulated in the host and associated with enhanced defence responses [7]. It has been demonstrated that the activation of plant defence responses is cost-intensive and causes a decrease in photosynthesis and an increased demand for assimilates, energy, reducing equivalents and carbon skeleton components that are provided by the primary metabolic pathways [33]. Simultaneously, the pathogen will put pressure on the host carbohydrate metabolism for its own needs, which will further increase the demand for assimilates.

To date, little is known of the molecular mechanisms in S. tritici that promote an early symptomless, biotrophic phase and a late highly destructive necrotrophic stage or those that mediate the transition between the two. Previous investigations focused on the secreted protein effectors for facilitating the initial symptomless growth and triggering PCD [10]. We observed strong, differential and growth stage-specific expression of transcripts encoding several small secreted cysteine-rich proteins with unknown function at the two stages, a hydrophobin 2 during the biotrophic stage and two LysM proteins, a necrosis-inducing protein and a cerato-platanin during the necrotrophic stage, which may have potential roles in pathogenesis or suppressing plant defence. This observation suggests that S. tritici secretes distinct classes of proteins that may first suppress plant defence responses and later induce host necrosis. It has been reported that some fungal effectors, for instance, P. infestans Avr3a and SNE1 and P. sojae Avr1b, can suppress PCD or host defences during biotrophic stages [3436]. On the other hand, the P. infestans PiNPP1.1 gene encoding a NLP effector, which is only highly expressed during necrotrophy, induces necrosis and is suppressed by the SNE1 effector [37]. Both LysM proteins identified in the present study have been functionally characterised, one of which was proposed to be an effector to facilitate the symptomless leaf colonisation of S. tritici through its activity in suppressing chitin-mediated plant defences [8]. Surprisingly, LysM transcripts were highly expressed during the necrotrophic stage here, which was associated with increased fungal chitin/biomass and plant defence responses. Thus, our data suggest that S. tritici manipulates the suppression of chitin-mediated wheat defences during the entire infection process and may be important for disease development due to the increased level of biomass or chitin. Expression of a necrosis-inducing protein has also been found in the oomycete Phytophthora sojae during transition from biotrophy to necrotrophy in soybean, which suggests that this protein facilitates the plant colonisation during the necrotrophic growth [38]. Three secreted necrosis-inducing protein effectors have been identified in barley pathogen Rhynchosporium commune and shown their contribution to disease development [39]. Cerato-platanin in phytopathogenic fungi acts as a phytotoxin and a pathogen-associated molecular pattern (PAMP) due to the induction of defence responses including phytoalexin synthesis and plant cell death [40]. Their strong expression at the transition and the necrotrophic stages, which was associated with the enhanced host defence responses, indicates a potential role in induction of the disease transition and/or facilitation of the necrotrophic growth. Fungal secreted small cysteine-rich hydrophobins have been implied in pathogenesis through acting as toxins or attaching fungal structures to the host surface [41]. Eight hydrophobin-like proteins are predicted in the annotated S. tritici genome, four of which are secreted [12]. In addition to one highly expressed hydrophobin shown in Table 2, four additional hydrophobins were identified with one containing a signal peptide here (Table S1). We speculate that hydrophobins play a role in facilitating the initial symptomless growth through its ability of mediating fungal hypha to attack the host cells.

The coverage of the transcripts encoding the proteins involved in transport, signalling and stress responses among the identified fungal transcripts was remarkable at both stages, indicating their important roles in facilitating both biotrophic and necrotrophic growth. Strongly expressed transcripts involved in defence-related signalling included G proteins, serine/threonine protein kinases (likely involved in MAPK pathways or Ca2+ signalling), PAS-domain proteins functioning as sensors of oxygen, redox potential, light and other stimuli [42] and EGF-like proteins known for a role in the immune system [43]. Some of these identifications differed from the 74 transcripts involved in four MAPK pathways, the cAMP-dependent pathways and G-signalling identified in a previous in planta EST sequencing study [5]. Disruption of S. tritici signalling genes can significantly affect and regulate the invasive growth, penetration or pathogenicity [44]. With respect to highly expressed transport-related transcripts, different classes corresponding to two growth stages were identified. The nutrient uptake and transport facilitation appeared to be active particularly during the biotrophic stage in contrast to the primary metabolism-related transport during the necrotrophic growth. As described above, ROS accumulated in the host from 3 to 14 dai (Figure 1) and is considered as one of the primary responses to S. tritici, starting at the infection site during the biotrophic stage and spreading to the entire tissue at the late necrotrophic stage [3]. ROS are effective in stopping growth of biotrophic pathogens. However, the roles of ROS in interactions with necrotrophic fungi are still ambiguous. It can both benefit and inhibit infection by necrotrophs, which may be able to produce ROS themselves and stimulate the ROS production from the hosts as well [7,45]. To survive in the harsh oxidative environments and colonise host cells, S. tritici, therefore, has to develop mechanisms to scavenge ROS and protect against ROS-induced damage, particularly during the necrotrophic stage when massive ROS accumulation occurs [3]. This was evident in the present study by increased expression levels of several fungal transcripts during infection encoding ROS-scavenging proteins including peroxidases, Cu2+/Zn2+ superoxide dismutases, glutathione S-transferases, hydroperoxide reductases, thioredoxins and peroxiredoxins as well as other stress-response proteins like cyclophilins and heat shock proteins (Table S1). Furthermore, excessive ROS production in the fungal cells after expose to the stress can give rise to the accumulation of aldehydes and alcohols [46]. S. tritici may rescue itself by increased expression of aldehyde dehydrogenase and alcohol dehydrogenase for detoxification of alcohols and aldehydes (Table 2).

During the entire infection cycle, S. tritici grows intercellularly throughout the leaf mesophyll cell layer and the host cell walls appear not to be breached [2,3], although the fungus can express a series of genes encoding CWDEs in planta, in particular a large number during the very late necrotrophic stage when the foliar tissue is fully necrotic, reflecting that stage of pathogenesis [4,5]. CWDEs are known to play a role in nutrient acquisition, host colonisation and pathogenicity in phytopathogenic fungi. However, scanning the S. tritici genome indeed revealed a relatively small number of genes annotated as CWDEs compared to some saprophytes and phytopathogenic fungi like Fusarium graminearum, Magnaporthe oryzae and Stagonospora nodorum [9,12]. Little is understood about the biological explanations for the strictly apoplastic colonising hyphae of S. tritici which do not breach host cell walls. In the present study, we identified relatively few transcripts encoding CWDEs expressed during the biotrophic and initial necrotrophic stages out of the total identifications and none of them was of high abundance, indicating that degradation of host cell walls may not play an essential role for nutrient acquisition and explaining the reduced capacity for breaching living host cell walls. Correspondingly, Rohel et al. [47] found that during the transition to necrotrophic growth where sporulation started, the fungus actually starved, indicating deficiency in obtaining sufficient nutrients. However, the fungus has to survive on the nutrients in the apoplast full of proteins, metabolites, sugars and ions [47]. Here, the identification of several transcripts annotated as peptidases, lipases and amylases with high abundance at two stages (Table 1, S2) suggests that degradation of host cell proteins, lipids and sugars can be an important source of fungal nutrients during colonisation. This is in agreement with the genomic mapping analysis revealing that several protease encoding gene families are expanded and 39 proteins have functions related to protein degradation [12].

In summary, the disease transition was clearly defined by characterising accumulation of H2O2, PR-proteins and fungal biomass in order to carry out the transcriptome analysis at the biotrophic and necrotrophic stages of the compatible interaction between wheat and S. tritici. RNA-Seq enabled the identification of considerable numbers of fungal and wheat transcripts in response to the infection at both stages. The results suggest that successful leaf colonisation is achieved by suppressing host defence responses during the symptomless biotrophic stage, activation of signalling, transport and antioxidative mechanisms, apoplastic nutrient acquisition and secretion of distinct classes of stage-specific molecules including potential effectors. The necrotrophic growth triggers enhanced host responses including anti-stress and defences, energy metabolism, signalling and transport as well as decreased photosynthesis. This transcriptome study paves the way for more detailed work on S. tritici pathogenicity and wheat resistance to S. tritici, which will be important for future crop protection.

Supporting Information

Figure S1.

Macroscopic appearance of susceptible wheat leaves after treatment with water (controls) or S. tritici. The leaves were photographed from 3 to 14 dai.

https://doi.org/10.1371/journal.pone.0081606.s001

(TIF)

Figure S2.

Schematic representation of RNA-Seq analysis of the interaction transcriptome between wheat and S. tritici. Total RNA was extracted from control and S. tritici-infected samples at 4, 10 and 13 dai. Illumina RNA-Seq was performed on cDNA libraries generated from each sample. All pre-processed reads of controls as well as the infected samples were combined and assembled. Subsequently, the transcripts from all the samples were combined, mapped to the fungal genome to filter out fungal transcripts prior to annotation. Finally, the clean reads of each sample were mapped to wheat and fungal transcripts. Bioinformatics analysis was conducted to obtain expression levels of both wheat and fungal transcripts in each sample and define differentially expressed wheat transcripts at each time point.

https://doi.org/10.1371/journal.pone.0081606.s002

(TIF)

Figure S3.

Expression profiles of S. tritici transcripts encoding potential effector proteins identified from S. tritici-infected wheat at 4, 10 and 13 dai. Expression levels of the transcripts were calculated by log2FPKM. (A) Ten expression clusters of 313 fungal transcripts encoding proteins containing signal peptides. (B) Ten expression clusters of 246 fungal transcripts which have been identified in the present study and in silico predicted as effectors by do Amaral et al. [12].

https://doi.org/10.1371/journal.pone.0081606.s003

(TIF)

Figure S4.

Functional category distribution in the ten expression clusters of differentially expressed wheat transcripts. Differentially expressed plant transcripts in response to S. tritici at 4, 10 and 13 dai are presented by expression ratios of the transcripts in infected (I) versus control (C) samples in log2-scale. Histogram representation of the category distribution is expressed as percentage of the amount of transcripts belonging to the cluster. Transcripts coding for unknown products were included in the analysis.

https://doi.org/10.1371/journal.pone.0081606.s004

(TIFF)

Table S1.

List of identified S. tritici transcripts during infection of wheat leaves at 4, 10 and 13 dai.

https://doi.org/10.1371/journal.pone.0081606.s005

(XLSX)

Table S2.

List of differentially expressed wheat transcripts in response to S. tritici at 4, 10 and 13 dai.

https://doi.org/10.1371/journal.pone.0081606.s006

(XLSB)

Table S3.

List of identified wheat transcripts with no change in abundance during S. tritici infection at 4, 10 and 13 dai.

https://doi.org/10.1371/journal.pone.0081606.s007

(XLSB)

Acknowledgments

We thank Daniel Amby (University of Copenhagen) for valuable scientific discussion and Haiyan Song (BGI) for the project coordination.

Author Contributions

Conceived and designed the experiments: FY. Performed the experiments: FY WL HJLJ. Analyzed the data: FY WL HJLJ. Contributed reagents/materials/analysis tools: FY WL HJLJ. Wrote the manuscript: FY WL.

References

  1. 1. Eyal Z (1999) The Septoria tritici and Stagonospora nodorum blotch diseases of wheat. Eur J Plant Pathol 105: 629–641. doi:https://doi.org/10.1023/A:1008716812259.
  2. 2. Kema GHJ, Yu DZ, Rijkenberg FHJ, Shaw MW, Baayen RP (1996) Histology of the pathogenesis of Mycosphaerella graminicola in wheat. Phytopathology 86: 777–786. doi:https://doi.org/10.1094/Phyto-86-777.
  3. 3. Shetty NP, Kristensen BK, Newman MA, Møller K, Gregersen PL, et al. (2003) Association of hydrogen peroxide with restriction of Septoria tritici in resistant wheat. Physiol Mol Plant P 62: 333–346.
  4. 4. Keon J, Antoniw J, Carzaniga R, Deller S, Ward JL et al. (2007) Transcriptional adaptation of Mycosphaerella graminicola to programmed cell death (PCD) of its susceptible wheat host. Mol Plant Microbe Interact 20: 178–193. doi:https://doi.org/10.1094/MPMI-20-2-0178. PubMed: 17313169.
  5. 5. Kema GH, van der Lee TA, Mendes O, Verstappen EC, Lankhorst RK et al. (2008) Large-scale gene discovery in the septoria tritici blotch fungus Mycosphaerella graminicola with a focus on in planta expression. Mol Plant Microbe Interact 21: 1249–1260. doi:https://doi.org/10.1094/MPMI-21-9-1249. PubMed: 18700829.
  6. 6. Rudd JJ, Keon J, Hammond-Kosack KE (2008) The Wheat mitogen-activated protein kinases TaMPK3 and TaMPK6 are differentially regulated at multiple levels during compatible disease interactions with Mycosphaerella graminicola. Plant Physiol 147: 802–815. doi:https://doi.org/10.1104/pp.108.119511. PubMed: 18441220.
  7. 7. Shetty NP, Mehrabi R, Lütken H, Haldrup A, Kema GHJ et al. (2007) Role of hydrogen peroxide during the interaction between the hemibiotrophic fungal pathogen Septoria tritici and wheat. New Phytol 174: 637–647. doi:https://doi.org/10.1111/j.1469-8137.2007.02026.x. PubMed: 17447918.
  8. 8. Marshall R, Kombrink A, Motteram J, Loza-Reyes E, Lucas J et al. (2011) Analysis of two in planta expressed LysM effector homologs from the fungus Mycosphaerella graminicola reveals novel functional properties and varying contributions to virulence on wheat. Plant Physiol 156: 756–769. doi:https://doi.org/10.1104/pp.111.176347. PubMed: 21467214.
  9. 9. Goodwin SB, M’Barek SB, Dhillon B, Wittenberg AHJ, Crane CF et al. (2011) Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis. PLoS Genet 7: e1002070. PubMed: 21695235.
  10. 10. Deller S, Hammond-Kosack KE, Rudd JJ (2011) The complex interactions between host immunity and non-biotrophic fungal pathogens of wheat leaves. J Plant Physiol 168: 63–71. doi:https://doi.org/10.1016/j.jplph.2010.05.024. PubMed: 20688416.
  11. 11. Rudd JJ, Antoniw J, Marshall R, Motteram J, Fraaije B et al. (2010) Identification and characterisation of Mycosphaerella graminicola secreted or surface-associated proteins with variable intragenic coding repeats. Fungal Genet Biol 47: 19–32. doi:https://doi.org/10.1016/j.fgb.2009.10.009. PubMed: 19887112.
  12. 12. do Amaral AM, Antoniw J, Rudd JJ, Hammond-Kosack KE (2012) Defining the predicted protein secretome of the fungal wheat leaf pathogen Mycosphaerella graminicola. PLOS ONE 7: e49904. doi:https://doi.org/10.1371/journal.pone.0049904. PubMed: 23236356.
  13. 13. Keon J, Rudd JJ, Antoniw J, Skinner W, Hargreaves J et al. (2005) Metabolic and stress adaptation by Mycosphaerella graminicola during sporulation in its host revealed through microarray transcription profiling. Mol Plant Pathol 6: 527–540. doi:https://doi.org/10.1111/j.1364-3703.2005.00304.x. PubMed: 20565677.
  14. 14. Ray S, Anderson JM, Urmeev FI, Goodwin SB (2003) Rapid induction of a protein disulfide isomerase and defense-related genes in wheat in response to the hemibiotrophic fungal pathogen Mycosphaerella graminicola. Plant Mol Biol 53: 741–754. doi:https://doi.org/10.1023/B:PLAN.0000019120.74610.52. PubMed: 15010608.
  15. 15. Shetty NP, Jensen JD, Knudsen A, Finnie C, Geshi N et al. (2009) Effects of β-1,3-glucan from Septoria tritici on structural defence responses in wheat. J Exp Bot 60: 4287–4300. doi:https://doi.org/10.1093/jxb/erp269. PubMed: 19880540.
  16. 16. Yang F, Melo-Braga MN, Larsen MR, Jørgensen HJL, Palmisano G (2013) Battle through signaling between wheat and the fungal pathogen Septoria tritici revealed by proteomics and phosphoproteomics. Mol Cell Proteomics 12: 2497–2508. doi:https://doi.org/10.1074/mcp.M113.027532. PubMed: 23722186.
  17. 17. Wang Z, Gerstein M, Snyder M (2009) RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 10: 57–63. doi:https://doi.org/10.1038/nrg2484. PubMed: 19015660.
  18. 18. Waalwijk C, Mendes O, Verstappen EC, de Waard MA, Kema GH (2002) Isolation and characterization of the mating-type idiomorphs from the wheat septoria leaf blotch fungus Mycosphaerella graminicola. Fungal Genet Biol 35: 277–286. doi:https://doi.org/10.1006/fgbi.2001.1322. PubMed: 11929216.
  19. 19. Yang F, Jensen JD, Spliid NH, Svensson B, Jacobsen S et al. (2010) Investigation of the effect of nitrogen on severity of Fusarium head blight in barley. J Proteomics 73: 743–752. doi:https://doi.org/10.1016/j.jprot.2009.10.010. PubMed: 19895910.
  20. 20. Yang F, Jensen JD, Svensson B, Jørgensen HJL, Collinge DB et al. (2010) Analysis of early events in the interaction between Fusarium graminearum and the susceptible barley (Hordeum vulgare) cultivar Scarlett. Proteomics 10: 3748–3755. doi:https://doi.org/10.1002/pmic.201000243. PubMed: 20925056.
  21. 21. Yang F, Jensen JD, Svensson B, Jørgensen HJL, Collinge DB et al. (2012) Secretomics identifies Fusarium graminearum proteins involved in the interaction with barley and wheat. Mol Plant Pathol 13: 445–453. doi:https://doi.org/10.1111/j.1364-3703.2011.00759.x. PubMed: 22044785.
  22. 22. Wang L, Si Y, Dedow LK, Shao Y, Liu P et al. (2011) A low-cost library construction protocol and data analysis pipeline for Illumina-based strand-specific multiplex RNA-Seq. PLOS ONE 6: e26426. doi:https://doi.org/10.1371/journal.pone.0026426. PubMed: 22039485.
  23. 23. Smeds L, Künstner A (2011) ConDeTri–a content dependent read trimmer for Illumina Data. PLOS ONE 6: e26314. doi:https://doi.org/10.1371/journal.pone.0026314. PubMed: 22039460.
  24. 24. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA et al. (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29: 644–652. doi:https://doi.org/10.1038/nbt.1883. PubMed: 21572440.
  25. 25. Pertea G, Huang XQ, Liang F, Antonescu V, Sultana R et al. (2003) TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics 19: 651–652. doi:https://doi.org/10.1093/bioinformatics/btg034. PubMed: 12651724.
  26. 26. Robinson MD, Oshlack A (2010) A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol 11: R25. doi:https://doi.org/10.1186/gb-2010-11-3-r25. PubMed: 20196867.
  27. 27. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5: 621–628. doi:https://doi.org/10.1038/nmeth.1226. PubMed: 18516045.
  28. 28. Chen S, Yang PC, Jiang F, Wei YY, Ma ZY et al. (2010) De novo analysis of transcriptome dynamics in the migratory locust during the development of phase traits. PLOS ONE 5: e15633. doi:https://doi.org/10.1371/journal.pone.0015633. PubMed: 21209894.
  29. 29. Audic S, Claverie JM (1997) The significance of digital gene expression profiles. Genome Res 7: 986–995. PubMed: 9331369.
  30. 30. Vargas WA, Martín JM, Rech GE, Rivera LP, Benito EP et al. (2012) Plant defense mechanisms are activated during biotrophic and necrotrophic development of Colletotricum graminicola in maize. Plant Physiol 158: 1342–1358. doi:https://doi.org/10.1104/pp.111.190397. PubMed: 22247271.
  31. 31. Koornneef A, Pieterse CMJ (2008) Cross talk in defense signaling. Plant Physiol 146: 839–844. doi:https://doi.org/10.1104/pp.107.112029. PubMed: 18316638.
  32. 32. Mahesh SK, Liu H, Qiu D (2012) The role of radical burst in plant defense responses to necrotrophic fungi. J Integr Agricult 11: 1305–1312. doi:https://doi.org/10.1016/S2095-3119(12)60127-0.
  33. 33. Bolton MD (2009) Primary metabolism and plant defense-fuel for the fire. Mol Plant Microbe Interact 22: 487–497. doi:https://doi.org/10.1094/MPMI-22-5-0487. PubMed: 19348567.
  34. 34. Bos JIB, Kanneganti TD, Young C, Cakir C, Huitema E et al. (2006) The C-terminal half of Phytophthora infestans RXLR effector AVR3a is sufficient to trigger R3a-mediated hypersensitivity and suppress INF1- induced cell death in Nicotiana benthamiana. Plant J 48: 165–176. doi:https://doi.org/10.1111/j.1365-313X.2006.02866.x. PubMed: 16965554.
  35. 35. Dou D, Kale SD, Wang X, Chen Y, Wang Q et al. (2008) Conserved C-terminal motifs required for avirulence and suppression of cell death by Phytophthora sojae effector Avr1b. Plant Cell 20: 1118–1133. doi:https://doi.org/10.1105/tpc.107.057067. PubMed: 18390593.
  36. 36. Kelley BS, Lee SJ, Damasceno CMB, Chakravarthy S, Kim BD et al. (2010) A secreted effector protein (SNE1) from Phytophthora infestans is a broadly acting suppressor of programmed cell death. Plant J 62: 357–366. doi:https://doi.org/10.1111/j.1365-313X.2010.04160.x. PubMed: 20128886.
  37. 37. Lee SJ, Rose JKC (2010) Mediation of the transition from biotrophy to necrotrophy in hemibiotrophic plant pathogens by secreted effector proteins. Plant Signal Behav 5: 769–772. doi:https://doi.org/10.4161/psb.5.6.11778. PubMed: 20400849.
  38. 38. Qutob D, Kamoun S, Gijzen M (2002) Expression of a Phytophthora sojae necrosis-inducing protein occurs during transition from biotrophy to necrotrophy. Plant J 32: 361–373. doi:https://doi.org/10.1046/j.1365-313X.2002.01439.x. PubMed: 12410814.
  39. 39. Kirsten S, Navarro-Quezada A, Penselin D, Wenzel C, Matern A et al. (2012) Necrosis-inducing proteins of Rhynchosporium commune, effectors in quantitative disease resistance. Mol Plant Microbe Interact 25: 1314–1325. doi:https://doi.org/10.1094/MPMI-03-12-0065-R. PubMed: 22712509.
  40. 40. de Oliveira AL, Gallo M, Pazzagli L, Benedetti CE, Cappugi G et al. (2011) The structure of the elicitor cerato-platanin (CP), the first member of the CP fungal protein family, reveals a double ψβ-barrel fold and carbohydrate binding. J Biol Chem 286: 17560–17568. doi:https://doi.org/10.1074/jbc.M111.223644. PubMed: 21454637.
  41. 41. Wösten HAB (2001) Hydrophobins: Multipurpose proteins. Annu Rev Microbiol 55: 625–646. doi:https://doi.org/10.1146/annurev.micro.55.1.625. PubMed: 11544369.
  42. 42. Taylor BL, Zhulin IB (1999) PAS domains: internal sensors of oxygen, redox potential, and light. Microbiol Mol Biol Rev 63: 479–506. PubMed: 10357859.
  43. 43. Black CG, Wang L, Wu T, Coppel RL (2003) Apical location of a novel EGF-like domain-containing protein of Plasmodium falciparum. Mol Biochem Parasitol 127: 59–68. doi:https://doi.org/10.1016/S0166-6851(02)00308-0. PubMed: 12615336.
  44. 44. Mehrabi R (2006) Signaling pathways involved in pathogenicity and development of the fungal wheat pathogen Mycosphaerella graminicola. PhD thesis. , Wageningen University.
  45. 45. Shetty NP, Jørgensen HJL, Jensen JD, Collinge DB, Shetty HS (2008) Roles of reactive oxygen species in interactions between plants and pathogens. Eur J Plant Pathol 121: 267–280. doi:https://doi.org/10.1007/s10658-008-9302-5.
  46. 46. Asiimwe T, Krause K, Schlunk I, Kothe E (2012) Modulation of ethanol stress tolerance by aldehyde dehydrogenase in the mycorrhizal fungus Tricholoma vaccinum. Mycorrhiza 22: 471–484. doi:https://doi.org/10.1007/s00572-011-0424-9. PubMed: 22159964.
  47. 47. Rohel EA, Payne AC, Fraaije BA, Hollomon DH (2001) Exploring infection of wheat and carbohydrate metabolism in Mycosphaerella graminicola transformants with differentially regulated green fluorescent protein expression. Mol Plant Microbe Interact 14: 156–163. doi:https://doi.org/10.1094/MPMI.2001.14.2.156. PubMed: 11204778.