Identification and characterization of genes frequently responsive to Xanthomonas oryzae pv. oryzae and Magnaporthe oryzae infections in rice

Disease resistance is an important factor that impacts rice production. However, the mechanisms underlying rice disease resistance remain to be elucidated. Here, we show that a robust set of genes has been defined in rice response to the infections of Xanthomonas oryzae pv. oryzae (Xoo) and Magnaporthe oryzae (Mor). We conducted a comprehensive analysis of the available microarray data from a variety of rice samples with inoculation of Xoo and Mor. A set of 12,932 genes was identified to be regulated by Xoo and another set of 2709 Mor-regulated genes was determined. GO enrichment analysis of the regulated genes by Xoo or Mor suggested mitochondrion may be an arena for the up-regulated genes and chloroplast be another for the down-regulated genes by Xoo or Mor. Cytokinin-related processes were most frequently repressed by Xoo, while processes relevant to jasmonic acid and abscisic acid were most frequently activated by Xoo and Mor. Among genes responsive to Xoo and Mor, defense responses and diverse signaling pathways were the most frequently enriched resistance mechanisms. InterPro annotation showed the zinc finger domain family, WRKY proteins, and Myb domain proteins were the most significant transcription factors regulated by Xoo and Mor. KEGG analysis demonstrated pathways including ‘phenylpropanoid biosynthesis’, ‘biosynthesis of antibiotics’, ‘phenylalanine metabolism’, and ‘biosynthesis of secondary metabolites’ were most frequently triggered by Xoo and Mor, whereas ‘circadian rhythm-plant’ was the most frequent pathway repressed by Xoo and Mor. The genes identified here represent a robust set of genes responsive to the infections of Xoo and Mor, which provides an overview of transcriptional reprogramming during rice defense against Xoo and Mor infections. Our study would be helpful in understanding the mechanisms of rice disease resistance.

hypersensitive reaction (HR), a specific programmed cell death (PCD), is often observed in plants. Many components have been demonstrated essential in the PTI and ETI processes in Arabidopsis thaliana [2,6,7].
Bacterial leaf blight is the most significant bacterial disease of rice. Its causal agent Xanthomonas oryzae pv. oryzae (Xoo), is a member of the gamma subdivision of the proteobacteria. Another popular rice disease is known as rice blast, which is caused by a filamentous ascomycete fungus, named Magnaporthe oryzae (Mor). Although disease resistance in rice has been extensively studied, there is still a little knowledge of the rice response to pathogen infections. It has been demonstrated that the invasions of Xoo and Mor on rice plants are mainly mediated through altering rice gene expression at the transcriptional level [8][9][10][11]. Hence, uncovering the transcriptional changes of rice genes during the infections of Xoo and Mor is of particular significance.
Up to date, plenty of microarray data of rice infected by Xoo and Mor have been produced, and these data give an opportunity to elucidate the mechanisms of rice response to the infections of the two pathogens. Previous microarray-based studies, however, focused on only a limited of samples. Although numerous rice genes responsive to Xoo and Mor infections have been identified, which led to insights into the rice resistance/susceptibility mechanisms, similar or specific results were usually obtained in different studies for the sake of specific experimental conditions.
The aim of this study is to determine a robust set of rice genes in response to Xoo and Mor infections, defining genes that are frequently regulated in diverse conditions through analysis of the publicly available rice microarray data sets associated with the infections of Xoo and Mor. To find out which mechanisms may be more common in rice response to the infections by these two pathogens, the gene sets determined were next analyzed for enrichment of Gene Ontology (GO) terms and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. Furthermore, the distributions of the significant enrichment of the GO and KEGG as well as the InterPro annotation were investigated.

