Next Article in Journal
Phytohormones Regulate the Development of Arbuscular Mycorrhizal Symbiosis
Next Article in Special Issue
Identification of Salt Stress Responding Genes Using Transcriptome Analysis in Green Alga Chlamydomonas reinhardtii
Previous Article in Journal
Effects of Panax Notoginseng Saponins on Esterases Responsible for Aspirin Hydrolysis In Vitro
Previous Article in Special Issue
Overexpression of a Novel ROP Gene from the Banana (MaROP5g) Confers Increased Salt Stress Tolerance
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 10
10.3390/ijms19103145
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Genome-Wide Association Study Reveals Candidate Genes Related to Salt Tolerance in Rice (Oryza sativa) at the Germination Stage

by
Jie Yu
1,2,†,
Weiguo Zhao
1,3,†,
Wei Tong
1,2,
Qiang He
1,4,
Min-Young Yoon
1,5,
Feng-Peng Li
1,6,
Buung Choi
1,7,
Eun-Beom Heo
1,8,
Kyu-Won Kim
9,* and
Yong-Jin Park
1,9,*
1
Department of Plant Resources, College of Industrial Sciences, Kongju National University, Yesan 32439, Korea
2
State Key Laboratory of Tea Plant Biology and Utilization, Anhui Agricultural University, Hefei 230036, China
3
School of Biotechnology, Jiangsu University of Science and Technology, Sibaidu, Zhenjiang, Jiangsu 212018, China
4
National Key Facility for Crop Resources and Genetic Improvement, Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing 100081, China
5
Leader of Eco. Energy & Bio (LEEBCOR), 190-26 Hwangyeonggongwon-ro, Asan-si, Chungcheongnam-do 31529, Korea
6
Suzhou GENEWIZ Biotechnology Co. LTD, C3 218 Xinghu Road Suzhou Industrial Park, Suzhou 215123, China
7
Chemical Safety Division, National Institute of Agricultural Sciences (NIAS), Wanju 55365, Korea
8
Breeding & Research Institute, Koregon Co. LTD, Anseong Center 60-34, Gokcheon-gil, Bogae-Myeon, Anseong-Si, Gyeonggi-Do 17509, Korea
9
Center of Crop Breeding on Omics and Artificial Intelligence, Kongju National University, Yesan 32439, Korea
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2018, 19(10), 3145; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19103145
Submission received: 21 September 2018 / Revised: 9 October 2018 / Accepted: 9 October 2018 / Published: 12 October 2018
(This article belongs to the Special Issue Salinity Tolerance in Plants)
Browse Figures
Versions Notes

Abstract

:
Salt toxicity is the major factor limiting crop productivity in saline soils. In this paper, 295 accessions including a heuristic core set (137 accessions) and 158 bred varieties were re-sequenced and ~1.65 million SNPs/indels were used to perform a genome-wide association study (GWAS) of salt-tolerance-related phenotypes in rice during the germination stage. A total of 12 associated peaks distributed on seven chromosomes using a compressed mixed linear model were detected. Determined by linkage disequilibrium (LD) blocks analysis, we finally obtained a total of 79 candidate genes. By detecting the highly associated variations located inside the genic region that overlapped with the results of LD block analysis, we characterized 17 genes that may contribute to salt tolerance during the seed germination stage. At the same time, we conducted a haplotype analysis of the genes with functional variations together with phenotypic correlation and orthologous sequence analyses. Among these genes, OsMADS31, which is a MADS-box family transcription factor, had a down-regulated expression under the salt condition and it was predicted to be involved in the salt tolerance at the rice germination stage. Our study revealed some novel candidate genes and their substantial natural variations in the rice genome at the germination stage. The GWAS in rice at the germination stage would provide important resources for molecular breeding and functional analysis of the salt tolerance during rice germination.

1. Introduction

Salt stress undesirably affects plant growth during all developmental stages. Therefore, it is a major threat to crop productivity [1] and this situation has lasted in some parts of the world for over 3000 years with growth every year. As a monocotyledonous model plants, rice feeds more than one half of the world’s population [2]. However, rice is also sensitive to salt stress and is currently listed as the most salt-sensitive cereal crop, which results in most cultivated varieties having a salinity threshold of 3 dSm−1 [3].
Seed germination is usually a very important stage in the seedling stable stand establishment and determines the success of crop production [4]. The effects of salt stress on seed germination are extremely complex involving various physical and biochemical cues. Generally, salt stress is negatively correlated with seed germination and seedling growth [5] in most plants such as Oryza sativa [6], Zea mays [7], Helianthus annuus [8], and Brassica spp. [9]. The impact of salt stress on seed germination is attributed to seed water uptake and ion toxic effect. During the seed germination, salinity alters the imbibition of water by reducing the osmotic potential of the germination medium [10], damages the ultrastructure of cells, tissue, and organs [11], changes the activity of enzymes [12], disturbs hormonal balance [13], alters protein metabolism [14], and reduces the use of seed reserves [15]. However, various environmental (external) and plant physiological (internal) factors affect seed germination under saline conditions including temperature, light, water and gasses, seed age, seed dormancy, nature of seed coat, seed morphology, and seedling vigor [16].
Recently, improving rice salt tolerance during the germination stage became more important because salinity may rapidly reduce the germination rate and percentage, which, in turn, may lead to a reduction of crop yields [17]. Many efforts have been made to improve seed germination and seedling vigor by optimizing the non-genetic factors [6]. However, success in improving salt tolerance in rice is made by identifying a major quantitative trait locus (QTL), which contributes to salt tolerance in rice. QTL analysis of seed germination have been reported in rice [18], soybean [19], wheat [20], Arabidopsis [21], and Brassica rapa [22]. However, it is difficult to develop rice elite varieties with a high level of salt tolerance due to a lack of understanding the mechanisms of salt tolerance during the seed germination stage. Moreover, QTLs conferring salt stress tolerance in rice were identified mainly at the seedling stage [23], but there are few reports on rice seed germination [18,24,25]. Wang et al. [18] detected 16 QTLs for the imbibition rate and germination percentage under control and salt stress with the recombinant inbred line (RIL) population derived from IR26/Jiucaiqing. Abe et al. [26] identified a candidate gene, OsGA20ox1, for a major QTL controlling seedling vigor in rice. Zheng et al. [27] identified 11 QTLs for salt tolerance at the germination and early seedling stage in japonica rice.
Genome-wide association study (GWAS) is an efficient method for detecting valuable natural variations in trait-associated loci as well as allelic variations in candidate genes underlying quantitative and complex traits [28,29]. Instead of SSR (simple sequence repeat) markers, which are commonly used for association mapping [30,31], SNPs (single-nucleotide polymorphisms) have become more popular for GWASs with the rapid development of NGS (next-generation sequencing) and the existence of high-density SNP markers by re-sequencing [32]. In rice, there are some successful reports to dissect genetic architecture of complex traits through GWAS [28,29,32,33]. However, limited studies have been carried out in rice to identify genes/QTLs for salt tolerance using GWAS. Kumar et al. [34] identified total 64 SNPs significantly associated with Na+/K+ ratio and other traits for reproductive stage salinity tolerance using GWAS. Yu et al. [35] identified 93 candidate genes significantly associated with salt tolerance at the rice seedling stage. Shi et al. [36] identified 22 significant salt tolerance associated SNPs based on the stress-susceptibility indices (SSIs) of vigor index (VI) and the mean germination time (MGT). Naveed et al. [37] identified 20 QTN for salinity tolerance at the germination and seedling stages in rice. Seed germination plays an important role in the cycle of plant growth. To our knowledge, there is little research until now on the identification of genes particularly for the germination stage salinity tolerance using GWAS in rice. Here, we applied GWAS mapping using ~1.65 million SNPs/indels covering all 12 rice chromosomes in a diverse rice collection to identify candidate genes and natural variation that may contribute to salt tolerance during the rice germination stage with the aim to guide breeding of salt-tolerant rice varieties.

2. Results

2.1. Phenotypic Screening and Evaluation

