Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 01 November 2024
Sec. Plant Abiotic Stress
This article is part of the Research Topic Advanced Breeding for Abiotic Stress Tolerance in Crops, Volume II View all 6 articles

GWAS combined with linkage analysis reveals major QTLs and candidate genes of salt tolerance in Japonica rice seedlings

Shanbin Xu&#x;Shanbin Xu1†Jie Zheng&#x;Jie Zheng1†Haoqiang Du&#x;Haoqiang Du1†Xiaodong DuXiaodong Du2Chong LiChong Li1Yuxuan DuanYuxuan Duan1Yanan CaiYanan Cai1Jingguo WangJingguo Wang1Hualong LiuHualong Liu1Luomiao YangLuomiao Yang1Wei XinWei Xin1Yan JiaYan Jia1Detang Zou*Detang Zou1*Hongliang Zheng*Hongliang Zheng1*
  • 1Key Laboratory of Germplasm Enhancement and Physiology & Ecology of Food Crop in Cold Region, Ministry of Education/College of Agriculture, Northeast Agricultural University, Harbin, China
  • 2Rice Research Institute, Heilongjiang Academy of Agricultural Sciences, Jiamusi, China

Background: Soil salinization is one of the significant factors limiting high crop yields and expansion of arable land, seriously affecting global agricultural production. Rice is an essential food crop throughout the world, and its seedlings are particularly susceptible to salt stress, which can directly affect the growth and development of rice and its final yield. We used the natural population as the material for genome-wide association study (GWAS) and the recombinant inbred line (RIL) population from CD (salt sensitive)/WD20342 (salt tolerant) hybridization as the material for linkage analysis.

Results: The degree of salt tolerance was evaluated by using the relative root length (RRL), relative root number (RRN), relative root fresh weight (RRFW), and relative root dry weight (RRDW) as indicators. Fifteen and six major quantitative trait loci (QTLs) were identified by GWAS and linkage analysis, respectively. Meanwhile, the GWAS identified the lead SNP (Chr2_22340368), which was located within qRRL2 and qRRDW2 identified by linkage analysis. GWAS, combined with linkage analysis, selected a 196-kb overlapping region on chromosome 2, including 22 candidate genes. LOC_Os02g36880 was discovered as the candidate gene involved in salt tolerance by haplotype analysis, qRT-PCR, and sequence analysis. The score of salinity toxicity (SST) and seedling survival rate (SSR) were determined for CRISPR/Cas9 mutants (CR-1 and CR-15) and wild-type (ZH11), respectively.

Conclusion: The phenotypic validation indicated that LOC_Os02g36880 negatively regulated the salt tolerance at the seedling stage. This study provides resources for breeding Japonica rice to improve its response to salt stress.

1 Introduction

Rice is an essential crop globally, and with the growing global population, the demand for higher rice yields is increasingly urgent. However, various abiotic factors including drought, salinity, alkalinity, and cold, significant impact crop yields (Radha et al., 2023; Raj and Nadarajah, 2023; Sarma et al., 2023). Among these, salt damage is one of the most critical factors reducing rice productivity, as rice is highly sensitive to saline conditions (Ruan et al., 2011; Chen et al., 2022). Soil salinity affects a substantial portion of arable land, and the shrinking availability of agricultural land further constrains global agricultural productivity. High salt concentrations lead to ionic toxicity, hyperosmotic stress, and secondary stresses, which together disrupt essential rice processes throughout its life cycle, from germination to senescence (Hanin et al., 2016; Zhu, 2016). Rice is particularly vulnerable to the adverse effects of salt stress during the seedling stages. Thus, understanding the genetic mechanisms and identifying relevant QTLs or genes is crucial for developing salt-tolerant rice cultivars, which forms the theoretical and practical foundation for breeding varieties.

Salt tolerance in rice is a typical quantitative trait jointly controlled by multiple genes and the integrated expression of various mechanisms (Liu et al., 2022; Phosuwan et al., 2024). Significant progress has been made by employing forward genetics strategies, such as GWAS and linkage analysis, to identify numerous QTLs and genes associated with salt tolerance. For instance, Kumar et al. (2015) evaluated twelve salt-related agronomic traits using a natural population of 220 rice accessions. By combining these traits with a gene microarray containing 6,000 single nucleotide polymorphisms (SNPs), they identified sixty-four loci, which explained 5% to 18% of the phenotypic variance. Similarly, Al-Tamimi et al. (2016) conducted a GWAS on 553 rice varieties identified 12 efficient QTLs associated with salt tolerance based on indicators such as relative growth rate, transpiration rate, and transpiration efficiency. Batayeva et al. (2018) identified 26 QTLs linked to salt tolerance in 203 rice accessions, with eleven of these QTLs overlapping with previously known loci. Rohila et al. (2019) used a hydroponic system to subject 162 rice accessions to salt stress during the seedling stage, leading to the identification of nine genomic regions associated with salt tolerance by analyzing over three million SNPs. Additionally, they discovered six new loci and identified sixteen potential genes. Yuan et al. (2020) utilized 664 cultivated rice materials for GWAS, identifying 21 QTLs associated with salt tolerance.