Results
Identification of genes responsive to Xoo and Mor infections in the rice microarray data We employed 69 pairs (including control and treatment) of microarray samples (consisting of 51 pairs of Xoo-infected samples and 18 pairs of Mor-infected samples) from 15 series of experiments to query Xooand Mor-induced gene expression changes, all of which were conducted by the use of the Affymetrix rice wholegenome arrays platform (GPL2025) (Additional file 4: Table S1). The data from the same platform were retrieved for the analysis to avoid the variance of different platforms, and the GEO2R tool was used to process all the samples uniformly to eliminate the technical variance of data transformation. Moreover, poor quality arrays with no match or matching multiple loci were discarded. Further, we only considered the samples with no less than 989 differentially expressed genes (DEGs) (P ≤ 0.05). In the end, we identified the DEGs from 29 pairs of Xoo-infected samples and 6 pairs of Mor-infected samples (Additional file 5: Table S2). The number of DEGs in these samples varied from 989 to 9769 genes (Additional file 5: Table S2). Totally, 12,932 DEGs (occurring at least three pairs of array samples) were identified in the Xoo-infected rice microarray data ( Table 1). Of the DEGs, 7452 genes were up-regulated and 5480 genes were down-regulated (Table 1). Also, 2709 DEGs (occurring at least three pairs of array samples) were identified in the Mor-infected rice microarray data ( Table 1). Out of these genes, 1615 were upregulated and 1094 were down-regulated by Mor (Table 1).
Comparing the two groups of the identified DEGs, we found that 11,075 genes were expressed differentially in common between Xooand Morinfected rice array samples, with 5580 DEGs being up-regulated and 5495 DEGs being down-regulated (Fig. 1). If only considering the DEGs that were present in at least three pairs of samples infected by Xoo or in at least two pairs of samples infected by Mor, it could be found that 3831 DEGs were shared by Xooand Morinfected array samples, with 2140 DEGs being up-regulated and 1691 DEGs being down-regulated (Additional file 1: Figure S1).
In a previous study, we showed that 882 rice genes contain pathogen-inducible cis-regulatory elements (PICEs) in their promoter regions [24]. Here we found that 389 DEGs contain the PICEs in their promoters (Additional file 3: Table S3). There is a 3.51% (389/11075) overlap or a 44.1% (389/882) overlap between the DEGs and the genes with PICEs in the promoters. And among the above 3831 DEGs, 304 genes were found to be overlapped with the prior set of genes carrying the PICEs in their promoters, which account for 7.93% of the DEGs (Additional file 7: Table S4). It seems that the PICEs in the promoters make genes accessible to be frequently regulated by pathogens.
GO enrichment analysis of the DEGs in the Xooand Morinfected rice samples GO enrichment analysis provides some detailed information on the potential functions of genes. However, the existing GO terms are too disorganized to succinctly describe the functional information of a large set of genes. Here, combining the routine GO terms, we tentatively use more broad terms to describe the GO enrichment results of the DEGs in the Xooand Mor-infected rice samples. GO enrichment analysis revealed that ribosome, snRNP complex and spliceosome, and diverse membranes were the most frequent cellular components in the up-regulated genes by Xoo (Fig. 2a and Additional file 8: Table S5). In contrast, chloroplast, membrane, and ribosome were the most frequent cellular components in the down-regulated genes by Xoo ( Fig. 2a and Additional file 9: Table S6).
Analysis on the molecular functions of the genes induced by Xoo showed that RNA binding, defense enzymes and other enzymes were significantly enriched in the up-regulated genes ( Fig. 2a and Additional file 12: Table S9). By contrast, defense enzymes and other enzymes, as well as transporters were the frequently enriched molecular function categories in the downregulated genes by Xoo ( Fig. 2a and Additional file 13: Table S10).
In the up-regulated genes by Mor, the analysis presented that the most frequently enriched cellular components were membrane, mitochondrion, proteasome, and ribosome ( Fig. 2b and Additional file 14: Table S11); whereas the most frequently enriched cellular components were chloroplast, snRNP complex and   Table S12).
On the biological processes that the induced genes by Mor involved in, enrichment analysis indicated that defense response, carbohydrate biosynthesis, transport and metabolism, and hormone regulation were the most frequent processes present in the up-regulated genes ( Fig. 2b and Additional file 16: Table S13); while transcription, post-transcriptional regulation, chloroplast and photosynthesis, and translation, post-translational regulation were the most frequent processes occurring in the down-regulated genes ( Fig. 2b and Additional file 17: Table S14).
Enrichment on the GO molecular functions of the regulated genes by Mor showed that defense enzymes, other enzymes, metal ion binding, and transporter were the most frequent categories in the up-regulated genes ( Fig. 2b and Additional file 18: Table S15); while in the down-regulated genes, nucleic acid binding, defense enzymes and other enzymes were included in the most frequent molecular functions ( Fig. 2b and Additional file 19: Table S16).
It has been established that most hormones are important in regulating rice disease resistance. Thus, we investigated the biological processes related to diverse hormones, which were activated or repressed after the Fig. 2 The distribution of the gene ontology annotation (GO) for the induced genes in the rice microarray data. a GO for the induced genes by Xoo; b GO for the induced genes by Mor infections of Xoo and Mor. We observed that some biological processes relative to hormones included the regulation of hormone-mediated signaling pathways, hormone biosynthetic/metabolic processes, and response to hormone (Fig. 3). Among the hormones, jasmonic acid and abscisic acid were prominent for the processes related to them were most frequently activated after infections by Xoo and Mor (Fig. 3). It is worth noting that the processes related to cytokinin (GO:0009735 and GO: 0009736) were most frequently repressed after Xoo infection (Fig. 3a). In samples infected by Mor, the processes relative to diverse hormones consisting of abscisic acid, auxin, cytokinin, jasmonic acid, and ethylene were evenly repressed (Fig. 3b).
Comparison of the most frequently observed GO enrichment results of the induced genes by Xoo and Mor revealed some interesting phenomena. Ribosome, membrane, mitochondria, and chloroplast were frequently observed, suggesting they were important cellular components where many genes were induced during the infections of Xoo or Mor. However, the occurring frequencies were different between the up-and downregulated patterns (Fig. 5). For example, the mitochondrion was a cellular component where the occurring frequency was higher in the up-regulated genes than the down-regulated genes by Xoo and Mor ( Fig. 2 and Fig.  5). Conversely, the occurring frequency of the chloroplast was higher in the down-regulated genes than the up-regulated genes by both of the pathogens ( Fig. 2 and Fig. 5). Therefore, when rice plants are subject to Xoo and Mor infections, mitochondria and chloroplast will turn into two important arenas, with the former being the one for up-regulated genes and the latter being another for down-regulated genes. Also, among the biological processes, defense response and hormone regulation were observed to be commonly and frequently activated by these two pathogens through up-regulating some related genes; and chloroplast organization and photosynthesis, development, and photomorphogenesis were commonly and frequently repressed by the two pathogens through down-regulating the relevant genes ( Fig. 2 and Fig. 5). Similarly, the molecular functions including defense enzymes, ribosome structural constituent and binding, and polysaccharide binding were commonly and frequently observed in the up-regulated genes by the two pathogens ( Fig. 2 and Fig. 5).
In addition, some frequent GO enrichment results were observed to be specific (Fig. 5). For example, Fig. 3 The distribution of the biological process (BP) relative to hormone regulation for the induced genes in the rice microarray data. a BP for the induced genes by Xoo; b BP for the induced genes by Mor 'snRNP complex and spliceosome' was only frequently occurring in the genes up-regulated by Xoo and downregulated by Mor; snoRNP was only frequently occurring in the up-regulated genes by Xoo ( Fig. 2 and Fig. 5). The genes involved in ribosome biogenesis and assembly process were only frequently observed among the genes with an up-regulated pattern of expression induced by Xoo and Mor ( Fig. 2 and Fig. 5). And the signal transduction process was only observed among the genes down-regulated by Mor ( Fig. 2b and Fig. 5). The genes with polysaccharide binding function were exclusively identified in those with an up-regulated pattern of expression induced by the two pathogens ( Fig. 2 and Fig. 5). Hence, some results obtained from the analysis of microarray data infected by Mor further supported that from the analysis of data infected by Xoo; on the other hand, comparison of the data infected by Mor to that by Xoo indicated that rice plants can use different mechanisms in response to their infections.