Individual value plots for GP with 0, 200 and 300 mM NaCl from a screening experiment using 12 randomly selected samples are shown in Figure 1a. According to the screening result, 200 mM NaCl fully exhibited their phenotypic variance, which resulted in the most diverse phenotypic distribution and facilitated discrimination of accessions with different salt tolerance levels. Thus, the treatment of 200 mM NaCl was chosen as the target salinity level for determining the salt tolerance of all accessions.
The following traits: GP, GE, GI, SL, and RL were examined under 0 and 200 mM NaCl salt stress during the rice germination stage. Descriptive statistics of the phenotypes related to salt tolerance during the germination stage of the current collection were presented in Table 1. Box plots of phenotypes including GE, GI, SL, and RL in the presence of 0 and 200 mM NaCl are shown in Figure 1b–d. The data suggested that seed germination traits were negatively influenced by salt stress. Salt stress inhibits shoot and root elongation dramatically but GE and RL were more affected than SL [18]. These findings indicated that most germination parameters under salt stress exhibited lower performance than under control conditions, which may restrict plant growth.
The correlation coefficients of phenotypes under control and salt-stress conditions were also evaluated (Table 2). RL was significantly and positively correlated with SL only under control conditions. Excepting RL and SL was also significant positively correlated with GI. GP, GE, and GI were all significantly and positively correlated with each other under control conditions while all phenotypes were significantly and positively correlated with each other in salt stress conditions. These results suggested that all phenotypes evaluated in this study could be used for GWAS and some overlapped results could be found among the phenotypes.

2.2. Principal Components Analysis (PCA)

PCA was performed with the 1.65 million high-quality SNPs/indels to mine the population structure in all rice accessions. Two components were suggested by the scree plot (Supplementary file 2: Figure S1a). Clear subpopulation structures were observed based on the first two PCs (PC1 and PC2), which resulted in two subpopulations, indica and japonica, with the admixture accessions located between the two groups (Figure 2a).
For the PCA using the phenotypes, we examined correlations between subspecies in salt tolerance levels using four main phenotypes that drove the differences among accessions. We used TASSEL to perform a PCA of R-GE, R-GI, R-RL, and R-SL in the rice collection. Most of the phenotypic variation (>91%) in the collection was explained by the first two PCs (Supplementary file 1: Figure S1b). Thus, we generated a PCA plot using PC1 and PC2. However, rice accessions in our study were not clustered into clearly defined groups (such as indica or japonica) based on the above four phenotypes (Figure 2b). This indicates that salt tolerance levels in rice (O. sativa) are not strongly correlated with the indica or japonica subgroups.

2.3. GWAS and Candidate Gene Identification

To generate the genotype dataset for GWAS, more than ~1.65 million SNPs/indels were identified across the accessions and subjected to GWAS applied with the CMLM [38]. The GWAS results were shown on Figure 3 and Supplementary file 1: Figure S2. We took associations held by the peaks with −log10 (p) value > 5 and adjusted p-value (FDR, false discovery rate) < 0.05 for further analysis since the cutoff of −log10 (p) value was five when the FDR ≤ 0.05. Under the salt stress condition, the significant signals were detected for GE, RL, SL, R/S, and R-R/S. In total, 10 SNPs were found significant for these traits and only one SNP (chr12_1628276) were found common for RL, R/S, and R-R/S traits. Excluding the common one, only 3, 2, 2, 1, and 1 SNPs were found uniquely associated for R-R/S, GE, RL, SL, and R/S, respectively. In addition, we also found two SNPs (chr02_1090174 and chr05_20164893) were the two strongest, significantly associated for GE in all observed traits (p < 10−9).
We further conducted a genome-wide LD analysis of the candidate peak regions and determined LD blocks harboring significant SNPs/indels that characterized in the last step as regions containing putative candidate genes. LD block analysis was detected in a 400 kb range centered on the highest −log10 (p) value (Figure 3c,d). Annotation of SNPs/indels from the 200 kb up-stream and down-stream ranges, together with the LD block analysis, resulted in the identification of 79 genes included in these peaks and some candidate genes have been reported previously to contribute to salt tolerance (Supplementary file 2: Table S2). Among the known genes, seven were associated with chromosome 1 (all in RL), 22 with chromosome 2 (7 in GE, 6 in R/S and 9 in R-R/S), 10 with chromosome 3 with R-R/S, 15 with chromosome 4 (3 in SL and 12 in RL), 6 with chromosome 5 (all in GE), 6 with chromosome 11 (all in R-R/S), and 13 with chromosome 12 (the peak region was identical in RL, R/S, and R-R/S). In the LD block analysis, most highly associated SNPs/indels were located in small or large LD blocks, which indicated that they were in significant linkage disequilibrium. Thus, these candidate genes may contribute to salt stress independently or co-operatively with other variations in other genes harboring these SNPs/indels. Simultaneously, we screened candidate genes containing many highly associated SNPs/indels in the genic region as well as some highly associated signals not located in known genes but suggesting that these unknown genes may also be related to salt tolerance (Supplementary file 4: Table S3). Some of those SNPs/indels were located in the coding region of the unknown genes rather than in the surrounding 200 kb regions. These genes could also be important determinants of salt tolerance in rice.

2.4. Natural Variations in Candidate Genes and Sequence Analysis

Based on the associated peaks identified in the GWAS and by determining the LD blocks test, we identified many candidate genes associated with salt tolerance during the rice germination stage. To mine functional and novel candidate genes, we investigated 17 final candidate genes (Table 3) that contained highly associated SNPs/indels within the coding region. Many of these SNPs/indels have been reported to play a role in the salt stress in rice such as OsAGO2 (Os04g0615700) [39], OsZIFL13 (Os12g0133300) [40], and OsHAK11 (Os04g0613900) [41], which are related to salt stress in rice. These genes are involved in the salt tolerance in rice by different pathways [42].
Natural variations of these 17 genes were mined and then functional variations were screened after checking the positions of the variations in genes and the corresponding amino acid change. Among the genes, OsMADS31, which is involved in floral organ specification and implicated in plant growth and development, was identified and predicted to be involved in salt tolerance [43]. As shown in Figure 4, one natural SNP substitution (T/A) was detected and caused an F/L amino acid change, which is presented by type 1 (reference sequence) and type 2 (variation) (Figure 4a). Furthermore, we generated a haplotype network of the whole collection, which was dominated by two common haplotypes including primarily the japonica type (type 1) and the indica type (type 2), respectively (Figure 4e). A phenotypic difference was observed in type 1 with 236 accessions and an average RL of 0.6293 and type 2 with 58 accessions and an RL of 0.9327 (Figure 4b). We conducted further orthologue alignment of OsMADS31 in several rice groups and other species (Figure 4c). Type 2 (candidate SNP) showed an F/L amino acid change compared to other rice groups and species (type 1, Oryza brachyantha, Oryza rufipogon, Oryza punctate, Hordeum vulgare, Triticum aestivum, Aegilops tauschii, and Triticum urartu). However, type 2 shared this F/L with three other rice species (Oryza glaberrima, Oryza barthii and Oryza glumaepatula). Oryza glaberrima and Oryza barthii are African rice and its wild type have higher salt tolerance than Oryza sativa species. Oryza glumaepatula is a wild rice found in South America usually in deep and sometimes flowing water, which may also have salt tolerance characteristics based on the presence of related genes. Four salt tolerant accessions with type 2 haplotype and four salt sensitive accessions without the haplotype were used for the real-time expression analysis. Generally, the relative RNA expression level of OsMADS31 was higher in type 1 than in type 2, which indicates that the gene expression is down-regulated in salt conditions when compared to the control (Figure 4d).
We also found several other functional SNPs/indels in the 17 candidate genes that were correlated with a phenotypic difference (Supplementary file 5: Table S4). These candidate genes may be related to rice salt tolerance, according to both previous reports and the natural variation mining in the current study. Novel polymorphisms of those genes may also contribute to salt tolerance that make the rice resistant to salinity.

