Knocking Out OsNAC050 Expression Causes Low Temperature Tolerance in Rice by Regulating Photosynthesis and the Sucrose Metabolic Pathway

1. Introduction

Low temperature is one type of environmental stress that has a strong effect on plant growth (Zhang et al. 2022). When the ambient temperatures are continuously lower than the lowest temperature at which plants can grow normally for a period of time, low temperature stress will occur (Yu, Liu, et al. 2022). The abnormal growth of plants caused by low temperatures above 0°C is called chilling injury. Many plants prefer warm temperatures; examples are rice, maize (Zea mays), kidney beans (Phaseolus vulgaris), and tropical and subtropical fruit crops such as bananas and pineapples, which are highly susceptible to cold damage(Guo et al. 2022; Yang et al. 2021). Under normal conditions, when the temperature is below 10°C, chilling damage will occur, which mainly causes discoloration of the plants, necrosis, and formation of lesions on the surfaces of leaves. In the early growth stage of rice seedlings, low temperatures and cold air currents are likely to increase the incidence of bacterial wilt that result in seedling mortality(Wang et al. 2021). However, in the late stage of rice panicle differentiation at the flowering stage, low temperature exposure for an extended period of time will also affect grain yield due to the abnormal development of pollen grains(Liu et al. 2021).

The NAC TF family are known to be related to abiotic stress (Sun et al. 2022). Transgenic rice lines overexpressing SNAC1, SNAC2, and SNAC3 show improved tolerance to salt, while RNAi lines in which these genes are silenced display increased plant susceptibility to these stresses(Zhang et al. 2021; Liu et al. 2014). The NAC-type transcription factor OsNAC2 regulates the expression of ABA-dependent genes and abiotic stress tolerance in rice(Shen et al. 2017). It has been reported that the suppressor line (ONAC095-SRDX) was reduced than wild type (Huang et al. 2016). Overexpressing the NAC transcription factor gene OsY37 promotes leaf senescence and causes early heading in rice (El Mannai et al. 2017). Knockdown of OsNAC5 renders RNAi plants more susceptible to oxidative damage(Takasaki et al. 2010).

With the maturity of gene knockout technology, CRISPR-Cas9 genome editing has rapidly become an important tool in genomic research due to its advantages of conceptually straightforward operation, precise targeting, and relatively high editing efficiency. Therefore, in this study, a member of the NAC transcription factor family, OsNAC050, was identified and knockout mutants were created using CRISPR-Cas9 gene editing. Our results will provide additional information on the functions of a member of the rice NAC transcription factor family.

2. Materials and Methods

2.1. Plant Growth Conditions and Stress Treatments

Plump seeds were chosen from ‘Nipponbare’ (Oryza sativa L. ssp. japonica) and the mutants. The hulls were removed and the seeds were surface-sterilized in 1% NaClO (bleach) for 15 min, followed by washing six times in sterilized distilled water. Fifty seeds were sterilized and placed on MS medium for germination in a light incubator at 28^oC under a 16 h light (3,000 Lux)/8 h dark photoperiod. WT and osnac050 mutant rice plants grown for four weeks (28 d) were selected for low temperature stress treatment, and the seedlings were transferred to a refrigerator at 4°C on day 29. Physiological indicators of the treated plants were assayed at the same times. The plants were sampled at the designated time points with three biological replicates.

2.2. Quantification of relative gene expression by using qRT-PCR

Total RNA was isolated from these samples using the modified Trizol Reagent (Invitrogen, Waltham, MA, USA) following the manufacturer’s protocol. The RNA samples were then subjected to qRT-PCR to quantify OsNAC050 expression in WT plants, using the primers shown in Supplemental Table 1. The rice Actin gene was used as the internal control for normalization of gene expression. The 2–DDCt method was used to calculate relative expression levels. To study the expression rule of OsNAC050 gene expression under different treatments, carry out abiotic stress on it, and take samples at the 0, 1, 3, 6, 9, 12 hour.

2.3. Subcellular localization of OsNAC050

The full length of the gene was cloned into the positioning vector by means of enzyme digestion and connection, and transformed into the rice protoplast(Yoo, Cho and Sheen 2007).：The subcellular localization vector: pBWA(V)HS-ccdb-GLosgfp. Nikon C2-ER laser confocal microscope was used for observation, in which the information of thew2 nuclear marker: Fusion fluorescent protein: mkate, excitation light 561nm, emission light 580nm. Nuclear localization signal protein: nls Amino acid sequence: MDPKKKRKV Fusion fluorescent protein: mKATE excited light at 561nm and emitted light at 580nm(Zhao et al. 2017).

