-
PDF
- Split View
-
Views
-
Cite
Cite
Thanikarn Udomchalothorn, Kitiporn Plaimas, Siriporn Sripinyowanich, Chutamas Boonchai, Thammaporn Kojonna, Panita Chutimanukul, Luca Comai, Teerapong Buaboocha, Supachitra Chadchawan, OsNucleolin1-L Expression in Arabidopsis Enhances Photosynthesis via Transcriptome Modification under Salt Stress Conditions, Plant and Cell Physiology, Volume 58, Issue 4, April 2017, Pages 717–734, https://doi.org/10.1093/pcp/pcx024
- Share Icon Share
OsNUC1 encodes rice nucleolin, which has been shown to be involved in salt stress responses. Expression of the full-length OsNUC1 gene in Arabidopsis resulted in hypersensitivity to ABA during germination. Transcriptome analysis of the transgenic lines, in comparison with the wild type, revealed that the RNA abundance of >1,900 genes was significantly changed under normal growth conditions, while under salt stress conditions the RNAs of 999 genes were found to be significantly regulated. Gene enrichment analysis showed that under normal conditions OsNUC1 resulted in repression of genes involved in photosynthesis, while in salt stress conditions OsNUC1 increased expression of the genes involved in the light-harvesting complex. Correspondingly, the net rate of photosynthesis of the transgenic lines was increased under salt stress. Transgenic rice lines with overexpression of the OsNUC1-L gene were generated and tested for photosynthetic performance under salt stress conditions. The transgenic rice lines treated with salt stress at the booting stage had a higher photosynthetic rate and stomatal conductance in flag leaves and second leaves than the wild type. Moreover, higher contents of Chl a and carotenoids were found in flag leaves of the transgenic rice. These results suggest a role for OsNUC1 in the modification of the transcriptome, especially the gene transcripts responsible for photosynthesis, leading to stabilization of photosynthesis under salt stress conditions.
Introduction
Environmental stress such as drought and saline soils can reduce various activities essential for respiration (Criddle et al. 1989) and photosynthesis (Yeo and Flower 1982), leading to reduction of crop yield. Exposure of plants to these stresses induces expression changes in genes with widely different roles (Yamaguchi-Shinozaki and Shinozaki 2006), presumably in order to cope with the deleterious effects of stress. The plant hormone ABA mediates many plant responses to physical stress, especially the response to water status such as drought and salinity. ABA regulates stomatal closure, plant development, such as seed maturation, germination, seedling growth and also adaptive responses to osmotic stress (Zhang et al. 2006).
RNA-binding proteins play diverse roles in plant growth and development (Lu and Fedoroff 2000). Several reports described the involvement of RNA-binding proteins in the ABA response and also in plant response to extreme environments such as those involving regulation of HYL1, ABH1 and CBP20 (Lu and Fedoroff 2000, Hugouvieux et al. 2001, Papp et al. 2004). Nucleolin (NUC) is an RNA-binding protein predominant in the nucleolus particularly around fibrillar centers (FCs) and the dense fibrillar component (DFC) (Tong et al. 1997, Pontvianne et al. 2007, Stepinski 2012). NUC is composed of four domains: the nuclear localization signal (NLS); the N-terminal acidic region; two RNA-recognition motifs (RRMs); and a glycine–arginine rich (GAR) region. Arabidopsis Nucleolin1, AtNUC-L1, the predominant nucleolin in Arabidopsis, is reported to be involved in sugar-inducible expression of ribosomal proteins. The disruption of the AtNUC-L1 gene resulted in phenotypes similar to those of ribosomal protein gene mutants (Kojima et al. 2007).
Two NUC genes are found in the rice genome, OsNUC1 and OsNUC2. Only OsNUC1 responds to salt stress conditions. Rice lines with higher salt tolerance display a higher level of OsNUC1 gene expression. In a previous study, we discovered by Northern blot analysis that OsNUC1 had two mRNA forms: OsNUC1-S and OsNUC1-L, 1.9 and 2.4 kb respectively. Overexpression of the OsNUC1-S gene in rice and its expression in Arabidopsis conferred salt tolerance in both species. Enhancement of the AtSOS1 and AtP5CS1 transcripts in Arabidopsis expressing OsNUC1-S suggests that the RRM and GAR regions of OsNUC1 function in salt resistance (Sripinyowanich et al. 2013).
Here, we describe the effect of expressing the OsNUC1-L isoform in Arabidopsis thaliana. We used transcriptome analysis to determine the impact of OsNUC1-L on global mRNAs, connecting the molecular changes to the salt tolerance phenotype of the transgenic lines. The overexpression of the OsNUC1-L isoform in rice was used to demonstrate that the phenotype in Arabidopsis is not due to artificial gain of function in the heterologous system.
Results
OsNUC1-L improved salt tolerance in Arabidopsis

The relative level of OsNUC1-L transcripts in transgenic Arabidopsis with the 35SCaMV:OsNUC1-L construct (T14, T22 and T24) (A). The phenotype of 16-day-old transgenic Arabidopsis expressing OsNUC1-L and the wild type, when grown on MS medium or MS medium supplemented with 150 mM NaCl (B). The relative growth rate of wild-type (WT) and transgenic Arabidopsis on MS medium supplemented with 150 mM NaCl (C) was determined using a completely randomized design (CRD) with 10 replicates, each of which consisted of at least 10 seedlings. Means with a different letter above the bar are significantly different (P < 0.05, Duncan’s Multiple Range Test).
All three transgenic lines expressing OsNUC1-L displayed increased salt tolerance, maintaining Chl content after growth on MS medium supplemented with 150 mM NaCl for 3 weeks (Fig. 1B). Interestingly, expression of OsNUC1-L did not prevent dry weight loss in Arabidopsis under salt stress conditions. The expression of OsNUC1-L only slightly enhanced relative growth rate during salt stress (Fig. 1C). This is in contrast to transgenic Arabidopsis expressing OsNUC1-S, which showed enhanced plant growth under stress conditions (Sripinyowanich et al. 2013).
OsNUC1-L altered the response to ABA