3. Discussion

3.1. Salt Tolerance at Rice Germination Stage

The seed germination is one of the most critical steps in the life cycle of a crop. Seed germination begins with water uptake while salinity prevents water imbibition, which inhibits seed germination [15]. Experiments have shown that increased salinity delays the initiation of germination, which leads to a reduced germination percentage. However, salt tolerance during the early growth stages is not always correlated during subsequent growth stages [44,45]. The seeds of crops in different genotypes may germinate adequately under salt stress. However, the seedling may not become fully established later. We observed differential inhibition of the root length and shoot length in our study, which suggested that salinity can influence the germination quality of the seed.
By using the optimized salinity (200 mM NaCl) for discriminating accessions with different salt tolerance levels, we characterized the salt-tolerance-related phenotypes in a collection comprising 295 rice varieties. Phenotypic differences between the control and 200 mM NaCl salinity conditions suggested that rice growth during the seed germination stage can be markedly inhibited by salt stress, which results in very low germination energy and index (GE and GI) as well as reduced root and shoot lengths. This may suppress rice seed germination especially in some direct-sowing areas and decreases plant density and yield markedly. Therefore, development of rice varieties with salt-tolerant seeds would prevent salinity-mediated plant and yield loss during the early growth stage.

3.2. Salt Tolerance Is Not Strongly Correlated with Rice Subgroups

According to Lee et al. [46], the salt tolerance of indica rice was higher than that of japonica rice at the seedling stage, which was determined by measuring shoot Na+ and K+ absorption. However, as revealed in a recent study of the salt tolerance of 115 O. sativa and O. glaberrima accessions, salt tolerance was not strongly correlated with O. sativa cultivar groups [35]. Most of the japonica types were salt sensitive, but accessions from the indica group and O. glaberrima showed a wide range of sensitivities [47]. In our study, we performed a PCA of all rice germplasm using both genotype and phenotype data. Inconsistent with the genotypic PCA, which separated the collection into clear groups, phenotypic PCA using germination-related phenotypes showed no clear grouping (Figure 2b). This indicated that salt tolerance levels during the seed germination stage are not well correlated with the rice (O. sativa) subgroup.

3.3. GWAS and Candidate Gene Identification

In some direct-sowing areas, salt tolerance in rice during the seed germination stage is particularly important. To improve rice productivity in such areas, novel genes and alleles associated with complex quantitative salinity tolerance traits must be identified in diverse rice accessions and salt-tolerant varieties bred. An alternate and complementary approach is GWAS, which takes advantage of historical recombination events and, thus, enables a high-resolution genome wide mapping for the identification of target genomic regions in response to complex quantitative traits in rice [29]. In this study, we used a core set of rice collections and multiple bred varieties to investigate candidate loci and genes that regulate important phenotypes under salt stress in rice at the germination stage. Twelve GWAS peaks representing new QTLs on chromosomes 1, 2, 3, 4, 5, 11, and 12 during the rice germination stage were identified. The current association mapping can serve as source of novel salt tolerance genes and alleles. Thus, we found abundant candidate regions with high association peaks in five traits and were distributed on seven chromosomes. Now many QTL analysis of rice salt tolerance have been reported, but it is difficult to directly compare the chromosomal location of marker–trait associations detected in this study with previously reported QTLs because different materials at different stages, descriptive traits, and molecular maps have been used. Wang et al. [18] detected 16 QTLs for the imbibition rate and the germination percentage. Kumar et al. [34] identified 64 SNPs (loci) significantly associated with salt stress-related traits by GWAS. Leon et al. [48] identified 85 additive QTLs for seedling salinity tolerance by GBS. Yu et al. [35] identified 25 SNPs (loci) significantly associated with salt stress-related traits by GWAS. Shi et al. [36] identified 22 SNPs based on SSIs of VI and MGT by GWAS. In this study, we also found that some SNPs associated with salt-tolerance traits overlapped or located in similar or proximal regions such as a significant SNP (chr04_31168058) near qRTL4.10 identified by Leon et al. [48] and SNPs (chr04_31164404) identified by Yu et al. (2017). This SNP is also located near the SNPs (chr04_34164920 and chr04_ 34292214) identified by Kumar et al. [34] associated with Na+/K+ ratio. Additionally, two QTLs for salt tolerance and potassium concentration were mapped just prior to this region, respectively, by Lin et al. [49] and Cai and Morishima [50]. The above results also indicated that chromosome 4 including many candidate genes in this region was found to be important for salt tolerance.
So far, about 70 salt tolerance QTLs had been located in rice using biparental mapping populations, but fine mapping and narrowing down reports are limited [34]. Driven by LD blocks to define the genomic regions for searching candidate genes has advantages over the fixed-window approach in which a certain distance from a significant SNP is considered to be the region containing candidate genes [51] by eliminating falsely included or excluded genes [52]. The wide candidate regions ranged from <1 kb to >1 Mb depending on the chromosomal position, which suggests that the resolution of the association mapping is highly dependent on the LD of the neighboring regions of the significant SNPs [34]. Since some of the LD blocks harboring significant SNPs did not contain an annotated gene, this method might have produced some false negatives or the identified region may have contained important DNA-binding or gene regulation sites, in which case, the causal gene was not detected in the LD block [53]. In this work, we used a 400 kb range of the strongest signal to locate the candidate genes, which is in line with previous studies [54]. From an LD block analysis, we obtained 79 candidate genes that had significant SNP/indel associations in LD block regions. Therefore, these regions and candidate genes have a statistically and genetically supported background and, therefore, may be important for the salt tolerance of rice during the germination stage. Apart from SNPs that had an association with previous known QTLs for salinity tolerance, there were a few SNPs, which hit specific genes that were known or functionally characterized for salt stress. Among 79 candidate genes including seven protein kinases (PK) (1 Serine/threonine protein kinase and 1 OsCDPK26), six ion exchanger and transporter related genes, five transcription factors (TFs), two electron carrier (peroxidase, Os01g0172600, and oxidoreductase, Os05g0411200), and two major facilitator superfamily proteins (Os12g0133100, OsZIFL12; Os12g0133300, OsZIFL13). In addition, we also found one stress-associated protein 18 (SAPs) (OsSAP18, Os02g0121600), one vacuolar ATPase assembly integral membrane protein (Os04g0612900), two argonaute family proteins (AGOs) (OsAGO2, Os04g0615700, and OsAGO3, Os04g0615800), one chloroplast precursor (Ferritin1, Os11g0106700), one calmodulin-like protein 3 (OsCML3, Os12g0132300), one Auxin efflux carrier protein (OsPIN1d, Os12g0133800), one Glycoside hydrolase (Os02g0532900), and one Glycosyltransferase (ALG3, Os01g0172000). The above results indicated that the candidate genes may play an essential role in salt tolerance mechanisms [34], which also indicated that salt tolerance genes are involved in ion pumps, calcium, the salt overly sensitive (SOS) pathway, mitogen-activated protein kinases (MAPK), glycine betaine, proline, and the reactive oxygen species pathways in a high salinity environment [42].
Furthermore, 17 candidate genes with high −log10 (p) value-associated signals inside the coding region were also mined and may play an important role in salt tolerance. OsHAKs are candidates for high-affinity K+ uptake transporters in the rice root. The transcription of OsHAK11 (Os04g0613900) is significantly induced by salt stress and K+ starvation, respectively [55]. AGOs (Os04g0615700) play important roles in the regulation of development and stress responses, antiviral immune response, transposons, and the regulation of chromatin structure and can affect the growth and development as well as the response to abiotic and biotic stress [56]. OsPIN (Os12g0133800), which encodes a member of the auxin efflux carrier proteins, is involved in the root elongation growth and lateral root formation patterns via the regulation of auxin distribution in rice [57]. The Germin family protein (Os02g0532500) had been revealed to be connected with a plant cell defense and diseases and to be highly resistant to sodium dodecyl sulfate (SDS) and proteases and important for early plant development and germination in plants [58]. SAP (Os12g0133700) is the A20/AN1 zinc-finger containing proteins, which can regulate the stress signaling in plants [59]. The Zinc-induced facilitator-like (ZIFL) family genes (Os12g0133100, OsZIFL12, Os12g0133300, OsZIFL13) are up-regulated under stress conditions [40].
Based on these regions and candidate genes, it may be possible to mine the natural variations of rice in response to salt stress in some tolerant accessions and apply those alleles to sensitive accessions via breeding methods. To the best of our knowledge, this is the first large-scale GWAS focusing on salt stress during the rice germination stage. These candidate regions and genes will facilitate the development of salt-tolerant rice varieties.