2.3. Targeted Mutagenesis of OsNAC050

OsNAC050 mutants were generated using the CRISPR-Cas9 knock out vector pZHY988. CRISPR-P (http://crispr.hzau.edu.cn/CRISPR2) was used to design the sgRNA (small guide RNA) and to predict potential off-target sites. After the sgRNA was synthesized, the vector was constructed using recombinant DNA methods. Ten positive clones subjected to Sanger sequencing. After detection and DNA sequencing identified the desired vector, rice callus was transformed using Agrobacterium tumefaciens EHA105(Wang et al. 2020). The sequences of primers used to detect off-targets are given in Supplementary Table 1 (Tang et al. 2016; Zhong et al. 2018). Double chain insertion T named osnac050 (T/T). Double chain deletion 5-bp deletion (ATGAG) named osnac050 (-5/-5).

2.4. Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) Analysis

qRT-PCR assays were used to quantify OsNAC050 expression profile analysis and ten pairs of gene-specific primers (Supplementary Table 1) were designed to generate the amplicons used to verify gene expression levels determined using RNA-seq data. We sampled roots, stems, and leaves at the seedling stage (4 weeks), and roots, stems, sheaths, leaves, and ears at the heading stage (10 weeks). After the seeds are vernalized at 4℃, they were planted in soil and grown under specified greenhouse conditions and watered regularly. During sampling, we selected three plants per line at the same growth stage for RNA extraction. To examine OsNAC050 expression levels under various stress and hormone treatments, we subjected four-week-old ‘Nipponbare’ seedlings to cold (4°C), heat (42°C), PEG 6000 (20%, w/v), NaCl (200 mM), H₂O₂ (1%), IAA (100 μM), ABA (100 μM), and GA₃ (100 μM) treatments. Also, qRT-PCR was used to validate the accuracy of the RNA-seq data using 10 randomly selected genes with significantly different expression profiles.Plants were sampled at different time periods, with three replicates for each sampling, and then RNA was extracted. Error bars show the standard deviation based on three independent biological replicates. We used the SYBR Green qPCR kit (Takara) according to the manufacturer's instructions to perform quantitative real-time PCR (qRT-PCR) on an iQ5.0 instrument (Bio Rad, USA). The relative gene expression levels were calculated using the 2^–ΔΔCt method. Quantification was achieved by normalizing the expression to that of the actin and GAPDH genes.

2.5. Physiological Measurements

Seeds of the WT and OsNAC050 mutants were grown in pots to the 4-week seedling stage, followed by the induction of cold stress. Five days later, we measured the physiological indices using our previously described methods [14].

2.6. Transcriptome analysis and data analysis

WT and osnac050(T/T), osnac050(-5/-5) plants were used for RNA-seq analysis. Plants were grown for 10 weeks under normal conditions or for eight weeks under normal conditions and then cold-treated for two weeks. Mixed samples were made by combining three T₂ -generation seedlings. The RNA samples were sequenced and analyzed at the Beijing Institute of Genomics (Shenzhen, China). Genes with p-value < 0.05 and fold change >2 are advised to define differentially expressed gene (DEGs).

3. Results

3.1. OsNAC050 is a cold-inducible transcription factor gene

Our experimental results showed that the expression of the OsNAC050 gene in leaves was significantly higher than in other tissues both in seedlings and heading stage. At heading stage, OsNAC050 showed significant expression in leaves; the relative expression in leaves was 143.95-fold that of roots, 80-fold that of stems, and 77.94-fold that of sheaths (Figure 1A). The subcellular localization results showed that OsNAC050 was mainly expressed in the nucleus (Figure 1B). In response to low temperature (4°C) treatment, the expression of OsNAC050 increased to 3.3-fold the initial expression level after about 9 h, and the expression of OsNAC050 changed significantly under low temperature stress. From the analyses of OsNAC050 expression in the different treatments, we found that OsNAC050 expression was induced by various stimuli, and it was particularly sensitive to PEG, NaCl, low temperature, and GA₃. Therefore, we hypothesized that OsNAC050 is involved in the response pathways for drought, high salt, low temperature, and GA₃ stresses in rice (Figure 1C).

3.2. Targeted mutagenesis of OsNAC050

