Identification of drought stress-responsive genes in rice (Oryza sativa) by meta-analysis of microarray data

Meta-analysis provides a systematic access to the previously studied microarray datasets that can recognize several common signatures of stresses. Three different datasets of abiotic stresses on rice were used for meta-analysis. These microarray datasets were normalized to regulate data for technical variation, as opposed to biological differences between the samples. A t-test was performed to recognize the differentially-expressed genes (DEGs) between stressed and normal samples. Gene ontology enrichment analysis revealed the functional distribution of DEGs in different stressed conditions. Further analysis was carried out using software RICE NET DB and divided into three different categories: biological process (homoiothermy and protein amino acid phosphorylation), cellular component (nucleus and membrane), and molecular function (zinc ion binding ad DNA binding). The study revealed that 5686 genes were constantly expressed differentially in Oryza sativa (2089 upregulated and 3597 downregulated). The lowest P value (P = 0.003756) among upregulated DEGs was observed for naringenin, 2-oxoglutrate 3-dioxygenase protein. The lowest P value (P = 0.002866816) among the downregulated DEGs was also recorded for retrotransposon protein. The network constructed from 48 genes revealed 10 hub genes that are connected with topological genes. These hub genes are stress responsive genes that may also be regarded as the marker genes for drought stress response. Our study reported a new set of hub genes (reference genes) that have potentially significant role in development of stress tolerant rice