Next Article in Journal
Genome Size Diversity and Its Impact on the Evolution of Land Plants
Previous Article in Journal
Characterization of TTN Novex Splicing Variants across Species and the Role of RBM20 in Novex-Specific Exon Splicing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GWAS Uncovers Differential Genetic Bases for Drought and Salt Tolerances in Sesame at the Germination Stage

1
Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, No.2 Xudong 2nd Road, Wuhan 430062, China
2
Centre d’Etudes Régional Pour l’Amélioration de l’Adaptation à la Sécheresse (CERAAS), Route de Khombole, Thiès, BP 3320, Senegal
3
College of Life and Environmental Sciences, Shanghai Normal University, Shanghai 200234, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 7 December 2017 / Revised: 25 January 2018 / Accepted: 8 February 2018 / Published: 14 February 2018
(This article belongs to the Section Plant Genetics and Genomics)

Abstract

:
Sesame has great potential as an industrial crop but its production is challenged by drought and salt stresses. To unravel the genetic variants leading to salinity and drought tolerances at the germination stage, genome-wide association studies of stress tolerance indexes related to NaCl-salt and polyethylene glycol-drought induced stresses were performed with a diversity panel of 490 sesame accessions. An extensive variation was observed for drought and salt responses in the population and most of the accessions were moderately tolerant to both stresses. A total of 132 and 120 significant Single Nucleotide Polymorphisms (SNPs) resolved to nine and 15 Quantitative trait loci (QTLs) were detected for drought and salt stresses, respectively. Only two common QTLs for drought and salt responses were found located on linkage groups 5 and 7, respectively. This indicates that the genetic bases for drought and salt responses in sesame are different. A total of 13 and 27 potential candidate genes were uncovered for drought and salt tolerance indexes, respectively, encoding transcription factors, antioxidative enzymes, osmoprotectants and involved in hormonal biosynthesis, signal transduction or ion sequestration. The identified SNPs and potential candidate genes represent valuable resources for future functional characterization towards the enhancement of sesame cultivars for drought and salt tolerances.

Graphical Abstract

1. Introduction

Sesame (Sesamum indicum L.) is a traditional crop considered as one of the oldest oilseed crops known by mankind [1]. It plays a very significant role in preserving food and nutritional security as well as livelihood improvement in developing regions of the world. Over the last decade, the production of sesame seeds has doubled, showing the increasing interest on this crop [2]. Sesame is especially significant in the arid and semi-arid regions where the cultivation of several major crops including maize, cotton, etc., becomes problematic. In fact, salinity and water deficit have a strong influence on agricultural development in these regions [3]. Because of the low rainfall coupled with an intense heat, drought frequently occurs and greatly impairs crop productivity. On the other hand, intense use of irrigation and applied fertilizers are the major causes of soil salinity which inhibits seed germination and plant growth, resulting in a drastic reduction of yields [4]. In contrast to most of the oilseed crops, sesame is a fairly resilient crop reported to be relatively tolerant to drought [5] and also rated moderately salt tolerant [6]. These features raise sesame as an attractive crop tailored for the challenging environmental conditions of the arid and semi-arid regions.
Despite sesame’s relative drought tolerance, severe water deficit constitutes a major impediment for its production. Especially, drought stress occurring at germination and anthesis stages has the most damaging effects on seed germination, seedling development, yield related traits and seed quality [7,8,9,10,11,12]. Similarly, detrimental effects of severe salt stress especially NaCl, on sesame seed germination, seedling growth and yield were documented by several authors [6,13,14,15,16]. Exposures to drought, salt as well as other abiotic and biotic factors are considerably decreasing sesame yield which is actually around 300–400 kg/ha [5] far from its full potential of 2500 kg/ha [17]. Therefore, enhancing sesame genotypes for drought and salinity tolerances constitutes one of the priorities in the current sesame research [2].
Several studies were performed on the screening for sources of tolerance in sesame germplasm in regard to drought and salt stresses [7,9,14,15,16,18,19,20,21,22,23,24,25,26]. It is reported that sesame accessions harbor different degrees of tolerance to these abiotic stresses by triggering various response mechanisms including the control of ion transfer from roots to leaves, the accumulation of osmotic regulators, a strong induction of antioxidative enzymes, the regulation of photosynthesis and phytohormones, etc. [18,19,20,21,22,23,24,25,26]. The diversity in sesame responses to salt and drought stresses are controlled by genetic factors. Therefore, understanding the genetic components for salt and drought tolerances by identifying the genetic loci and the candidate genes associated with these traits would be an economical, feasible and efficient way to accelerate the progress of abiotic tolerance breeding in sesame [2]. Seed germination is one the most sensitive stages of sesame growth cycle [9]. Accordingly, detecting some genetic variants associated with sesame tolerance to drought and salinity at the germination stage will be a major asset for sesame cultivation in arid and semi-arid regions.
Genome Wide Association Study (GWAS) approach is widely employed to identify relationships between molecular markers or candidate genes and traits in a given population based on linkage disequilibrium. It offers advantages over traditional linkage analyses, such as more accurate positioning and mapping, simultaneous assessments of multiple alleles at a locus and no requirement for linkage group construction [27,28]. It has been successfully applied in several crops such as alfalfa, Arabidopsis thaliana, soybean, maize, rice, asparagus bean, Brassica napus, etc., to dissect complex traits including salt and drought tolerance [29,30,31,32,33,34,35,36,37,38,39]. GWAS was also successfully applied in sesame to unravel the genetic basis of its high oil production and quality, and some other key agronomic traits by fully sequencing a population composed of 705 worldwide accessions [40]. The release of the genotypic data of this panel provided inestimable genomic information to conduct further GWAS on other important traits in sesame. This study aimed at identifying the genetic bases of salt (NaCl-induced) and drought (polyethylene glycol (PEG)-induced) tolerances at the seed germination stage in 490 diverse sesame accessions using GWAS.

2. Materials and Methods

2.1. Plant Materials

Seeds of 490 sesame accessions were obtained from the China National Gene Bank, Oil Crops Research Institute, Chinese Academy of Agricultural Sciences (Wuhan, China). The 490 accessions used for the association analysis in this study were from 33 countries of Asia, Africa, America, and Europe (Table S1), aiming to capture a high geographical diversity. All accessions were self-pollinated for four generations by growing under natural growing conditions in Sanya, Hainan province, China (109.187° E, 18.38° N, altitude 11 m). Seed lots from 2013 harvest were mainly used in this study. However, for the following accessions G656, G041, G051, G123, G128, G162, G175, G176, G189, G204, G208, G213, G215, G224, G228, G233, G242, G262, G264, G299, G301, G315, G352, G391, G400, G401, G407, G413 with limited seed quantities, seed lots from 2013 were used for controls while seed lots from 2014 were used for stress treatments.

2.2. Salt and Polyethylene Glycol Stress Treatment

A pilot experiment was conducted in order to identify the suitable concentrations of NaCl and PEG 6000 which could induce most stress conditions and could easily distinguish sesame germplasm according to their tolerance levels. For the NaCl, different concentrations were tested including 0 mM, 40 mM, 60 mM, 80 mM and 100 mM on 15 diverse accessions (Table S2). For each accession, 50 seeds were germinated in petri dishes (90 mm diameter) containing two layers of filter paper (Whatman™, Malaga, WA, Australia). Two treatments were applied: a control condition in which seeds were soaked with 10 mL deionized water; a stress condition in which seeds were soaked with 10 mL of the different concentrations of NaCl. Petri dishes were maintained in the dark into a climatic chamber (SANYO, Gallenkamp PLC, Loughborough, UK) set at 28 °C for five days. Every day, deionized water and the NaCl solutions were added into the petri dishes to compensate for evaporation in the control and stress treatments, respectively. The experiment was arranged in a completely randomized design with three replications for each accession in both treatments. In regard to the PEG 6000, similar procedures were followed as for NaCl. Different concentrations of 0%, 5%, 10%, 15%, 20%, 25% and 30% PEG 6000 were tested on 15 accessions. Then, the numbers of germinated seeds were counted and the germination percentage was estimated through the formula: (number of germinated seeds/number of plated seeds (50)) × 100.
For the 15 accessions tested, a mean germination percentage was estimated. According to the suitable concentrations identified for NaCl and PEG 6000 treatments in the pilot experiment, the 490 accessions were submitted to 60 mM NaCl and 15% PEG 6000 following experimental procedures described above. The whole experiment was repeated three times. After five days, the number of germinated seeds (GR) and the fresh weight (FW (g)) of the seedlings were recorded in both control and stress conditions.

2.3. Statistical Analysis