The Cas9 sgRNA was designed to target a sequence near the 5’ end of the first exon of OsNAC050. The design site and sequence are shown in Figure 2A. Following transformation of ‘Nipponbare’, a total of 15 T₀-generation transgene-positive individual plants were obtained, and amplify the sequence near the mutation site for sequencing (Figure 2B). The results of sequencing and comparison showed that four of the 15 transgene-positive plants assayed carried mutations in OsNAC050; the mutation rate was 26.7% and the biallelic mutation rate was 6.7%. The sequence analysis results for the 15 transgenic positive plants tested are shown in Figure 2C. Among them, only OsNAC050-sgRNA01-09 was a homozygous biallelic mutant, with a single-base T insertion at the same position in both alleles, named osnac050 (T/T). OsNAC050-sgRNA01-01 has a 1bp insertion (A). OsNAC050-sgRNA01-04 has a 1-bp deletion (G) osnac050 and OsNAC050-sgRNA01-11 has a 5-bp deletion (ATGAG) in one allele. After T2 generation separation, screening homozygous mutants named osnac050 (-1/-1) and osnac050 (-5/-5). The mutants were screened and identified, and the non-transgenic and homozygous mutant plants were identified and propagated for subsequent experiments. All the following experiments were performed with T₂-generation homozygous mutant plants.

Upstream (5’) and downstream primers (3’) were designed to amplify the 10 predicted most likely off-target sites, and the target fragments were amplified for Sanger sequencing analysis (Supplemental Figure 1). Among the 10 most likely off-target sites, no off-target mutations were identified.

3.3. Loss of function of OsNAC050 increases cold tolerance in rice seedlings

After five days of low temperature treatment, the mutants were showed that grow better than the wild type (Figure 3A). The experiment was repeated several times and the results were consistent, which indicated that knocking out OsNAC050 expression resulted in improved tolerance to low temperature in rice seedlings. The results of NBT and DAB staining for the WT and osnac050 mutant plants in the treatments showed that more superoxide ion (O^2-) accumulated in the leaves of the WT lines compared to the mutants in the low temperature treatment, and the amount of blue precipitate was significantly higher in WT than in the mutant lines (Figure 3B). The relative water contents, chlorophyll contents, SOD (superoxide dismutase), POD (peroxidase), and CAT (catalase) activities were higher in the mutants, which also had less accumulation of MDA, compared to WT (Figure 3C-H). The O^2- content in the WT plants under low temperature stress was 4-4.5 times that of the mutant plants, and the accumulation of H₂O₂ in the WT plants was 3.5-4 times that of the mutants, indicating that the mutants experienced less low-temperature-induced oxidative damage (Figure 3I-J). In conclusion, the rice OsNAC050 mutant lines were better at dealing with reactive oxygen species (ROS) caused by low temperature stress, which is consistent with the increased tolerance to low temperature in the mutant plants.

3.4. Knocking out OsNAC050 expression changed global gene expression in rice.

According to the genome-wide expression analysis of osnac050 mutants, the gene expression of the whole genome was greatly changed after knocking out the OsNAC050 gene (Figure 4 A, C). The number of differentially expressed genes (DEGs) in the samples was counted. Under normal conditions, a comparison of WT and osnac050(T/T) identified 6,605 DEGs, of which 4,294 were up-regulated and 2,311 were down-regulated. The WT vs osnac050(-5/-5) comparison had a total of 2,714 DEGs, in which 1,772 showed up-regulated expression and 942 were down-regulated (Figure 4B). In the low temperature treatment, there were 3,851 DEGs between WT and osnac050(T/T), including 1,425 up-regulated genes and 2,426 down-regulated genes. The comparison of WT and osnac050(-5/-5) under low temperature had 7,000 differentially expressed genes, of which 4,200 were up-regulated and 2,800 were down-regulated (Figure 4D). Ten genes were randomly selected to verify the reliability of our transcriptome data (Supplemental figure 2).