3.4. Novel Natural Variations of Candidate Genes

Investigation of new natural variations in focal traits can extend the tolerant varieties’ functional alleles to other non-tolerant varieties. Breeding methods can then be used to transfer them to elite lines to produce tolerant varieties. Using the results of the GWAS and LD analysis, the haplotypes of candidate genes can be targeted and the functional alleles involved in responses can be identified. New alleles in rice have been reported [60,61] and provide insight for researchers and breeders into the underlying mechanisms, which facilitates the breeding of improved varieties. According to Arora et al. [42], OsMADS31 expression was low and not markedly affected by salt and cold stress. However, the expression was relatively down-regulated in seedlings under a salt condition. Nevertheless, OsMADS31 expression was higher in seeds than during the panicle stage. In this study, we found that OsMADS31 was associated with salt tolerance in rice at the seed germination stage with a down-regulated expression in the salt condition (Figure 4d). The contribution of MADS-box genes to flower organ specification is well developed in eudicots, but not very well in rice. Therefore, the roles of MADS genes and other candidate genes identified here using GWAS at the seed germination stage in response to salt stress should be investigated further. Moreover, by adapting functional studies (such as those performed using TALEN and CRISPR/Cas 9), the functions of genes and gene variations can be determined. Natural variations that have functional signals could be a good starting point for the exploration of gene-based assays of phenotypically different individuals such as salt tolerant vs. sensitive, drought resistant vs. susceptible, and more.
Overall, we investigated the genetics architecture of natural variation in rice salt-tolerance-related traits at the germination stage by GWAS mapping in 295 rice accessions. A total of 79 candidate genes were determined by LD blocks analysis. In addition, by detecting the highly associated variations located inside the genic region that overlapped with the results of LD block analysis, we finally characterized 17 genes that may contribute to salt tolerance during the seed germination stage. The salt tolerance related novel candidate genes would provide important resources for molecular breeding and functional analysis of the salt tolerance during the rice germination.

4. Materials and Methods

4.1. Materials

A core set of 137 rice accessions and 158 bred varieties from the National Gene Bank of the Rural Development Administration (RDA-Genebank, Korea) [62,63] was re-sequenced in the current study (Supplementary file 1: Table S1). We conducted a field experiment during the rice-growing season at the Kongju National University experimental farm and young leaves from a single plant were collected and immediately kept at −80 °C prior to genomic DNA extraction using the DNeasy Plant Mini Kit (Qiagen, Hilden, Germany). Qualified DNA was sent for the whole genome re-sequencing.

4.2. Whole Genome Re-Sequencing and Variation Detection

The genomes of all 295 rice accessions were sequenced with an average coverage of approximately 7.8× on an Illumina HiSeq 2000 or 2500 Sequencing Systems Platform (Illumina Inc., San Diego, CA, USA). Raw reads were aligned against the rice reference genome (IRGSP 1.0) [64] for genotypes calling and only SNPs/indels without the missing value and a minor allele frequency (MAF) > 0.05 and containing genotype calls for all 295 accessions that were used. Lastly, ~1.65 million high-quality SNPs/Indels were obtained and used for the further GWAS [65].

4.3. Evaluation of Salt Stress and Phenotyping

We first carried out the pre-screening experiment using 12 randomly selected samples to determine the optimum level of NaCl concentration for the evaluation of salt stress during the germination stage. Seed germination were initially screened by germinating 30 seeds per genotype in petri dishes with two layers of filter papers soaked in two different NaCl concentrations: 200 and 300 mM NaCl. The germination percentage was recorded daily for 10 days. At the concentration of 300 mM NaCl, seeds hardly germinated and the seedlings did not grow out enough to be able to measure root and shoot length. Therefore, in this study, we used 0 mM NaCl (non-stress) and 200 mM NaCl (salt stress) for phenotyping all 295 accessions.
The following experiments were performed in petri dishes containing two-layered filter paper. Thirty seeds were first washed in water, then sterilized in 1% sodium hypochlorite solution for 10 min, and washed three times in deionized distilled water. Thereafter, seeds of each accession were soaked in petri dishes and then incubated at 30 °C with 40% relative humidity. Petri dishes were randomized in an incubator and three replicates of each accession under control and salt conditions (200 mM) were adopted. The solution was replaced every two days to maintain the NaCl concentration and the distilled water volume, respectively. The daily germination seed was measured and filter papers were replaced as necessary. Plumule emergence was taken as an index of germination. The length reached about 2 mm. At the end of day 10, we measured the root length (RL) and shoot length (SL) of the seedlings and the R/S (root/shoot ratio) was also calculated. Based on these experiments, several germination stage-related phenotypes (list below) were calculated and subjected to a GWAS. The mean value of the three biological replicates was calculated and used in the further analysis.
Germination Percentage (GP)
GP was recorded daily for 10 days and was calculated using the formula below.
GP = Number of germinated seeds at 10 days/Total number of seeds tested × 100%
Germination energy (GE)
GE was recorded daily for four days and was calculated using the formula below.
GE = Number of germinated seeds at four days/Total number of seeds tested × 100%
Germination index (GI)
GI was calculated using the formula below.
GI=Σ(Gt/t), where Gt is the number of seeds that germinated on day t (Alvarado et al. 1987, and Ruan et al. 2002).
Relative germination energy (R-GE)
R-GE was calculated by using the formula below.
R-GE = GE200/GEcontrol.
Relative germination index (R-GI)
R-GI was calculated by using the formula below.
R-GI = GI200/GIcontrol.
Relative root length (R-RL)
R-RLwas calculated by using the formula below R-RL = RL200/RLcontrol.
Relative shoot length (R-SL)
R-SLwas calculated by using the formula below R-SL = SL200/SLcontrol.
Relative R/S (R-R/S)
R-R/S was calculated by using the formula below R-R/S = R/S200/R/Scontrol.

4.4. Principal Components and GWAS Analysis

Principal components analysis (PCA) of the genotype with ~1.65 million SNPs/indels and four main salt-tolerance-related phenotypes: R-GE, R-GI, R-RL, and R-SL was conducted using GAPIT and Trait Analysis by Association, Evolution and Linkage (TASSEL) 5.0 [66]. Principal component analyses (PCA) in genotypic and phenotypic were also performed using GAPIT and TASSEL 5 [66].
GWAS was performed in the GAPIT package (Genome Association and Prediction Integrated Tool) in which an advanced kinship clustering algorithm was implemented [38]. Only SNPs with adjusted p-values < 0.05 were considered significantly associated. Gene loci containing the SNPs with significantly associated peaks in the Manhattan plot of the GWAS result were considered to be candidate genes related to salt tolerance.

4.5. Linkage Disequilibrium (LD) Block, Haplotype Analysis, and Expression Analysis

