Skip to main content

Comparative transcriptomic analysis of the super hybrid rice Chaoyouqianhao under salt stress

Abstract

Background

Soil salinization is a threat to food security. China is rich in saline land resources for potential and current utilization. The cultivation and promotion of salt-tolerant rice varieties can greatly improve the utilization of this saline land. The super hybrid rice Chaoyouqianhao (CY1000) is one of the most salt-tolerant rice varieties and is widely used, but the molecular mechanism underlying its salt tolerance is not clear.

Results

In this study, the characteristics of CY1000 and its parents were evaluated in the field and laboratory. The results showed that aboveground parts of CY1000 were barely influenced by salt stress, while the roots were less affected than those of its parents. A comparative transcriptomic strategy was used to analyze the differences in the response to salt stress among the male and female parents of CY1000 at the seedling stage and the model indica rice 93–11. We found that the salt tolerance of CY1000 was mainly inherited from its male parent R900, and its female parent GX24S showed hardly any salt tolerance. To adapt to salt stress, CY1000 and R900 upregulated the expression of genes associated with soluble component synthesis and cell wall synthesis and other related genes and downregulated the expression of most genes related to growth material acquisition and consumption. In CY1000 and R900, the expression of genes encoding some novel key proteins in the ubiquitination pathway was significantly upregulated. After treatment with MG-132, the salt tolerance of CY1000 and R900 was significantly decreased and was almost the same as that of the wild type after salt stress treatment, indicating that ubiquitination played an important role in the salt tolerance mechanism of CY1000. At the same time, we found that some transcription factors were also involved in the salt stress response, with some transcription factors responding only in hybrid CY1000, suggesting that salt tolerance heterosis might be regulated by transcription factors in rice.

Conclusion

Our results revealed that the ubiquitination pathway is important for salt tolerance in rice, and several novel candidate genes were identified to reveal a novel salt tolerance regulation network. Additionally, our work will help clarify the mechanism of heterosis in rice. Further exploration of the molecular mechanism underlying the salt tolerance of CY1000 can provide a theoretical basis for breeding new salt-tolerant rice varieties.

Peer Review reports

Introduction

Soil salinization and secondary salinization have become global ecological problems affecting cultivated land and threatening world food security [1]. China’s inland saline-alkali land area is nearly 100 million hectares, and the coastal beach area is nearly 2.333 million hectares [2]. Eighty percent of saline-alkali land is in a state of desolation. In total, 13.33 million hectares of this saline-alkali land is rich in rainwater resources and can be used for rice cultivation in the near future, with great potential for comprehensive utilization. This will be the most economical and effective way to develop and utilize coastal beaches and inland saline-alkali land and guarantee land area for the cultivation and popularization of saline-alkali-tolerant rice varieties [3].

Under the leadership of academician Yuan Longping, the National Hybrid Rice Research Centre has selected a number of salt-tolerant rice varieties with adaptation to salinity levels of 0.3–0.6% and yields of over 4500 kg/ha. Among them, the super hybrid rice variety “Chaoyouqianhao” (CY1000) is the best and is widely used [4]. CY1000, belonging to the fifth generation of super hybrid rice, is an indica two-line hybrid rice variety. It is a cross between the super rice variety R900 and photothermosensitive two-line male sterile line Guangxiang24S (GX24S). CY1000 passed the national examination and approval (national rice examination 20,170,053) in 2017. Although the salt tolerance advantage of CY1000 is being gradually utilized, its molecular mechanism of salt tolerance remains unclear.

To date, studies have found that the main mechanisms of salt resistance are ion balance regulation [5, 6], osmotic balance regulation [7, 8], reactive oxygen species scavenging [9] and nutrient balance regulation [10, 11]. These regulatory effects are mainly mediated by the salt overly sensitive (SOS) pathway and plant hormone signal transduction [12]. When rice is exposed to high salinity, a large amount of Na+ enters the cells through univalent cation channels, which stimulates a rapid increase in the intracellular Ca2+ concentration [13]. Then, the calcium-binding protein OsSOS3/CBL4 senses this signal, interacts with and phosphorylates OsSOS2/OsCIPK24, a member of the SnRK3 family, and rapidly activates the Na+/H+ antiporter OsSOS1 on the plasma membrane, which promotes Na+ excretion from cells and maintains ion balance in cells [14, 15]. This is the primary mechanism underlying the response of rice to salt stress. Subsequently, plant hormones (such as ABA and ethylene) and epigenetic regulation also respond to salt stress signals and play an important role in plant salt tolerance [11, 16, 17].

To explore the salt tolerance mechanism of CY1000, we used comparative transcriptomics to examine the response of CY1000 to salt stress at the seedling stage, analyzed the differences in gene expression among male parents, female parents and the model indica rice Line 93–11, and analyzed the molecular mechanism of CY1000 salt tolerance to provide a theoretical basis for breeding new salt-tolerant rice varieties.

Results

The field planting results of CY1000 showed that it had good salt tolerance and high yield in high-salt fields

