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 28oC 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), H2O2 (1%), IAA (100 μM), ABA (100 μM), and GA3 (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 T2 -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).
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
2 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
3) 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
2 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
2O
2), 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.
Figure 1.
OsNAC050 expression profile analysis. (A) The expression of OsNACO50 in different tissues and organs of rice was detected and quantified using qRT-PCR. Samples of roots, sheaths, and leaves were collected from 4-week old seedlings, and roots, stems, sheaths, leaves, and ears were sampled during the reproductive growth period from plants prior to heading. The error bars represents the standard error (SE) based on three biological replicates. (B) Nuclear localization of the OsNAC050 protein in rice protoplasts. NLS, nuclear location signal. Scale bars=10 µm. (C) OsNAC050 expression levels in response to various abiotic stresses and hormone treatments. Four week old seedlings were treated with cold treatment (4 oC), heat treatment (42oC), PEG 6000 (20%, w/v), NaCl (200 mm), H2O2 (1%), IAA (100 µm), ABA (100 μm) And GA3 (100 µm). The relative expression levels of OsNAC050 were determined using qRT-PCR at the specified times. The error bars indicate that the SE is based on three independent biological replicates.
Figure 1.
OsNAC050 expression profile analysis. (A) The expression of OsNACO50 in different tissues and organs of rice was detected and quantified using qRT-PCR. Samples of roots, sheaths, and leaves were collected from 4-week old seedlings, and roots, stems, sheaths, leaves, and ears were sampled during the reproductive growth period from plants prior to heading. The error bars represents the standard error (SE) based on three biological replicates. (B) Nuclear localization of the OsNAC050 protein in rice protoplasts. NLS, nuclear location signal. Scale bars=10 µm. (C) OsNAC050 expression levels in response to various abiotic stresses and hormone treatments. Four week old seedlings were treated with cold treatment (4 oC), heat treatment (42oC), PEG 6000 (20%, w/v), NaCl (200 mm), H2O2 (1%), IAA (100 µm), ABA (100 μm) And GA3 (100 µm). The relative expression levels of OsNAC050 were determined using qRT-PCR at the specified times. The error bars indicate that the SE is based on three independent biological replicates.
Figure 2.
Targeted mutagenesis using CRISPR-Cas9 gene editing. (A) sgRNA site design in the first exon of OsNAC050. The target sequence is shown in black, and the PAM (protospacer adjacent motif) site is shown in red. (B) Gel electrophoresis of amplified DNA fragments from 15 independent OsNAC050-sgRNA01 T0 -generation lines for sanger sequencing to identify the mutant genotype. WT, wild type ‘Nipponbare’. (C) Sanger sequencing of the targets site of OsNAC050-sgRNA01 T0 -generation lines. (D) Phenotypic analysis of OsNAC050 mutant lines under germinating period. Experiments showthatthegerminationspeedofmutantsdecreases.
Figure 2.
Targeted mutagenesis using CRISPR-Cas9 gene editing. (A) sgRNA site design in the first exon of OsNAC050. The target sequence is shown in black, and the PAM (protospacer adjacent motif) site is shown in red. (B) Gel electrophoresis of amplified DNA fragments from 15 independent OsNAC050-sgRNA01 T0 -generation lines for sanger sequencing to identify the mutant genotype. WT, wild type ‘Nipponbare’. (C) Sanger sequencing of the targets site of OsNAC050-sgRNA01 T0 -generation lines. (D) Phenotypic analysis of OsNAC050 mutant lines under germinating period. Experiments showthatthegerminationspeedofmutantsdecreases.
Figure 3.
Cold-tolerance in the osnac050 mutants. (A) Phenotypes of osnac050 T2 -generation mutant plants exposed to cold stress. (B) Levels of O2- and H2O2 in WT and osnac050 T2 -generation mutant lines subjected to cold stress. Rice leaf samples were stained separately with nitroblue tetrazolium (NBT) and diaminobenzidine (DAB). (C) Determination of relative water content and (D) chlorophyll content after 5-days of low temperature stress. (E-H) Malondialdehyde (MDA) content (E), Superoxide dismutase (SOD) activity (F), Peroxidase (POD) activity (G), and Catalase (CAT) activity after low temperature stress (H). O2- production rate and (J) H2O2 content after low temperature stress. Bars represent the mean ± SE of three independent experiments. ∗ and ∗ ∗ represent significant differences at p.
Figure 3.
Cold-tolerance in the osnac050 mutants. (A) Phenotypes of osnac050 T2 -generation mutant plants exposed to cold stress. (B) Levels of O2- and H2O2 in WT and osnac050 T2 -generation mutant lines subjected to cold stress. Rice leaf samples were stained separately with nitroblue tetrazolium (NBT) and diaminobenzidine (DAB). (C) Determination of relative water content and (D) chlorophyll content after 5-days of low temperature stress. (E-H) Malondialdehyde (MDA) content (E), Superoxide dismutase (SOD) activity (F), Peroxidase (POD) activity (G), and Catalase (CAT) activity after low temperature stress (H). O2- production rate and (J) H2O2 content after low temperature stress. Bars represent the mean ± SE of three independent experiments. ∗ and ∗ ∗ represent significant differences at p.
Figure 4.
Global gene expression changes in OsNAC050-knockout rice plants. (A) Cluster analysis of the most important differentially expressed genes (DEGs) between the WT and osnac050 T2 -generation mutant lines. Targeted knockout of OsNAC050 expression resulted in changes in gene expression, physiology, growth, and development compared with ‘Nipponbare’ and the control without low temperature stress treatment. Relative expression corresponds to the log2 (FPKM) value of each gene. WT, osnac050(T/T), and osnac050(-5/-5) T2 mutant line, expression profiles were obtained from RNA-seq data. Total RNA was extracted from a mixture of three different plants. Genes with p-value < 0.05 and foldchange >2 are advised to define differentially expressed genes (DEGs). (B) DEGs in the WT vs. osnac050(T/T) and WT vs. osnac050(-5/-5) comparisons without cold stress. (C) Cluster analysis of the most important DEGs between the WT and osnac050 T2 mutant lines after cold stress. (D) DEGs in the WT vs. osnac050(T/T) and WT vs. osnac050(-5/-5) lines after cold stress. (E) Gene ontology (GO) classification of DEGs between the WT and two osnac050 mutants under normal (red) and cold stress conditions (blue) into the three main GO classes CC (Cellular Component), MF (Molecular Function), and BP (Biological Process). The x-axis displays the GO terms selected by the user, and the y-axis displays the percentage of genes (specific genes divided by total genes.
Figure 4.
Global gene expression changes in OsNAC050-knockout rice plants. (A) Cluster analysis of the most important differentially expressed genes (DEGs) between the WT and osnac050 T2 -generation mutant lines. Targeted knockout of OsNAC050 expression resulted in changes in gene expression, physiology, growth, and development compared with ‘Nipponbare’ and the control without low temperature stress treatment. Relative expression corresponds to the log2 (FPKM) value of each gene. WT, osnac050(T/T), and osnac050(-5/-5) T2 mutant line, expression profiles were obtained from RNA-seq data. Total RNA was extracted from a mixture of three different plants. Genes with p-value < 0.05 and foldchange >2 are advised to define differentially expressed genes (DEGs). (B) DEGs in the WT vs. osnac050(T/T) and WT vs. osnac050(-5/-5) comparisons without cold stress. (C) Cluster analysis of the most important DEGs between the WT and osnac050 T2 mutant lines after cold stress. (D) DEGs in the WT vs. osnac050(T/T) and WT vs. osnac050(-5/-5) lines after cold stress. (E) Gene ontology (GO) classification of DEGs between the WT and two osnac050 mutants under normal (red) and cold stress conditions (blue) into the three main GO classes CC (Cellular Component), MF (Molecular Function), and BP (Biological Process). The x-axis displays the GO terms selected by the user, and the y-axis displays the percentage of genes (specific genes divided by total genes.
Figure 5.
Transcriptomic analysis of system regulatory genes in the WT and osnac050 mutant lines under low temperature stress. (A and B) Heatmaps showing changes in gene expression in the DEGs in two major pathways between WT and two osnac050 mutants with and without cold treatment. (A) Photosynthesis. (B) Starch and sucrose metabolism. Log2 fold-change (FC) values for DEGs in WT and osnac050(T/T) and osnac050(-5/-5) mutant lines are shown before (cold-) and after (cold+) low temperature treatment. Genes with p-value < 0.05 and fold change >2 are advised to define differentially expressed genes.
Figure 5.
Transcriptomic analysis of system regulatory genes in the WT and osnac050 mutant lines under low temperature stress. (A and B) Heatmaps showing changes in gene expression in the DEGs in two major pathways between WT and two osnac050 mutants with and without cold treatment. (A) Photosynthesis. (B) Starch and sucrose metabolism. Log2 fold-change (FC) values for DEGs in WT and osnac050(T/T) and osnac050(-5/-5) mutant lines are shown before (cold-) and after (cold+) low temperature treatment. Genes with p-value < 0.05 and fold change >2 are advised to define differentially expressed genes.