LD analysis. LD analysis was calculated using TASSEL 5 [66] based on the high-quality variations (with neither missing genotype calls over all accessions nor MAF < 0.05) in a 400 kb range determined by the most closely associated SNP/indel. An LD block was recognized when the top 95% confidence intervals of the D’ value exceeded 0.98 and the lower bounds exceeded 0.70 [67]. Loci with significant variations harbored by LD blocks were then defined as the candidate genes.
Haplotype analysis. With about 7.3× depth of genome coverage, we constructed the haplotyping of the identified candidate genes. Nucleotide polymorphisms on the target genes were captured according to the rice reference genome (IRGSP 1.0). The orthologous genes of the target candidate genes in several other plants were provided by Ensembl Plants (http://plants.ensembl.org). Alignments of orthologous gene sequences were conducted using Geneious (http://www.geneious.com) [68]. In addition, the TCS [69] haplotype network was conducted by PopART v 1.7 [70].
Gene expression analysis by qRT-PCR. Germinated seeds with shoot and root after 10 days in control (H2O) and salt (200 mM NaCl) conditions were collected and used for expression analysis. Total RNA was prepared using an RNA extraction kit (Qiagen, Hilden, Germany). cDNA was synthesized according to the manufacturer’s instructions using the PrimeScriptTM RT reagent Kit (TaKaRa, Shiga, Japan). Real-time PCR was carried out using the SYBR Green method with the primers of OsMADS31 (MADS-F: TGGCTTCACTGACTCTGCAA, MADS-R: TACATACCCGGCTGTGCATC). Relative expression levels were calculated using the 2−ΔΔCT method [71] with Ubiquitin 5 (UBQ 5) as the internal control [72] under three replicated tests.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/19/10/3145/s1.

Author Contributions

Conceptualization, Y.-J.P. Formal analysis, Q.H. Funding acquisition, Y.-J.P. Investigation, M.-Y.Y. and F.L. Methodology, M.-Y.Y. and B.C. Project administration, K.-W.K. and Y.-J.P. Software, W.T. and E.-B.H. Writing–original draft, J.Y. and W.Z.

Funding

National Research Foundation of Korea (NRF), Rural Development Administration, Republic of Korea (RDA).

Acknowledgments

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (NRF-2017R1A2B3011208). This work was carried out with the support of “Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ013405)” Rural Development Administration, Republic of Korea.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

GWASgenome-wide association study
LDlinkage disequilibrium
NGSnext-generation sequencing
SNPsingle-nucleotide polymorphism
INDELinsertion and deletion
MAFminor allele frequency
RLroot length
SLshoot length
RSroot/shoot ratio
GPGermination percentage
GEGermination energy
GIGermination index
CMLMcompressed mixed linear model
PCAPrincipal component analysis

