Background
Rice is one of the most important cereal crops in the world but is susceptible to cold stress (CS). In this study, we carried out parallel transcriptomic analysis at the reproductive stage on the anthers of two Japonica rice varieties with contrasting CS resistance: cold susceptible Longjing11 (LJ11) and cold resistant Longjing25 (LJ25).
Conclusion
Collectively, our findings provide comprehensive insights into complex molecular mechanisms of CS response and can aid in CS resistant molecular breeding of rice in cold regions.
Results
According to the obtained results, a total of 16,762 differentially expressed genes (DEGs) were identified under CS, including 7,050 and 14,531 DEGs in LJ25 and LJ11, respectively. Examining gene ontology (GO) enrichment identified 35 up- and 39 down-regulated biological process BP GO terms were significantly enriched in the two varieties, with 'response to heat' and 'response to cold' being the most enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 33 significantly enriched pathways. Only the carbon metabolism and amino acid biosynthesis pathways with down-regulated DEGs were enriched considerably in LJ11, while the plant hormone signal transduction pathway (containing 153 DEGs) was dramatically improved. Eight kinds of plant hormones were detected in the pathway, while auxin, abscisic acid (ABA), salicylic acid (SA), and ethylene (ETH) signaling pathways were found to be the top four pathways with the most DEGs. Furthermore, the protein-protein interaction (PPI) network analysis identified ten hub genes (co-expressed gene number ≥ 30), including six ABA-related genes. Various DEGs (such as OsDREB1A, OsICE1, OsMYB2, OsABF1, OsbZIP23, OsCATC, and so on) revealed distinct expression patterns among rice types when the DEGs between LJ11 and LJ25 were compared, indicating that they are likely responsible for CS resistance of rice in cold region.