Germination response of transgenic OsNUC1 (T14, T22 and T24, dashed lines) and wild-type Arabidopsis (solid line) on MS medium (A, B) and MS medium supplemented with various concentrations of ABA (C, D).
In order to understand the enhanced response to ABA, the expression of genes encoding ABA synthesis and ABA-responsive genes was investigated via transcriptome analysis of the T22 line. We found that OsNUC1-L expression reduced the expression of the NCED3 gene, encoding 9-cis-epoxycarotenoid dioxygenase, which catalyzes the rate-limiting step in ABA biosynthesis, conversion of 9-cis-violaxanthin to xanthoxin (Qin and Zeevaart 2002) (Supplementary Fig. S2). At the same time, expression of OsNUC1-L induced the expression of AAO3 and NCED6 (Table 1). Moreover, expression of the gene encoding CYP707A4, an enzyme mediating ABA degradation, was 5-fold higher than in the wild type. The expression of some ABA-responsive genes was also changed in the transgenic lines. We found a significant reduction of ABA-responsive gene expression, such as genes encoding the late embryogenesis abundant protein RD29A, and dehydrin RAB18-like protein (Table 2).
mRNA abundance of ABA biosynthesis genes in 35SCaMV:OsNUC1 and wild-type Arabidopsis grown on MS medium with and without 150 mM NaCl
Gene (locus No.) . | Control . | Salt stress (150 mM NaCl) . | ||||
---|---|---|---|---|---|---|
. | Average read (after normalization) . | Fold change . | Average read (after normalization) . | Fold change . | ||
. | WT . | 35SCaMV:OsNUC1-L . | . | WT . | 35SCaMV:OsNUC1-L . | . |
ABA1/ZEP (AT5G67030) | 74.96 | 82.60 | 1.10 | 121.67 | 106.96 | 0.88 |
ABA4 (AT1G67080) | 0.64 | 0.60 | 0.93 | 1.64 | 1.01 | 0.61 |
NCED2 (AT4G18350) | 197.04 | 195.78 | 0.99 | 203.54 | 209.30 | 1.03 |
NCED3 (AT3G14440) | 9,811.30 | 4,342.56 | 0.44* | 6,010.93 | 3,845.76 | 0.64* |
NCED5 (AT1G30100) | 207.46 | 225.80 | 1.09 | 235.47 | 281.34 | 1.19 |
NCED6 (AT3G24220) | 0.32 | 1.58 | 4.91 | 2.68 | 1.65 | 0.62 |
NCED9 (AT1G78390) | 159.55 | 177.61 | 1.11 | 224.41 | 213.82 | 0.95 |
ABA2 (AT1G52340) | 21.90 | 30.74 | 1.40 | 56.09 | 44.96 | 0.80 |
AAO3 (AT2G27150) | 67.89 | 183.65 | 2.70 | 340.96 | 160.49 | 0.47* |
ABA3 (AT1G16540) | 20.59 | 22.82 | 1.11 | 33.29 | 14.06 | 0.42* |
CYP707A1 (AT4G19230) | – | – | – | – | – | – |
CYP707A2 (AT2G29090) | 2.39 | 2.91 | 1.22 | 8.56 | 3.66 | 0.43 |
CYP707A3 (AT5G45340) | 63.44 | 74.46 | 1.17 | 317.68 | 256.78 | 0.81 |
CYP707A4 (AT3G19270) | 21.92 | 122.50 | 5.59* | 76.42 | 75.99 | 0.99 |
BG1 (AT1G52400) | 501.21 | 585.94 | 1.17 | 519.16 | 546.56 | 1.05 |
Gene (locus No.) . | Control . | Salt stress (150 mM NaCl) . | ||||
---|---|---|---|---|---|---|
. | Average read (after normalization) . | Fold change . | Average read (after normalization) . | Fold change . | ||
. | WT . | 35SCaMV:OsNUC1-L . | . | WT . | 35SCaMV:OsNUC1-L . | . |
ABA1/ZEP (AT5G67030) | 74.96 | 82.60 | 1.10 | 121.67 | 106.96 | 0.88 |
ABA4 (AT1G67080) | 0.64 | 0.60 | 0.93 | 1.64 | 1.01 | 0.61 |
NCED2 (AT4G18350) | 197.04 | 195.78 | 0.99 | 203.54 | 209.30 | 1.03 |
NCED3 (AT3G14440) | 9,811.30 | 4,342.56 | 0.44* | 6,010.93 | 3,845.76 | 0.64* |
NCED5 (AT1G30100) | 207.46 | 225.80 | 1.09 | 235.47 | 281.34 | 1.19 |
NCED6 (AT3G24220) | 0.32 | 1.58 | 4.91 | 2.68 | 1.65 | 0.62 |
NCED9 (AT1G78390) | 159.55 | 177.61 | 1.11 | 224.41 | 213.82 | 0.95 |
ABA2 (AT1G52340) | 21.90 | 30.74 | 1.40 | 56.09 | 44.96 | 0.80 |
AAO3 (AT2G27150) | 67.89 | 183.65 | 2.70 | 340.96 | 160.49 | 0.47* |
ABA3 (AT1G16540) | 20.59 | 22.82 | 1.11 | 33.29 | 14.06 | 0.42* |
CYP707A1 (AT4G19230) | – | – | – | – | – | – |
CYP707A2 (AT2G29090) | 2.39 | 2.91 | 1.22 | 8.56 | 3.66 | 0.43 |
CYP707A3 (AT5G45340) | 63.44 | 74.46 | 1.17 | 317.68 | 256.78 | 0.81 |
CYP707A4 (AT3G19270) | 21.92 | 122.50 | 5.59* | 76.42 | 75.99 | 0.99 |
BG1 (AT1G52400) | 501.21 | 585.94 | 1.17 | 519.16 | 546.56 | 1.05 |
Gene expression was compared using DESeq (Anders and Huber 2010) and adjusting the false discovery rate using the Benjamini–Hochberg method (FDA < 0.05). Fold changes crossing the threshold of P = 0.05 are marked in red and with an asterisk.
mRNA abundance of ABA biosynthesis genes in 35SCaMV:OsNUC1 and wild-type Arabidopsis grown on MS medium with and without 150 mM NaCl
Gene (locus No.) . | Control . | Salt stress (150 mM NaCl) . | ||||
---|---|---|---|---|---|---|
. | Average read (after normalization) . | Fold change . | Average read (after normalization) . | Fold change . | ||
. | WT . | 35SCaMV:OsNUC1-L . | . | WT . | 35SCaMV:OsNUC1-L . | . |
ABA1/ZEP (AT5G67030) | 74.96 | 82.60 | 1.10 | 121.67 | 106.96 | 0.88 |
ABA4 (AT1G67080) | 0.64 | 0.60 | 0.93 | 1.64 | 1.01 | 0.61 |
NCED2 (AT4G18350) | 197.04 | 195.78 | 0.99 | 203.54 | 209.30 | 1.03 |
NCED3 (AT3G14440) | 9,811.30 | 4,342.56 | 0.44* | 6,010.93 | 3,845.76 | 0.64* |
NCED5 (AT1G30100) | 207.46 | 225.80 | 1.09 | 235.47 | 281.34 | 1.19 |
NCED6 (AT3G24220) | 0.32 | 1.58 | 4.91 | 2.68 | 1.65 | 0.62 |
NCED9 (AT1G78390) | 159.55 | 177.61 | 1.11 | 224.41 | 213.82 | 0.95 |
ABA2 (AT1G52340) | 21.90 | 30.74 | 1.40 | 56.09 | 44.96 | 0.80 |
AAO3 (AT2G27150) | 67.89 | 183.65 | 2.70 | 340.96 | 160.49 | 0.47* |
ABA3 (AT1G16540) | 20.59 | 22.82 | 1.11 | 33.29 | 14.06 | 0.42* |
CYP707A1 (AT4G19230) | – | – | – | – | – | – |
CYP707A2 (AT2G29090) | 2.39 | 2.91 | 1.22 | 8.56 | 3.66 | 0.43 |
CYP707A3 (AT5G45340) | 63.44 | 74.46 | 1.17 | 317.68 | 256.78 | 0.81 |
CYP707A4 (AT3G19270) | 21.92 | 122.50 | 5.59* | 76.42 | 75.99 | 0.99 |
BG1 (AT1G52400) | 501.21 | 585.94 | 1.17 | 519.16 | 546.56 | 1.05 |
Gene (locus No.) . | Control . | Salt stress (150 mM NaCl) . | ||||
---|---|---|---|---|---|---|
. | Average read (after normalization) . | Fold change . | Average read (after normalization) . | Fold change . | ||
. | WT . | 35SCaMV:OsNUC1-L . | . | WT . | 35SCaMV:OsNUC1-L . | . |
ABA1/ZEP (AT5G67030) | 74.96 | 82.60 | 1.10 | 121.67 | 106.96 | 0.88 |
ABA4 (AT1G67080) | 0.64 | 0.60 | 0.93 | 1.64 | 1.01 | 0.61 |
NCED2 (AT4G18350) | 197.04 | 195.78 | 0.99 | 203.54 | 209.30 | 1.03 |
NCED3 (AT3G14440) | 9,811.30 | 4,342.56 | 0.44* | 6,010.93 | 3,845.76 | 0.64* |
NCED5 (AT1G30100) | 207.46 | 225.80 | 1.09 | 235.47 | 281.34 | 1.19 |
NCED6 (AT3G24220) | 0.32 | 1.58 | 4.91 | 2.68 | 1.65 | 0.62 |
NCED9 (AT1G78390) | 159.55 | 177.61 | 1.11 | 224.41 | 213.82 | 0.95 |
ABA2 (AT1G52340) | 21.90 | 30.74 | 1.40 | 56.09 | 44.96 | 0.80 |
AAO3 (AT2G27150) | 67.89 | 183.65 | 2.70 | 340.96 | 160.49 | 0.47* |
ABA3 (AT1G16540) | 20.59 | 22.82 | 1.11 | 33.29 | 14.06 | 0.42* |
CYP707A1 (AT4G19230) | – | – | – | – | – | – |
CYP707A2 (AT2G29090) | 2.39 | 2.91 | 1.22 | 8.56 | 3.66 | 0.43 |
CYP707A3 (AT5G45340) | 63.44 | 74.46 | 1.17 | 317.68 | 256.78 | 0.81 |
CYP707A4 (AT3G19270) | 21.92 | 122.50 | 5.59* | 76.42 | 75.99 | 0.99 |
BG1 (AT1G52400) | 501.21 | 585.94 | 1.17 | 519.16 | 546.56 | 1.05 |
Gene expression was compared using DESeq (Anders and Huber 2010) and adjusting the false discovery rate using the Benjamini–Hochberg method (FDA < 0.05). Fold changes crossing the threshold of P = 0.05 are marked in red and with an asterisk.
Expression levels of known ABA-inducible genes (Hoth et al. 2002) in 35SCaMV:OsNUC1-L transgenics compared with wild-type Arabidopsis
Gene (locus no.) . | Control . | Salt stress . | ||||
---|---|---|---|---|---|---|
. | Average reads (after normalization) . | Log2fold change . | Average reads (after normalization) . | Log2fold change . | ||
. | WT . | 35SCaMV:OsNUC1-L . | . | WT . | 35SCaMV:OsNUC1-L . | . |
Late embryogenesis abundant protein (At5g06760) | 6,703.68 | 4,298.66 | 0.64 | 12,235.58 | 7,962.26 | 0.65* |
Low-temperature-induced protein 78, RD29A (At5g52310) | 2,101.38 | 1,589.91 | 0.76 | 6,300.14 | 4,290.21 | 0.68* |
Homeobox-leucine zipper protein AtHB-12 (At3g61890) | 180.11 | 249.01 | 1.38 | 187.25 | 177.67 | 0.95 |
AtACP5, acid phosphatase type 5 (At3g17790) | 254.86 | 306.23 | 1.20 | 619.01 | 528.03 | 0.85 |
AtKIN1 (At5g15960) | 614.63 | 363.04 | 0.59* | 1,466.41 | 910.95 | 0.62* |
Dehydrin RAB18-like protein (At5g66400) | 60.79 | 124.36 | 2.05* | 462.96 | 358.55 | 0.77 |
Homeodomain transcription factor AtHB-(At2g46680) | 532.97 | 739.37 | 1.39 | 1,953.16 | 1,729.46 | 0.89 |
AtP5C1, delta-1-pyrroline 5-carboxylase synthetase (At2g39800) | 320.81 | 565.34 | 1.76* | 1,053.88 | 900.09 | 0.85 |
AtERD10, hypothetical protein (At1g20450) | 1,262.55 | 1,959.37 | 1.55* | 4,221.04 | 4,231.15 | 1.00 |
AtRD20 (At2g33380) | 11,536.06 | 9,246.12 | 0.80 | 7,888.81 | 7,661.10 | 0.97 |
AtABI1 (At4g26080) | 180.31 | 214.81 | 1.19 | 819.79 | 586.52 | 0.72* |
AtABF3 (At4g34000) | 599.62 | 633.01 | 1.06 | 989.19 | 774.77 | 0.78* |
AtABI2 (At5g57050) | 88.21 | 122.78 | 1.39 | 284.56 | 227.19 | 0.80 |
Cold-regulated protein cor15b precursor (At2g42530) | 730.31 | 1,025.37 | 1.40 | 1,488.52 | 1,276.51 | 0.86 |
Alcohol dehydrogenase (At1g77120) | 282.40 | 320.51 | 1.13 | 817.78 | 706.51 | 0.86 |
Putative receptor-like protein kinase RPK1 (At1g69270) | 15.02 | 15.20 | 1.01 | 166.09 | 106.02 | 0.64 |
Hypothetical protein, AtCOR47 (At1g20440) | 107.39 | 98.87 | 0.92 | 185.89 | 161.51 | 0.87 |
Cold-regulated protein COR6.6, KIN2 (At5g15970) | 100.40 | 143.60 | 1.43 | 3,285.65 | 3,414.42 | 1.04 |
Cold-regulated protein cor15a precursor (At2g42540) | 2,841.94 | 2,864.95 | 1.01 | 13,720.20 | 11,956.36 | 0.87 |
ABA-responsive protein-like (At5g13200) | 1,603.03 | 1,162.40 | 0.73 | 7,848.29 | 6,324.54 | 0.81 |
Protein phosphatase 2C (AtP2C-HA) (At1g72770) | 76.42 | 114.33 | 1.50 | 236.49 | 212.89 | 0.90 |
PIP1B aquaporin (At2g45960) | 126.37 | 117.83 | 0.93 | 303.26 | 241.84 | 0.80 |
Gene (locus no.) . | Control . | Salt stress . | ||||
---|---|---|---|---|---|---|
. | Average reads (after normalization) . | Log2fold change . | Average reads (after normalization) . | Log2fold change . | ||
. | WT . | 35SCaMV:OsNUC1-L . | . | WT . | 35SCaMV:OsNUC1-L . | . |
Late embryogenesis abundant protein (At5g06760) | 6,703.68 | 4,298.66 | 0.64 | 12,235.58 | 7,962.26 | 0.65* |
Low-temperature-induced protein 78, RD29A (At5g52310) | 2,101.38 | 1,589.91 | 0.76 | 6,300.14 | 4,290.21 | 0.68* |
Homeobox-leucine zipper protein AtHB-12 (At3g61890) | 180.11 | 249.01 | 1.38 | 187.25 | 177.67 | 0.95 |
AtACP5, acid phosphatase type 5 (At3g17790) | 254.86 | 306.23 | 1.20 | 619.01 | 528.03 | 0.85 |
AtKIN1 (At5g15960) | 614.63 | 363.04 | 0.59* | 1,466.41 | 910.95 | 0.62* |
Dehydrin RAB18-like protein (At5g66400) | 60.79 | 124.36 | 2.05* | 462.96 | 358.55 | 0.77 |
Homeodomain transcription factor AtHB-(At2g46680) | 532.97 | 739.37 | 1.39 | 1,953.16 | 1,729.46 | 0.89 |
AtP5C1, delta-1-pyrroline 5-carboxylase synthetase (At2g39800) | 320.81 | 565.34 | 1.76* | 1,053.88 | 900.09 | 0.85 |
AtERD10, hypothetical protein (At1g20450) | 1,262.55 | 1,959.37 | 1.55* | 4,221.04 | 4,231.15 | 1.00 |
AtRD20 (At2g33380) | 11,536.06 | 9,246.12 | 0.80 | 7,888.81 | 7,661.10 | 0.97 |
AtABI1 (At4g26080) | 180.31 | 214.81 | 1.19 | 819.79 | 586.52 | 0.72* |
AtABF3 (At4g34000) | 599.62 | 633.01 | 1.06 | 989.19 | 774.77 | 0.78* |
AtABI2 (At5g57050) | 88.21 | 122.78 | 1.39 | 284.56 | 227.19 | 0.80 |
Cold-regulated protein cor15b precursor (At2g42530) | 730.31 | 1,025.37 | 1.40 | 1,488.52 | 1,276.51 | 0.86 |
Alcohol dehydrogenase (At1g77120) | 282.40 | 320.51 | 1.13 | 817.78 | 706.51 | 0.86 |
Putative receptor-like protein kinase RPK1 (At1g69270) | 15.02 | 15.20 | 1.01 | 166.09 | 106.02 | 0.64 |
Hypothetical protein, AtCOR47 (At1g20440) | 107.39 | 98.87 | 0.92 | 185.89 | 161.51 | 0.87 |
Cold-regulated protein COR6.6, KIN2 (At5g15970) | 100.40 | 143.60 | 1.43 | 3,285.65 | 3,414.42 | 1.04 |
Cold-regulated protein cor15a precursor (At2g42540) | 2,841.94 | 2,864.95 | 1.01 | 13,720.20 | 11,956.36 | 0.87 |
ABA-responsive protein-like (At5g13200) | 1,603.03 | 1,162.40 | 0.73 | 7,848.29 | 6,324.54 | 0.81 |
Protein phosphatase 2C (AtP2C-HA) (At1g72770) | 76.42 | 114.33 | 1.50 | 236.49 | 212.89 | 0.90 |
PIP1B aquaporin (At2g45960) | 126.37 | 117.83 | 0.93 | 303.26 | 241.84 | 0.80 |
Gene expression was compared using DESeq (Anders and Huber 2010) and adjusting the false discovery rate using the Benjamini–Hochberg method (FDA < 0.05). Fold changes crossing the threshold of P = 0.05 are marked in red and with an asterisk.
Expression levels of known ABA-inducible genes (Hoth et al. 2002) in 35SCaMV:OsNUC1-L transgenics compared with wild-type Arabidopsis
Gene (locus no.) . | Control . | Salt stress . | ||||
---|---|---|---|---|---|---|
. | Average reads (after normalization) . | Log2fold change . | Average reads (after normalization) . | Log2fold change . | ||
. | WT . | 35SCaMV:OsNUC1-L . | . | WT . | 35SCaMV:OsNUC1-L . | . |
Late embryogenesis abundant protein (At5g06760) | 6,703.68 | 4,298.66 | 0.64 | 12,235.58 | 7,962.26 | 0.65* |
Low-temperature-induced protein 78, RD29A (At5g52310) | 2,101.38 | 1,589.91 | 0.76 | 6,300.14 | 4,290.21 | 0.68* |
Homeobox-leucine zipper protein AtHB-12 (At3g61890) | 180.11 | 249.01 | 1.38 | 187.25 | 177.67 | 0.95 |
AtACP5, acid phosphatase type 5 (At3g17790) | 254.86 | 306.23 | 1.20 | 619.01 | 528.03 | 0.85 |
AtKIN1 (At5g15960) | 614.63 | 363.04 | 0.59* | 1,466.41 | 910.95 | 0.62* |
Dehydrin RAB18-like protein (At5g66400) | 60.79 | 124.36 | 2.05* | 462.96 | 358.55 | 0.77 |
Homeodomain transcription factor AtHB-(At2g46680) | 532.97 | 739.37 | 1.39 | 1,953.16 | 1,729.46 | 0.89 |
AtP5C1, delta-1-pyrroline 5-carboxylase synthetase (At2g39800) | 320.81 | 565.34 | 1.76* | 1,053.88 | 900.09 | 0.85 |
AtERD10, hypothetical protein (At1g20450) | 1,262.55 | 1,959.37 | 1.55* | 4,221.04 | 4,231.15 | 1.00 |
AtRD20 (At2g33380) | 11,536.06 | 9,246.12 | 0.80 | 7,888.81 | 7,661.10 | 0.97 |
AtABI1 (At4g26080) | 180.31 | 214.81 | 1.19 | 819.79 | 586.52 | 0.72* |
AtABF3 (At4g34000) | 599.62 | 633.01 | 1.06 | 989.19 | 774.77 | 0.78* |
AtABI2 (At5g57050) | 88.21 | 122.78 | 1.39 | 284.56 | 227.19 | 0.80 |
Cold-regulated protein cor15b precursor (At2g42530) | 730.31 | 1,025.37 | 1.40 | 1,488.52 | 1,276.51 | 0.86 |
Alcohol dehydrogenase (At1g77120) | 282.40 | 320.51 | 1.13 | 817.78 | 706.51 | 0.86 |
Putative receptor-like protein kinase RPK1 (At1g69270) | 15.02 | 15.20 | 1.01 | 166.09 | 106.02 | 0.64 |
Hypothetical protein, AtCOR47 (At1g20440) | 107.39 | 98.87 | 0.92 | 185.89 | 161.51 | 0.87 |
Cold-regulated protein COR6.6, KIN2 (At5g15970) | 100.40 | 143.60 | 1.43 | 3,285.65 | 3,414.42 | 1.04 |
Cold-regulated protein cor15a precursor (At2g42540) | 2,841.94 | 2,864.95 | 1.01 | 13,720.20 | 11,956.36 | 0.87 |
ABA-responsive protein-like (At5g13200) | 1,603.03 | 1,162.40 | 0.73 | 7,848.29 | 6,324.54 | 0.81 |
Protein phosphatase 2C (AtP2C-HA) (At1g72770) | 76.42 | 114.33 | 1.50 | 236.49 | 212.89 | 0.90 |
PIP1B aquaporin (At2g45960) | 126.37 | 117.83 | 0.93 | 303.26 | 241.84 | 0.80 |
Gene (locus no.) . | Control . | Salt stress . | ||||
---|---|---|---|---|---|---|
. | Average reads (after normalization) . | Log2fold change . | Average reads (after normalization) . | Log2fold change . | ||
. | WT . | 35SCaMV:OsNUC1-L . | . | WT . | 35SCaMV:OsNUC1-L . | . |
Late embryogenesis abundant protein (At5g06760) | 6,703.68 | 4,298.66 | 0.64 | 12,235.58 | 7,962.26 | 0.65* |
Low-temperature-induced protein 78, RD29A (At5g52310) | 2,101.38 | 1,589.91 | 0.76 | 6,300.14 | 4,290.21 | 0.68* |
Homeobox-leucine zipper protein AtHB-12 (At3g61890) | 180.11 | 249.01 | 1.38 | 187.25 | 177.67 | 0.95 |
AtACP5, acid phosphatase type 5 (At3g17790) | 254.86 | 306.23 | 1.20 | 619.01 | 528.03 | 0.85 |
AtKIN1 (At5g15960) | 614.63 | 363.04 | 0.59* | 1,466.41 | 910.95 | 0.62* |
Dehydrin RAB18-like protein (At5g66400) | 60.79 | 124.36 | 2.05* | 462.96 | 358.55 | 0.77 |
Homeodomain transcription factor AtHB-(At2g46680) | 532.97 | 739.37 | 1.39 | 1,953.16 | 1,729.46 | 0.89 |
AtP5C1, delta-1-pyrroline 5-carboxylase synthetase (At2g39800) | 320.81 | 565.34 | 1.76* | 1,053.88 | 900.09 | 0.85 |
AtERD10, hypothetical protein (At1g20450) | 1,262.55 | 1,959.37 | 1.55* | 4,221.04 | 4,231.15 | 1.00 |
AtRD20 (At2g33380) | 11,536.06 | 9,246.12 | 0.80 | 7,888.81 | 7,661.10 | 0.97 |
AtABI1 (At4g26080) | 180.31 | 214.81 | 1.19 | 819.79 | 586.52 | 0.72* |
AtABF3 (At4g34000) | 599.62 | 633.01 | 1.06 | 989.19 | 774.77 | 0.78* |
AtABI2 (At5g57050) | 88.21 | 122.78 | 1.39 | 284.56 | 227.19 | 0.80 |
Cold-regulated protein cor15b precursor (At2g42530) | 730.31 | 1,025.37 | 1.40 | 1,488.52 | 1,276.51 | 0.86 |
Alcohol dehydrogenase (At1g77120) | 282.40 | 320.51 | 1.13 | 817.78 | 706.51 | 0.86 |
Putative receptor-like protein kinase RPK1 (At1g69270) | 15.02 | 15.20 | 1.01 | 166.09 | 106.02 | 0.64 |
Hypothetical protein, AtCOR47 (At1g20440) | 107.39 | 98.87 | 0.92 | 185.89 | 161.51 | 0.87 |
Cold-regulated protein COR6.6, KIN2 (At5g15970) | 100.40 | 143.60 | 1.43 | 3,285.65 | 3,414.42 | 1.04 |
Cold-regulated protein cor15a precursor (At2g42540) | 2,841.94 | 2,864.95 | 1.01 | 13,720.20 | 11,956.36 | 0.87 |
ABA-responsive protein-like (At5g13200) | 1,603.03 | 1,162.40 | 0.73 | 7,848.29 | 6,324.54 | 0.81 |
Protein phosphatase 2C (AtP2C-HA) (At1g72770) | 76.42 | 114.33 | 1.50 | 236.49 | 212.89 | 0.90 |
PIP1B aquaporin (At2g45960) | 126.37 | 117.83 | 0.93 | 303.26 | 241.84 | 0.80 |
Gene expression was compared using DESeq (Anders and Huber 2010) and adjusting the false discovery rate using the Benjamini–Hochberg method (FDA < 0.05). Fold changes crossing the threshold of P = 0.05 are marked in red and with an asterisk.
OsNUC1-L globally changed the gene expression

