Skip to main content

Genome-wide association study of seedling stage salinity tolerance in temperate japonica rice germplasm



Salinity has a significant impact on rice production in coastal, arid and semi-arid areas in many countries, including countries growing temperate rice, such as Kazakhstan. Recently, the complete genomes of 3000 rice accessions were sequenced through the 3 K rice genome project, and this set included 203 temperate japonica rice accessions. To identify salinity-tolerant germplasm and related genes for developing new salinity-tolerant breeding lines for the temperate japonica rice growing regions, we evaluated the seedling stage salinity tolerance of these sequenced temperate japonica rice accessions, and conducted genome-wide association studies (GWAS) for a series of salinity tolerance related traits.


There were 27 accessions performed well (SES < 5.0) under moderate salinity stress (EC12), and 5 accessions were tolerant under both EC12 and EC18. A total of 26 QTLs were identified for 9 measured traits. Eleven of these QTLs were co-located with known salinity tolerance genes. QTL/gene clusters were observed on chromosome 1, 2, 3, 6, 8 and 9. Six candidate genes were identified for five promising QTLs. The alleles of major QTL Saltol and gene O S HKT1;5 (SKC1) for Na+/K+ ratio identified in indica rice accessions were different from those in the temperate japonica rice accessions used in this study.


Salinity tolerant temperate japonica rice accessions were identified in this study, these accessions are important resources for breeding programs. SNPs located in the promising QTLs and candidate genes could be used for future gene validation and marker assisted selection. This study provided useful information for future studies on genetics and breeding of salinity tolerance in temperate japonica rice.


Rice is grown worldwide in areas where temperatures are suitable and irrigation water is available. This includes many areas with soil salinity problems. Rice is rated as a salt-sensitive crop, especially at early seedling and reproductive stages [1, 2]. The threshold for salt stress in rice is 3 dS/m, with a 12% reduction in yield for each unit increase in dS/m beyond this value [3]. The narrow genetic pool of japonica germplasm for salinity tolerance and the extent of the salinity challenge are common to many rice-growing regions. Rice production is limited by soil salinity problems in many tropical countries such as Bangladesh and India, and temperate countries such as Korea, Japan, Iran and Kazakhstan [4]. For example, the rice growing regions in Kazakhstan are located in the most northern regions of rice cultivation on earth (from Kyzyl-Orda 44°48’N to Balkhash 46°00’N) with total area of about 100,000 ha and paddy production of 300,000–500,000 tons per year. Soil salinity and low temperatures at the early growing stage are the main environmental problems for rice production in Kazakhstan [5]. However, most of the previous studies on rice salinity tolerance have been conducted in tropical regions using indica accessions. Some salinity-tolerant rice accessions (Nona Bokra, Pokkali, PSBRc50, FL478) have been identified and used in breeding programs, and many salinity-tolerant indica rice varieties have been developed. IRRI released a significant number of salt-tolerant varieties from 2011 to 2013. These include three varieties released in Bangladesh (Bina dhan 10, BRRI dhan 55 and 61), one variety in India (Luna Sankhi), eight ‘Salinas’ varieties in the Philippines and one variety in Gambia [6]. However, very few salinity-tolerant japonica rice accessions have been identified, such as Harra (Spain), Agami (Egypt) and Daeyabyeo (Korea) [7]. Thus, it is important to identify salinity-tolerant japonica rice germplasm and use them to develop salinity-tolerant japonica varieties for temperate rice growing regions.

Rice is very sensitive to salinity at the seedling stage. Soil salinity usually causes stunting or even death of seedlings [8]. Based on the ability of seedlings to grow in salinized nutrient solution, a screening protocol was developed and widely used for germplasm screening, breeding selection and genetic studies [9]. The typical mechanism of salinity tolerance in rice is the Na+ exclusion or uptake reduction, and increased absorption of K+ to maintain a suitable Na+-K+ balance in the shoot [10]. Thus, shoot Na+, K+ content and Na+/K+ ratio are also used as valid criterion in measuring salinity tolerance in rice [9].

Salinity tolerance is a complex trait controlled by quantitative trait loci (QTLs). Many QTLs for salinity tolerance have been identified [11,12,13,14,15,16,17,18,19,20,21,22,23,24]. However, most of the QTLs were identified in bi-parental populations derived from indica accessions such as Pokkali and Nona Bokra, very few QTLs for salinity tolerance have been identified in temperate japonica accessions [15, 25]. A recent study showed that significant differences in phenotypic response to salinity exist within the japonica accessions of European Rice Core collection (ERCC) [26]. More genetic studies should be carried out to understand the salinity tolerance QTLs in temperate japonica accessions and to facilitate the subsequent use of these QTLs in breeding programs.

Recent development of high density markers such as SNPs enables identification of trait-marker association through association mapping, and GWAS for many important agronomic traits have been reported [27,28,29]. With the collaboration among the Chinese Academy of Agricultural Sciences (CAAS), the Beijing Genomics Institute (BGI) and the International Rice Research Institute (IRRI), 3000 rice accessions of known diversity were systematically sequenced [30]. After alignment with the reference genome (Nipponbare), around 18.9 million SNPs were identified. This database provided valuable information for understanding the genetic mechanism of useful traits in these diversity accessions.

