Next Article in Journal
Identification and Expression Analysis of Gretchen Hagen 3 (GH3) in Kiwifruit (Actinidia chinensis) During Postharvest Process
Next Article in Special Issue
Genetic Diversity, Population Structure, and Marker-Trait Association for Drought Tolerance in US Rice Germplasm
Previous Article in Journal / Special Issue
Assessment of the Genetic Diversity of Rice Germplasms Characterized by Black-Purple and Red Pericarp Color Using Simple Sequence Repeat Markers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 8
Issue 11
10.3390/plants8110472
plants-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

Identification of Superior Alleles for Seedling Stage Salt Tolerance in the USDA Rice Mini-Core Collection

by
Jai S. Rohila
1,*,
Jeremy D. Edwards
1,
Gioi D. Tran
2,
Aaron K. Jackson
1 and
Anna M. McClung
1
1
Dale Bumpers National Rice Research Center, United States Department of Agriculture (USDA)—Agricultural Research Service (ARS), Stuttgart, AR 72160, USA
2
Genetics and Plant Breeding Division, Cuu Long Delta Rice Research Institute, Tan Thanh, Thoi Lai, Can Tho 905660, Vietnam
*
Author to whom correspondence should be addressed.
Plants 2019, 8(11), 472; https://doi.org/10.3390/plants8110472
Submission received: 1 October 2019 / Revised: 30 October 2019 / Accepted: 31 October 2019 / Published: 5 November 2019
(This article belongs to the Special Issue Genomics for Plant Breeding)
Browse Figures
Versions Notes

Abstract

:
Salt stress is a major constraint to rice acreage and production worldwide. The purpose of this study was to evaluate the natural genetic variation available in the United States Department of Agriculture (USDA) rice mini-core collection (URMC) for early vigor traits under salt stress and identify quantitative trait loci (QTLs) for seedling-stage salt tolerance via a genome-wide association study (GWAS). Using a hydroponic system, the seedlings of 162 accessions were subjected to electrical conductivity (EC) 6.0 dS m−1 salt stress at the three-to-four leaf stage. After completion of the study, 59.4% of the accessions were identified as sensitive, 23.9% were identified as moderately tolerant, and 16.7% were identified as highly tolerant. Pokkali was the most tolerant variety, while Nerica-6 was the most sensitive. Adapting standard International Rice Research Institute (IRRI) protocols, eight variables associated with salt tolerance were determined. The GWAS of the URMC, using over three million single-nucleotide polymorphisms (SNPs), identified nine genomic regions associated with salt tolerance that were mapped to five different chromosomes. Of these, none were in the known Saltol QTL region, suggesting different probable genes and mechanisms responsible for salt tolerance in the URMC. The study uncovered genetic loci that explained a large portion of the variation in salt tolerance at the seedling stage. Fourteen highly salt-tolerant accessions, six novel loci, and 16 candidate genes in their vicinity were identified that may be useful in breeding for salt stress tolerance. Identified QTLs can be targeted for fine mapping, candidate gene verification, and marker-assisted breeding in future studies.

1. Introduction

Rice is a staple food for billions of people, and it is the only food crop that provides more than 50% of daily calories by direct consumption to more than half of the world’s population [1]. According to recent projections, the worldwide demand for rice is expected to increase [2], but the long-term environmental sustainability of rice production is under question due to its large fresh-water irrigation requirements. Rice is generally grown in flooded paddies and surface water is often supplemented with water pumped from underground wells, which sometimes is of poor quality due to the presence of high amounts of dissolved salts, especially in depleted aquifers [3,4]. Moreover, a rice crop can be exposed to salt stress due to natural soil salinity or by salt-water intrusion in irrigation resources and fields in low-lying coastal areas. Salt-water intrusion in rice fields has occurred due to sea-level rise and inundations during tidal surges from tropical storms or tsunamis along the Gulf Coast in the USA, Japan, Vietnam, and elsewhere [5,6]. Rice is susceptible to the quantity and the quality of irrigation water, as it directly affects rice yield and grain quality [7]. Salt stress is one of the most important yield-limiting abiotic stresses in flooded lowland rice production areas and because there are limited opportunities for the mitigation of salt-water intrusion, efforts toward crop adaptation through breeding are needed. In general, rice is affected by saline conditions above 2 dS m−1 at the seedling stage and at reproductive stages [8,9,10]. Salinity conditions in the field during early plant growth stages affect stand establishment, and later in the growth cycle can influence many morphological (e.g., root length and diameter, photosynthetically active canopy area, plant height), physiological (e.g., net photosynthesis, stomatal conductance), and biochemical processes (e.g., accumulation of reactive oxygen species, osmolytes, ion homeostasis), which eventually contribute to grain yield reduction [9,11].
Multiple studies have identified several basic molecular mechanisms by which plants respond to and confer resistance to salinity stress [12,13,14,15]. These include: transporters, detoxification of Na+-induced reactive oxygen species, storage (or compartmentalization) of salts in the apoplast and vacuole, secretion of salt directly out of the plant system via salt glands, an increase of osmolytes, and many more [16,17,18]. Knowledge of the genes responsible for these various mechanisms is critical for providing effective molecular tools for rice-breeding programs aimed at improving various aspects of salt tolerance.
Molecular markers have been identified that are linked to easily measurable phenotypic traits that convey salt tolerance such as shoot and root dry weight, plant height, root length, and Na+/K+ accumulation/exclusion in the shoot [19]. Research on the genetics of salt tolerance, including identification of the Saltol quantitative trait loci (QTL) [20], has aided in the adaptation of rice plants to salt stress, and has been found to primarily control shoot Na+/K+ homeostasis. However, considering the numerous processes that salt stress impacts in plants, combining multiple QTLs is needed to achieve higher levels of tolerance throughout the plant growth cycle [19]. Prior to the rice genome being sequenced, salt–stress tolerance mutants and major genes were mapped with the use of low-throughput markers such as: RFLP, RAPD, AFLP and SSR [20,21]. Genome sequencing has enabled the high-resolution mapping of quantitative traits using large high-throughput panels of single-nucleotide polymorphisms (SNPs) [22,23,24,25]. These molecular markers coupled with marker-assisted selection (MAS)-based breeding programs offer a means of incorporating salt tolerance into recently developed high-yielding germplasm [26].
However, the continued development of salt tolerant lines through breeding will require the identification of germplasm possessing genetic variability across a wide spectrum of salt tolerance-related traits. Compared to biparental mapping populations, diversity panels provide a much larger pool of genetic resources for the discovery of beneficial genetic variants (or alleles) from landraces and germplasm that can be examined for traits that are important components of salt stress tolerance and their associated QTLs [19,27]. The dissection of these traits by phenotyping and genotyping will reveal new insights into the biological mechanisms underlying these phenotypes. Numerous diversity panels have been developed and extensively genotyped using various platforms [28,29,30,31,32,33,34]. The United States Department of Agriculture (USDA) rice mini-core collection (URMC) is a subset of 217 accessions containing similar phenotypic, genotypic, and geographical diversity as present in the USDA rice core collection of 1794 accessions representing rice-growing countries around the globe [35,36]. The URMC has been successfully used by rice researchers for identifying sources of germplasm and novel QTLs for several important traits such as: grain yield [35,37,38], grain yield components and harvest index [39], silica content in hulls [40], grain quality [23], grain protein content [41], grain starch quality [42], sheath blight resistance [43,44], and cold tolerance at seed germination and the seedling stage [45]. These evaluations have provided insights into the beneficial phenotypic and genotypic diversity available in the URMC for rice breeding programs.
The purpose of this study was to identify novel sources of tolerance to salt stress at the seedling stage across the diversity of Oryza sativa represented by the URMC, to compare the range of phenotypic diversity for salt tolerance in the URMC to that of current international salt tolerant varieties, and to identify salt tolerance-associated genomic regions to enable the introduction of superior salt tolerance alleles into new rice varieties.

2. Results

2.1. Phenotypic Variation for Seedling Salt Tolerance Traits in USDA Rice Mini-Core

We evaluated a panel of 162 diverse rice (Oryza sativa L.) accessions from URMC that represented four rice subpopulations: Tropical japonica (TRJ), Temperate japonica (TEJ), Indica (IND), and Aus (AUS) [35]. Six cultivars were included as salt-sensitive (IR29, Nerica 6, A69-1) and salt-tolerant (Pokkali, IR45427-2B-2-2B-1-1, IR65196-3B-5-2-2) checks [17,46,47]. A detailed view of all the traits measured in this study is provided in Table 1a,b. This study successfully differentiated salt-tolerant genotypes from salt-sensitive genotypes as evident by the Least Squares Means (LSMeans) (Table 1a) and the range values (Table 1b) of salt stress injury (SSI) scores. The day 16 (d16) LSMeans and range values in the sensitive checks were 7.53 and 5.2–9.00, respectively compared to 2.32 and 1.00–4.60 for the tolerant checks. The accessions in URMC showed a wide variation for SSI scores on d10 and d16, which ranged from a value of 1 to 9 with an LSMean of 6.39 on d16. Among the four subpopulations in the URMC, TRJ was identified as the most sensitive, with a mean SSI score of 5.63 on d10, and 7.55 on d16. The IND subpopulation (and two accessions in the TEJ) was the most tolerant. On average, the sensitive checks only increased in plant height by 8.25 cm following 10 d of salt exposure, whereas the tolerant checks grew 22.44 cm (Table 1a). The URMC and breeding lines from Vietnam were intermediate to the checks in plant height (PHT), and there were no significant differences among the four subpopulations at either d10 or d16. A similar trend was observed for the change in number of green leaves following 14 d of salt stress, with the URMC and breeding lines from Vietnam being intermediate to the checks. The subpopulations overlapped in number of green leaves, but the TRJ group was more senescent than the IND. A similar trend was observed for total, shoot, and root biomass. For the other traits that were measured, subpopulations with low SSI scores were associated with higher values for biomass, Δ green leaf number, and, although not significant, Δ PHT (Table 1a). As shown in Table 1b, a wide range in variation was observed for all these traits in URMC and among the lines from Vietnam; however, except for the SSI d10 score, there were no accessions as tolerant as the tolerant checks.
A two-way dendrogram was constructed based on the SSI score at d16, ΔPHT, Δ green leaf number, and total biomass plant−1 (Figure 1). The results demonstrate clear differences between the tolerant (e.g., Pokkali) and the sensitive check cultivars (e.g., IR 29). Although the majority of URMC accessions were clustered as susceptible, there were a number of accessions that were deemed highly or moderately tolerant.
For a better understanding of the relationships among the measured traits in these 162 accessions, we computed Pearson’s correlation coefficient between each trait combination (Table 2). All traits were highly correlated with each other and contribute to the growth and senescence factors that influence the SSI scores. The Δ green leaf number had a strong negative correlation with SSI scores at both d10 (r = −0.80, p < 0.0001) and d16 (r = −0.86, p < 0.0001) and a strong positive correlation with the other measured traits such as shoot biomass (r = 0.71, p < 0.0001) and root biomass (r = 0.72, p < 0.0001).

