Meta-analysis of microarray and RNAseq data reveal OsbZIP52 to mediate salt stress responses in sensitive, tolerant and halophyte rice varieties
CABI Agriculture and Bioscience volume 4, Article number: 31 (2023)
The development of salt-tolerant rice has become urgent due to climate change and rising global rice consumption. A large-scale analysis using different but related platforms has become imperative to filter out candidate genes responsible for salinity tolerance and salinity stress-responsive pathways. Such genes can be used to find prospective candidate salt resistance genes in donor rice genotypes and transfer them to high-yielding rice varieties. We performed a meta-analysis to screen out candidate genes using stress-related three microarray and one RNASeq datasets from NCBI. As different genotypes of rice and different salinity stress conditions were considered in our analysis, the sensitivity of the results is expected to be multi-fold higher. Our analysis revealed the differentially expressed genes (DEGs) OsbZIP52 and OsLTP2.5 to be common between leaf and root tissues. These genes were further compared with those of the wild halophytic rice Oryza coarctata expression data in stress conditions to understand the significance of these genes. The OsbZIP52 gene homolog of Oryza coarctata was the only one found to be differentially expressed. The expression level of OsbZIP52 was quantified using RT-qPCR and observed downregulated expression in salt stress in root and leaf tissues of four rice cultivars (2 salt-tolerant and 2 salt-sensitive). Promoter and motif analysis revealed a high number of variations in promoter and motif regions of the gene in IR29 salt-sensitive rice. Expression correlation analysis and Gene Ontology study suggested that OsbZIP52 interacts with genes that are engaged in stress response and participate in stress-responsive pathways. Collectively this study increases our understanding of the differential gene expression in various stress conditions in root and leaf tissues. It also helped identify a critical regulatory transcription factor in assisting the plant in combating salinity stress.
Rice is one of the most important food crops globally, with 154 million hectares under cultivation and humans consuming 85% of the total rice produced. Rice is the primary diet of Asia, where 70% of the world’s 1.3 billion people live (Anon 1997). Around 1.0 billion hectares of salt-affected soils are predicted to exist globally, which significantly hampers rice productivity (Singh 2022). While rice (Oryza sativa L.) is by far the most important crop and the most practical alternative for starting crop production in saline soil (Ismail et al. 2007; Singh et al. 2010), it is projected to suffer the most from salt stress, especially in nations with extensive coastlines (Swaminathan and Kesavan 2012). Rice is one of the most vulnerable crop species to salt stress (Steduto et al. 2012), and often a small amount of 50 mM NaCl can have a negative impact (Hairmansis et al. 2014; Munns and Tester 2008; Yeo and Flowers 1986).
By 2025, the global total rice consumption would be 800 mT, whereas rice demand in Asia is predicted to rise further by about 70% in the next 30 years due to population expansion (Kubo and Purevdorj 2004; Muthayya et al. 2014). This is likely to produce a double blow to producers and consumers alike with the expected rise in salinity stress and a loss in production per unit area. Salinity remediation has its set of constraints, and so rice genetic modification for higher salinity tolerance is critical (Goff 2002).
Salt stress has two critical effects: ionic imbalance and osmotic stress. Membrane damage, increased lipid peroxidation, and the formation of reactive oxygen species are all consequences of these effects (Kumar et al. 2013). Excess Na+ accumulation results in nutritional imbalance, which ultimately reduces the growth and development of plants (Li et al. 2017). However, some plants have evolved methods to cope with the salty environment. The accumulation of suitable osmolytes in the cytosol, the reduction of osmotic potential to maintain water intake, and the increased activity of enzymes that detoxify reactive oxygen species (ROS) (Schmidt et al. 2013) are some of the key mechanisms. The type of downstream response to abscisic acid (Kurotani et al. 2015), Ca2+ binding proteins (Garg et al. 2015), photosynthesis apparatus damage prevention and preservation (Gill et al. 2013), and energy production and compartmentalization also influence the tolerance or sensitivity of a rice genotype. In stress-tolerant cultivars, several catalytic, DNA binding, transcription regulator activities (Hirayama and Shinozaki 2010; Yamaguchi-Shinozaki and Shinozaki 2006), carbon fixation, glycolysis, and metabolic pathways have been reported to be enriched (Walia et al. 2005), and genes encoding enzymes related to these systems are promising targets to improve stress tolerance in economically important crops (Pardo 2010; Umezawa et al. 2006).
Gene response and the relative importance of those in combating salinity stress must be prioritized to perform the candidate gene approach method. Over a thousand research publications on salinity stress in rice and other crops have been published worldwide in the previous decade, and data on differential gene expression studies have been stored in several public databases. This availability of data opens up a plethora of opportunities to screen out candidate gene list responses to salinity stress and use them further for gene editing or modification. There are no environmental, food, or feed safety concerns with the candidate gene strategy because the donor and receivers are both from the rice gene pool. Through meta-analysis, a group of possible candidate salinity tolerance-linked genes can be identified, validated, and used in rice improvement strategies.
A meta-analysis combines multiple scientific studies followed by statistical validation (Normand 1999; Raudenbush et al. 1991). The challenges of identifying essential salinity tolerance-linked genes and prioritizing them can be resolved using a meta-analysis approach. While meta-analysis regarding abiotic stress such as drought (Khowaja et al. 2009; Shaik and Ramakrishna 2013; Swamy et al. 2011; Trijatmiko et al. 2014), cold or heat stress (Raza et al. 2020), and biotic stresses like rice blast (Ballini et al. 2008) and bacterial blight (Shaik and Ramakrishna 2013) of rice has been implemented, a meta-analysis of salinity stress of rice has not. A platform for rice meta-analysis is now available (McLaren et al. 2005), and thus the available datasets need to be investigated appropriately to find out the candidate genes responsible for salinity tolerance. There are not many studies that performed a wide range of meta-analyses considering different tissue types and different stress conditions for rice. This research used an integrated meta-analysis approach and discovered crucial genes implicated in salinity response. Both microarray and RNASeq analysis methods were implemented, which made the analysis more robust and sensitive. This study considered different tolerant and sensitive rice varieties and various salt stresses to pinpoint the most potent genes expressed in leaf and root tissues responding to salt stress. Various bioinformatics analyses were performed to analyze the effect of potential genes in response to salt stress. The result of this novel study will aid in understanding the salt stress tolerance mechanism by looking at the different rice varieties and stress conditions and coming to a common ground regarding the genes that are likely responsible for stress tolerance.
Materials and methods
Salinity stress data were retrieved from GEO: Gene Expression Omnibus (Barrett et al. 2012) (www.ncbi.nlm.nih.gov/geo/). GSE76613 (Kong et al. 2019), GSE58603 (Wang et al. 2016), and GSE21651 (Mishra et al. 2018) datasets of salt stress experiments were chosen for Microarray data analysis. The RNASeq data (Illumina and 454) of rice salinity were retrieved from Sequence Retrieval Archive (SRA) at NCBI, and their accessions were GSE109341 (Formentin et al. 2018). The count data obtained from NCBI were all already preprocessed and directly utilized for differential expression analysis. These datasets were employed in a rice salinity meta-analysis where independent differtential expression analysis was performed on each dataset.
Datasets were categorized into two types: leaf and root. Root data includes 3 japonica varieties (TCN1: salt-sensitive, Baldo: salt-tolerant, and Vialone Nano: salt-sensitive) and 3 indica varieties (TNG67: salt-tolerant, PL177: salt-tolerant, and IR64: salt-sensitive), whereas leaf data encompasses of 4 indica varieties (PL177: salt-tolerant, IR64: salt-sensitive, CSR11: salt-tolerant, and VSR156: salt-sensitive) and 2 japonica varieties (Baldo: salt-tolerant, and Vialone Nano: salt-sensitive). The datasets, along with genotypes, treatment, and stress concentration, were recorded in Table 1.
Identification of differentially expressed genes (DEGs) using microarray datasets
Microarray data analysis for leaf tissue was performed by GEO2R, an available interactive web tool at NCBI. Using the Bioconductor project’s GEOquery and limma R packages (Ritchie et al. 2015), GEO2R compares original submitter-supplied processed data tables. Bioconductor is an open-source software project of R programming language that provides tools for analyzing high-throughput genomic data. GEOquery parses GEO data into R data that other R tools can use. The limma (Linear Models for Microarray Analysis) package is one of the most widely used statistical tests for identifying differentially expressed genes. It can handle a wide range of experimental designs and data sources, as well as apply multiple-testing corrections to P-values to assist in eliminating false-positive results. GEO2R analysis was performed using default parameters where Benjamini–Hochberg procedure controls false discovery rate (FDR) and GEO2R automatically performs log2-transformation on the values in log space and quantile normalization on expression data to have an identical distribution of all samples. The rank of differentially expressed genes was based on log fold change > 1 or < − 1 and with a threshold P-value < 0.05. So, genes with a log fold change greater than 1 or less than − 1 were considered differentially expressed genes (Hong and Breitling 2008). Independent expression analysis was performed on each dataset and different numbers of DEGs were identified in each experiment of the microarray data. They were combined to identify the most significant genes between all microarray data experiments.
Identification of differentially expressed genes (DEGs) using RNASeq datasets
The raw read count data was collected for the RNASeq dataset from Sequence Retrieval Archive (SRA) at NCBI. Differentially expressed genes (DEGs) identification analyses were then performed using the DESeq2 method in R programming language by comparing the read counts of the transcripts of the control and salt-treated samples. DESeq2 analysis performs three consecutive steps to identify differentially expressed genes, such as size factor estimation, dispersion estimation, and FDR statistics. DESeq2 uses the median of ratios method to normalize the count data and employ Benjamini–Hochberg FDR (Benjamini and Hochberg 1995) statistics to screen out differentially expressed genes. DESeq2 analysis was performed in a default setting where default normalization method, size factor, wald test and a negative binomial generalized linear model (GLM) were used to assess differential expression. For differentially expressed genes, a log fold change greater than 1 or less than − 1 and adjusted P-value < 0.05 were considered as cut-off values to preserve consistency over the whole experiment.
Differentially expressed genes in leaf and downstream analysis
Microarray and RNASeq analysis under leaf category were performed, and the genes of the two analyses were compared. Genes from both microarray and RNASeq analysis that passed the cut-off value were considered as significantly differentially expressed genes. The shared genes between the two analyses were selected for further analysis. To predict the relationship between differentially expressed genes, protein–protein interactions were analyzed in Cytoscape software using the STRING v11.0 database (Shannon 2003; Szklarczyk et al. 2019). The RAP (Rice Annotation Project) ID of DEGs in leaf was submitted, and a PPI network was constructed with the default confidence level of 0.4 utilizing STRING database. Gene ontology enrichment study was performed to identify the enriched “Biological process” in response to salt stress. The analysis was also performed using the STRING web tool which utilizes rank-based enrichment detection algorithms and visualized through the R programming language and Cytoscape software.
Differentially expressed genes in root and downstream analysis
Differentially expressed genes from microarray and RNASeq analysis of root were compared, and the shared genes between the two studies were filtered out. These shared genes were then subjected to a protein–protein interaction network (PPI) and gene ontology (GO) enrichment analysis. Cytoscape software was used to perform PPI network analysis by using the STRING database (https://version-11-0b.string-db.org/), and the GO Biological process was also observed from the STRING v11.0 web tool. A default condifence score of 0.4 was used for GO analysis. For clear visualization, R and Cytoscape tools were used.
DEGs between leaf and root, and their significance
Resultant differentially expressed genes from leaf and root were examined to distinguish the shared genes between the two tissues. In addition, log fold changes of common genes in different rice varieties from the experiments, as well as their involvement in protein networks and biological processes through the STRING database, were observed.
Comparison of DEGs with halophyte Oryza coarctata
Oryza coarctata (Oryza coarctata) is a wild halophyte that has adapted to high salinity (20–40 dS ml−1) and can withstand long periods of complete submersion in saline water. The differential genes for Oryza were obtained from a previously performed transcriptomic (Garg et al. 2014). The study provides DEGs from different stress conditions with at least twofold changes. A comparison between the differential expression of genes obtained from the above analysis and the differential expression of genes from Oryza coarctata was performed.
Real-time qPCR and data anlaysis
RT-qPCR analysis was carried out to confirm that the bZIP52 gene is differently expressed between salt-stressed and control samples in both salt-sensitive and salt-tolerant genotypes. Two salt-sensitive cultivars (BR28 and IR29) and two salt-tolerant cultivars (Horkuch and Pokkali) were selected for this purpose. The experiment was performed in the net house of the Plant Biotechnology Laboratory, University of Dhaka. Seeds from these cultivars were incubated at 50 °C for 5 days to break their dormancy. Seeds soaked in distilled water were then placed in a 37 °C incubator for germination for 3 days. The seeds were then transferred to Styrofoam sheets floating in two trays containing Yoshida’s solution. Yoshida’s solution was changed every 2d. After 14 days, 60 mM salt stress was applied to one tray; no salt was applied to an identical control tray. Then the salt concentration of the salt stress treatment tray was increased by 20 mM each day until it reached 100 mM salt concentration. Seedlings were randomly taken from both trays after 24 h of 100 mM salt application. Samples were washed and dried with tissue paper. Shoot and root tissues were collected separately in liquid nitrogen. Therefore, the total number of samples was [4 cultivars × 2 treatments × 3 biological replicates × 2 tissues = 48].
Total RNA was isolated from the collected tissues using TRIZOL following the manufacturer’s protocol (Invitrogen, USA) and quantified using Qubit™ 4 Fluorometer (Thermo Fisher Scientific Inc.). Total RNA purity and degradation were again evaluated on 1% agarose gels before proceeding to the RT-qPCR step. 1.5 μg samples of purified total RNA from all four cultivars were reverse transcribed into cDNA using SuperScript™ III First-Strand Synthesis System for RT-qPCR (Invitrogen, USA). RT-qPCR was performed in triplicate using SYBRTM Select Master Mix for CFX on CFX Opus 96 Real-Time PCR System (Bio-Rad) according to the manufacturer’s recommendations. The primers were designed using the Primer3Plus program and are listed in Additional file 6: Table S6. The thermal cycling conditions were as follows: 95 °C for 3 min followed by 40 cycles of 95 °C for 10 s and 62 °C for the 30 s. The specificity of amplification was verified at the end of each PCR run using the melting curve data. To determine the relative fold differences for each sample, the Ct value for each gene was normalized to the Ct value of the sno gene, and the difference between cultivars was calculated using the 2 − ΔΔCt method (Livak and Schmittgen 2001). Data analysis was performed in R using multiple T-test statistical analysis where P-value < 0.05 was considered significant against control.
Genetic variant analysis
The variant analysis was performed on the promoter and motif regions of the OsbZIP52 gene between salt-tolerant (Horkuch and Pokkali) and salt-sensitive varieties (IR29). A Meme suite server was utilized to find the transcription factor binding regions within the promoter of the gene and examined the variations between tolerant and sensitive varieties. It utilizes an expectation–maximization algorithm to iteratively build position-specific probability matrices (PSPMs) representing the motifs. The MEME algorithm also incorporates background models to estimate motif significance (Bailey et al. 2015). For this analysis, promoter sequence of OsbZIP2 was entered into the Meme suite server and as a result, the server provided with transcription factor binding regions within the promoter. Later, multiple sequence alignment was used to those particular binding sites from both tolerant and sensitive varieties in order to see how the transcription factor binding sequences varied. Moreover, motif regions from coding sequences of the gene from tolerant and sensitive varities were sreached using the Tomtom algorithm and subjected to multiple sequence alignment between varieties. Tomtom employs a similarity measure, such as the Pearson correlation coefficient or the Kullback–Leibler divergence, to assess the similarity between motifs and identify potential matches. (Gupta et al. 2007). Tomtom allows to compare a query motif against a database of known motifs to find potential matches and infer functional annotations.
The significant differentially expressed gene was then subjected to co-expression analysis. The co-expression analysis was performed in the genevestigator software. Genevestigator (Hruz et al. 2008) application contains gene expression databases and tools for analyzing expression data. The co-expression analysis was conducted by analyzing samples of stress-related microarray datasets for tolerant plants with a maximum threshold value of 0.9 and P-value < 0.05. Using the euclidian distance approach, the Hierarchical clustering algorithm depicted pattern similarities of co-expressed genes in different stress expression data based on log fold change in transcript abundance with the largest difference from control.
Expression analysis in different plants
The expression pattern of the most significant gene was then observed in different plants, such as Arabidopsis thaliana (Accession number: GSE27548, GSE27550, and GSE53990), Medicago truncatula (Accession number: E-MTAB-2681, and GSE13907), Glycine max (Accession number: PRJNA516324, E-MTAB-4352, and PRJNA306380), and Zea mays (Accession number: PRJNA556806, PRJNA657262, and E-MTAB-4258), under different stress conditions using genevestigator software. The software selects samples from all available microarray or RNASeq datasets from different databases, including NCBI, array express, and many others, for each plant based on different stress conditions and performs expression analysis using a log fold change cut-off of 1.
The analysis focused on identifying genes consistently regulated across different studies, providing insights into the shared molecular mechanisms underlying salt stress response in rice. The complete structure of the research is represented in Fig. 1.
Microarray and RNASeq analysis result of leaf tissue
The expression analysis was performed on two leaf microarray dataset and obtained 767 statistically significant DEGs. One leaf RNASeq analysis revealed 650 statistically significant DEGs. The comparison between microarray and RNASeq analysis results provided 106 DEGs (Additional file 1: Table S1) common in both leaf tissue analyses in different rice varieties (Fig. 2A). The expression pattern of these 106 genes was investigated through a heatmap (Fig. 2B), and two clusters of genes were observed across the pattern expressing in the opposite manner. The general differential expression pattern appeared to be quite similar between tolerant and sensitive plants, but the heat map revealed a small variation that could be attributable to the stress treatment time, as each data set employed varied treatment periods such as 3 h, 24 h, 3 days, and so on. Also, some salt stress genes were found expressed differently in tolerant and sensitive plants, including OsTPS2 (improve tolerance under cold and salinity stress), OsLEA14A (improve tolerance against dehydration, salinity, CuSO4, and HgCl2) (Hu et al. 2019), OsMST6 (enhances membrane integrity and stability under salt stress) (Monfared et al. 2020), OsHsfB2b (acts as a negative regulator in response to salinity stress) (Xiang et al. 2013), and many others.
Protein–Protein network and gene ontology analysis of leaf DEGs
The interactions between the predicted 106 differentially expressed leaf proteins were investigated through the STRING web tool. 33 proteins were depicted in the network that showed interaction with one or more proteins (Fig. 2C). OsJ_04024 (putative heat shock protein 70 kDa) protein was identified to interact with the highest number of proteins (8 proteins) in the network. The network also includes some other salt stress-related genes, including OsPYL4, Os01g0337500, OsBZ8, OsTPS2, OsSRFP1, and many others. Gene ontology analysis to explore enriched biological processes revealed that these DEGs were significantly enriched in 15 biological processes (Additional file 2: Table S2) terms (FDR < 0.04) (Fig. 2D), including terpenoid metabolic process (GO:0006721)(4 genes), primary metabolic process (GO:0044238)(16 genes), cellular metabolic process (GO:0044237)(15 genes), and organic substance biosynthetic process (GO:1901576)(10 genes).
Microarray and RNASeq analysis result of root tissue
The GEO2R analysis of two root microarray datasets provided 82, and DESeq2 analysis of one root RNASeq studies unraveled 754 DEGs. Further observation of the results of these two analyses identified 29 common (Fig. 3A) DEGs (Additional file 3: Table S3). The expression pattern of these 29 genes was observed via heatmap, and the dendrogram divided all the genes into 2 clusters (Fig. 3B). The overall differential expression pattern looked quite similar between tolerant and sensitive, but the subtle difference that was observed from the heat map are due to stress response genes and different stress treatment procedure as each data used different treatment periods and conditions (Table 1). Various salt stress genes were discovered from the analysis that expressed differentially in root tissue, such as OsHOX24 (enhances susceptibility to abiotic stresses through fine-tuning salt-responsive gene expression in rice) (Bhattacharjee et al. 2017), OsHOX22 (functions as a negative regulator in salt and drought tolerance in rice) (Zhang et al. 2012b), OsPP2C30 (regulates abiotic stress via ABA signaling pathways) (Singh et al. 2015), and Os03g0757600 (response to abiotic stress through glycosylation of flavonoids) (Dong et al. 2020).
protein–protein network and gene ontology analysis of root DEGs
The relationship between predicted DEGs in root tissue, protein–protein interaction network was constructed using the STRING database. Only 5 of the genes were depicted in the network that interacts with one another (Fig. 3C), and OsJ_009875 (probable protein phosphatase 2C 30) interacts with the highest (2 proteins) number of proteins in the network. Gene ontology analysis showed that the predicted DEGs were enriched in 29 biological process terms (FDR < 0.04) (Fig. 3D), among them, regulation of cellular process (GO:0050794) (5 genes), cellular response to acid chemical (GO:0071229) (2 genes), carbohydrate derivative metabolic process (GO:1901135) (3 genes), cellular response to chemical stimulus (GO:0070887) (3 genes), and nitrogen compound metabolic process (GO:0006807) (9 genes) were included (Additional file 4: Table S4).
Significant DEGs in leaf and root tissues and downstream analysis
Differential expression analysis revealed 106 and 29 dysregulated genes, respectively, from leaf and root experiments. Among these genes, two genes (OsbZIP52: b-ZIP transcription factor 52 and OsLTP2.5: non-specific lipid transfer protein 2.5) were shared between the leaf and root that were differentially expressed. Log fold change of OsbZIP52 gene in different rice varieties across the analysis was observed (Fig. 4A). The dysregulation of OsbZIP52 was found to vary based on the plant's geographical distribution and salt response in different tissues. STRING network analysis showed that OsbZIP52 forms a protein–protein network composed of 11 significant proteins (Fig. 4B). These proteins were enriched in 13 critical biological processes (Additional file 5: Table S5) terms (FDR < 0.03) (Fig. 4C), including signal transduction (GO:0007165) (2 genes), cellular response to an organic substance (GO:0071310) (2 genes), response to abscisic acid (GO:0009737) (2 genes), response to stress (GO:0006950) (3 genes), and response to oxygen-containing compound (GO:1901700) (4 genes). On the other hand, OsLTP2.5 was also found to be dysregulated based on the regional distribution of plants and tissues (Fig. 5A). OsLTP2.5 formed a network of 9 proteins (Fig. 5B), but no significant GO terms have been found. However, recent studies found that the involvement of nsLTP in abiotic stress responses may aid plant adaptation to change environmental circumstances (Liu et al. 2015). Therefore OsbZIP52 and OsLTP2.5 were of great significance in tolerating stress conditions.
Analysis of resultant DEGs with Oryza coarctata expression analysis
The significant genes were then compared with Oryza differential expression analysis in different stress conditions. Transcriptomic analysis in different stress conditions of O. coarctata has been performed in past experiments (Garg et al. 2014). The expression data of that experiment revealed that the bZIP9 gene of Oryza with the ‘Arabidopsis’ annotation of basic leucine zipper 9 and a homolog to the OsbZIP52 rice gene (HomoloGene:116482) was found to be differentially expressed in submergence plus salt stress condition. Such a response of a specific gene to salt stress in a highly salt-tolerant plant clearly depicted the significance of the gene.
Expression analysis of OsbZIP52 genes by RT-qPCR
The expression profile of the OsbZIP52 gene in the root and shoot of tolerant and sensitive plants was studied to determine the physiological and functional importance of the gene during salt stress. It was observed that the gene was significantly downregulated under salt stress compared to control conditions in roots of both salt-tolerant (Horkuch, Pokkali) and salt-sensitive cultivars (BR28, IR29) (Fig. 6A). A similar result was found for shoot tissue but the decrease in expression of the gene was only significant for Horkuch but not in Pokkali, IR29, and BR28 (Fig. 6B). So there appeared to be no difference in the response of OsbZIP52 to salt stress between salt sensitive and tolerant geneotypes.
Variant analysis between tolerant and sensitive varities
Promoter and gene sequence of OsbZIP52 gene of two salt-tolerant and one salt-sensitive genotypes were obtained from Plant Biotechnology Laboratory, University of Dhaka. After analyzing the OsbZIP52 promoter regions of two salt-tolerant verities (Horkuch and Pokkali) and one salt-sensitive verity (IR29), a greater number of variations was observed in IR29 where different abiotic stress-responsive transcription factors bind which act by increasing expression of genes involved in ion channel (AP2-ERF family) (Xie et al. 2019), by regulating alternative splicing (ARF family) (Ye et al. 2020), and by controlling Na+ and K+ transport systems (bHLH and WRKY) (Rajappa et al. 2020) (Fig. 7 and Additional file 7: Table S7). The differences in sequence variation between tolerant and sensitive varieties was also visualized in Fig. 7 As a result, regulation in the expression of the OsbZIP52 gene in salt-sensitive varieties may be altered due to these variations in the promoter regions during abiotic stress including salt stress.
Using the MEME suit 7 motif regions were predicted for the salt-tolerant varieties. These motifs correspond to a number of domains of different transcription factor families which play a role in response to abiotic stresses. There were many variations in the genomic motif regions of salt-sensitive verity (IR29) found by performing multiple sequence alignment using horkuch, pokkali and IR29 varieties, which may indicate a lack in sensitivity of the functionality of OsbZIP52 TF during abiotic stresses (Additional file 8: Table S8).
Co-expression analysis of the significant gene
To further assess the significance of the OsbZIP52 gene, a co-expression analysis using the genevestigator software was performed. Tolerant plants (FL478, PL177, and Pokkali) expression data were used from GSE13735, GSE58603, and GSE14403 experiments to perform co-expression analysis. Co-expressional analysis showed different salt stress-response genes were negatively co-related with OsbZIP52 (LOC_Os06g45140) gene, including OsLEA3:Late Embryogenesis Protein (LOC_Os05g46480) (Duan and Cai 2012), RAB16A protein (LOC_Os11g26790) (Ganguly et al. 2012), glutathione-S-transferase (LOC_Os01g49710) (Sharma et al. 2014), copper/zinc superoxide dismutase like protein (LOC_Os07g46990) (Prashanth et al. 2008), and DEAD-box helicase ATP-binding protein (LOC_Os11g38670) (Macovei et al. 2012). The expression pattern of these genes, along with OsbZIP52, was visualized using the Hierarchical clustering method (Fig. 4D). OsbZIP52 (LOC Os06g45140) was found to be down-regulated, whereas the other salt response genes were found to be up-regulated, according to the heatmap. Thus, downregulation of OsbZIP52 in tolerant plants may lead to increased expression of other salt response genes, resulting in salt tolerance.
Expression pattern observation in different plants
The expression of the OsbZIP52 gene in rice and other plants was studied under various stress conditions. BZIP9 gene from Arabidopsis thaliana, GLYMA_17G188500 gene from Glycine max, MTR_1g008990 gene from Medicago truncatula, and Zm00001d003529 gene from Zea mays were selected for expression analysis as they are similar to OsbZIP52 gene. Genevestigator software was used in order to analyze the expression pattern of these genes, and the software used available stress-related microarray experiments for each plant. The analysis revealed that genes from each plant were differentially expressed with a log fold change greater than 1. Therefore, each gene responds to abiotic stresses and may be engineered to observe its role in the corresponding plants (Fig. 8).
Salinity stress is a major agricultural limitation that has a global impact on the food supply. Much research has been carried out to decipher the complex systems that operate under stress, but more work is required due to the differences and complexities of the species, genotype, and type of tissue involved. In this study, a meta-analysis was performed using available microarray and RNASeq datasets of root and leaf tissues to find the association of various genes or factors that responds to salt stress. The objective was to find common genes between these two tissues, which may enable us to work on both tissues and overall enhance tolerance. This research considered different tolerant (TNG67, PL177, CSR11, and Baldo) and sensitive (TCN1, IR64, VSR156, and Vialone Nano) rice species with varying salt stresses, which broadened the scope of possible solutions to the salinity problem globally. 106 and 29 differentially expressed genes were discovered in leaf and root tissues, respectively, from the analysis.
In the network of leaf DEGs, OsJ_04024 (heat-shock protein 70, HSP70) was found to have the highest interactions in the network of leaf DEGs, which is a key regulator in response to heat stress (Hu et al. 2009). This protein also responds to cold stress (Lee et al. 2009), dehydration stress (Choudhary et al. 2009), salt stress (Hoang et al. 2015), and some more abiotic stresses (Chankova et al. 2014; Qi et al. 2011). The HSPs family are essential in cell homeostasis, transport of synthesized proteins through cell organelles, and folding, preventing misfolded, denatured, and aggregated proteins induced by stress (Balchin et al. 2016; Ratajczak et al. 2009; Tyedmers et al. 2010). Some essential genes were also discovered from the leaf DEGs. OsPYL4 (abscisic receptor PYL4) is a PYLs family protein that promotes the activation of ABA-regulated genes and results in the acquisition of abiotic stress resistance (Ma et al. 2009; Melcher et al. 2009; Tian et al. 2015). OsBZ8 (Abscisic acid Responsive Element (ABRE)-binding factor) is an important trans-acting factor that regulates NaCl-stress-induced gene expression. This protein strongly interacts with an ABRE-based promoter and regulates abiotic stress-inducible genes (Mukherjee et al. 2006). Ring-type E3-Ubiquitin Ligases (such as OsSRFP1) were reported to interact with the transketolase enzyme to induce NADPH production under salinity stress, which is important in combating ROS-induced damage under abiotic stresses. OsSRFP1 is an unusual ligase protein that acts as a negative regulator under cold stress and a positive regulator under heat stress (Melo et al. 2021; Tunc-Ozdemir et al. 2009). Monosaccharide transporters (OsMST6) are proteins that transport xylose, glucose, and galactose across the hydrophobic membrane and contribute to abiotic stress tolerance (Monfared et al. 2020; Wang et al. 2008). Germin-like proteins (OsGLP8-12) have been reported to function as auxin receptors (Yin et al. 2009) and superoxide dismutase activity which defends against various stresses (Breen and Bellgard 2010). GLPs have the ability to convert superoxide to H2O2 and CO2 while also strengthening the cell wall through protein coupling and glycosylation (Li et al. 2016; Rietz et al. 2012). Therefore differential expression analysis of leaf tissue unravels several crucial genes that actively respond to abiotic stresses and regulate various stress-related pathways. Leaf DEGs were also enriched in important GO biological process terms, such as GO:0006721 (terpenoid metabolic process) (Qiu et al. 2008), GO:1901576 (organic substance metabolic process) (Li et al. 2020), GO:0016998 (cell wall macromolecule catabolic process) (Kim et al. 2015), and many others that are responsive to stress.
OsJ_009875 (probable protein phosphatase 2C 30; OsPP2C30) was the most interacted protein in the Root DEGs network, and a previous study revealed that the protein involves in the ABA-signaling pathway as well as response to abiotic stresses (Singh et al. 2015). Several root DEGs were found that play a crucial role in salinity stress and stress-responsive pathways. Nitrate transporters (OsNRT2.3) form a bridge between hormones such as abscisic acid, ethylene, cytokines, etc., which respond to stress conditions (Krouk 2016). These transporters are also involved with nitrate and k+ or cl− shuttling, and loss of its function prevents plant growth under salinity stress (Drechsler et al. 2015; Fan et al. 2017; Hsu and Tsay 2013; Taochy et al. 2015). OsHox22 and OsHox24 genes that were found to be differentially expressed in root tissue have been reported to have a crucial role in response to abiotic stresses (Bhattacharjee et al. 2016; Jain and Khurana 2009; Zhang et al. 2012b). The network proteins were also seen enriched in salt-responsive GO biological process terms, including GO:0005975 (carbohydrate metabolic process), GO:0050794 (regulation of cellular process), GO:0006807 (nitrogen compound metabolic process) (Jadamba et al. 2020; Wang et al. 2012; Yan et al. 2005; Zhang et al. 2012a). By observing GO terms of both leaf and root tissues, 4 GO Biological processes were found common, such as the primary metabolic process (GO:0044238), organic substance metabolic process (GO:0071704), cellular metabolic process (GO:0044237), and organic substance biosynthetic process (GO:1901576).
After analyzing differentially expressed genes from leaf and root tissues, OsbZIP52 and OsLTP2.5 genes were common. To assess the significance of these genes, downstream bioinformatics analyses were undertaken. The OsbZIP52 logFC observation revealed that gene expression might vary based on tissue types and geographic distribution. The action mechanism of this gene, response to salt, may differ between tissues as we observed relative difference in expression between tissues from Fig. 6, and the regional distribution of the plants (such as Baldo and Vialone Nano are small-grained risotto rice from Europe (Cirillo et al. 2009), whereas CSR11 and VSR156 rice genotypes are from Asia (Mishra et al. 2018) may also contribute to this variation of expression. PPI network analysis of the OsbZIP52 gene revealed the association with essential genes (RAB25, BZIP50, Os03g0158200) that are highly responsive to stress, and the proteins of the network were enriched in GO Biological pathways such as response to stress (GO:0006950), response to abscisic acid (GO:0009737), response to hormone (GO:0009725), signal transduction (GO:0007165), and many other stress-responsive biological processes. OsbZIP52 gene was greatly induced by low temperatures and responded to drought, heat, and other abiotic stresses (Liu and Chu 2015; Liu et al. 2012; Zhang et al. 2013). A PPI network was also constructed for OsLTP2.5, a non-specific lipid transfer protein, but no enriched GO terms were found. A previous study suggested that non-specific lipid transfer proteins were crucial in resisting biotic and abiotic stresses (Choi et al. 2008; Jang et al. 2002; Jung et al. 2003; Kielbowicz-Matuk et al. 2008; Liu et al. 2012; Molina and Garcia-Olmedo 1993). They are also vital to plant growth and development. Therefore, OsLTP2.5 may play an essential role in response to salt stress and become a target for further genetic research to resist salinity.
Oryza coarctata expression analysis with resultant DEGs revealed that the basic leucine zipper 9 annotated gene, a homolog of the OsbZIP52 gene, was differentially expressed in submergence and stress conditions. Oryza coarctata is a plant of high salinity tolerance which can tolerate up to 40 dS ml−1 salt stress (Garg et al. 2014). Differential expression of basic leucine zipper 9 in Oryza coarctata under stress conditions explains the significance of the gene to tolerate salinity stress. Consequently, the OsbZIP52 gene is anticipated to play an important role in response to salt stress. The dysregulation of the OsbZIP52 gene in different rice varieties considered by the analysis was examined, and it was observed that the up and down-regulation varied based on species and tissue types.
bZIP transcription factors are key regulators in plant salt stress responses. They perceive and signal salt stress, regulating the expression of stress-responsive genes. These transcription factors control ion homeostasis by modulating ion transporters and channels (Liu et al. 2023). They also regulate osmotic adjustment by influencing the synthesis and accumulation of compatible solutes (Yang et al. 2019). Additionally, bZIP factors play a role in antioxidant defense by regulating genes involved in scavenging reactive oxygen species (Dvořák et al. 2021). They interact with hormone signaling pathways, particularly abscisic acid, affecting stomatal closure and other hormonal responses (Schlögl et al. 2008). There isn’t any clear indication about the regulation mechanism of bzip transcription factor under salt stress and so further analysis is required to understand the complex mechanism of this family protein under salt stress. The OsbZIP52 gene expression and promoter sequence differences between tolerant and sensitive types were found to be significantly altered in this study, which highlights the need for deeper characterization and understanding of proteome interactions for salt-tolerant crop development.
OsbZIP52 expression was found to be downregulated in root and shoot tissues of salt-tolerant and salt-sensitive verities under salt stress in the RT-qPCR analysis. There were no differences in the response between tolerant and sensitive varieties in terms of gene expression. Since the RT-qPCR analysis was done for only a one-time point of 24 h, it is not enough to get a comprehensive understanding. In addition, the variant analysis showed a significant alteration in the promoter and motif regions of the OsbZIP52 gene in salt-sensitive varieties. Although there wasn’t a significant difference observed between salt-tolerant and salt-sensitive varieties in a time point of 24 h, extensive analysis may be performed to perceive the behaviour of the gene expression under salt stress at different time points in order to understand the role of the OsbZIP52 gene in defending against salt stress in tolerant genotypes.
The co-expression analysis of OsbZIP52 genes was performed to identify the co-expressing genes associated with salt stress response using tolerant plant data in genevestigator software. In comparison to the downregulation of OsbZIP52, multiple salt response genes were up-regulated, such as OsLEA3 and RAB16A, genes that enhance salt tolerance, copper/zinc superoxide dismutase, glutathione S transferase, and DEAD-box helicase ATP-binding protein are such salt response proteins that play a critical role in response to salinity stress. This result was also found consistent with the previously performed test (Liu et al. 2012), where it was observed that OsbZIP52 over-expression cell lines showed a significant decrease in the expression of several salt-stress genes, including OsLEA3, LOC_Os05g39250, Rab25, etc. Therefore, OsbZIP52 down-regulation may enhance salinity tolerance.
Expression pattern observation in Oryza sativa, Arabidopsis thaliana, Medicago truncatula, Glycine max, and Zea mays plants for OsbZIP52 revealed that the gene was differentially expressed in each plant in different tissues and stress conditions. As this gene was found to respond to different abiotic stresses in various plant species, this gene may become a common point to increase the tolerance of a wide range of plant species.
Following substantial expression investigation employing Microarray and RNASeq analytical techniques, the OsbZIP52 gene was found to be the most crucial gene, having important interactions with other essential genes and participating in salt-responsive pathways. Compared to OsbZIP52, which is downregulated, the co-expression analysis found that salt stress-responsive genes up-regulated in tolerant plants. The study considered a wide range of rice plants along with different salt concentrations, which helped to understand the significance of the OsbZIP52 gene better. Altering the expression of this gene may control the growth of different rice varieties in any saline condition. Gene-specific study of the OsbZIP52 gene may establish potential evidence behind its role against salt stress. The findings of the study can become a stronghold in future research studies and assist in the development of the salt-tolerant transgenic breed.
The research performed meta-analyses on several salinity-related microarray and RNASeq datasets and filtered out the OsbZIP52 gene to be the most significant gene for salinity response shared between root and leaf. Differences in the response of OsbZIP52 to 24 h of salt stress were not observed among the salt-tolerant and sensitive genotypes. However, there is further scope in fine-tuning this response at different times and levels of stress. The research also showed that OsbZIP52 responds to other abiotic stresses and is differentially expressed in other plant species. Therefore, this gene can be used as a potential target for genetic engineering or mutant analysis for future investigation and will serve as a valuable genetic resource for the development of salt-tolerant rice.
Anon. Human Development Report 1997. 1997.
Bailey TL, Johnson J, Grant CE, Noble WS. The MEME suite. Nucleic Acids Res. 2015;43(W1):W39-49. https://doi.org/10.1093/NAR/GKV416.
Balchin D, Hayer-Hartl M, Hartl FU. In vivo aspects of protein folding and quality control. Science. 2016;353(6294):aac4354. https://doi.org/10.1126/science.aac4354.
Ballini E, Morel JB, Droc G, Price A, Courtois B, Notteghem J-L, Tharreau D. A genome-wide meta-analysis of rice blast resistance genes and quantitative trait loci provides new insights into partial and complete resistance. Mol Plant-Microbe Interact. 2008;21(7):859–68. https://doi.org/10.1094/MPMI-21-7-0859.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, Yefanov A, Lee H, Zhang N, Robertson CL, Serova N, Davis S, Soboleva A. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2012;41(D1):991–5. https://doi.org/10.1093/nar/gks1193.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B (methodol). 1995;57(1):289–300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x.
Bhattacharjee A, Khurana JP, Jain M. Characterization of rice homeobox genes, OsHOX22 and OsHOX24, and over-expression of OsHOX24 in transgenic arabidopsis suggest their role in abiotic stress response. Front Plant Sci. 2016. https://doi.org/10.3389/fpls.2016.00627.
Bhattacharjee A, Sharma R, Jain M. Over-expression of OsHOX24 confers enhanced susceptibility to abiotic stresses in transgenic rice via modulating stress-responsive gene expression. Front Plant Sci. 2017;8:628. https://doi.org/10.3389/FPLS.2017.00628.
Breen J, Bellgard M. Germin-like Proteins (GLPs) in cereal genomes: gene clustering and dynamic roles in plant defence. Funct Integr Genomics. 2010;10(4):463–76. https://doi.org/10.1007/s10142-010-0184-1.
Chankova SG, Dimova EG, Mitrovska Z, Miteva D, Mokerova DV, Yonova PA, Yurina NP. Antioxidant and HSP70B responses in chlamydomonas reinhardtii genotypes with different resistance to oxidative stress. Ecotoxicol Environ Saf. 2014;101:131–7. https://doi.org/10.1016/j.ecoenv.2013.11.015.
Choi AM, Lee SB, Cho SH, Hwang I, Hur C-G, Suh MC. Isolation and characterization of multiple abundant lipid transfer protein isoforms in developing sesame (Sesamum indicum L.) seeds. Plant Physiol Biochem. 2008;46(2):127–39. https://doi.org/10.1016/j.plaphy.2007.10.003.
Choudhary MK, Basu D, Datta A, Chakraborty N, Chakraborty S. Dehydration-responsive nuclear proteome of rice (Oryza sativa L.) illustrates protein network, novel regulators of cellular adaptation, and evolutionary perspective. Mol Cell Proteomics. 2009;8(7):1579–98. https://doi.org/10.1074/mcp.M800601-MCP200.
Cirillo A, Del Gaudio S, Di Bernardo G, Galderisi U, Cascino A, Cipollaro M. Molecular characterization of Italian rice cultivars. Eur Food Res Technol. 2009;228(6):875–81. https://doi.org/10.1007/S00217-008-1000-1.
Dong NQ, Sun Y, Guo T, Shi CL, Zhang YM, Kan Y, Xiang YH, Zhang H, Yang YB, Li YC, Zhao HY, Yu HX, Lu ZQ, Wang Y, Ye WW, Shan JX, Lin HX. UDP-glucosyltransferase regulates grain size and abiotic stress tolerance associated with metabolic flux redirection in rice. Nat Commun. 2020;11(1):1–16. https://doi.org/10.1038/s41467-020-16403-5.
Drechsler N, Zheng Y, Bohner A, Nobmann B, von Wirén N, Kunze R, Rausch C. Nitrate-dependent control of shoot K homeostasis by NPF7.3/NRT1.5 and SKOR in Arabidopsis. Plant Physiol. 2015. https://doi.org/10.1104/pp.15.01152.
Duan J, Cai W. OsLEA3–2, an abiotic stress induced gene of rice plays a key role in salt and drought tolerance. PLoS ONE. 2012;7(9):e45117. https://doi.org/10.1371/JOURNAL.PONE.0045117.
Dvořák P, Krasylenko Y, Zeiner A, Šamaj J, Takáč T. Signaling toward reactive oxygen species-scavenging enzymes in plants. Front Plant Sci. 2021;11:618835. https://doi.org/10.3389/FPLS.2020.618835/FULL.
Fan X, Naz M, Fan X, Xuan W, Miller AJ, Xu G. Plant nitrate transporters: from gene function to application. J Exp Bot. 2017;68(10):2463–75. https://doi.org/10.1093/jxb/erx011.
Formentin E, Barizza E, Stevanato P, Falda M, Massa F, Tarkowskà D, Novák O, Lo Schiavo F. Fast regulation of hormone metabolism contributes to salt tolerance in rice (Oryza sativa Spp. Japonica, L.) by inducing specific morpho-physiological responses. Plants (basel, Switzerland). 2018;7(3):75. https://doi.org/10.3390/PLANTS7030075.
Ganguly M, Datta K, Roychoudhury A, Gayen D, Sengupta DN, Datta SK. Overexpression of Rab16A gene in indica rice variety for generating enhanced salt tolerance. Plant Signal Behav. 2012;7(4):502–9. https://doi.org/10.4161/PSB.19646.
Garg R, Verma M, Agrawal S, Shankar R, Majee M, Jain M. Deep transcriptome sequencing of wild halophyte rice, porteresia coarctata, provides novel insights into the salinity and submergence tolerance factors. DNA Res Int J Rapid Publ Rep Genes Genomes. 2014;21(1):69–84. https://doi.org/10.1093/dnares/dst042.
Garg R, Narayana Chevala VV, Shankar R, Jain M. Divergent DNA methylation patterns associated with gene expression in rice cultivars with contrasting drought and salinity stress response. Sci Rep. 2015. https://doi.org/10.1038/srep14922.
Gill SS, Tajrishi M, Madan M, Tuteja N. A DESD-box helicase functions in salinity stress tolerance by improving photosynthesis and antioxidant machinery in rice (Oryza sativa L. Cv. PB1). Plant Mol Biol. 2013;82(1–2):1–22. https://doi.org/10.1007/s11103-013-0031-6.
Goff SA. A draft sequence of the rice genome (Oryza sativa L. Ssp. Japonica). Science. 2002;296(5565):92–100. https://doi.org/10.1126/science.1068275.
Gupta S, Stamatoyannopoulos JA, Bailey TL, Noble WS. Quantifying similarity between motifs. Genome Biol. 2007;8(2):1–9. https://doi.org/10.1186/GB-2007-8-2-R24/TABLES/3.
Hairmansis A, Berger B, Tester M, John Roy S. Image-based phenotyping for non-destructive screening of different salinity tolerance traits in rice. Rice. 2014;7(1):1–10. https://doi.org/10.1186/S12284-014-0016-3.
Hirayama T, Shinozaki K. Research on plant abiotic stress responses in the post-genome era: past, present and future. Plant J. 2010;61(6):1041–52. https://doi.org/10.1111/j.1365-313X.2010.04124.x.
Hoang TM, Moghaddam L, Williams B, Khanna H, Dale J, Mundree SG. Development of salinity tolerance in rice by constitutive-overexpression of genes involved in the regulation of programmed cell death. Front Plant Sci. 2015;6:175. https://doi.org/10.3389/fpls.2015.00175.
Hong F, Breitling R. A comparison of meta-analysis methods for detecting differentially expressed genes in microarray experiments. Bioinformatics. 2008;24(3):374–82. https://doi.org/10.1093/bioinformatics/btm620.
Hruz T, Laule O, Szabo G, Wessendorp F, Bleuler S, Oertle L, Widmayer P, Gruissem W, Zimmermann P. Genevestigator V3: a reference expression database for the meta-analysis of transcriptomes. Adv Bioinform. 2008. https://doi.org/10.1155/2008/420747.
Hsu P-K, Tsay Y-F. Two phloem nitrate transporters, NRT1.11 and NRT1.12, are important for redistributing xylem-borne nitrate to enhance plant growth. Plant Physiol. 2013;163(2):844–56. https://doi.org/10.1104/pp.113.226563.
Hu W, Guocheng H, Han B. Genome-wide survey and expression profiling of heat shock proteins and heat shock factors revealed overlapped and stress specific response under abiotic stresses in rice. Plant Sci. 2009;176(4):583–90. https://doi.org/10.1016/J.PLANTSCI.2009.01.016.
Hu T, Liu Y, Zhu S, Qin J, Li W, Zhou N. Overexpression of OsLea14-A improves the tolerance of rice and increases Hg accumulation under diverse stresses. Environ Sci Pollut Res. 2019;26(11):10537–51. https://doi.org/10.1007/S11356-019-04464-Z.
Ismail AM, Heuer S, Thomson MJ, Wissuwa M. Genetic and genomic approaches to develop rice germplasm for problem soils. Plant Mol Biol. 2007;65(4):547–70. https://doi.org/10.1007/S11103-007-9215-2.
Jadamba C, Kang K, Paek N-C, Lee SI, Yoo S-C. Overexpression of Rice Expansin7 (Osexpa7) confers enhanced tolerance to salt stress in rice. Int J Mol Sci. 2020;21(2):454. https://doi.org/10.3390/ijms21020454.
Jain M, Khurana JP. Transcript profiling reveals diverse roles of auxin-responsive genes during reproductive development and abiotic stress in rice. FEBS J. 2009;276(11):3148–62. https://doi.org/10.1111/j.1742-4658.2009.07033.x.
Jang C, Kim D, Bu S, Kim J, Lee S, Kim J, Johnson J, Seo Y. Isolation and characterization of lipid transfer protein (LTP) genes from a wheat-rye translocation line. Plant Cell Rep. 2002;20(10):961–6. https://doi.org/10.1007/s00299-001-0424-x.
Jung HW, Kim W, Hwang BK. Three pathogen-inducible genes encoding lipid transfer protein from pepper are differentially activated by pathogens, abiotic, and environmental stresses*. Plant Cell Environ. 2003;26(6):915–28. https://doi.org/10.1046/j.1365-3040.2003.01024.x.
Khowaja FS, Norton GJ, Courtois B, Price AH. Improved resolution in the position of drought-related QTLs in a single mapping population of rice by meta-analysis. BMC Genom. 2009. https://doi.org/10.1186/1471-2164-10-276.
Kielbowicz-Matuk A, Rey P, Rorat T. The organ-dependent abundance of a solanum lipid transfer protein is up-regulated upon osmotic constraints and associated with cold acclimation ability. J Exp Bot. 2008;59(8):2191–203. https://doi.org/10.1093/jxb/ern088.
Kim S-W, Jeong H-J, Jung K-H. Integrating omics analysis of salt stress-responsive genes in rice. Genes Genom. 2015;37(8):645–55. https://doi.org/10.1007/s13258-015-0293-2.
Kong W, Zhong H, Gong Z, Fang X, Sun T, Deng X, Li Y. Meta-analysis of salt stress transcriptome responses in different rice genotypes at the seedling stage. Plants. 2019;8(3):64. https://doi.org/10.3390/PLANTS8030064.
Krouk G. Hormones and nitrate: a two-way connection. Plant Mol Biol. 2016;91(6):599–606. https://doi.org/10.1007/s11103-016-0463-x.
Kubo M, Purevdorj M. The future of rice production and consumption. J Food Distrib Res. 2004;35(1):128–42. https://doi.org/10.22004/AG.ECON.27145.
Kumar K, Kumar M, Kim S-R, Ryu H, Cho Y-G. Insights into genomics of salt stress response in rice. Rice. 2013. https://doi.org/10.1186/1939-8433-6-27.
Kurotani KI, Hayashi K, Hatanaka S, Toda Y, Ogawa D, Ichikawa H, Ishimaru Y, Tashita R, Suzuki T, Ueda M, Hattori T, Takeda S. Elevated levels of CYP94 family gene expression alleviate the jasmonate response and enhance salt tolerance in rice. Plant Cell Physiol. 2015;56(4):779–89. https://doi.org/10.1093/pcp/pcv006.
Lee DG, Ahsan N, Lee SH, Lee JJ, Bahk JD, Kang KY, Lee BH. Chilling stress-induced proteomic changes in rice roots. J Plant Physiol. 2009;166(1):1–11. https://doi.org/10.1016/J.JPLPH.2008.02.001.
Li L, Xu X, Chen C, Shen Z. Genome-wide characterization and expression analysis of the germin-like protein family in rice and arabidopsis. Int J Mol Sci. 2016;17(10):1622. https://doi.org/10.3390/ijms17101622.
Li Q, Yang A, Zhang WH. Comparative studies on tolerance of rice genotypes differing in their tolerance to moderate salt stress. BMC Plant Biol. 2017. https://doi.org/10.1186/s12870-017-1089-0.
Li Q, Ma C, Tai H, Qiu H, Yang A. Comparative transcriptome analysis of two rice genotypes differing in their tolerance to saline-alkaline stress. PLoS ONE. 2020;15(12):e0243112. https://doi.org/10.1371/journal.pone.0243112.
Liu X, Chu Z. Genome-wide evolutionary characterization and analysis of BZIP transcription factors and their expression profiles in response to multiple abiotic stresses in brachypodium distachyon. BMC Genom. 2015. https://doi.org/10.1186/s12864-015-1457-9.
Liu C, Wu Y, Wang X. BZIP transcription factor OsbZIP52/RISBZ5: a potential negative regulator of cold and drought stress response in rice. Planta. 2012;235(6):1157–69. https://doi.org/10.1007/s00425-011-1564-z.
Liu F, Zhang X, Lu C, Zeng X, Li Y, Fu D, Wu G. Non-specific lipid transfer proteins in plants: presenting new advances and an integrated functional analysis. J Exp Bot. 2015;66(19):5663–81. https://doi.org/10.1093/jxb/erv313.
Liu H, Tang X, Zhang N, Li S, Si H. Role of BZIP transcription factors in plant salt stress. Int J Mol Sci. 2023;24(9):7893. https://doi.org/10.3390/IJMS24097893.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods (san Diego, Calif). 2001;25(4):402–8. https://doi.org/10.1006/METH.2001.1262.
Ma Y, Szostkiewicz I, Korte A, Moes D, Yang Y, Christmann A, Grill E. Regulators of PP2C phosphatase activity function as abscisic acid sensors. Science. 2009. https://doi.org/10.1126/science.1172408.
Macovei A, Vaid N, Tula S, Tuteja N. A new DEAD-box helicase ATP-binding protein (OsABP) from rice is responsive to abiotic stress. Plant Signal Behav. 2012;7(9):1138. https://doi.org/10.4161/PSB.21343.
McLaren CG, Bruskiewich RM, Portugal AM, Cosico AB. The international rice information system. A platform for meta-analysis of rice crop data. Plant Physiol. 2005;139(2):637–42. https://doi.org/10.1104/pp.105.063438.
Melcher K, Ng LM, Zhou XE, Soon FF, Xu Y, Suino-Powell KM, Park SY, Weiner JJ, Fujii H, Chinnusamy V, Kovach A, Li J, Wang Y, Li J, Peterson FC, Jensen DR, Yong E-L, Volkman BF, Cutler SR, Zhu J-K, Xu HE. A gate–latch–lock mechanism for hormone signalling by abscisic acid receptors. Nature. 2009;462(7273):602–8. https://doi.org/10.1038/nature08613.
Melo FV, Oliveira MM, Saibo NJM, Lourenço TF. Modulation of abiotic stress responses in rice by E3-ubiquitin ligases: a promising way to develop stress-tolerant crops. Front Plant Sci. 2021. https://doi.org/10.3389/fpls.2021.640193.
Mishra P, Singh N, Jain A, Jain N, Mishra V, Pushplatha G, Sandhya KP, Singh NK, Rai V. Identification of cis-regulatory elements associated with salinity and drought stress tolerance in rice from co-expressed gene interaction networks. Bioinformation. 2018;14(3):123–31. https://doi.org/10.6026/97320630014123.
Molina A, Garcia-Olmedo F. Developmental and pathogen-induced expression of three barley genes encoding lipid transfer proteins. Plant J. 1993;4(6):983–91. https://doi.org/10.1046/j.1365-313X.1993.04060983.x.
Monfared HH, Chew JK, Azizi P, Xue G-P, Ee S-F, Kadkhodaei S, Hedayati P, Ismail I, Zainal Z. Overexpression of a rice monosaccharide transporter gene (OsMST6) confers enhanced tolerance to drought and salinity stress in Arabidopsis thaliana. Plant Mol Biol Report. 2020;38(1):151–64. https://doi.org/10.1007/s11105-019-01186-x.
Mukherjee K, Choudhury A, Gupta B, Gupta S, Sengupta DN. An ABRE-binding factor, OSBZ8, is highly expressed in salt tolerant cultivars than in salt sensitive cultivars of indica rice. BMC Plant Biol. 2006. https://doi.org/10.1186/1471-2229-6-18.
Munns R, Tester M. Mechanisms of salinity tolerance. Annu Rev Plant Biol. 2008;59:651–81. https://doi.org/10.1146/ANNUREV.ARPLANT.59.032607.092911.
Muthayya S, Sugimoto JD, Montgomery S, Maberly GF. An overview of global rice production, supply, trade, and consumption. Ann N Y Acad Sci. 2014;1324(1):7–14. https://doi.org/10.1111/nyas.12540.
Normand ST. Meta-analysis: formulating, evaluating, combining, and reporting. Stat Med. 1999;18(3):321–59.
Pardo JM. Biotechnology of water and salinity stress tolerance. Curr Opin Biotechnol. 2010;21(2):185–96. https://doi.org/10.1016/j.copbio.2010.02.005.
Prashanth SR, Sadhasivam V, Parida A. Over expression of cytosolic copper/zinc superoxide dismutase from a mangrove plant Avicennia Marina in indica Rice var Pusa Basmati-1 confers abiotic stress tolerance. Transgenic Res. 2008;17(2):281–91. https://doi.org/10.1007/S11248-007-9099-6.
Qi Y, Wang H, Zou Y, Liu C, Liu Y, Wang Y, Zhang W. Over-expression of mitochondrial heat shock protein 70 suppresses programmed cell death in rice. FEBS Lett. 2011;585(1):231–9. https://doi.org/10.1016/j.febslet.2010.11.051.
Qiu D, Xiao J, Xie W, Liu H, Li X, Xiong L, Wang S. Rice Gene Network inferred from expression profiling of plants overexpressing OsWRKY13, a positive regulator of disease resistance. Mol Plant. 2008;1(3):538–51. https://doi.org/10.1093/mplant/ssn012.
Rajappa S, Krishnamurthy P, Kumar PP. Regulation of AtKUP2 expression by BHLH and WRKY transcription factors helps to confer increased salt tolerance to Arabidopsis thaliana plants. Front Plant Sci. 2020;11:1311. https://doi.org/10.3389/FPLS.2020.01311/BIBTEX.
Ratajczak E, Ziętkiewicz S, Liberek K. Distinct activities of Escherichia coli small heat shock proteins IbpA and IbpB promote efficient protein disaggregation. J Mol Biol. 2009;386(1):178–89. https://doi.org/10.1016/j.jmb.2008.12.009.
Raudenbush SW, Hunter JE, Schmidt FL. Methods of meta-analysis: correcting error and bias in research findings. J Am Stat Assoc. 1991;86(413):242. https://doi.org/10.2307/2289738.
Raza Q, Riaz A, Bashir K, Sabar M. Reproductive tissues-specific meta-QTLs and candidate genes for development of heat-tolerant rice cultivars. Plant Mol Biol. 2020;104(1–2):97–112. https://doi.org/10.1007/s11103-020-01027-6.
Rietz S, Bernsdorff FEM, Cai D. Members of the germin-like protein family in Brassica napus are candidates for the initiation of an oxidative burst that impedes pathogenesis of Sclerotinia sclerotiorum. J Exp Bot. 2012;63(15):5507–19. https://doi.org/10.1093/jxb/ers203.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucl Acids Res. 2015;43(7):e47. https://doi.org/10.1093/nar/gkv007.
Schlögl PS, Fábio TS, Nogueira RD, Felix JM, De Rosa VE, Vicentini R, Leite A, Ulian EC, Menossi M. Identification of new ABA- and MEJA-activated sugarcane BZIP genes by data mining in the SUCEST database. Plant Cell Rep. 2008;27(2):335–45. https://doi.org/10.1007/S00299-007-0468-7/FIGURES/4.
Schmidt R, Mieulet D, Hubberten HM, Obata T, Hoefgen R, Fernie AR, Fisahn J, San Segundo B, Guiderdoni E, Schippers JHM, Mueller-Roeber B. SALT-RESPONSIVE ERF1 regulates reactive oxygen species-dependent signaling during the initial response to salt stress in rice. Plant Cell. 2013;25(6):2115–31. https://doi.org/10.1105/tpc.113.113068.
Shaik R, Ramakrishna W. Genes and co-expression modules common to drought and bacterial stress responses in Arabidopsis and rice. PLoS ONE. 2013;8(10):e77261. https://doi.org/10.1371/journal.pone.0077261.
Shannon P. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Sharma R, Sahoo A, Devendran R, Jain M. Over-expression of a rice tau class glutathione S-transferase gene improves tolerance to salinity and oxidative stresses in Arabidopsis. PLoS ONE. 2014;9(3):e92900. https://doi.org/10.1371/JOURNAL.PONE.0092900.
Singh A. Soil salinity: a global threat to sustainable development. Soil Use Manag. 2022;38(1):39–67. https://doi.org/10.1111/SUM.12772.
Singh RK, Redoña E, Refuerzo L. Varietal improvement for abiotic stress tolerance in crop plants: special reference to salinity in rice. In: Abiotic stress adaptation in plants: physiological, molecular and genomic foundation. 2010:387–415. https://doi.org/10.1007/978-90-481-3112-9_18.
Singh A, Jha SK, Bagri J, Pandey GK. ABA inducible rice protein phosphatase 2C confers ABA insensitivity and abiotic stress tolerance in Arabidopsis. PLoS ONE. 2015;10(4):e0125168. https://doi.org/10.1371/journal.pone.0125168.
Steduto P, Hsiao TC, Fereres E, Raes D. Crop yield response to water. FAO irrigation and drainage paper (66). 2012.
Swaminathan MS, Kesavan PC. Agricultural research in an era of climate change. Agric Res. 2012;1(1):3–11. https://doi.org/10.1007/S40003-011-0009-Z.
Swamy BPM, Vikram P, Dixit S, Ahmed HU, Kumar A. Meta-analysis of grain Yield QTL identified during agricultural drought in grasses showed consensus. BMC Genom. 2011. https://doi.org/10.1186/1471-2164-12-319.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, Jensen LJ, von Mering C. STRING V11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucl Acids Res. 2019;47(D1):D607–13. https://doi.org/10.1093/nar/gky1131.
Taochy C, Gaillard I, Ipotesi E, Oomen R, Leonhardt N, Zimmermann S, Peltier J-B, Szponarski W, Simonneau T, Sentenac H, Gibrat R, Boyer J-C. The Arabidopsis Root Stele Transporter NPF2.3 contributes to nitrate translocation to shoots under salt stress. Plant J. 2015;83(3):466–79. https://doi.org/10.1111/tpj.12901.
Tian X, Wang Z, Li X, Lv T, Liu H, Wang L, Niu H, Bu Q. Characterization and functional analysis of pyrabactin resistance-like abscisic acid receptor family in rice. Rice. 2015. https://doi.org/10.1186/s12284-015-0061-6.
Trijatmiko KR, Supriyanta, Prasetiyono J, Thomson MJ, Vera Cruz CM, Moeljopawiro S, Pereira A. Meta-analysis of quantitative trait loci for grain yield and component traits under reproductive-stage drought stress in an upland rice population. Mol Breed. 2014;34(2):283–95. https://doi.org/10.1007/s11032-013-0012-0.
Tunc-Ozdemir M, Miller G, Song L, Kim J, Sodek A, Koussevitzky S, Misra AN, Mittler R, Shintani D. Thiamin confers enhanced tolerance to oxidative stress in Arabidopsis. Plant Physiol. 2009;151(1):421–32. https://doi.org/10.1104/pp.109.140046.
Tyedmers J, Mogk A, Bukau B. Cellular strategies for controlling protein aggregation. Nat Rev Mol Cell Biol. 2010;11(11):777–88. https://doi.org/10.1038/nrm2993.
Umezawa T, Fujita M, Fujita Y, Yamaguchi-Shinozaki K, Shinozaki K. Engineering drought tolerance in plants: discovering and tailoring genes to unlock the future. Curr Opin Biotechnol. 2006;17(2):113–22. https://doi.org/10.1016/j.copbio.2006.02.002.
Walia H, Wilson C, Condamine P, Liu X, Ismail AM, Zeng L, Wanamaker SI, Mandal J, Xu J, Cui X, Close TJ. Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage. Plant Physiol. 2005;139(2):822–35. https://doi.org/10.1104/pp.105.065961.
Wang Y, Xiao Y, Zhang Y, Chai C, Wei G, Wei X, Xu H, Wang M, Ouwerkerk PBF, Zhu Z. Molecular cloning, functional characterization and expression analysis of a novel monosaccharide transporter gene OsMST6 from rice (Oryza sativa L.). Planta. 2008;228(4):525. https://doi.org/10.1007/s00425-008-0755-8.
Wang H, Zhang M, Guo R, Shi D, Liu B, Lin X, Yang C. Effects of salt stress on ion balance and nitrogen metabolism of old and young leaves in rice (Oryza sativa L.). BMC Plant Biol. 2012. https://doi.org/10.1186/1471-2229-12-194.
Wang WS, Zhao XQ, Li M, Huang LY, Jian Long Xu, Zhang F, Cui YR, Bin Ying Fu, Li ZK. Complex molecular mechanisms underlying seedling salt tolerance in rice revealed by comparative transcriptome and metabolomic profiling. J Exp Bot. 2016;67(1):405–19. https://doi.org/10.1093/JXB/ERV476.
Xiang J, Ran J, Zou J, Zhou X, Liu A, Zhang X, Peng Y, Tang N, Luo G, Chen X. Heat shock factor OsHsfB2b negatively regulates drought and salt tolerance in rice. Plant Cell Rep. 2013;32(11):1795–806. https://doi.org/10.1007/s00299-013-1492-4.
Xie Z, Nolan TM, Jiang H, Yin Y. AP2/ERF transcription factor regulatory networks in hormone and abiotic stress responses in Arabidopsis. Front Plant Sci. 2019;10:228. https://doi.org/10.3389/FPLS.2019.00228/BIBTEX.
Yamaguchi-Shinozaki K, Shinozaki K. TRANSCRIPTIONAL REGULATORY NETWORKS IN CELLULAR RESPONSES AND TOLERANCE TO DEHYDRATION AND COLD STRESSES. Ann Rev Plant Biol. 2006;57(1):781–803. https://doi.org/10.1146/annurev.arplant.57.032905.105444.
Yan S, Tang Z, Su W, Sun W. Proteomic analysis of salt stress-responsive proteins in rice root. Proteomics. 2005;5(1):235–44. https://doi.org/10.1002/pmic.200400853.
Yang S, Xu K, Chen S, Li T, Xia H, Chen L, Liu H, Luo L. A stress-responsive BZIP transcription factor OsbZIP62 improves drought and oxidative tolerance in rice. BMC Plant Biol. 2019. https://doi.org/10.1186/S12870-019-1872-1.
Ye Y, Wang J, Wang W, Xu LA. ARF family identification in tamarix chinensis reveals the salt responsive expression of TcARF6 targeted by MiR167. PeerJ. 2020. https://doi.org/10.7717/PEERJ.8829/SUPP-1.
Yeo AR, Flowers TJ. Salinity resistance in rice (Oryza sativa L.) and a pyramiding approach to breeding varieties for saline soils. Aust J Plant Physiol. 1986;13(1):161–73. https://doi.org/10.1071/PP9860161.
Yin K, Han X, Xu Z, Xue H. Arabidopsis GLP4 is localized to the golgi and binds auxin in vitro. Acta Biochim Biophys Sin. 2009;41(6):478–87. https://doi.org/10.1093/abbs/gmp036.
Zhang Z, Liu Q, Song H, Rong X, Ismail AM. Responses of different rice (Oryza sativa L.) genotypes to salt stress and relation to carbohydrate metabolism and chlorophyll content. Afr J Agric Res. 2012a;7(1):19–27. https://doi.org/10.5897/AJAR.9000605.
Zhang S, Haider I, Kohlen W, Jiang L, Bouwmeester H, Meijer AH, Schluepmann H, Liu C-M, Ouwerkerk PBF. Function of the HD-Zip I gene Oshox22 in ABA-mediated drought and salt tolerances in rice. Plant Mol Biol. 2012b;80(6):571–85. https://doi.org/10.1007/s11103-012-9967-1.
Zhang X, Rerksiri W, Liu A, Zhou X, Xiong H, Xiang J, Chen X, Xiong X. Transcriptome profile reveals heat response mechanism at molecular and metabolic levels in rice flag leaf. Gene. 2013;530(2):185–92. https://doi.org/10.1016/j.gene.2013.08.048.
The authors would like to thank Sabrina M Elias from Independent University, Bangladesh for providing the sequence of OsbZIP52 gene of Horkuch, Pokkali and IR29 rice varieties.
This project was supported by a grant from the Biotechnology Research Center, University of Dhaka.
Ethics approval and consent to participate
The authors have no competing interests to declare that are relevant to the content of this article.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Differentially expressed genes in leaf tissue.
List of enriched GO terms from leaf DEGs.
Differentially expressed genes in root tissue.
List of enriched GO terms from root DEGs.
List of enriched GO terms from OsbZIP52 protein network.
Primer used in real-time qPCR.
Promoter region variation analysis by comparing salt-tolerant (Horkuch and Pokkali) and salt-sensitive (IR29) rice varieties.
Motif region variation analysis by comparing salt-tolerant (Horkuch and Pokkali) and salt-sensitive (IR29) rice varieties.
About this article
Cite this article
Chatterjee, D., Shohan, M.U.S., Tamanna, N. et al. Meta-analysis of microarray and RNAseq data reveal OsbZIP52 to mediate salt stress responses in sensitive, tolerant and halophyte rice varieties. CABI Agric Biosci 4, 31 (2023). https://doi.org/10.1186/s43170-023-00173-3