The objectives of this study were to screen a set of temperate japonica accessions which had been fully sequenced, to identify SNPs associated with seedling stage salinity tolerance, and to identify potential candidate genes for the promising QTLs. We selected the subset of temperate japonica accessions from the 3000 sequenced accessions, evaluated the seedling stage salinity tolerance of these accessions using 16 different traits (indexes), and conducted genome-wide association studies between these traits and core SNP markers.


Plant materials

The 3 K Rice Genomes Project sequenced 3000 rice genomes with an average sequencing depth of 14×. The dataset includes publicly available genome sequences derived from 3000 accessions of rice with global representation (from 89 countries) of genetic and functional diversity [30]. Phylogenetic analyses based on SNP data confirmed differentiation of the O. sativa gene pool into 5 varietal groups – indica, aus/boro, basmati/sadri, tropical japonica and temperate japonica [31]. There were 203 accessions in the temperate japonica group (Additional file 1: Figure S1). We requested all the sequenced temperate japonica accessions, however, seed availability was limited to 191 accessions, and these accessions were used in this study (Additional file 2: Table S1).

Evaluation of salinity tolerance at seedling stage

The seeds of selected accessions were treated at 50 °C for 5 days in an oven to break seed dormancy. Two pre-germinated seeds were sown in each well on a styrofoam seedling float that was placed on a tray filled with distilled water. Each accession was sown into four wells, and each seeding float contained 20 accessions plus one tolerant check variety (FL478) and one susceptible check variety (NSIC Rc222). The experiment was conducted in the phytotron glasshouse maintained at 29/21 °C day/night temperature and minimum relative humidity of 70%. The experimental design was a randomized complete block design (RCBD) with 2 treatments (control and salt treatment), and 3 replications for each treatment. For the salt treatment, four days after seeding (DAS), we replaced the distilled water with salinized Peter’s Professional™ nutrient solution (20 N-20P-20 K) as hydroponic solution at a rate of 1 g per liter water and with 400 mg of ferrous sulfate, the electrical conductivity was adjusted to EC = 6 dSm−1 by adding NaCl to the nutrient solution, and the pH was adjusted to 5.1 daily. At 7 DAS, the hydroponic solution was changed to maintain at EC = 12 dSm−1. After two weeks in EC 12 dSm−1, the visual reactions of plant to salinity stress were evaluated using percentage of damage (PD) and the Standard Evaluation Score (SES) for salinity tolerance (Additional file 3: Table S2) [32]. After scoring, one seedling of each replicate was sampled for measurement of shoot Na+ (SNa) and K+ (SK) content. Then the electrical conductivity of the hydroponic solution was changed to maintain at EC = 18 dSm−1. A final scoring of PD and SES was done after 7 days exposure to EC = 18 dSm−1. The control was always maintained in distilled water (EC0), and all samplings and measurements were the same as salt treatment.

Measurement of Na+ and K+ content

One seedling from each replicate was sampled and washed 3 times with distilled water. The seedling was put in paper bags and dried in an oven at 70 °C for at least 3 days. Dry weight of the seedling (SDW) was measured, and then ten milligrams of dried leaf were cut into small crumbs (<1 cm) and placed in 10 ml Falcon tubes, 10 ml 0.1 N acetic acid (CH3COOH) was added to each tube, then heated in a waterbath at 90 °C for 2 h. The solution was filtered through Whatman No1 filter paper placed in a glass funnel, and the filtration was collected in a 10 ml flask as stock solution. The stock solution was then diluted 10 times (1 ml of stock solution in 9 ml of nanopure water), then sodium and potassium content were measured using a flame spectrophotometer (Sherwood Model 420).

Genome sequence data

Approximately 18.9 million single nucleotide polymorphisms (SNPs) were discovered in 3000 rice genomes when aligned to the reference genome of the temperate japonica accession Nipponbare [31]. The SNPs and allele information were organized into a SNP-Seek system [33, 34]. The core SNP V2.1 of the selected 191 temperate japonica rice accessions were downloaded from the SNP-seek system and used for data analysis.

Data analysis

For the phenotypic data of salinity tolerance related traits, the mean value of three replications was calculated and used for genome-wide association study. The relative reduction rate of each trait was calculated as (control-treatment)×100/control. Basic statistical information of the traits is shown in Additional file 4: Table S3. The correlations among different traits were calculated by using MINITAB V14.0 (Minitab Inc.).

The core SNP dataset of the selected 191 temperate japonica accessions included 365,710 SNPs. There were a large proportion of missing calls along with many heterozygous SNPs detected in the dataset. Fifteen accessions with high heterozygotes (>1%) and high missing sites (>25%) were removed. The Trait Analysis by Association Evolution and Linkage (TASSEL) program version 5.2.18 [35], was used to filter the sites at a maximum count of 158 of the 176 remaining accessions, which accounts for sites in which 90% of the accessions have a call and a minimum frequency of 0.05 for the minor allele. The above criteria resulted in 68,786 filtered sites. Finally 176 accessions with 68,786 SNPs were used for making cladogram tree with neighbor-joining and generating kinship matrix with centered IBS (default) [36]. Principal components analysis was done using these filtered 68,786 SNPs with default settings. A united data file with the genotype and phenotype of the accessions was created by using union join. The united file along with kinship matrix was used to analyze marker-trait associations using a mixed linear model (MLM) [37]. The compression level was set to optimum level, and variance component estimation was set to P3D. The significant threshold was set at p < 0.0001 (−log10 p-value >4.0) and/or SNP marker R2 ≥ 0.1 [38]. The QQ plot is shown in Additional file 5: Figure S2. The identified QTLs were named using the CGSNL nomenclature [39], and mapped on the rice genome using Mapchart 2.30 [40]. A promising QTL was considered when many SNPs were lined up near the peak of the QTL. The chromosomes and QTL regions of the promising QTLs were selected and re-analyzed using general linear model (GLM) with 1000 times permutation. Candidate genes (near peak SNPs) for the promising QTLs were searched on the rice genome browser in the rice SNP seek database [34] and literature reports.