2.2. Identification of SNPs and QTLs Associated with Early Vigor and Tolerance for Salt Injury under Salt Stress

To identify genomic regions that were associated with salt tolerance in rice, we performed a genome-wide association study (GWAS) using 3.2 million SNP markers available for the URMC and mapped eight traits that were highly correlated with salt tolerance at the seedling stage. Manhattan plots of GWAS results and the identified segments were investigated, and nine genomic regions of interest were selected based on the significance of association across the traits, number of significant SNPs, and allele frequencies. The nine selected genomic regions were represented by nine unique SNPs, which were distributed on five different chromosomes of rice and were associated with one or more of the eight studied traits. Based on proximity to the centromere and lack of recombination, two selected regions on chromosome 3 were merged into a single region leaving a total of eight selected regions, seven of which were significantly associated with the SSI score and one that was associated with plant height. Figure 2 shows the approximate locations of the seven SSI d16 score-associated representative SNPs (Table 3a,b). Manhattan plots and respective quantile–quantile (Q-Q) distributions are presented as Figure S1. Table 3a shows the selected SNP chromosome positions (in base pairs, bp) with reference and alternate alleles, frequencies for the reference allele in different subpopulations (based on 3000 sequenced rice accessions [48,49]), and the subpopulations with significant SNP–phenotype associations. Additionally, information on all the significant regions was identified, and overlapping significant regions across all subpopulations and all traits are provided in Tables S1 and S2.
The effects of each of the selected alternate alleles (relative to the reference allele) were calculated using the associated phenotypic traits in the subpopulation(s) where the significant association was detected (Table 3b). Each of the selected SNPs, except for SNP 1_33398135, were found to be associated with SSI score, Δ green leaf number, and root and shoot biomass. Three SNPs were significantly associated with Δ PHT, SSI score, Δ green leaf number, and root and shoot biomass. SNP 1_33398135 was associated with only one measured trait, Δ PHT. The highest observed allele effect for Δ PHT was 3.998 for SNP 1_33398135. The total phenotypic variance (Rsquare) of the associated trait explained by the identified selected SNP ranged from 6% to 71% with a median value of 33%.
To understand the potential of these SNPs in providing salt tolerance, putative candidate genes within the eight selected chromosomal regions were identified. In each chromosomal region, the SNP with lowest p-value was considered as the peak SNP. For identifying putative genes, we took an approach of selecting the genes close to the peak SNP and considering only those that were annotated by the rice scientific community (Oryzabase, MSU, RAP-DB), and published reports were used as evidence of the possible roles of the identified SNPs. Our analysis resulted in 16 candidate genes (Table 4), which are physically shown on each chromosome in Figure 3. A number of genes shown in this table have been found to encode transcription factors, including functional proteins such as Apetala2 (AP2) and Oligosaccharyltransferase3 (OST3), respectively. Additionally, there were several genes identified in this study that have not been implicated as having a direct role in salt tolerance, but are shown to have broader roles in overall abiotic stress tolerance. These genes may serve as novel targets in improving our understanding of salt tolerance in rice at the seedling stage.

3. Discussion