Meanwhile, researchers have also found several QTLs using linkage analysis. De Leon et al. (2016) utilized a RIL population derived from a cross between the Bengal (salt-sensitive) and the Pokkali (salt-tolerant) to identify 85 additive QTLs for nine salt-related traits, including salt injury score, shoot length, root length. Nakhla et al. (2021) identified 23 and 46 QTLs controlling salt tolerance at the germination and seedling stages, respectively, through linkage mapping. Zeng et al. (2021) used the a BC1F2 population derived from a cross between WJZ and Nipponbare to identify the major salt tolerance QTL qGR6.2, during the germination stage and predicted 11 candidate genes. Geng et al. (2023) employed a RIL population derived from a cross between Milyang 23 (salt-tolerant) and Jileng 1 (salt-sensitive) to identify 12 QTLs through linkage mapping.

In this study, we explored the genetic foundation of salt tolerance in rice using GWAS and linkage analysis with RRL, RRN, RRFW, and RRDW during the rice seedling salt treatment. Both GWAS and linkage analysis are powerful tools for detecting QTLs associated with complex traits, and their combination can significantly enhance the accuracy of QTL detection. Our GWAS identified a lead SNP (Chr2_22340368), within the QTLs qRRL2 and qRRDW2, which were also identified through linkage analysis. A 196-kb region on chromosome 2 was selected for further examination based on the overlapping regions identified through both GWAS and linkage mapping. Knockout mutants were subsequently used for the functional verification, leading to the identification of LOC_Os02g36880 as a gene that contributes to salinity resistance.

2 Materials and methods

2.1 Plant materials

The natural population was conducted using the 295 Japonica rice materials gathered from Chinese provinces such as Heilongjiang, Ningxia, and Jilin, along with foreign varieties mainly from Japan, Russia, and Korea (Li et al., 2019). The linkage mapping population consisted of 189 RILs, a cross between the CD (salt-sensitive) and WD20342 (salt-tolerant) (Li et al., 2023). Supplementary Figure S1 shows CD and WD20342 salt tolerance phenotypes. In Heilongjiang Province, China, at the Northeast Agricultural University Experimental Station (45°52′N, 127°04′E), all materials were planted.

2.2 Phenotype data

All materials were baked for 48 hours at 55°C, after choosing and disinfecting the seeds for 20 minutes with a 3% NaClO solution. The seeds were cultured hydroponically at 30°C for 4 days. Using the Yoshida nutrient solution, sixty-four evenly germinated seedlings of each type were divided into salt-treated and control groups for hydroponic cultivation. The germinated seeds were incubated at 26 and 24°C, with a light and dark cycle of 14 and 10 hours, respectively. Once the seedlings reached two leaves and one heart, the control group received standard nutrient solution application; the treatment group received salt stress treatment. Pretreatment of the seedlings was 50 mM NaCl, which was applied for 3 days, and formally treated with 120 mM NaCl for 7 days. Five plants were randomly chosen for the measurement in each material, including the root length, the root number, the root fresh weight, and the root dry weight under control and salt treatment. Further calculations were performed to determine the RRL, RRN, RRFW, and RRDW: the relative values represent the ratio of phenotypic values during salt treatments to phenotypic values under control conditions. Three replications of the experiment were set up.

2.3 GWAS