Diversity and population structure of selected temperate japonica rice accessions

There were 203 temperate japonica accessions included in the 3000 sequenced genomes. However, the seeds of 12 accessions were not available at the time of request, and 15 accessions showed high heterozygous SNP sites and a high rate of missing data. Thus, 176 temperate japonica accessions from 32 countries were used for analysis. The principal component analysis of the 176 accessions showed that all the accessions were closely linked in one group, and no significant sub-group was classified (Fig. 1). When the origin of the accessions were considered, accessions from most of the countries were randomly distributed, except that some of the Chinese and Korean accessions were closely linked (Additional file 6: Figure S3). Thus, population structure was not considered in the following GWAS analysis.

Fig. 1

Principal component analysis (PCA) of 176 temperate japonica rice varieties using 68,786 SNPs

Salinity tolerance of selected temperate japonica rice accessions

There was a wide range of variation in salinity tolerance at seedling stage among the 176 temperate japonica accessions. However, based on the visual evaluation, only 3 of the evaluated accessions (Nep Ngau from Vietnam, Bai Mang Ai Zhong from China and Shinchiku Iku 97 from Taiwan) were as tolerant as the tolerant check variety FL478 at EC12 (SES = 3), and none of them was as tolerant as FL478 at EC18 (Table 1). Twenty-seven salinity tolerant accessions were listed in Additional file 7: Table S4. Among them, only 7 accessions were tolerant (SES ≤ 5) at both EC12 and EC18 including Bai Mang Ai Zhong and Shinchiku Iku 97.

Table 1 Number of accessions in different SES categories

Among the evaluated traits, shoot dry weight under normal conditions (control) was not correlated with other traits under salt treatment except shoot dry weight (Table 2). Since SES scores and percentage of damage (PD) are based on visual evaluation of the injury of the seedlings, SES was correlated with PD. Shoot Na+ and K+ content under normal conditions (control) were not significantly correlated with PD and SES score, however leaf Na+ content and Na+/K+ ratio were correlated with PD and SES under salt stress. SES12 and SES18 are significantly correlated, scoring once at EC12 treatment is enough for evaluating the salinity tolerance of a variety. Shoot dry weight (include its relative deduction compared to control) and shoot Na+ content at EC12 treatment are significantly correlated with PD and SES, thus SDW12, SDWD and SNa12 are good indices for evaluating salinity tolerance of rice varieties. For EC12 treatment, shoot dry weight (SDW12) and its relative deduction (SDWD) are significantly correlated with shoot Na+ (SNa12), K+ (SK12) content and Na+/K+ ratio (SNa/K12), the change of Na+, K+ content and their ratio in the shoot directly affected the growth of the plant.

Table 2 Pearson correlation among salinity tolerance related traits

Genome-wide association study of traits related to seedling stage salinity tolerance

Percentage of damage and standard evaluation score

For the percentage of damage of the seedlings after salinity treatment at EC12 for 2 weeks, no significant QTLs was identified. However, three QTLs on chromosomes 2, 4 and 11 were identified for percentage of damage of the seedlings after salinity treatment at EC18 (Table 3, Additional file 8: Figure S4).

Table 3 Identified QTLs for the assessed salinity responsive traits

Same as PD, no QTL was identified for SES of EC12 treatment, but three QTLs were identified on chromosomes 2, 4 and 11 for EC18 treatment. Since the SES and PD were significantly correlated, QTLs identified for SES and PD were also located in the same regions on the chromosomes.

Shoot dry weight

There was a wide range of variation of shoot dry weight among the accessions. Under normal conditions (control), two QTLs related to shoot dry weight were identified on chromosomes 2 and 9. However, no QTL was identified for shoot dry weight after salinity treatment and the relative reduction of shoot dry weight.

Shoot Na+ and K+ content and Na+/K+ ratio

QTLs for shoot Na+ content were identified on chromosomes 2, 6, 8 and 9 under normal conditions (control) and on chromosomes 4 and 10 under salinity treatment (EC12). After salinity treatment, the leaf Na+ content increased comparing to control, and QTLs for this increment were identified on chromosomes 5 and 6.

For shoot K+ content, no QTL was identified for control and salinity treatment. After salinity treatment, the leaf K+ content increased, and QTLs for this increment were identified on chromosomes 1, 2, 3, 4 and 12.

QTLs for shoot Na+/K+ ratio were identified on chromosomes 3 and 11 under normal conditions (control), but no QTL was identified for salinity treatment (EC12). After salinity treatment, the leaf Na+/K+ ratio varied among accessions, with some showing higher ratios and some lower. One QTL for the change of Na+/K+ ratio were identified on chromosomes 3.

Some of the identified QTLs on chromosomes 2 (qPD18_2.1, qSES18_2.1 and qSKD_2.1), 3 (qSNa/K0_3.1 and qSKD_3.1), 4 (qPD18_4.1 and qSES18_4.1), 9 (qSNa0_9.1 and qSDW0_9.1) and 11 (qPD18_11.1 and qSES18_11.1) are located in the same position or very closely linked (Table 3, Additional file 9: Figure S5).