Based on the recorded data, the stress tolerance (ST) index was estimated as the ratio of the number of germinated seeds and fresh weight under stress conditions to the same traits under stress free conditions [41]. Therefore, ST-SGR (%) (Stress tolerance for germination percentage under salt) and ST-SFW (%) (Stress tolerance for fresh weight under salt) represent the salt tolerance index. While, ST-DGR (%) (Stress tolerance for germination percentage under PEG 6000) and ST-DFW (%) (Stress tolerance for fresh weight under PEG 6000) represent drought tolerance index.
For each trait, data were analyzed using analysis of variance (ANOVA) based on General Linear Model procedure in R3.2.0 software with the packages “Ade4” [42] and “Agricolae” [43] considering all effects as fixed. ANOVA results were considered significant at p < 0.05 and the mean comparisons were performed using the Tukey honest significant difference (HSD) test. Additionally, the least square means, standard deviation, variance, and descriptive statistics such as the coefficient of variation, range, skewness and kurtosis were estimated. Correlation coefficients among the stress indexes of PEG and NaCl treatments were calculated by Pearson’s method at a significance level of p < 0.05 using the “Corrplot” package [44]. The variation of the different traits in the control and stress conditions was represented employing the “ggplot2” package [45]. Pearson dissimilarity coefficients of stress tolerance indexes were used for hierarchical cluster analysis and the dendrogram based on Ward’s method was constructed in Minitab software (Pennsylvania State University, PA, USA).

2.4. Phenotype-Genotype Association Analysis

The association panel used in this study along with other 215 sesame accessions was previously fully re-sequenced in works of Wei et al. [40]. A total of 1,005,413 common Single Nucleotide Polymorphisms (SNPs) covering the whole genome with minor allele frequency >0.03 were used for genome wide association analysis in this study. Average Linkage Disequilibrium (LD) region was estimated to 88 kb in the whole genome of sesame [40]. The relative kinship analysis was implemented using the package SPAGeDi (Université Libre de Bruxelles, Brussels, Belgium) [46]. To get insight into the population structure of the association panel, a dendrogram based on simple matching coefficient was constructed using the software package MEGA7 (Temple University, Philadelphia, PA, USA) [47]. In addition, the program STRUCTURE 2.3.4 (Stanford University, Stanford, CA, USA) [48] was used for a Bayesian clustering analysis with the admixture model and correlated allele frequencies. Five runs were performed for each subpopulation K (1 to 10) [49]. The burn-in time and iterations for each run were set to 30,000 and 70,000, respectively and the true K was determined according to the method described by Evanno et al. [50]. Association analysis was performed for ST-SGR, ST-SFW, ST-DGR and ST-DFW using the EMMAX model [51] based on the Mixed Model. The matrix of pair-wise genetic distance derived from simple matching coefficients was used as the variance–covariance matrix of the random effect. Significance was defined at a uniform threshold of p < 5.54 × 10−7 (−log10 (p) > 6). The value r2 derived from linear regressions were calculated to examine the phenotypic variance explained (PVE) of each peak SNPs using the Minitab software (Pennsylvania State University, PA, USA). Before fitting the model, each marker was coded with the value 0 used for the reference allele and the value 1 for the alternative allele.

2.5. Mining of Potential Candidate Genes

A region of 88 kb (LD [40]) including at least three significant SNPs was considered a Quantitative trait locus (QTL). To identify the potential candidate genes around the significant peak SNPs associated with each trait, the full gene list in the QTL was searched using the sesame reference genome information [52]. The corresponding putative homologs of all genes in A. thaliana were retrieved from the database SesameFG [53] with cut off E-value ≤1 × 10−40. We retained the potential candidate genes based on two criteria: (a). the gene is annotated; (b). the gene annotation is related to abiotic stress response and/or the homolog in Arabidopsis is characterized as an abiotic stress responsive gene. Gene Ontology (GO) annotation of the potential candidate genes was carried out using Blast2GO tool v.3.1.3 (Biobam, Valencia, Spain) [54] and plotted with WEGO tool (BGI, Shenzhen, China) [55].

3. Results

3.1. Identification of Suitable Concentrations for NaCl and Polyethylene Glycol Stress Induction

From 30% to 25% PEG 6000 concentration applications, no seed could germinate, while at 20% concentration, only few seeds (12%) germinated. Therefore, 15% PEG 6000 was defined as the appropriate concentration for examining drought tolerance in the whole germplasm (Figure 1A). Concerning the NaCl concentrations, results showed that 100 mM NaCl was detrimental to the seeds. At 80 mM, only few seeds (3%) germinated; therefore, 60 mM NaCl was selected as the suitable concentration for salt stress tolerance screening in the association panel (Figure 1B). The Figures S1 and S2 present the phenotypes of seedlings from tolerant and sensitive sesame genotypes in responses to PEG 6000 (0% vs. 15%) and NaCl (0 vs. 60 mM), respectively. The growth of the tolerant genotypes was less affected under stress as compared to the sensitive genotypes.

3.2. Phenotypic Variation for Salt and Drought Tolerances in the Sesame Germplasm

An extensive phenotypic variation was displayed for drought and salt tolerance indexes in the sesame association panel (Figure 2). Normal or nearly normal distributions were observed for all the stress indexes in the mapping population. ST-DFW ranged from 7% to 276% with most of the accessions displaying 50% tolerance. Similarly, ST-DGR ranged from 2% to 266% and most accessions showed 60% tolerance. Concerning the salt stress, a similar trend was observed. ST-SFW values were between 2% and 294% with the majority of the accessions exhibiting 60% tolerance. For the ST-SGR values, most of the accessions displayed a good tolerance (mean ST-SGR value around 70%). Overall, the fresh weight of the seedlings was found to be more affected than the number of germinated seeds, although both traits were significantly (p < 0.001) reduced under drought and salt stresses (Table 1). It is worth to mention that for some accessions, the fact that different seed lots were used for controls and stressed conditions, may explain the outperformance of stressed samples over the controls.
Figure 3A,B presents the clustering patterns of the 490 accessions according to their tolerance levels to salt and drought. In total, three main clusters were observed for both stresses: the cluster of tolerant accessions is colored in blue, the one of moderately tolerant accessions is colored in black and finally the cluster of sensitive accessions is marked in carmine color. In general, the tolerant accessions and sensitive ones were less represented than the moderately tolerant accessions for drought and salt stresses (Table 2). In addition, the responses of sesame accessions to PEG and NaCl are mostly inconsistent (Figure 3C). For instance, only 27 accessions were commonly tolerant to drought and salt stresses (TT), 82 commonly sensitive to both stresses (SS) and 90 accessions identified as commonly moderately tolerant to both stresses (MM). The remaining 291 accessions displayed differential responses to drought and salt stresses (Figure 3C). Similarly, the correlation analysis of the drought and salt tolerance indexes indicated that traits from the same stress were strongly and significantly correlated whereas the correlation of traits between drought and salt were significantly weak (p < 0.001) (Figure 4). This indicates that PEG-induced drought stress and NaCl-induced salt stress lead to differential responses in sesame and are probably underlined by different genetic factors.

3.3. Genome-Wide Association Studies for Drought and Salt Tolerance Indexes

A total of 1M SNPs were employed for the GWAS analysis of drought and salt tolerance traits, resulting in a high marker density of 2700 SNPs per Mb. STRUCTURE and dendrogram showing genetic relationships revealed two recognizable groups with a very weak differentiation between them (Fst = 0.041) (Figure S3 and Table S1).
To get insight into the genetic variants associated with drought and salt tolerances at the germination stage in sesame, a GWAS was conducted for stress tolerance index on four traits (Figure 5). In regard to the PEG-induced drought stress, a total of 132 significant SNPs (−log10 (p) > 6) were uncovered, with 73 and 59 significant SNPs for ST-DFW and ST-DGR, respectively (Figure 5A,B). These SNPs were located on the linkage groups (LGs) 1, 4, 5, 7, 8, 10 and 11. The phenotypic variance explained (PVE) values of the peak SNPs were between 4.93% (SNP10174676) and 10.45% (SNP765526). Because of the dense coverage of SNPs in the whole genome, we assumed that the peak significant SNPs from a true association must be in a strong LD with the surrounding SNPs. Hence, this has led us to retain five QTLs (qDFW1.1, qDFW5.1, qDFW7.1, qDFW8.1 and qDFW10.1) and four QTLs (qDGR1.1, qDGR4.1, qDGR7.1 and qDGR11.1) for ST-DFW and ST-DGR, respectively. Between both indexes, we found two common QTLs regions (qDFW1.1 similar to qDGR1.1 and qDFW7.1 similar to qDGR7.1) located on LG1 and LG7, led by the SNPs SNP1657071 and SNP10174187, respectively.
Concerning the NaCl-induced salt stress, we identified in total 120 significant SNPs (−log10 (p) > 6) including 90 and 30 significant SNPs for ST-SFW and ST-SGR, respectively (Figure 5C,D). Their PVE values ranged from 2.92% (SNP1848856) to 7.7% (SNP11548217). These SNPs were distributed on all the LGs except for the LGs9, 13 and 14. Following the same assumption described above, we identified 12 QTLs (qSFW1.1, qSFW2.1, qSFW3.1, qSFW3.2, qSFW4.1, qSFW4.2, qSFW5.1, qSFW6.1, qSFW7.1, qSFW10.1, qSFW11.1 and qSFW12.1) for ST-SFW located on LGs 1, 2, 3, 4, 5, 6, 7, 10, 11 and 12 (Figure 5C). Similarly, three QTLs (qSGR1.1, qSGR6.1 and qSGR16.1) were retained for ST-SGR located on LGs 1, 6 and 16 (Figure 5D). In contrast to drought stress indexes, no common QTL was found for the two salt tolerance indexes.
By comparing the results from drought and salt GWAS, we found only two common QTLs (qSFW5.1 similar to qDFW5.1 and qSFW7.1 similar to qDGR7.1) associated with the stress indexes. These QTLs were linked to the SNPs SNP1848856 and SNP10174187 located on the LG5 and LG7, respectively.

