Comparative transcriptomic profiling of the two-stage response of rice to Xanthomonas oryzae pv. oryzicola interaction with two different pathogenic strains

Abstract Background Two-tiered plant immune responses involve cross-talk among defense-responsive (DR) genes involved in pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI), effector-triggered immunity (ETI) and effector-triggered susceptibility (ETS). Bacterial leaf streak (BLS), caused by Xanthomonas oryzae pv. oryzicola (Xoc) is an important bacterial disease that causes serious threats to rice yield and quality. Transcriptomic profiling provides an effective approach for the comprehensive and large-scale detection of DR genes that participate in the interactions between rice and Xoc. Results In this study, we used RNA-seq to analyze the differentially expressed genes (DEGs) in susceptible rice after inoculation with two naturally pathogenic Xoc strains, a hypervirulent strain, HGA4, and a relatively hypovirulent strain, RS105. First, bacterial growth curve and biomass quantification revealed that differential growth occurred beginning at 1 day post inoculation (dpi) and became more significant at 3 dpi. Additionally, we analyzed the DEGs at 12 h and 3 days post inoculation with two strains, representing the DR genes involved in the PTI and ETI/ETS responses, respectively. Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed on the common DEGs, which included 4380 upregulated and 4019 downregulated genes and 930 upregulated and 1383 downregulated genes identified for the two strains at 12 h post inoculation (hpi) and 3 dpi, respectively. Compared to those at 12 hpi, at 3 dpi the number of common DEGs decreased, while the degree of differential expression was intensified. In addition, more disease-related GO pathways were enriched, and more transcription activator-like effector (TALE) putative target genes were upregulated in plants inoculated with HGA4 than in those inoculated with RS105 at 3 dpi. Then, four DRs were randomly selected for the BLS resistance assay. We found that CDP3.10, LOC_Os11g03820, and OsDSR2 positively regulated rice resistance to Xoc, while OsSPX3 negatively regulated rice resistance. Conclusions By using an enrichment method for RNA-seq, we identified a group of DEGs related to the two stages of response to the Xoc strain, which included four functionally identified DR genes.