Under the same treatment conditions, the DEGs in the intersection of the WT vs osnac050(T/T) vs osnac050(-5/-5) comparison are more representative of the differential gene expression changes resulting from knocking out OsNAC050 (Supplemental figure 3). GO (Gene Ontology) annotation was performed on some of the DEGs. A total of 73 sub-categories were annotated, among which the largest number of DEGs was concentrated in the Molecular Function (MF) category. This provides some evidence that knocking out OsNAC050 expression may affect the association between cells. The terms ‘organic cyclic compound binding’ and ‘heterocyclic compound binding’ were greatly affected, indicating that OsNAC050 may affect the synthesis of organic compounds. The significant changes in the Cellular Component (CC) category are mainly in ‘membrane structure’ (membranes and membrane parts), so we can speculate that knocking out OsNAC050 will lead to changes in the main components of the cell membrane and organelle membranes. In the Biological Process (BP) GO category, the proportion of DEGs in ‘response to stimulus’ was high, which indicates that OsNAC050 could affect the stress response in plants (Figure 4E).

The results of the DEG annotation based on the KEGG (Kyoto Encyclopedia of Genes and Genomes) database are as follows: in the no treatment group, cell cycle, DNA replication and phenylpropanoid biosynthesis were the three most affected signaling pathways. The DEGs are also concentrated in the metabolic processes of ‘mismatch repair’, ‘purine metabolism’, and ‘base excision repair’ (Supplemental Figure 4). In the low temperature group, the DEGs were mainly concentrated in the ‘metabolic pathways’ and the ‘biosynthesis of secondary metabolites. In addition, there are DEGs enriched in ‘drug metabolism cytochrome P450’ and ‘photosynthesis’, both of which can directly affect light reception rice and thus affect the photosynthesis rate. At the same time, the two pathways of ‘starch and sucrose metabolism’ and ‘amino sugar and nucleotide sugar metabolism’ were partially enriched (Supplemental Figure 5). These results may suggest that osnac050 mutant plants can withstand low temperatures by increasing the biological processes involved in soluble sugar content and photosynthetic efficiency.

Analysis of the DEGs identified possible transcription-level responses to low temperature

Because the DEGs in the KEGG analysis were annotated to the ‘photosynthesis’ and ‘starch and sucrose anabolic pathways’, and both play an important role in plant growth and development and the low temperature response. The significantly differential gene expression in the two signaling pathways was analyzed (Figure 5). There were 13 main genes in which the expression levels changed in ‘photosynthesis’. LOC4348378 participates in the photosynthetic electron transport process (GO: 0045156), transfers electrons in the cyclic electron transport pathway of photosynthesis, affects the formation of the chloroplast thylakoid membrane (GO: 0009535), and directly affects photosystem II (GO: 0009523). LOC4342404 is the main gene encoding the psaK subunit in the photosystem I reaction center in the chloroplast, and mainly affects photosystem I. LOC4333650, LOC4342395, LOC4333624, and LOC112938704 participate in the photosynthetic electron transport chain (GO: 0009767), affecting the photosystem II reaction center (GO: 0009539). LOC9270637 is involved in chlorophyll synthesis (GO:0009507) (Figure 5A). Taken together, these results suggest that OsNAC050 may regulate the response of genes that encode proteins in the photosystem I and photosystem II reaction centers to low temperature.

The anabolic pathway of starch and sucrose is a signaling mechanism that directly affects the response of plants to low temperature. The higher the soluble sugar content, the stronger the cold resistance of the plant. The expression of 26 genes in ‘starch and sucrose biological processes’ changed significantly in OsNAC05-knockout plants. Among them, the expression of LOC4344497 was up-regulated, so it may promote the fructokinase activity encoded by this gene (GO: 0008865). LOC4333062 regulates synthetic transferase activity (GO: 0016740), sucrose synthase activity (GO: 0016157), and transglycosylation activity (GO: 0016757), which directly affects sucrose metabolism (GO: 0005985); this gene is up-regulated at low temperatures and its expression promotes the above biological processes. LOC4334029 is involved in many pathways, and the protein it encodes regulates both cellulase activity (GO: 0008810) and hydrolase activity (GO: 0004553), which directly affects the process of polysaccharide catabolism (GO: 0000272) and the cellulose catabolism pathway (GO: 0030245). LOC4333841 and LOC4336142 are involved in a variety of carbohydrate metabolism and synthesis pathways, affecting hydrolase activity (GO: 0004553), β-galactosidase activity (GO: 0004565), β-mannosidase activity (GO: 0004567), β-mannosidase activity, β-Glucosidase activity (GO: 0008422), endoglucanase 1,3-β-D-glucosidase activity (GO: 0042973), pectin β-glucosidase activity (GO: 0050224), and cellobiose glucosidase activity (GO: 0080079). LOC4332041 is mainly involved in the starch catabolism pathway (GO: 0005983) and sucrose catabolism process (GO: 0005987). In summary, this study shows that knocking out OsNAC050 significantly changes the sucrose and starch metabolic pathways, and it is very likely that knocking out OsNAC050 expression can promote the accumulation of soluble sugars through the sucrose and starch metabolic synthesis pathway, thereby regulating the low temperature response in rice (Figure 5B).