Enriched InterPro annotation of the DEGs in rice induced by Xoo and Mor infections
To further explore the possible functions of the DEGs in rice induced by Xoo and Mor infections, InterPro annotation enrichment analysis was conducted. The most frequent and significant InterPro annotations associated with the up-regulated DEGs by Xoo were diverse kinds of enzymes, e.g., dehydrogenase, synthase/synthetase, transferase, kinase, and glycoside hydrolase, followed by ribosomal proteins, translation proteins, and transcription factors; while among the down-regulated DEGs by Xoo, the most frequent and significant InterPro annotations were diverse enzymes including transferase, synthase/synthetase, reductase, dehydrogenase, and peptidase, etc., followed by transporter, transcription factors and ribosomal proteins ( Table 2, Additional file 20: Table S17 and Additional file 21: S18).
In the up-regulated DEGs by Mor, the analysis showed that diverse kinds of enzymes (e.g., dehydrogenase, kinase, transferase, hydrolase, and synthase/ synthetase), transporters, and transcription factors were the top three enriched InterPro annotations; and in the down-regulated DEGs by Mor, it was demonstrated that various enzymes including kinase, transferase, synthase/synthetase, and so on, were annotated to be the most frequent products associated with the DEGs, followed by transcription factors, heat shock proteins and WD40 domain proteins (Table 3,  Additional file 22: Table S19 and Additional file 23: Table S20.). Therefore, a number of enzymes were induced in rice response to Xoo and Mor infections, with some being up-regulated and others down-regulated. Among the enzymes, kinase, transferase, and synthase/synthetase were frequently activated or repressed during the rice response to Xoo and Mor infections. According to the above GO molecular function analysis of the DEGs, some activated enzymes belong to defense enzymes.