Using Plink 2.0 software (Chang et al., 2015), 788,396 SNPs were identified with a minor allele frequency (MAF) of ≥5% and a missing rate of ≤20%. The 295 materials were analysed by GWAS using 788,396 SNPs. The calculated population structure (Q) and kinship between individuals were used to adjust the population structure of GWAS (Li et al., 2019). The mixed linear model (MLM) included in TASSEL 5.0 (Bradbury et al., 2007) was used for GWAS. The threshold for the identification of SNPs significantly associated with traits was set to –log10(P) > 5.26, determined by the genetic type 1 error calculator (GEC; http://statgenpro.psychiatry.hku.hk/gec/), which calculates the effective number of independent markers. Generate the Manhattan map using the R program ‘qqman’. To determine the leading SNPs with the lowest P-value, redundant SNPs were filtered using a least-distance interval. Utilizing the LDBlockShow software, the pairwise R2 value between any two SNPs within the interval of leading SNPs ± 2 Mb was determined. The LD attenuation interval of leading SNPs between 1.5 and 2.0 Mb was the average of the top 10% R2 values plus 0.2 (Dong et al., 2020).

2.4 Linkage analysis

Using 10 K Array targeted sequencing technology from MOLBREEDING Biotechnology Company, the genotypes of the RIL lines were determined. Following biparental polymorphism analysis and de-redundant analysis and by the IciMapping-Bin project, there was a grouping of specific SNP markers with the same genotype into one “BIN”. Using 189 RILs from CD and WD20342 and 978 bin markers, a genetic linkage mapping was produced. With an average distance of 2.52 cM between markers, the map contains 2465.33 cm of the genome (Supplementary Figure S2). The inclusive composite interval mapping method used in ICIMapping Ver.4.2 was utilized for linkage analysis, with LOD > 2.5 as the threshold value (Lei et al., 2021).

2.5 Haplotype analysis

Non-synonymous mutant SNPs in the promoter (1.5 kb before ATG) and exonic regions of genes were retrieved using the RiceSNPSeekDatabase website (https://snp-seek.irri.org/_snp.zul). Haplotype analysis was done by the DnaSP software application. The number of materials needs to be ten or more to be classified as different haplotypes.

2.6 qRT-PCR and sequence analysis

Using qRT-PCR analysis in both standard and salinized conditions, the expression levels of the candidate genes associated with both CD and WD20342 were confirmed. Roche LightCycler96 was used to carry out the qRT-PCR analysis. Supplementary Table S3 contains all of the primer sequences. The promoter and exonic regions of genes were amplified using PCR.

2.7 LOC_Os02g36880 mutant plants

By utilizing CRISPR/Cas9 technology, the mutants were created, and the T1 generation LOC_Os02g36880 knockout mutant seeds in the Zhonghua11 (ZH11) background were made by BIOGLE GENETECH company (http://www.biogle.cn/). PCR amplification was used to identify two homozygous mutants (CR-1, CR-15) in the T1 generation. These plants were then selfed to produce the T3 generation. The experiment was conducted with three replications using the wild-type (ZH11) and T3 generation homozygous mutants as the experimental materials. The experimental materials were separated into two groups. One group got a regular nutrient solution when they reached two leaves and one heart, while the other group received a nutrient solution with 120 mM NaCl. The therapy lasted 7 days, followed by a recovery period of 10 days. The treatment lasted for 7 days, and the recovery period lasted for 10 days. Salt tolerance was determined by measuring the SST and SSR.

3 Results

3.1 Phenotypes of GWAS and RIL population

Data from 295 Japonica rice materials and RIL lines were analyzed statistically to determine four salt tolerance phenotypes: RRL, RRN, RRFW, and RRDW. Under the salt stress during the seedling stage, the RRL, RRN, RRFW, and RRDW among the 295 materials varied in range from 0.39-0.88, 0.41-0.88, 0.38-0.87, and 0.39-0.87, respectively, with corresponding coefficients of variation of 13.65%, 12.79%, 14.74%, and 14.71% (Supplementary Table S1). The ranges of variation for RRL, RRN, RRFW, and RRDW in the RIL lines were 0.37-0.87, 0.44-0.88, 0.40-0.89, and 0.39-0.89, with the coefficients of variation of 17.05%, 13.89%, 17.22%, and 17.78%, respectively (Supplementary Table S1). The phenotypic data demonstrate a wide range of phenotypic diversity to Japonica rice seedling salt tolerance among the 295 materials and RIL lines. The phenotypic values of four traits significantly differed between WD20342 and CD, revealing that WD20342 was more salinity-tolerant than CD (Supplementary Table S1). The RRL, RRN, RRFW, and RRDW phenotypic values of the 295 materials and RIL lines were typically distributed and showed continuous phenotypic variation for each trait, suggesting that these traits are quantitative and affected by multiple factors (Figure 1).

Figure 1
www.frontiersin.org

Figure 1. Phenotypic variation in the RRL, RRN, RRFW, and RRDW. (A-D) The RRL, RRN, RRFW, and RRDW of 295 materials. (E-H) The RRL, RRN, RRFW, and RRDW of RIL lines.

3.2 GWAS of salt-related traits in the seedling stage

GWAS was performed using the 788,396 SNPs acquired from previous research (Duan et al., 2022; Li et al., 2022). Fifteen lead SNPs were significantly related to RRL, RRN, RRFW, and RRDW by GWAS (Table 1; Figure 2). Four QTLs linked to RRL were identified on chromosomes 2, 3, 6, and 11, with R2 values varying between 9.52% and 11.30%. Three QTLs linked to RRN were identified on chromosomes 6, 10, and 11, with R2 values ranging from 8.35% to 11.12%. Five QTLs linked to RRFW were identified on chromosomes 4, 5, 6, 7, and 12, with R2 values ranging from 10.14% to 12.21%. Meanwhile, three QTLs linked to RRDW were identified on chromosomes 4, 6, and 9, with R2 values ranging from 8.42% to 10.27%.

Table 1
www.frontiersin.org

Table 1. Lead SNPs associated with RRL, RRN, RRFW, and RRDW were identified by GWAS.

Figure 2
www.frontiersin.org

Figure 2. Manhattan plots and Q-Q plots. (A-D) Manhattan plot for the RRL, RRN, RRFW, and RRDW. (E-H) Q-Q plot for the RRL, RRN, RRFW, and RRDW.

3.3 QTL mapping of RIL

Six major QTLs were found on chromosomes 1, 2, 6, 10, and 11 in the RIL lines. The explained phenotypic variance varied from 4.82% to 17.26%, while the LOD values ranged from 2.55 to 7.57. There was one QTL associated with RRL (qRRL2), one QTL associated with RRN (qRRN6), two QTLs associated with RRFW (qRRFW1 and qRRFW10), and two QTLs associated with RRDW (qRRDW2 and qRRDW11) (Table 2). Furthermore, qRRL2 and qRRDW2 were defined as the same QTL because of the same interval, located in the physical region between markers C2_21864234 and C2_24239570 and explaining 9.38-17.26% of the phenotypic variance.

Table 2
www.frontiersin.org

Table 2. QTLs for RRL, RRN, RRFW, and RRDW were identified by linkage analysis.

Significantly, the GWAS discovered the lead SNP related to RRL (Chr2_22340368), which was located within qRRL2 and qRRDW2. On chromosome 2, the LD block region was estimated to be 22.241-22.437 Mb (196 kb). Using GWAS and linkage analysis, a 196-kb overlap region was identified (Figure 3). Therefore, we identified the 196-kb overlap region as the significant QTL for candidate gene identification.

Figure 3
www.frontiersin.org

Figure 3. Co-location results of GWAS and linkage analysis. (A) The local Manhattan plot and LD heatmap surround the lead SNP. (B) Mapped to the interval between markers C2_21864234 and C2_24239570 by linkage analysis. (C) The lead SNP (C2_22340368) was identified by the GWAS. (D) The 196-kb region.

3.4 Haplotype analysis of candidate genes

The Rice Genome Annotation Project confirms that the 196-kb region contained 22 genes (Supplementary Table S2). We conducted haplotype analysis on 22 genes and identified 4 genes (LOC_Os02g36880, LOC_Os02g36950, LOC_Os02g37000, LOC_Os02g37080) in the overlapping region. LOC_Os02g36880 was classified into three haplotypes based on non-synonymous mutant SNPs in the exon region. LOC_Os02g36950, LOC_Os02g37000, and LOC_Os02g37080 were separated into two haplotypes due to non-synonymous mutation SNPs in the exon region (Figure 4).

Figure 4
www.frontiersin.org

Figure 4. Haplotype analysis of candidate genes. (A-D) The haplotype analysis of LOC_Os02g36880, LOC_Os02g36950, LOC_Os02g37000, and LOC_Os02g37080. (The * and ** suggest significance of ANOVA at P < 0.05 and P < 0.01, respectively).

3.5 Identification of candidate genes

Both parents, CD and WD20342, were treated with 120 mM NaCl for 0, 1, 3, 6, 12, and 24 h, respectively. Four genes (LOC_Os02g36880, LOC_Os02g36950, LOC_Os02g37000, LOC_Os02g37080) were determined by qRT-PCR analysis. LOC_Os02g36950 was no induction at all under salt stress, and expression of the other genes (LOC_Os02g36880, LOC_Os02g37000, and LOC_Os02g37080) were up-regulated under stress. Interestingly, only LOC_Os02g36880 showed significant upregulation of CD and WD20342 expression under salt stress, with significant differences in expression levels. In the 24h salt stress, the expression level of LOC_Os02g36880 is more than four times higher in CD than WD20342 (Figure 5). Meanwhile, qRT-PCR analysis of other functionally annotated genes within the 196-kb region showed that the eight genes expressed between the two parents were not differentially expressed (Supplementary Figure S3).

Figure 5
www.frontiersin.org

Figure 5. The candidate gene expression patterns were confirmed. (A-D) The expression of LOC_Os02g36880, LOC_Os02g36950, LOC_Os02g37000, and LOC_Os02g37080 genes under control conditions and salt stress was analyzed. (** P < 0.01, *** P < 0.001, Students’t-test).

The sequencing of promotor and exonic regions of the genes (LOC_Os02g36880, LOC_Os02g36950, LOC_Os02g37000, and LOC_Os02g37080) was conducted for the two parents, CD and WD20342, after analyzing the initial findings. LOC_Os02g36880 in the CD (salt-sensitive) matches the reference sequence from Nipponbare, as shown by the findings. Contrarily, LOC_Os02g36880 in WD20342 (salt-tolerant) had one SNP mutation in the promoter region, 6-bp insertion mutations in the first exon region, one SNP mutation in the second exon region, and two SNP mutations in the third exon region (Supplementary Figure S4). The three genes (LOC_Os02g36950, LOC_Os02g37000, and LOC_Os02g37080) in CD and WD20342 was not differ regarding mutations. Consequently, LOC_Os02g36880 was found to be the functional gene in salt tolerance, and its function was further verified.

3.6 Validation of the LOC_Os02g36880 mutant

We selected two homozygous CRISPR/Cas9 mutants, CR-1 with a one-base (A) insertion and CR-15 with a one-base (T) deletion (Figure 6) to investigate whether LOC_Os02g36880 was associated with salt tolerance. These plants were propagated to the T3 generation for further analysis. The analysis found no significant difference in wild-type and CRISPR/Cas9 mutant seedlings under normal conditions. While treated with salt, mutants CR-1 and CR-15 showed higher salt tolerance than wild-type, with average SST values of 3.4 and 4.2, and SSR values of 77.0% and 72.0%, respectively. In contrast, the wild-type was an average SST of 7.8 and an average SSR of 27.0% (Figure 7). The findings suggest that the knockout mutant of LOC_Os02g36880 enhances salt tolerance in seedlings.

Figure 6
www.frontiersin.org

Figure 6. Target sequences of WT (ZH11) and knockout mutants (CR-1, CR-15).

Figure 7
www.frontiersin.org

Figure 7. Salt tolerance phenotype of WT and knockout mutants. (A) Phenotype of control conditions; (B) Phenotype of the salt stress conditions; (C) The SST under salt stress; (D) The SSR under salt stress. (** P < 0.01; Students’ t-test).

4 Discussion

The combination of GWAS and linkage analysis to locate gene loci for traits can effectively improve the accuracy and precision of association loci and obtain reliable target intervals. In recent years, this method has been widely applied in genetic studies of rice traits, such as salt tolerance (Kong et al., 2021), alkalinity tolerance (Li et al., 2020), disease resistance (Xiao et al., 2019), grain size (Kang et al., 2020), flag leaf trait (Wang et al., 2022). In this study, we combined these two approaches to characterize seedling salt tolerance in 295 materials and 189 RIL lines and obtained fifteen and six major QTLs, respectively. Comparing these genes to the findings of earlier research, they were near or overlapped with the loci of several known QTLs/genes. Ammar et al. (2009) found a salt tolerance-associated QTL (qNaLV-3.1), and the lead SNP (Chr3_25974535), detected by GWAS, were located within qNaLV-3.1. A QTL (qRFW-6) associated with the radicle fresh weight under salt stress was discovered by Mardani et al. (2014) on chromosome 6, and it contained the lead SNPs Chr6_19195187, Chr6_12954487, and Chr6_15555901 that were identified by GWAS. Atefeh et al. (2013) found a QTL (qRL-10b) on chromosome 10 associated with root length under osmotic stress, where the lead SNP (Chr10_4241399) identified by GWAS was located. The lead SNP (Chr12_19562318) detected by GWAS was located within the QTL (qSNC-12) for the concentrations of Na+ in shoots under salt stress detected by Zheng et al. (2015) using linkage analysis, and was close to the gene regulating iron distribution, OsbHLH133 (Wang et al., 2013). The QTL qRRN6 contained known genes OsCYP19-4 (Yoon et al., 2016), which has peptidyl-prolyl cis-trans isomerase activity, and, OsCYP19-4 overexpression plants enhance cold tolerance. Meanwhile, the GWAS identified the lead SNP (Chr2_22340368), which is located within qRRL2 and qRRDW2 identified by linkage analysis, both within the QTL (qNAK-2) associated with salt tolerance identified by Yao et al. (2005).

We found that LOC_Os02g36880, a NAC transcription factor, contains an apical meristem structural domain. NAC transcription factors have a role in regulating abiotic stress tolerance, leaf senescence, lateral root growth, programmed cell death, etc (Huysmans et al., 2018; Xia et al., 2018; Zhang et al., 2018; Singh et al., 2021; Wang et al., 2021). For example, the OsNAC106 gene is a NAC transcription factor that negatively regulates salt tolerance, and the osnac106 deletion mutant is more salinity tolerant than normal plants (Sakuraba et al., 2015). Nakashima et al. (2007) found that transgenic overexpressing OsNAC6 had increased tolerance to high salt and drought stresses. Hong et al. (2016) identified ONAC022 as a transcriptional activator that positively regulates salt tolerance by modulating abscisic acid signaling, and overexpression of OsNAC022 transgenic plants were more salt-tolerant, accumulating less Na+ in roots and shoots than wild-type. It was shown that the NAC family of transcription factors is crucial in the study of salt tolerance.

In the previous study, Fang et al (Fang et al., 2014). found that the expression level of LOC_Os02g36880 was significantly reduced under drought and increased under high salt stress. Wang et al. (2023) combined metabolomic and transcriptomic profiles and found that LOC_Os02g36880 may be associated with low-nitrogen tolerance. The gene LOC_Os02g36880 was determined to be the most likely functional gene linked to salt tolerance. Seedlings of the CRISPR/Cas9 mutant of the LOC_Os02g36880 gene were substantially more saline-tolerant than wild-type, demonstrating that LOC_Os02g36880 negatively regulates salt tolerance. We combined with haplotype analysis and found that CD belongs to Hap1 (AC) and WD20342 belongs to Hap2 (GA), as WD20342 is more salt-tolerant than CD, indicating that Hap2 belongs to the dominant salt tolerant haplotype. This provides a theoretical basis for breeding new salt-tolerant rice varieties and is essential for improving the utilization efficiency of saline land.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

SX: Formal analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing. JZ: Formal analysis, Investigation, Validation, Writing – original draft, Writing – review & editing. HD: Formal analysis, Investigation, Validation, Writing – original draft, Writing – review & editing. XD: Formal analysis, Writing – review & editing. CL: Validation, Writing – review & editing. YD: Validation, Writing – review & editing. YC: Writing – review & editing. JW: Writing – review & editing. HL: Writing – review & editing. LY: Writing – review & editing. WX: Writing – review & editing. YJ: Writing – review & editing. DZ: Conceptualization, Supervision, Writing – review & editing. HZ: Conceptualization, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by Key Research and Development Projects of Heilongjiang Province (grant No. 2022ZX02B04), the National Natural Science Foundation of China (grant No. U20A2025), the “Academic Backbone” Project of Northeast Agricultural University (grant No. 20XG24), and the Mining of Salinity Tolerant Germplasm Resources and Variety Selection and Breeding of Coldland Early Japonica Rice (grant No. CZKYF2022-1-C005).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1462856/full#supplementary-material

Supplementary Figure 1 | Phenotypes of CD and WD20342 seedlings under control and salt stress. (A) Under control conditions. (B) Under salt stress. Bar = 5 cm.

Supplementary Figure 2 | Linkage analysis results of RIL lines derived from CD/WD20342.

Supplementary Figure 3 | Expression patterns of the other eight genes under normal growth and salt stress.

Supplementary Figure 4 | The gene structure and sequence difference analysis of LOC_Os02g36880. (A) The gene structure of LOC_Os02g36880. (B) The sequence difference analysis of LOC_Os02g36880 in WD20342, CD, and Nipponbare.

References

Al-Tamimi, N., Brien, C., Oakey, H., Berger, B., Saade, S., Ho, Y. S., et al. (2016). Salinity tolerance loci revealed in rice using high-throughput non-invasive phenotyping. Nat. Commun. 7, 11. doi: 10.1038/ncomms13342

Crossref Full Text | Google Scholar

Ammar, M. H. M., Pandit, A., Singh, R. K., Sameena, S., Chauhan, M. S., Singh, A. K., et al. (2009). Mapping of QTLs controlling Na+, K+ and Cl- ion concentrations in salt tolerant indica rice variety CSR27. J. Plant Biochem. Biotechnol. 18, 139–150. doi: 10.1007/bf03263312

Crossref Full Text | Google Scholar

Atefeh, Sabouri, Hossein, Sabouri, Marjorie, De, et al. (2013). Genetic analysis seedling 333 vigour under osmotic stress in rice by QTL mapping. Russ. Agric. Sci.

Google Scholar

Batayeva, D., Labaco, B., Ye, C. R., Li, X. L., Usenbekov, B., Rysbekova, A., et al. (2018). genome-wide association study of seedling stage salinity tolerance in temperate japonica rice germplasm. BMC Genet. 19, 11. doi: 10.1186/s12863-017-0590-7

PubMed Abstract | Crossref Full Text | Google Scholar

Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308

PubMed Abstract | Crossref Full Text | Google Scholar

Chang, C. C., Chow, C. C., Tellier, L., Vattikuti, S., Purcell, S. M., Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4, 16. doi: 10.1186/s13742-015-0047-8

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, G., Hu, K. M., Zhao, J. H., Guo, F. F., Shan, W. F., Jiang, Q. Q., et al. (2022). genome-wide association analysis for salt-induced phenotypic and physiologic responses in rice at seedling and reproductive stages. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.822618

Crossref Full Text | Google Scholar

De Leon, T. B., Linscombe, S., Subudhi, P. K. (2016). Molecular dissection of seedling salinity tolerance in rice (Oryza sativa L.) using a high-density GBS-based snp linkage map. Rice 9, 22. doi: 10.1186/s12284-016-0125-2

PubMed Abstract | Crossref Full Text | Google Scholar

Dong, S. S., He, W. M., Ji, J. J., Zhang, C., Yang, T. L. (2020). LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files. Briefings Bioinf. 22, 6. doi: 10.1101/2020.06.14.151332

Crossref Full Text | Google Scholar

Duan, Y. X., Zheng, H. L., Wen, H. R., Qu, D., Cui, J. N., Li, C., et al. (2022). Identification of candidate genes for salt tolerance at the germination stage in japonica rice by genome-wide association analysis. Agriculture-Basel 12, 15. doi: 10.3390/agriculture12101588

Crossref Full Text | Google Scholar

Fang, Y. J., Xie, K. B., Xiong, L. Z. (2014). Conserved miR164-targeted NAC genes negatively regulate drought resistance in rice. J. Exp. Bot. 65, 2119–2135. doi: 10.1093/jxb/eru072

PubMed Abstract | Crossref Full Text | Google Scholar

Geng, L. Y., Zhang, W., Zou, T., Du, Q., Ma, X. D., Cui, D., et al. (2023). Integrating linkage mapping and comparative transcriptome analysis for discovering candidate genes associated with salt tolerance in rice. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1065334

Crossref Full Text | Google Scholar

Hanin, M., Ebel, C., Ngom, M., Laplaze, L., Masmoudi, K. (2016). New insights on plant salt tolerance mechanisms and their potential use for breeding. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.01787

PubMed Abstract | Crossref Full Text | Google Scholar

Hong, Y. B., Zhang, H. J., Huang, L., Li, D. Y., Song, F. M. (2016). Overexpression of a stress-responsive NAC transcription factor gene ONACO22 improves drought and salt tolerance in rice. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.00004

Crossref Full Text | Google Scholar

Huysmans, M., Buono, R. A., Skorzinski, N., Radio, M. C., De Winter, F., Parizot, B., et al. (2018). NAC transcription factors ANAC087 and ANAC046 control distinct aspects of programmed cell death in the arabidopsis columella and lateral root cap. Plant Cell 30, 2197–2213. doi: 10.1105/tpc.18.00293

PubMed Abstract | Crossref Full Text | Google Scholar

Kang, J. W., Kabange, N. R., Phyo, Z., Park, S. Y., Lee, S. M., Lee, J. Y., et al. (2020). Combined linkage mapping and genome-wide association study identified QTLs associated with grain shape and weight in rice (Oryza sativa L.). Agronomy-Basel 10, 22. doi: 10.3390/agronomy10101532

Crossref Full Text | Google Scholar

Kong, W. L., Zhang, C. H., Zhang, S. C., Qiang, Y. L., Zhang, Y., Zhong, H., et al. (2021). Uncovering the novel QTLs and candidate genes of salt tolerance in rice with linkage mapping, RTM-GWAS, and RNA-seq. Rice 14, 12. doi: 10.1186/s12284-021-00535-3

PubMed Abstract | Crossref Full Text | Google Scholar

Kumar, V., Singh, A., Mithra, S. V. A., Krishnamurthy, S. L., Parida, S. K., Jain, S., et al. (2015). Genome-wide association mapping of salinity tolerance in rice (Oryza sativa). DNA Res. 22, 133–145. doi: 10.1093/dnares/dsu046

PubMed Abstract | Crossref Full Text | Google Scholar

Lei, L., Han, Z. H., Cui, B. W., Yang, L. M., Liu, H. L., Wang, J. G., et al. (2021). Mapping of a major QTL for salinity tolerance at the bud burst stage in rice (Oryza sativa L) using a high-density genetic map. Euphytica 217, 8. doi: 10.1007/s10681-021-02901-0

Crossref Full Text | Google Scholar

Li, J., Xin, W., Wang, W. P., Zhao, S. J., Xu, L., Jiang, X. D., et al. (2022). Mapping of candidate genes in response to low nitrogen in rice seedlings. Rice 15, 16. doi: 10.1186/s12284-022-00597-x

PubMed Abstract | Crossref Full Text | Google Scholar

Li, N., Zheng, H. L., Cui, J. N., Wang, J. G., Liu, H. L., Sun, J., et al. (2019). Genome-wide association study and candidate gene analysis of alkalinity tolerance in japonica rice germplasm at the seedling stage. Rice 12, 12. doi: 10.1186/s12284-019-0285-y

PubMed Abstract | Crossref Full Text | Google Scholar

Li, S. S., Xu, S. B., Zheng, J., Du, H. Q., Li, C., Shen, S., et al. (2023). Joint QTL mapping and transcriptome sequencing analysis reveal candidate genes for salinity tolerance in oryza sativa L. ssp. japonica seedlings. Int. J. Mol. Sci. 24, 15. doi: 10.3390/ijms242417591

Crossref Full Text | Google Scholar

Li, X. W., Zheng, H. L., Wu, W. S., Liu, H. L., Wang, J. G., Jia, Y., et al. (2020). QTL mapping and candidate gene analysis for alkali tolerance in japonica rice at the bud stage based on linkage mapping and genome-wide association study. Rice 13, 11. doi: 10.1186/s12284-020-00412-5

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, C. T., Mao, B. G., Yuan, D. Y., Chu, C. C., Duan, M. J. (2022). Salt tolerance in rice: physiological responses and molecular mechanisms. Crop J. 10, 13–25. doi: 10.1016/j.cj.2021.02.010

Crossref Full Text | Google Scholar

Mardani, Z., Rabiei, B., Sabouri, H., Sabouri, A. (2014). Identification of molecular markers linked to salt-tolerant genes at germination stage of rice. Plant Breed. 133, 196–202. doi: 10.1111/pbr.12136

Crossref Full Text | Google Scholar

Nakashima, K., Tran, L. S. P., Nguyen, D. V., Fujita, M., Maruyama, K., Todaka, D., et al. (2007). Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice. Plant J. 51, 617–630. doi: 10.1111/j.1365-313X.2007.03168.x

PubMed Abstract | Crossref Full Text | Google Scholar

Nakhla, W. R., Sun, W. Q., Fan, K., Yang, K., Zhang, C. P., Yu, S. B. (2021). Identification of QTLs for salt tolerance at the germination and seedling stages in rice. Plants-Basel 10, 12. doi: 10.3390/plants10030428

Crossref Full Text | Google Scholar

Phosuwan, S., Nounjan, N., Theerakulpisut, P., Siangliw, M., Charoensawan, V. (2024). Comparative quantitative trait loci analysis framework reveals relationships between salt stress responsive phenotypes and pathways. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1264909

PubMed Abstract | Crossref Full Text | Google Scholar

Qiu, X. J., Yuan, Z. H., Liu, H., Xiang, X. J., Yang, L. W., He, W. J., et al. (2015). Identification of salt tolerance-improving quantitative trait loci alleles from a salt-susceptible rice breeding line by introgression breeding. Plant Breed. 134, 653–660. doi: 10.1111/pbr.12321

Crossref Full Text | Google Scholar

Radha, B., Sunitha, N. C., Sah, R. P., Azharudheen, M. T. P., Krishna, G. K., Umesh, D. K., et al. (2023). Physiological and molecular implications of multiple abiotic stresses on yield and quality of rice. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.996514

PubMed Abstract | Crossref Full Text | Google Scholar

Raj, S. R. G., Nadarajah, K. (2023). QTL and candidate genes: techniques and advancement in abiotic stress resistance breeding of major cereals. Int. J. Mol. Sci. 24, 35. doi: 10.3390/ijms24010006

Crossref Full Text | Google Scholar

Rohila, J. S., Edwards, J. D., Tran, G. D., Jackson, A. K., McClung, A. M. (2019). Identification of superior alleles for seedling stage salt tolerance in the USDA rice mini-core collection. Plants-Basel 8, 23. doi: 10.3390/plants8110472

PubMed Abstract | Crossref Full Text | Google Scholar

Ruan, S. L., Ma, H. S., Wang, S. H., Fu, Y. P., Xin, Y., Liu, W. Z., et al. (2011). Proteomic identification of OsCYP2, a rice cyclophilin that confers salt tolerance in rice (Oryza sativa L.) seedlings when overexpressed. BMC Plant Biol. 11, 15. doi: 10.1186/1471-2229-11-34

PubMed Abstract | Crossref Full Text | Google Scholar

Sakuraba, Y., Piao, W., Lim, J. H., Han, S. H., Kim, Y. S., An, G., et al. (2015). Rice ONAC106 inhibits leaf senescence and increases salt tolerance and tiller angle. Plant Cell Physiol. 56, 2325–2339. doi: 10.1093/pcp/pcv144

PubMed Abstract | Crossref Full Text | Google Scholar

Sarma, B., Kashtoh, H., Tamang, T. L., Bhattacharyya, P. N., Mohanta, Y. K., Baek, K. H., et al. (2023). Abiotic stress in rice: visiting the physiological response and its tolerance mechanisms. Plants-Basel 12, 41. doi: 10.3390/plants12233948

Crossref Full Text | Google Scholar

Singh, S., Koyama, H., Bhati, K. K., Alok, A. (2021). The biotechnological importance of the plant-specific NAC transcription factor family in crop improvement. J. Plant Res. 134, 475–495. doi: 10.1007/s10265-021-01270-y

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Z. F., Cheng, J. P., Chen, Z. W., Huang, J., Bao, Y. M., Wang, J. F., et al. (2012). Identification of QTLs with main, epistatic and QTL x environment interaction effects for salt tolerance in rice seedlings under different salinity conditions. Theor. Appl. Genet. 125, 807–815. doi: 10.1007/s00122-012-1873-z

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y. R., Jiang, J., Qian, Y. K., Miao, S. Y., Wang, W. S., Xu, J. L., et al. (2023). Multi-Omics analysis reveals the regulatory and metabolic mechanisms underlying low-nitrogen tolerance at the flowering stage in rice. Agronomy-Basel 13, 16. doi: 10.3390/agronomy13020578

Crossref Full Text | Google Scholar

Wang, J. X., Wang, T., Wang, Q., Tang, X. D., Ren, Y., Zheng, H. Y., et al. (2022). QTL mapping and candidate gene mining of flag leaf size traits in japonica rice based on linkage mapping and genome-wide association study. Mol. Biol. Rep. 49, 63–71. doi: 10.1007/s11033-021-06842-8

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, L., Ying, Y. H., Narsai, R., Ye, L. X., Zheng, L. Q., Tian, J. L., et al. (2013). Identification of OsbHLH133 as a regulator of iron distribution between roots and shoots in oryza sativa. Plant Cell Environ. 36, 224–236. doi: 10.1111/j.1365-3040.2012.02569.x

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, H. L., Zhang, Y., Wang, T., Yang, Q., Yang, Y. L., Li, Z., et al. (2021). An alternative splicing variant of PtRD26 delays leaf senescence by regulating multiple NAC transcription factors in Populus. Plant Cell 33, 1594–1614. doi: 10.1093/plcell/koab046

PubMed Abstract | Crossref Full Text | Google Scholar

Xia, J. N., Zhao, Y. J., Burks, P., Pauly, M., Brown, P. J. (2018). A sorghum NAC gene is associated with variation in biomass properties and yield potential. Plant Direct 2, 11. doi: 10.1002/pld3.70

Crossref Full Text | Google Scholar

Xiao, S. Z., Wang, B. X., Liu, Y. Q., Miao, T. H., Zhang, H. L., Wen, P. Z., et al. (2019). Genome-wide association study and linkage analysis on resistance to rice black-streaked dwarf virus disease. Mol. Breed. 39, 10. doi: 10.1007/s11032-019-0980-9

Crossref Full Text | Google Scholar

Yao, M. T., Wang, J. F., Chen, H. Y., Zhai, H. Q., Zhang, H. S. (2005). Inheritance and QTL mapping of salt tolerance in rice. Rice Sci. 12, 25–32.

Google Scholar

Yoon, D. H., Lee, S. S., Park, H. J., Lyu, J. I., Chong, W. S., Liu, J. R., et al. (2016). Overexpression of OsCYP19-4 increases tolerance to cold stress and enhances grain yield in rice (Oryza sativa). J. Exp. Bot. 67, 69–82. doi: 10.1093/jxb/erv421

PubMed Abstract | Crossref Full Text | Google Scholar

Yuan, J., Wang, X. Q., Zhao, Y., Khan, N. U., Zhao, Z. Q., Zhang, Y. H., et al. (2020). Genetic basis and identification of candidate genes for salt tolerance in rice by GWAS. Sci. Rep. 10, 9. doi: 10.1038/s41598-020-66604-7

PubMed Abstract | Crossref Full Text | Google Scholar

Zeng, P., Zhu, P. W., Qian, L. F., Qian, X. M., Mi, Y. X., Lin, Z. F., et al. (2021). Identification and fine mapping of qGR6.2, a novel locus controlling rice seed germination under salt stress. BMC Plant Biol. 21, 14. doi: 10.1186/s12870-020-02820-7

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, L., Yao, L., Zhang, N., Yang, J. W., Zhu, X., Tang, X., et al. (2018). Lateral root development in potato is mediated by stu-mi164 regulation of NAC transcription factor. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00383

Crossref Full Text | Google Scholar

Zheng, H. L., Zhao, H. W., Liu, H. L., Wang, J. G., Zou, D. T. (2015). QTL analysis of Na+ and K+ concentrations in shoots and roots under NaCl stress based on linkage and association analysis in japonica rice. Euphytica 201, 109–121. doi: 10.1007/s10681-014-1192-3

Crossref Full Text | Google Scholar

Zhu, J. K. (2016). Abiotic stress signaling and responses in plants. Cell 167, 313–324. doi: 10.1016/j.cell.2016.08.029

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Japonica rice seedling, salt tolerance, QTL, candidate gene, GWAS, linkage analysis

Citation: Xu S, Zheng J, Du H, Du X, Li C, Duan Y, Cai Y, Wang J, Liu H, Yang L, Xin W, Jia Y, Zou D and Zheng H (2024) GWAS combined with linkage analysis reveals major QTLs and candidate genes of salt tolerance in Japonica rice seedlings. Front. Plant Sci. 15:1462856. doi: 10.3389/fpls.2024.1462856

Received: 10 July 2024; Accepted: 25 September 2024;
Published: 01 November 2024.

Edited by:

Meng Jiang, Zhejiang University, China

Reviewed by:

Biaolin Hu, Jiangxi Academy of Agricultural Sciences, China
Xiaodeng Zhan, Chinese Academy of Agricultural Sciences, China

Copyright © 2024 Xu, Zheng, Du, Du, Li, Duan, Cai, Wang, Liu, Yang, Xin, Jia, Zou and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hongliang Zheng, zhenghongliang008@126.com; Detang Zou, zoudtneau@126.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.