Here, we present the results from a study with 162 rice accessions coming from the URMC diversity panel (N = 123) that included four subpopulations of Oryza sativa L., along with 24 entries from a Vietnam breeding program, seven entries from US breeding germplasm, and six varieties that have been documented to display a range in susceptibility for salt exposure. Previously, it had been reported that seedling vigor under salt stress is an effective method of selection in rice breeding programs [63]. Thus, we grew the accessions hydroponically and, at the 3–4 leaf stage, rice seedlings were exposed to a constant EC 6.0 dS m−1 salt stress for 16 days. Genomic and phenotypic data were analyzed using 118 of the accessions that had been resequenced [24] to identify salt tolerant genotypes, QTLs, and SNPs using GWAS.
We demonstrated that four traits (SSI score, ΔPHT, Δ green leaf number, and total biomass plant−1) were effective at differentiating salt-tolerant from salt-sensitive genotypes, and we observed significant variation in the URMC for the tolerance traits (Table 1). While increased plant height, shoot biomass, and root biomass are important for salt tolerance at the seedling stage, the ability to produce new leaves and maintain the functions of pre-existing leaves is also important for rice seedlings to remain physiologically and biochemically sustained under salt stress. Although none of the accessions were more tolerant than Pokkali, a known salt-tolerant genotype, 14 URMC accessions were statistically clustered as highly tolerant (Figure 1). In addition, about half of the advanced lines from Vietnam’s salt tolerance breeding program were also scored as highly tolerant (data not shown). Due to the concern of salt-water intrusion in many lowland rice-growing areas in Vietnam, they have effectively deployed the Saltol allele from Pokkali in a number of these varieties. However, there were several accessions in the URMC that performed better than the advanced salt-tolerant Vietnamese varieties but lacked the Saltol allele, indicating that there is the opportunity to pyramid diverse alleles to make additional gains in salt tolerance. Moreover, gains that come from multiple genes may act through synergistic mechanisms beyond what might be predicted from the estimated additive effects of each allele alone, and that is where trait-based modeling for salt tolerance in rice could support rice-breeding programs [64].
SSI scores, taken at 10 and 16 d post exposure to salt, were highly correlated (r = 0.90, p < 0.0001) and, as expected, the plants received a lower SSI score (less senescence) at the d10 rating (Table 1). This indicates that evaluations of germplasm for salt tolerance can be effectively determined at this earlier time point, which will increase the speed of selection. A strong correlation was determined between Δ green leaf number and root biomass (r = 0.7286, p < 0.0001). This suggests that maintaining a photosynthetically active leaf area under salt stress conditions is an adaptive tolerance mechanism [65]. Further, a stable or increased number of green leaves and biomass in tolerant germplasm under salt conditions may allow more salt to be sequestered in the vacuoles of roots and shoots, as excessive concentrations of Na+ in cytoplasm is detrimental to the physiological and biochemical processes of rice plants. The vacuolar sequestration of Na+ is known to support better seedling survival and plant growth under such stressed conditions [66,67].
Lee et al. [68] reported that tolerance to salinity was greater in Indica than Japonica subspecies. Our GWAS analysis was conducted with 118 URMC accessions that were represented by about twice as many Indica accessions than either the Tropical japonica or Aus subpopulations, affording the opportunity of identifying novel salt-tolerant alleles from Indica that would benefit US breeding programs, which are predominantly TRJ-based. GWAS analysis utilizing over three million SNPs identified nine unique SNPs that were highly associated with eight desirable salt tolerant traits.
In GWAS, the ability to detect a phenotypic association with a SNP is dependent on the allele frequency of the SNP and linkage disequilibrium with the functional variant. There is low statistical power to detect associations with rare alleles, and SNPs with low minor allele frequency are filtered out in GWAS pipelines. When SNP–phenotype associations were found in only some subpopulations, this was often a consequence of little or no SNP variation in the subpopulation(s) lacking the association (Figure S2). For some SNPs, there were moderate allele frequencies present in some subpopulations where associations were not detected, and other subpopulations with significant associations. This may suggest that the SNP was not in linkage disequilibrium with the functional variant in all the subpopulations or that the genetic background of the subpopulation may be influencing the genetic pathways involved in salt tolerance.
Considering the complex nature of salt stress tolerance, each identified SNP may have a small individual effect; thus, considering the aggregate effects of multiple SNPs using MAS can be a viable option for improving rice varieties [69]. This variation in phenotypic traits among diverse germplasm and between the subpopulations is an indication that genetic variation is available in specific subpopulations and could be utilized in a wide diversity of rice-breeding programs. However, the IND subpopulation had more desirable salt tolerant traits (Table 1) and desirable alleles at the identified SNPs than AUS or TRJ (TEJ was only represented by two accessions, and hence was not analyzed independently).
SNP 1_2594296 was identified with the alternate allele associated with a lower SSI score and increased leaf number and biomass in the IND group, although this allele is also found in moderate abundance in the TEJ and AUS subpopulations (Table 3a,b). The salt-tolerant alternate allele is predicted in 20.3% of IND suggesting this could be selected for improving salt tolerance when breeding within IND. It should be noted that in this study, the alternate allele does not appear to be associated with the phenotype in the AUS subpopulation, so for using this SNP in other subpopulations via introgression and MAS, it will be necessary to verify that the introgression is from an IND source that is likely to be linked with the functional variant for salt tolerance and not from an AUS source. SNP S1_33398135 has an alternate allele that is found almost exclusively in IND and is associated with an increased value for the Δ PHT trait under salt conditions. Likewise, SNP 2_2053382 should be useful in rice breeding because it has an allele that was found almost exclusively in IND and was associated with a reduced SSI score and increased Δ PHT, Δ green leaf number, and the biomass traits.
A study by Shi et al. [70] identified a QTL for vigor index (a measure of seed germination and shoot elongation during salt stress) on rice chromosome 2 that occurred 1.38 Mb proximal to the SNP S2_2053382 region identified in our study. The Shi et al. study utilized 478 entries [305 Indica, and 85 Japonica, 65 Aus, 16 Basmati (ARO), and seven intermediate (admixture) accessions] with 6,361,920 SNPs from the 3000 genome data set and performed a GWAS for seed germination under 60 mM of salt (NaCl) stress. This locus may be another useful target for MAS to improve salt tolerance using diversity from IND. The Os02g0452 gene encodes a G-protein γ subunit and is reported to be a signal transducer during salt stress, and is a likely candidate for salt stress tolerance in this region [52,71].
SNP S3_17374343 is located within a broad GWAS peak spanning the centromere of chromosome 3 (18.9–19.9 Mb). TRJ and TEJ groups have almost exclusively the reference allele. Variability at the SNP is limited primarily to the AUS and IND groups where the alternate allele is associated with an increased SSI score, and reduced Δ green leaf number and biomass. Thus, this is an alternate allele to be selected against in rice breeding for improving salt tolerance. However, MAS may be complicated in this area due to the recombination suppression that occurs near the centromeric region, which also contains several important yield-related genes [72]. SNP S3_36149293 has an alternate allele found in the AUS group that is associated with a reduced SSI score and increased Δ green leaf number and biomass. This salt tolerance-associated allele from AUS may be a useful target for MAS to improve salt tolerance by introgression to the rice-breeding lines of other subpopulations. Gibberellin 20 oxidase, GA20OX1, is a likely candidate gene for this locus. GA20OX1 is a paralog of the semi-dwarf Sd-1 “green revolution” gene that has been found to increase cytokinin activity in rice [73]. The alternate allele of SNP S4_31361839 is the most frequent allele in all subpopulations except in TEJ, and is associated with an increased SSI score, as well as reduced Δ green leaf number and biomass. The reference allele is associated with higher salt tolerance and is found in most TEJ accessions (81.1%) and a moderate number of IND accessions (26.5%). Either TEJ or IND could be used as a source of salt tolerance alleles in MAS for this locus to improve the germplasm in other subpopulations.
A GWAS study performed by Patishtan et al. [30] examined 306 rice accessions utilizing a 700 K high-density SNP chip after three durations of 50-mM salt treatments: short (6 h), medium (7 d), and long (30 d) on 15-day-old seedlings. Eight of their identified SNPs occurred within 1.5 Mb or less of the representative SNPs identified in this study. They reported two SNPs associated with d7 shoot potassium level on chromosome 1 and on chromosome 4, which are in proximity to the S1_2594296 and S4_31361839 SNPs of our study. Four SNPs on chromosome 10 associated with 6-h shoot sodium level were identified near SNP S10_11743519, which was identified in this study. A SNP each for the 6-h shoot sodium level and d30 shoot sodium level were identified on chromosome 10 near the SNP S10_18801757. GLP4-1, an auxin-binding gene, is a likely candidate gene for this locus [74]. The SNP 10_11743519 alternate allele was found primarily in the TEJ (41.9%) and TRJ (22.8%) groups, and was associated with a reduced SSI score, and increased Δ PHT, Δ green leaf number, and biomass.
We identified two SNPs, S10_11743519 and SNP 10_18801757, on chromosome 10. Genomic regions in the proximity of these SNPs were also identified in a GWA study by Kumar et al. [75] that examined salt stress using a 6K SNP chip and 220 rice accessions. The authors identified a SNP within 1.3 MB of S10_11743519 that was associated with spikelet fertility and productive tillers under season-long EC 10 dS m−1 salt (7 NaCl: 1 Na2SO4: 2 CaCl2) stress in field conditions. The salt tolerance-associated alternate allele would be a good target for MAS to improve salt tolerance within TEJ or through the introgression of the TEJ or TRJ allele into IND or AUS. This locus is an example of the diversity from Japonica contributing to salt tolerance, whereas other reported loci derive tolerance from IND or AUS sources. ERF51, a likely candidate gene located in this region (Table 4), is also reported to be associated with cold tolerance [76]. The SNP 10_18801757 alternate allele is found primarily in AUS and IND subpopulations, and is associated with lower SSI scores, and increased Δ PHT, Δ green leaf number, and biomass. The salt tolerance-associated allele from IND or AUS would be a good target for MAS to improve salt tolerance within these groups or to introgress into other subpopulations such as TRJ or TEJ. As discussed above, six SNPs (1_2594296, 2_2053382, S3_36149293, S4_31361839, 10_11743519, and 10_18801757) in this study were identified as being useful for MAS.
In total, the study identified six potential novel loci that are promising targets for MAS with tolerant alleles coming from the Indica subspecies in five cases and the Japonica subspecies in one case. Pyramiding multiple salt tolerance alleles into a single genetic background is possible through molecular breeding, and may result in highly salt tolerant germplasm that could be a useful parent for rice-breeding programs. This study presents previously untapped genetic variation and the identification of novel rice germplasm that can be used to improve salt tolerance in rice, which in turn can help increase rice production to meet the demands of the world’s increasing population.

4. Materials and Methods

4.1. Plant Material, Growth Conditions, Salt Stress Treatment

