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

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.

Whole-transcriptome profiling and identification of cold tolerance-related ceRNA networks in japonica rice varieties.

Introduction: Low-temperature stress negatively impacts rice yield, posing a significant risk to food security. While previous studies have explored the physiological and linear gene expression alterations in rice under low-temperature conditions, the changes in competing endogenous RNA (ceRNA) networks remain largely unexamined. Methods: We conducted RNA sequencing on two japonica rice varieties with differing cold-tolerance capabilities to establish ceRNA networks. This enabled us to investigate the transcriptional regulatory network and molecular mechanisms that rice employs in response to low-temperature stress. Results: We identified 364 differentially expressed circular RNAs (circRNAs), 224 differentially expressed microRNAs (miRNAs), and 12,183 differentially expressed messenger RNAs (mRNAs). WRKY family was the most prominent transcription factor family involved in cold tolerance. Based on the expression patterns and targeted relationships of these differentially expressed RNAs, we discerned five potential ceRNA networks related to low-temperature stress in rice: osa-miR166j-5p from the miR166 family was associated with cold tolerance; osa-miR528-3p and osa-miR156j-3p were linked to stress response; and osa-miR156j-3p was involved in the antioxidant system. In addition, Os03g0152000 in the antioxidant system, as well as Os12g0491800 and Os05g0381400, correlated with the corresponding stress response and circRNAs in the network. A gene sequence difference analysis and phenotypic validation of Os11g0685700 (OsWRKY61) within the WRKY family suggested its potential role in regulating cold tolerance in rice. Discussion and conclusion: We identified Os11g0685700 (OsWRKY61) as a promising candidate gene for enhancing cold tolerance in japonica rice. The candidate miRNAs, mRNAs, and circRNAs uncovered in this study are valuable targets for researchers and breeders. Our findings will facilitate the development of cold-tolerant rice varieties from multiple angles and provide critical directions for future research into the functions of cold-tolerance-related miRNAs, mRNAs, and circRNAs in rice.