From 2018 to 2020, field experiments were carried out in Sanya city, Hainan Province (109.17°E, 18.36°N), for three consecutive years (Fig. 1). CY1000 exhibited advantages in seed setting rate, yield under 0.3% salt stress and yield reduction under 0.3% salt stress (Fig. 1B). In 2020, 13.3 ha of CY1000 was planted and subjected to 0.3% salinity, and the average yield was 7.6 T/ha (Fig. 1A).

Fig. 1
figure 1

Yield performance of CY1000 under salt stress. A The yield of different varieties under 0.3% salt stress. B Salt stress (0.3%) reduced the yield of different varieties

CY1000 and its male parent R900 had high salt tolerance, but its female parent GX24S showed almost no salt tolerance

To further understand the salt resistance mechanism of CY1000, the responses of the parents and the model species 93–11 to salt stress at the seedling stage were analyzed. When 7-day-old seedlings of CY1000, the parents R900 and GX24S, and 93–11 were planted in 0.6% NaCl, the growth of the four varieties was inhibited to varying degrees after 7 d. CY1000 showed the strongest salt tolerance (Fig. 2A-B); the survival rate was still 98% after 7 d, and the stem length was barely different from that of untreated CY1000 (Fig. 2B-C). The growth of the female parent GX24S and that of 93–11 was significantly inhibited; the survival rate decreased to below 40%, and the growth of the stems and roots was significantly affected, showing a state of stopped growth. In addition, the male parent R900 also showed strong salt tolerance; the survival rate reached 85%, and the length of the stem showed no significant difference compared to that of the untreated materials (Fig. 2B&D). Interestingly, the roots of CY1000 and R900 were also influenced by salt stress, but they also exhibited no growth inhibition. After careful observation, it was found that there were obvious salt crystals at the stem nodes of CY1000 seedlings, which indicates that the balance of the Na+ emission concentration is one of the important factors underlying salt resistance.

Fig. 2
figure 2

Results of salt tolerance treatment of seedlings of CY1000 and its parents. A Photos of 7 d salt-treated and control rice seedlings (scale bar, 1 cm). B Survival rate of CY1000 and its parents after 7 d of salt treatment. C The length of roots and stems of three rice seedlings after 7 d of salt treatment. Each value was verified by more than 20 samples. The p value was tested by Student’s t test

Comparative omics revealed the molecular mechanism by which CY1000 and its parents respond to salt stress

In total, 6.3 Gb of data was obtained for each replicate of each sample, with Q20 over 95.8% and Q30 over 92%. The total amount of data was more than 163 Gb. More than 106.8 Gb of clean reads was obtained. When using parametric analysis, the proportion of the total mapped genes in the genome was more than 91% in most samples; only the proportion of R_Un_1 was 83.66%. Further statistics on the distribution of reads on the genome showed that 88.45% of the reads were mapped to genes, and 96.27% of the reads were mapped to exons. The PCA and Pearson correlation results are shown in Fig. S1. Although a few replicates were not most correlated, they all clustered separately. RNA sequencing (RNA-Seq) showed that the number of differentially expressed genes in 93–11 and GX24S was quite different (Fig. S1, Table S1). GX24S had the fewest (2940 DEGs, including 1685 downregulated and 1255 upregulated genes) differentially expressed genes, while 93–11 had the highest number (5809 DEGs, including 2793 downregulated and 3016 upregulated genes) of differentially expressed genes; CY1000 had 5261 DEGs, including 2828 downregulated and 2433 upregulated genes; and R900 had 4114 DEGs, including 2292 downregulated and 1822 upregulated genes (Fig. 3A, Fig. S1B-E).

Fig. 3
figure 3

Differences in CY1000 and other varieties under salt treatment. A Heatmap of all expressed genes from samples. B GO and KEGG functional analyses of the three varieties under salt treatment. C Venn diagram of GO terms among samples. D Observation of the salt tolerance and phenotype of CY1000. Salt crystals were observed in the nodes of CY1000 stems

After clustering the expression patterns, all samples were divided into two branches according to treatment and no treatment, indicating that there were common mitigation patterns in the salt stress response among different varieties at different levels (Tables S2, S3, S4 and S5). The response patterns of GX24S were the most similar, indicating that they were insensitive to salt stress, which may be the reason for the poor salt resistance. CY1000 and R900 had more similar expression spectra, which indicates that the salt resistance of CY1000 was inherited from R900.

GO and KEGG functional enrichment analyses also showed that CY1000 and R900 had more similar response patterns, including “phenylpropanoid biosynthesis”, “carbon fixation in photosynthetic organization”, “response to oxidation” and other pathways, which were not enriched in 93–11 and showed no difference in GX24S, indicating that they play a certain role in salt tolerance (Figs. 3B-C, 4B).

Fig. 4
figure 4

Ubiquitination is the main mechanism underlying the salt tolerance of CY1000. A Venn diagram of DEGs among different rice varieties. B-C The GO and KEGG functional annotation of 371 DEGs specifically expressed in CY1000 and R900. D. Photos of rice seedlings under salt stress after MG-132 treatment