Disease resistance/susceptibility-related (DRR/DSR) genes associated with Xoo and Mor infections in rice
To identify the DRR/DSR genes, we performed a more detailed InterPro analysis on the Xooand Mor-regulated genes, combining with literature mining. Among the upregulated genes by Xoo, the analysis indicated that genes from 106 InterPro annotated items were identified as DRR/DSR genes and the most frequently up-regulated genes by Xoo encode NAD(P)-binding domain proteins (IPR016040), which involve in reactive oxygen species (ROS) and SA signaling (Additional file 24: Table S21). Genes encoding chaperonin Cpn60/TCP-1 (IPR002423) and GroEL-like apical domain proteins (IPR027409) were also frequently activated, which are related to PCD and defense response, respectively (Additional file 24: Table S21). Among the down-regulated genes by Xoo, the analysis showed that genes from 75 InterPro items were identified as DRR/DSR genes (Additional file 25: Table S22). Three groups of genes were most frequently down-regulated by Xoo, and they encode NAD(P)-binding domain proteins (IPR016040) (functioning in ROS and SA signaling), major facilitator superfamily (IPR020846 and IPR011701) (acting as defense proteins), and NAF/FISL domain proteins (IPR004041 and IPR018451) (involved in PTI), respectively (Additional file 25: Table S22). Among genes responsive to Mor, the results demonstrated that genes from 39 InterPro annotated items were identified to be DRR/DSR genes, with up-regulated expression pattern (Additional file 26: Table S23). Interestingly, the most frequently activated genes by Mor also encode NAD(P)-binding domain proteins (IPR016040), and function in ROS and SA signaling (Additional file 26: Table S23). Among the downregulated genes by Mor, the putative DRR/DSR genes were confirmed to be distributed in 19 InterPro annotated items (Additional file 27: Table S24). Three groups of genes, which encode PB1 domain (Phox/Bem1p) proteins (IPR000270) (involved in defense response), DnaJ domain proteins (IPR001623) (implicated in cell death), and CCT domain proteins (IPR010402) (related to defense response), respectively, were most frequently repressed by Mor infection (Additional file 27: Table S24).
We further investigated the mechanisms on disease resistance of the identified DRR/DSR genes in rice response to Xoo and Mor infections through consulting a great number of papers. The results showed that various disease resistance mechanisms were conferred by the up-regulated genes by Xoo (Fig. 4a). For example, 29% (264) of the up-regulated genes by Xoo played roles during defense responses or as defense proteins. One hundred ninety six genes (22%) were probably involved in diverse signaling pathways, including SA, JA, ET, MAPK, receptor kinase and so on, among the up-regulated genes by Xoo. 12% (105) of the up-regulated genes by Xoo were associated with PCD, HR or other cell death. A group of 101 up-regulated genes by Xoo (11%) was related to basal and innate immunity including PTI, ETI. 7% (61) were implicated in ROS/oxidative stress, and 3% (30) were found as transcription factors, among the upregulated genes by Xoo (Fig. 4a). A similar result was observed among the down-regulated genes by Xoo (Fig.  4b). For instance, 22% (161) of the down-regulated genes by Xoo, the greatest group of genes, played a part in defense responses or as defense proteins. A group of 146 genes (20%) was related to diverse signaling pathways mediated by SA, JA, ET, ABA, MAPK, and receptor kinase, among the down-regulated genes by Xoo. 10% (76) of the down-regulated genes by Xoo were in association with PCD, HR and other cell death. 9% (65) were involved in ROS, oxidative stress or antioxidant related protein. Another group of 51 down-regulated genes by Xoo (7%) was implicated in basal and innate immunity including PTI, ETI. In addition, 48 down-regulated genes by Xoo (6%) were present, which encoded transcription factors (Fig. 4b).
Analysis of the genes regulated by Mor revealed some disease resistance mechanisms (Additional file 1: Figure  S1). Defense responses (49 genes, 27%) were observed to be most frequently induced in rice after infection by Mor. A group of 32 genes (18%) activated by Mor was implicated in signaling mediated by SA, ET, MAPK, etc. Some induced genes encoded antifungal proteins or were related to the synthesis of phytoalexin. Induced genes involving in PCD, HR and other cell death, and genes related to ROS accounted for 10% (19) and 8% (15), respectively (Additional file 2: Figure S2). Among the down-regulated genes by Mor, 54% (27) encoded transcription factors; ten genes (20%) were associated with cell death; four genes (8%) were present in defense responses; additionally, several genes were related to auxin signaling, ethylene signaling, and NPR1 degradation, respectively (Additional file 3: Figure S3). Top regulated genes after Xoo and Mor infections We next determined the top regulated genes by Xoo and Mor in rice. We found that nine genes, which were occurring in at least 18 pairs of samples, were up-regulated after Xoo infection (Table 4). Os04g0650800, which encodes the D-3-phosphoglycerate dehydrogenase 3 in chloroplast, was found to be most frequently upregulated, with occurring in 22 pairs of Xoo-infected samples ( Table 4). Four of the up-regulated genes are involved in defense responses (Os03g0235000, Os09g0 491772, Os08g0508800, and Os09g0484200) ( Table 4).
Another nine genes were observed to be frequently down-regulated by Xoo, with happening in at least 18 pairs of samples (Table 5). Among the genes, one encodes the cinnamoyl-CoA reductase 1 (Os09g0491820), and was most frequently down-regulated by Xoo, with occurring in 21 pairs of samples; two genes (Os05g0204600 and Os06g0571800) encode different transcription factors (Table 5).
In Mor-infected samples, 16 genes were observed to be up-regulated in all the six pairs of samples analyzed (Table 6). Among these genes, five belong to defense genes (Os01g0713200, Os01g0963000, Os02g0569900, Os02g0570700, and Os10g0542900), one encodes a WRKY transcription factor (Os01g0584900), and another gene (Os04g0647900) encodes an LRR receptor-like serine/threonine-protein kinase ( Table 6). In the samples, ten genes were found to be down-regulated, with occurring in five pairs of samples ( Table 7). Out of the genes, Os01g0719100 and Os03g0607700 encode two different zinc finger proteins, and Os09g0536400 encodes a defense protein (thaumatin-like protein 1b) ( Table 7).
Of the up-regulated transcription factors in the Morinfected samples, the WRKY proteins (IPR003657) were the maximum of the induced transcription factors, followed by the Myb domain factors (IPR006447, IPR017930) (Additional file 29: Table S26). Among the down-regulated transcription factors in these samples, various kinds of zinc finger domain proteins were also the most frequently induced transcription factors, and the Myb domain proteins (IPR006447, IPR017930, and IPR001005) were the second frequently induced (Additional file 29: Table S26). Therefore, the zinc finger domain family, WRKY proteins, and Myb domain proteins were the most significant proteins among the differentially expressed transcription factors in most of the Xooand Mor-infected samples (Fig. 5).