A group of 162 rice accessions (30 AUS, 93 IND, two TEJ, and 37 TRJ) representing 55 rice-growing countries was selected from the URMC along with 24 modern varieties from a Vietnamese breeding program and seven accessions with known salt tolerance or salt susceptibility along with three salt-tolerant checks (Pokkali, IR45427-2B-2-2B-1-1, and IR65196-3B-5-2-2) and three salt-sensitive checks (IR29, Nerica 6, and A69-1) were used to evaluate salt tolerance at the seedling stage (Table S3). Of these, DNA sequence data was available for 118 URMC accessions [30 AUS, 58 IND, and 30 Japonicas (two TEJ, and 28 TRJ)] [23,24] for genomic analysis. Hydroponic experiments were conducted under greenhouse conditions at the USDA-ARS Dale Bumpers National Rice Research Center, Stuttgart, Arkansas, USA. Seeds were planted according to the protocol of Gregorio et al. [68,77] with minor modifications. Briefly, each experimental unit was represented by two healthy seeds placed into each of 5 wells (1.5-cm diameter) that had been created in a 30 cm (L), 27 cm (W), 1.5 cm (D) Styrofoam floater with a fiberglass mesh attached to the bottom to support the rice seeds in the nutrient solution. Each floater was put in an individual 11.4-L black plastic tub that was 31.8 cm (L), 29 cm (W), 12.2 cm (D) in dimension. Twenty entries, including two repetitive check entries, Pokkali (resistant check) and IR 29 (susceptible check), were planted in each floater in a tub. The experiment was conducted as a randomized complete block design with three replications. Once seedlings were approximately 2 cm tall, the seedlings were thinned to one plant per well. When the seedlings were at the two-leaf stage, airsoft 6-mm pellets (Soft Air, Dallas, TX, USA) were added to each well to support the plant. Aluminum foil was wrapped around each tub and floater to prevent sunlight from entering the nutrient solution to help prevent algae formation. When >80% of seedlings were at the three-to-four leaf stage (generally about two weeks after germination), salt treatment was started (designated as day 0 or d_0) by adding NaCl (salt) to the nutrient solution. Electrical conductivity (EC) was kept constant at 6.0 ± 0.2. The pH (5 ± 0.2 and EC (6.0 ± 0.2) were adjusted every other day using 1 M of NaOH and 1 M of HCl. Nutrient solution was changed once a week to ensure the availability of sodium ions (Na+) for toxicity induction in rice seedlings. A side-by-side control treatment (without salt) was run on the check entries for quality control measures of the salt treatments.

4.2. Phenotypic Data Collection, Processing, and Analysis

Phenotypic symptoms of sodium toxicity in the form of salt stress injury (SSI) scores were evaluated on the 10th day (d10) and 16th day (d16) after subjecting the seedlings to salt treatment. SSI scores were recorded on a scale of 1 (tolerant) to 9 (susceptible), as suggested by Gregorio et al. [77] and shown in Figure S3. In this investigation, we also measured the following seven early vigor traits that were related to plant growth and development [78]: plant height on the day of initiating salt stress (PHT_d0), plant height on day 10 of salt stress (PHT_d10), plant height on day 16 of salt stress (PHT_d16), green leaf number on the day of initiating salt stress (GLN_d0), green leaf number on day 14 of salt stress (GLN_d14), shoot biomass on day 16 of salt stress (SB_d16), and root biomass on day 16 of salt stress (RB_d16). On completion of the 16th day of salt treatment, the rice seedlings were harvested by cutting them into two portions at the crown and dried individually at 60 °C for 96 h to determine the dry weights of: shoot biomass (SB) and root biomass (RB) separately.
The following four values for traits of interest were calculated as below:
  • Growth in plant height under 10 days of salt stress (Δ PHT-d10): PHT_d10 – PHT_d0
  • Growth in plant height under 16 days of salt stress (Δ PHT-d16): PHT_d16 – PHT_d0
  • Green leaf number (GLN) produced under 14 days salt stress (Δ green leaf#): GLN_d14 – GLN_d0
  • Total biomass after 16 days of salt stress (TB): SB_d16 + RB_d16
The best linear unbiased prediction (BLUP) estimates of the means for all traits were used for genomic analysis as explained by Eizenga et al. [79]. Briefly, to estimate the BLUP values for each measured and calculated trait, the variation among rice accessions was considered as a random effect. In the model, a normal distribution was used for all the studied traits. The values for each parameter were the means of three independent replications where each replication had five plants. An ANOVA was performed on growth variables using the standard least squares procedure of JMP13.2 (SAS Institute, Cary, NC, USA). Tukey’s HSD (Honestly Significant Difference) analysis was used for the comparisons of the growth parameters measured among the genotypes. Pearson’s correlation coefficients were determined between pairs of traits using JMP14.3.0 (SAS Institute, Cary, NC, USA).

4.3. Genome-Wide Association Mapping, Identification of Single-Nucleotide Polymorphisms (SNPs) and Candidate Genes

The USDA rice mini-core collection (URMC) was sequenced to an average read depth of 1.5x by Wang et al. [24], and the raw reads were later recalled as described in Huggins et al. [23]. A mixed linear model (MLM) in Tassel version 5 [80,81] was used to perform a genome-wide association study (GWAS) on 118 accessions in the URMC. Of these 118 accessions, six panels were used to capture the genetic diversity present that is specific to each subpopulation/subspecies. The first panel, “All”, was comprised of 118 accessions with 3,029,968 SNPs and contained four rice subpopulations present in this study; IND, AUS, TEJ, and TRJ. No aromatic accessions from the URMC were used in this study. The “IND” panel contained 58 indica subpopulation accessions and 1,883,603 SNPs. The “AUS” panel had 30 aus accessions and 1,753,773 SNPs. The “TRJ” panel had 28 tropical japonica accessions and 1,067,722 SNPs. Only two temperate japonica accessions represented the Japonica subspecies.
For each panel, SNPs were filtered to a minimum allele frequency of 0.05 to remove rare alleles, and SNPs with 20% or greater missing calls were removed. Tassel version 5 [80] was used to generate a centered identity by state kinship matrix and principal component analysis (PCA) for each filtered panel. The kinship matrix and first three principal components were used in the MLM and the options of no compression of the kinship matrix and the variance components estimated for each marker-trait combination were used.
The R script qqman [82] was used to generate Manhattan and Q-Q plots. For exploratory purposes, SNPs with −log10 of p values >5 were considered significant, and increased stringency was applied during data interpretation as needed. A Perl script was used to define the associated chromosome regions from individual SNPs or groups of physically linked SNPs. Chromosome regions encompassed 50 Kb in both directions around each individual significant SNP, and were extended to include nearby significant SNPs occurring within 200 Kb. The Perl script designated a ‘Peak SNP’ for each region, which corresponded to the SNP with the most significant p-value found within the region.
A Perl script was used to help identify candidate gene(s) within 250 Kb of the significant regions, as described in Huggins et al. [23]. This Perl script used gene annotation from MSU7 [83] and the Rice Annotation Project (RAP1; http://rapdb.dna.affrc.go.jp/, accessed 14 September 2018) [84] plus two new annotation files that were created to aid in candidate gene identification. The first file was based on the results of the candidate genes identified in Cohen and Leach [22], which examined general biotic and abiotic stress response genes in rice. The second file merged all MSU7 annotations [83] (accessed 30 May 2019) with data in Oryzabase (OrzbaseGeneListEn_20190424010057; https://shigen.nig.ac.jp/rice/oryzabase/download/gene, accessed 24 April 2019) that shared an RAP ID. As described in Huggins et al. [23], an additional Perl script was used to examine the regions that overlapped across the different traits and subpopulation/subspecies panels.
SNPs were grouped into chromosomal regions by processing GWAS results with custom Perl scripts using the methods reported by Huggins et al. [23]. Based on p-values and allele effects, chromosomal regions and corresponding peak SNPs that are highly predictive of SSI scores and related phenotypes were selected for detailed analysis. Allele effects, significance, and R-squared values for each of the selected SNPs were calculated separately for each subpopulation by ANOVA in JMP 14.0.0. To determine allele distribution among a larger panel of rice lines than the URMC alone, we further examined the representative significant SNPs in public databases. The allele frequency information was computed by using RiceVarMap, which contains genotype information on 4726 rice accessions [49]. If allele frequency information was not present in RiceVarMap, then it was computed from the SNP-seek database, which contains genotype information on 3000 rice accessions [48,85].

Supplementary Materials

The following are available online at https://www.mdpi.com/2223-7747/8/11/472/s1. Figure S1: Manhattan and quantile–quantile (Q-Q) plots for the entire GWA analysis conducted in this study, Figure S2: Violin plots showing effect of each identified significant SNP/allele (x axis) on the salt stress injury (SSI) score on d16 (y axis) of salt stress in each rice subpopulations (x axis), Figure S3: Representative pictures of criteria used for salt stress injury (SSI) scoring for salt tolerance, Table S1: Significant regions identified across all subpopulations/panels for all measured traits, Table S2: Overlapping significant regions identified across all subpopulations/panels for all measured traits, Table S3: A complete list of the rice accessions evaluated for salt tolerance at three-four leaf seedling stage.

Author Contributions

J.S.R., J.D.E., A.M.M., G.D.T. conceived and designed the experiments. G.D.T. conducted the experiments under J.S.R.’s supervision. J.S.R., J.D.E., A.K.J., A.M.M., G.D.T. analyzed and interpreted the data, J.R. wrote the first draft of the manuscript. All authors revised and approved the final submitted version.

Funding

This research received no external funding.

Acknowledgments

This research utilized facilities and assistance provided by the USDA Agricultural Research Service Dale Bumpers National Rice Research Center. Rice (Oryza sativa) accessions were obtained from the USDA-ARS, Dale Bumpers National Rice Research Center, Stuttgart, Arkansas, Genetic Stocks Oryza Collection (www.ars.usda.gov/GSOR). The authors acknowledge help received from Ms. Tiffany Sookaserm in performing the hydroponic experiment. Financial support from the Cuu Long Delta Rice Research Institute (CLRRI), Vietnam in the form of a six-month scholarship to TDG is also acknowledged. The authors also thank Drs. Trevis Huggins and Santosh Sharma for their critical peer review. The USDA is an equal opportunity provider and employer.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McNally, K.L.; Childs, K.L.; Bohnert, R.; Davidson, R.M.; Zhao, K.; Ulat, V.J.; Zeller, G.; Clark, R.M.; Hoen, D.R.; Bureau, T.E.; et al. Genomewide SNP variation reveals relationships among landraces and modern varieties of rice. Proc. Natl. Acad. Sci. USA 2009, 106, 12273–12278. [Google Scholar] [CrossRef] [Green Version]
  2. Nigatu, G.; Hansen, J.; Childs, N.; Seeley, R. Sub-Saharan Africa Is Projected to Be the Leader in Global Rice Imports; Amber Waves: 2017, United States Department of Agriculture, Economic Research Service. Available online: https://www.ers.usda.gov/amber-waves/2017/october/sub-saharan-africa-is-projected-to-be-the-leader-in-global-rice-imports/ (accessed on 30 September 2019).
  3. Kang, M.; Jackson, R.B. Salinity of deep groundwater in California: Water quantity, quality, and protection. Proc. Natl. Acad. Sci. USA 2016, 113, 7768–7773. [Google Scholar] [CrossRef] [Green Version]
  4. Verma, A.K.; Gupta, S.K.; Isaac, R.K. Use of saline water for irrigation in monsoon climate and deep water table regions: Simulation modeling with SWAP. Agric. Water Manag. 2012, 115, 186–193. [Google Scholar] [CrossRef]
  5. Wassmann, R.; Hien, N.X.; Hoanh, C.T.; Tuong, T.P. Sea level rise affecting the Vietnamese Mekong delta: Water elevation in the flood season and implications for rice production. Clim. Chang. 2004, 66, 89–107. [Google Scholar] [CrossRef]
  6. Takagi, H.; Tamiru, M.; Abe, A.; Yoshida, K.; Uemura, A.; Yaegashi, H.; Obara, T.; Oikawa, K.; Utsushi, H.; Kanzaki, E.; et al. MutMap accelerates breeding of a salt-tolerant rice cultivar. Nat. Biotechnol. 2015, 33, 445. [Google Scholar] [CrossRef]
  7. Reba, M.L.; Daniels, M.; Chen, Y.; Sharpley, A.; Bouldin, J.; Teague, T.G.; Daniel, P.; Henry, C.G. A statewide network for monitoring agricultural water quality and water quantity in Arkansas. J. Soil Water Conserv. 2013, 68, 45A–49A. [Google Scholar] [CrossRef] [Green Version]
  8. Aguilar, M.; Fernández-Ramírez, J.L.; Aguilar-Blanes, M.; Ortiz-Romero, C. Rice sensitivity to saline irrigation in Southern Spain. Agric. Water Manag. 2017, 188, 21–28. [Google Scholar] [CrossRef]
  9. Grattan, S.R.; Zeng, L.; Shannon, M.C.; Roberts, S.R. Rice is more sensitive to salinity than previously thought. Calif. Agric. 2002, 56, 189–198. [Google Scholar] [CrossRef] [Green Version]
  10. 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]
  11. Zeng, L.; Shannon, M.C. Salinity effects on seedling growth and yield components of rice. Crop Sci. 2000, 40, 996–1003. [Google Scholar] [CrossRef]
  12. Morton, M.J.L.; Awlia, M.; Al-Tamimi, N.; Saade, S.; Pailles, Y.; Negrão, S.; Tester, M. Salt stress under the scalpel—Dissecting the genetics of salt tolerance. Plant J. 2019, 97, 148–163. [Google Scholar] [CrossRef] [PubMed]
  13. Munns, R. Genes and salt tolerance: Bringing them together. New Phytol. 2005, 167, 645–663. [Google Scholar] [CrossRef] [PubMed]
  14. 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]
  15. Li, Y.-F.; Zheng, Y.; Vemireddy, L.R.; Panda, S.K.; Jose, S.; Ranjan, A.; Panda, P.; Govindan, G.; Cui, J.; Wei, K.; et al. Comparative transcriptome and translatome analysis in contrasting rice genotypes reveals differential mRNA translation in salt-tolerant Pokkali under salt stress. BMC Genom. 2018, 19, 935. [Google Scholar] [CrossRef] [PubMed]
  16. Deinlein, U.; Stephan, A.B.; Horie, T.; Luo, W.; Xu, G.; Schroeder, J.I. Plant salt-tolerance mechanisms. Trends Plant Sci. 2014, 19, 371–379. [Google Scholar] [CrossRef] [Green Version]
  17. Negrão, S.; Courtois, B.; Ahmadi, N.; Abreu, I.; Saibo, N.; Oliveira, M.M. Recent updates on salinity stress in rice: From physiological to molecular responses. Crit. Rev. Plant Sci. 2011, 30, 329–377. [Google Scholar] [CrossRef]
  18. Yuan, F.; Leng, B.; Wang, B. Progress in studying salt secretion from the salt glands in recretohalophytes: How do plants secrete salt? Front. Plant Sci. 2016, 7, 977. [Google Scholar] [CrossRef]
  19. Thomson, M.J.; De Ocampo, M.; Egdane, J.; Rahman, M.A.; Sajise, A.G.; Adorada, D.L.; Tumimbang-Raiz, E.; Blumwald, E.; Seraj, Z.I.; Singh, R.K.; et al. Characterizing the Saltol quantitative trait locus for salinity tolerance in rice. Rice 2010, 3, 148–160. [Google Scholar] [CrossRef]
  20. Bonilla, P.; Dvorak, J.; Mackell, 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]
  21. Zhang, G.-Y.; Guo, Y.; Chen, S.-L.; Chen, S.-Y. RFLP tagging of a salt tolerance gene in rice. Plant Sci. 1995, 110, 227–234. [Google Scholar] [CrossRef]
  22. Cohen, S.P.; Leach, J.E. Abiotic and biotic stresses induce a core transcriptome response in rice. Sci. Rep. 2019, 9, 6273. [Google Scholar] [CrossRef] [PubMed]
  23. Huggins, T.D.; Chen, M.H.; Fjellstrom, R.G.; Jackson, A.K.; McClung, A.M.; Edwards, J.D. Association analysis of three diverse rice (Oryza sativa L.) germplasm collections for loci regulating grain quality traits. Plant Genome 2019, 12, 170085. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, H.; Xu, X.; Vieira, F.G.; Xiao, Y.; Li, Z.; Wang, J.; Nielsen, R.; Chu, C. The power of inbreeding: NGS-based GWAS of rice reveals convergent evolution during rice domestication. Mol. Plant 2016, 9, 975–985. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, W.; Mauleon, R.; Hu, Z.; Chebotarov, D.; Tai, S.; Wu, Z.; Li, M.; Zheng, T.; Fuentes, R.R.; Zhang, F.; et al. Genomic variation in 3010 diverse accessions of Asian cultivated rice. Nature 2018, 557, 43–49. [Google Scholar] [CrossRef]
  26. Rana, M.M.; Takamatsu, T.; Baslam, M.; Kaneko, K.; Itoh, K.; Harada, N.; Sugiyama, T.; Ohnishi, T.; Kinoshita, T.; Takagi, H.; et al. Salt Tolerance Improvement in Rice through Efficient SNP Marker-Assisted Selection Coupled with Speed-Breeding. Int. J. Mol. Sci. 2019, 20, 2585. [Google Scholar] [CrossRef]
  27. Singh, A.K.; Ansari, M.W.; Pareek, A.; Singla-Pareek, S.L. Raising salinity tolerant rice: Recent progress and future perspectives. Physiol. Mol. Biol. Plants 2008, 14, 137–154. [Google Scholar] [CrossRef]
  28. McCouch, S.R.; Wright, M.H.; Tung, C.W.; Maron, L.G.; McNally, K.L.; Fitzgerald, M.; Singh, N.; DeClerck, G.; Agosto-Perez, F.; Korniliev, P.; et al. Open access resources for genome-wide association mapping in rice. Nat. Commun. 2016, 7, 10532. [Google Scholar] [CrossRef]
  29. Al-Tamimi, N.; Brien, C.; Oakey, H.; Berger, B.; Saade, S.; Ho, Y.S.; Schmöckel, S.M.; Tester, M.; Negrão, S. Salinity tolerance loci revealed in rice using high-throughput non-invasive phenotyping. Nat. Commun. 2016, 7, 13342. [Google Scholar] [CrossRef] [Green Version]
  30. Patishtan, J.; Hartley, T.N.; Fonseca de Carvalho, R.; Maathuis, F.J.M. Genome-wide association studies to identify rice salt-tolerance markers. Plant Cell Environ. 2018, 41, 970–982. [Google Scholar] [CrossRef]
  31. 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]
  32. Ahmadi, N.; Negrão, 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, 881–895. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, H.; Zhang, D.; Wang, M.; Sun, J.; Qi, Y.; Li, J.; Wei, X.; Han, L.; Qiu, Z.; Tang, S.; et al. A core collection and mini core collection of Oryza sativa L. in China. Theor. Appl. Genet. 2011, 122, 49–61. [Google Scholar] [CrossRef] [PubMed]
  34. Tiwari, K.K.; Singh, A.; Pattnaik, S.; Sandhu, M.; Kaur, S.; Jain, S.; Tiwari, S.; Mehrotra, S.; Anumalla, M.; Samal, R.; et al. Identification of a diverse mini-core panel of Indian rice germplasm based on genotyping using microsatellite markers. Plant Breed. 2015, 134, 164–171. [Google Scholar] [CrossRef]
  35. Agrama, H.A.; Yan, W.; Lee, F.; Fjellstrom, R.; Chen, M.-H.; Jia, M.; McClung, A. Genetic assessment of a mini-core subset developed from the USDA rice genebank. Crop Sci. 2009, 49, 1336–1346. [Google Scholar] [CrossRef]
  36. Yan, W.; Rutger, J.N.; Bryant, R.J.; Bockelman, H.E.; Fjellstrom, R.G.; Chen, M.-H.; Tai, T.H.; McClung, A.M. Development and evaluation of a core subset of the USDA rice germplasm collection. Crop Sci. 2007, 47, 869–878. [Google Scholar] [CrossRef]
  37. Li, X.; Yan, W.; Agrama, H.; Hu, B.; Jia, L.; Jia, M.; Jackson, A.; Moldenhauer, K.; McClung, A.; Wu, D. Genotypic and phenotypic characterization of genetic differentiation and diversity in the USDA rice mini-core collection. Genetica 2010, 138, 1221–1230. [Google Scholar] [CrossRef]
  38. Li, X.; Yan, W.; Agrama, H.; Jia, L.; Shen, X.; Jackson, A.; Moldenhauer, K.; Yeater, K.; McClung, A.; Wu, D. Mapping QTLs for improving grain yield using the USDA rice mini-core collection. Planta 2011, 234, 347–361. [Google Scholar] [CrossRef]
  39. Li, X.; Yan, W.; Agrama, H.; Jia, L.; Jackson, A.; Moldenhauer, K.; Yeater, K.; McClung, A.; Wu, D. Unraveling the complex trait of harvest index with association mapping in rice (Oryza sativa L.). PLoS ONE 2012, 7, e29350. [Google Scholar] [CrossRef]
  40. Bryant, R.; Proctor, A.; Hawkridge, M.; Jackson, A.; Yeater, K.; Counce, P.; Yan, W.; McClung, A.; Fjellstrom, R. Genetic variation and association mapping of silica concentration in rice hulls using a germplasm collection. Genetica 2011, 139, 1383–1398. [Google Scholar] [CrossRef]
  41. Bryant, R.J.; Jackson, A.K.; Yeater, K.M.; Yan, W.G.; McClung, A.M.; Fjellstrom, R.G. Genetic variation and association mapping of protein concentration in brown rice using a diverse rice germplasm collection. Cereal Chem. 2013, 90, 445–452. [Google Scholar] [CrossRef]
  42. Li, K.; Bao, J.; Corke, H.; Sun, M. Association Analysis of Markers Derived from Starch Biosynthesis Related Genes with Starch Physicochemical Properties in the USDA Rice Mini-Core Collection. Front. Plant Sci. 2017, 8, 424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Jia, L.; Yan, W.; Agrama, H.A.; Yeater, K.; Li, X.; Hu, B.; Moldenhauer, K.; McClung, A.; Wu, D. Searching for germplasm resistant to sheath blight from the USDA rice core collection. Crop Sci. 2011, 51, 1507–1517. [Google Scholar] [CrossRef]
  44. Jia, L.; Yan, W.; Zhu, C.; Agrama, H.A.; Jackson, A.; Yeater, K.; Li, X.; Huang, B.; Hu, B.; McClung, A.; et al. Allelic analysis of sheath blight resistance with association mapping in rice. PLoS ONE 2012, 7, e32703. [Google Scholar] [CrossRef] [PubMed]
  45. Schläppi, M.R.; Jackson, A.K.; Eizenga, G.C.; Wang, A.; Chu, C.; Shi, Y.; Shimoyama, N.; Boykin, D.L. Assessment of five chilling tolerance traits and GWAS mapping in rice using the USDA mini-core collection. Front. Plant Sci. 2017, 8, 957. [Google Scholar] [CrossRef] [PubMed]
  46. Chowdhury, A.D.; Haritha, G.; Sunitha, T.; Krishnamurthy, S.L.; Divya, B.; Padmavathi, G.; Ram, T.; Sarla, N. Haplotyping of rice genotypes using simple sequence repeat markers associated with salt tolerance. Rice Sci. 2016, 23, 317–325. [Google Scholar] [CrossRef]
  47. Yamamoto, A.; Sawada, H.; Shim, I.S.; Usui, K.; Fujihara, S. Effect of salt stress on physiological response and leaf polyamine content in NERICA rice seedlings. Plant Soil Environ. 2011, 57, 571–576. [Google Scholar] [CrossRef] [Green Version]
  48. Alexandrov, N.; Tai, S.; Wang, W.; Mansueto, L.; Palis, K.; Fuentes, R.R.; Ulat, V.J.; Chebotarov, D.; Zhang, G.; Li, Z.; et al. SNP-Seek database of SNPs derived from 3000 rice genomes. Nucleic Acids Res. 2015, 43, D1023–D1027. [Google Scholar] [CrossRef]
  49. Zhao, H.; Yao, W.; Ouyang, Y.; Yang, W.; Wang, G.; Lian, X.; Xing, Y.; Chen, L.; Xie, W. RiceVarMap: A comprehensive database of rice genomic variations. Nucleic Acids Res. 2014, 43, D1018–D1022. [Google Scholar] [CrossRef]
  50. Duan, Y.-B.; Li, J.; Qin, R.-Y.; Xu, R.-F.; Li, H.; Yang, Y.-C.; Ma, H.; Li, L.; Wei, P.-C.; Yang, J.-B. Identification of a regulatory element responsible for salt induction of rice OsRAV2 through ex situ and in situ promoter analysis. Plant Mol. Biol. 2016, 90, 49–62. [Google Scholar] [CrossRef]
  51. Kong, W.; Zhong, H.; Deng, X.; Gautam, M.; Gong, Z.; Zhang, Y.; Zhao, G.; Liu, C.; Li, Y. Evolutionary Analysis of GH3 Genes in Six Oryza Species/Subspecies and Their Expression under Salinity Stress in Oryza sativa ssp. japonica. Plants 2019, 8, 30. [Google Scholar] [CrossRef]
  52. Yadav, D.K.; Islam, S.M.S.; Tuteja, N. Rice heterotrimeric G-protein gamma subunits (RGG1 and RGG2) are differentially regulated under abiotic stress. Plant Signal. Behav. 2012, 7, 733–740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Jia, Q.; Zheng, C.; Sun, S.; Amjad, H.; Liang, K.; Lin, W. The role of plant cation/proton antiporter gene family in salt tolerance. Biol. Plant. 2018, 62, 617–629. [Google Scholar] [CrossRef]
  54. Xu, F.; Fang, J.; Ou, S.; Gao, S.; Zhang, F.; Du, L.; Xiao, Y.; Wang, H.; Sun, X.; Chu, J.; et al. Variations in CYP78A13 coding region influence grain size and yield in rice. Plant Cell Environ. 2015, 38, 800–811. [Google Scholar] [CrossRef] [PubMed]
  55. Han, R.; Jiang, W.; Tang, S.; Wan, J.; Long, Q.; Pan, X.; Shi, Q.; Wu, Z. Xanthine Dehydrogenase Involves in the Response of Photosystem and Reactive Oxygen Metabolism to Drought Stress in Rice. Russ. J. Plant Physiol. 2018, 65, 404–411. [Google Scholar] [CrossRef]
  56. Gu, J.-F.; Qiu, M.; Yang, J.-C. Enhanced tolerance to drought in transgenic rice plants overexpressing C4 photosynthesis enzymes. Crop J. 2013, 1, 105–114. [Google Scholar] [CrossRef] [Green Version]
  57. Agrawal, G.K.; Jwa, N.-S.; Shibato, J.; Han, O.; Iwahashi, H.; Rakwal, R. Diverse environmental cues transiently regulate OsOPR1 of the “octadecanoid pathway” revealing its importance in rice defense/stress and development. Biochem. Biophys. Res. Commun. 2003, 310, 1073–1082. [Google Scholar] [CrossRef]
  58. Farid, A.; Malinovsky, F.G.; Veit, C.; Schoberer, J.; Zipfel, C.; Strasser, R. Specialized Roles of the Conserved Subunit OST3/6 of the Oligosaccharyltransferase Complex in Innate Immunity and Tolerance to Abiotic Stresses. Plant Physiol. 2013, 162, 24–38. [Google Scholar] [CrossRef] [Green Version]
  59. Chen, G.; Zou, Y.; Hu, J.; Ding, Y. Genome-wide analysis of the rice PPR gene family and their expression profiles under different stress treatments. BMC Genom. 2018, 19, 720. [Google Scholar] [CrossRef]
  60. Li, L.; Xu, X.; Chen, C.; Shen, Z. Genome-Wide Characterization and Expression Analysis of the Germin-Like Protein Family in Rice and Arabidopsis. Int. J. Mol. Sci. 2016, 17, 1622. [Google Scholar] [CrossRef]
  61. Rashid, M.; Guangyuan, H.; Guangxiao, Y.; Hussain, J.; Xu, Y. AP2/ERF Transcription Factor in Rice: Genome-Wide Canvas and Syntenic Relationships between Monocots and Eudicots. Evol. Bioinform. 2012, 8, EBO-S9369. [Google Scholar] [CrossRef]
  62. Kim, Y.-O.; Kang, H.; Ahn, S.-J. Overexpression of phytochelatin synthase AtPCS2 enhances salt tolerance in Arabidopsis thaliana. J. Plant Physiol. 2019, 240, 153011. [Google Scholar] [CrossRef] [PubMed]
  63. Sahi, C.; Singh, A.; Kumar, K.; Blumwald, E.; Grover, A. Salt stress response in rice: Genetics, molecular biology, and comparative genomics. Funct. Integr. Genomics 2006, 6, 263–284. [Google Scholar] [CrossRef] [PubMed]
  64. Paleari, L.; Movedi, E.; Confalonieri, R. Trait-based model development to support breeding programs. A case study for salt tolerance and rice. Sci. Rep. 2017, 7, 4352. [Google Scholar] [CrossRef] [PubMed]
  65. Acosta-Motos, J.R.; Ortuño, M.F.; Bernal-Vicente, A.; Diaz-Vivancos, P.; Sanchez-Blanco, M.J.; Hernandez, J.A. Plant responses to salt stress: Adaptive mechanisms. Agronomy 2017, 7, 18. [Google Scholar] [CrossRef]
  66. Böhm, J.; Messerer, M.; Müller, H.M.; Scholz-Starke, J.; Gradogna, A.; Scherzer, S.; Maierhofer, T.; Bazihizina, N.; Zhang, H.; Stigloher, C.; et al. Understanding the molecular basis of salt sequestration in epidermal bladder cells of Chenopodium quinoa. Curr. Biol. 2018, 28, 3075–3085. [Google Scholar] [CrossRef]
  67. Wu, H.; Shabala, L.; Zhou, M.; Su, N.; Wu, Q.; Ul-Haq, T.; Zhu, J.; Mancuso, S.; Azzarello, E.; Shabala, S. Root vacuolar Na+ sequestration but not exclusion from uptake correlates with barley salt tolerance. Plant J. 2019, 100, 55–67. [Google Scholar] [CrossRef]
  68. 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]
  69. De Leon, T.B.; Linscombe, S.; Subudhi, P.K. 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]
  70. 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]
  71. Swain, D.M.; Sahoo, R.K.; Srivastava, V.K.; Tripathy, B.C.; Tuteja, R.; Tuteja, N. Function of heterotrimeric G-protein γ subunit RGG1 in providing salinity stress tolerance in rice by elevating detoxification of ROS. Planta 2017, 245, 367–383. [Google Scholar] [CrossRef]
  72. Blary, A.; Jenczewski, E. Manipulation of crossover frequency and distribution for plant breeding. Theor. Appl. Genet. 2019, 132, 575–592. [Google Scholar] [CrossRef] [PubMed]
  73. Wu, Y.; Wang, Y.; Mi, X.-F.; Shan, J.-X.; Li, X.-M.; Xu, J.-L.; Lin, H.-X. The QTL GNP1 encodes GA20ox1, which increases grain number and yield by increasing cytokinin activity in rice panicle meristems. PLoS Genet. 2016, 12, e1006386. [Google Scholar] [CrossRef] [PubMed]
  74. Davidson, R.M.; Manosalva, P.M.; Snelling, J.; Bruce, M.; Leung, H.; Leach, J.E. Rice Germin-Like Proteins: Allelic diversity and relationships to early stress responses. Rice 2010, 3, 43–55. [Google Scholar] [CrossRef]
  75. Kumar, V.; Singh, A.; Mithra, S.V.A.; 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]
  76. Wang, H.; Lee, A.-R.; Park, S.-Y.; Jin, S.-H.; Lee, J.; Ham, T.-H.; Park, Y.; Zhao, W.-G.; Kwon, S.-W. Genome-wide association study reveals candidate genes related to low temperature tolerance in rice (Oryza sativa) during germination. 3 Biotech 2018, 8, 235. [Google Scholar] [CrossRef]
  77. Gregorio, G.B.; Senadhira, D.; Mendoza, R.D. Screening Rice for Salinity Tolerance; IRRI Discussion Paper Series No. 22; International Rice Research Institute (IRRI): Manila, Philippines, 1997; p. 30. [Google Scholar]
  78. Negrão, S.; Schmöckel, S.M.; Tester, M. Evaluating physiological responses of plants to salinity stress. Ann. Bot. 2017, 119, 1–11. [Google Scholar] [CrossRef]
  79. Eizenga, G.C.; Jia, M.H.; Jackson, A.K.; Boykin, D.L.; Ali, M.L.; Shakiba, E.; Tran, N.T.; McCouch, S.R.; Edwards, J.D. Validation of yield component traits identified by genome-wide association mapping in a tropical japonica × tropical japonica rice biparental mapping population. Plant Genome 2019, 12, 180021. [Google Scholar] [CrossRef]
  80. 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]
  81. Zhang, Z.; Ersoz, E.; Lai, C.-Q.; Todhunter, R.J.; Tiwari, H.K.; Gore, M.A.; Bradbury, P.J.; Yu, J.; Arnett, D.K.; Ordovas, J.M.; et al. Mixed linear model approach adapted for genome-wide association studies. Nat. Genet. 2010, 42, 355. [Google Scholar] [CrossRef]
  82. Turner, S.D. Qqman: An R package for visualizing GWAS results using Q-Q and manhattan plots. Biorxiv 2014, 005165. [Google Scholar] [CrossRef]
  83. 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]
  84. Sakai, H.; Lee, S.S.; Tanaka, T.; Numa, H.; Kim, J.; Kawahara, Y.; Wakimoto, H.; Yang, C.-C.; Iwamoto, M.; Abe, T.; et al. Rice Annotation Project Database (RAP-DB): An Integrative and Interactive Database for Rice Genomics. Plant Cell Physiol. 2013, 54, e6. [Google Scholar] [CrossRef] [PubMed]
  85. Mansueto, L.; Fuentes, R.R.; Borja, F.N.; Detras, J.; Abriol-Santos, J.M.; Chebotarov, D.; Sanciangco, M.; Palis, K.; Copetti, D.; Poliakov, A.; et al. Rice SNP-seek database update: New SNPs, indels, and queries. Nucleic Acids Res. 2017, 45, D1075–D1081. [Google Scholar] [CrossRef] [PubMed]