After clustering the upregulated and downregulated genes, it was found that for salt tolerance, the upregulated genes of CY1000 and R900 were involved in oxidative stress, soluble component synthesis, and cell wall synthesis, while the downregulated genes were involved in the acquisition and consumption of most growth substances (Tables S2 and S3). The above results indicate that oxidative stress, cell wall synthesis, soluble sugars and other metabolic pathways also play a role in the salt stress resistance of CY1000 and R900, which is consistent with other research results [18,19,20,21]. Then, physiological indices, including CAT, POD and soluble sugar and protein, were detected (Fig. 3E). According to the results, after salt treatment, all four phenotypes showed increases in these activities, but CY1000 and R900 had higher rates of increase than 93–11 and GX24S, suggesting that these two salt-tolerant phenotypes had more sensitive response mechanisms with regard to these activities.

Interestingly, GO and KEGG functional enrichment analyses showed that CY1000 was most similar to 93–11, followed by R900, and most of the functional modules that differed were the same in the two comparisons, but there was only one overlap with GX24S, indicating that the salt resistance of CY1000 was inherited from R900 but had some “atavism”, which may be part of the reason for heterosis (Fig. 3B).

Ion balance and ubiquitination are key factors affecting the salt tolerance of CY1000

To further reveal the unique salt resistance mechanism of CY1000 and identify the most critical regulatory nodes, we analyzed the four varieties by Venn diagram analysis. A total of 371 genes were specific to CY1000 and R900, which means that these genes are the most critical genes for salt resistance (Fig. 4A, Table S6). Through enrichment analysis of these 371 genes, it was found that these genes were still enriched in pathways associated with ion concentration regulation (ion response and transport, mainly MFS domain-containing proteins), redox activity (flavone and glutathione metabolism), soluble substances and so on.

Among the downregulated genes, those encoding MFS domain-containing proteins, protein kinase domain-containing proteins and peroxidases were the most important. Among the upregulated genes, those encoding ubiquitin-like domain-containing proteins (Ubxs) and F-box domain-containing proteins, which were all previously unidentified, accounted for the majority; this indicates that enhanced ubiquitination is another reason for the strong salt resistance of CY1000 and R900 (Table 1, Fig. 5, Table S6).

Table 1 Salt-responsive ubiquitination genes in CY1000 and R900
Fig. 5
figure 5

qRT-PCR of selected Ubiquitination genes in treated four phenotypes. Each samples has three bio-replicates, and it was tested by student’s t-test. * represents p-value < 0.05 and ** represents p-value< 0.01

Three Ubx proteins and one F-box domain-containing protein were significantly upregulated (143.9-, 114.3-, 16.6- and 56.1-fold, respectively) after salt stress according to transcriptome data. Additionally, in the qRT–PCR results, these four genes showed a high increase in CY1000, as they had a significant increase in R900. In addition, six ring-type domain-containing proteins were also significantly induced by salt stress (three of them were upregulated more than 5-fold, and the other three were downregulated more than 2-fold), indicating that ubiquitination played an important role in the salt tolerance mechanism of CY1000. Interestingly, according to qRT–PCR, two ring-type domain-containing protein genes, BGIOSGA020610 and BGIOSGA011087, did not show the same decreasing pattern in R900 as in CY1000, so the two genes need to be further functionally studied to reveal their significance in the salt tolerance response.

Then, after the seedlings of CY1000 and its parents were treated with the ubiquitin inhibitor MG-132, the salt tolerance of CY1000 and R900 decreased significantly, and there was almost no difference between the phenotypes of CY1000 and R900 after salt stress treatment (Fig. 4D&E). These results indicate that ubiquitination-mediated regulation is the key pathway of salt tolerance in CY1000.

Role of transcription factors in the salt tolerance of CY1000

Transcription factors are the main regulatory factors. Among 371 salt resistance genes, 24 transcription factors were found, including 3 ERFs, 3 NACs, 5 bHLHs, and 2 WRKYs (Table S8). Four were upregulated, while 20 were downregulated, indicating that CY1000 and R900 inhibit the acquisition and metabolism of growth substances by downregulating the expression of these transcription factors. Among them, BGIOSGA012597 (ERF) and BGIOSGA018654 (NAC) had higher expression abundance in CY1000, and the difference factor after salt treatment was similar to that observed for R900, but the abundance increased by more than 2.5–3-fold. This finding indicated that this ERF may have a more active regulatory role in the hybrid CY1000, which may explain why CY1000 is more salt resistant. The abundances of BGIOSGA0209257 (NAC), BGIOSGA02290 (B3), BGIOSGA23706 (WRKY) and BGIOSGA02706 (NAC) were low, and these transcription factors were more strongly downregulated after salt treatment. Except for BGIOSGA02090 (B3), the expression abundance of these transcription factors in GX24S was similar to that in R900, indicating that the local expression level and salt stress response of these genes developed after the formation of hybrids.