Gene enrichment of the significant gene expression difference between the transgenic Arabidopsis with expression of the OsNUC1-L gene and the wild type under normal and salt stress conditions. The enrichment analysis was performed for compartmentalization (A), molecular function (B) and biological function (C) categories. The red bar indicates the higher gene count under salt stress conditions in the transgenic line found in those categories.
To study the global effect of OsNUC1-L gene expression on the transcriptome of Arabidopsis, we used the cytoscape plugin ClueGo (Bindea et al. 2009) to generate and visualize gene set enrichment under normal and salt stress conditions. The significantly enriched terms with more up- and down-regulated genes were visualized in red and blue color, respectively. The node size represents the degree of significance and the color density represents the number of genes in each node.

The gene set enrichment analysis of significantly expressed genes in transgenic OsNUC1-L under normal (A) and salt stress (B) conditions. Nodes in the network show the biological processes significantly enriched, with up- and down-regulated genes shown in red and blue, respectively. The node size represents the degree of significance and the color density the number of genes in each node. The enrichment network was generated by ClueGo (Bindea et al. 2009) with Benjamini and Hochberg correction (corrected P < 0.05).
Under salt stress conditions, the transcriptome analysis still showed the up-regulation in stress-responsive genes and the cellular alkene metabolic process. Moreover, some genes assigned to the auxin metabolic process were also enhanced (Fig. 4B). Although overall genes responding to stress were up-regulated, genes responding to arsenic-containing substances and oxygen level were decreased in the transgenic lines. The down-regulated genes were found in negative regulation of signal transduction, seed maturation, circadian rhythm, positive regulation of post-embryonic development and carbohydrate transport categories (Fig. 4B).
Light-harvesting system protein and photosynthetic membrane components that displayed significant changes in transgenic OsNUC1-L Arabidopsis during normal and salt stress conditions
Experiment . | Location . | Locus name . | Gene name . | Description . | Fold change . |
---|---|---|---|---|---|
Control | PSII | ATCG00020 | PSBA | PSII reaction center protein a | 0.428 |
ATCG00270 | PSBD | PSII reaction center protein d | 0.464 | ||
ATCG00280 | CP43 | CP43 subunit of the PSII reaction center | 0.437 | ||
ATCG00340 | PSAB | D1 subunit of PSI reaction center | 0.566 | ||
ATCG00350 | PSAA | psaA protein comprising the reaction center for PSI | 0.518 | ||
ATCG00630 | PSAJ | Subunit J of PSI | 0.532 | ||
ATCG00680 | PSBB | PSII reaction center protein b | 0.449 | ||
ATCG00710 | PSBH | PSII reaction center protein h | 0.457 | ||
ATCG01060 | PSAC | PsaC subunit of PSI | 0.394 | ||
Cyt b6f complex | AT5G60920 | COB | Glycosylphosphatidylinositol-anchored protein | 0.652 | |
ATCG00540 | PETA | Photosynthetic electron transfer a | 0.640 | ||
ATCG00720 | PETB | Photosynthetic electron transfer b | 0.401 | ||
ATCG00730 | PETD | Photosynthetic electron transfer d | 0.490 | ||
ATMG00160 | COX2 | Cyt c oxidase subunit 2 | 0.534 | ||
ATMG01360 | COX1 | Cyt c oxidase subunit 1 | 0.475 | ||
Salt stress | Light-harvesting complex | AT1G29910 | CAB3 | Chl a/b-binding protein 3 | 1.403 |
AT1G29920 | CAB2 | Chl a/b-binding protein 2 | 1.914 | ||
AT2G05070 | LHCB2.2 | PSII light-harvesting complex gene 2.2 | 1.682 | ||
AT2G05100 | LHCB2.1 | PSII light-harvesting complex gene 2.1 | 1.665 | ||
AT3G27690 | LHCB2.3 | PSII light-harvesting complex gene 2.3 | 1.881 | ||
AT4G31390 | ACDO1 | ABC1-like kinase related to Chl degradation and oxidative stress 1 | 0.683 | ||
AT5G54270 | LHCB3 | Light-harvesting Chl b-binding protein 3 | 1.390 | ||
ATCG00020 | PSBA | PSII reaction center protein A | 0.709 |
Experiment . | Location . | Locus name . | Gene name . | Description . | Fold change . |
---|---|---|---|---|---|
Control | PSII | ATCG00020 | PSBA | PSII reaction center protein a | 0.428 |
ATCG00270 | PSBD | PSII reaction center protein d | 0.464 | ||
ATCG00280 | CP43 | CP43 subunit of the PSII reaction center | 0.437 | ||
ATCG00340 | PSAB | D1 subunit of PSI reaction center | 0.566 | ||
ATCG00350 | PSAA | psaA protein comprising the reaction center for PSI | 0.518 | ||
ATCG00630 | PSAJ | Subunit J of PSI | 0.532 | ||
ATCG00680 | PSBB | PSII reaction center protein b | 0.449 | ||
ATCG00710 | PSBH | PSII reaction center protein h | 0.457 | ||
ATCG01060 | PSAC | PsaC subunit of PSI | 0.394 | ||
Cyt b6f complex | AT5G60920 | COB | Glycosylphosphatidylinositol-anchored protein | 0.652 | |
ATCG00540 | PETA | Photosynthetic electron transfer a | 0.640 | ||
ATCG00720 | PETB | Photosynthetic electron transfer b | 0.401 | ||
ATCG00730 | PETD | Photosynthetic electron transfer d | 0.490 | ||
ATMG00160 | COX2 | Cyt c oxidase subunit 2 | 0.534 | ||
ATMG01360 | COX1 | Cyt c oxidase subunit 1 | 0.475 | ||
Salt stress | Light-harvesting complex | AT1G29910 | CAB3 | Chl a/b-binding protein 3 | 1.403 |
AT1G29920 | CAB2 | Chl a/b-binding protein 2 | 1.914 | ||
AT2G05070 | LHCB2.2 | PSII light-harvesting complex gene 2.2 | 1.682 | ||
AT2G05100 | LHCB2.1 | PSII light-harvesting complex gene 2.1 | 1.665 | ||
AT3G27690 | LHCB2.3 | PSII light-harvesting complex gene 2.3 | 1.881 | ||
AT4G31390 | ACDO1 | ABC1-like kinase related to Chl degradation and oxidative stress 1 | 0.683 | ||
AT5G54270 | LHCB3 | Light-harvesting Chl b-binding protein 3 | 1.390 | ||
ATCG00020 | PSBA | PSII reaction center protein A | 0.709 |
Light-harvesting system protein and photosynthetic membrane components that displayed significant changes in transgenic OsNUC1-L Arabidopsis during normal and salt stress conditions
Experiment . | Location . | Locus name . | Gene name . | Description . | Fold change . |
---|---|---|---|---|---|
Control | PSII | ATCG00020 | PSBA | PSII reaction center protein a | 0.428 |
ATCG00270 | PSBD | PSII reaction center protein d | 0.464 | ||
ATCG00280 | CP43 | CP43 subunit of the PSII reaction center | 0.437 | ||
ATCG00340 | PSAB | D1 subunit of PSI reaction center | 0.566 | ||
ATCG00350 | PSAA | psaA protein comprising the reaction center for PSI | 0.518 | ||
ATCG00630 | PSAJ | Subunit J of PSI | 0.532 | ||
ATCG00680 | PSBB | PSII reaction center protein b | 0.449 | ||
ATCG00710 | PSBH | PSII reaction center protein h | 0.457 | ||
ATCG01060 | PSAC | PsaC subunit of PSI | 0.394 | ||
Cyt b6f complex | AT5G60920 | COB | Glycosylphosphatidylinositol-anchored protein | 0.652 | |
ATCG00540 | PETA | Photosynthetic electron transfer a | 0.640 | ||
ATCG00720 | PETB | Photosynthetic electron transfer b | 0.401 | ||
ATCG00730 | PETD | Photosynthetic electron transfer d | 0.490 | ||
ATMG00160 | COX2 | Cyt c oxidase subunit 2 | 0.534 | ||
ATMG01360 | COX1 | Cyt c oxidase subunit 1 | 0.475 | ||
Salt stress | Light-harvesting complex | AT1G29910 | CAB3 | Chl a/b-binding protein 3 | 1.403 |
AT1G29920 | CAB2 | Chl a/b-binding protein 2 | 1.914 | ||
AT2G05070 | LHCB2.2 | PSII light-harvesting complex gene 2.2 | 1.682 | ||
AT2G05100 | LHCB2.1 | PSII light-harvesting complex gene 2.1 | 1.665 | ||
AT3G27690 | LHCB2.3 | PSII light-harvesting complex gene 2.3 | 1.881 | ||
AT4G31390 | ACDO1 | ABC1-like kinase related to Chl degradation and oxidative stress 1 | 0.683 | ||
AT5G54270 | LHCB3 | Light-harvesting Chl b-binding protein 3 | 1.390 | ||
ATCG00020 | PSBA | PSII reaction center protein A | 0.709 |
Experiment . | Location . | Locus name . | Gene name . | Description . | Fold change . |
---|---|---|---|---|---|
Control | PSII | ATCG00020 | PSBA | PSII reaction center protein a | 0.428 |
ATCG00270 | PSBD | PSII reaction center protein d | 0.464 | ||
ATCG00280 | CP43 | CP43 subunit of the PSII reaction center | 0.437 | ||
ATCG00340 | PSAB | D1 subunit of PSI reaction center | 0.566 | ||
ATCG00350 | PSAA | psaA protein comprising the reaction center for PSI | 0.518 | ||
ATCG00630 | PSAJ | Subunit J of PSI | 0.532 | ||
ATCG00680 | PSBB | PSII reaction center protein b | 0.449 | ||
ATCG00710 | PSBH | PSII reaction center protein h | 0.457 | ||
ATCG01060 | PSAC | PsaC subunit of PSI | 0.394 | ||
Cyt b6f complex | AT5G60920 | COB | Glycosylphosphatidylinositol-anchored protein | 0.652 | |
ATCG00540 | PETA | Photosynthetic electron transfer a | 0.640 | ||
ATCG00720 | PETB | Photosynthetic electron transfer b | 0.401 | ||
ATCG00730 | PETD | Photosynthetic electron transfer d | 0.490 | ||
ATMG00160 | COX2 | Cyt c oxidase subunit 2 | 0.534 | ||
ATMG01360 | COX1 | Cyt c oxidase subunit 1 | 0.475 | ||
Salt stress | Light-harvesting complex | AT1G29910 | CAB3 | Chl a/b-binding protein 3 | 1.403 |
AT1G29920 | CAB2 | Chl a/b-binding protein 2 | 1.914 | ||
AT2G05070 | LHCB2.2 | PSII light-harvesting complex gene 2.2 | 1.682 | ||
AT2G05100 | LHCB2.1 | PSII light-harvesting complex gene 2.1 | 1.665 | ||
AT3G27690 | LHCB2.3 | PSII light-harvesting complex gene 2.3 | 1.881 | ||
AT4G31390 | ACDO1 | ABC1-like kinase related to Chl degradation and oxidative stress 1 | 0.683 | ||
AT5G54270 | LHCB3 | Light-harvesting Chl b-binding protein 3 | 1.390 | ||
ATCG00020 | PSBA | PSII reaction center protein A | 0.709 |