References

  1. Sakadevan, K.; Nguyen, M.L. Extent, impact, and response to soil and water salinity in arid and semiarid regions. Adv. Agron. 2010, 109, 55. [Google Scholar]
  2. Mather, K.A.; Caicedo, A.L.; Polato, N.R.; Olsen, K.M.; McCouch, S.; Purugganan, M.D. The extent of linkage disequilibrium in rice (Oryza sativa L.). Genetics 2007, 177, 2223–2232. [Google Scholar] [CrossRef] [PubMed]
  3. 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 (Oryza sativa L.). PLoS ONE 2012, 7, e51202. [Google Scholar] [CrossRef] [PubMed]
  4. Almansouri, M.; Kinet, J.M.; Lutts, S. Effect of salt and osmotic stresses on germination in durum wheat (Triticum durum Desf.). Plant and Soil 2001, 231, 243–254. [Google Scholar] [CrossRef]
  5. Rehman, S.; Harris, P.J.; Bourne, W.F.; Wilkin, J. The relationship between ions, vigour and salinity tolerance of acacia seeds. Plant Soil 2000, 220, 229–233. [Google Scholar] [CrossRef]
  6. Xu, S.; Hu, B.; He, Z.; Ma, F.; Feng, J.; Shen, W.; Yang, J. Enhancement of salinity tolerance during rice seed germination by presoaking with hemoglobin. Int. J. Mol. Sci. 2011, 12, 2488–2501. [Google Scholar] [CrossRef] [PubMed]
  7. Khodarahmpour, Z.; Ifar, M.; Motamedi, M. Effects of NaCl salinity on maize (Zea mays L.) at germination and early seedling stage. Afr. J. Biotechnol. 2014, 11, 298–304. [Google Scholar] [CrossRef]
  8. Mutlu, F.; Bozcuk, S. Salinity-induced changes of free and bound polyamine levels in sunflower (Helianthus annuus L.) roots differing in salt tolerance. Pak. J. Bot. 2007, 39, 1097–1102. [Google Scholar]
  9. Ulfat, M.; Athar, H.U.R.; Ashraf, M.; Akram, N.A.; Jamil, A. Appraisal of physiological and biochemical selection criteria for evaluation of salt tolerance in canola (Brassica napus L.). Pak. J. Bot. 2007, 39, 1593–1608. [Google Scholar]
  10. Khan, M.A.; Weber, D.J. Ecophysiology of High Salinity Tolerant Plants; Springer Science & Business Media: New York, NY, USA, 2006. [Google Scholar]
  11. Koyro, H.W. Ultrastructural Effects of Salinity in Higher Plants. In Salinity, Environment-Plants-Molecules; Springer: New York, NY, USA, 2002; ISBN 978-1-4020-0492-6. [Google Scholar]
  12. Gomes-Filho, E.; Lima, C.R.F.M.; Costa, J.H.; da Silva, A.C.M.; Lima, M.d.G.S.; de Lacerda, C.F.; Prisco, J.T. Cowpea ribonuclease, properties and effect of NaCl-salinity on its activation during seed germination and seedling establishment. Plant Cell Rep. 2008, 27, 147–157. [Google Scholar] [CrossRef] [PubMed]
  13. Khan, M.A.; Rizvi, Y. Effect of salinity, temperature, and growth regulators on the germination and early seedling growth of Atriplex griffithii var. stocksii. Can. J. Bot. 1994, 72, 475–479. [Google Scholar] [CrossRef]
  14. Yupsanis, T.; Moustakas, M.; Eleftheriou, P.; Damianidou, K. Protein phosphorylation-dephosphorylation in alfalfa seeds germinating under salt stress. J. Plant Physiol. 1994, 143, 234–240. [Google Scholar] [CrossRef]
  15. Othman, Y.; Al-Karaki, G.; Al-Tawaha, A.; Al-Horani, A. Variation in germination and ion uptake in barley genotypes under salinity conditions. World J. Agric. Sci. 2006, 2, 11–15. [Google Scholar]
  16. Wahid, A.; Rasul, E.; Rao, A.U.R. Germination of Seeds and Propagules under Salt Stress; Handbook of Plant and Crop Stress: Boca Raton, FL, USA, 1999; Volume 2, pp. 153–167. [Google Scholar]
  17. Foolad, M.; Hyman, J.; Lin, G. Relationships between cold-and salt-tolerance during seed germination in tomato, Analysis of response and correlated response to selection. Plant Breed. 1999, 118, 49–52. [Google Scholar] [CrossRef]
  18. Wang, Z.F.; Wang, J.F.; Bao, Y.M.; Wu, Y.Y.; Zhang, H.S. Quantitative trait loci controlling rice seed germination under salt stress. Euphytica 2011, 178, 297–307. [Google Scholar] [CrossRef]
  19. Csanádi, G.; Vollmann, J.; Stift, G.; Lelley, T. Seed quality QTLs identified in a molecular map of early maturing soybean. Theor. Appl. Genet. 2001, 103, 912–919. [Google Scholar] [CrossRef]
  20. Bai, C.; Liang, Y.; Hawkesford, M.J. Identification of QTLs associated with seedling root traits and their correlation with plant height in wheat. J. Exp. Bot. 2013, 64, 1745–1753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. DeRose-Wilson, L.; Gaut, B.S. Mapping salinity tolerance during Arabidopsis thaliana germination and seedling growth. PLoS ONE 2011, 6, e22832. [Google Scholar] [CrossRef] [PubMed]
  22. Basnet, R.K.; Duwal, A.; Tiwari, D.N.; Xiao, D.; Monakhos, S.; Bucher, J.; Visser, R.G.F.; Groot, S.P.C.; Bonnema, G.; Maliepaard, C. Quantitative trait locus analysis of seed germination and seedling vigor in Brassica rapa reveals QTL hotspots and epistatic interactions. Front. Plant Sci. 2015, 6, 1032. [Google Scholar] [CrossRef] [PubMed]
  23. Hoang, T.M.L.; Tran, T.N.; Nguyen, T.K.T.; Williams, B.; Wurm, P.; Bellairs, S.; Mundree, S. Improvement of salinity stress tolerance in rice, challenges and opportunities. Agronomy 2016, 6, 54. [Google Scholar] [CrossRef]
  24. Wang, Z.; Wang, J.; Bao, Y.; Wang, F.; Zhang, H. Quantitative trait loci analysis for rice seed vigor during the germination stage. Zhejiang Univ. Sci. B (Biomed. & Biotechnol.) 2010, 11, 958–964. [Google Scholar] [Green Version]
  25. Cheng, J.; He, Y.; Yang, B.; Lai, Y.; Wang, Z.; Zhang, H. Association mapping of seed germination and seedling growth at three conditions in indica rice (Oryza sativa L.). Euphytica 2015, 206, 103–115. [Google Scholar] [CrossRef]
  26. Abe, A.; Takagi, H.; Fujibe, T.; Aya, K.; Kojima, M.; Sakakibara, H. OsGA20ox1.; a candidate gene for a major QTL controlling seedling vigor in rice. Theor. Appl. Genet. 2012, 125, 647–657. [Google Scholar] [CrossRef] [PubMed]
  27. Zheng, H.; Liu, B.; Zhao, H.; Wang, J.; Liu, H.; Sun, J.; Xing, J.; Zou, D. Identification of QTLs for salt tolerance at the germination and early seedling stage using linkage and association analysis in japonica rice. Chin. J. Rice Sci. 2014, 28, 358–366. [Google Scholar]
  28. Huang, X.; Wei, X.; Sang, T.; Zhao, Q.; Feng, Q.; Zhao, Y.; Li, C.; Zhu, C.; Lu, T.; Zhang, Z.; et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 2010, 42, 961–967. [Google Scholar] [CrossRef] [PubMed]
  29. Zhao, K.; Tung, C.W.; Eizenga, G.C.; Wright, M.H.; Ali, M.L.; Price, A.H.; Norton, G.J.; Islam, M.R.; Reynolds, A.; Mezey, J.; et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat. Commun. 2011, 2, 467. [Google Scholar] [CrossRef] [PubMed]
  30. Zhao, W.; Park, E.J.; Chung, J.W.; Park, Y.J.; Chung, I.M.; Ahn, J.K.; Kim, G.H. Association analysis of the amino acid contents in rice. J. Integr. Plant Biol. 2009, 51, 1126–1137. [Google Scholar] [CrossRef] [PubMed]
  31. Li, G.; Na, Y.W.; Kwon, S.W.; Park, Y.J. Association analysis of seed longevity in rice under conventional and high-temperature germination conditions. Plant Syst. Evol. 2014, 300, 389–402. [Google Scholar] [CrossRef]
  32. 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–39. [Google Scholar] [CrossRef] [PubMed]
  33. Ma, X.; Feng, F.; Wei, H.; Mei, H.; Xu, K.; Chen, S.; Li, T.; Liang, X.; Liu, H.; Luo, L. Genome-wide association study for plant height and grain yield in rice under contrasting moisture regimes. Front. Plant Sci. 2016, 7, 1801. [Google Scholar] [CrossRef] [PubMed]
  34. Kumar, V.; Singh, A.; Mithra, S.V.; Krishnamurthy, S.L.; Parida, S.K.; Jain, S.; Tiwari, K.K.; Kumar, P.; Rao, A.R.; Sharma, S.K.; et al. Genome wide association mapping of salinity tolerance in rice (Oryza sativa). DNA Res. 2015, 22, 133–145. [Google Scholar] [CrossRef] [PubMed]
  35. Yu, J.; Zhao, W.; He, Q.; Kim, T.S.; Park, Y.J. Genome-wide association study and gene set analysis for understanding candidate genes involved in salt tolerance at the rice seedling stage. Mol. Genet. Genomics 2017, 292, 1391–1403. [Google Scholar] [CrossRef] [PubMed]
  36. Shi, Y.; Gao, L.; Wu, Z.; Zhang, X.; Wang, M.; Zhang, C.; Zhang, F.; Zhou, Y.; Li, Z. Genome-wide association study of salt tolerance at the seed germination stage in rice. BMC Plant Biol. 2017, 17, 92. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Naveed, S.A.; Zhang, F.; Zhang, J.; Zheng, T.Q.; Meng, L.J.; Pang, Y.L.; Xu, J.L.; Li, Z.K. Identification of QTN and candidate genes for salinity tolerance at the germination and seedling stages in rice by genome-wide association analyses. Sci. Rep. 2018, 8, 6505. [Google Scholar] [CrossRef] [PubMed]
  38. Lipka, A.E.; Tian, F.; Wang, Q.; Peiffer, J.; Li, M.; Bradbury, P.J.; Gore, M.A.; Buckler, E.S.; Zhang, Z. GAPIT, genome association and prediction integrated tool. Bioinformatics 2012, 28, 2397–2399. [Google Scholar] [CrossRef] [PubMed]
  39. Kapoor, M.; Arora, R.; Lama, T.; Nijhawan, A.; Khurana, J.P.; Tyagi, A.K.; Kapoor, S. Genome-wide identification, organization and phylogenetic analysis of Dicer-like, argonaute and RNA-dependent RNA polymerase gene families and their expression analysis during reproductive development and stress in rice. BMC Genom. 2008, 9, 451. [Google Scholar] [CrossRef] [PubMed]
  40. Ricachenevsky, F.K.; Sperotto, R.A.; Menguer, P.K.; Sperb, E.R.; Lopes, K.L.; Fett, J.P. ZINC-INDUCED FACILITATOR-LIKE family in plants, lineage-specific expansion in monocotyledons and conserved genomic and expression features among rice (Oryza sativa) paralogs. BMC Plant Biol. 2011, 11, 20. [Google Scholar] [CrossRef] [PubMed]
  41. Okada, T.; Nakayama, H.; Shinmyo, A.; Yoshida, K. Expression of OsHAK genes encoding potassium ion transporters in rice. Plant Biotechnol. 2008, 25, 241–245. [Google Scholar] [CrossRef] [Green Version]
  42. Wang, J.; Chen, L.; Wang, Y.; Zhang, J.; Liang, Y.; Xu, D. A computational systems biology study for understanding salt tolerance mechanism in rice. PLoS ONE 2013, 8, e64929. [Google Scholar] [CrossRef] [PubMed]
  43. Arora, R.; Agarwal, P.; Ray, S.; Singh, A.K.; Singh, V.P.; Tyagi, A.K.; Kapoor, S. MADS-box gene family in rice, genome-wide identification, organization and expression profiling during reproductive development and stress. BMC Genom. 2007, 8, 242. [Google Scholar] [CrossRef] [PubMed]
  44. Zeng, L.; Shannon, M.; Grieve, C. Evaluation of salt tolerance in rice genotypes by multiple agronomic parameters. Euphytica 2002, 127, 235–245. [Google Scholar] [CrossRef]
  45. Ferdose, J.; Kawasaki, M.; Taniguchi, M.; Miyake, H. Differential sensitivity of rice cultivars to salinity and its relation to ion accumulation and root tip structure. Plant Prod. Sci. 2009, 12, 453–461. [Google Scholar] [CrossRef]
  46. Lee, K.S.; Choi, W.Y.; Ko, J.C.; Kim, T.S.; Gregorio, G.B. Salinity tolerance of japonica and indica rice (Oryza sativa L.) at the seedling stage. Planta. 2003, 216, 1043–1046. [Google Scholar] [CrossRef] [PubMed]
  47. Platten, J.D.; Egdane, J.A.; Ismail, A.M. 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. [Google Scholar] [CrossRef] [PubMed]
  48. Leon, T.B.; Steven, L.; Prasanta, K.S. Molecular dissection of seedling salinity tolerance in rice (Oryza sativa L.) using a high-density GBS-based SNP linkage map. Rice 2016, 9, 52. [Google Scholar] [CrossRef] [PubMed]
  49. Lin, H.X.; Yanagihara, S.; Zhuang, J.Y.; Senboku, T.; Zheng, K.L.; Yashima, S. Identification of QTLs for salt tolerance in rice via molecular markers. Chin. J. Rice Sci. 1998, 12, 72–78. [Google Scholar]
  50. Cai, H.W.; Morishima, H. QTL clusters reflect character associations in wild and cultivated rice. Theor. Appl. Genet. 2002, 104, 1217–1228. [Google Scholar] [PubMed]
  51. Courtois, B.; Audebert, A.; Dardou, A.; Roques, S.; Ghneim-Herrera, T.; Droc, G.; Frouin, J.; Rouan, L.; Goze, E.; Kilian, A.; et al. Genome-wide association mapping of root traits in a japonica rice panel. PLoS ONE 2013, 8, e78037. [Google Scholar] [CrossRef] [PubMed]
  52. Chen, C.; DeClerck, G.; Tian, F.; Spooner, W.; McCouch, S.; Buckler, E. PICARA.; an analytical pipeline providing probabilistic inference about a priori candidates genes underlying genome-wide association QTL in plants. PLoS ONE 2012, 7, e46596. [Google Scholar] [CrossRef] [PubMed]
  53. Sur, I.; Tuupanen, S.; Whitington, T.; Aaltonen, L.A.; Taipale, J. Lessons from functional analysis of genome-wide association studies. Cancer Res. 2013, 73, 4180–4184. [Google Scholar] [CrossRef] [PubMed]
  54. Xu, X.; Liu, X.; Ge, S.; Jensen, J.D.; Hu, F.; Li, X.; Dong, Y.; Gutenkunst, R.N.; Fang, L.; Huang, L.; et al. Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes. Nat. Biotechnol. 2012, 30, 105–111. [Google Scholar] [CrossRef] [PubMed]
  55. Reddy, I.N.B.L.; Kim, B.K.; Yoon, I.S.; Kim, K.H.; Kwon, T.R. Salt tolerance in rice, focus on mechanisms and approaches. Rice Sci. 2017, 24, 123–144. [Google Scholar] [CrossRef]
  56. Yang, Y.; Zhong, J.; Ouyang, Y.D.; Yao, J. The integrative expression and co-expression analysis of the AGO gene family in rice. Gene 2013, 528, 221–235. [Google Scholar] [CrossRef] [PubMed]
  57. Inahashi, H.; Shelley, I.J.; Yamauchi, T.; Nishiuchi, S.; Takahashi-Nosaka, M.; Matsunami, M.; Ogawa, A.; Noda, Y.; Inukai, Y. OsPIN2, which encodes a member of the auxin efflux carrier proteins, is involved in root elongation growth and lateral root formation patterns via the regulation of auxin distribution in rice. Physiol. Plantarum 2018. [Google Scholar] [CrossRef] [PubMed]
  58. Rebecca, M.D.; Patrick, A.R.; Patricia, M.M.; Jan, E.L. Germins, a diverse protein family important for crop improvement. Plant Sci. 2009, 177, 499–510. [Google Scholar]
  59. Kothari, K.S.; Dansana, P.K.; Giri, J.; Tyagi, A.K. Rice stress associated protein 1 (OsSAP1) interacts with aminotransferase (OsAMTR1) and pathogenesis-related 1a protein (OsSCP) and regulates abiotic stress responses. Front. Plant Sci. 2016, 7, 1057. [Google Scholar] [CrossRef] [PubMed]
  60. Konishi, S.; Izawa, T.; Lin, SY.; Ebana, K.; Fukuta, Y.; Sasaki, T.; Yano, M. An SNP caused loss of seed shattering during rice domestication. Science 2006, 312, 1392–1396. [Google Scholar] [CrossRef] [PubMed]
  61. Hu, B.; Wang, W.; Ou, S.; Tang, J.; Li, H.; Che, R.; Zhang, Z.; Chai, X.; Wang, H.; Wang, Y. Variation in NRT1.1B contributes to nitrate-use divergence between rice subspecies. Nat. Genet. 2015, 47, 834–838. [Google Scholar] [CrossRef] [PubMed]
  62. Kim, K.W.; Chung, H.K.; Cho, G.T.; Ma, K.H.; Chandrabalan, D.; Gwag, J.G.; Kim, T.S.; Cho, E.G.; Park, Y.J. PowerCore, a program applying the advanced M strategy with a heuristic search for establishing core sets. Bioinformatics 2007, 23, 2155–2162. [Google Scholar] [CrossRef] [PubMed]
  63. Zhao, W.; Cho, G.T.; Ma, K.H.; Chung, J.W.; Gwag, J.G.; Park, Y.J. Development of an allele-mining set in rice using a heuristic algorithm and SSR genotype data with least redundancy for the post-genomic era. Mol. Breed. 2010, 26, 639–651. [Google Scholar] [CrossRef]
  64. Kawahara, Y.; de la Bastide, M.; Hamilton, J.P.; Kanamori, H.; McCombie, W.R.; Ouyang, S.; Schwartz, D.C.; Tanaka, T.; Wu, J.; Zhou, S.; et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice 2013, 6, 4. [Google Scholar] [CrossRef] [PubMed]
  65. Kim, T.S.; He, Q.; Kim, K.W.; Yoon, M.Y.; Ra, W.H.; Li, F.P.; Tong, W.; Yu, J.; Oo, W.H.; Choi, B.; et al. Genome-wide resequencing of KRICE_CORE reveals their potential for future breeding.; as well as functional and evolutionary studies in the post-genomic era. BMC Genom. 2016, 17, 408. [Google Scholar] [CrossRef] [PubMed]
  66. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL, software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef] [PubMed]
  67. Gabriel, S.B.; Schaffner, S.F.; Nguyen, H.; Moore, J.M.; Roy, J.; Blumenstiel, B.; Higgins, J.; DeFelice, M.; Lochner, A.; Faggart, M. The structure of haplotype blocks in the human genome. Science 2002, 296, 2225–2229. [Google Scholar] [CrossRef] [PubMed]
  68. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C. Geneious basic, an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  69. Clement, M.; Posada, D.; Crandall, K.A. TCS, a computer program to estimate gene genealogies. Mol. Ecol. 2000, 9, 1657–1659. [Google Scholar] [CrossRef] [PubMed]
  70. Leigh, J.W.; Bryant, D. PopART, Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  71. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  72. Jain, M.; Nijhawan, A.; Tyagi, A.K.; Khurana, J.P. Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem. Biophys. Res. Commun. 2006, 345, 646–651. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 03145 g001 550
