Identification of rice OsCBP60s genes
The protein sequence of all the AtCBP60s genes (Reddy et al., 2002) were extracted using Phytozome11 (https://phytozome.jgi.doe.gov/pz/portal.html). BLASTP searches were performed using the complete protein sequence of each AtCBP60 protein separately as the query sequence and total list of OsCBP60s were prepared. The duplicate entries and apparent incomplete sequences were removed manually. Fifteen candidate OsCBP60s genes were obtained (Table S1). Further, protein sequences of 15 OsCBP60s and 8 AtCBP60s were used to create a phylogenetic tree using neighbor joining (N-J) tree with bootstrapping (500 reps) in MEGA7 software (http://www.megasoftware.net/) (Fig. 1). Two major clusters were obtained, which were further divided into different sub cluster as shown in Fig. 1. Majority of the AtCBP60s were grouped in cluster I except SARD1 and AtCBP60g. Both SARD1 and AtCBP60g were previously shown to be involved in systemic acquired resistance (SAR) in Arabidopsis (Zhang et al. 2010). However, OsCBP60s were dispersed in both the cluster. Interestingly, LOC_Os01g04280 (OsCBP60_14) and SARD1 (AtCB60h) were fell within the same clade.
Protein domain search
The Pfam analysis showed that like AtCBP60s all the identified OsCBP60s contains only one domain i.e. calmodulin binding protein (CBP) domain (Fig. 2). However, the size and position of CBP domain varies in Arabidopsis and rice CBP60s. The size of CBP domain varies from 283 to 346 aa.
Subcellular localization
In order to predict the subcellular localization, the protein sequences of AtCBP60s and OsCBP60s were subjected to WoLF PSORT (https://www.genscript.com/wolf-psort.html) analysis. The subcellular localization studies showed that AtCBP60s and OsCBP60s localizes to various cell organelles like nucleus, mitochondria, cytoskeleton, vacuoles, peroxisome and cytoplasm (Supplementary Table S1). Interestingly, eight OsCBP60s (OsCBP60_1, OsCBP60_3, OsCBP60_4, OsCBP60_6, OsCBP60_8, OsCBP60_9, OsCBP60_11, OsCBP60_12) and five AtCBP60s (At5g57580, At2g18750, At4g25800, At2g24300, At4g31000) were maximally predicted to localize to nucleus. Six OsCBP60s (OsCBP60_13, OsCBP60_14, OsCBP60_15, OsCBP60_5, OsCBP60_7, OsCBP60_10 were predicted to localize to chloroplast (Table S1). Majority of OsCBP60s gene were found on Chr 12 followed by Chr 3.
Identification of transcription factor binding sits (TFBSs)
The upstream sequences (1000 bp) of OsCBP60s genes were analysed using PlantPan VERSION. 2/. Several TFBSs involved in stress responses were found to be enriched in the promoter sequence of OsCBP60s (Fig. S1). In total 1529 stress responsive TFBSs were observed in positive or negative strand of OsCBP60s (Fig. 3 and 4).
A number of TFBSs associated with light stress include ATCT-motif, TCCC-motif, TCT-motif, GT1-motif, G-box, GATA-motif, AE-box, GTGGC-motif, Sp1, Box II, ACE, MRE are present in OsCBP60s (Supplementary Table S3). Apart from this other stress motifs MYC, STRE, TC-rich repeats, CCAAT-box, LTR, MBS, GC-motif, ARE, WRE3, W box, CARE, MYB-like sequence, MYB recognition site, MYB, Myb, AC-I, WUN-motif, CCAAT-box, as-1 (Supplementary Table S3). Oxidative stress responsive elements like ARE (AAACCA), ethylene-responsive ERE (ATTTTAAA), AREs (Anaerobic responsive elements), GC-motif (CCCCCG) presents in some of members of OsCBP60s genes with consequences TGACG (Table S3).
TFBSs involved in hormonal regulation like methyl-jasmonate CGTCA and TGACG-motif , AuxRR-core, TGA-element (auxin responsiveness), AAGAA-motif, DRE core, DRE1 and ABRE (abscisic acid), TCA element (salicylic acid), TATC-box, GARE-motif and P-box (gibberellin-responsiveness) were found to be enriched in OsCBP60s (Supplementary Table S3). These motifs are present in studied OsCBP60s genes (Table S3).
Pathogenesis genes play qualitative role in cellular development and this can be clearly seen by presence of motifs involved in cellular development viz. CAT-box (meristem expression), RY-element (seed-specific regulation), circadian, GCN4_motif (endosperm expression), O2-site (zein metabolism regulation), MSA-like (cell cycle regulation), CCGTCC-box, NON (CAACGGCCACG), MBSI (flavonoid biosynthetic) and dOCT (CACGGATC) are important for tissue specific expression (Table S3).
Expression profiling of OsCBP60s after pathogens treatment
To gain further information on the biological function of OsCBP60s, Quantitative real-time reverse-transcriptase (RT) polymerase chain reaction (qPCR) analysis were performed in rice seedlings treated with two most devastating rice pathogens (M. oryzae and X. oryzae) at 12, 24 and 48 h after treatment. In present investigation, differential expression of OsCBP60s was observed in rice seedlings treated with pathogens (Fig. 5 and 6). OsPR1a was also found to be upregulated in all the time points studied after treatment with M. oryzae and X. oryzae (Fig. 5 and 6). The upregulation of OsPR1a in all the time point studied indicates onset of defense mechanism in rice against both the pathogens
Treatment with M. oryzae leads to the up regulation of all OsCBP60s with different magnitude in all the time point studied. Five OsCBP60s (OsCBP60_5, OsCBP60_6, OsCBP60_7, OsCBP60_11 and OsCBP60_14) showed high level of upregulation at early time point i.e. 12 h. Interestingly, six OsCBP60s namely OsCBP60_4, OsCBP60_5, OsCBP60_9, OsCBP60_10, OsCBP60_12 and OsCBP60_15 showed consistently up regulated at 12 h, 24 h and 48 h (Fig. 5). Three OsCBP60s (OsCBP60_5, OsCBP60_6 and OsCBP60_11) showed high level of upregulation at 24 h after treatment with M. oryzae.
Treatment with X. oryzae caused differential expression of OsCBP60s (Fig. 6). qRT-PCR analysis showed that induced expression of 4 OsCBP60s (OsCBP60_6, OsCBP60_7, OsCBP60_8 and OsCBP60_9) occurred as early as 12 h after treatment with X. oryzae but downregulates at 24 h and 48 h. Interestingly four genes (OsCBP60_4, OsCBP60_5, OsCBP60_14 and OsCBP60_15) showed consistently up regulation at 12 h, 24 h and 48 h (Fig. 6) after treatment with X. oryzae. Taken together, OsCBP60_5, OsCBP60_10 and OsCBP60_15) were found to be commonly upregulated in M. oryzae and X. oryzae treatments.
Responses of OsCBP60s to defence related phytohormones treatment
Typically salicylic acid (SA) and jasmonic acid (JA) are considered to play a key role in plant responses to biotic stresses (Bari and Jones 2009; Prasad et al., 2009; Robert-Seilaniantz et al., 2011). The growing evidences also indicates the implication of brassinosteroid (BR) in defence responses in plant (Bari and Jones 2009; Sahni et al., 2016). To assess possible effects of these defence related phytohormones on the expression of OsCBP60s, qRT-PCR analysis were performed in rice seedlings treated with SA, JA and 24-epibrassinosteroid (EBR). The expression of OsPR1a was also analyzed to see the effectiveness of phytohormones treatment.
SA treatments showed the differential expression of OsCBP60s (Fig. 7). Induced expression of 4 OsCBP60s (OsCBP60_5, OsCBP60_13, OsCBP60_14 and OsCBP60_15) occurred as early as 6 h after treatment with SA (Fig. 7). Nine OsCBP60s (OsCBP60_1, OsCBP60_2, OsCBP60_5, OsCBP60_6, OsCBP60_7, OsCBP60_9, OsCBP60_10, OsCBP60_11 and OsCBP60_14) showed higher level of upregulation at 24h compared to other time point studied. Two OsCBP60s (OsCBP60_13 and OsCBP60_15) showed upregulation at 12 h, followed by down regulation at 24 h and 48 h of SA treatment (Fig. 7).
Most studies have identified the antagonistic interactions between the SA and JA mediated signaling pathways (Takahashi et al., 2004). The down regulation of OsPR1a at 24 h and 48 h after JA treatment confirms the effectiveness of treatment. Five OsCBP60s (OsCBP60_1, OsCBP60_5, OsCBP60_7, OsCBP60_13, and OsCBP60_14) showed downregulation in all the time studied after JA treatment (Fig. 8). However, the expression of two OsCBP60s (OsCBP60_4 and OsCBP60_15) were found to be upregulated at 12 h followed by downregulation at 24 h and 48 h of JA treatment (Fig. 8). Interestingly, 4 OsCBP60s (OsCBP60_7, OsCBP60_8, OsCBP60_9 and OsCBP60_12) were upregulated at all the time point studied after JA treatment (Fig. 8).
In EBR treated rice samples all the 15 OsCBP60s genes were found to be upregulated (Fig. 9). The maximum upregulation was observed in OsCBP60_6 followed by OsCBP60_7.
Taken together, OsCBP60_5 and OsCBP60_9 were found to be commonly upregulated in EBR and JA treatments. Similarly OsCBP60_1, OsCBP60_2, OsCBP60_3, OsCBP60_5, OsCBP60_13 and OsCBP60_14 were found to be upregulated in both EBR and SA treatments. OsCBP60_8 was found to be commonly upregulated in SA, JA, EBR phytohormones. These data helps to interpret that EBR mediates biotic stress responses in rice by modulating the both SA and JA pathways.