Candidate genes for promising QTLs

Based on the significant SNPs in the QTL regions, six QTLs on chromosomes 2, 3 and 4 were promising (Fig. 2). By searching previous reported genes related to salinity tolerance, six candidate genes located at or near the peak SNPs were identified for five promising QTLs (Table 4). There are 1–5 SNPs within these genes, these SNP markers could be used for validate the candidate genes.

Fig. 2

Manhattan plots of promising QTLs. Manhattan plots of selected QTLs were generated by using GLM with 1000 times of permutation

Table 4 Candidate genes for promising QTLs


It is difficult to introduce salinity tolerance of indica accessions into japonica accessions that have high grain quality and high yield [41]. Inter sub-specific crossing barriers between indica and japonica background also limit the ready introgression of salinity tolerance QTLs [42]. Breeding methods have been limited by the lack of salinity-tolerant japonica accessions. To date, very few salinity-tolerant rice accessions have been identified and used in breeding programs, and of those which have been identified, most are indica accessions. In this study, we identified some salinity-tolerant temperate japonica rice accessions such as, Shinchiku Iku 97, Lomellino, Norin 21, 81A32 and 36,037–1, showing good tolerance at both EC12 and EC18. These accessions could be important sources of salinity tolerance for future temperate rice breeding programs.

Many salinity tolerance QTLs have been identified in different studies [11,12,13,14,15,16,17,18,19,20,21,22,23,24]. However, most of them were mapped using different populations and markers. It is difficult to compare the positions of those QTLs and the QTLs identified in our current study because different markers were used. Fortunately, many genes related to rice salinity tolerance have been published in the QTL Annotation Rice Online database [43], and those genes can be tracked on the genome browser of rice SNP-seek database [34]. We compared the positions of the QTLs identified in this study and the salinity tolerance related genes in the genome browser. The results showed many QTL/gene clusters on chromosomes 1, 2, 3, 4, 6, 8, 9 and 11. Some of the QTLs identified in this study were co-located with salinity tolerance genes on chromosomes 1, 2, 3, 6, 8 and 9 (Additional file 9: Figure S5). These QTL/gene clusters might be important loci for salinity tolerance in temperate japonica rice varies, and related SNP markers could be used for QTL validation and marker assisted selection in breeding programs.

By comparing the SNPs in the QTL regions, six QTLs with a large number of SNPs showing low p values were identified as very promising (Fig. 2). Six candidate genes located at or near the peak SNPs were identified for five QTLs (Table 4). In the QTL region of qSDW0_2.1, there is no known gene controlling shoot dry weight, only a homeobox domain containing protein gene (Os02g13310) was down-regulated under desiccation and salt stress [44], however, no SNP was found in this gene.

In the QTL region of qPD18_2.1 and qSES18_2.1, gene Os02g40530 is an MYB transcription factor which was up-regulated under drought and/or salt stresses [45]. Another gene Os02g40730 is a Ammonium transporter (OsAMT1.3) which was significantly up-regulated under salt treatment of rice varieties FL478 and IR29, since high Cl concentrations can inhibit NO3 transport, up-regulation of ammonium uptake would be beneficial to withstand salt stress [46].

In the QTL region of qSKD_3.1 and qSNa/K0_3.1, gene Os03g07450 is a homeobox associated leucine zipper gene, which was down-regulated under desiccation and salt stress [44]. Gene Os03g07480 encodes a sucrose transporter (SUT). Sucrose allocation between tissues is a fundamental process in all multicellular organisms. Thus, sucrose transporter genes play an essential role in phloem loading and assimilate partitioning [47]. In sweet sorghum, salinity-tolerant species accumulate more sucrose by enhancing the synthesis and reducing the decomposition of sucrose under salt stress, while salt-sensitive species show enhanced decomposition of sucrose [48]. There are nine sucrose transporter genes in the rice genome [49, 50]. Sucrose transporter genes are essential for long distance sucrose transport and for osmo-protectant activities during drought and salinity stresses. Cultivars with higher expression of SUTs should be able to tolerate drought and salinity stresses better than those with lower expression [51].

In the QTL region of qSNa12_4.1, genes Os04g33640 and Os04g33720 belong to gene family of glycosyl hydrolases (GH). Glycosyl hydrolase functions in both the biosynthetic and hydrolytic pathways of raffinose metabolism, especially under certain abiotic stress conditions such as drought, high salinity or high temperature [52, 53]. Gene Os04g33640 (Glycosyl hydrolase family 17) showed the highest alternate transcripts among 579 genes in a cold-stressed leaf library [54]. Os04g33720 was specifically expressed in lowland rice, but not expressed in upland rice [55]. Os04g33720 was also shown to play an important signalling role in leaves subject to biotic and abiotic stresses [56, 57].

Although previous studies showed that some of the candidate genes were directly or indirectly related to salinity tolerance, the function of many of these genes are still unknown. Further study is needed to validate the expression and effect of these genes under salt stress.

