Whole-Transcriptome Profiling and Functional Prediction of Long Non-Coding RNAs Associated with Cold Tolerance in Japonica Rice Varieties

Abstract

Low-temperature chilling is a major abiotic stress leading to reduced rice yield and is a significant environmental threat to food security. Low-temperature chilling studies have focused on physiological changes or coding genes. However, the competitive endogenous RNA mechanism in rice at low temperatures has not been reported. Therefore, in this study, antioxidant physiological indices were combined with whole-transcriptome data through weighted correlation network analysis, which found that the gene modules had the highest correlation with the key antioxidant enzymes superoxide dismutase and peroxidase. The hub genes of the superoxide dismutase-related module included the UDP-glucosyltransferase family protein, sesquiterpene synthase and indole-3-glycerophosphatase gene. The hub genes of the peroxidase-related module included the WRKY transcription factor, abscisic acid signal transduction pathway-related gene plasma membrane hydrogen-ATPase and receptor-like kinase. Therefore, we selected the modular hub genes and significantly enriched the metabolic pathway genes to construct the key competitive endogenous RNA networks, resulting in three competitive endogenous RNA networks of seven long non-coding RNAs regulating three co-expressed messenger RNAs via four microRNAs. Finally, the negative regulatory function of the WRKY transcription factor OsWRKY61 was determined via subcellular localization and validation of the physiological indices in the mutant.