Figure 1. Determination of the optimum NaCl concentration and main phenotypes under salt stress and control conditions. (a) Individual value plot of germination percentage in the presence of 0, 200, and 300 mM NaCl (each dot represents an individual). (b–d) Box plots for phenotypic values in the presence of 0 and 200 mM NaCl (the asterisks are extreme outliers). GP: germination percentage, GE: germination energy, GI: Germination index, SL: shoot length, RL: root length.
Figure 1. Determination of the optimum NaCl concentration and main phenotypes under salt stress and control conditions. (a) Individual value plot of germination percentage in the presence of 0, 200, and 300 mM NaCl (each dot represents an individual). (b–d) Box plots for phenotypic values in the presence of 0 and 200 mM NaCl (the asterisks are extreme outliers). GP: germination percentage, GE: germination energy, GI: Germination index, SL: shoot length, RL: root length.
Ijms 19 03145 g001
Ijms 19 03145 g002 550
Figure 2. Principal components analysis (PCA) using genotype and phenotype data. (a). For genotype data, 295 accessions were divided into ind (indica) and jap (japonica) based on PC1 and PC2 along with an admixture group. (b). For phenotype data, no clear grouping was observed.
Figure 2. Principal components analysis (PCA) using genotype and phenotype data. (a). For genotype data, 295 accessions were divided into ind (indica) and jap (japonica) based on PC1 and PC2 along with an admixture group. (b). For phenotype data, no clear grouping was observed.
Ijms 19 03145 g002
Ijms 19 03145 g003 550
Figure 3. Genome-wide association mapping and LD block analysis for root length (RL) under salt stress (200 mM) conditions. (a) Manhattan plot from association mapping using the CMLM. (b) QQ plot of expected and observed P values. (c) The peak region on chromosome 4 along with the LD blocks. (d) The peak region on chromosome 12 along with the LD blocks. In (c,d) pair-wise LD between SNPs is indicated as D’ values: red indicates a value of 1 and yellow indicates 0. The LD region was 200 kb upstream and downstream of the top −log (p) value in the peak range.
Figure 3. Genome-wide association mapping and LD block analysis for root length (RL) under salt stress (200 mM) conditions. (a) Manhattan plot from association mapping using the CMLM. (b) QQ plot of expected and observed P values. (c) The peak region on chromosome 4 along with the LD blocks. (d) The peak region on chromosome 12 along with the LD blocks. In (c,d) pair-wise LD between SNPs is indicated as D’ values: red indicates a value of 1 and yellow indicates 0. The LD region was 200 kb upstream and downstream of the top −log (p) value in the peak range.
Ijms 19 03145 g003
Ijms 19 03145 g004 550
Figure 4. Haplotyping and sequence analysis of Os04g0614100, which was correlated with a phenotypic difference. (a) One functional SNP of the candidate gene in the CDS region. Type 1 is the reference and type 2 is the SNP. (b) The phenotypic difference based of the functional SNP. (c) Amino acid sequence alignment using several orthologues in various rice subgroups and species. Red box indicates the target amino acid change caused by the functional SNP. (d) RNA expression levels in rice accessions with type 1 and type 2. (e) Haplotype network analysis. Circle size is proportional to the number of samples within a given haplotype. Lines between haplotypes represent mutational steps between alleles. Colors denote rice designation.
Figure 4. Haplotyping and sequence analysis of Os04g0614100, which was correlated with a phenotypic difference. (a) One functional SNP of the candidate gene in the CDS region. Type 1 is the reference and type 2 is the SNP. (b) The phenotypic difference based of the functional SNP. (c) Amino acid sequence alignment using several orthologues in various rice subgroups and species. Red box indicates the target amino acid change caused by the functional SNP. (d) RNA expression levels in rice accessions with type 1 and type 2. (e) Haplotype network analysis. Circle size is proportional to the number of samples within a given haplotype. Lines between haplotypes represent mutational steps between alleles. Colors denote rice designation.
Ijms 19 03145 g004
Table
Table 1. Descriptive statistics for the traits in the control and salt-treated (200 mM NaCl) rice accessions.
Table 1. Descriptive statistics for the traits in the control and salt-treated (200 mM NaCl) rice accessions.
TraitSalinity Level (NaCl/mM)Mean ± SD aRangeMedianIQR b
GP00.97 ± 0.060.47–1.000.990.97–1.00
2000.87 ± 0.180.14–1.000.940.83–0.98
GE00.95 ± 0.080.37–1.000.980.94–0.99
2000.25 ± 0.290–0.970.130.03–0.39
GI047.66 ± 6.7118.33–57.8747.5543.40–53.44
20020.66 ± 8.531.51–43.7221.0515.48–26.03
SL06.07 ± 1.383.25–12.625.825.08–6.90
2001.45 ± 0.630.38–3.871.291.03–1.76
RL011.28 ± 2.894.28–18.8911.119.46–12.81
2000.69 ± 0.490.013–3.440.550.39–0.87
a Standard deviation. b Interquartile range. GP: Germination percentage. GE: Germination energy. GI: Germination index. SL: shoot length. RL: root length.
Table
Table 2. Pearson correlation coefficients among traits under control and salt stress (200 mM NaCl) conditions.
Table 2. Pearson correlation coefficients among traits under control and salt stress (200 mM NaCl) conditions.
TraitGPGEGISLRL
ControlGP
GE0.946 ***
GI0.624 ***0.710 ***
SL−0.001 ns0.037 ns0.168 **
RL0.073 ns0.087 ns0.110 ns0.254 ***
200 mMGP
GE0.394 ***
GI0.758 ***0.849 ***
SL0.439 ***0.738 ***0.725 ***
RL0.357 ***0.618 ***0.594 ***0.712 ***
GP: germination percentage. GE: germination energy. GI: germination index. SL: shoot length. RL: root length. *, **, ***, ns: significant at the 0.05, 0.01, and 0.001 probability level and not significant, respectively.
Table
Table 3. Candidate genes with highly associated signals in the coding region that overlapped with the GWAS and LD analysis.
Table 3. Candidate genes with highly associated signals in the coding region that overlapped with the GWAS and LD analysis.
Chr_Pos aTraitp-ValueFDR bGene IDDescription
chr02_19605493R/S1.72 × 10−70.00779Os02g0532500Germin family protein, Germin-like protein 2-4
Os02g0532900Glycoside hydrolase family 17 protein
Os02g0533300Carbonic anhydrase, CAH1-like domain, containing protein
Os02g0533800Similar to ATPase inhibitor
chr04_31168058RL2.42 × 10−70.00859Os04g0612900Vacuolar ATPase assembly integral membrane protein VMA21-like domain-containing protein
Os04g0613900Similar to Potassium transporter 18, OsHAK11
Os04g0614000Similar to Peroxisomal 2,4-dienoyl-CoA reductase
Os04g0614100MADS-box domain-containing protein, OsMADS31
Os04g0614600Similar to Viroid RNA-binding protein, aminotransferase
Os04g0614500Pyridoxal phosphate-dependent transferase, major region, subdomain 1 domain-containing protein
Os04g0615100Similar to Lecithine cholesterol acyltransferase-like protein
Os04g0615700Protein argonaute 2, OsAGO2
chr12_1628276RL, R/S, R-R/S2.02 × 10−70.00859Os12g0133100Major facilitator superfamily protein, OsZIFL12
Os12g0133300zinc-induced facilitator-like 13, OsZIFL13
Os12g01334004′-phosphopantetheinyl transferase domain-containing protein
Os12g0133700Stress-activated protein kinase pathway-regulating phosphatase 1
Os12g0133800Similar to Auxin efflux carrier protein, OsPIN1d
a The position was based on the annotation data on Os-Nipponbare-Reference-IRGSP-1.0 (RAP-DB, http://rapdb.dna.affrc.go.jp/). b FDR: False discovery rate. FDR Adjusted p values were calculated by GAPIT applying the Benjamini-Hochberg (1995) FDR-controlling procedure. RL: root length. R/S: root/shoot ratio. R-R/S: relative root/shoot ratio.

Share and Cite

MDPI and ACS Style

Yu, J.; Zhao, W.; Tong, W.; He, Q.; Yoon, M.-Y.; Li, F.-P.; Choi, B.; Heo, E.-B.; Kim, K.-W.; Park, Y.-J. A Genome-Wide Association Study Reveals Candidate Genes Related to Salt Tolerance in Rice (Oryza sativa) at the Germination Stage. Int. J. Mol. Sci. 2018, 19, 3145. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19103145

AMA Style

Yu J, Zhao W, Tong W, He Q, Yoon M-Y, Li F-P, Choi B, Heo E-B, Kim K-W, Park Y-J. A Genome-Wide Association Study Reveals Candidate Genes Related to Salt Tolerance in Rice (Oryza sativa) at the Germination Stage. International Journal of Molecular Sciences. 2018; 19(10):3145. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19103145

Chicago/Turabian Style

Yu, Jie, Weiguo Zhao, Wei Tong, Qiang He, Min-Young Yoon, Feng-Peng Li, Buung Choi, Eun-Beom Heo, Kyu-Won Kim, and Yong-Jin Park. 2018. "A Genome-Wide Association Study Reveals Candidate Genes Related to Salt Tolerance in Rice (Oryza sativa) at the Germination Stage" International Journal of Molecular Sciences 19, no. 10: 3145. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19103145

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop