Distribution of ST genes in 550 Indica/Xian rice accessions
To identify the distribution of ST genes in the rice population, as shown in Fig. 1, we selected nine ST-related SNPs in seven reported genes, including STG5I12S, SKC1P140A, SKC1R184H, OsHKT2;3I77T, OsHKT1;1L94K, OsWRKY53A173G, OsSTL1P289S, RST1A530G and RST1E611G, which have been identified as key sites conferring ST in rice (Ren et al. 2005; Mishra et al. 2016; Wang et al. 2018; Yuan et al. 2020; Deng et al. 2022; Lv et al. 2022; Yu et al. 2023; Wei et al. 2024). We subsequently investigated the proportion of allelic SNPs (tolerance/sensitive, T/S) in a resequencing population of 550 Indica/Xian rice accessions (Meng et al. 2021). The findings indicated that six T-SNPs were highly conserved in this population, including SKC1140A (79.27%), RST1611G (93.82%), OsWRKY53173G (92.55%), STG512S (92.91%), OsSTL1289S (100%), and OsHKT2;377T (88.55%). The remaining three elite alleles, SKC1184H (31.63%), RST1530G (6.36%), and OsHKT1;194K (19.82%), accounted for relatively rare proportions (Table 1). Furthermore, we divided all 550 accessions into twenty-one haplotypes on the basis of the T/S types of the nine SNPs, among which Hap1 and Hap2 were two major haplotypes consisting of 202 and 105 cultivars containing 6 and 7 ST-related SNPs, respectively (Fig. 1). Interestingly, rice is highly sensitive to salt stress (Yu et al. 2023; Wei et al. 2024); however, each Indica/Xian rice accession was found to carry at least four out of seven detected ST genes, and cultivars of Hap3 and Hap5 carried all seven reported ST genes (Fig. 1).
Table 1
Distribution of nine SNPs in the 550 Indica/Xian accessions.
SNP alleles | Base Type (S/T)1 | Homozygous(S) | Heterozygous | Homozygous(T) | Percentages of T (%) |
STG5I12S | A/C | 30 | 9 | 511 | 92.91 |
SKC1P140A | C/G | 95 | 19 | 436 | 79.27 |
SKC1R184H | G/A | 358 | 18 | 174 | 31.63 |
OsHKT1;1L94K | G/A | 428 | 13 | 109 | 19.82 |
OsHKT2;3I77T | T/C | 56 | 7 | 487 | 88.55 |
OsWRKY53A173G | C/G | 37 | 4 | 509 | 92.55 |
OsSTL1P298S | G/A | 0 | 0 | 550 | 100 |
RST1A530G | C/G | 504 | 11 | 35 | 6.36 |
RST1E611G | A/G | 32 | 2 | 516 | 93.82 |
1 S/T represents the salt-sensitive and salt-tolerant type respectively. |
Assessment of the functions of three OsHKT family genes
SKC1 and OsHKT1;1 are both HKT family genes that have been verified to facilitate the export of Na+ to maintain the cellular ion balance during saltwater stress (Ren et al. 2005; Wang et al. 2015; Takagi et al. 2015; Lv et al. 2022). To evaluate the functions of SKC1140A, SKC1184H and OsHKT1;194K, we compared ST at the seeding stage between two haplotypes whose T-SNPs were identical except for the target haplotype. Hap1-Hap6 commonly contained five T-SNPs, namely, STG512S, OsHKT2;377T, OsWRKY53173G, OsSTL1289S and RST1611G. Compared with Hap4 rice accessions, which carry only the above five common T-SNPs, Hap1 and Hap6 rice accessions carry additional T-SNPs of SKC1140A and OsHKT1;194K, respectively. Compared with Hap1 rice accessions, Hap2 and Hap5 rice accessions carry one more T-SNP of SKC1184H or OsHKT1;194K, respectively, and Hap3 rice accessions carry both T-SNPs of SKC1184H and OsHKT1;194K (Fig. 1).
Eight plants selected from six random cultivars of each haplotype were subjected to ST comparison treatment (sup Table 1). After 150 mM NaCl treatment for 7 d and rehydration for 2 d, obvious damage was observed in the tested groups compared with the control group (Fig. 2A, B). The relative shoot fresh weight (rSFW) was subsequently used as the ST index. Compared with that of the Hap4 plants, the increased rSFW was approximately 23.8% greater for those from Hap1, representing the ST contribution of SKC1140A. Consistent with the results of the comparison between the Hap1 and Hap4 groups, the rSFW of Hap5 was approximately 27.5% greater than that of Hap6 (Fig. 2C). These results suggested that the haplotype SKC1140A could increase the ST in all the tested genetic backgrounds. Similar to the results of the analysis of rSFW for the other groups, Hap2 and Hap3 improved by approximately 18.9% and 19.0%, respectively, over those of Hap1 and Hap5 (Fig. 2C). On the basis of the previously predicted structure of SKC1 (Ren et al. 2005), the amino acid variations of SKC1140A and SKC1184H were both located in a cytoplasm-exposed loop (sup Fig. 1), indicating that SKC1184H likely contributed to ST through synergistic effects with SKC1140A. However, no significant differences were detected between the two comparison haplotypes of groups Hap6 and Hap4, Hap5 and Hap1, or Hap3 and Hap2, indicating that the rice varieties that carried OsHKT1;194K did not contribute to increased ST under these treatment conditions (Fig. 2C). Because OsHKT1;194K has been shown to enhance ST in several previous reports (Wang et al. 2015; Alfatih et al. 2023; Wei et al. 2024), an extended treatment was performed between Hap2 and Hap3 (Fig. 2D). After 14 d of treatment, a 37.9% improvement in rSFW and a 17.4% decrease in the dead leaf rate (DLR) were detected in Hap3 compared with those in Hap2 (Fig. 2E, F).
Several SNPs located in OsHKT2;3 were revealed to be associated with ST traits via GWAS, among which OsHKT2;3I77T was detected twice in different studies (Mishra et al. 2016; Yu et al. 2023). Therefore, the ST of OsHKT2;377T might be worthy of evaluation. According to the haplotype groups, Hap3 and Hap9 were found to differ only in OsHKT2;377T (Fig. 1). Six cultivars of Hap3 (OsHKT2;377T) and five cultivars of Hap9 (OsHKT2;377I) were selected for rSFW comparison (sup Table 1). Compared with the Hap3 plants, the Hap9 plants significantly improved after 150 mM NaCl treatment, suggesting the importance of OsHKT2;377T for the maintenance of ST in Indica/Xian cultivars (Fig. 3A, B).
Functional evaluation of RST1
RST1 is a major ST gene that encodes a negative regulator of OsARF18 and affects NH4+ metabolism. Evaluation was performed with RST1611G and RST1530G, two key T-SNPs detected via GWAS (Deng et al. 2022). Our haplotype analysis revealed that Hap2, Hap1, and Hap9 were distinguished from Hap15, Hap8 and Hap10 only at RST1611G, respectively (Fig. 1), and Hap12 was distinguished from Hap11 at RST1530G (Fig. 1). Thus, a comparison analysis of ST was performed for each of the four pairs of groups. As shown in Fig. 4, significant ST was observed on Hap2, Hap1, and Hap9 plants compared with that of Hap15, Hap8 and Hap10 plants (Fig. 4A, C, E), with 88.6%, 34.2% and 41.1% improvements in rSFW (sup Table 1), respectively, after 150 mM NaCl treatment for 7 d and rehydration for another 2 d (Fig. 4B, D, F). To our surprise, Hap12 plants that carried the T-SNP RST1530G were more sensitive to salt treatment than Hap11 plants that carried RST1530A (Fig. 4G). The statistical analysis of rSFW also revealed that the Hap12 haplotype was lower than the Hap11 haplotype (Fig. 4H; sup Fig. 2; sup Table 1). These findings suggest that the ST of Hap11 plants is associated with RST1530A rather than RST1530G. Alternatively, we identified 504 and 516 varieties from 550 accessions that carried RST1530A and RST1611G, respectively (Fig. 1). A total of 478 out of 504 accessions of RST1530A simultaneously carried the T-SNP of RST1611G, suggesting that they collaboratively contributed to ST and were tightly linked during the development of ST.
Evaluation of salt tolerance of OsWRKY53173G
The natural variation of OsWRKY53173G confers ST, which was shown to repress the expression of SKC1 and OsMPKK10.2 in tolerant rice varieties (Yu et al. 2023). For the functional evaluation of OsWRKY53173G, salt treatment was performed as described above. Three groups of haplotypes contained the elite allele of OsWRKY53173G, including Hap2, Hap1 and Hap4, each comprising 3 rice varieties, which were separately compared with Hap18, Hap21 and Hap17, which contained the other identical ST genes and the sensitive allele of OsWRKY53173G (Fig. 5A-F, sup Table 1). After 7 d of treatment with 150 mM NaCl followed by rehydration for 2 d, all the varieties from Hap2 (Fig. 5A) and Hap1 (Fig. 5C) were more tolerant to salt stress than were the control varieties from the Hap18 and Hap21 groups, with 38.6% and 33.7% improvements in the rSFW in each group compared with those in the control group, respectively (Fig. 5B, D). However, there was no obvious difference for all the varieties in the comparison of Hap4 and Hap17 (Fig. 5E, F; sup Table 1). Comparing the ST genes in all three groups, in addition to OsWRKY53173G, we found that additional SKC1140A and/or SKC1184H genes were present in the comparison groups Hap2 to Hap18 and Hap1 to Hap21 but not in Hap4 to Hap17. Therefore, OsWRKY53173G may play a role in ST that is genetically dependent on SKC1140A and/or SKC1184H.
OsWRKY53 negatively regulates ST and functions as a trans-repressor of SKC1 (Yu et al. 2023). The natural variation of OsWRK53173G increased ST and was genetically dependent on SKC1140A and/or SKC1184H (Yu et al., 2023; Fig. 5A-F). However, whether natural variation in the ST alleles of OsWRK53173G can repress the expression of SKC1 remains unclear. Here, we selected two rice varieties from each of the six haplotypes shown in Fig. 5A-F and performed a quantitative analysis of the relative expression of OsWRKY53 and SKC1 in roots without or with 150 mM NaCl for 8 h. As shown in Fig. 5G, 10 out of the 12 rice varieties that contained OsWRK53173G or not presented upregulated expression of SKC1, and only two varieties, Hap4#-1 and Hap17#2, presented downregulated expression after salt treatment (Fig. 5G). Moreover, 8, 1 and 3 of the 12 rice varieties presented upregulated, downregulated and no significant changes in the expression of OsWRKY53, respectively, after salt treatment (Fig. 5H). Interestingly, the salt-regulated expression of SKC1 was correlated with that of OsWRKY53 (R2 = 0.8020, P < 0.01). Overall, we did not observe the repressed expression of SKC1 caused by either OsWRKY53173G or OsWRKY53173A in natural variations under salt treatment.