Among the previously reported salinity tolerance QTLs, a major QTL for shoot K+ concentration (qSKC1 or O S HKT1;5) explaining 40.1% of the total phenotypic variance was identified on chromosome 1 and subsequently investigated by map-based cloning. It was shown to encode a HTK1 Na+ transporter and is involved in Na+ and K+ homeostasis [12, 58, 59]. Another QTL controlling shoot Na+/K+ ratio, SalTol, was originally identified using a RIL population derived from the cross IR29/Pokkali, explained 64% of the phenotypic variation [9], and was validated to the same region of SKC1 [11]. Haq et al. [14] also identified a QTL for Na+/K+ ratio between 11.1 and 14.6 Mb on chromosome 1 from a tropical japonica variety Moroberekan. A recent association mapping effort using varieties from the japonica cultivar group has also identified the SalTol genomic region as controlling important aspects of salinity tolerance [26]. It has been suggested that O S HKT1;5 may be the causal gene underlying these QTLs and functionally effective for salt tolerance across both indica and japonica accessions. However, Platten et al. (2013) identified seven major and three minor alleles of O S HKT1;5 gene from Oryza sativa and Oryza glaberrima, the aromatic allele conferred the highest leaf Na+ exclusion, and the japonica allele the least [60]. A salinity-tolerant allele of the O S HKT1;5 gene was also found in the wild rice species O. rufipogon and O. nivara [61]. In our current study, no QTL for Na+/K+ ratio was identified in the chromosomal region of these QTLs/genes. There are 6 SNPs in the O S HKT1;5 gene Oso1g20160, and 2 of them were polymorphic between the salinity-tolerant indica rice varieties in which qSKC1 and SalTol were identified, Nona Bokra and Pokkali. The haplotype of the 176 temperate japonica rice accessions were not the same as Nona Bokna and Pokkali, suggesting that the allele of O S HKT1;5 gene was not functional in the 176 temperate japonica rice accessions used in this study (Additional file 10: Table S5). A previous study showed that SNP mutations may cause the destabilization of a transmembrane domain in O S HKT1;5 and increase the probability of the O S HKT1;5 phosphorylation [62]. Rice lines containing Leu in position 395 of O S HKT1;5 gene exhibited higher shoot Na+ concentration than those containing Val [63]. The change of Asn to Asp in HKT1-type transporters also established altered cation selectivity and uptake dynamics, which enabled some crucifer species to acquire improved salt tolerance [64]. Thus, the SNP haplotype variation in the O S HKT1;5 alleles may render this gene non-functional in the temperate japonica accessions used in this study.

Under moderate salinity stress (EC12), some temperate japonica rice accessions achieved the same level of salinity tolerance as the check variety FL478, despite not having the same alleles at the SKC1 loci associated with Na+/K+ ratio. This suggests differences between indica and japonica subspecies in the effect of QTLs and genes involved in salinity tolerance. It is important to identify those temperate japonica specific salinity tolerance QTLs/genes that can significantly improve the salinity tolerance of future temperate japonica rice breeding lines.


In this study, we evaluated the seedling stage salinity tolerance of a subset of sequenced temperate japonica rice from 3000 sequenced rice accessions, and did genome-wide association study for 16 salinity responsive traits. Twenty-seven salinity-tolerant accessions were identified, and 26 QTLs related to 9 salinity tolerance traits were mapped on chromosomes 1, 2, 3, 4, 6, 8, 9 and 11. Six candidate genes were identified for five promising QTLs. The salinity tolerant varieties and QTLs could be used in related breeding programs. This study provided useful information for future studies on genetics and breeding of salinity tolerance in temperate japonica rice.



Genome-wide association study


Percentage of damage


Quantitative trait locus


Quantitative trait loci


Shoot dry weight


Standard evaluation score