Plants 08 00472 g001 550
Figure 1. Two-way phenotypic dendrogram showing three clusters for variation in salt tolerance among 162 rice accessions including 123 URMC accessions. Four early vigor traits associated with salt stress tolerance at the seedling stage were used for constructing the dendrogram and are denoted by numbers: 1: SSI score d16; 2: Δ PHT d16; 3: Δ green leaf#; and 4: Total biomass plant−1 on the x-axis of the heatmap. The blue color illustrates low values, while the red color illustrates high values for the measured trait.
Figure 1. Two-way phenotypic dendrogram showing three clusters for variation in salt tolerance among 162 rice accessions including 123 URMC accessions. Four early vigor traits associated with salt stress tolerance at the seedling stage were used for constructing the dendrogram and are denoted by numbers: 1: SSI score d16; 2: Δ PHT d16; 3: Δ green leaf#; and 4: Total biomass plant−1 on the x-axis of the heatmap. The blue color illustrates low values, while the red color illustrates high values for the measured trait.
Plants 08 00472 g001
Plants 08 00472 g002 550
Figure 2. Genome-wide association analysis for salt stress injury score (SSI) on day 16 of salt stress treatment using single-nucleotide polymorphisms (SNPs) from the USDA rice mini-core. Manhattan plots are shown for three rice subpopulations: (a) Tropical japonica (TRJ), (b) Aus (AUS), and (c) Indica (IND). The SNP positions are shown across the 12 rice chromosomes (X axis) with a −log10 (p) value for each SNP (Y axis). The horizontal black line at −log10(p) = 5 is the significance threshold used in the analysis. Black arrows indicate approximate regions for the representative significant SNPs.
Figure 2. Genome-wide association analysis for salt stress injury score (SSI) on day 16 of salt stress treatment using single-nucleotide polymorphisms (SNPs) from the USDA rice mini-core. Manhattan plots are shown for three rice subpopulations: (a) Tropical japonica (TRJ), (b) Aus (AUS), and (c) Indica (IND). The SNP positions are shown across the 12 rice chromosomes (X axis) with a −log10 (p) value for each SNP (Y axis). The horizontal black line at −log10(p) = 5 is the significance threshold used in the analysis. Black arrows indicate approximate regions for the representative significant SNPs.
Plants 08 00472 g002
Plants 08 00472 g003 550
Figure 3. Physical location of candidate genes on five rice chromosomes based on identified SNPs associated with various salt tolerant traits at the seedling stage. Saltol, a major quantitative trait loci (QTL) [20], is shown for reference. SNPs are shown in rectangular boxes on the left side of the chromosomes. The candidate genes are shown in italics on the right side of the chromosomes.
Figure 3. Physical location of candidate genes on five rice chromosomes based on identified SNPs associated with various salt tolerant traits at the seedling stage. Saltol, a major quantitative trait loci (QTL) [20], is shown for reference. SNPs are shown in rectangular boxes on the left side of the chromosomes. The candidate genes are shown in italics on the right side of the chromosomes.
Plants 08 00472 g003
Table 1. Comparison of United States Department of Agriculture (USDA) rice mini-core and its four subpopulations, breeding lines from Vietnam, and salt tolerant and sensitive controls for seedling salt stress: (a) Least Squares Mean (LSMean) values and Tukey’s HSD analysis; (b) Range [minimum (Min) – maximum (Max)] values. SSI = salt stress injury, PHT = plant height, URMC = USDA rice mini-core collection, AUS = Aus, IND = Indica, TEJ = Temperate japonica, TRJ = Tropical japonica.
Table
(a)
(a)
ClassificationGenotypeSSI Score-d10SSI Score-d16Δ PHT-d10 (cm)Δ PHT-d16 (cm)Δ Green Leaf#-d14Total Biomass Plant−1 (g)Shoot Biomass Plant−1 (g)Root Biomass Plant−1 (g)
ISensitive checks5.45 a7.53 a8.25 c8.29 d−2.97 d0.16 d0.13 d0.03 d
Tolerant checks1.38 d2.32 d22.44 a29.43 a0.64 a1.30 a1.08 a0.22 a
URMC *4.61 b6.39 b12.98 b14.07 b−2.06 c0.28 c0.23 c0.05 c
Vietnamese lines *3.52 c4.94 c10.45 c12.06 c−1.24 b0.35 b0.29 b0.07 b
II *AUS4.46 B6.50 B12.96 A13.94 A−2.08 A,B0.27 B0.22 B0.05 B
IND3.95 C5.48 C12.61 A14.00 A−1.76 A0.34 A0.27 A0.07 A
TEJ3.20 C4.93 B,C14.32 A15.45 A−1.33 A,B0.44 A0.36 A0.08 A
TRJ5.63 A7.55 A12.04 A12.89 A−2.24 B0.20 C0.16 C0.04 C
ALL δ4.43 ± 0.076.15 ± 0.0912.57 ± 0.2013.75 ± 0.23−1.93 ± 0.060.29 ± 0.010.24 ± 0.010.06 ± 0.001
* Checks (tolerant and sensitive) were excluded from calculations for the LSMeans and Tukey’s HSD analysis. Column values not followed by the same letter within a classification are significantly different (α = 0.050). Note: Lowercase letters are used for comparisons in I, and uppercase letters are used for comparisons in II. δ ALL in this table denotes all 162 entries.
Table
(b)
(b)
GenotypeMin/MaxSSI Score-d10SSI Score-d16Δ PHT-d10 (cm)Δ PHT-d16 (cm)Δ Green Leaf#-d14Total Biomass Plant−1 (g)Shoot Biomass Plant−1 (g)Root Biomass Plant−1 (g)
Sensitive checkMin3.805.202.602.60−4.000.110.010.01
Max7.409.0011.3011.30−2.000.270.210.06
Tolerant checkMin1.001.006.008.80−1.000.420.350.07
Max3.804.6030.2047.702.001.991.660.33
URMCMin1.001.803.604.00−4.000.030.030.01
Max9.009.0027.7030.401.000.990.810.18
Vietnamese linesMin1.402.000.004.40−5.000.140.110.02
Max6.608.2018.5020.701.000.810.670.15
AUS *Min1.302.005.005.00−4.000.090.070.02
Max7.409.0020.5025.901.000.810.660.14
IND *Min1.001.800.004.00−5.000.120.090.02
Max7.409.0027.7030.401.000.990.810.18
TEJ *Min1.602.0011.5011.50−3.000.210.170.04
Max4.607.2017.1019.101.000.710.580.13
TRJ *Min3.003.204.704.70−4.000.030.030.01
Max9.009.0025.5030.400.000.480.390.09
ALL *Min1.001.800.004.00−5.000.030.030.01
Max9.009.0027.7030.401.000.990.810.18
* Checks (tolerant and sensitive) were excluded from calculations for these values.
Table
Table 2. Correlation coefficients among seedling salt tolerance traits determined under hydroponic conditions (N = 162).
Table 2. Correlation coefficients among seedling salt tolerance traits determined under hydroponic conditions (N = 162).
SSI Score-d10SSI Score-d16Δ PHT-d10Δ PHT-d16Δ Green Leaf#-d14Total BiomassShoot BiomassRoot Biomass
SSI score-d101
SSI score-d160.9041 *1
Δ PHT-d10−0.4270 *−0.4309 *1
Δ PHT-d16−0.5366 *−0.5708 *0.9367 *1
Δ green leaf#-d14−0.7985 *−0.8577 *0.4030 *0.5393 *1
Total biomass−0.7218 *−0.7495 *0.5838 *0.7177 *0.7170 *1
Shoot biomass−0.7177 *−0.7417 *0.5885 *0.7191 *0.7112 *0.9991 *1
Root biomass−0.7267 *−0.7696 *0.5447 *0.6918 *0.7286 *0.9809 *0.9722 *1
* = p < 0.0001.
Table 3. GWAS analysis of USDA rice mini-core accessions: (a) Representative SNP locations with reference (Ref.) or alternate (Alt.) allele associated with salt-tolerant traits in specific subpopulations of the USDA rice mini-core collection, and subpopulation SNP allele frequencies across 3000 rice genomes; (b) Significant SNP effects in each subpopulation as estimated by ANOVA.
Table
(a)
(a)
SNP ID Chr.Position (bp)Ref. AlleleAlt. AlleleSub-Pop with Alt. Allele Associated with TraitsTraits Associated with the Alt. AlleleFrequencies for the Ref. Allele
TRJTEJAROAUSIND
S1_259429612,594,296GCINDSSI score, Δ green leaf#, biomass95.7%50.2%85.9%64.1%79.7%
S1_33398135133,398,135CTINDΔ PHT100.0%99.9%100.0%98.9%43.2%
S2_205338222,053,382ACINDSSI score, Δ PHT, Δ green leaf#, biomass99.9%100.0%100.0%99.2%65.7%
S3_17374343317,374,343CTAUS, INDSSI score, Δ green leaf#, biomass97.2%99.7%19.8%91.4%10.3%
S3_21582523321,582,523CTAUS, INDSSI score, Δ green leaf#, biomass90.1%98.0%5.2%41.3%44.2%
S3_36149293336,149,293ATAUSSSI score, Δ green leaf#, biomass99.8%100.0%100.0%57.2%98.4%
S4_31361839431,361,839GAINDSSI score, Δ green leaf#, biomass2.6%81.1%1.0%0.7%26.5%
S10_117435191011,743,519CTTRJSSI score, Δ PHT, Δ green leaf#, biomass77.2%58.1%100.0%100.0%99.6%
S10_188017571018,801,757ACAUS, INDSSI score, Δ PHT, Δ green leaf#, biomass96.6%50.2%95.8%63.9%86.1%
* All the frequencies for the reference alleles were computed using the SNP-seek database of the 3000 rice genomes [48] and RiceVarMap [49]. TRJ = Tropical japonica; TEJ = Temperate japonica; ARO = aromatic; AUS = Aus; IND = Indica. Chr. denotes chromosome number.
Table
(b)
(b)
SNPSub-PopTraitsProb>|t|RsquareAlternate Allele Effect
1_2594296INDd10 SSI score0.00070.32−1.247
d16 SSI score<0.00010.47−2.152
Δ green leaf#-d14<0.00010.431.043
Shoot biomass d16<0.00010.450.126
Root biomass d16<0.00010.650.035
1_33398135INDΔ PHT-d10<0.00010.483.543
Δ PHT-d16<0.00010.443.998
2_2053382INDd10 SSI score0.00200.22−0.904
d16 SSI score0.00020.30−1.659
Δ PHT-d100.09350.071.536
Δ PHT-d160.04260.102.215
Δ green leaf#-d140.00020.300.822
Shoot biomass d16<0.00010.390.118
Root biomass d16<0.00010.530.036
3_17374343AUSd10 SSI score0.01100.281.200
d16 SSI score0.09740.131.024
Shoot biomass d160.01960.24−0.079
Root biomass d160.02320.23−0.017
INDd10 SSI score<0.00010.311.708
d16 SSI score<0.00010.322.422
Δ green leaf#-d140.00010.29−1.117
Shoot biomass d160.00140.21−0.123
Root biomass d16<0.00010.31−0.038
3_21582523AUSd10 SSI score0.00080.511.336
d16 SSI score0.00520.401.509
Δ green leaf#-d140.02000.29−0.649
Shoot biomass d160.01160.34−0.079
Root biomass d160.00420.41−0.018
INDd10 SSI score0.00100.330.987
d16 SSI score0.00060.361.460
Δ green leaf#-d140.00060.36−0.740
Shoot biomass d160.02310.18−0.069
Root biomass d160.01350.21−0.018
3_36149293AUSd10 SSI score0.01240.25−0.800
d16 SSI score0.00080.41−1.426
Δ green leaf#-d140.00110.390.686
Shoot biomass d160.03310.190.051
Root biomass d160.00890.270.013
4_31361839INDd10 SSI score0.01960.401.113
d16 SSI score0.01510.431.654
Δ green leaf#-d140.01710.42−0.816
Shoot biomass d160.00030.71−0.164
Root biomass d160.00040.69−0.040
10_11743519TRJd10 SSI score<0.00010.58−1.248
d16 SSI score<0.00010.66−1.603
Δ PHT-d100.04160.212.016
Δ PHT-d160.00860.333.237
Δ green leaf#-d140.00040.510.710
Shoot biomass d160.00070.480.077
Root biomass d160.00330.390.015
10_18801757AUSd10 SSI score0.00010.46−1.185
d16 SSI score<0.00010.53−1.647
Δ green leaf#-d14<0.00010.460.819
Shoot biomass d160.01030.240.048
Root biomass d160.00280.300.013
INDd10 SSI score0.01220.13−0.817
d16 SSI score0.00070.22−1.579
Δ PHT-d160.09180.061.948
Δ green leaf#-d140.00130.200.723
Shoot biomass d160.00060.230.100
Root biomass d160.00050.230.026
Table
Table 4. Candidate genes in the chromosomal region of representative SNPs associated with salt stress tolerance at the seedling stage in rice.
Table 4. Candidate genes in the chromosomal region of representative SNPs associated with salt stress tolerance at the seedling stage in rice.
SNPMSU7 Locus IDPositionGene Name Putative FunctionReference
S1_2594296LOC_Os01g04800Chr01:2202264 - 2203860AP2/EREBP129Transcriptional reprogramming during salt stress [50]
S1_33398135LOC_Os01g57610Chr01:33308448 - 33311391OsGH3.1Auxin-responsive role in salinity tolerance[51]
S2_2053382LOC_Os02g04520Chr02:2007232 - 2009961G-protein γ subunitSignal transducer during salt stress[52]
LOC_Os02g04630Chr02:2070492 - 2075369Vacuolar cation/proton exchangerProtects primary cell mechanisms mediated by ion homeostasis[53]
S3_17374343 (and S3_21582523)LOC_Os03g30420Chr03:17340415 - 17340601GL3.2Cytochrome P450, stress tolerance[22,54]
LOC_Os03g31550Chr03:17985563 - 17998498Aldehyde oxidase putativeAbiotic stress response[22,55]
LOC_Os03g31750Chr03:18153143 - 18160127Pyruvate orthophosphate dikinase (PPDK)Improved plant physiology under abiotic stress [56]
LOC_Os03g32314Chr03:18485606 - 18488371AOC1, allene oxide cyclase Response to salt and other abiotic stresses[57]
LOC_Os03g33590Chr03:19197601 - 19204090Interferon-related developmental regulatorSalt sensitivity in Arabidopsis[22]
LOC_Os03g36730Chr03:20363507 - 20370047OST3/OST6 family proteinHypersensitivity to osmotic and salt stress[58]
LOC_Os03g37260Chr03:20650253 - 20653294Pentatricopeptide repeatPost-transcriptional gene regulation under abiotic stresses[59]
S3_36149293LOC_Os03g63970Chr03:36152044 - 36152517GA20OX1, gibberellin 20 oxidase 1Growth and development, downregulated under abiotic stress[22]
S4_31361839LOC_Os04g52720Chr04:31392876 - 31395185GLP4-1, germin-like protein 4-1Salt and osmotic stress tolerance [60]
S10_11743519LOC_Os10g22600Chr10:11731592 - 11732917ERF51, ethylene response factor 51Osmotic stress tolerance[61]
S10_18801757LOC_Os10g35150Chr10:18754181 - 18758096Expressed proteinDownregulated during abiotic stress[22]
LOC_Os10g35460Chr10:18976009 - 18979110Phytochelatin synthase Regulation of osmolytes, Na+/K+ ratio, sequestration in the vacuole for cellular redox homeostasis[22,62]