3.4. Allelic Effects of the Associated SNPs on the Salt and Drought Tolerance Indexes in Sesame

The distributions of the stress tolerance indexes were further examined in individuals that carried each allele of the peak SNPs (Figure 6). We defined as favorable, the allele at the peak SNP which increases stress tolerance indexes. In general, the non-reference alleles at the peak SNPs were found to be the favorable alleles. Only the reference allele cytosine “C” at the SNP SNP23040151 significantly associated with ST-SGR, was identified as the favorable allele over the variant thymine “T” allele. In addition, there was no significant difference between the groups of accessions with the favorable alleles and the groups of accessions harboring the unfavorable alleles at the SNP SNP17546324, SNP14830572, SNP12613392, SNP4353089, SNP1848856 and SNP9099698 (Figure 6A–D).
The tolerant and sensitive groups of sesame accessions under PEG-induced drought stress (Table 2) were compared for variation and enrichment of favorable alleles at the peak SNPs associated with ST-DFW and ST-DGR. As shown in Figure 6E, the tolerant accessions harbored significantly more favorable alleles (ranging from one to seven, with a mean of three favorable alleles) than the sensitive group (ranging from zero to two, with a mean of zero favorable allele). The accession G440 harboring the highest number of favorable alleles (seven) has the best drought tolerance index (ST-GR and ST-FW).
Similarly, under NaCl-induced salt stress, there was a significant difference between the two groups (tolerant vs. sensitive). The tolerant accessions harbored more favorable alleles (ranging from one to five, with a mean of two favorable alleles) than the sensitive accessions (ranging from zero to two, with a mean of zero favorable allele). The accession G430 holds five favorable alleles and was one of the most salt tolerant genotype based on values of ST-SGR and ST-SFW.

3.5. Assigning Significant SNPs Associated with Drought and Salt Tolerance to Potential Candidate Genes

To assess the putative candidate genes associated with the significant SNPs for salt and drought tolerances in sesame, we retrieved all genes in the 88 kb window (LD region) around each peak SNP. A total of 241 and 151 genes were found for salt and drought stresses, respectively (Table 3, Table 4, Tables S3 and S4). The gene number around the peak SNPs ranged from seven to 33 for the drought tolerance index and from six to 32 for the salt stress index.
Based on the annotation information of the retrieved genes and the functions described for their homologs in A. thaliana, we subsequently retained 13 putative candidate genes associated with PEG-induced drought tolerance (Table 3 and Table S5) and 27 potential candidate genes for NaCl-induced salt tolerance in sesame (Table 4 and Table S6). These genes encode transcription factors, antioxidative enzymes, osmoprotectants and are involved in hormonal biosynthesis or ion sequestration. Six common genes (SIN_1007701, SIN_1007708, SIN_1007698, SIN_1008841, SIN_1008842 and SIN_1009337) were identified for drought and salt stresses. Gene ontology analysis of the potential candidate genes established that they are related to stress responses (Figure 7). Both drought and salt potential candidate genes were predominantly enriched in metabolic process, response to abiotic or biotic stresses and response to stress (Figure 7A,B).

4. Discussion

4.1. Drought and Salt Responses Are Governed by Different Genetic Components in Sesame

Drought and salinity stresses determine the primary cause of worldwide crop loss [56]. The germination stage is one of the crucial stages of the sesame crop growth cycle especially when it is confronted with abiotic stresses including waterlogging, drought and salt [9]. Although some QTLs and candidate genes were previously reported for waterlogging tolerance in sesame [57,58], very limited researches have been conducted concerning drought and salt stresses [2]. In this study, we observed that most of the sesame accessions were moderately tolerant to drought and salinity stresses as previously reported by several authors [5,6,14,19,21,23]. We noticed that the responses of sesame accessions to drought and salinity stresses were quite different based on stress tolerance indexes. This is further supported by the weak correlation between salt and drought tolerance indexes. Patade et al. [59,60] and Lokhande et al. [61] reported similar observations in sugarcane and Sesuvium portulacastrum, respectively. They demonstrated that PEG-induced osmotic stress and NaCl-induced salinity stress trigger differential biochemical and physiological responses. In fact, plants exposed to salt stress accumulate saline ions as an osmoticum stored in the vacuoles to prevent toxicity [62]. Additionally, they also synthesize compatible solutes such as proline, polyols, glycine betaine etc., for cellular osmotic adjustment [63]. In contrast, plants exposed to dehydration stress solely rely on the synthesis of compatible solutes for maintenance of the cell turgor (Lokhande et al. [61]).
Based on the phenotypic observations, we speculate that the genetic bases of drought and salt responses in sesame are different. Interestingly, the GWAS results corroborated well our hypothesis since, only two of the 20 identified QTLs were common to both abiotic stresses. In Arabidopsis, Kreps et al. [64] also showed that the transcriptomes involved in several abiotic stresses including cold, salt and drought were principally stimulus specific. Our results suggested that the genes controlling drought and salt stresses belong to fundamentally different genetic pathways in sesame.

4.2. GWAS Is an Effective Approach to Identify Functional SNPs and Candidate Genes for Drought and Salt Tolerances in Sesame

The modest LD decay rate in sesame and the high marker density coupled with the low genetic differentiation of our association panel are advantageous for GWAS implementation in this study. Most of the significant associations occurred in clusters of SNPs, delineating some important genomic regions governing drought and salt responses. In total, we identified 20 QTLs with modest effects for drought and salt tolerance indexes at the germination stage in sesame. The number of significant associated SNPs and their phenotypic contributions are similar to reports of salt and drought tolerances in Brassica napus, alfalfa and soybean [31,36,38]. The QTLs are distributed on 12 out of the 16 LGs of the sesame genome, suggesting that salt and drought responses are complex traits which implicate several genomic regions [56]. Moreover, these significant QTL regions harbored a very limited number of genes which is crucial to rapidly pinpoint the causative genes. Through classical bi-parental QTL mapping for plant height [65] and for yield related traits [66] in sesame, QTLs expanding on very large genomic regions and harboring hundreds of genes were discovered. Spotting the causative genes that govern these complex agronomic traits seems to be a daunting challenge. Therefore, GWAS is a more effective approach to dissect complex traits as compared to the classical bi-parental QTL mapping in sesame.
The favorable alleles at the peak SNPs improve tolerance under drought and salt treatments. We found more favorable alleles in the drought-tolerant accessions than in the salt-tolerant accessions, which may explain why sesame is naturally more tolerant to drought than salinity stress [16]. In addition, the accessions with the highest number of favorable alleles at the peak SNPs were the most tolerant. Knowing that marker-based allele pyramiding is very effective in crop improvement [67,68], we propose to confirm in future studies, whether the favorable alleles would have a positive pyramiding effect in improving drought and salt tolerances in sesame.

4.3. Discovering New Functional Genes for the Enhancement of Drought and Salt Tolerances in Sesame