Analysis of 1121 differentially expressed genes specific to CY1000 revealed heterosis

We next analyzed the specific salt stress response genes of CY1000 and found that they were mainly enriched in ribosome, C/N metabolism, chlorophyll synthesis, etc., wherein ribosome and other pathways were significantly upregulated, while chlorophyll and other pathways were significantly downregulated. These results indicate that CY1000 not only inherited the salt tolerance mechanisms of R900 but also developed further enhanced salt tolerance characteristics, that is, it underwent heterosis by regulating these mechanisms.

To explore the mechanism underlying the heterosis of CY1000, we studied the transcription factors among the 1121 genes (Table S7). There were 68 transcription factors, including 10 bHLHs, 6 ERFs, 7 NACs, 5 bZIPs and 5 WRKYs. Thirty-two of them were downregulated, and 36 were upregulated. Among them, 6 NACs, 6 bHLHs, 4 ERFs and 4 WRKYs were downregulated, and 4 bZIPs were upregulated. In particular, bHLH BGIOSGA027780 was upregulated 65.5-fold under salt treatment in CY1000. Additionally, three WRKYs (BGIOSGA014175, BGIOSGA036756 and BGIOSGA034613) were downregulated at least 20-fold. One NF-YA (BGIOSGA005501) and two bHLH (BGIOSGA000097 and BGIOSGA006236) were downregulated approximately 10-fold. All these specific genes in hybrids need to be further studied to clarify their relationships with salt tolerance and heterosis.

Discussion

The salt resistance of CY1000 hybrids was mainly inherited from R900, and oxidative stress, balance of ion concentration and increase in soluble matter were the reasons for salt resistance

It has been widely reported that the basic model of rice salt resistance involves oxidative stress, ion concentration balance, an increase in soluble matter increases and so on [10,11,12]. In our results and reported studies, the 93–11 line also has such a response mechanism, but it did not show the corresponding salt resistance, which indicates that in addition to these, it has a stronger salt resistance mechanism (Fig. 3) [22]. In addition, R900 showed the same pattern, and GX24S did not respond to salt stress. This also indicated that the salt resistance of CY1000 was inherited from R900, but CY1000 showed stronger salt resistance (Fig. 2B), reflecting a certain degree of heterosis. However, GX24S had no salt tolerance, which indicated that the heterosis was not caused by superposition. Through multiomic comparative analysis, 371 genes were found to be specific to CY1000 and R900, which indicated that these genes were the main functional genes for salt resistance in the two varieties. These genes were mainly involved in the balance of ion concentration (ion transport, preservation, etc.) and oxidative stress, and included, for example, genes encoding MFS proteins and those involved in the synthesis of antioxidant substances, indicating that these components constitute the main pathways of salt resistance (Fig. 4B&C, Table S6).

Increased ubiquitination is another important reason for salt resistance

Moreover, we found that Ubx proteins (sharing 100% similarity with ubiquitin; accession numbers a2z1x7, b8bcu7 and a2xeq8) and many F-box domain-containing proteins, including 10 unidentified proteins, were significantly different among the 371 genes, which shows that ubiquitination plays a key role in this process. Previous studies have shown that ubiquitination plays an important role in salt tolerance in plants [23,24,25,26,27]. For example, expression of the ospub15 gene of rice (Oryza sativa) was induced under salt stress, and the root and stem development of the ospub15 mutant was slow [28]. Moreover, the RNA interference (RNAi)-generated ospub15 mutant plant was extremely sensitive to salt stress. Further studies have shown that ospub15 can reduce reactive oxygen species (ROS) bursts and cell mortality, thus positively regulating the salt stress response [28]. For example, in A. thaliana, the E2-regulated enzymes ubiquitin carrier proteins (UBCs) 1 and 2 regulate the ubiquitination of H2B and then regulate MYB42 by inducing mpk4 expression. This leads to enhanced phosphorylation of MYB42 under salt stress to maintain the stability of MYB42, followed by regulation of the tolerance of plants to salt stress [29].

According to our results, ubiquitination was increased (16-fold for one and more than 114-fold for the other two in CY1000 and more than 8-fold and 62-fold in R900), indicating that ubiquitination plays an important role in the salt resistance of CY1000 and R900. We used MG-132, an inhibitor of ubiquitin degradation, to treat CY1000 and R900 and then observed the changes under salt stress. The salt tolerance of CY1000 and R900 decreased, indicating that ubiquitination played an important role.

In addition, we found that ubiquitination-related genes in CY1000 had a stronger response pattern than those in R900, which may be one of the reasons for the heterosis of CY1000 (Table 1, Table S6).

CY1000 heterosis occurs due to its more active regulatory mode, in which some transcription factors play a role