The gene set enrichment analysis of significantly regulated genes in transgenic OsNUC1-L under normal (A) and salt stress (B) conditions. Nodes in the network show the specific cellular comparment, where the significantly enriched terms with more up- and down-regulated genes are shown in red and blue, respectively. The node size represents the degree of significance and the color gradient represents the proportion of genes in each node. The enrichment network was generated by ClueGo (Bindea et al. 2009) with Benjamini and Hochberg correction (corrected P < 0.05).
OsNUC1-S has strong effects on gene expression during salt stress

Gene enrichment of the significant gene expression difference between the transgenic Arabidopsis with expression of the OsNUC1-S gene and the wild type under normal and salt stress conditions. The enrichment analysis was performed for compartmentalization (A), molecular function (B) and biological function (C) categories.

The gene set enrichment analysis of significantly expressed genes in transgenic OsNUC1-S under salt stress conditions. The enrichment network was generated by ClueGo (Bindea et al. 2009) with Benjamini and Hochberg correction (corrected P < 0.01). The node size represents the number of gene in each node, and the darker colors represent higher siginificant levels.

The gene set enrichment analysis of significantly regulated genes in transgenic OsNUC1-S under salt stress conditions, showing cellular compartment (A) and molecular function (B). Nodes in the network show the specific cellular comparment. The node sizes represent the degree of significance.The enrichment network was generated by ClueGo (Bindea et al. 2009) with Benjamini and Hochberg correction (P < 0.01).
OsNUC1-L improved photosynthetic activity during salt stress by regulating the expression of genes in the thylakoid membrane