To date, only few candidate genes are available for drought tolerance enhancement in sesame [10,11,12,69,70]. In contrast to drought, no gene or functional molecular marker is available for the improvement to salinity tolerance in sesame [2]. Here, we unraveled 13 and 27 potential candidate genes associated with the stress tolerance indexes for PEG-induced drought and NaCl-induced salinity, respectively. Most of the potential candidate genes are the homologs of well-described genes related to abiotic stress responses or plant development processes in the model species Arabidopsis. For example, the gene SiOPR3 (SIN_1024693) detected for drought stress in sesame is the homolog of the gene AT2G06050 (ATOPR3) in Arabidopsis. ATOPR3, an essential component of the jasmonic acid biosynthesis, is needed for the increase of abscisic acid in desiccating Arabidopsis roots which lead to a high drought tolerance [71]. Another interesting gene identified for drought is SiWRKY69 (SIN_1019661). This gene is around the SNP SNP765526 with the highest PVE (10.45%). The homolog of SIN_1019661 in Arabidopsis is AT3G58710, described as a non-memory gene functioning in response to dehydration stress. AT3G58710 was constitutively and significantly up-regulated (>6-fold the expression level in the control plants) under multiple drought stresses [72]).
Between ST-DGR and ST-DFW, we detected a common QTL (qDGR1.1 similar to qDFW1.1) located on the LG1. By examining the potential candidate genes in this genomic region, we retrieved the gene SiCCD8 (SIN_1021566), the homolog of the Arabidopsis gene AT4G32810 (ATCCD8). ATCCD8 functions as a carotenoid cleavage dioxygenase and is involved in plant growth especially the hypocotyl development under several stress conditions [73]. The peak associated SNP SNP1657071 was located in an intron of SiCCD8 with the allele ‘‘G’’ identified as the favorable allele. We infer that this intronic SNP may affect alternative splicing of the SiCCD8 mRNA, its expression level and ultimately may alter the seed germination and the seedling growth potential under drought stress in sesame.
Concerning the salt stress, we also discovered several promising genes. In the QTL qSFW2.1 led by the SNP SNP15050812, eight potential candidate genes which function in salt stress responses were uncovered. The gene SiMLP31 (SIN_1021337) is annotated as a MLP-like protein which is the homolog of Arabidopsis gene MLP31 (AT1G70840). According to Rajjou et al. [74], MLP31 is implicated in salicylic acid synthesis which markedly improves seed germination and seedling vigor under salt stress. Based on these results, we speculate that the gene SiMLP31 may hold similar attributes in sesame. Another interesting gene detected in this QTL region is SiANTH (SIN_1021330), the homolog of the gene ANTH (AT2G01600) in Arabidopsis. ANTH is a phosphatidic acid-binding protein which is recruited to membranes in response to salt stress in Arabidopsis roots [75]. Furthermore, a tandem array of several sesame Peroxidase encoding genes (SiPOD) was found in the same QTL region. These genes (SIN_1021330, SIN_1021327, SIN_1021326, SIN_1021325, SIN_1021324, SIN_1021323 and SIN_1021322) may have resulted from a whole genome duplication event [76] and are likely to play the same function since all of them have the same homolog in Arabidopsis (AT1G14550). AT1G14550 is reported to be involved in several abiotic stresses including drought, hypoxia and salt stress [77,78]. Plant peroxidases have been ascribed a variety of biological functions, including hydrogen peroxidase detoxification, lignin biosynthesis, hormonal signaling and stress response [79]. Under abiotic stresses such as salt, drought, cold, heat etc., plant accumulates excessive amount of reactive oxygen species (ROS) in cells which damages many cellular components. Peroxidase (POD) plays key roles in cellular ROS detoxification and is therefore a cardinal element of the plant antioxidant defense [80]. Based on these results, we deduce that this genomic region in the LG2 might be an important QTL for salt tolerance in sesame.
Identifying some QTLs associated with multiple desirable traits is highly sought in plant improvement strategies. Here, we found two QTLs (qSFW5.1 similar to qDFW5.1 and qSFW7.1 similar to qDGR7.1) commonly associated with drought and salt tolerance indexes which may be useful for the development of sesame varieties tolerant to both stresses. In the first QTL region located on the LG5, we uncovered three promising candidate genes. Within them, the gene SiHKT1 (SIN_1007698) is a sodium transporter crucial for plant survival under salt stress. Mounting evidence has highlighted this gene associated with ion content as an essential component in salinity tolerance in Arabidopsis and other plant species [81,82,83]. It is worth mentioning that by using GR and FW data under control condition for GWAS analysis, we were unable to detect the QTL region harboring SiHKT1 (Data not shown). Therefore, identifying SiHKT1 with the stress tolerance indexes data suggests that SiHKT1 is not involved in sesame germination but in sesame response to salt stress at the germination stage. Nonetheless, whether this gene also functions in drought response in sesame or there is another gene in this QTL which is involved in drought response, is still unclear and will need additional analysis. In the second QTL located on LG7 common to both stresses, we also detected three potential candidate genes. Within these genes, SiDREB2A (SIN_1009337) that is the homolog of the gene AtDREB2A (AT5G18450) in Arabidopsis is well characterized for its involvement in drought and salt tolerances in several species [84,85,86,87].
To capture more effective candidate genes, we will conduct in future studies, RNA-seq analyses under PEG and NaCl stresses. Also, transgenic experiments and gene-based association analysis will help to shed more light on the functions of the candidate genes as well as the genetic variants altering their expression levels. Additionally, we will transform the significant peak SNPs into allele-specific markers such as Cleaved Amplified Polymorphic Sequences or Kompetitive allele specific PCR, in order to examine their potency as diagnostic molecular markers for sesame breeding programs.

5. Conclusions

In this study, stress tolerance indexes for four traits related to PEG-induced drought stress and NaCl-induced salt stress were assessed using a worldwide panel of sesame accessions at the germinating stage. Most of the accessions were moderately tolerant to both stresses. In addition, we found that sesame accessions respond differently to drought and salt stresses. GWAS revealed 20 QTLs including nine QTLs for drought and 15 QTLs for salt, with modest phenotypic contributions. Only two common QTLs were identified for both stresses, suggesting that the genes controlling these two abiotic stresses belong to fundamentally different genetic pathways. In total, 13 potential candidate genes associated with drought responses and 27 potential candidate genes for salt response were uncovered. These genes encode transcription factors, antioxidative enzymes, osmoprotectants and are involved in hormonal biosynthesis or ion sequestration. Interestingly, most the potential candidate genes are described to be involved in responses to abiotic stresses. Altogether, the phenotyped biomaterials, the discovered functional SNPs and the potential candidate genes for drought and salt responses will be helpful for marker-assisted breeding programs aimed at enhancing salt and drought tolerances in sesame cultivars.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/9/2/87/s1, Figure S1: Responses of tolerant and sensitive sesame genotypes to PEG 6000 at the concentration 15%; Figure S2: Responses of tolerant and sensitive sesame genotypes to NaCl at the concentration 60 mM; Figure S3: Population genetics of the association panel. (A) Neighbor-joining tree of all the varieties calculated from whole-genome SNPs revealed two recognizable groups colored in red and purple; (B,C) Estimated population structure in the association assessed by STRUCTURE. Each individual is represented by a thin vertical bar, partitioned into k = 2 colored segments representing two subpopulations; Table S1: Full list of the 490 sesame accessions used in this study, their origins, breeding status, sequencing information, dendrogram grouping, salt and drought tolerance status; Table S2: List of the 15 accessions used for preliminary salt and drought screening experiments; Table S3: Full list of the genes detected around the peak SNPs for PEG-induced drought stress in the sesame association panel; Table S4: Full list of the genes detected around the peak SNPs for NaCl-induced salt stress in the sesame association panel; Table S5: List of the significant peak SNPs and annotation of associated potential candidate genes for PEG-induced drought stress; Table S6: List of the significant peak SNPs and annotation of associated potential candidate genes for NaCl-induced salt stress.

Acknowledgments

This work is financially supported by the National Science-technology Support Plan Project the National Science-technology Support Plan Project (2013BAD01B03-08), the Opening Project of the Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, P. R. China (2017007), the China Agriculture Research System (CARS-15), the Agricultural Science and Technology Innovation Project of Chinese Academy of Agricultural Sciences (CAAS-ASTIP-2013-OCRI), the National Natural Science Foundation of China (31401412) and the Fundamental Research Funds for Central Non-profit Scientific Institution (1610172014003). K.D. acknowledges the fellowship offered by the Chinese Scholarship Council (2015GXY934).

Author Contributions

X.W., Y.Z., L.W. and X.Z. conceived and designed the experiments; D.L., Y.Z., A.L., Y.Z., R.Z. and K.D. performed the experiments; K.D., L.D. and X.W. analyzed the data; K.D., A.L., L.D., Y.Z., R.Z. and X.Z. contributed reagents/materials/analysis tools; K.D. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ashri, A. Sesame Breeding. Plant Breed. Rev. 1998, 16, 179–228. [Google Scholar]
  2. Dossa, K.; Diouf, D.; Wang, L.; Wei, X.; Zhang, Y.; Niang, M.; Fonceka, D.; Yu, J.; Mmadi, M.A.; Yehouessi, L.W.; et al. The emerging oilseed crop Sesamum indicum enters the “Omics” era. Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef] [PubMed]
  3. Ramirez, R.; Gutierrez, D.; Villafane, R.; Lizaso, J.I. Salt tolerance of sesame genotypes at germination, vegetative and maturity stages. Commun. Soil Sci. Plant Anal. 2005, 36, 2405–2419. [Google Scholar] [CrossRef]
  4. Radhakrishnan, R.; Lee, I.-J. Penicillium–sesame interactions: A remedy for mitigating high salinity stress effects on primary and defense metabolites in plants. Environ. Exp. Bot. 2015, 116, 47–60. [Google Scholar] [CrossRef]
  5. Islam, F.; Gill, R.A.; Ali, B.; Farooq, M.A.; Xu, L.; Najeeb, U.; Zhou, W. Sesame. In Breeding Oilseed Crop for Sustainable Production: Opportunities and Constraints; Gupta, S.K., Ed.; Academic Press: Cambridge, MA, USA, 2016; pp. 135–147. [Google Scholar]
  6. Yousif, H.Y.; Bingham, F.T.; Yermason, D.M. Growth, mineral composition and seed oil of sesame (Sesamum indicum L.) as affected by NaCl. Soil Sci. Soc. Am. Proc. 1972, 36, 450–453. [Google Scholar] [CrossRef]
  7. Hassanzadeh, M.; Asghari, A.; Jamaati-e-Somarin, S.; Saeidi, M.; Zabihi-e-Mahmoodabad, R.; Hokmalipour, S. Effects of water deficit on drought tolerance indices of sesame (Sesamum indicum L.) genotypes in Moghan Region. Res. J. Environ. Sci. 2009, 3, 116–121. [Google Scholar] [CrossRef]
  8. Sun, J.; Rao, Y.; Le, M.; Yan, T.; Yan, X.; Zhou, H. Effects of drought stress on sesame growth and yield characteristics and comprehensive evaluation of drought tolerance. Chin. J. Oil Crop Sci. 2010, 32, 525–533. (In Chinese) [Google Scholar]
  9. Boureima, S.; Eyletters, M.; Diouf, M.; Diop, T.A.; Van Damme, P. Sensitivity of seed germination and seedling radicle growth to drought stress in sesame (Sesamum indicum L.). Res. J. Environ. Sci. 2011, 5, 557–564. [Google Scholar] [CrossRef]
  10. Dossa, K.; Niang, M.; Assogbadjo, A.E.; Cisse, N.; Diouf, D. Whole genome homology-based identification of candidate genes for drought resistance in (Sesamum indicum L.). Afr. J. Biotechnol. 2016, 15, 1464–1475. [Google Scholar]
  11. Dossa, K.; Wei, X.; Li, D.; Fonceka, D.; Zhang, Y.; Wang, L.; Yu, J.; Boshou, L.; Diouf, D.; Cissé, N.; et al. Insight into the AP2/ERF transcription factor superfamily in sesame (Sesamum indicum) and expression profiling of the DREB subfamily under drought stress. BMC Plant Biol. 2016, 16, 171. [Google Scholar] [CrossRef] [PubMed]
  12. Dossa, K.; Diouf, D.; Cissé, N. Genome-wide investigation of Hsf genes in sesame reveals their segmental duplication expansion and their active role in drought stress response. Front. Plant Sci. 2016, 7, 1522. [Google Scholar] [CrossRef] [PubMed]
  13. Cerda, A.; Bingham, F.T.; Hoffman, G.J. Interactive effect of salinity and phosphorus on sesame. Soil Sci. Soc. Am. J. 1977, 41, 915–918. [Google Scholar] [CrossRef]
  14. Koca, H.; Bor, M.; Ozdemir, F.; Turkan, I. The effect of salt stress on lipid peroxidation, antioxidative enzymes and proline content of sesame cultivars. Environ. Exp. Bot. 2007, 60, 344–351. [Google Scholar] [CrossRef]
  15. Yahya, A. Selectivity and partitioning of potassium and sodium in sesame. J. Plant Nutr. 2010, 33, 670–683. [Google Scholar] [CrossRef]
  16. Bazrafshan, A.H.; Ehsanzadeh, P. Growth, photosynthesis and ion balance of sesame (Sesamum indicum L.) genotypes in response to NaCl concentration in hydroponic solutions. Photosynthetica 2014, 52, 134–147. [Google Scholar] [CrossRef]
  17. Weiss, E.A. Oil Seed Crop, 2nd ed.; Black Well Science: Malden, MA, USA, 2000. [Google Scholar]
  18. Betram, K.; Janssens, M.J.J.; Abdalwahab, A. Breeding for drought tolerance in sesame (Sesamum indicum). In Proceedings of the Conference on Technological and Institutional Innovations for Sustainable Rural Development, Gottingen, Germany, 8–10 October 2003; p. 135. [Google Scholar]
  19. Gehlot, H.S.; Purohit, A.; Shekhawat, N.S. Metabolic changes and protein patterns associated with adaptation to salinity in Sesamum indicum cultivars. J. Cell. Mol. Biol. 2005, 4, 31–39. [Google Scholar]
  20. Boureima, S.; Oukarroum, A.; Diouf, M.; Cissé, N.; Van Damme, P. Screening for drought tolerance in mutant germplasm of sesame (Sesamum indicum) probing by chlorophyll a fluorescence. Environ. Exp. Bot. 2012, 81, 37–43. [Google Scholar] [CrossRef]
  21. Boureima, S.; Diouf, M.; Amoukou, A.I.; Damme, V.P. Screening for sources of tolerance to drought in sesame induced mutants: Assessment of indirect selection criteria for seed yield. Int. J. Pure Appl. Biosci. 2016, 4, 45–60. [Google Scholar] [CrossRef]
  22. Bahrami, H.; Razmjoo, J.; Jafari, A.O. Effect of drought stress on germination and seedling growth of sesame cultivars (Sesamum indicum L.). Int. J. AgriSci. 2012, 2, 423–428. [Google Scholar]
  23. Kadkhodaie, A.; Zahedi, M.; Razmjoo, J.; Pessarakli, M. Changes in some anti-oxidative enzymes and physiological indices among sesame genotypes (Sesamum indicum L.) in response to soil water deficits under field conditions. Acta Physiol. Plant. 2013, 36, 641–650. [Google Scholar] [CrossRef]
  24. Kadkhodaie, A.; Razmjoo, J.; Zahedi, M.; Pessarakli, M. Selecting sesame genotypes for drought tolerance based on some physiochemical traits. Agron. J. 2014, 106, 111–118. [Google Scholar] [CrossRef]
  25. Bazrafshan, A.H.; Ehsanzadeh, P. Evidence for differential lipid peroxidation and antioxidant enzyme activities in Sesamum indicum L. genotypes under NaCl salinity. J. Agric. Sci. Technol. 2016, 18, 207–222. [Google Scholar]
  26. Hussein, Y.; Amin, G.; Azab, A.; Gahin, H. Induction of drought stress resistance in (Sesamum indicum L.) plant by salicylic acid and kinetin. J. Plant Sci. 2015, 10, 128–141. [Google Scholar] [CrossRef]
  27. Flint-Garcia, S.A.; Thuillet, A.C.; Yu, J.; Pressoir, G.; Romero, S.M.; Mitchell, S.E.; Doebley, J.; Kresovich, S.; Goodman, M.M.; Buckler, E.S. Maize association population: A high-resolution platform for quantitative trait locus dissection. Plant J. 2005, 44, 1054–1064. [Google Scholar] [CrossRef] [PubMed]
  28. Huang, X.; Han, B. Natural variations and genome-wide association studies in crop plants. Annu. Rev. Plant Biol. 2014, 65, 531–551. [Google Scholar] [CrossRef] [PubMed]
  29. Mao, H.; Wang, H.; Liu, S.; Li, Z.; Yang, X.; Yan, J.; Li, J.; Tran, L.P.; Qin, F. A transposable element in a NAC gene is associated with drought tolerance in maize seedlings. Nat. Commun. 2015, 6, 8326. [Google Scholar] [CrossRef] [PubMed]
  30. 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] [PubMed]
  31. Zhang, T.; Yu, L.; Zheng, P.; Li, Y.; Rivera, M.; Main, D.; Greene, S.L. Identification of loci associated with drought resistance traits in heterozygous autotetraploid Alfalfa (Medicago sativa L.) using genome-wide association studies with genotyping by sequencing. PLoS ONE 2015, 10, e0138931. [Google Scholar] [CrossRef] [PubMed]
  32. Clauw, P.; Coppens, F.; Korte, A.; Herman, D.; Slabbinck, B.; Dhondt, S.; Van Daele, T.; De Milde, L.; Vermeersch, M.; Maleux, K.; et al. Leaf growth response to mild drought: Natural variation in Arabidopsis sheds light on trait architecture. Plant Cell 2016, 28, 2417–2434. [Google Scholar] [CrossRef] [PubMed]
  33. Julkowska, M.M.; Klei, K.; Fokkens, L.; Haring, M.A.; Schranz, E.M.; Testerink, C. Natural variation in rosette size under salt stress conditions corresponds to developmental differences between Arabidopsis accessions and allelic variation in the LRR-KISS gene. J. Exp. Bot. 2016, 67, 2127–2138. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, X.; Wang, H.; Liu, S.; Ferjani, A.; Li, J.; Yan, J.; Yang, X.; Qin, F. Genetic variation in ZmVPP1 contributes to drought tolerance in maize seedlings. Nat. Genet. 2016, 48, 1233–1241. [Google Scholar] [CrossRef] [PubMed]
  35. Xu, L.; Hu, K.; Zhang, Z.; Guan, C.; Chen, S.; Hua, W.; Li, J.; Wen, J.; Yi, B.; Shen, J.; et al. Genome-wide association study reveals the genetic architecture of flowering time in rapeseed (Brassica napus L.). DNA Res. 2016, 23, 43–52. [Google Scholar] [CrossRef] [PubMed]
  36. Zeng, A.; Chen, P.; Korth, K.; Hancock, F.; Pereira, A.; Brye, K.; Wu, C.; Shi, A. Genome-wide association study (GWAS) of salt tolerance in worldwide soybean germplasm lines. Mol. Breed. 2017, 37, 30. [Google Scholar] [CrossRef]
  37. Patishtan, J.; Hartley, T.N.; de Carvalho, R.F.F.; Maathuis, J.M. Genome wide association studies to identify rice salt-tolerance markers. Plant Cell Environ. 2017. [Google Scholar] [CrossRef] [PubMed]
  38. Tan, M.; Liao, F.; Hou, L.; Wang, J.; Wei, L.; Jian, H.; Xu, X.; Li, J.; Liu, L. Genome-wide association analysis of seed germination percentage and germination index in Brassica napus L. under salt and drought stresses. Euphytica 2017, 213, 40. [Google Scholar] [CrossRef]
  39. Liu, X.; Yu, L. Genome-wide association mapping of loci associated with plant growth and forage production under salt stress in Alfalfa (Medicago sativa L.). Front. Plant Sci. 2017, 8, 853. [Google Scholar] [CrossRef] [PubMed]
  40. Wei, X.; Liu, K.; Zhang, Y.; Feng, Q.; Wang, L.; Zhao, Y.; Li, D.; Zhao, Q.; Zhu, X.; Zhu, X.; et al. Genetic discovery for oil production and quality in sesame. Nat. Commun. 2015, 6, 8609. [Google Scholar] [CrossRef] [PubMed]
  41. Long, N.V.; Dolstra, O.; Malosetti, M.; Kilian, B.; Graner, A.; Visser, R.G.F.; van der Linden, C.G. Association mapping of salt tolerance in barley (Hordeum vulgare L.). Theor. Appl. Genet. 2013, 126, 2335–2351. [Google Scholar] [CrossRef] [PubMed]
  42. Dray, S.; Dufour, A.B.; Chessel, D. The ade4 package-II: Two-table and K-table methods. R News 2007, 7, 47–52. [Google Scholar]
  43. De Mendiburu, F. agricolae: Statistical Procedures for Agricultural Research, R Package Version 1.2-1. 2014. Available online: http://CRAN.R-project.org/package=agricolae (accessed on 15 August 2017).
  44. Wei, T.; Simko, V. corrplot: Visualization of a Correlation Matrix, R Package Version 0.77. 2016. Available online: http://CRAN.R-project.org/package=corrplot, (accessed on 15 August 2017).
  45. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2009. [Google Scholar]
  46. Hardy, O.J.; Vekemans, X. SPAGEDi: A versatile computer program to analyze spatial genetic structure at the individual or population levels. Mol. Ecol. Notes 2002, 2, 618–620. [Google Scholar] [CrossRef]
  47. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis Version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  48. Falush, D.; Stephens, M.; Pritchard, J.K. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 2003, 164, 1567–1587. [Google Scholar] [PubMed]
  49. Dossa, K.; Wei, X.; Zhang, Y.; Fonceka, D.; Yang, W.; Diouf, D.; Liao, B.; Cissé, N.; Zhang, X. Analysis of genetic diversity and population structure of sesame accessions from Africa and Asia as major centers of its cultivation. Genes 2016, 7, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [PubMed]
  51. Kang, H.; Sul, J.; Service, S.; Zaitlen, N.; Kong, S.; Freimer, N.; Sabatti, C.; Eskin, E. Variance component model to account for sample structure in genome-wide association studies. Nat. Genet. 2010, 42, 348–354. [Google Scholar] [CrossRef] [PubMed]
  52. Wang, L.; Yu, S.; Tong, C.; Zhao, Y.; Liu, Y.; Song, C.; Zhang, Y.; Zhang, X.; Wang, H.; Wang, X.; et al. Genome sequencing of the high oil crop sesame provides insight into oil biosynthesis. Genome Biol. 2014, 15, R39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wei, X.; Gong, H.; Yu, J.; Liu, P.; Wang, L.; Zhang, Y.; Zhang, X. Sesame FG: An integrated database for the functional genomics of sesame. Sci. Rep. 2017, 7, 2342. [Google Scholar] [CrossRef] [PubMed]
  54. Götz, S.; García-Gómez, J.M.; Terol, J.; Williams, T.D.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talón, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef] [PubMed]
  55. Ye, J.; Fang, L.; Zheng, H.; Zhang, Y.; Chen, J.; Zhang, Z.; Wang, J.; Li, S.; Li, R.; Bolund, L.; et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, W293–W297. [Google Scholar] [CrossRef] [PubMed]
  56. Mondini, L.; Pagnotta, M.A. Drought and salt stress in cereal. In Sustainable Agriculture Reviews; Lichtfouse, E., Goyal, A., Eds.; Springer: Cham, Switzerland, 2015; p. 16. [Google Scholar]
  57. Zhang, Y.; Wang, L.; Li, D.; Gao, Y.; Lu, H.; Zhang, X. Mapping of sesame waterlogging tolerance QTL and identification of excellent waterlogging tolerant germplasm. Sci. Agric. Sin. 2014, 47, 422–430. (In Chinese) [Google Scholar]
  58. Wang, L.; Li, D.; Zhang, Y.; Gao, Y.; Yu, J.; Wei, X.; Zhang, X. Tolerant and susceptible sesame genotypes reveal waterlogging stress response patterns. PLoS ONE 2016, 11, e0149912. [Google Scholar] [CrossRef] [PubMed]
  59. Patade, V.Y.; Bhargava, S.; Suprasanna, P. Salt and drought tolerance of sugarcane and antioxidant defense. J. Plant Interact. 2011, 6, 275–282. [Google Scholar] [CrossRef]
  60. Patade, V.Y.; Bhargava, S.; Suprasanna, P. Effects of NaCl and iso-osmotic PEG stress on growth, osmolytes accumulation and antioxidant defense in cultured sugarcane cells. Plant Cell Tissue Organ. Cult. 2012, 108, 279–286. [Google Scholar] [CrossRef]
  61. Lokhande, V.H.; Nikam, T.D.; Penna, S. Differential osmotic adjustment to iso-osmotic NaCl and PEG stress in the in vitro cultures of Sesuvium portulacastrum (L.) L. J. Crop Sci. Biotechnol. 2010, 13, 251–256. [Google Scholar] [CrossRef]
  62. Munns, R.; Tester, M. Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 2008, 59, 651–681. [Google Scholar] [CrossRef] [PubMed]
  63. De Olivera, A.B.; Alencar, N.L.M.; Gomes-Filho, E. Comparison between the Water and Salt Stress Effects on Plant Growth and Development; Intech: Rijeka, Croatia, 2013; pp. 67–94. [Google Scholar]
  64. Kreps, J.A.; Wu, Y.; Chang, H.S.; Zhu, T.; Wang, X.; Harper, J.F. Transcriptome changes for Arabidopsis in response to salt, osmotic and cold stress. Plant Physiol. 2002, 130, 2129–2141. [Google Scholar] [CrossRef] [PubMed]
  65. Wang, L.; Xia, Q.; Zhang, Y.; Zhu, X.; Zhu, X.; Li, D.; Ni, X.; Gao, Y.; Xiang, H.; Wei, X.; et al. Updated sesame genome assembly and fine mapping of plant height and seed coat color QTLs using a new high-density genetic map. BMC Genom. 2016, 17, 31. [Google Scholar] [CrossRef] [PubMed]
  66. Wu, K.; Liu, H.; Yang, M.; Tao, Y.; Ma, H.; Wu, W.; Zuo, Y.; Zhao, Y. High-density genetic map construction and QTLs analysis of grain yield-related traits in sesame (Sesamum indicum L.) based on RAD-Seq technology. BMC Plant Biol. 2014, 14, 274. [Google Scholar] [CrossRef] [PubMed]
  67. Sacco, A.; Di, M.A.; Lombardi, N.; Trotta, N.; Punzo, B.; Mari, A.; Barone, A. Quantitative trait loci pyramiding for fruit quality traits in tomato. Mol. Breed. 2013, 31, 217–222. [Google Scholar] [CrossRef] [PubMed]
  68. Zhang, B.; Li, W.; Chang, X.; Li, R.; Jing, R. Effects of favorable alleles for water soluble carbohydrates at grain filling on grain weight under drought and heat stresses in wheat. PLoS ONE 2014, 29, e102917. [Google Scholar] [CrossRef] [PubMed]
  69. Dossa, K.; Li, D.; Wang, L.; Zheng, X.; Yua, J.; Wei, X.; Fonceka, D.; Diouf, D.; Liao, B.; Cisse, N.; et al. Dynamic transcriptome landscape of Sesame (Sesamum indicum L.) under progressive drought and after rewatering. Genom. Data 2017, 11, 122–124. [Google Scholar] [CrossRef] [PubMed]
  70. Dossa, K.; Li, D.; Wang, L.; Zheng, X.; Liu, A.; Yu, J.; Wei, X.; Zhou, R.; Fonceka, D.; Diouf, D.; et al. Transcriptomic, biochemical and physio-anatomical investigations shed more light on responses to drought stress in two contrasting sesame genotypes. Sci. Rep. 2017, 7. [Google Scholar] [CrossRef] [PubMed]
  71. De Ollas, C.; Arbona, V.; Gomez-Cadenas, A. Jasmonoyl isoleucine accumulation is needed for abscisic acid build-up in roots of Arabidopsis under water stress conditions. Plant Cell Environ. 2015, 38, 2157–2170. [Google Scholar] [CrossRef] [PubMed]
  72. Ding, Y.; Liu, N.; Virlouvet, L.; Riethoven, J.; Fromm, M.; Avramova, Z. Four distinct types of dehydration stress memory genes in Arabidopsis thaliana. BMC Plant Biol. 2013, 13, 229. [Google Scholar] [CrossRef] [PubMed]
  73. Marzec, M.; Muszynska, A.; Gruszka, D. The role of strigolactones in nutrient-stress responses in plants. Int. J. Mol. Sci. 2013, 14, 9286–9304. [Google Scholar] [CrossRef] [PubMed]
  74. Rajjou, L.; Belghazi, M.; Huguet, R.; Robin, C.; Moreau, A.; Job, C.; Job, D. Proteomic investigation of the effect of salicylic acid on Arabidopsis seed germination and establishment of early defense mechanisms. Plant Physiol. 2006, 141, 910–923. [Google Scholar] [CrossRef] [PubMed]
  75. McLoughlin, F.; Arisz, S.A.; Dekker, H.L.; Kramer, G.; de Koster, C.G.; Haring, M.A.; Munnik, T.; Testerink, C. Identification of novel candidate phosphatidic acid-binding proteins involved in the salt-stress response of Arabidopsis thaliana roots. Biochem. J. 2013, 450, 573–581. [Google Scholar] [CrossRef] [PubMed]
  76. Yu, J.; Ke, T.; Tehrim, S.; Sun, F.; Liao, B.; Hua, W. PTGBase: An integrated database to study tandem duplicated genes in plants. Database 2015, 2015, bav017. [Google Scholar] [CrossRef] [PubMed]
  77. Ma, S.; Gong, Q.; Bohnert, H.J. Dissecting salt stress pathways. J. Exp. Bot. 2006, 57, 1097–1107. [Google Scholar] [CrossRef] [PubMed]
  78. Yang, C.; Hsu, F.; Li, J.; Wang, N.; Shih, M. The AP2/ERF transcription factor AtERF73/HRE1 modulates ethylene responses during hypoxia in Arabidopsis. Plant Physiol. 2011, 156, 202–212. [Google Scholar] [CrossRef] [PubMed]
  79. Gao, C.; Wang, Y.; Liu, G.; Wang, C.; Jiang, J.; Yang, C. Cloning of ten peroxidase (POD) genes from Tamrix hispida and characterization of their responses to abiotic stress. Plant Mol. Biol. Rep. 2010, 28, 77. [Google Scholar] [CrossRef]
  80. Farooq, M.; Wahid, A.; Kobayashi, N.; Fujita, D.; Basra, S.M.A. Plant drought stress: Effects, mechanisms and management. Agron. Sustain. Dev. 2009, 29, 185–212. [Google Scholar] [CrossRef]
  81. Shkolnik-Inbar, D.; Adler, G.; Bar-Zvi, D. ABI4 down regulates expression of the sodium transporter HKT1;1 in Arabidopsis roots and affects salt tolerance. Plant J. 2013, 73, 993–1005. [Google Scholar] [CrossRef] [PubMed]
  82. Ali, A.; Park, H.C.; Aman, R.; Ali, Z.; Yun, D. Role of HKT1 in Thellungiella salsuginea, a model extremophile plant. Plant Signal. Behav. 2013, 8, e25196. [Google Scholar] [CrossRef] [PubMed]
  83. Ali, A.; Raddatz, N.; Aman, R.; Kim, S.; Park, H.C.; Jan, M.; Baek, D.; Khan, I.U.; Oh, D.; Lee, S.Y.; et al. A single amino-acid substitution in the sodium transporter HKT1 associated with plant salt tolerance. Plant Physiol. 2016, 171, 2112–2126. [Google Scholar] [CrossRef] [PubMed]
  84. Dubouzet, J.G.; Sakuma, Y.; Ito, Y.; Kasuga, M.; Dubouzet, E.G.; Miura, S.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt- and cold-responsive gene expression. Plant J. 2003, 33, 751–763. [Google Scholar] [CrossRef] [PubMed]
  85. Qin, F.; Sakuma, Y.; Tran, L.P.; Maruyama, K.; Kidokoro, S.; Fujita, M.; Umezawa, T.; Sawano, Y.; Miyazono, K.; Tanokura, M.; et al. Arabidopsis DREB2A-interacting proteins function as RING E3 ligases and negatively regulate plant drought stress-responsive gene expression. Plant Cell 2008, 20, 1693–1707. [Google Scholar] [CrossRef] [PubMed]
  86. Sakuma, Y.; Maruyama, K.; Osakabe, Y.; Qin, F.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. Plant Cell 2006, 18, 1292–1309. [Google Scholar] [CrossRef] [PubMed]
  87. Schmidt, R.; Mieulet, D.; Hubberten, H.; Obata, T.; Hoefgen, R.; Fernie, A.R.; Fernie, A.R.; Fisahn, J.; San Segundo, B.; Guiderdoni, E.; et al. SALT-RESPONSIVE ERF1 regulates reactive oxygen species-dependent signaling during the initial response to salt stress in rice. Plant Cell 2013, 25, 2115–2131. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Pilot experiment showing the mean germination percentage of 15 sesame accessions under various concentrations of NaCl and PEG 6000. (A) PEG 6000; (B) NaCl. PEG: polyethylene glycol.