4. Discussion

NAC TFs have been reported to be involved in the regulation of plant responses to various stresses, including cold, drought, and salinity (Yarra and Wei 2021; Ren et al. 2021). After the genomes of Arabidopsis, rice and other plant species were sequenced, NAC genes were systematically identified and shown to encode an important family of stress-related transcription factors (Mao et al. 2017; Liu et al. 2020). More than 100 genes have been classified as NAC family genes in rice by multiple research groups. Although several OsNAC genes have been functionally identified, the functions of many have yet to be determined. Therefore, why OsNAC TFs have so many functional homologues and how the functions of many OsNAC proteins are regulated through complex networks remains unclear. In addition, chilling injury has always been an environmental factor that limits the growth of rice, and exposure to low temperatures at the seedling stage can significantly inhibit rice development, resulting in a large reduction in grain yield (Chen et al. 2022; Kang et al. 2022). In our study, we showed that osnac050 mutant rice seedlings were tolerant to low temperatures after knocking out OsNAC050 expression. OsNAC050 expression is specifically induced in seedlings in response to cold stress conditions, but showed no obvious response to ABA, IAA, and hydrogen peroxide stresses. Under normal conditions, OsNAC050 is mainly expressed in leaves, and the relative expression is up to 100-fold higher than in other tissues such as roots, stems, and panicles. In addition, we found that the initial expression level decreased between one and three hours after exposure to low temperature, and then increased rapidly between six and twelve hours. Taken together, these results show that OsNAC050 may not be a transcription factor gene that responds initially and rapidly to cold stress, but rather may cause secondary and long-term responses that affect cold stress tolerance in rice (Figure 1).

After transforming ‘Nipponbare’ with the CRISPR/Cas9 vector designed to target a region of the first exon of OsNAC050, we identified only three mutants in the 15 transgenic plants. T₂ generation plants were isolated for the three mutations and the homozygous mutant plants were selected for propagation and then used in the experiment. Our results showed that the seeds germination slower after knocking out OsNAC050 expression. Because we found that gibberellic acid (GA₃) can inhibit expression of OsNAC050, we speculated that it might affect the gibberellin pathway, which in turn affects germination, making germination of osnac050 mutant seeds slower than in WT (Figure 2D).

Photosynthesis is probably the most basic and complex physiological process that occurs in plants, but it can be severely affected by adverse environmental conditions such as drought, salinity, and high temperature(Joo et al. 2022; Xiong et al. 2021). Photosynthesis involves sensitive pigments and photosystems as well as CO₂ reduction pathways, which, along with chloroplast ultrastructure and electron transport systems, may be disrupted by abiotic stresses, which ultimately affect photosynthetic capacity(Gong et al. 2014). The accumulation of photosynthetic pigments such as chlorophyll, which is involved in capturing solar radiation to drive the photosynthetic machinery, is a potential biochemical indicator of drought tolerance in different crops(Liu, Quan and Bartels 2022). In this study, the osnac050 mutant plants exhibited higher chlorophyll contents (Figure 3D), which is very likely the reason for the increased tolerance to low temperatures (Bao et al. 2016). ROS, such as superoxide radicals (O2-), hydrogen peroxide (H₂O₂), monooxyl radicals (O21), alkoxyl radicals (RO), and hydroxyl radicals (OH-), which are highly reactive towards proteins, lipids, carbohydrates, and DNA, are produced in plants in response to various abiotic stresses including low temperature (Ramanjulu and Bartels 2002). MDA levels are a well-known indicator used to determine the degree of lipid peroxidation caused by ROS (reactive oxygen species) damage. In the present study, a higher MDA content was detected in WT plants after exposure to low temperature stress, indicating that these plants suffered more extensive oxidative damage to cell membranes. In response to ROS toxicity, plants have evolved efficient antioxidant mechanisms to adapt(Mohamadi Esboei et al. 2022). These mechanisms mainly include non-enzymatic scavengers as well as antioxidant enzymes (Jiang et al. 2021). Non-enzymatic antioxidants include ascorbic acid and reduced glutathione, and enzymatic antioxidants include superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD), which are important in the drought response. Rapid accumulation under stress conditions can minimize oxidative damage(Toczewska et al. 2022). In the present study, mutant plants exhibited higher SOD, CAT, and POD activities than WT plants after cold stress (Figure 3). OsNAC050 enhances low temperature tolerance by stimulating the antioxidant defense system as measured by MDA content and SOD, CAT, and POD activities.