Changes in photosynthesis parameters: (A) net photosynthesis, (B) stomatal conductivity, (C) transpiration rate, (D) intracellular CO2 concentration and (E) water use efficiency, measured before and after 1 d of salt treatment (150 mM NaCl). Means were compared using one-way analysis of variance. All parameters ae shown as the mean ± SE. * indicates significant difference from the wild type (WT) (P < 0.05, Duncan’s Multiple Range Test). Nine plants were used for the analysis with a completely randomized design (CRD).
Transcriptome analysis of Arabidopsis with OsNUC1-S expression also suggested changes in the chloroplast gene function (Fig. 8A, B) and photosynthetic process (Fig. 7). Unfortunately, we could not evaluate the photosynthetic responses in Arabidopsis with OsNUC1-S expression due to irregular inheritance of the transgene in the line selected for advancement.
Overexpression of OsNUC1-L enhances photosynthesis under salt stress conditions in rice

OsNUC1 relative gene expression in wild-type (WT) (A) and OsNUC1-L-overexpressing transgenic lines, FL1 (B), FL2 (C) and FL3 (D) under normal and salt stress conditions. The experiment was performed with three biological replicates. Relative gene expression was calculated in comparison with the expression of the wild type on Day 0.

Photosynthesis rate (A, B) and stomatal conductance (C and D) of flag leaves (A and C) and second leaves (B and D) of wild-type and transgenic rice, FL1, FL2 and FL3 at the booting stage grown in normal conditions or treated with salt stress for 9 d. The experiment was designed in random complete block design with four replications. Three individuals of each homozygous T3 transgenic line and the wild type were used for each replication. Means with a different letter above the bar are significantly different (P < 0.05, Duncan’s Multiple Range Test).