KEGG enrichment analysis of the DEGs induced by Xoo and Mor in rice
To search the enriched pathways potentially targeted by the DEGs in rice induced by Xoo and Mor infections, KEGG annotation was performed. The analysis indicated Similarly the terms in light blue down-pentagon and those in both the sides of the down-pentagon demonstrate the terms that were down-regulated in the infections by both pathogens and specific to one infection respectively. Terms in light yellow rectangles were found both in up-and down-regulated gene sets in response to one infection respectively. Terms in brown rectangle were found both in up-and down-regulated gene sets that genes involved in the following pathways were most frequently enriched in the up-regulated genes by Xoo: 'phenylpropanoid biosynthesis' (osa00940), 'biosynthesis of amino acids' (osa01230), 'biosynthesis of antibiotics' (osa01130), 'ribosome biogenesis in eukaryotes' (osa03008), and 'biosynthesis of secondary metabolites' (osa01110) pathways (Additional file 30: Table S27). In the down-regulated genes by Xoo, the most frequent pathways targeted were 'biosynthesis of secondary metabolites' (osa01110), 'metabolic pathways' (osa01100), 'biosynthesis of antibiotics' (osa01130), and 'carbon metabolism' (osa01200) (Additional file 31: Table S28).
To better outline the enrichment results of the KEGG targeted by the DEGs in rice induced by Xoo and Mor, all the pathways enriched were divided into eight categories ( Table 8, Table 10). As shown in Table 8, disease resistance related pathways, carbohydrate metabolism, and genetic information processing were the top three frequently enriched categories in the up-regulated genes by Xoo; and in the down-regulated genes by Xoo, carbohydrate metabolism, disease resistance related pathways, and biosynthesis and metabolism of amino acids were observed to be the most frequently enriched pathways. Comparison of the pathways targeted by the genes induced by Xoo indicates that the pathways related to genetic information processing and lipid metabolism were more frequently enriched in the up-regulated genes, while the pathways involving metabolism of cofactors and vitamins were more frequently enriched in the down-regulated genes (Table 8).
We next analyzed the KEGG enrichment of the genes induced by Mor. This analysis led to a similar KEGG enrichment of the up-regulated genes by Mor to that of the genes activated by Xoo. As listed in Additional file 32: Table S29, 'biosynthesis of secondary metabolites' (osa01110), 'biosynthesis of amino acids' (osa01230), 'phenylpropanoid biosynthesis' (osa00940), and 'biosynthesis of antibiotics' (osa01130) were also frequently enriched in Mor-up-regulated genes. However, in the Mor-down-regulated genes, the frequent KEGG enrichment annotations were 'plant hormone signal transduction' (osa04075), 'circadian rhythm-plant' (osa04712), and 'porphyrin and chlorophyll metabolism' (osa00860)  (Additional file 33: Table S30). Broad KEGG enrichment results indicated that disease resistance related pathways, biosynthesis and metabolism of amino acids, and carbohydrate metabolism categories were most frequently enriched in the up-regulated genes by Mor; while in the genes repressed by Mor, the most frequent enrichment categories of KEGG were disease resistance related pathways, genetic information processing, and biosynthesis and metabolism of amino acids (Table 10). Further analysis showed that the most frequently activated disease resistance related pathway by Mor in rice was 'phenylalanine, tyrosine and tryptophan biosynthesis' (osa00400), followed by 'phenylalanine metabolism' (osa00360), 'phenylpropanoid biosynthesis' (osa00940), 'biosynthesis of antibiotics' (osa01130), and 'plant-pathogen interaction' (osa04626) ( Table 11). In contrast, 'plant hormone signal transduction' (osa04075), and 'circadian rhythm-plant' (osa04712) were the most frequently repressed pathways related to disease resistance by Mor in rice (Table 11).

Discussion
Transcriptomic analysis is a powerful tool to reveal the interactions between host and pathogen [61][62][63][64]. In this study, we analyzed 35 pairs of samples from 69 pairs of Xooor Mor-infected rice microarray samples deposited in the GEO database. Some rice genes were identified to be frequently regulated by Xoo and Mor in various conditions. During the interactions between rice and these two pathogens, the pathogen attacks coincided with rice defenses. We separated the up-regulated genes and down-regulated genes by pathogens in rice, in other words, there was no gene overlap between them. GO, InterPro and KEGG annotations and enrichment analyses were performed respectively on the up-and down-regulated genes in each rice sample infected by Xoo or Mor. We further revealed the frequent GO, Inter-Pro and KEGG annotations and enrichment results of rice genes regulated by Xoo and Mor.
Some studies have been carried out to investigate the responsive mechanisms in rice to the infections of Xoo and Mor in specific conditions. The roles of resistance genes were focused on in early studies on rice response  to Xoo infection [12]. Recently, analyses of transcriptomic data revealed the detailed responses to Xoo and Mor infections in rice. It has been demonstrated that the transcriptional changes during the compatible interaction between rice and Mor were very similar to the changes of incompatible interaction; nevertheless, the changes were higher during the incompatible interaction [8]. In plants including rice, transcriptional reprogramming is a universal process during the infections of pathogens [9,65,66]. Our results indicate that the frequently regulated genes involved in basic biological processes from transcription to translation, and their regulation. A number of researches reported the expression profiles of rice defense response genes [9,67]. For example, chitinase activities were detected to play roles in rice resistance to Mor [68,69]. Our results show that the expected defense enzymes were the most frequent products of the regulated genes by Xoo and Mor (Fig. 2, Fig.  5, Additional file 12: Table S9, Additional file 11: S8. Additional file 16: S13. and Additional file 17: S14).
Genes related to secondary metabolism were observed to be up-regulated and significantly enriched in rice infected by Mor [8]. In Arabidopsis, metabolic pathways were observed to be frequently influenced by diverse pathogens [70]. The resistance to Mor was shown in wild rice plants to be involved in lipid metabolism, phenylpropanoid and diterpenoid metabolism [69]. During the interaction of a pair of rice near-isogenic lines (NILs) with Xoo, phenylpropanoid biosynthesis was identified as the most conspicuous pathway [11]. In our analysis, KEGG enrichment indicated that 'phenylpropanoid biosynthesis' (osa00940) and 'biosynthesis of secondary metabolites' (osa01110) were observed to be two of the most frequently enriched pathways in the up-regulated genes by Xoo and Mor (Fig. 5, Additional file 30: Table  S27 and Additional file 32: Table S29). However, 'biosynthesis of secondary metabolites' (osa01110) was also the most frequent pathway targeted by the down-regulated genes by Xoo (Additional file 31: Table S28). Our analysis further showed 'biosynthesis of antibiotics' (osa01130) and 'biosynthesis of amino acids' (osa01230) were another two of the pathways most frequently enriched by the up-regulated genes by Xoo and Mor (Fig. 5, Additional file 30: Table S27 and Additional file 32: Table S29).
Analysis of the DRR/DSR genes responsive to Xoo and Mor infections suggests that various disease resistance mechanisms were activated or repressed in rice. Among the mechanisms, defense response, signaling pathways, PCD, HR and other cell death, basal and innate immunity were the most eminent arenas that fight between rice and these two pathogens (Fig. 4, Additional file 1: Figure  S1 and Additional file 2: Figure S2). It is worth noting that genes encoding antifungal proteins and related to the synthesis of phytoalexin were activated in rice after Mor but Xoo infection, suggesting that different mechanisms were triggered ( Fig. 5 and Additional file 2: Figure S2).
WRKY transcription factors were again and again observed to be implicated in the rice resistance to Mor [8,68,69]. In our analysis, WRKY proteins (IPR003657) were the most frequently activated transcription factors by Mor (Table S26). The WRKY factors were also frequently up-regulated by Xoo (Fig. 5, Additional file 8: Table S5). In addition, the zinc finger domain proteins, and the Myb domain proteins were demonstrated to be the most significant factors with differential expressions in most samples infected by Xoo and Mor (Fig. 5, Additional file 8: Table S5 and Additional file 27: Table S24).
Signalings mediated by abscisic acid and cytokinin were revealed to be successively activated by Mor infection in rice [22]. Abscisic acid was also suggested being important in the rice (carrying Xa7) response to Xoo infection and high temperature stress [10]. However, our results indicate that the processes relative to abscisic acid may not only be frequently activated by the infections of Xoo and Mor (GO:0009737~response to abscisic acid; GO:0009738~abscisic acid-activated signaling pathway) (Fig. 3 and Fig. 5), but also be repressed by Mor infection (GO:0009738~abscisic acid-activated signaling pathway; GO:0009788~negative regulation of abscisic acid-activated signaling pathway) (Fig. 3b). In our results, processes with relation to cytokinin were strikingly and frequently repressed by Xoo and Mor infections (GO:0009735~response to cytokinin; GO:0009736~cytokinin-activated signaling pathway) ( Fig. 3 and Fig. 5). Two pathways mediated by jasmonic acid and ethylene were demonstrated to be activated by Xoo and Mor infections [11,68,69]. Our results show that the processes relative to jasmonic acid were frequently activated in the Xooand Mor-infected rice (GO: 2000022~regulation of jasmonic acid mediated signaling pathway; GO:0009695~jasmonic acid biosynthetic process; GO:0009753~response to jasmonic acid) ( Fig. 3 and Fig. 5), and they may be repressed by Mor infection (Fig. 3b and Fig. 5). The processes relative to ethylene (GO:0009693~ethylene biosynthetic process; GO: 0009723~response to ethylene; GO:0009873~ethyleneactivated signaling pathway) were also repressed in Morinfected rice (Fig. 3b).
Previous studies have shown that the mechanisms by which rice plants respond to the infections of Xoo and Mor are not quite identical or even completely opposite under different experimental systems and experimental conditions. Our results reveal which mechanisms are occurring with a relatively high probability.
We show that Os04g0650800 was most frequently activated by Xoo. This gene encodes the D-3-phosphoglycerate dehydrogenase 3 in chloroplast. It has been established that the phosphorylated pathway in plastidial glycolysis is one of three different serine biosynthesis pathways in plants, and 3-phosphoglycerate dehydrogenase (PGDH) is the first committed enzyme [71]. PGDH was revealed to be quantitatively the most important enzyme in maintaining serine homeostasis at the whole plant level [72]. Our result suggests that Os04g0650800 is likely a key gene in rice defense response and the activation of serine biosynthesis dependent on a phosphorylated pathway may be important in rice response to Xoo invasion.
As shown above, the genes responsive to Xoo and Mor infections may be up-regulated or down-regulated. Further experiments of molecular genetics are necessary to determine that the up-or down-regulated expressions of these genes are active or passive, and to reveal the significance of the genes with highly frequent expressions in response to the infections of these two pathogens. In this study, only six pairs of Mor-infected samples were retrieved and analyzed, therefore more extensive analyses are needed on the genes and mechanisms responsive to Mor infection in rice.

Conclusions
A robust set of genes has been defined in rice response to the infections of Xoo and Mor. Mitochondrion and chloroplast may be important organelles for rice response to Xoo and Mor infections. Processes with relation to cytokinin, jasmonic acid, and abscisic acid were most frequently operated by Xoo and Mor. A great number of enzymes were in rice favored to be manipulated by Xoo and Mor. Defense responses and diverse signaling pathways were indispensable among the responsive mechanisms of rice to Xoo and Mor. Pathways including phenylpropanoid biosynthesis, biosynthesis of antibiotics, phenylalanine metabolism, and biosynthesis of secondary metabolites were most frequently triggered by Xoo and Mor. Circadian rhythm-plant was the most frequent pathway repressed by Xoo and Mor.

Materials and methods
Retrieval and analysis of the microarray datasets of rice infected by Xoo and Mor The Xooand Mor-infected microarray data were collected from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) [73,74]. Sixty-nine pairs (including control and treatment) of raw microarray data (.CEL) of rice in response to the infections of Xoo and Mor were obtained by SRA Toolkit. The microarray data included 15 series data sets consisting of 51 pairs of Xoo-infected rice samples and 18 pairs of Mor-infected rice samples. All the microarray experiments were conducted on the same Affymetrix platform (Rice Genome Array GPL2025). The Bioconductor package Simpleaffy was used for quality assessment [75]. GEOquery and limma R packages from GEO2R were used to identify the differential expression of all samples [76]. Poor quality samples were discarded and all the rest of samples (35 pairs of samples) were used to analyze the differentially expressed genes.