According to the above results, we know that the salt resistance of CY1000 was basically derived from R900. However, the difference in the salt resistance phenotype should be due to differences in expression levels. Our transcriptome data also showed that the difference in the expression of 371 genes with common salt resistance was greater in CY1000 (Fig. 4, Table S6). Particularly, among the 371 genes, a series of TFs and functional genes were expressed at higher levels in the hybrid CY1000 than in the parental line R900. For example, two novel TFs, BGIOSGA012597 (ERF) and BGIOSGA018654 (NAC), were 2.5 to 3 times abundant in the hybrid. BGIOSGA018654 was most similar to AtNAC075 (also called SND4), which is considered involved in the negative regulation of flowering and secondary cell wall biosynthesis [30]. BGIOSGA012597 (ERF) was closest to Arabidopsis RAP2–3, which is comprehensively involved in stress responses, such as heat, H2O2 and osmotic stresses [31]. According to these results, it can be concluded that these TFs might be valuable for alleviating salt stress, and CY1000 showed better salt tolerance than R900. The conditions for the high expression of these TFs and genes in heteroalacles need to be further clarified to explain heterosis.

Moreover, there were 1121 DEGs specifically expressed in CY1000 but not in R900 (Fig. 4A, Table S7). More responsive genes are a main cause of heterosis, which has already been confirmed in many reports, revealing a more active regulatory mode [22].

Transcription factors are the main regulators in plants, probably the main reason for the active regulation mode in CY1000. In our data, we found that 24 transcription factors may play an important role in the salt tolerance response, and 5 of them showed a relatively strong response, which may play a role in heterosis (Table S8). However, most of the transcription factors did not show obvious expression and response advantages. However, 68 of the 1121 differentially expressed genes specific to CY1000 were not differentially expressed in the parents, indicating that this transcription regulatory network did play an important role in the more active regulation of heterosis. Most of these transcription factor families, including the bHLH, WRKY and NAC families, are involved in salt stress resistance, and some of them have been proven to be salt stress-related factors, such as BGIOSGA011061, which is most similar to DREB26 in Arabidopsis thaliana and is related to salt stress [32]. BGIOSGA036921 (NAC domain-containing protein 77, onac132, onac300) also responded to salt stress [33]. Some transcription factors only play a role in the response of CY1000 to salt stress, indicating that the regulation of these transcription factors may be the mechanism of heterosis. Moreover, not only stress-responsive TFs but also many developmental TFs were expressed specifically and extremely differentially in salt-treated CY1000. For instance, BGIOSGA027780 (bHLH), which is similar to HEC3, which is responsible for transmitting tract and stigma development in Arabidopsis, was upregulated 65-fold [34]. These results indicated that heterosis was caused by the change in the degree of regulation of some genes after crossing between the two parents, resulting in a more sensitive response to stress.

Conclusions

Soil salinization is a threat to food security, and China is rich in saline land resources for potential and current utilization. The super hybrid rice variety CY1000 is one of the best salt-tolerant rice varieties and is widely used. Therefore, the analysis of the molecular mechanism of salt tolerance of CY1000 can provide a theoretical basis for breeding new salt-tolerant rice varieties. Therefore, an analysis of a comparative transcriptomic dataset of seedlings of CY1000 and its parental lines R900 and GX24S along with 93–11 under salinity stress and controlled conditions was conducted. CY1000 inherited its salt tolerance from the R900 parent. A total of 371 DEGs were found to be specific to only the salt tolerant phenotypes CY1000 and R900, suggesting that these genes are associated with the salt response. Among these genes, MFS, especially ubiquitination-related genes, exhibited the most responsive patterns. MG-132 pretreatment repressed the salt tolerance of CY1000 and R900, suggesting that increased ubiquitination is important for salt tolerance. Additionally, some candidate ubiquitination-related genes and TFs were found and considered important for the salt response. Additionally, CY1000 also showed a better salt tolerance than the parent R900, indicating that heterosis exists in CY1000. Among the 371 genes, a series of TFs and genes were expressed at higher levels in CY1000 than in R900, supporting the occurrence of heterosis. Moreover, a total of 1121 DEGs specific to CY1000 were found, which could be the other reason for heterosis. Concludingly, our results provide a picture of how CY1000 responds to salt stress and heterosis.

Materials and methods

Plant materials, seed germination and salt stress treatment