Photosynthetic pigment content, Chl a (A and B), Chl b (C and D) and carotenoid (E and F), of flag leaves (A, C, E) and second leaves (B, D, F) of wild-type and transgenic rice with OsNUC1-L overexpression, when rice plants were treated with salt stress at 8–10 dS m–1 for 9 d at the booting stage. The experiment was designed in a random complete block design with four replications. Three individuals of each homozygous T3 transgenic line and the wild type were used for each replication. Means with a different letter above the bar are significantly different (P < 0.05, Duncan’s Multiple Range Test).

PSII efficiency determined by Chl fluorescence (Fv/Fm) after 30 min of dark adaptation in flag leaves (A) and second leaves (B) of wild-type (WT) and transgenic rice lines, FL1, FL2 and FL3 grown in normal conditions or treated with salt stress for 9 d. The experiment was designed in a random complete block design with four replications. Three individuals of each homozygous T3 transgenic line and the wild type were used for each replication. Means with a different letter above the bar are significantly different (P < 0.05, Duncan’s Multiple Range Test).
To explore the correlation between the changes in the photosynthesis response and gene expression in the OsNUC1-L overexpression rice lines, we performed quantitative RT-PCR of the light reaction-associated genes that are homologous to the Arabidopsis genes displaying significant differences upon OsNUC1-L overexpression. Four rice genes, LOC_Os01g52240.1, encoding a CAB, LOCOS01g57942.1 encoding Cyt b6, LOC_Os03g39610.1 encoding a CAB and LOC_Os08g35420.1 encoding a photosynthetic reaction center protein, which respectively are homologous to AT1G29920 (CAB2), ATCG00720 (PETB), AT2G05070 (LHCB2.2) and ATCG00020 (PSBA), were examined in the flag leaf and second leaf of transgenic rice in comparison with the wild type.