Single nucleotide polymorphism


  1. 1.

    Zeng L, Shannon M. Salinity effects on the seedling growth and yield components of rice. Crop Sci. 2000;40:996–1003.

    Article  Google Scholar 

  2. 2.

    Rao P, Mishra B, Gupta S, Rathore A. Reproductive stage tolerance to salinity and alkalinity stresses in rice genotypes. Plant Breed. 2008;127:256–61.

    Article  Google Scholar 

  3. 3.

    Maas E. Agricultural salinity assessment and management. In: Tanji K, editor. Asce manualsandreportsonengineering. New York: American Society of Civil Engineers; 1990. p. 262–304.

    Google Scholar 

  4. 4.

    ESCAP. State of the environment in asia and pacific. New York: United Nations; 1995.

    Google Scholar 

  5. 5.

    Mamonov L, Taranov O, Sedlovskiy A. Biotechnological approaches for rice plant breeding in republic of kazakhstan. Biotechnol Theory Pract. 2010;2:62–6.

    Google Scholar 

  6. 6.

    Gregorio G, Singh R, Islam R, Sajise A, Vergara G. Salt-tolerant varieties developed to increase Rice self-sufficiency in salt-affected farmlands. In: 4th international rice congress: 27 October-1 November 2014; Bangkok, Thailand. Los Banos: International Rice Research Institute; 2014.

    Google Scholar 

  7. 7.

    Negrão S, Courtois B, Ahmadi N, Abreu I, Saibo N, Oliveira M. Recent updates on salinity stress in rice: from physiological to molecular responses. Crit Rev Plant Sci. 2011;30(4):329–77.

    Article  Google Scholar 

  8. 8.

    Munns R, Tester M. Mechanisms of salinity tolerance. Annu Rev Plant Biol. 2008;59:651–81.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Gregorio G, Senadhira D, Mendoza R. In: Institute IRR, editor. Screening rice for salinity tolerance. Manila: IRRI; 1997.

    Google Scholar 

  10. 10.

    Lee KS, Choi WY, Ko JC, Kim TS, Gregorio GB. Salinity tolerance of japonica and indica rice (oryza sativa l.) at the seedling stage. Planta. 2003;216(6):1043–6.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Bonilla P, Dvorak J, Mackill D, Deal K, Gregorio G. Rflp and sslp mapping of salinity tolerance genes in chromosome 1 of rice (oryza sativa l.) using recombinant inbred lines. Philipp Agric Sci. 2002;85:68–76.

    Google Scholar 

  12. 12.

    Lin H, Zhu M, Yano M, Gao J, Liang Z, Su W, Hu X, Ren Z, Chao D. Qtls for na+ and k+ uptake of the shoots and roots controlling rice salt tolerance. Theor Appl Genet. 2004;108(2):253–60.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Thomson M, Ocampo MD, Egdane J, Rahman M, Sajise A, Adorada D, Tumimbang-Raiz E, Blumwald E, Seraj Z, Singh R, et al. Characterizing the saltol quantitative trait locus for salinity tolerance in rice. Rice. 2010;3:148–60.

    Article  Google Scholar 

  14. 14.

    Haq U, Gorham J, Akhtar J, Akhtar N, Steele K. Dynamic quantitative trait loci for salt stress components on chromosome 1 of rice. Funct Plant Biol. 2010;37(7):634–45.

    Article  Google Scholar 

  15. 15.

    Lee S, Ahn J, Cha Y, Yun D, Lee M, Ko J, Lee K, Eun M. Mapping qtls related to salinity tolerance of rice at the young seedling stage. Plant Breed. 2007;126:43–6.

    Article  Google Scholar 

  16. 16.

    Prasad S, Bagali P, Hittalmani S, Shashidhar H. Molecular mapping of quantitative trait loci associated with seedling tolerance to salt stress in rice (oryza sativa l.). Curr Sci. 2000;78:162–4.

    CAS  Google Scholar 

  17. 17.

    Sabouri H, Rezai A, Moumeni A, Kavousi A, Katouzi M, Sabouri A. Qtls mapping of physiological traits related to salt tolerance in young rice seedlings. Biol Plant. 2009;53(4):657–62.

    Article  Google Scholar 

  18. 18.

    Takehisa H, Shimodate T, Fukuta Y, Ueda T, Yano M, Yamaya T, Kameya T, Sato T. Identification of quantitative trait loci for plant growth of rice in paddy field flooded with salt water. Field Crop Res. 2004;89:85–95.

    Article  Google Scholar 

  19. 19.

    Koyama M, Levesley A, Koebner R, Flowers T, Yeo A. Quantitative trait loci for component physiological traits determining salt tolerance in rice. Plant Physiol. 2001;125:406–22.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Islam M, Salam M, Hassan L, Collard B, Singh R, Gregorio G. Qtl mapping for salinity tolerance at seedling stage in rice. Emir J Food Agric. 2011;23:137–46.

    Article  Google Scholar 

  21. 21.

    Ammar M, Pandit A, Singh R, Sameena S, Chauhan M, Singh A, Sharma P, Gaikwad K, Sharma T, Mohapatra T, et al. Mapping of qtls controlling na+, k+ and cl- ion concentrations in salt tolerant indica rice variety csr27. J Plant Biochem Biotechnol. 2009;18:139–50.

    CAS  Article  Google Scholar 

  22. 22.

    Cheng L, Wang Y, Meng L, Hu X, Cui Y, Sun Y, Zhu L, Ali J, Xu J, Li Z. Identification of salt-tolerant qtls with strong genetic background effect using two sets of reciprocal introgression lines in rice. Genome. 2011;5(55):45–55.

    Google Scholar 

  23. 23.

    Pandit A, Rai V, Bal S, Sinha S, Kumar V, Chauhan M, Gautam R, Singh R, Sharma P, Singh A. Combining qtl mapping and transcriptome profiling of bulked rils for identification of functional polymorphism for salt tolerance genes in rice. Mol Gen Genomics. 2010;284:121–36.

    CAS  Article  Google Scholar 

  24. 24.

    Wang Z, Chen Z, Cheng J, Lai Y, Wang J, Bao Y, Huang J, Zhang H. Qtl analysis of na? And k? Concentrations in roots and shoots under different levels of nacl stress in rice. PLoS One. 2012;7:e51202.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Zheng H, Zhao H, Liu H, Wang J, Zou D. Qtl analysis of na+ and k+ concentrations in shoots and roots under nacl stress based on linkage and association analysis in japonica rice. Euphytica. 2015;201:109–21.

    CAS  Article  Google Scholar 

  26. 26.

    Ahmadi N, Negrao S, Katsantonis D, Frouin J, Ploux J, Letourmy P, Droc G, Babo P, Trindade H, Bruschi G, et al. Targeted association analysis identified japonica rice varieties achieving na+/k+ homeostasis without the allelic make-up of the salt tolerant indica variety nona bokra. Theor Appl Genet. 2011;123(6):881–95.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Zhao K, Tung C, Eizenga G, Wright M, Ali M. Genome-wide association mapping reveals a rich genetic architecture of complex traits in oryza sativa. Nat Commun. 2011;2:467.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Huang X, Wei X, Sang T, Zhao Q, Feng Q, Zhao Y, Li C, Zhu C, Lu T, Zhang Z, et al. Genomewide studies of 14 agronomic traits in rice landraces. Nat Genet. 2010;42(11):961–7.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Huang X, Zhao Y, Wei X, Li C, Wang A, Zhao Q, Li W, Guo Y, Deng L, Zhu C, et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2012;44:32–42.

    Article  Google Scholar 

  30. 30.

    Li J, Wang J, Zeigler R. The 3,000 rice genomes project: new opportunities and challenges for future rice research. GigaScience. 2014;3:8.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    The3KRGP. The 3000 rice genomes project. GigaScience. 2014;3:7.

    Article  Google Scholar 

  32. 32.

    IRRI. Standard evaluation system for rice. 4th ed. Manila: IRRI; 1996.

    Google Scholar 

  33. 33.

    Alexandrov N, Tai S, Wang W, Mansueto L, Palis K, Fuentes RR, Ulat VJ, Chebotarov D, Zhang G, Li Z, et al. Snp-seek database of snps derived from 3000 rice genomes. Nucleic Acids Res. 2015;43:1023–7.

    Article  Google Scholar 

  34. 34.

    Rice snp seek database. International rice research institute. In: http://oryzasnporg/iric-portal/indexzul. Accessed 1 May 2017.

  35. 35.

    Bradbury P, Zhang Z, Kroon D, Casstevens T, Ramdoss Y, Buckler E. Tassel: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23:2633–5.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Endelman J, Jannink J. Shrinkage estimation of the realized relationship matrix. G3. 2012;2:1405–13.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Zhang Z, Ersoz E, Lai C, Todhunter R, Tiwari H, Gore M, Bradbury P, Yu J, Arnett D, Ordovas J, et al. Mixed linear model approach adapted for genome-wide association studies. Nat Genet. 2010;42:355–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Ye C, Tenorio F, Argayoso M, Laza M, Koh H, Redona E, Jagadish K, Gregorio G. Identifying and confirming quantitative trait loci associated with heat tolerance at flowering stage in different rice populations. BMC Genet. 2015;16:41.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    McCouch S. Gene nomenclature system for rice. Rice. 2008;1:72–84.

    Article  Google Scholar 

  40. 40.

    Voorrips R. Mapchart: software for the graphical presentation of linkage maps and qtls. J Hered. 2002;93(1):77–8.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Yen C, Lin J. Screening, inheritance and linkage marker analyses of salt tolerance in mutated scented japonica rice. Plant Prod Sci. 2011;14(3):260–9.

    CAS  Article  Google Scholar 

  42. 42.

    Chen J, Ding J, Ouyang Y, Du H, Yang J, Cheng K, Zhao J, Qiu S, Zhang X, Yao J, et al. A triallelic system of s5 is a major regulator of the reproductive barrier and compatibility of indica–japonica hybrids in rice. PNAS. 2008;105(32):11436–41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Qtl annotation rice online database. National institute of agrobiological sciences. In: http://qtaroabraffrcgojp/. Accessed 1 May 2017.

  44. 44.

    Jain M, Tyagi A, Khurana J. Genome-wide identification, classification, evolutionary expansion and expression analyses of homeobox genes in rice. FEBS J. 2008;275:2845–61.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Chantarachot T, Buaboocha T, Gu H, Chadchawan S. Putative calmodulin-binding r2r3-myb transcription factors in rice. Thai J Bot. 2012;4:101–12.

    Google Scholar 

  46. 46.

    Senadheera P, Singh R, Maathuis FM. Differentially expressed membrane transporters in rice roots may contribute to cultivar dependent salt tolerance. J Exp Bot. 2009;60:2553–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Siahpoosh M. Sugar partitioning and sink-source modification in plants. J Rice Res. 2014;2:e106.

    Google Scholar 

  48. 48.

    Sui N, Yang Z, Liu M, Wang B. Identification and transcriptomic profiling of genes involved in increasing sugar content during salt stress in sweet sorghum leaves. BMC Genomics. 2015;16:534.

    Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Aoki N, Hirose T, Cofield G, Whitfeld P, Furbank R. The sucrose transporter gene family in rice. Plant Cell Physiol. 2003;44:223–32.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Shiratake K. Genetics of sucrose transporter in plants. Genes Genomes Genomics. 2010;1:73–80.

    Google Scholar 

  51. 51.

    Ibraheem O, Dealtry G, Roux S, Bradley G. The effect of drought and salinity on the expressional levels of sucrose transporters in rice cultivar plants. Plant Omics J. 2011;4:68–74.

    CAS  Google Scholar 

  52. 52.

    Nishizawa A, Yabuta Y, Shigeoka S. Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 2008;147:1251–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Maruyama K, Takeda M, Kidokoro S, Yamada K, Sakuma Y, Urano K, Fujita M, Yoshiwara K, Matsukura S, Morishita Y. Metabolic pathways involved in cold acclimation identified by integrated analysis of metabolites and transcripts regulated by dreb1a and dreb2a. Plant Physiol. 2009;150:1972–80.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Venu R, Sreerekha M, Sheshu M, Nobuta K, Mohan K, Chen S, Jia Y, Meyers B, Wang G. Deep transcriptome sequencing reveals the expression of key functional and regulatory genes involved in the abiotic stress signaling pathways in rice. J Plant Biol. 2013;56:216–31.

    CAS  Article  Google Scholar 

  55. 55.

    Gao F, Zhang H, Wang H, Gao H, LI Z. Comparative transcriptional profiling under drought stress between upland and lowland rice using cdna-aflp. Chin Sci Bull. 2009;54:3555–71.

    CAS  Google Scholar 

  56. 56.

    Roitsch T, Balibrea M, Hofmann M, Proels R, Sinha A. Extracellular invertase: key metabolic enzyme and pr protein. J Exp Bot. 2003;54:513–24.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Roitsch T, Gonzalez M. Function and regulation of plant invertases: sweet sensations. Trends Plant Sci. 2004;9:606–13.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Ren Z, Gao J, Li L, Cai X, Huang W, Chao D, Zhu M, Wang Z, Luan S, Lin H. A rice quantitative trait locus for salt tolerance encodes a sodium transporter. Nat Genet. 2005;37:1141–6.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Platten D, Cotsaftis O, Berthomieu P, Bohnert H, Davenport R, Fairbairn D, Horie T, Leigh R, Lin H, Luan S, et al. Nomenclature for hkt transporters, key determinants of plant salinity tolerance. Trends Plant Sci. 2006;11:372–4.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Platten J, Egdane J, Ismail A. Salinity tolerance, na+exclusion and allele mining of hkt1; 5 in oryza sativa and o. Glaberrima: many sources, many genes, one mechanism? BMC Plant Biol. 2013;13:32.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Mishra S, Singh B, Panda K, Singh B, Singh N, Misra P, Rai V, Singh N. Association of snp haplotypes of hkt family genes with salt tolerance in indian wild rice germplasm. Rice. 2016;9:15.

    Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Negrão S, Almadanim C, Pires I, Abreu I, Maroco J, Courtois B, Gregorio G, McNally K, Oliveira M. New allelic variants found in key rice salt-tolerance genes: an association study. Plant Biotechnol J. 2013;11:87–100.

    Article  PubMed  Google Scholar 

  63. 63.

    Cotsaftis O, Plett D, Shirley N, Tester M, Hrmova M. A two-staged model of na+ exclusion in rice explained by 3d modeling of hkt transporters and alternative splicing. PLoS One. 2012;7:e39865.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Ali A, Raddatz N, Aman R, Kim S, Park H, Jan M, Baek D, Khan I, Oh D, Lee S, et al. A single amino-acid substitution in the sodium transporter hkt1 associated with plant salt tolerance. Plant Physiol. 2016;171:2112–26.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


