Comparative transcriptome profiling of a rice line carrying Xa39 and its parents triggered by Xanthomonas oryzae pv. oryzae provides novel insights into the broad-spectrum hypersensitive response
- Fan Zhang†1,
- Li-Yu Huang†1, 2,
- Fan Zhang1,
- Jauhar Ali3,
- Casiana Vera Cruz3,
- Da-Long Zhuo1,
- Zheng-Lin Du4,
- Zhi-kang Li1Email author and
- Yong-Li Zhou1Email author
© Zhang et al.; licensee BioMed Central. 2015
Received: 25 September 2014
Accepted: 6 February 2015
Published: 21 February 2015
Bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), is a devastating rice disease worldwide. Xa39 is a resistance (R) gene with a broad-spectrum hypersensitive response (BSHR) to Xoo. Nevertheless, the molecular mechanisms of resistance mediated by Xa39 remain unclear. In this study, the transcriptome profiling of a rice line carrying Xa39 and its parents at the early stage of Xoo infection were investigated.
A rice introgression line H471 carrying Xa39 exhibited a typical local hypersensitive response phenotype, accompanied by programmed cell death after inoculation with the Xoo Philippines’ race 9b. Transcriptome profiling of H471 and its parents at 1 and 2 days post-inoculation was performed using RNA sequencing. In total, 306 differentially expressed genes (DEGs) were identified in H471 compared with its recurrent parent Huang-Hua-Zhan after inoculation with Xoo. Among them, 121 (39.5%) genes, with functional enrichments that were related to defense response, protein amino acid phosphorylation, and apoptosis, were found to be constitutively expressed. The other 185 (60.5%) genes, with GO terms that belonged to defense response, were significantly responsive to Xoo infection in H471. Ten up-regulated and 12 down-regulated genes encoding intracellular immune receptors were identified in H471 compared with Huang-Hua-Zhan. LOC_Os11g37759, which was located in the fine-mapping region harboring Xa39, is a Xa39 candidate gene. The putative BSHR-related co-regulatory networks were constructed using 33 DEGs from four functional groups, including gibberellic acid receptors and brassinosteroid regulators, which were differentially co-expressed with LOC_Os11g37759 in infected H471. Our results indicated that there might be cross-talk between the Xa39-mediated signal transduction cascades and the GA/BR signaling pathway, and that the defense mechanism was related to diverse kinases, transcription factors, post-translational regulation, and R genes.
The present study provides the comprehensive transcriptome profile of a rice introgression line carrying Xa39 and its parents, and identifies a set of DEGs involved in BSHR mediated by Xa39. These data provide novel insights into the regulatory networks of plant disease resistance mediated by R genes, and the identified DEGs will serve as candidates for Xa39 cloning and for further understanding the molecular mechanism of BSHR.
KeywordsRice Xanthomonas oryzae pv. oryzae R gene Transcriptome sequencing
Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo), is a devastating rice disease worldwide . The development and deployment of resistant cultivars carrying major resistance (R) genes has been the most effective approach for BB management . However, rapid resistance loss in rice varieties carrying a single R gene is a problem for breeders [3-6]. Sustainable control measures for BB require a better understanding of resistance mechanisms in rice .
During the co-evolution between hosts and microorganisms, plants evolved a repertoire of R genes to defend themselves against pathogens by mounting effective, fine-tuned immune responses . Like Arabidopsis, rice has evolved a two-layered innate immune system that includes pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) and effector-triggered immunity (ETI) . PTI, the first layer of defense, is governed by pattern recognition receptors (PRRs) that recognize highly conserved PAMPs, triggering a relatively weak immune response that restricts colonization by invading organisms. To circumvent PTI, adapted pathogens can deliver effector molecules directly into the plant cell. Through co-evolution with pathogens, plants have accordingly developed intracellular immune receptors (R proteins) that can recognize the presence of certain pathogen effector molecules and trigger ETI . In contrast to PTI, ETI, the second layer of defense, is a rapid and robust response, usually associated with localized programmed cell death (PCD), referred to as the hypersensitive response (HR), which is defined as a localized and rapid cell death response at sites of pathogen attack .
Functional genomic surveys of pathogen effectors have indicated that these proteins are highly diverse in sequence, as well as in molecular function [12-14]. However, the cognate R proteins in plants are structurally conserved. The numerous R proteins identified in Arabidopsis thaliana and in rice (Oryza sativa) have typically consisted of a variable amino terminus, a nucleotide-binding site (NBS) domain in the middle, and a leucine-rich repeat (LRR) domain at the carboxyl terminus. NBS-LRR type R genes in monocots, such as rice, usually carry a putative coiled-coil (CC-NBS-LRR) domain at the N-terminus .
Serving as a model system to elucidate the interactions between pathogens and monocotyledon plants, 39 R genes, 28 dominant and 11 recessive, conferring resistance to BB have been registered (http://www.shigen.nig.ac.jp/rice/oryzabase/gene/list) and identified , and six genes, Xa1, xa5, xa13, Xa21, Xa26/Xa3, and Xa27, have been cloned (http://www.shigen.nig.ac.jp/rice/oryzabase/gene/list). The first class of R genes contains receptor kinases, including the cloned Xa21 and Xa26 [17,18], but the molecules produced by Xoo that are recognized by XA21 and XA26 have not yet been characterized. Xa1 represents the second major class of R genes, the NBS-LRR group. The gene for the elicitor signal from the pathogen has not been identified for Xa1. However, some genes encode diverse proteins that are expressed and trigger ETI after the recognition of the pathogen-delivered transcription activator-like (TAL) effector genes avrXa27 [19-21]. Recently, there have been great advances in understanding the recognition between the R protein in rice and effector in Xoo including Xa27 and AvrXa27 [22-24], Xa10 and AvrXa10 (EBEAvrXa10) , and Xa23 and TalC9b/AvrXa23  however, research on how the R proteins trigger signal transduction cascades leading to HR in rice is still obscure.
Xa39 is a novel dominant R gene with a broad-spectrum hypersensitive response (BSHR) to Xoo . In this study, we compared transcriptome profiling by RNA sequencing (RNA-Seq) and identified differentially expressed genes (DEGs) in a rice introgression line (IL) H471 compared with its parents. Our results elucidated some interesting molecular mechanisms underlying BSHR mediated by R proteins in rice.
Phenotype and DNA ladder detection of H471-Xa39 in response to Xoo
DNA fragmentation rates of HHZ and H471 were detected after inoculation. DNA fragmentation was observed in both HHZ and H471 at 1 dpi (Figure 1C), and more fragmented DNA was observed in H471 than HHZ at 3 dpi. Notably, the DNA in HHZ was degraded at 5 dpi because of necrosis, whereas the DNA in H471 showed no changes, owing to the cessation of cell death.
Gene expression profiling of three rice genotypes under inoculation and non-inoculation conditions
Mapping results of RNA sequencing reads of the rice introgression line H471 carrying Xa39 and parents Huang-Hua-Zhan (HHZ) and PSBRC28 (P28) at 1 and 2 days post-inoculation (dpi) with Xanthomonas oryzae pv. oryzae and under non-inoculation conditions (CK)
Biological replication I
Biological replication II
Total mapped unique genes
Total filtered pair- end reads
Total mapped reads (%)
Total filtered pair- end reads
Total mapped reads (%)
2 × 10,068,639
2 × 10,083,264
2 × 8,894,670
2 × 9,326,860
2 × 10,384,097
2 × 10,326,754
2 × 5,232,832
2 × 11,822,492
2 × 10,296,887
2 × 6,174,204
2 × 12,505,956
2 × 11,484,590
2 × 10,902,112
2 × 11,010,814
2 × 5,873,150
2 × 5,453,981
2 × 7,462,058
2 × 10,400,231
The expressed genes in all samples were subjected to a cluster analysis. As shown in Additional file 2A: Figure S2, the three genotypes at 1 and 2 dpi, as well as the control, were separated from each other. The three control samples clustered into a group, and the inoculated samples clustered into another group, suggesting that most of the expressed genes had similar expression patterns in response to PXO349 infection, even in different genotypes. The inoculated H471 and HHZ clustered into a group, consistent with their similar genetic backgrounds, and samples at 1 and 2 dpi were in different subgroups, indicating that the gene expression responses to pathogen infection were dynamic.
To investigate the intrinsic differences in gene expression between the resistant and susceptible genotypes, the expression levels of genes in H471 were compared with those in HHZ and P28 under the non-inoculation condition. There were 1,161 DEGs identified between H471 and P28, and 343 DEGs between H471 and HHZ (Additional file 2B: Figure S2 and Additional file 3: Table S1). Between HHZ and P28, there were 1,286 DEGs because of the larger genetic difference between these two lines than between H471 and its parents. The genes with higher basal expression levels in H471 compared with HHZ were mainly functionally enriched in oxidoreductase activity, lyase activity, carboxylic acid metabolic process, response to stress, and cofactor binding (Additional file 4: Table S2).
DEGs between the resistant and susceptible genotypes after inoculation with Xoo
These 306 DEGs could be classified into two groups based on their expression patterns before and after the Xoo inoculation of H471. The first group contained 121 (39.5%) DEGs, which were differentially expressed under both non-inoculation and inoculation conditions. Among them, 64 genes were differentially up-regulated in H471 compared with HHZ after inoculation, whereas 57 genes were differentially down-regulated after inoculation (Additional file 5: Table S3). The first group is functionally enriched in defense response, protein amino acid phosphorylation, and apoptosis-associated proteins by GO analysis (Additional file 6: Table S4). The second group contained 185 DEGs that were responsive to Xoo infection at a minimum of one time point after inoculation. Among them, only one significant GO term (GO: 0006952, defense response), characterizing nine DEGs, was identified. These nine disease R genes included three up-regulated and six down-regulated genes responsive to Xoo infection (Additional file 5: Table S3).
DEGs related to protein kinases in H471 and HHZ respond to Xoo infection
Protein kinases play central roles in signal recognition and the subsequent activation of plant defense mechanisms during pathogen infection. In this study, 31 genes encoding protein kinases were differentially expressed in H471 compared with HHZ, and 18 were responsive to Xoo infection (Additional file 7: Table S5). Among these DEGs, three types of receptor-like protein kinases (RLKs) were identified. The first group was composed of five genes encoding cysteine-rich receptor kinases (CRKs), which were differentially expressed in H471 compared with HHZ at 1 dpi (Additional file 7: Table S5). Among them, two genes encoding cysteine-rich RLK 10 (LOC_Os07g43560 and LOC_Os07g43570) were up-regulated, whereas the other three genes (LOC_Os04g30030, LOC_Os11g28104, and LOC_Os04g30040) were down-regulated. The second group was composed of two genes encoding the LRR protein kinases (LRKs) (LOC_Os08g14950 and LOC_Os11g36190), which were up-regulated in H471 compared with HHZ after inoculation with Xoo, and six genes encoding LRKs (LOC_Os03g56270, LOC_Os11g35500, LOC_Os08g14940, LOC_Os11g36140, LOC_Os11g36150, and LOC_Os03g28270), which were down-regulated in H471 compared with HHZ after Xoo inoculation (Additional file 7: Table S5). Of these, the expression of LOC_Os08g14950 was induced by Xoo infection, while the expression levels of the other four genes (LOC_Os08g14940, LOC_Os11g36140, LOC_Os11g36150, and LOC_Os03g28270) were inhibited in H471 at 1 dpi. The final group of RLKs consisted of seven genes encoding wall-associated kinases (WAKs), which were differentially expressed in H471 compared with HHZ at 1 dpi. Of the differentially expressed WAKs, OsWAK38 (LOC_Os04g29680) was up-regulated in H471 compared with HHZ at 1 dpi, while the other six genes (LOC_Os04g30160, LOC_Os04g30110, LOC_Os04g29770, LOC_Os04g30250, LOC_Os04g29790, and LOC_Os04g29810) were down-regulated. These genes were all located in the introgressed fragments of chromosome 4 in H471 (Additional file 7: Table S5).
In addition to RLKs, two genes encoding calcium/calmodulin-dependent protein kinases (CPKs) were differentially expressed in H471 at 1 dpi. One gene (LOC_Os03g50330) was up-regulated, and the other (LOC_Os03g43440) was down-regulated (Additional file 7: Table S5).
DEGs involved in phytohormone signaling pathways in response to Xoo infection
Comparisons of the expression levels of representative hormone genes involved in rice innate immunity between the rice introgression line H471 carrying Xa39 and the recurrent parent Huang-Hua-Zhan (HHZ)
H471- CK vs HHZ- CK
H471- 1dpi vs HHZ- 1dpi
H471-2dpi vs HHZ-2dpi
H471- 1dpi vs H471- CK
H471- 2dpi vs H471- CK
SA pathway related gene
SA pathway related gene
JA responsive gene
ABA/JA related gene
JA/Ethylene related gene
Ethylene related gene
Differential expression of transcription factors (TFs) and genes possibly related to post-transcription regulation in H471 and HHZ
Differential expression of transcription process components in the rice introgression line H471 compared with the recurrent parent Huang-Hua-Zhan (HHZ) after inoculation with Xanthomonas oryzae pv. oryzae
H471- 1dpi vs HHZ- 1dpi
H471- 2dpi vs HHZ- 2dpi
H471- 1dpi vs H471- CK
H471- 2dpi vs H471- CK
HRT-like/transcription repressor HOTR
CO-like/OsCO3/CCT/B-box zinc finger protein
CO-like/CCT/B-box zinc finger protein
PPR repeat domain containing protein
Pentatricopeptide repeat (PPR) genes target effectors to the correct site on the correct transcripts, and are thus involved in many post-transcriptional processes . In this study, three genes encoding PPRs were significantly expressed in infected H471 compared with HHZ (Table 3).
Defense-related DEGs in the resistant and susceptible genotypes
The 10 R genes up-regulated in H471 compared with HHZ included four genes belonging to the NB-ARC domain-containing protein family (LOC_Os03g40194, LOC_Os04g30930, LOC_Os08g14810 and LOC_Os11g39320), three genes encoding disease resistance proteins containing LRR/NBS-LRR domains (LOC_Os11g35850, LOC_Os05g33520 and LOC_Os11g37759), and three genes encoding pathogenesis-related proteins (PRs) (LOC_Os12g36880, LOC_Os12g36830 and LOC_Os12g36850) (Figure 3 and Additional file 10: Table S7). A further analysis indicated that four (LOC_Os04g30930, LOC_Os08g14810, LOC_Os11g39320, and LOC_Os11g37759) of the up-regulated R genes in H471 were co-localized in the introgressed regions (Figure 3 and Additional file 10: Table S7), implying roles in Xoo resistance.
PRs, which are involved in plant resistance to pathogens, are commonly induced in plants resistant to pathogens of viral, fungal, and bacterial origin . For example, PR10/Bet v1, which is induced in hot pepper (Capsicum annuum) by incompatible interactions with Xanthomonas campestris pv. vesicatoria, may function as a ribonuclease. Subsequent phosphorylation of CaPR10/Bet v1 increases its ribonucleolytic activity to cleave invading viral RNAs, and this activity is important to its in vivo antiviral pathway . Besides the three up-regulated genes encoding PRs, the expression of three genes (LOC_Os12g36830, LOC_Os12g36850, and LOC_Os12g36880) encoding the pathogenesis-related Bet v I family protein were inhibited in Xoo-infected H471.
Notably, two NBS-LRR class genes (LOC_Os05g33520 and LOC_Os11g37759) were induced in H471 after inoculation with Xoo (Additional file 8: Figure S3). LOC_Os05g33520 encoded a homolog resistant to phytophthora 1 in Arabidopsis, which positively regulates the Phytophthora brassicae-induced oxidative burst and enhances expression of defense-related genes . LOC_Os11g37759, encoding a CC-NBS-LRR disease resistance protein, was expressed at a higher level in H471 compared with its parents before Xoo inoculation and was only induced by Xoo infection in H471 (Additional file 8: Figure S3). Notably, LOC_Os11g37759 was the only DEG out of 10 up-regulated R genes located in the 97.4-kb region (from 22233214 to 22330619 bp) on chromosome 11 harboring the Xa39 gene based on our previous fine mapping , suggesting it was the most likely candidate gene of Xa39.
Putative immunity mechanism of broad-spectrum resistance mediated by Xa39
The BSHR phenotype of rice IL carrying Xa39 is different from that of its parents. It is possible that a rare unequal rearrangement event, insertion, or deletion has occurred, which causes the novel resistance . In our previous study, two genes encoding CC-NBS-LRR disease resistance proteins, LOC_Os11g37740 and LOC_Os11g37759, were mapped in the 97.4-kb region harboring the Xa39 gene on rice chromosome 11 . In this study, LOC_Os11g37759 was differentially up-regulated in H471 compared with its parents, suggesting that the hypothetical mechanism of Xa39 might involve this gene.
Among 39 R genes resistant to BB, Xa39 and Xa23 exhibit similar BSHR, with LLs shorter than 1 cm [16,35]. Xa21, which is discussed as a putative PRR, and the signaling network mediated by Xa21 have been studied extensively . Xa21-containing rice lines also have broader resistance spectra, but the LLs were longer than 5 mm when inoculated with Xoo strains possessing strong virulence, including PXO349. This is in accordance with PTI playing a pivotal role in the defense against a broad spectrum of potential pathogens, but triggering a relatively weak immune response in contrast to ETI [9,19]. To date, it has been reported that Xa23 can be recognized by the TAL effector talC9b/avrXa23 and triggers resistance to Xoo ; however, the signaling network leading to HR still remains to be elucidated. In this study, the comparative transcriptome profiling of H471 and its parents provided several interesting insights into the molecular mechanisms of BSHR in rice.
First, the resistance mediated by Xa39 might not be related to ion fluxes and oxidative bursts. The signal-specific activation of plant PRRs by PAMPs leads to a series of cellular processes, including ion fluxes, oxidative bursts by the production of reactive oxygen species (ROS), activation of downstream MAPK cascades, transcriptional changes, and the production of antimicrobial compounds, such as PR and phytoalexins . In this study, six genes encoding PRs were differentially expressed; however, there were no significant differences in the expression levels of the representative PRRs of PTI, genes related to oxidative burst, or downstream genes related to cell death due to ion fluxes and ROS in H471 compared with HHZ after inoculation with Xoo (Additional file 8: Figure S3 and Additional file 12: Table S9). This suggests that hypersensitive cell death at sites of Xoo infection did not result from the signal-specific activation of PRRs in H471.
Second, GA and BR signaling pathways may be involved in BSHR mediated by Xa39. It is generally accepted that SA plays a major role in activating defenses against biotrophic pathogens, whereas JA is usually associated with defenses against necrotrophic pathogen attacks . Grewal et al.  reported that JA and ethylene signaling were involved in the early defense of rice cells, and four overlapping genes related to JA or ethylene were up-regulated or down-regulated in H471 after inoculation with Xoo (Table 2). However, none of them showed a significantly different expression level between H471 and HHZ after Xoo inoculation (Table 2). Interestingly, several genes related to GA and BR signaling were differentially expressed in H471 compared with its parents. In rice, it was inferred that the GA regulator OsSERK2 positively regulates immunity mediated by XA21 and XA3 [19,39]. LOC_Os11g37759, encoding a CC-NBS-LRR protein, is structurally related to XA21 and XA3, and OsSERK1/OsBAK1 is the closest paralog to OsSERK2 in rice. Interestingly, we found that there was a strong co-expression of LOC_Os11g37759 and OsSERK1/OsBAK1 (Figure 4). In rice, SLR1 (DELLA family), a key regulator of GA signaling, functions in disease resistance to Xoo as a positive regulator through cross-talk with the JA signaling pathway via the GA receptor GID1 . In this study, there were no significant differences in the expression levels of SLR1 and GID1 between H471 and HHZ; however, the GA receptor GID1L2 was differentially co-expressed with LOC_Os11g37759. The BR signaling-related gene OsBAK1 (LOC_Os08g07760) and BR biosynthesis-related gene LOC_Os03g40540 were expressed at higher levels in H471 compared with HHZ after Xoo infection (Table 2). The activation of BR signaling inhibits PTI by the BAK1-independent recognition of the fungal PAMP chitin . Although the type of immunity activated by Xa39 remains a mystery, the differential expression levels of these hormone regulators suggests that there might be cross-talk between the Xa39-mediated signal transduction cascades and the GA/BR signaling pathway.
Finally, diverse categories of DEGs related to signal transduction and regulation were identified. WAKs span the plasma membrane, allowing cells to recognize and respond to their extracellular environment . CPKs can regulate downstream components of complex immune and stress signaling networks as positive or negative regulators . For example, OsCPK12 is involved in both salt-stress tolerance and blast disease resistance in rice . Several CRKs are reportedly associated with resistance to pathogen infection and activation of PCD in Arabidopsis [44-46], barley , and wheat . In this study, DEGs encoding LRKs, WAKs, CPKs, and CRKs were identified at 1 and 2 dpi, indicating that a series of signaling events was regulated by complex communications at the early stage of infection in H471. TFs in several families such as WRKY, MADs box, and NAC are involved in biotic and abiotic stresses, as well as pathogen invasion . In this study, besides the two TFs in the WRKY family, TFs in B3, bHLH/OsIRO2, M-type, HRT-like, and CO-like families were differentially expressed in H471 compared with their parents, suggesting that functions of diverse TFs in plant resistance remain to be elucidated.
HR accompanying PCD at the site of attack is an effective defense strategy against pathogens and nematodes that feed on live plant cells . However, unlike the nucleotide binding oligomerization domain (NOD)-like receptors (NLRs) and LRR-containing proteins in animals, which are known to trigger cell apoptosis through the activity of caspases that activate proinflammatory cytokines , plants lack the homologous caspases. The hypothesis that R-protein-mediated PCD might involve the perturbation of multiple cellular processes came from a report showing that in A. thaliana, the resistance conferred by the R protein RPP4 against the obligate biotroph Hyaloperonospora arabidopsidis is not mediated by a single gene, but rather by multiple downstream genes that encode proteins (mostly enzymes with very diverse functions) . The DEGs identified in this study, which might be involved in perception, signal transduction, transcriptional regulation, and defense-related genes, will be studied to determine their functions and to understand the molecular networks involved in the BSHR mediated by R genes in rice.
The BB R gene Xa39, which is able to trigger a resistance response to Xoo characterized by PCD in rice plants, activates DEGs encoding proteins with diverse functions leading to BRHR. These results provide novel insights into the regulatory networks of plant resistance mediated by R proteins.
Plant materials and artificial inoculation
Three rice genotypes were used in this study. HHZ, with a high yield and good quality, is a widely used indica inbred rice line in southern China. P28 is an indica rice variety from the Philippines. H471, carrying Xa39, is a BC1F6 IL with a few chromosomal fragments introgressed from the donor parent P28 into the background of the recurrent parent HHZ. A genome-wide single nucleotide polymorphism analysis by re-sequencing showed that H471 differs from HHZ at 26 genomic segments, which originate from P28, with sizes in the range of 57–6,057 kb (Additional file 13: Figure S4). The recombination events were judged as described in Huang et al. , using the genome-resequencing data of HHZ, P28, and H471 (unpublished data). Seeds of H471, HHZ, and P28 were sown in a seedling nursery and 30-day-old seedlings were transplanted into a screenhouse of the Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, China. There were nine plants in each row, with a spacing of 20 × 17 cm. The Philippine’s representative strain of Xoo, PXO349, was used to artificially inoculate H471 and its parents. The strain was incubated on peptone sucrose agar at 30°C for 2 days, and the inoculum was prepared by suspending the bacterial mass in sterile water at a concentration of 108 cells mL−1. Five plants of each line were inoculated with PXO349 in the four to five uppermost leaves of each plant using the leaf-clipping method  in three replications at the tillering stage (plant age of 65 days). LLs were measured on all inoculated leaves at 14 dpi when the lesions were stable.
DNA fragmentation of HHZ and H471 leaves infected by PXO349
Leaves of HHZ and H471 inoculated by PXO349 at 1, 3, and 5 dpi were collected for a DNA fragmentation analysis. A 5-mm leaf fragment under the incision was collected and extracted using the DNeasy Plant Mini Kit (Cat. No.69106) following the manufacturer’s manual.
Preparation of samples for RNA-Seq
Leaves of plants were clipped with sterile water and PXO349, respectively. 1-cm-long leaf tips from H471, HHZ, and P28 were dissected at 0 day after treatment with water, and at 1 day and 2 days after inoculation with PXO349. Two replicate leaf samples were taken for each rice line treatment at each time point. Samples of H471, HHZ, and P28 collected at 0 dpi were used as control and named as H471-CK, HHZ-CK, and P28-CK, respectively. Samples of H471, HHZ, and P28 collected at 1 dpi were named as H471-1dpi, HHZ-1dpi, and P28-1dpi, respectively; and samples collected at 2 dpi were named as H471-2dpi, HHZ-2dpi, and P28-2dpi, respectively. All samples were immediately frozen in liquid nitrogen after collection, and stored at −80°C.
RNA extraction, mRNA-Seq library construction, and sequencing
For all RNA-Seq samples, total RNA was extracted with TRIzol Reagent (Invitrogen, Carlsbad, CA, USA), and quantified using a Qubit RNA assay kit (Applied Biosystems, Foster City, CA, USA). RNA integrity was checked using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Total mRNAs were isolated by oligo(dT) selection using Dynal® magnetic beads (Invitrogen). The paired-end fragment library was constructed following the TruSeq RNA Sample Preparation kit (Illumina, San Diego, CA, USA) protocol with minor modifications. Cluster generation of the produced libraries was performed using cBot, and sequenced on a HiSeq 2000 platform (Illumina) with paired-end 100-bp reads. Primary data analysis and base calling were performed using the Illumina instrument software.
Transcriptome data processing and analysis
All primary sequencing data are available at the GEO database under the accession number GSE62488 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE62488). Low-quality nucleotides (< Q20) were trimmed from raw sequences for each sample, and then pair-end reads with lengths less than 30 bp were removed using an in-house Perl script. Retained high-quality pair-end reads of rice for each sample were mapped to the rice genome of the Rice Genome Annotation Project (RGAP) at MSU  using TopHat , and assembled with Cufflinks  to construct unique transcript sequences using the parameter –g –b –u –o. Cuffcompare was used to compare the assembled transfrags of each sample to the reference annotation, and to build a non-redundant transcript dataset among the samples . The number of mapped clean reads for each gene was counted and normalized into the reads per kilobase per million value . Next, Cuffdiff was used to identify DEGs, and genes with P values ≤ 0.001 were marked as significantly different between the two samples . GO terms of rice DEGs were identified in accordance to the methodology described by Du et al. . A Venn diagram was built using software available online (http://bioinfogp.cnb.csic.es/tools/venny/index.html) .
Confirmation of the DEGs by quantitative real-time RT-PCR (qRT-PCR)
To validate the Illumina sequencing results, a subset of DEGs were verified by qRT-PCR (n = 108). An independent set of samples was collected at 6, 12, 24, 36, 48, and 60 h after inoculation with PXO349 to use for the expression analysis of some important DEGs. The sequence for each rice gene was obtained from the MSU rice database . Sequences from each gene were used to design primers using Primer 5 software (http://frodo.wi.mit.edu/; Additional file 14: Table S10). Independent biological repetitions of each experiment were performed in triplicate. Expression levels of 17 rice genes were tested in 20-μL reactions using the SYBR® Green PCR Master Mix kit (Applied Biosystems) following the manufacturer’s protocol via an ABI Prism 7900 Sequence Detection System (Applied Biosystems).
Availability of supporting data
The sequencing raw data of this article are available through a GEO database at the NCBI (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE62488). And the data sets supporting the results of this article are included within the article and its additional files.
This research was supported by the National Natural Science Foundation of China (Grant Nos. 31161140349 and U1201211), the National High-tech Program of China (No. 2014AA10A603), and the CAAS Innovative Team Award.
- Mew TW. Current status and future prospects of research on bacterial blight of rice. Annu Rev Phytopathol. 1987;25:359–82.View ArticleGoogle Scholar
- Ogawa T. Methods and strategy for monitoring race distribution and identification of resistance genes to bacterial leaf blight (Xanthomonas campetris pv. oryzae) in rice. Jarq-Jpn Agr Res Q. 1993;27:71–80.Google Scholar
- Lee SW, Choi SH, Han SS, Lee DG, Lee BY. Distribution of Xanthomonas oryzae pv. oryzae Strains Virulent to Xa21 in Korea. Phytopathology. 1999;89(10):928–33.View ArticlePubMedGoogle Scholar
- Marella LS, George MLC, Vera Cruz CM, Bernardo MA, Nelson RJ, Leung H. Identification of resistance genes effective against rice bacterial blight pathogen in eastern India. Plant Disease. 2001;85(5):506–12.View ArticleGoogle Scholar
- Zeng LX, Huang SH, Wu SZ. The resistance of IRBB21 (Xa21) against 5 races of Guangdong province. Acta Phytopathologica Sinica. 2002;29(2):97–100.Google Scholar
- Zhang Q. Highlights in identification and application of resistance genes to bacterial blight. Chinese Journal of Rice Science. 2005;19(5):453–9.Google Scholar
- Nino-Liu DO, Ronald PC, Bogdanove AJ. Xanthomonas oryzae pathovars: model pathogens of a model crop. Mol Plant Pathol. 2006;7(5):303–24.View ArticlePubMedGoogle Scholar
- Bai S, Liu J, Chang C, Zhang L, Maekawa T, Wang Q, et al. Structure-function analysis of barley NLR immune receptor MLA10 reveals its cell compartment specific activity in cell death and disease resistance. Plos Pathog. 2012;8(6):e1002752.View ArticlePubMed CentralPubMedGoogle Scholar
- Jones JD, Dangl JL. The plant immune system. Nature. 2006;444(7117):323–9.View ArticlePubMedGoogle Scholar
- Boller T, He SY. Innate immunity in plants: an arms race between pattern recognition receptors in plants and effectors in microbial pathogens. Science. 2009;324(5928):742–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Tsuda K, Katagiri F. Comparing signaling mechanisms engaged in pattern-triggered and effector-triggered immunity. Curr Opin Plant Biol. 2010;13(4):459–65.View ArticlePubMedGoogle Scholar
- Petnicki-Ocwieja T, Schneider DJ, Tam VC, Chancey ST, Shan L, Jamir Y, et al. Genomewide identification of proteins secreted by the Hrp type III protein secretion system of Pseudomonas syringae pv. tomato DC3000. Proc Natl Acad Sci U S A. 2002;99(11):7652–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Salzberg SL, Sommer DD, Schatz MC, Phillippy AM, Rabinowicz PD, Tsuge S, et al. Genome sequence and rapid evolution of the rice pathogen Xanthomonas oryzae pv. oryzae PXO99A. BMC Genomics. 2008;9:204.View ArticlePubMed CentralPubMedGoogle Scholar
- Cabral A, Stassen JH, Seidl MF, Bautor J, Parker JE, Van den Ackerveken G. Identification of Hyaloperonospora arabidopsidis transcript sequences expressed during infection reveals isolate-specific effectors. PLoS One. 2011;6(5):e19328.View ArticlePubMed CentralPubMedGoogle Scholar
- Pan Q, Wendel J, Fluhr R. Divergent evolution of plant NBS-LRR resistance gene homologues in dicot and cereal genomes. J Mol Evol. 2000;50(3):203–13.PubMedGoogle Scholar
- Zhang F, Zhuo DL, Huang LY, Wang WS, Xu JL, Vera Cruz C, et al. Xa39, a novel dominant gene conferring broad-spectrum resistance to Xanthomonas oryzae pv. oryzae in rice. Plant Pathol 2014; doi:10.1111/ppa.12283.Google Scholar
- Song WY, Wang GL, Chen LL, Kim HS, Pi LY, Holsten T, et al. A receptor kinase-like protein encoded by the rice disease resistance gene, Xa21. Science. 1995;270(5243):1804–6.View ArticlePubMedGoogle Scholar
- Sun X, Cao Y, Yang Z, Xu C, Li X, Wang S, et al. Xa26, a gene conferring resistance to Xanthomonas oryzae pv. oryzae in rice, encodes an LRR receptor kinase-like protein. Plant J. 2004;37(4):517–27.View ArticlePubMedGoogle Scholar
- Liu W, Liu J, Triplett L, Leach JE, Wang GL. Novel insights into rice innate immunity against bacterial and fungal pathogens. Annu Rev Phytopathol. 2014;52:213–41.View ArticlePubMedGoogle Scholar
- Wang C, Zhang X, Fan Y, Gao Y, Zhu Q, Zheng C, et al. XA23 is an executor R protein and confers broad-spectrum disease resistance in rice. Mol Plant 2014 Nov 9. pii: ssu132. [Epub ahead of print].Google Scholar
- Tian D, Wang J, Zeng X, Gu K, Qiu C, Yang X, et al. The rice TAL effector-dependent resistance protein XA10 triggers cell death and calcium depletion in the endoplasmic reticulum. Plant Cell. 2014;26(1):497–515.View ArticlePubMed CentralPubMedGoogle Scholar
- Gu K, Yang B, Tian D, Wu L, Wang D, Sreekala C, et al. R gene expression induced by a type-III effector triggers disease resistance in rice. Nature. 2005;435(7045):1122–5.View ArticlePubMedGoogle Scholar
- Yang B, Sugio A, White FF. Os8N3 is a host disease-susceptibility gene for bacterial blight of rice. Proc Natl Acad Sci U S A. 2006;103(27):10503–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Gu K, Tian D, Qiu C, Yin Z. Transcription activator-like type III effector AvrXa27 depends on OsTFIIAγ5 for the activation of Xa27 transcription in rice that triggers disease resistance to Xanthomonas oryzae pv. oryzae. Mol Plant Pathol. 2009;10(6):829–35.View ArticlePubMedGoogle Scholar
- Abbruscato P, Nepusz T, Mizzi L, Del Corvo M, Morandini P, Fumasoni I, et al. OsWRKY22, a monocot WRKY gene, plays a role in the resistance response to blast. Mol Plant Pathol. 2012;13(8):828–41.View ArticlePubMedGoogle Scholar
- Liu X, Bai X, Wang X, Chu C. OsWRKY71, a rice transcription factor, is involved in rice defense response. J Plant Physiol. 2007;164(8):969–79.View ArticlePubMedGoogle Scholar
- Peng Y, Bartley LE, Chen X, Dardick C, Chern M, Ruan R, et al. OsWRKY62 is a negative regulator of basal and Xa21-mediated defense against Xanthomonas oryzae pv. oryzae in rice. Mol Plant. 2008;1(3):446–58.View ArticlePubMedGoogle Scholar
- Barkan A, Small I. Pentatricopeptide repeat proteins in plants. Annu Rev Plant Biol. 2014;65:415–42.View ArticlePubMedGoogle Scholar
- Edreva A. Pathogenesis-related proteins: research progress in the last 15 years. Gen Appl Plant Physiology. 2005;31:105–24.Google Scholar
- Park CJ, Kim KJ, Shin R, Park JM, Shin YC, Paek KH. Pathogenesis-related protein 10 isolated from hot pepper functions as a ribonuclease in an antiviral pathway. Plant J. 2004;37(2):186–98.View ArticlePubMedGoogle Scholar
- Belhaj K, Lin B, Mauch F. The chloroplast protein RPH1 plays a role in the immune response of Arabidopsis to Phytophthora brassicae. Plant J. 2009;58(2):287–98.View ArticlePubMedGoogle Scholar
- Cao P, Jung KH, Choi D, Hwang D, Zhu J, Ronald PC. The Rice Oligonucleotide Array Database: an atlas of rice gene expression. Rice. 2012;5(1):17.View ArticlePubMedGoogle Scholar
- Chinchilla D, Zipfel C, Robatzek S, Kemmerling B, Nurnberger T, Jones JD, et al. A flagellin-induced complex of the receptor FLS2 and BAK1 initiates plant defence. Nature. 2007;448(7152):497–500.View ArticlePubMedGoogle Scholar
- Hayashi N, Inoue H, Kato T, Funao T, Shirota M, Shimizu T, et al. Durable panicle blast-resistance gene Pb1 encodes an atypical CC-NBS-LRR protein and was generated by acquiring a promoter through local genome duplication. Plant J. 2010;64(3):498–510.View ArticlePubMedGoogle Scholar
- Zhang Q, Lin SC, Zhao BY, Wang CL, Yang WC, Zhou YL, et al. Identification and tagging a new gene for resistance to bacterial blight (Xanthomonas oryzae pv. oryzae) from O. rufipogon. Rice Genetics Newsletter. 1998;15:138–42.Google Scholar
- Spoel SH, Dong X. How do plants achieve immunity? Defence without specialized immune cells. Nat Rev Immunol. 2012;12(2):89–100.View ArticlePubMedGoogle Scholar
- Adie BA, Perez-Perez J, Perez-Perez MM, Godoy M, Sanchez-Serrano JJ, Schmelz EA, et al. ABA is an essential signal for plant resistance to pathogens affecting JA biosynthesis and the activation of defenses in Arabidopsis. Plant Cell. 2007;19(5):1665–81.View ArticlePubMed CentralPubMedGoogle Scholar
- Grewal RK, Gupta S, Das S. Xanthomonas oryzae pv oryzae triggers immediate transcriptomic modulations in rice. BMC Genomics. 2012;13:49.View ArticlePubMed CentralPubMedGoogle Scholar
- Chen X, Zuo S, Schwessinger B, Chern M, Canlas PE, Ruan D, et al. An XA21-associated kinase (OsSERK2) regulates immunity mediated by the XA21 and XA3 immune receptors. Mol Plant. 2014;7(5):874–92.View ArticlePubMed CentralPubMedGoogle Scholar
- Albrecht C, Boutrot F, Segonzac C, Schwessinger B, Gimenez-Ibanez S, Chinchilla D, et al. Brassinosteroids inhibit pathogen-associated molecular pattern-triggered immune signaling independent of the receptor kinase BAK1. Proc Natl Acad Sci U S A. 2012;109(1):303–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Fritz-Laylin LK, Krishnamurthy N, Tor M, Sjolander KV, Jones JD. Phylogenomic analysis of the receptor-like proteins of rice and Arabidopsis. Plant Physiol. 2005;138(2):611–23.View ArticlePubMed CentralPubMedGoogle Scholar
- Asano T, Tanaka N, Yang G, Hayashi N, Komatsu S. Genome-wide identification of the rice calcium-dependent protein kinase and its closely related kinase gene families: comprehensive analysis of the CDPKs gene family in rice. Plant Cell Physiol. 2005;46(2):356–66.View ArticlePubMedGoogle Scholar
- Asano T, Hayashi N, Kobayashi M, Aoki N, Miyao A, Mitsuhara I, et al. A rice calcium-dependent protein kinase OsCPK12 oppositely modulates salt-stress tolerance and blast disease resistance. Plant J. 2012;69(1):26–36.View ArticlePubMedGoogle Scholar
- Zhang X, Han X, Shi R, Yang G, Qi L, Wang R, et al. Arabidopsis cysteine-rich receptor-like kinase 45 positively regulates disease resistance to Pseudomonas syringae. Plant physiology and biochemistry: PPB/Societe francaise de physiologie vegetale. 2013;73:383–91.View ArticleGoogle Scholar
- Chen K, Fan B, Du L, Chen Z. Activation of hypersensitive cell death by pathogen-induced receptor-like protein kinases from Arabidopsis. Plant Mol Biol. 2004;56(2):271–83.View ArticlePubMedGoogle Scholar
- Ederli L, Madeo L, Calderini O, Gehring C, Moretti C, Buonaurio R, et al. The Arabidopsis thaliana cysteine-rich receptor-like kinase CRK20 modulates host responses to Pseudomonas syringae pv. tomato DC3000 infection. J Plant Physiol. 2011;168(15):1784–94.View ArticlePubMedGoogle Scholar
- Rayapuram C, Jensen MK, Maiser F, Shanir JV, Hornshoj H, Rung JH, et al. Regulation of basal resistance by a powdery mildew-induced cysteine-rich receptor-like protein kinase in barley. Mol Plant Pathol. 2012;13(2):135–47.View ArticlePubMedGoogle Scholar
- Yang K, Rong W, Qi L, Li J, Wei X, Zhang Z. Isolation and characterization of a novel wheat cysteine-rich receptor-like kinase gene induced by Rhizoctonia cerealis. Sci Rep. 2013;3:3021.PubMed CentralPubMedGoogle Scholar
- Ting JP, Willingham SB, Bergstralh DT. NLRs at the intersection of cell death and immunity. Nat Rev Immunol. 2008;8(5):372–9.View ArticlePubMedGoogle Scholar
- Wang W, Barnaby JY, Tada Y, Li H, Tor M, Caldelari D, et al. Timing of plant immune responses by a central circadian regulator. Nature. 2011;470(7332):110–4.View ArticlePubMedGoogle Scholar
- Huang X, Feng Q, Qian Q, Zhao Q, Wang L, Wang A, et al. High-throughput genotyping by whole-genome resequencing. Genome Res. 2009;19(6):1068–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Kauffman HE, Reddy APK, Hsieh SPY, Merca SD. A improved technique for evaluation of resistance of rice varieties to Xanthomonas oryzea. Plant Dis Rep. 1973;57:537–41.Google Scholar
- Kawahara Y, de la Bastide M, Hamilton JP, Kanamori H, McCombie WR, Ouyang S, et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013;6(1):4.View ArticlePubMedGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.View ArticlePubMed CentralPubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.View ArticlePubMedGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38(Web Server issue):W64–70.View ArticlePubMed CentralPubMedGoogle Scholar
- Oliveros JC. VENNY. An interactive tool for comparing lists with Venn Diagrams. 2007. http://bioinfogp.cnb.csic.es/tools/venny/index.html.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.