Global analysis of differentially expressed genes between two Japonica rice varieties induced by low temperature during the booting stage by RNA-Seq

As one of the abiotic stresses, low temperature severely threatens rice production during its entire growth period, especially during the booting stage. In the present study, transcriptome analysis was performed comparing Longjing (LJ) 25 (chilling-tolerant) and LJ 11 (chilling-sensitive) rice varieties to identify genes associated with chilling tolerance in rice spikelets. A total of 23 845 expressed genes and 13 205 differentially expressed genes (DEGs) were identified, respectively. Gene ontology (GO) enrichment analyses revealed ‘response to cold’ (containing 180 DEGs) as the only category enriched in both varieties during the entire cold treatment period. Through MapMan analysis, we identified nine and six DEGs related to the Calvin cycle and antioxidant enzymes, respectively, including OsRBCS3, OsRBCS2, OsRBCS4, OsAPX2 and OsCATC, that under chilling stress were markedly downregulated in LJ11 compared with LJ25. Furthermore, we predicted their protein–protein interaction (PPI) network and identified nine hub genes (the threshold of co-expressed gene number ≥ 11) in Cytoscape, including three RuBisCO-related genes with 14 co-expressed genes. Under chilling stress, antioxidant enzyme activities (peroxidase (POD) and catalase (CAT)) were downregulated in LJ11 compared with LJ25. However, the content of malondialdehyde (MDA) was higher in LJ11 compared with LJ25. Collectively, our findings identify low temperature responsive genes that can be effectively used as candidate genes for molecular breeding programmes to increase the chilling tolerance of rice.