To investigate the molecular changes in the osnac50 mutant plants in more depth, we sequenced their transcriptomes and compared them to that of the wild-type ‘Nipponbare’ (Figure 4). Identification of differentially expressed genes between the mutant and WT transcriptomes showed that photosynthesis and the synthesis of starch and sucrose were the two metabolic pathways that were significantly altered in the mutants, and these are the main driving force for plant growth and can provide the necessary energy for plants to synthesize organic compounds(Di et al. 2022). The results of this study found that knocking out OsNAC050 mainly affected the synthesis of chloroplast thylakoid membranes and the transport of photosynthetic electron chains in photosystem II (Wang et al. 2022). The chloroplast is the main reaction center of photosynthesis. Chloroplasts have an inner and outer membrane, which provides a place for the attachment of enzymes and other proteins and electron transport(Li et al. 2021). Among the membrane structures, thylakoids are stacked membranous disks found within the chloroplast matrix and are where the light reactions occur, so the thylakoid membranes are essential for photosynthesis(Szechynska-Hebda et al. 2015; Lu et al. 2022). The photosynthetic electron transport chain is a series of systems responsible for connecting electrons on the photosynthetic membrane. This system can carry the main substrates of photosynthesis such as water or oxygen while transferring electrons(Gai et al. 2020; Yamori et al. 2011). Due to the obvious changes in these two signaling pathways, we can speculate that knocking out OsNAC050 expression may promote the light response in photosynthesis, thereby enhancing low temperature tolerance in the osnac050 mutants.

The synthesis of starch and sucrose are also important products of photosynthesis (Lopez-Gonzalez et al. 2022). Low temperature stress will lead to a reduction in enzyme activity, affecting the ability of plant cells to synthesize starch and sucrose (Apriyanto et al. 2022). At the same time, carbohydrate accumulation can help plants resist damage caused by low temperatures, and the sucrose component in the soluble sugars can also increase plant tolerance to low temperature exposure (Yu, Xu, et al. 2022; Xie et al. 2022). Our comparative transcriptomic analysis found that after knocking out OsNAC050, most of the genes encoding enzymes such as sucrose synthase, cellulase, hydrolase, and others were down-regulated, and these genes are mainly involved in various catabolic pathways, such as fructose catabolism, polysaccharide decomposition, cellulose hydrolysis, glucose decomposition, and starch catabolism (Shi et al. 2022; Nilholm et al. 2022). This implies that the osnac050 mutant plants respond to low temperature because the synthesis rate of enzymes in the catabolism processes is reduced, resulting in a decrease in the decomposition rate, thus leading to the accumulation of sucrose and starch, which in turn increases low temperature tolerance in the mutant plants (Figure 5)(Roth et al. 2022). This part of our study analyzed the phenotypic and transcriptome data of the osnac050 mutants in a deep exploration of the possible molecular mechanism of low temperature tolerance in rice, which will have important applications in rice breeding and research.

5. Conclusions

The above results show that the expression of the rice transcription factor gene OsNAC050 is induced by low temperature. After knocking out the expression of OsNAC050, the mutant plants displayed a low temperature-tolerant phenotype when exposed to low temperatures. A comparative analysis of the transcriptomes of ‘Nipponbare’ and the osnac50 mutants by RNA-seq was used to determine the possible mechanisms, and it was found that the ‘photosynthesis’ and ‘sucrose and glucose’ synthesis pathways were affected after knocking out the expression of OsNAC050. The differentially expressed genes in ‘photosynthesis’ included several photosynthetic pathways involving genes of the photosystem I and photosystem II active centers. In the sucrose and glucose anabolic pathways, we focused on genes involved in the decomposition of cellulose and sucrose, and we found that most of the metabolic pathways involved in the processes of decomposition and hydrolysis were changed. In conclusion, we analyzed the expression characteristics and gene function of the transcription factor gene OsNAC050, and determined its possible mechanism of action and regulatory target sites. The results of our study provide an important theoretical basis for breeding for low temperature tolerance in rice.