Share and Cite

MDPI and ACS Style

Rohila, J.S.; Edwards, J.D.; Tran, G.D.; Jackson, A.K.; McClung, A.M. Identification of Superior Alleles for Seedling Stage Salt Tolerance in the USDA Rice Mini-Core Collection. Plants 2019, 8, 472. https://doi.org/10.3390/plants8110472

AMA Style

Rohila JS, Edwards JD, Tran GD, Jackson AK, McClung AM. Identification of Superior Alleles for Seedling Stage Salt Tolerance in the USDA Rice Mini-Core Collection. Plants. 2019; 8(11):472. https://doi.org/10.3390/plants8110472

Chicago/Turabian Style

Rohila, Jai S., Jeremy D. Edwards, Gioi D. Tran, Aaron K. Jackson, and Anna M. McClung. 2019. "Identification of Superior Alleles for Seedling Stage Salt Tolerance in the USDA Rice Mini-Core Collection" Plants 8, no. 11: 472. https://doi.org/10.3390/plants8110472

APA Style

Rohila, J. S., Edwards, J. D., Tran, G. D., Jackson, A. K., & McClung, A. M. (2019). Identification of Superior Alleles for Seedling Stage Salt Tolerance in the USDA Rice Mini-Core Collection. Plants, 8(11), 472. https://doi.org/10.3390/plants8110472

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