CY1000 is a super hybrid rice variety cultivated by the Hunan Hybrid Rice Research Center and approved by the state in 2017, and its parents are GX24S and R900; the rice line 93–11 was used as a control. Variety information can be queried at the China Rice Data Center (https://www.ricedata.cn/variety/varis/616642.htm). All rice seeds came from long-term reproduction in our laboratory, and there are no intellectual property disputes related to these materials. The seeds used in this study were harvested in the autumn of 2020, and batches of fresh seeds were placed on wet filter paper for germination, followed by transplantation into seedling pots for hydroponic seedling cultivation for 7 d. Then, some of the seedlings were hydroponically grown in 0.6% (100 mM) NaCl solution, and some of the seedlings were transferred to fresh water solution for culture. After 7 d of growth, the phenotypes were observed, and whole plants were taken for omics sequencing. Each phenotype had three replicates for RNA-seq.

Total RNA extraction, mRNA purification and cDNA library construction

Total RNA was extracted by TRIzol [35], the gel was removed, the absorbance value was measured, and the quality was detected by an Agilent 2100 Bioanalyzer.

The mRNA containing the poly-A structure was enriched by magnetic beads, and the fragment length was 200–300 bp. Then, the cDNA library (380 bp) was constructed by using 6-base random primers.

High-throughput RNA sequencing and analysis

Based on the Illumina sequencing platform, these libraries were sequenced by paired-end (PE) (2 * 150 bp) sequencing.

First, the raw data were filtered, and cutadapt [36] was used to remove the sequences with 3′ end connectors. Reads with average mass fractions lower than Q20 were removed. The filtered clean data were aligned with HISAT2 (http://ccb.jhu.edu/software/hisat2/index.shtml) [37]. This software compares filtered reads with a reference genome. HISAT2 uses an improved BWT algorithm. The reference genome in Ensembl was Oryza_ indica.ASM465v1.dna.toplevel.fa.

We used the union scheme of htseq to compare the read count value of each gene, which was used as the original gene expression level. The counts of reads were positively correlated with the real expression levels of genes and the length and sequencing depth of genes. To make the gene expression levels of different genes and different samples comparable, we used fragments per kilobase per million fragments mapped (FPKM) values to standardize the expression levels. FPKM values indicate the number of fragments from a gene per thousand base length per million fragments. For PE sequencing, each fragment has two reads, and the FPKM value is calculated based on only the number of fragments of two reads that can be compared to the same transcript. In the transcriptome, we generally considered genes with FPKM > 1 to be expressed.

Identification of differentially expressed genes

We used DESeq [38] to analyze the differences in gene expression and used the conditions for differential gene expression as follows: expression difference multiple | log2fold change | > 1, p value < 0.05.

Enrichment analysis

topGO was used for GO enrichment analysis. During the analysis, the gene list and gene number of each term were calculated by using the different genes annotated by GO terms. Then, the p value (p value < 0.05) was calculated by the hypergeometric distribution method to determine the GO terms with significant enrichment of different genes compared with the whole-genome background to determine the main biological functions of different genes.

MG-132 treatment

After germination and seedling development, the seeds were treated with 10 μM MG-132 for 6 h before salt treatment and then transferred to salt water for 7 d. The phenotype was observed, and the survival rate was analyzed by ANOVA.

Physiological index detection and qRT–PCR

Four phenotypes, 93–11, GX64S, R900 and CY1000, were treated with 0.6% NaCl solution, and the whole plant was sampled to detect the CAT and POD activities and the contents of soluble sugar and soluble protein. Additionally, the expression levels of the ten ubiquitination-related genes were tested by qRT–PCR. Actin was used as the reference gene, and all the primers are listed in Table S9. All experiments were conducted with three replicates.

Availability of data and materials

The raw reads produced in this study were deposited in the NCBI SRA under bioProject number PRJNA785275. The clean data were submitted to SRA under SRP349079.

References

  1. Butcher K, Wick AF, DeSutter T, Chatterjee A, Harmon J. Soil salinity: a threat to global food security. Agron J. 2016;108(6):2189–200.

    Article  CAS  Google Scholar 

  2. Wang W-J, He H-S, Zu Y-G, Guan Y, Liu Z-G, Zhang Z-H, et al. Addition of HPMA affects seed germination, plant growth and properties of heavy saline-alkali soil in northeastern China: comparison with other agents and determination of the mechanism. Plant Soil. 2010;339(1–2):177–91.

    Google Scholar 

  3. Qin H, Li Y, Huang R. Advances and challenges in the breeding of salt-tolerant Rice. Int J Mol Sci. 2020;21(21):8385.

    Article  CAS  Google Scholar 

  4. Huang Z, Lv Q, Xin Y, Yuan L, Fu X, Zhu L, et al. Heterotic performance of the main yield traits in different types of Indica hybrid rice. Food Energy Secur. 2020;9(3):e210.

    Google Scholar 

  5. Wang H, Zhang M, Guo R, Shi D, Liu B, Lin X, et al. Effects of salt stress on ion balance and nitrogen metabolism of old and young leaves in rice (Oryza sativa L.). BMC Plant Biol. 2012;12:194.

    Article  Google Scholar 

  6. Zhang Y, Fang J, Wu X, Dong L. Na(+)/K(+) Balance and transport regulatory mechanisms in weedy and cultivated Rice (Oryza sativa L.) under salt stress. BMC Plant Biol. 2018;18(1):375.

    Article  CAS  Google Scholar 

  7. Huang SC, Xin SC, Xie GQ, Han J, Liu ZL, Wang B, et al. Mutagenesis reveals that the rice OsMPT3 gene is an important osmotic regulatory factor. Crop J. 2020;8(3):465–79.

    Article  Google Scholar 

  8. Nounjan N, Chansongkrow P, Charoensawan V, Siangliw JL, Toojinda T, Chadchawan S, et al. High performance of photosynthesis and osmotic adjustment are associated with salt tolerance ability in Rice carrying drought tolerance QTL: physiological and co-expression network analysis. Front Plant Sci. 2018;9:1135.

    Article  Google Scholar 

  9. Vaidyanathan H, Sivakumar P, Chakrabarty R, Thomas G. Scavenging of reactive oxygen species in NaCl-stressed rice (Oryza sativa L.)—differential response in salt-tolerant and sensitive varieties. Plant Sci. 2003;165(6):1411–8.

    Article  CAS  Google Scholar 

  10. Liu C, Mao B, Yuan D, Chu C, Duan M. Salt tolerance in rice: physiological responses and molecular mechanisms. Crop J. 2021;10:3–25.

  11. Kumar K, Kumar M, Kim SR, Ryu H, Cho YG. Insights into genomics of salt stress response in rice. Rice. 2013;6:1–5.

    Article  Google Scholar 

  12. van Zelm E, Zhang Y, Testerink C. Salt tolerance mechanisms of plants. Annu Rev Plant Biol. 2020;71:403–33.

    Article  Google Scholar 

  13. Wu GQ, Wang SM. Calcium regulates K+/Na+ homeostasis in rice (Oryza sativa L.) under saline conditions. Plant Soil Environ. 2012;58(3):121–7.

    Article  CAS  Google Scholar 

  14. Wu SJ, Ding L, Zhu JK. SOS1, a genetic locus essential for salt tolerance and potassium acquisition. Plant Cell. 1996;8(4):617–27.

    Article  CAS  Google Scholar 

  15. Soni P, Kumar G, Soda N, Singla-Pareek SL, Pareek A. Salt overly sensitive pathway members are influenced by diurnal rhythm in rice. Plant Signal Behav. 2013;8(7):e24738.

    Article  Google Scholar 

  16. Zhang J, Jia W, Yang J, Ismail AM. Role of ABA in integrating plant responses to drought and salt stresses. Field Crop Res. 2006;97(1):111–9.

    Article  Google Scholar 

  17. Cao WH, Liu J, He XJ, Mu RL, Zhou HL, Chen SY, et al. Modulation of ethylene responses affects plant salt-stress responses. Plant Physiol. 2007;143(2):707–19.

    Article  CAS  Google Scholar 

  18. Couee I, Sulmon C, Gouesbet G, El Amrani A. Involvement of soluble sugars in reactive oxygen species balance and responses to oxidative stress in plants. J Exp Bot. 2006;57(3):449–59.

    Article  CAS  Google Scholar 

  19. Chaves MM, Flexas J, Pinheiro C. Photosynthesis under drought and salt stress: regulation mechanisms from whole plant to cell. Ann Bot. 2009;103(4):551–60.

    Article  CAS  Google Scholar 

  20. Parvaiz A, Satyawati S. Salt stress and phyto-biochemical responses of plants – a review. Plant Soil Environ. 2008;54(3):89–99.

    Article  CAS  Google Scholar 

  21. Oliveira DM, Mota TR, Salatta FV, Sinzker RC, Koncitikova R, Kopecny D, et al. Cell wall remodeling under salt stress: insights into changes in polysaccharides, feruloylation, lignification, and phenolic metabolism in maize. Plant Cell Environ. 2020;43(9):2172–91.

    Article  CAS  Google Scholar 

  22. Jahan N, Lv Y, Song M, Zhang Y, Shang L, Lu Y, et al. Transcriptomic analysis of short-term salt-stress response in mega hybrid Rice seedlings. Agronomy. 2021;11(7):1328.

    Article  CAS  Google Scholar 

  23. Zhang Z, Li J, Liu H, Chong K, Xu Y. Roles of ubiquitination-mediated protein degradation in plant responses to abiotic stresses. Environ Exp Bot. 2015;114:92–103.

    Article  CAS  Google Scholar 

  24. Lou L, Yu F, Tian M, Liu G, Wu Y, Wu Y, et al. ESCRT-I component VPS23A sustains salt tolerance by strengthening the SOS module in Arabidopsis. Mol Plant. 2020;13(8):1134–48.

    Article  CAS  Google Scholar 

  25. Wang W, Wang W, Wu Y, Li Q, Zhang G, Shi R, et al. The involvement of wheat U-box E3 ubiquitin ligase TaPUB1 in salt stress tolerance. J Integr Plant Biol. 2020;62(5):631–51.

    Article  CAS  Google Scholar 

  26. Tian M, Lou L, Liu L, Yu F, Zhao Q, Zhang H, et al. The RING finger E3 ligase STRF1 is involved in membrane trafficking and modulates salt-stress response in Arabidopsis thaliana. Plant J. 2015;82(1):81–92.

    Article  Google Scholar 

  27. Wang J, Qin H, Zhou S, Wei P, Zhang H, Zhou Y, et al. The ubiquitin-binding protein OsDSK2a mediates seedling growth and salt responses by regulating gibberellin metabolism in Rice. Plant Cell. 2020;32(2):414–28.

    Article  CAS  Google Scholar 

  28. Park JJ, Yi J, Yoon J, Cho LH, Ping J, Jeong HJ, et al. OsPUB15, an E3 ubiquitin ligase, functions to reduce cellular oxidative stress during seedling establishment. Plant J. 2011;65(2):194–205.

    Article  CAS  Google Scholar 

  29. Sun Y, Zhao J, Li X, Li Y. E2 conjugases UBC1 and UBC2 regulate MYB42-mediated SOS pathway in response to salt stress in Arabidopsis. New Phytol. 2020;227(2):455–72.

    Article  CAS  Google Scholar 

  30. Fujiwara S, Mitsuda N. ANAC075, a putative regulator of VASCULAR-RELATED NAC-DOMAIN7, is a repressor of flowering. Plant Biotechnol (Tokyo). 2016;33(4):255–65.

    Article  CAS  Google Scholar 

  31. Papdi C, Perez-Salamo I, Joseph MP, Giuntoli B, Bogre L, Koncz C, et al. The low oxygen, oxidative and osmotic stress responses synergistically act through the ethylene response factor VII genes RAP2.12, RAP2.2 and RAP2.3. Plant J. 2015;82(5):772–84.

    Article  CAS  Google Scholar 

  32. Kazama D, Itakura M, Kurusu T, Mitsuda N, Ohme-Takagi M, Tada Y. Identification of chimeric repressors that confer salt and osmotic stress tolerance in Arabidopsis. Plants (Basel). 2013;2(4):769–85.

    Article  Google Scholar 

  33. Fang Y, You J, Xie K, Xie W, Xiong L. Systematic sequence analysis and identification of tissue-specific or stress-responsive genes of NAC transcription factor family in rice. Mol Gen Genomics. 2008;280(6):547–63.

    Article  CAS  Google Scholar 

  34. Ballester P, Martinez-Godoy MA, Ezquerro M, Navarrete-Gomez M, Trigueros M, Rodriguez-Concepcion M, et al. A transcriptional complex of NGATHA and bHLH transcription factors directs stigma development in Arabidopsis. Plant Cell. 2021;33(12):3645–57.

    Article  Google Scholar 

  35. Zhu MD, Zhang M, Gao DJ, Zhou K, Tang SJ, Zhou B, et al. Rice OsHSFA3 gene improves drought tolerance by modulating polyamine biosynthesis depending on Abscisic acid and ROS levels. Int J Mol Sci. 2020;21(5):1857.

    Article  CAS  Google Scholar 

  36. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17(1):10–2.

    Google Scholar 

  37. Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37(8):907–15.

    Article  CAS  Google Scholar 

  38. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

This work was supported by the 2020 Research Program of Sanya Yazhou Bay Science and Technology City (SKJC-2020-02-006) and the Hainan Major Science and Technology Projects (ZDKJ202001). The funding bodies did not play a role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript but only provided financial support.

Funding

Supported by the 2020 Research Program of Sanya Yazhou Bay Science and Technology City (SKJC-2020-02-006) and the Hainan Major Science and Technology Projects (ZDKJ202001).

Author information

Authors and Affiliations

Authors

Contributions

Data curation, Z. M and Z. B; Funding acquisition, G. X-Y and A. Z-Y; Methodology, Z. M-D; Project administration, D. H-F; Resources, L. J-R and A. Z-Y; Software, Z. M; Validation, L. J-W and W. Z-W; Writing – original draft, G. X-Y and Z. M; Writing – review & editing, D.H-F. All authors reviewed the manuscript. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Zhou Bin, Ai Zhi-Yong or Deng Hua-Feng.

Ethics declarations

Consent for publication

All authors agree to the publication of this manuscript.

Competing interests

The authors declare that they have no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Information on transcriptomic data. A. PCA of all samples. B-E. Gene expression levels of CY1000, R900, GX24S and 93–11. F. Correlation analysis of all samples.

Additional file 2: Table S1.

Expression levels of all identified genes in all samples. C, R, G and Y are short for CY1000, R900, GX24S and 93–11.

Additional file 3: Table S2.

DEGs of CY1000 under salt stress treatment.

Additional file 4: Table S3.

DEGs of R900 under salt stress treatment.

Additional file 5: Table S4.

T DEGs of GX24S under salt stress treatment.

Additional file 6: Table S5.

T DEGs of 93–11 under salt stress treatment.

Additional file 7: Table S6.

The 371 DEGs specifically associated with the salt response.

Additional file 8: Table S7.

The 1121 DEGs specific to CY1000.

Additional file 9: Table S8.

24 TFs specifically responded to salt stress in CY1000 and R900.

Additional file 10: Table S9.

Primers for the qRT-PCR.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xia-Yu, G., Meng, Z., Ming-Dong, Z. et al. Comparative transcriptomic analysis of the super hybrid rice Chaoyouqianhao under salt stress. BMC Plant Biol 22, 233 (2022). https://doi.org/10.1186/s12870-022-03586-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-022-03586-w

Keywords