This study was partially supported by Kazakhstan government project 2113/GF4 “Biotechnology and development of high-yield and salt-tolerant rice breeding lines”, the Bill & Melinda Gates Foundation project Cereal Systems Initiative for South Asia (CSISA), and the high-end talent development program (2016HE006) of Yunnan province (China).


Kazakhstan government project 2113/GF4 and the high-end talent development program (2016HE006) of Yunnan province (China).

Availability of data and materials

The sequence and SNP data of all accessions are available in the SNP-Seek system ( The seeds of all the used accessions are available in the International Rice genebank ( List of temperate japonica accessions used are available in Additional file 2: Table S1. Manhattan plots of all traits are available in Additional file 8: Figure S3.

Author information




DB, CY, BU, GD, GV, RR and HL designed the experiments. DB, BL and CY conducted the experiments. CY, DB, and LX analyzed the data and wrote the paper. All authors read and approved the final manuscript for publishing.

Corresponding author

Correspondence to Changrong Ye.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Figure S1. Selection of 203 temperate japonica from the 3000 sequenced rice accessions. Left is the archaeopteryx tree of 3000 sequenced rice accessions using 365,710 SNPs, and the right is the archaeopteryx tree of 176 temperate japonica accessions using 68,786 SNPs. (DOCX 133 kb)

Additional file 2:

Table S1. List of temperate japonica accessions used. (DOCX 28 kb)

Additional file 3:

Table S2. Standard evaluation score (SES) of visual salt injury at seedling stage. (DOCX 12 kb)

Additional file 4:

Table S3. Statistics of the measured traits. (DOCX 14 kb)

Additional file 5:

Figure S2. QQ plot of the measured traits. (DOCX 57 kb)

Additional file 6:

Figure S3. Archaeopteryx tree of 176 temperate japonica rice accessions using 68,786 SNPs. The origin country of the accession was shown after the sequence ID. (DOCX 697 kb)

Additional file 7:

Table S4. List of salinity tolerant accessions at seedling stage based on SES12. (DOCX 16 kb)

Additional file 8:

Figure S4. Manhattan plots of different traits related to salinity tolerance at seedling stage. The genotyping and phenotyping data along with kinship matrix were analyzed using mixed linear model (MLM). (DOCX 1043 kb)

Additional file 9:

Figure S5. Distribution of identified QTLs on chromosomes. The Os genes are genes related to salinity tolerance from QTARO database ( and rice SNP-seek database ( The number on the left is the position of the QTL/gene in Mb. (DOCX 711 kb)

Additional file 10:

Table S5. Hyplotype of SKC1 gene (Oso1g20160, OsHKT1 Na+ transporter) in salinity-tolerant varieties (Nona Bokra and Pokkali) and 191 temperate japonica accessions. (DOCX 35 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Batayeva, D., Labaco, B., Ye, C. et al. Genome-wide association study of seedling stage salinity tolerance in temperate japonica rice germplasm. BMC Genet 19, 2 (2018).

Download citation


  • Genome-wide association studies
  • Seedling stage
  • Salinity tolerance
  • Candidate gene
  • Temperate japonica