Relative gene expression of four rice genes, LOC_Os01g52240.1(A), encoding a CAB, LOCOS01g57942.1 (B) encoding Cyt b6, LOC_Os03g39610.1 (C) encoding a CAB, and LOC_Os08g35420.1 (D) encoding a photosynthetic reaction center protein in transgenic rice (FL3) at the flowering stage under normal conditions or 9 d after salt stress. The experiment was performed using the T3 transgenic line (FL3), with three biological replicates of three technical replicates used for the quantiative PCR. The non-transgenic wild type (WT) was used as controls. Means with a different letter above the bar are significantly different (P < 0.05, Duncan’s Multiple Range Test).
Discussion
A considerable amount of information is available on the biochemistry and function of NUC. By ultrastructural cell analysis, NUC is mostly found in the nucleolus, particularly around FCs and the DFC (Tong et al. 1997, Pontvianne et al. 2007, Stepinski 2012), where it is involved in rRNA transcription and pre-ribosomal formation (Medina et al. 2010). Yeast and Arabidopsis mutants lacking NUC function display overproduction of unprocessed rRNA (Lee et al. 1992. Petricka and Nelson 2007. Pontvianne et al. 2007). At the same time, NUC also affects the 40S/60S ribosomal subunit ratio. Yeast defective in NUC gene function displayed a dramatic decrease in the 40S/60S ribosomal subunit ratio (Kondo and Inouye 1992). Thus, NUC is associated with rRNA transcription and maturation, as well as ribosome biogenesis. Consistent with this association, the level of NUC expression is directly related to cell activity and protein translation (Stepinski 2012). In alfalfa, NUC expression was absent in stationary phase cells, but high in dividing cells (Bogre et al. 1996).
The two NUC splice isoforms appear to have at least some separate functions. Transgenic expression of the large isoform of rice nucleolin in A. thaliana, OsNUC1-L, did not prevent dry weight loss under salt stress, in contrast to the effect found when overexpressing OsNUC1-S (Sripinyowanich et al. 2013). The subcellular localization of green fluorescent protein (GFP)–OsNUC1 fusion protein showed that both OsNUC1-L and OsNUC1-S could be imported into the nucleus. The difference in the phenotype of the OsNUC1-L and OsNUC1-S expression lines may result from differential binding of RNA targets (Ginisty et al. 1999). However, we cannot rule out the possibility that OsNUC1-L and OsNUC1-S can regulate transcription itself. Surprisingly, OsNUC1-L was more stable in transgenic Arabidopsis, while OsNUC1-S inheritance appeared compromised, limiting a complete examination of the physiological response to Arabidopsis OsNUC1-S transgenics.
The expression of the OsNUC1-L gene in A. thaliana resulted in both up- and down-regulation of multiple mRNAs. The mechanism through which NUC regulates the mRNA level of target genes is still unclear. NUC is a multifunctional protein involved in multiple regulatory pathways. The well-documented RNA-binding activity of NUC might affect RNA turnover and thus abundance. Binding of NUC to DNA has also been documented (Olson et al. 1983, Bouche et al. 1984, Sapp et al. 1986) as well as its ability to interact with histones and induce chromatin decondensation (Erard et al. 1988, Kharrat et al. 1991, Gaume et al. 2011, Goldstein et al. 2013). A chromatin-modifying role for NUC has been well documented in rDNA transcription (reviewed by Durut and Sáez-Vásquez 2015). Therefore, OsNUC1 may regulate the transcription of target genes by modifying local chromatin properties. Moreover, transgenic expression of OsNUC1-L during salt stress up-regulated a large group of genes involved in ribosome biogenesis (Fig. 5), suggesting that OsNUC1-L stabilized growth under saline conditions by stimulating the translation machinery, a role consistent with known NUC properties (Ginisty et al. 1999).
An interesting effect of OsNUC1-L expression in transgenic Arabidopsis was hypersensitivity to the hormone ABA. Under an ABA concentration that delays germination of wild-type seed, the OsNUC1-L transgenic lines were completely inhibited. This phenomenon is similar to that displayed by mutants of genes encoding proteins with RNA binding activities (Kuhn and Schroeder 2003). Increased responsivenes to ABA, such as delayed seed germination, root growth inhibition and rapid stomatal closure upon ABA application, link this class of proteins to ABA signal transduction. One such gene, ABA-HYPERSENSITIVE 1 (ABH1), encodes the large subunit of the nuclear cap-binding protein complex that was showed to participate in several steps of RNA processing, nuclear export and mRNA decay (Makarov et al. 2002). The abh1 knock-out mutant, identified in a screen for altered ABA responses, displays delayed seed germination, rapid stomatal closure and a reduced wilting phenotype during drought stress (Hugouvieux et al. 2001) as well as altered post-transcriptional RNA processing (Kuhn and Schroeder 2003). Another example of an RNA-binding protein affecting the ABA response and plant stress is provided by the glycine-rich RNA-binding protein, ATGRP7-1, whose mutation in A. thaliana results in increased expressoin of ABA-induced genes and in other phenotypes related to those displayed by NUC mutants. Overexpression of ATGRP7-1 increases cold tolerance (Cao et al. 2006, Kim et al. 2008, Kim et al. 2010). Complementation of the atgrp7 mutant with OsGRP1 and OSGRP4 restored mRNA export from the nucleus to the cytosol, supporting broad conservation of NUC and a relationship between RNA metabolism, hormone response and stress adaptation.
Chloroplasts and ribosomes are most directly affected by salt stress. Salt stress and high light intensity induce photoinhibition (Takahashi and Murata 2008), probably resulting from the balance between the rate of photodamage to and rate of repair of PSII (Allakhverdiev et al. 2005). Hypersaline conditions suppress multiple pathways: CO2 fixation by inhibiting the activity of ribulose-1,5-bisphosphate carboxylase oxygenase (Solomon et al. 1994), protein translation by reducing ribosome stability (Brady et al. 1984) and the repair mechanism of PSII by inhibiting D1 protein synthesis (Allakhverdiev et al.2005). The NUC gene is strongly induced after exposure to light (Tong et al. 1997, Reichler et al. 2001) consistent with a role for NUC in the high light response. Gene enrichment analysis of salt-treated Arabidopsis OsNUC1-L transgenics revealed that genes with increased expression encode proteins important in chloroplast and ribosome function (Fig. 5B). OsNUC1-L up-regulation of these genes may contribute to the stability of these organelles.
During salt stress conditions, the net photosynthesis rate in OsNUC1-L transgenic lines was significantly increased compared with wild-type Arabidopsis (Fig. 9A). Moreover, transcriptomic analysis revealed that during salt stress the OsNUC1-L transgenic lines increased expression of the genes encoding PSII proteins, such as the PSII light-harvesting complex genes (LHCB2.1, LHCB2.2 and LHCB2.3) genes, significantly enhancing light-harvesting system protein and photosynthetic membrane components (Table 3). Genes encoding CAB2 and CAB3 (AT1G29910 and AT1G29920) and proteins of the PSII light-harvesting complex were also up-regulated, as were genes encoding photosynthetic membrane components (Table 3). The CABs (also LHCB for light-harvesting CAB) are the main component in photosystems and the most abundant proteins in plant green tissues. The T-DNA knockouts of some LHCB genes increased stomatal opening and sensitivity to drought stress (Xu et al. 2012). Plausibly, increased expression of the LHCB genes could increase tolerance to stress. Consistent with this, salt tolerance in Indian mustard was associated with higher expression of proteins involved in stabilization of photosystems (Yousuf et al. 2015, Yousuf et al. 2016), and higher salt tolerance of ‘Pokkali’ rice was associated with a higher level of proteins involved in the light reaction of photosynthesis (Saleethong et al. 2016). Higher expression of OsNUC1 is associated with the higher salt tolerance of this rice (Sripinyowanich et al. 2013). To validate the role of OsNUC1-L, we produced and characterized transgenic rice lines overxpressing OsNUC1-L. Similar to transgenic Arabidopsis expressing OsNUC1-L, we detected a higher photosynthetic rate under salt stress in this transgenic rice at the flowering stage, although not in seedlings. Based on detailed gene expression analysis of the significantly regulated photosynthesis genes in the RNA-seq data, all tested genes were expressed at a higher level in the flag leaves of the transgenic lines in salt stress conditions. This supports the role of OsNUC1-L in enhancement of photosynthesis via the protection of light reaction components under salt stress. Decreased expression of the tested genes in the second leaf suggests that the salt stress effects may be strong in that tissue. Rice flag leaf is the major photosynthetic organ responsible for grain filling. More than half of the carbohydrates in rice seeds are from flag leaf photosynthesis (Li et al. 1998). Therefore, the differential protection from salt stress may occur between the flag and second leaf. These results indicate that the protection of photosynthesis under salt stress found in OsNUC1-L-expressing Arabidopsis is not a novel and spurious NUC function acquired in the Arabidopsis heterologous system, but rather it derives from a conserved role for OsNUC1-L in photosynthesis stabilization under salt stress conditions.
In conclusion, we explored NUC function in the context of environmental stress by expressing OsNUC1-L in the background of wild-type A. thaliana. Plants expressing the rice NUC displayed altered regulation of the genes encoding key enzymes in ABA biosynthesis and degradation such as NCED3, AAO3, ABA3 and CYP707A4. Consistent with a role for NUC in ABA signaling, they also displayed changes in the expression level of some ABA-responsive genes. These changes had phenotypic consequences as the rice NUC expressors displayed hypersensitivity to ABA during germination. A question remaining to be addressed is whether this response is due to the changes in ABA accumulation or ABA sensitivity. Our transcriptome analysis of rice NUC transgenic plants indicated that OsNUC1-L plays a role in RNA metabolism and specifically regulates the expression level of genes involved in stress responses, photosynthesis, hormone response and transport. The use of cellular compartment criteria connected OsNUC1-L expression to up-regulation of genes in photosynthesis, especially in the light-harvesting complex. Our results showed that during stress conditions OsNUC1-L increased the transcript level of LCHBs proteins, leading to maintenance of photosynthetic capacity, increasing WUE and growth. The higher capacity for photosynthesis maintenance under salt stress conditions was found in both Arabidopsis and rice with OsNUC1-L overexpression. We conclude that OsNUC1 is involved in photosystem stabilization during salt stress conditions.
Materials and Methods
Generation of transgenic OsNUC1-L- and OsNUC1-S-expressing Arabidopsis
Full-length OsNUC1-L and OsNUC1-S cDNA was PCR amplified, cloned into the pCR®- Blunt II-TOPO® vector (Invitrogen) and sequenced. Then it was subcloned into the plant gene transformation vector, pJim19 plasmid, to generate the expression construct under the Cauliflower mosaic virus 35S promoter (35SCaMV promoter). The recombinant plasmid was used to transfer the genes into wild-type Arabidopsis (Arabidopsis thaliana ecotype Columbia-0) by floral dipping (Clough and Bent 2008). The kanamycin resistance gene was used to select the successfully transformed plants. The T2 homozygous transgenics were used for further experiments.
Generation of OsNUC1-L-overexpressing transgenic rice
The OsNUC1-L cDNA was cloned into the pWM101 vector to generate the overexpression construct of OsNUC1-L under the ubiquitin promoter. Transformation was performed according to Sripinyowanich et al. (2013). Briefly, rice calli were regenerated from scutellum on MS basal medium containing 2.5 mg l–1 2,4-D and 3% (w/v) sucrose. Then, Agrobacterium transformation was done to transfer the required construct to rice calli. The transformed calli were regenerated to seedlings and selected on MS basal medium containing 50 mg l–1 hygromycin B. Three T3 homozygous transgenic rice lines, FL1, FL2 and FL3, were used for characterization.
Germination test
After surface sterilization, seeds were placed on half-strength MS medium (Murashige and Skoog 1962), stratified at 4°C for 2 d and then transferred to growth chambers at 25°C, with a light intensity of 100 µmol m–1 s–1 and 16 h days/8 h nights. To determine ABA sensitivity during germination, seeds of OsNUC1-L transgenic and wild-type Arabidopsis were germinated on half-strength MS medium with or without 0.2, 0.4 and 0.6 µM ABA (Sigma).
OsNUC1-L expression quantification
The gene expression of OsNUC1 in 14-day-old rice seedlings grown under normal or salt stress conditions for 6 d was detected with the method according to Sripinyowanich et al (2013). The OsEf-1α gene (GenBank accession No. GQ848074.1) was used as internal control. The relative expression was calculated with the above formula.
Transcriptome study using RNA-seq
Three biological replicates of the wild-type and transgenic Arabidopsis were used for transcriptome analysis. Two-week-old plants grown on MS medium were transferred to MS medium containing 100 mM NaCl for salt stress conditions and to freshly perpared MS medium for control conditions. After 5 h of treatment, total RNA was isolated with Plant RNA Reagent (Invitrogen). The mRNA was uniquely isolated with Dynabead mRNA purification (Invitrogen) and converted to cDNA with SuperScript III (Invitrogen). The cDNA was fragmented by dsDNA fragmentase (NEB). Fragment sizes of 300 bp were selected and purified using the AMPure protocol. The DNA fragments were enriched by PCR for 10 cycles using adaptors with unique DNA barcodes for each cDNA library. All six libraries were generated using the Genome Analyzer (Illumina), loading them together in the same sequencing lane of the flow cell.
The transcriptome analysis compared transcripts from the wild-type and transgenic plants under salt stress and control conditions. Briefly, all short-sequence reads produced from Genome Analyzer were processed to remove the DNA barcode and grouped according to the DNA barcode in the right category using the pipeline created by Missirian et al. (2012). The sequence reads were uniquely aligned and mapped to the Arabidopsis genome (TAIR10) using Bowtie aligner (Langmead et al. 2009) and TopHat (Trapnell et al. 2009). The DESeq program (Anders and Huber 2010) was used to identify differentially expressed genes. Significantly regulated genes were considered differentially expressed.
Determination of gas exchange parameters and Chl fluorescence
To examine the effect of salt stress on photosynthetic activity in OsNUC1-L-overexpressing transgenic lines, the seeds of the wild type and three 35SCaMV:OsNUC1-L transgenic lines were grown on soil mix for 4 weeks at 25°C, 180 µmol photon m–1 s–1 and 16 h days/8 h nights. The net photosynthetic rate, stomatal conductance and transpiration rate were determined using a Gas Analysis System (Li-Cor, LI-6400) with a 6400-15 extended reach 1 cm diameter chamber. The gas exchange measurement was performed at 25°C under an external light source of 650 µmol photon m–1 s–1, the reference carbon dioxide concentration of 350 µmol m–1 s–1 and relative humidity of 40–44%. For the salt treatment, wild-type and transgenic Arabidopsis were exposed to 25 ml of 100 mM NaCl per pot by irrigation. The WUE was calculated by the formula is given as: WUE = Pn/E, where Pn = net rate of photosynthesis (µmol CO2 m–1 s–1) and E = transpiration rate (mmol m–2 s–1)
The rate of photosynthesis and stomatal conductance were determined in rice at the booting stage. At this stage, rice leaves have been enlarged which is better for the measurement of photosyntheic behavior than seedling leaves. Rice plants were grown on soil in 30 cm diameter pots until flowering in the greenhouse under natural light conditions. For salt stress treatment, when rice plants reached the booting stage, 150 mM NaCl solution was poured to cover the soil and the soil electric conductivity was determined to reach 8–10 dS m–1 during salt stress treatment. The photosynthetic rate and stomatal conductance were determined a Gas Analysis System (Li-Cor, LI-6400) after 9 d of treatment. The gas exchange measurement was performed at 30–33°C under an external light source of 1,200 µmol photon m–1 s–1, the reference carbon dioxide concentration of 380 µmol m–1 s–1 and relative humidity of 56–67%.
Fv/Fm was determined in flag leaves and second leaves of normal and salt-stressed plants using Pocket PEA (Hansatech Instruments, Ltd.) after 30 min of dark adaptation.
Determination of the content of photosynthetic pigments
Photosynthetic pigment concentration was determined by the methods of Shabala et al. (1998) and Lichtenthaler (1987). The pigments were extracted with 80% acetone, and Chl a, Chl b and carotenoids were measured using a UV–visible spectrophotometer (Agilent Technology) at absorbance wavelengths of 662, 644 and 470 nm. The pigment contents were calculated with the following equations: [Chl a] = 9:784D662 – 0:99D644; [Chl b] = 21:42D664 – 4:65D662; and [carotenoid] = (1,000D470 – 1:90[Chl a] – 63:14[Chl b])/214.
Supplementary data
Supplementary data are available at PCP online.
Funding
This research was supported by the Ratchadaphiseksomphot Endowment Fund (2015–16) and Chulalongkorn University [RES560530061-FW, CU-58-012-FW, and, CUAASC; Ratchadaphiseksomphot Fund for Postdoctoral Fellowship; Chulalongkorn University [to T.U. and S.S]; and the Science Achievement Scholarship of Thailand [to C.B.].
Abbreviations
- AAO3
aldehyde oxidase 3
- ABH1
ABA hypersensitive 1
- AtEF-1α
Eukaryotic elongation factor1-α subunit
- CAB
Chl a/b-binding protein
- 35SCaMV promoter
Cauliflower mosaic virus 35S promoter
- CBP20
nuclear cap-binding protein subunit 2
- Ci
intracellular CO2 concentration
- CYP707A4
abscisic acid 8′-hydroxylase 4
- DFC
dense fibrillar component
- E
transpiration rate
- FC
fibrillar center
- GAR
glycine–arginine rich
- GO
Gene Ontology
- gs
stomatal conductance
- HYL1
Hyponastic leaves1 mutant
- LHCB
PSII light-harvesting complex
- MS
Murashige and Skoog
- NCED
9-cis-epoxycarotenoid dioxygenase
- NLS
nuclear localization signal
- NUC
nucleolin
- OsNUC1
Rice nucleolin 1
- Pn
net photosynthetic rate
- PSBA
PSII reaction center protein
- RRM
RNA-recognition motif
- RT-PCR
real-time PCR
- WUE
water use efficiency
Acknowledgments
We would like to thanks Dr. Isabelle M. Henry and Ms. Kathie Ngo, Genome Center, University of California Davis for their valuable help in bioinformatic analysis.
Disclosures
The authors have no conflicts of interest to declare.