Figure 1. Pilot experiment showing the mean germination percentage of 15 sesame accessions under various concentrations of NaCl and PEG 6000. (A) PEG 6000; (B) NaCl. PEG: polyethylene glycol.
Genes 09 00087 g001
Figure 2. Phenotypic variation of the stress tolerance indexes for the measured traits. (A) Stress tolerance for fresh weight under PEG 6000 (ST-DFW) (%); (B) Stress tolerance for germination percentage under PEG 6000 (ST-DGR) (%); (C) Stress tolerance for fresh weight under salt (ST-SFW) (%); (D) Stress tolerance for germination percentage under salt (ST-SGR) (%).
Figure 2. Phenotypic variation of the stress tolerance indexes for the measured traits. (A) Stress tolerance for fresh weight under PEG 6000 (ST-DFW) (%); (B) Stress tolerance for germination percentage under PEG 6000 (ST-DGR) (%); (C) Stress tolerance for fresh weight under salt (ST-SFW) (%); (D) Stress tolerance for germination percentage under salt (ST-SGR) (%).
Genes 09 00087 g002
Figure 3. Classification of 490 sesame accessions according to their levels of tolerance/sensitivity to PEG and NaCl stresses. (A) Dendrogram and heat map hierarchical clustering of the accessions according to their tolerance levels to NaCl. (B) Dendrogram and heat map hierarchical clustering of the accessions according to their tolerance levels to PEG. The cluster colored in carmine gathers the sensitive accessions, the blue one is for the tolerant accessions and the black cluster represents the moderately tolerant accessions. (C) Comparison of the responses of the accessions to PEG and NaCl. The bars showed the proportion of accessions according to their responses to both stresses. The bars colored in blue, black and carmine represent the accessions which were commonly tolerant, moderately tolerant and sensitive to both stresses, respectively. The bars colored in gray represent the genotypes displaying contrasting responses to PEG and NaCl treatments. T = Tolerant, M = Moderately tolerant, S = Sensitive, based on the classification of the dendrogram. The first letter refers to PEG while the second letter refers to NaCl.
Figure 3. Classification of 490 sesame accessions according to their levels of tolerance/sensitivity to PEG and NaCl stresses. (A) Dendrogram and heat map hierarchical clustering of the accessions according to their tolerance levels to NaCl. (B) Dendrogram and heat map hierarchical clustering of the accessions according to their tolerance levels to PEG. The cluster colored in carmine gathers the sensitive accessions, the blue one is for the tolerant accessions and the black cluster represents the moderately tolerant accessions. (C) Comparison of the responses of the accessions to PEG and NaCl. The bars showed the proportion of accessions according to their responses to both stresses. The bars colored in blue, black and carmine represent the accessions which were commonly tolerant, moderately tolerant and sensitive to both stresses, respectively. The bars colored in gray represent the genotypes displaying contrasting responses to PEG and NaCl treatments. T = Tolerant, M = Moderately tolerant, S = Sensitive, based on the classification of the dendrogram. The first letter refers to PEG while the second letter refers to NaCl.
Genes 09 00087 g003
Figure 4. Correlation analysis of the stress tolerance indexes. ***, ** Significant difference at p < 0.001 and p < 0.01, respectively.
Figure 4. Correlation analysis of the stress tolerance indexes. ***, ** Significant difference at p < 0.001 and p < 0.01, respectively.
Genes 09 00087 g004
Figure 5. Manhattan plots and QQ plots of genome-wide association studies using the Mixed model for drought and salt tolerance indexes. (A) ST-DFW; (B) ST-DGR; (C) ST-SFW; (D) ST-SGR. The significant trait-associated QTLs commonly identified for the two indexes of the same stress are highlighted with dotted red lines and those commonly identified for drought and salt stresses are highlighted with dotted purple lines. QQ plot: Quantile-Quantile plots.
Figure 5. Manhattan plots and QQ plots of genome-wide association studies using the Mixed model for drought and salt tolerance indexes. (A) ST-DFW; (B) ST-DGR; (C) ST-SFW; (D) ST-SGR. The significant trait-associated QTLs commonly identified for the two indexes of the same stress are highlighted with dotted red lines and those commonly identified for drought and salt stresses are highlighted with dotted purple lines. QQ plot: Quantile-Quantile plots.
Genes 09 00087 g005
Figure 6. Allelic effects of the peak SNPs on NaCl-induced salt stress and PEG-induced drought tolerance indexes. (A) ST-SFW (%); (B) ST-SGR; (C) ST-DFW (%); (D) ST-DGR; (E) Mean number of favorable alleles identified in the tolerant and sensitive accessions to drought and salt stresses. *, **, *** significant difference at p < 0.05, p < 0.01, p < 0.001, respectively. For each SNP, the first allele is the reference allele and the second allele is the variant. Gray bars represent the favorable allele while the green bars represent the unfavorable allele. Error bars represent standard deviation.
Figure 6. Allelic effects of the peak SNPs on NaCl-induced salt stress and PEG-induced drought tolerance indexes. (A) ST-SFW (%); (B) ST-SGR; (C) ST-DFW (%); (D) ST-DGR; (E) Mean number of favorable alleles identified in the tolerant and sensitive accessions to drought and salt stresses. *, **, *** significant difference at p < 0.05, p < 0.01, p < 0.001, respectively. For each SNP, the first allele is the reference allele and the second allele is the variant. Gray bars represent the favorable allele while the green bars represent the unfavorable allele. Error bars represent standard deviation.
Genes 09 00087 g006
Figure 7. Biological process enrichment of the putative candidate genes. (A) PEG-induced drought stress; (B) NaCl-induced salt stress. GO: Gene Ontology.
Figure 7. Biological process enrichment of the putative candidate genes. (A) PEG-induced drought stress; (B) NaCl-induced salt stress. GO: Gene Ontology.
Genes 09 00087 g007
Table 1. Phenotypic variation of GR and FW under control and stress conditions.
Table 1. Phenotypic variation of GR and FW under control and stress conditions.
TraitTreatment (T)MeanSD RangeVarianceSkewnessKurtosisA TA × T
GRDistilled water409.4033–5088.40−1.753.65*********
PEG stress25.244.600–50112.42−1.461.94
Salt stress24.233.460–50109.51−1.411.66
FW (g)Distilled water1.490.521.08–2.930.26−0.28−0.30*********
PEG stress0.660.180.0.5–1.610.08−0.270.00
Salt stress0.70.160.05–1.920.220.230.05
A = Accessions, T = Treatments, GR = number of germinated seeds and FW = Fresh weight of the seedlings. *** Significant difference at p < 0.001.
Table 2. Distribution of the accessions according to their tolerance levels to salt and drought stresses.
Table 2. Distribution of the accessions according to their tolerance levels to salt and drought stresses.
StressTolerance LevelsNumber of AccessionsPercentage (%)Mean ST-GRMean ST-FW
Sensitive195400.7666.28
SaltModerate230470.8781.83
Tolerant65131.17100.62
Sensitive176360.7850.31
PEGModerate207420.9456.75
Tolerant107221.0691.14
Table 3. SNPs significantly associated with PEG-induced drought tolerance indexes at the germination stage in sesame.
Table 3. SNPs significantly associated with PEG-induced drought tolerance indexes at the germination stage in sesame.
TraitLGQTLPeak SNP–log10 (p)Ref
Base
SNP
Base
MAFGenes in LD PVE (%)Candidate Gene IDGene Name
ST-DFWLG1qDFW1.1SNP16570718.63AG0.06247.23SIN_1021558SiCCD8
SIN_1021566SiEMF1
LG5qDFW5.1SNP18488568.55CT0.07246.89SIN_1007701SiGPAT3
SIN_1007708SiAGL37
SIN_1007698SiHKT1
LG7qDFW7.1SNP101746766.54GA0.04274.93SIN_1008841SiGDH2
SIN_1008842SiCYP76C7
SIN_1009337SiDREB
LG8qDFW8.1SNP76552612.88AG0.033310.45SIN_1019660SiRABA1D
SIN_1019661SiWRKY69
LG10qDFW10.1SNP106972008.85GA0.0677.67NANA
ST-DGRLG1qDGR1.1SNP16570718.63AG0.06247.23SIN_1021566SiEMF1
SIN_1021558SiCCD8
LG4qDGR4.1SNP60894198.66TG0.07107.22SIN_1001572SiGRV2
LG7qDGR7.1SNP101741879.01TC0.03266.61SIN_1008842SiCYP76C7
SIN_1009337SiDREB
SIN_1008841SiGDH2
LG11qDGR11.1SNP90996986.59GC0.06255.3SIN_1024695SiGRF5
SIN_1024693SiOPR3
PVE (%): Phenotypic variance explained, refbase: reference allele, SNPbase: non-reference allele, MAF: Minor Allele Frequency. LG: Linkage Group.
Table 4. SNPs significantly associated with NaCl-induced salt tolerance indexes at the germination stage in sesame.
Table 4. SNPs significantly associated with NaCl-induced salt tolerance indexes at the germination stage in sesame.
TraitLGQTLPeak SNP–log10 (p)REF
BASE
SNP
Base
MAFGenes in LD PVE (%)Candidate Gene IDGene Name
ST-SFWLG1qSFW1.1SNP11482556.16CT0.05203.93SIN_1021624SiLHCB6
LG2qSFW2.1SNP150508127.45GT0.05218.03SIN_1021337SiMLP31
SIN_1021330SiANTH
SIN_1021327SiPOD
SIN_1021326SiPOD
SIN_1021325SiPOD
SIN_1021324SiPOD
SIN_1021323SiPOD
SIN_1021322SiPOD
LG3qSFW3.1SNP175463246.66GA0.03164.04SIN_1015378SiHSFA1
LG3qSFW3.2SNP101735689.92TG0.08158.32SIN_1017475SiDUF538
LG4qSFW4.1SNP53851009.74TC0.0567.45SIN_1018894SiCC-NBS-LRR
LG4qSFW4.2SNP16620878.04CT0.0686.52SIN_1008463SiUDG
LG5qSFW5.1SNP18488568.55CT0.07242.92SIN_1007701SiGPAT3
SIN_1007708SiAGL37
SIN_1007698SiHKT1
LG6qSFW6.1SNP126133929.5TC0.0666.52SIN_1018616SiNAC43
LG7qSFW7.1SNP101098616.45CT0.0495.88SIN_1008841SiGDH2
SIN_1008842SiCYP76C7
SIN_1009337SiDREB
LG10qSFW10.1SNP148305726.95TC0.03234.91SIN_1026087SiCP24
LG11qSFW11.1SNP115482176.83AC0.04327.7SIN_1013032NA
LG12qSFW12.1SNP57550306.37TC0.06166.68SIN_1006749SiWRKY14
SIN_1006753SiLSD1
ST-SGRLG1qSGR1.1SNP53157067.33TC0.06105.91SIN_1026318SiXXT5
LG6qSGR6.1SNP230401518.18CT0.03248.45SIN_1022410SiXTH15
LG16qSGR16.1SNP43530898.98GC0.07107.45SIN_1003799SiG6PD1

Share and Cite

MDPI and ACS Style

Li, D.; Dossa, K.; Zhang, Y.; Wei, X.; Wang, L.; Zhang, Y.; Liu, A.; Zhou, R.; Zhang, X. GWAS Uncovers Differential Genetic Bases for Drought and Salt Tolerances in Sesame at the Germination Stage. Genes 2018, 9, 87. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9020087

AMA Style

Li D, Dossa K, Zhang Y, Wei X, Wang L, Zhang Y, Liu A, Zhou R, Zhang X. GWAS Uncovers Differential Genetic Bases for Drought and Salt Tolerances in Sesame at the Germination Stage. Genes. 2018; 9(2):87. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9020087

Chicago/Turabian Style

Li, Donghua, Komivi Dossa, Yanxin Zhang, Xin Wei, Linhai Wang, Yujuan Zhang, Aili Liu, Rong Zhou, and Xiurong Zhang. 2018. "GWAS Uncovers Differential Genetic Bases for Drought and Salt Tolerances in Sesame at the Germination Stage" Genes 9, no. 2: 87. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9020087

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