Identification of the differentially expressed genes (DEGs) induced by Xoo and Mor infections
A series of data generated by GEO2R were then evaluated and analyzed. The obtained p-values for multiple testing were corrected according to Benjamini and Hochberg procedure [77]. Probe IDs were converted to gene symbols according to the GPL2025 annotation file provided by GEO. To identify the DEGs, an adjusted pvalue was used to filter genes with no differential expressions. By setting the threshold of adj. P-value < 0.05, the genes with differential expression induced by Xoo and Mor were screened. Given that more than 30,000 genes were found in rice genome, and more than 1000 DEGs were indentified in many rice samples infected by pathogens, the infected samples with too few DEGs may be resulted from unknown experiment factors. Among the analyzed samples, 34 samples were found with more than 1000 DEGs, and a sample was found with nearly 1000 DEGs (989 DEGs). To avoid inaccurate results, samples with less than 989 DEGs were discarded. In the discarded 34 samples, no DEGs could be identified in 13 samples, DEGs with less than 200 genes were identified from nine samples, and the other samples were found with several hundred DEGs.
To reveal the potential function tendency of the above identified DEGs, GO enrichment was conducted through the DAVID tool (https://david.ncifcrf. gov/) [81,82]. And the enrichment analysis was deduced based on p-value ≤0.05. Pathway analysis was also carried out through the DAVID tool to unravel the functional annotation of the identified DEGs, and the pathways with p-value ≤0.05 were retrieved as significant enrichment [81,82].