Next Article in Journal
Involvement of Molecular Mechanisms between T/B Cells and IL-23: From Palmoplantar Pustulosis to Autoimmune Diseases
Next Article in Special Issue
Mechanism and Utilization of Ogura Cytoplasmic Male Sterility in Cruciferae Crops
Previous Article in Journal
Exploration of Diazaspiro Cores as Piperazine Bioisosteres in the Development of σ2 Receptor Ligands
Previous Article in Special Issue
Fine Mapping and Gene Analysis of restorer-of-fertility Gene CaRfHZ in Pepper (Capsicum annuum L.)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 15
10.3390/ijms23158260
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

A Large-Scale Genomic Association Analysis Identifies the Candidate Genes Regulating Salt Tolerance in Cucumber (Cucumis sativus L.) Seedlings

by
Dongrang Liu
†,
Shaoyun Dong
†,
Han Miao
,
Xiaoping Liu
,
Caixia Li
,
Jianan Han
,
Shengping Zhang
* and
Xingfang Gu
*
Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(15), 8260; https://doi.org/10.3390/ijms23158260
Submission received: 8 June 2022 / Revised: 19 July 2022 / Accepted: 19 July 2022 / Published: 27 July 2022
(This article belongs to the Special Issue Vegetable Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Salt stress seriously restricts plant growth and development, affects yield and quality, and thus becomes an urgent problem to be solved in cucumber stress resistance breeding. Mining salt tolerance genes and exploring the molecular mechanism of salt tolerance could accelerate the breeding of cucumber germplasm with excellent salt stress tolerance. In this study, 220 cucumber core accessions were used for Genome-Wide Association Studies (GWAS) and the identification of salt tolerance genes. The salinity injury index that was collected in two years showed significant differences among the core germplasm. A total of seven loci that were associated with salt tolerance in cucumber seedlings were repeatedly detected, which were located on Chr.2 (gST2.1), Chr.3 (gST3.1 and gST3.2), Chr.4 (gST4.1 and gST4.2), Chr.5 (gST5.1), and Chr.6 (gST6.1). Within these loci, 62 genes were analyzed, and 5 candidate genes (CsaV3_2G035120, CsaV3_3G023710, CsaV3_4G033150, CsaV3_5G023530, and CsaV3_6G009810) were predicted via the functional annotation of Arabidopsis homologous genes, haplotype of extreme salt-tolerant accessions, and qRT-PCR. These results provide a guide for further research on salt tolerance genes and molecular mechanisms of cucumber seedlings.

1. Introduction

Cucumber (Cucumis sativus L.) is one of the most important vegetables worldwide, and according to the FAOSTAT (https://www.fao.org/faostat/zh/#home (accessed on 17 January 2022)), the global total grown area and production of cucumber was 2.25 million hectares and 90.35 million tons, respectively, in 2020. However, cucumber is very sensitive to salt stress, which is one of the important limiting factors in cucumber cultivation [1,2].
In general, salt stress restricts plant physiological activity through osmotic and toxic effects, depending on the salt concentration [3]. Previous studies have shown that plants have evolved multiple complex mechanisms to adapt to salt stress in environment, including antioxidant systems, osmotic regulation, transcriptional regulation, and signal transduction, etc. [4,5]. Salt stress will cause significant damage to cucumber seedlings, which is mainly manifested in obvious symptoms of chlorosis and the decline of various growth indicators of the plant and root system, and the damage is significantly increased with the increase of salt concentration and treatment time [6,7]. Therefore, mining genes that are related to plant salt tolerance is an important and complex process. Researches have shown that cucumber salt tolerance is a quantitative trait that is controlled by multiple genes and is easily affected by environmental factors [7,8]. Liu, et al. [7] evaluated the salt tolerance of a RILs population in cucumber seedlings and detected a major QTL on chromosome 6 using salt injury index. As this is currently the only reported QTL controlling the salt tolerance of cucumber seedlings, it is urgent to identify more loci and genes to elucidate the mechanism of cucumber salt tolerance.
Recently, linkage mapping and genome-wide association studies have been widely used in mining genes that are related to vegetable phenotypic variation [9]. Hou, et al. [10] conducted GWAS analysis on 212 diverse melon (Cucumis melo L.) accession seedlings using chilling injury index and identified 11 low-temperature tolerance genes. Khanzada, et al. [11] identified the phenotypes of 228 Brassica accession seedlings under drought stress and 85 drought tolerance genes were identified by GWAS analysis in rapeseed (Brassica napus). Zhang, et al. [12] used GWAS analysis of 83 non-heading Chinese cabbage accessions to identify 17 loci that were associated with turnip mosaic virus (TuMV) resistance. In cucumber, a large number of loci and genes that are related to cucumber growth phenotypes and resistance traits have been unearthed through GWAS analysis. Wei, et al. [13] used the heat injury index of 86 cucumber core accessions to conduct GWAS analysis and detected seven GWAS loci (gHII4.1, gHII4.2, gHII5.1, gHII5.2, gHII6.1, gHII6.2, and gHII7.1). Wang, et al. [14] repeatedly detected four GWAS loci (gLTS1.1, gLTS3.1, gLTS4.1, and gLTS5.1) by using the chilling injury index of 172 cucumber core accessions seedlings. Lee, et al. [15] conducted GWAS analysis of cucumber powdery mildew resistance using 264 accessions from all over the world and predicted two novel candidate genes (Csa5G453160 and Csa5G471070). Liu, et al. [16,17] used the core germplasm to predict five downy mildew-related candidate genes (Csa1G575030, Csa2G060360, Csa4G064680, Csa5G606470, and Csa7G004020) and eight candidate genes that were associated with powdery mildew (Csa2G030010, Csa2G030020, Csa2G030030, Csa3G414050, Csa4G022350, Csa5G48880, Csa5G603950, and Csa5G603960), respectively. However, no study on the investigation of cucumber salt tolerance via GWAS has been reported.
Here, GWAS was performed using the salt injury index of 220 cucumber core accessions, the loci that were repeatedly detected in two years, that were further investigated, and candidate genes were mined through haplotype and qRT-PCR analyses. These genes could lay the foundation for the elaboration of cucumber salt stress response mechanism and the molecular improvement of excellent salt-tolerant cucumber varieties.

2. Results

2.1. Genetic Diversity of Salt Tolerance in Cucumber Germplasm

The 220 accessions were evaluated for salt tolerance two times, in June 2020 (Salt_2020) and May 2021 (Salt_2021). A total of three replicates were set for each experiment, each of which had 5 plants. At 14 days old, the seedlings were exposed to salt stress by 200 mmol·L−1 NaCl. After 7 days of salt treatment, the salt injury was divided into six scales, based on the chlorosis and wilting degree of the first true leaf and cotyledons (Figure 1a). For each accession, the salt injury index of each experiment was calculated by the formula on the basis of the six scales and shown in Table S2. The phenotypic data of the two treatments showed that the salt injury index of different accessions was significantly different. The average salt injury index was 31.18 and 28.80, with the coefficient of variation of 41.35% and 47.84%, respectively (Table 1). As shown in Table 1 and Figure 1c, the normal distribution test (p > 0.05) and the violin plot showed that the salt injury index follow the normal distribution pattern, and the salt injury index correlation among the two experiments was significant (Figure 1c,d).
According to difference in the salt injury index, 220 accessions were divided into 5 groups (I, II, III, IV, and V) by the Ward method (Figure 2). Group Ⅰ had 58 highly salt-tolerant accessions, Group Ⅱ had 45 salt-tolerant accessions, Group Ⅲ had 66 intermediate salt-tolerant accessions, Group Ⅳ had 31 salt-sensitive accessions, and Group Ⅴ had 20 highly salt-sensitive accessions (Table S3). The results of cluster analysis were further analyzed in combination with the ecotype of the germplasm. The 220 core accessions included 1 European-Asia ecotype, 5 American slicing ecotypes, 6 Xishuangbanna ecotypes, 12 Japan ecotypes, 16 American processing ecotypes, 17 India ecotypes, 18 North China-European hybrid ecotypes, 30 South China ecotypes, 31 European ecotypes, and 84 North China ecotypes. It was difficult to conclude the salt tolerance of the Xishuangbanna ecotype, American slicing ecotype, and the European-Asia ecotype, since only several lines were included. As such, we analyzed the ecotypes which contained more than 10 accessions. The results showed that most of the Japan ecotypes showed high salt-tolerance; the European, North China-European hybrid and North China ecotypes exhibited intermediate salt-tolerance; while the Indian ecotypes were more salt-sensitive (Figure 3).

2.2. Genome-Associated Analysis of Salt Tolerance in Seedlings

Using the salt injury index data of the two experiments in 2020 and 2021, genome-wide association analysis (GWAS) was performed and 7 genetic loci (gST2.1, gST3.1, gST3.2, gST4.1, gST4.2, gST5.1, and gST6.1) that were repeatedly detected. These 7 loci were distributed over five chromosomes (Figure 4), which were considered as novel loci for salt tolerance in cucumber seedlings.

2.3. Candidate Gene Analysis for the Novel Loci

In order to identify the potential candidate genes in these loci, the 100-kb region around each peak of SNPs was used for further analysis, and genes in each locus are shown in Table S4. To identify the candidate genes, functional annotation of the orthologs in Arabidopsis, haplotype of extreme salt-tolerant accessions, and qRT-PCR were analyzed. According to the results of the clustering analysis, Group Ⅰ (58 highly salt-tolerant accessions) and Group Ⅴ (20 highly salt-sensitive accessions) were used for haplotype analysis. Based on the phenotypic data and ecotype of accessions, salt-tolerant accessions (R45, R63, R78, and CG104) and salt-sensitive accessions (R59, R132, CG13, and CG14) were used for expression analysis. At 14 days old, the seedlings were exposed to salt stress by 200 mmol·L−1 NaCl, with a total of 4 time points, namely 0 h, 10 h, 48 h, and 96 h post-salt treatment. The first true leaf of each seedling was harvested as the sample for qRT-PCR. Through the above analysis, we preliminarily predicted five candidate genes that were related to salt tolerance in gST2.1, gST3.1, gST4.2, gST5.1, and gST6.1 (Table S5).
For gST2.1, the 100 kb (23,409,202~23,509,202 bp) region on Chr.2 was analyzed using pairwise LD correlations (Figure 5a). According to CuGenDB, 11 annotated genes are located in this region, and there was one gene (CsaV3_2G035120) that was related to salt stress (Figure 5b and Table S5). CsaV3_2G035120 (encoding an isocitrate lyase) had a SNP in the second exon region that led to the mutation of glycine to alanine (Figure 5c). In the highly salt-tolerant and highly salt-sensitive accessions, 49 out of the 58 highly salt-tolerant accessions (84.48%) carried the G haplotype, while 16 out of the 20 highly salt-sensitive (80%) accessions had the alternate C haplotype (Figure 5d). qRT-PCR analysis revealed that the relative expression level of this gene was significantly up-regulated in highly salt-tolerant accessions (R45 and R63) and highly salt-sensitive accessions (R132, CG13, and CG14) after salt treatment (Figure 5e).
For gST3.1, the 100 kb (21,496,729~21,596,729 bp) region on Chr.3 was analyzed using pairwise LD correlations (Figure 6a). According to CuGenDB, six annotated genes are located in this region, and there was one gene (CsaV3_3G023710) that was related to salt stress (Figure 6b and Table S5). CsaV3_3G023710 (encoding an ankyrin repeat family protein) had two SNPs in the first exon region that led to the mutation of lysine to glutamate and isoleucine to methionine (Figure 6c). In the highly salt-tolerant and highly salt-sensitive accessions, 57 out of the 58 highly salt-tolerant accessions (98.28%) carried the AA haplotype, while 15 out of the 20 highly salt-sensitive (75%) accessions had the alternate GG haplotype (Figure 6d). qRT-PCR analysis revealed that the relative expression level of this gene was significantly up-regulated in the highly salt-tolerant accessions (R45 and R63) and the highly salt-sensitive accessions (R132, CG13, and CG14) after salt treatment (Figure 6e).
For gST3.2 and gST4.1, we found two genes (CsaV3_3G030250 and CsaV3_4G027680) that were related to salt stress according to CuGenDB. There were nine SNPs on CsaV3_3G030250 and one SNP on CsaV3_4G027680, which were located in the CDS (coding sequence) region that led to amino acid variation (Figures S1c and S2c). But according to the haplotype analysis, all of these SNPs did not have any differential change between the highly salt-tolerant and highly salt-sensitive accessions (Figures S1d and S2d). As such, the candidate genes in loci gST3.2 and gST4.1 had not been predicted.
For gST4.2, the 100 kb (23,423,770~23,523,770 bp) region on Chr.4 was analyzed using pairwise LD correlations (Figure 7a). According to CuGenDB, 13 annotated genes are located in this region, and there was one gene (CsaV3_4G033150) that was related to salt stress (Figure 7b and Table S5). CsaV3_4G033150 (encoding a phosphoenolpyruvate carboxylase protein) had one SNP in the fourth exon region that led to the mutation of leucine to valine (Figure 7c). In the highly salt-tolerant and highly salt-sensitive accessions, 51 out of the 58 highly salt-tolerant accessions (87.93%) carried the C haplotype, while 10 out of the 20 highly salt-sensitive (50%) accessions had the alternate G haplotype (Figure 7d). qRT-PCR analysis revealed that after salt treatment, the relative expression level of this gene was significantly down-regulated in the highly salt-tolerant accessions (R45, R63, R78, and CG104), while was significantly up-regulated in the highly salt-sensitive accessions (R132 and CG14) (Figure 7e).
For gST5.1, the 100 kb (17,864,273~17,964,273 bp) region on Chr.5 was analyzed using pairwise LD correlations (Figure 8a). According to CuGenDB, eight annotated genes are located in this region, and there was one gene (CsaV3_5G023530) that was related to salt stress (Figure 8b and Table S5). CsaV3_5G023530 (encoding a PHD finger protein) had one SNP in the 19th exon region that led to the mutation of glutamate to aspartate (Figure 8c). In the highly salt-tolerant and highly salt-sensitive accessions, 24 out of the 58 highly salt-tolerant accessions (58.62%) carried the G haplotype, while all highly salt-sensitive (100%) accessions had the alternate C haplotype (Figure 8d). qRT-PCR analysis revealed that after salt treatment, the relative expression level of this gene was significantly up-regulated in the highly salt-tolerant accessions (R45 and R63), while was significantly down-regulated in the highly salt-sensitive accessions (R59 and CG14) (Figure 8e).
For gST6.1, the 100 kb (7,916,377~8,016,377 bp) region on Chr.6 was analyzed using pairwise LD correlations (Figure 9a). According to CuGenDB, 12 annotated genes are located in this region, and there was one gene (CsaV3_6G009810) that was related to salt stress (Figure 9b and Table S5). CsaV3_6G009810 (encoding a pesticidal crystal cry8Ba protein) had three SNPs in the CDS (coding sequence) region that led to the mutation of aspartate to glutamate, histidine to glutamine, and methionine to isoleucine (Figure 9c). In the highly salt-tolerant and highly salt-sensitive accessions, all the highly salt-tolerant accessions (100%) carried the GTT haplotype, while 13 out of the 20 highly salt-sensitive (65%) accessions had the alternate AGG haplotype (Figure 9d). qRT-PCR analysis revealed that after salt treatment, the relative expression level of this gene was significantly up-regulated in the highly salt-tolerant accessions (R45 and R63) and highly salt-sensitive accessions (R132, CG13, and CG14), while was significantly down-regulated in the highly salt-tolerant accessions (R78 and CG104) (Figure 9e).

3. Discussion

The acquisition of phenotypic data is the basis of GWAS. At present, there are many studies on the identification of salt tolerance phenotypic data of cucumber, involving morphological indicators, physiological and biochemical indicators, and cytological indicators [18]. We obtained effective phenotypic data by using the salt injury index classification survey [7]. It can efficiently obtain germplasm salt injury index data that can be used for subsequent analysis. As such, we identified the salt tolerance of different accessions by referring to this method (Figure 1a). Meanwhile, we conducted the experiments in two different environments (Shouguang in 2020 and Beijing in 2021) to ensure the accuracy of the results. The results showed that the two phenotypic data in this study have significant correlations, demonstrating the reliability of the data and the applicability of this method (Table 1 and Figure 1b).
The stress resistance of accessions in different ecotype groups was different. Cao, et al. [19] believed that the salt tolerance of South China ecotype cucumber germplasm was generally better than that of North China cucumber germplasm at the germination stage. Wei, et al. [13] identified the heat injury index of 86 cucumber core accession seedlings and found that the Indian ecotype and Xishuangbanna ecotype accessions were mostly highly heat-tolerant. Wang, et al. [14] used the chilling injury index of 87 cucumber core accession seedlings as an indicator and found that most of the North China ecotype and Japan ecotype accessions were tolerant to low temperature, while most of the European ecotype and India ecotype accessions were sensitive to low temperature. In this study, 220 core accessions were used to identify the salt tolerance of seedlings and we found that Japan ecotype and Indian ecotype accessions were generally highly salt-tolerant and highly salt-sensitive, respectively (Figure 3).
In recent years, great progress has been made in revealing the key genes of plant salt tolerance through GWAS in a variety of crops. Baxter, et al. [20] used 337 Arabidopsis accessions to conduct GWAS analysis of Na+ accumulation in leaves under salt stress and identified the key locus of Arabidopsis salt tolerance AtHKT1;1 via combined genetic complementation and gene expression studies. Patishtan, et al. [21] used 306 rice accessions as materials and used Na+ and K+ content in the tissues as indictors for GWAS analysis and predicted 70 candidate genes. Do, et al. [22] detected a GWAS locus that was associated with four traits, including leaf scorch degree, chlorophyll content ratio, leaf sodium content, and leaf chlorine content using 305 soybean accessions. Oyiga, et al. [23] used 150 wheat accessions to conduct GWAS analysis with Na+ and K+ contents in leaves and predicted 22 candidate genes for salt tolerance with combined gene expression analysis. Sun, et al. [24] measured the relative survival rate and salt tolerance of 713 cotton accessions under salt stress and predicted six candidate genes by combining GWAS analysis and expression analysis. In this study, gST2.1, gST3.1, gST3.2, gST4.1, gST4.2, gST5.1, and gST6.1 were new salt tolerance loci in cucumber through the identification of salt tolerance of core germplasm and GWAS analysis, which was also the first application of GWAS analysis in mining salt tolerance genes in cucumber (Figure 4 and Table S4).
The mechanism of cucumber response to salt stress is a complex and dynamic process. Through homologous gene cloning and expression analysis, it was found that the main mechanisms of salt stress regulation in Arabidopsis also play an important role in cucumber, such as SOS (salt-overly-sensitive) and ABA [18]. Currently, a QTL locus that is associated with salt tolerance in cucumber has been reported [7]. In addition, a large number of gene families that might be related to salt tolerance of cucumber have been identified in recent years, such as CDPK [25], SOD [26], AQP [27], bHLH [28], WRKY [29], GRX [30], and CAMTA [31]. In this study, five candidate genes that were associated with salt stress were predicted from these seven loci (Table S5). CsaV3_2G035120 encodes a glyoxylate cycle enzyme isocitrate lyase (ICL) that is involved in salt tolerance. Transcriptome analysis showed that rice isocitrate lyase gene (OsICL) was induced by salt stress, and the overexpression of rice calmodulin gene (OsCAM1-1) enhanced OsICL expression by salt stress [32]. The overexpression of OsICL in Arabidopsis AtICL knockout mutant by Yuenyong, et al. [33] further supports the role of ICL in plant salt tolerance through the glyoxylate cycle. Our study found the expression of CsaV3_2G035120 was also significantly up-regulated by salt stress in cucumber core accessions (R45, R63, R132, CG13, and CG14) (Figure 5e).
CsaV3_3G023710 encodes an ankyrin repeat family protein and the homologous gene in Arabidopsis is ITN1 (AT3G12360). ITN1 encodes a protein with an ankyrin motif and transmembrane domains that is involved in salt tolerance and may act through abscisic acid (ABA) signaling pathways and promote reactive oxygen species (ROS) production in Arabidopsis [34]. The expression of CsaV3_3G023710 was also significantly up-regulated by salt stress in highly salt-sensitive cucumber core accessions (R132, CG13, and CG14) (Figure 6e).
CsaV3_4G033150 encodes Arabidopsis phosphoenolpyruvate carboxylase (PEPC) proteins. Sanchez, et al. [35] found that PEPC genes show differential expression in Arabidopsis organs in response to salt and drought. Zhong, et al. [36] found the expression of the PEPC gene in cucumber leaves increased under NaCl stress. We also found that CsaV3_4G033150 expression was significantly up-regulated in two highly salt-sensitive accessions (R132 and CG14), while significantly down-regulated in all the highly salt-tolerant accessions (Figure 7e).
CsaV3_5G023530 encodes a plant homologous domain (PHD) finger family protein. PHD finger proteins are a family of important zinc finger transcription factors that are responsible for regulating the transcription and chromatin state and responding to various stresses, and this family of genes has been reported to be associated with salt tolerance in many plants, such as Arabidopsis thaliana [37], Moso bamboo [38], Brassica rapa [39], and Solanum tuberosum [40]. Zheng, et al. [41] found that the expression of CsaV3_5G023530 homologous gene PRE2 (AT1G77800) in Arabidopsis was significantly increased under salt treatment, which may be involved in ABA and salt responses in Arabidopsis. We also found that the expression level of CsaV3_5G023530 changed significantly in several cucumber accessions under salt stress, especially in the highly salt-sensitive accession CG13, which showed an extremely significant increase (Figure 8e).
CsaV3_6G009810 is predicted to encode a pesticidal crystal cry8Ba protein in CuGenDB, but more information in NCBI shows that it encodes a peptide of unknown function with no similarity to any known sequence to date. AtRASD1 (AT5G48310), a homologue gene of CsaV3_6G009810 in Arabidopsis, has been shown to activate abiotic stress-related responses such as salt and drought in Arabidopsis by encoding an unknown protein that is involved in ABA-dependent signal transduction pathways [42]. In this study, there were significant changes in the expression level of CsaV3_6G009810 in all the highly salt-tolerant and -sensitive cucumber accessions (Figure 9e).
In fact, the change trend of gene expression in different accessions may be variable. We found that the changes in the expression of five genes in these accessions were diverse, and the expression of the same gene was up-regulated in some of the highly salt-tolerant or -sensitive accessions, while down-regulated in others. For four highly salt-tolerant accessions, the expression of CsaV3_2G035120 was up-regulated in two accessions, while down-regulated in one accession; the expression of CsaV3_3G023710 was up-regulated in two accessions, while down-regulated in one accession; the expression of CsaV3_6G009810 was up-regulated in one accession, while down-regulated in two accessions. For four highly salt-sensitive accessions, the expression of CsaV3_2G035120 was significantly up-regulated by salt stress in three accessions, while significantly down-regulated in one accession; the expression of CsaV3_4G033150 was up-regulated in two accessions, while down-regulated in two accessions; the expression of CsaV3_5G023530 was up-regulated in one accession, while down-regulated in two accessions; the expression of CsaV3_6G009810 was up-regulated in three accessions, while down-regulated in one accession. In addition, we also found that the expression of CsaV3_5G023530 was down-regulated in two cucumber accessions under salt stress, while its homologous gene in Arabidopsis was significantly up-regulated. We think that the possible reason for this result is not only that the gene expression is variable, but also that the regulation mechanism of this gene in cucumber is different from that in Arabidopsis. In any case, the up-regulation or down-regulation of gene expression is strong evidence that these genes respond to salt stress in cucumber.
To sum up, CsaV3_2G035120, CsaV3_3G023710, CsaV3_4G033150, CsaV3_5G023530, and CsaV3_6G009810 are predicted to be involved in the salt stress response in cucumber. The next research will focus on the functional validation of these candidate genes via overexpression and CRISPR in cucumber. This work could provide some references for elucidating the mechanism of cucumber salt stress response and breeding salt-tolerant cucumber cultivars. In addition, the highly salt-tolerant and highly salt-sensitive accessions that were identified in this study can be used as experimental materials for the follow-up gene cloning and salt-tolerant cucumber cultivars breeding, such as selecting accessions to construct RILs populations for QTL mapping.

4. Materials and Methods

4.1. Plant Materials

The plant materials with a total of 220 accessions that were used in this study are from the cucumber research group of the Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, China. The code and ecotype information of 220 accessions can be found in Table S1.

4.2. Evaluation of Salt Tolerance Ability and Data Analysis

The salt treatment experiments were carried out in two different environments, including Shouguang (36°51′ N, 118°50′ E) in June 2020 and Beijing (39°9′ N, 116°3′ E) in May 2021. A total of three replicates were set for each experiment, each of which had 5 plants. The cucumbers were cultivated in 54 cm × 28 cm × 9 cm soil pots which had sufficient nursery substrate, with trays for water collection at the bottom. The nursery substrate that was used in each experiment was the same, ideal for cucumber growth and did not affect the salt stress responses. The average ambient temperature in the greenhouse during the experiments was 26.65 ℃ (Shouguang) and 25.05 ℃ (Beijing), respectively. At 14 days old, the seedlings were exposed to salt stress by adding 2 L of 200 mmol·L−1 NaCl into the trays regularly to ensure that each seedling was saturated with the same amount of salt solution. The first true leaf and cotyledon showed different degrees of chlorosis and wilting after 7 days of salt treatment. The salt treatment and the calculation of the salt injury index were carried as described by Liu, et al. [7]. The salt injury index was calculated as below: salt injury index = S × ∑(s × n)/N. Where “S” represents the highest salt injury rating scale, “s” represents the salt injury grade, “n” represents the number of plants in the salt injury grade, and “N” represents the total number of plants. For each of the cucumber seedlings that were grown under 200 mmol·L−1 NaCl, the phenotypic data were collected for the degree of salt injury. The salt symptoms within the 220 cucumber accessions were divided into six scales. WPS Office v.11.1 (Kingsoft, Beijing, China) was used to calculate salt injury index.
Variance analysis and cluster analysis of the phenotypic data were implemented using SPSS v.19.0 (International Business Machines Corporation, Armonk, NY, USA). The cluster analysis used the Ward method and the Euclidean distance method for measurement.

4.3. Genome-Wide Association Study (GWAS) and Linkage-Disequilibrium (LD) Analysis

The re-sequenced genome data and SNPs of these 220 core accessions are available in the NCBI Database (PRJNA171718 and PRJNA831637). A general linear model (GLM) was used for the association tests, with an estimated relatedness matrix as a covariate [43]. GWAS was conducted, and the genome-wide lowest p-value was recorded. The 5% lowest tail was taken from the 200 recorded minimal p-values as the threshold for genome-wide significance.
The software Plink [44] was used to calculate the LD coefficient (r2) between the pairwise high-quality SNPs and the results were used to estimate LD decay [16]. Manhattan plots and the LD heatmaps were drawn by qqman package in R environment [45].

4.4. Identification of QTLs and Candidate Gene Analysis

The regions of the significant trait-associated SNPs were identified as candidate regions of QTLs by repeated loci of two sets of the salt injury index. The cucumber genome database website CuGenDB (http://cucurbitgenomics.org/ (accessed on 2 August 2021)) was used to retrieve candidate resistance genes that were located within a 100-kb interval, from 50 kb upstream and downstream of each peak SNPs from the candidate regions [46]. Gene function was analyzed and predicted by BLAST in the CuGenDB, Swiss-Prot (https://www.uniprot.org/ (accessed on 2 August 2021)), TAIR (https://www.arabidopsis.org/index.jsp (accessed on 2 August 2021)) and Gene Ontology databases (http://amigo.geneontology.org/amigo/landing (accessed on 2 August 2021)) [47]. Furthermore, highly salt-tolerant and highly salt-sensitive core accessions that were obtained by cluster analysis were used as materials, and the genes were analyzed by combination of germplasm haplotype analysis and qRT-PCR to determine the candidate genes that were related to salt tolerance.

4.5. RNA Extraction and qRT-PCR Verification

Using four highly salt-tolerant accessions (R45, R63, R78, and CG104) and four highly salt-sensitive accessions (R59, R132, CG13, and CG14) as experimental materials, a total of 4 time points were set, namely 0 h, 10 h, 48 h, and 96 h after salt treatment. There were three replicates for each time point. The first true leaf of each seedling was harvested and immediately put into liquid nitrogen, and stored at −80 ℃. The TaKaRa MiniBEST Universal RNA Extraction Kit (TakaRa, Kusatsu, Japan) was used to extract the total RNA of each sample that was collected at the different time points. Reverse transcription of the extracted RNA into cDNA was performed with the UEIris II RT-PCR System (Biodee, Beijing, China). According to the CDS sequence of the candidate gene on the reference genome, primer 6.0 software was used to design the primers (Table S6). qRT-PCR was performed with ChamQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China), using the CFX96 Real-Time System (Bio-Rad, Hercules, CA, USA). Analysis of the relative expression data of the candidate genes was performed using the 2−ΔΔCt method [48].

5. Conclusions

We reported 58 highly salt-tolerant accessions and 20 highly salt-sensitive accessions from 220 cucumber core accessions. A total of seven loci, namely gST2.1, gST3.1, gST3.2, gST4.1, gST4.2, gST5.1, and gST6.1, that regulate salt tolerance in cucumber seedlings, were repeatedly detected. Furthermore, five candidate genes (CsaV3_2G035120, CsaV3_3G023710, CsaV3_4G033150, CsaV3_5G023530, and CsaV3_6G009810) within these loci that are involved in the salt stress response were predicted through gene analysis, haplotype, and qRT-PCR analyses. This study not only provides a reference for the subsequent excavation of candidate genes and the molecular mechanism of salt tolerance in cucumber, but also provides a theoretical basis for cultivating new high-yield and high-quality cucumber varieties with highly salt-tolerant, which will accelerate the cultivation of salt-tolerant cucumber varieties.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23158260/s1.

Author Contributions

Conceptualization, X.G., S.Z. and H.M; methodology, D.L., C.L. and J.H.; investigation, D.L. and S.D.; data curation, X.L.; writing—original draft preparation, D.L.; writing—review and editing, D.L. and S.D.; visualization, D.L.; supervision, X.G., S.Z. and H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Key-Area Research and Development Program of Shandong Province (2021LZGC016), Key-Area Research and Development Program of Guangdong Province (2020B020220001), the Earmarked Fund for Modern Agro-industry Technology Research System (CARS-23), Science and Technology Innovation Program of the Chinese Academy of Agricultural Science (CAAS-ASTIP-IVFCAAS), Central Public-Interest Scientific Institution Basal Research Fund (No.Y2017PT52), and the key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture and Rural Affairs, P.R. China.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors express their appreciation to the reviewers for their valuable comments on this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baysal, G.; Tipirdamaz, R. The effect of salinity on lipid peroxidation and some antioxidant enzyme activities in two cucumber cultivars. Acta Hortic. 2007, 729, 199–203. [Google Scholar] [CrossRef]
  2. Zhang, T.; Shi, Z.; Zhang, X.; Zheng, S.; Wang, J.; Mo, J. Alleviating effects of exogenous melatonin on salt stress in cucumber. Sci. Hortic. 2020, 262, 109070. [Google Scholar] [CrossRef]
  3. Shu, S.; Chen, L.; Lu, W.; Sun, J.; Guo, S.; Yuan, Y.; Li, J. Effects of exogenous spermidine on photosynthetic capacity and expression of Calvin cycle genes in salt-stressed cucumber seedlings. J. Plant Res. 2014, 127, 763–773. [Google Scholar] [CrossRef]
  4. Liang, W.; Ma, X.; Wan, P.; Liu, L. Plant salt-tolerance mechanism: A review. Biochem. Biophis. Res. Commun. 2018, 495, 286–291. [Google Scholar] [CrossRef]
  5. Gupta, B.; Huang, B.R. Mechanism of Salinity Tolerance in Plants: Physiological, Biochemical, and Molecular Characterization. Int. J. Genom. 2014, 2014, 701596. [Google Scholar] [CrossRef] [PubMed]
  6. Shen, J.; Jiang, J.; Tian, S.; Liu, X. Effects of NaCl Stress on Growth and Cell Membrane Permeability of Cucumber Seedlings. Guizhou Agric. Sci. 2016, 8, 19–24. [Google Scholar] [CrossRef]
  7. Liu, D.; Dong, S.; Bo, K.; Miao, H.; Li, C.; Zhang, Y.; Zhang, S.; Gu, X. Identification of QTLs Controlling Salt Tolerance in Cucumber (Cucumis sativus L.) Seedlings. Plants 2021, 10, 85. [Google Scholar] [CrossRef]
  8. Li, Y.; Si, L. Genetic analysis of salt-tolerance of cucumber at seedling stage. J. Northwest A F Univ. 2011, 39, 127–131. [Google Scholar] [CrossRef]
  9. Simko, I.; Jia, M.; Venkatesh, J.; Kang, B.-C.; Weng, Y.; Barcaccia, G.; Lanteri, S.; Bhattarai, G.; Foolad, M.R. Genomics and Marker-Assisted Improvement of Vegetable Crops. Crit. Rev. Plant Sci. 2021, 40, 303–365. [Google Scholar] [CrossRef]
  10. Hou, J.; Zhou, Y.; Gao, L.; Wang, Y.; Yang, L.; Zhu, H.; Wang, J.; Zhao, S.; Ma, C.; Sun, S.; et al. Dissecting the Genetic Architecture of Melon Chilling Tolerance at the Seedling Stage by Association Mapping and Identification of the Elite Alleles. Front. Plant Sci. 2018, 9, 1577. [Google Scholar] [CrossRef] [Green Version]
  11. Khanzada, H.; Wassan, G.M.; He, H.; Mason, A.S.; Keerio, A.A.; Khanzada, S.; Faheem, M.; Solangi, A.M.; Zhou, Q.; Fu, D.; et al. Differentially evolved drought stress indices determine the genetic variation of Brassica napus at seedling traits by genome-wide association mapping. J. Adv. Res. 2020, 24, 447–461. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, R.; Liu, C.; Song, X.; Sun, F.; Xiao, D.; Wei, Y.; Hou, X.; Zhang, C. Genome-wide association study of turnip mosaic virus resistance in non-heading Chinese cabbage. 3 Biotech 2020, 10, 363. [Google Scholar] [CrossRef] [PubMed]
  13. Wei, S.; Zhang, S.; Bo, K.; Wang, W.; Miao, H.; Dong, S.; Gu, X.; Zhang, S. Evaluation and Genome-wide Association Study (GWAS) of Seedling Thermotolerance in Cucumber Core Germplasm. J. Plant Genet. Resour. 2019, 20, 1223–1231. [Google Scholar] [CrossRef]
  14. Wang, W.; Song, Z.; Bo, K.; Dong, S.; Wei, S.; Miao, H.; Li, J.; Zhang, S.; Gu, X. Evaluation and Genome-wide Association Study (GWAS) of Low-temperature Tolerance at Seedling Stage in Cucumber Core Germplasm. J. Plant Genet. Resour. 2019, 20, 1606–1612. [Google Scholar] [CrossRef]
  15. Lee, H.Y.; Kim, J.G.; Kang, B.C.; Song, K. Assessment of the Genetic Diversity of the Breeding Lines and a Genome Wide Association Study of Three Horticultural Traits Using Worldwide Cucumber (Cucumis spp.) Germplasm Collection. Agronomy 2020, 10, 1736. [Google Scholar] [CrossRef]
  16. Liu, X.; Lu, H.; Liu, P.; Miao, H.; Bai, Y.; Gu, X.; Zhang, S. Identification of Novel Loci and Candidate Genes for Cucumber Downy Mildew Resistance Using GWAS. Plants 2020, 9, 1659. [Google Scholar] [CrossRef]
  17. Liu, X.; Gu, X.; Lu, H.; Liu, P.; Miao, H.; Bai, Y.; Zhang, S. Identification of Novel Loci and Candidate Genes for Resistance to Powdery Mildew in a Resequenced Cucumber Germplasm. Genes 2021, 12, 584. [Google Scholar] [CrossRef]
  18. Liu, D.; Dong, S.; Miao, H.; Bo, K.; Zhang, S.; Gu, X. Research Progress on Genetic Breeding of Cucumber Tolerance for Salt Stress. China Veg. 2021, 7, 14–23. [Google Scholar] [CrossRef]
  19. Cao, Q.; Li, L.; Kong, S.; Qiu, A.; Zhang, Y.; Sun, X. Evaluation of Salt Tolerance of Newly-bred Greenhouse Cucumber Cultivars at Germination Stage. J. Plant Genet. Resour. 2014, 15, 546–553. [Google Scholar] [CrossRef]
  20. Baxter, I.; Brazelton, J.N.; Yu, D.; Huang, Y.S.; Lahner, B.; Yakubova, E.; Li, Y.; Bergelson, J.; Borevitz, J.O.; Nordborg, M.; et al. A coastal cline in sodium accumulation in Arabidopsis thaliana is driven by natural variation of the sodium transporter AtHKT1;1. PLoS Genet. 2010, 6, e1001193. [Google Scholar] [CrossRef] [Green Version]
  21. 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] [PubMed]
  22. Do, T.D.; Vuong, T.D.; Dunn, D.; Clubb, M.; Valliyodan, B.; Patil, G.; Chen, P.; Xu, D.; Nguyen, H.T.; Shannon, J.G. Identification of new loci for salt tolerance in soybean by high-resolution genome-wide association mapping. BMC Genom. 2019, 20, 318. [Google Scholar] [CrossRef] [PubMed]
  23. Oyiga, B.C.; Sharma, R.C.; Baum, M.; Ogbonnaya, F.C.; Leon, J.; Ballvora, A. Allelic variations and differential expressions detected at quantitative trait loci for salt stress tolerance in wheat. Plant Cell Environ. 2018, 41, 919–935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Sun, Z.; Li, H.; Zhang, Y.; Li, Z.; Ke, H.; Wu, L.; Zhang, G.; Wang, X.; Ma, Z. Identification of SNPs and Candidate Genes Associated With Salt Tolerance at the Seedling Stage in Cotton (Gossypium hirsutum L.). Front. Plant Sci. 2018, 9, 1011. [Google Scholar] [CrossRef] [Green Version]
  25. Xu, X.; Liu, M.; Lu, L.; He, M.; Qu, W.; Xu, Q.; Qi, X.; Chen, X. Genome-wide analysis and expression of the calcium-dependent protein kinase gene family in cucumber. Mol. Genet. Genom. 2015, 290, 1403–1414. [Google Scholar] [CrossRef] [PubMed]
  26. Zhou, Y.; Hu, L.; Wu, H.; Jiang, L.; Liu, S. Genome-Wide Identification and Transcriptional Expression Analysis of Cucumber Superoxide Dismutase (SOD) Family in Response to Various Abiotic Stresses. Int. J. Genom. 2017, 2017, 7243973. [Google Scholar] [CrossRef] [Green Version]
  27. Zhu, Y.; Yang, L.; Liu, N.; Yang, J.; Zhou, X.; Xia, Y.; He, Y.; He, Y.; Gong, H.; Ma, D.; et al. Genome-wide identification, structure characterization, and expression pattern profiling of aquaporin gene family in cucumber. BMC Plant Biol. 2019, 19, 345. [Google Scholar] [CrossRef] [Green Version]
  28. Li, J.; Wang, T.; Han, J.; Ren, Z. Genome-wide identification and characterization of cucumber bHLH family genes and the functional characterization of CsbHLH041 in NaCl and ABA tolerance in Arabidopsis and cucumber. BMC Plant Biol. 2020, 20, 272. [Google Scholar] [CrossRef]
  29. Chen, C.; Chen, X.; Han, J.; Lu, W.; Ren, Z. Genome-wide analysis of the WRKY gene family in the cucumber genome and transcriptome-wide identification of WRKY transcription factors that respond to biotic and abiotic stresses. BMC Plant Biol. 2020, 20, 443. [Google Scholar] [CrossRef]
  30. Yang, Y.; Xue, W.; Chen, P.; Yuan, X.; Li, X.; Zhang, T.; Chen, S. Identification and expression analyzes of CC-type glutaredoxin in cucumber (Cucumis sativus L.) under abiotic stress. Sci. Hortic. 2021, 289, 110417. [Google Scholar] [CrossRef]
  31. Gao, R.; Luo, Y.; Yun, F.; Wu, X.; Wang, P.; Liao, W. Genome-Wide Identification, Expression Profile, and Alternative Splicing Analysis of CAMTA Family Genes in Cucumber (Cucumis sativus L.). Agronomy 2021, 11, 1827. [Google Scholar] [CrossRef]
  32. Yuenyong, W.; Chinpongpanich, A.; Comai, L.; Chadchawan, S.; Buaboocha, T. Downstream components of the calmodulin signaling pathway in the rice salt stress response revealed by transcriptome profiling and target identification. BMC Plant Biol. 2018, 18, 335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Yuenyong, W.; Sirikantaramas, S.; Qu, L.; Buaboocha, T. Isocitrate lyase plays important roles in plant salt tolerance. BMC Plant Biol. 2019, 19, 472. [Google Scholar] [CrossRef] [Green Version]
  34. Sakamoto, H.; Matsuda, O.; Iba, K. ITN1, a novel gene encoding an ankyrin-repeat protein that affects the ABA-mediated production of reactive oxygen species and is involved in salt-stress tolerance in Arabidopsis thaliana. Plant J. 2008, 56, 411–422. [Google Scholar] [CrossRef] [PubMed]
  35. Sanchez, R.; Flores, A.; Cejudo, F. Arabidopsis phosphoenolpyruvate carboxylase genes encode immunologically unrelated polypeptides and are differentially expressed in response to drought and salt stress. Planta 2006, 223, 901–909. [Google Scholar] [CrossRef] [PubMed]
  36. Zhong, M.; Yuan, Y.; Shu, S.; Sun, J.; Guo, S.; Yuan, R.; Tang, Y. Effects of exogenous putrescine on glycolysis and Krebs cycle metabolism in cucumber leaves subjected to salt stress. Plant Growth Regul. 2016, 79, 319–330. [Google Scholar] [CrossRef]
  37. Zang, D.; Li, H.; Xu, H.; Zhang, W.; Zhang, Y.; Shi, X.; Wang, Y. An Arabidopsis Zinc Finger Protein Increases Abiotic Stress Tolerance by Regulating Sodium and Potassium Homeostasis, Reactive Oxygen Species Scavenging and Osmotic Potential. Front. Plant Sci. 2016, 7, 1272. [Google Scholar] [CrossRef]
  38. Gao, Y.; Liu, H.; Wang, Y.; Li, F.; Xiang, Y. Genome-wide identification of PHD-finger genes and expression pattern analysis under various treatments in moso bamboo (Phyllostachys edulis). Plant Physiol. Biochem. 2018, 123, 378–391. [Google Scholar] [CrossRef]
  39. Alam, I.; Liu, C.; Ge, H.; Batool, K.; Yang, Y.; Lu, Y. Genome wide survey, evolution and expression analysis of PHD finger genes reveal their diverse roles during the development and abiotic stress responses in Brassica rapa L. BMC Genom. 2019, 20, 773. [Google Scholar] [CrossRef] [Green Version]
  40. Qin, M.; Luo, W.; Zheng, Y.; Guan, H.; Xie, X. Genome-wide identification and expression analysis of the PHD-finger gene family in Solanum tuberosum. PLoS ONE 2019, 14, e0226964. [Google Scholar] [CrossRef]
  41. Zheng, K.; Wang, Y.; Wang, S. The non-DNA binding bHLH transcription factor Paclobutrazol Resistances are involved in the regulation of ABA and salt responses in Arabidopsis. Plant Physiol. Biochem. 2019, 139, 239–245. [Google Scholar] [CrossRef] [PubMed]
  42. Podia, V.; Milioni, D.; Martzikou, M.; Haralampidis, K. The role of Arabidopsis thaliana RASD1 gene in ABA-dependent abiotic stress response. Plant Biol. 2018, 20, 307–317. [Google Scholar] [CrossRef] [PubMed]
  43. Bo, K.; Wei, S.; Wang, W.; Miao, H.; Dong, S.; Zhang, S.; Gu, X. QTL mapping and genome-wide association study reveal two novel loci associated with green flesh color in cucumber. BMC Plant Biol. 2019, 19, 243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  45. Turner, S.D. qqman: An R package for visualizing GWAS results using Q-Q and manhattan plots. bioRxiv 2014. [Google Scholar] [CrossRef]
  46. Yao, L.; Li, Y.; Ma, C.; Tong, L.; Du, F.; Xu, M. Combined genome-wide association study and transcriptome analysis reveal candidate genes for resistance to Fusarium ear rot in maize. J. Integr. Plant Biol. 2020, 62, 1535–1551. [Google Scholar] [CrossRef]
  47. Shi, L.; Yang, Y.; Xie, Q.; Miao, H.; Bo, K.; Song, Z.; Wang, Y.; Xie, B.; Zhang, S.; Gu, X. Inheritance and QTL mapping of cucumber mosaic virus resistance in cucumber (Cucumis Sativus L.). PLoS ONE 2018, 13, e0200571. [Google Scholar] [CrossRef]
  48. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Ijms 23 08260 g001 550
Figure 1. Phenotypic characterization of salt tolerance in 14 day old cucumber seedlings that were treated with 200 mmol·L−1 NaCl for seven days. A total of three replicates were set for each experiment, with a total of 15 replicates per accession. (a) Salt injury of cucumber seedlings was divided into 6 grades. 0, the cotyledons and the first true leaf show no symptoms of salt injury; 1, the cotyledons are slightly yellowed; 3, the cotyledons and the first true leaf are significantly yellowed; 5, the cotyledons are yellowed and wilted, and the first true leaf is significantly yellowed and wilted; 7, the cotyledons are dehydrated, the large areas of first true leaf are yellowed and wilted; 9, both the cotyledons and the first true leaf are yellowed and dehydrated. (b) Violin plot depicting the phenotypic distribution of salt injury index in the two treatments. (c) Frequency distribution and Pearson rank correlations of the mean salt injury index in two experiments. *** indicate significance at p < 0.001. (d) Heatmap depicting the salt injury index in two treatments.
Figure 1. Phenotypic characterization of salt tolerance in 14 day old cucumber seedlings that were treated with 200 mmol·L−1 NaCl for seven days. A total of three replicates were set for each experiment, with a total of 15 replicates per accession. (a) Salt injury of cucumber seedlings was divided into 6 grades. 0, the cotyledons and the first true leaf show no symptoms of salt injury; 1, the cotyledons are slightly yellowed; 3, the cotyledons and the first true leaf are significantly yellowed; 5, the cotyledons are yellowed and wilted, and the first true leaf is significantly yellowed and wilted; 7, the cotyledons are dehydrated, the large areas of first true leaf are yellowed and wilted; 9, both the cotyledons and the first true leaf are yellowed and dehydrated. (b) Violin plot depicting the phenotypic distribution of salt injury index in the two treatments. (c) Frequency distribution and Pearson rank correlations of the mean salt injury index in two experiments. *** indicate significance at p < 0.001. (d) Heatmap depicting the salt injury index in two treatments.
Ijms 23 08260 g001
Ijms 23 08260 g002 550
Figure 2. Cluster analysis of salt tolerance of core accession seedlings based on salt injury index. The I, II, III, IV, and V indicates Group I, Group II, Group III, Group IV, and Group V, respectively.
Figure 2. Cluster analysis of salt tolerance of core accession seedlings based on salt injury index. The I, II, III, IV, and V indicates Group I, Group II, Group III, Group IV, and Group V, respectively.
Ijms 23 08260 g002
Ijms 23 08260 g003 550
Figure 3. The box plot depicting the salt injury index of different ecotypes. The I, II, III, IV, and V indicates Group I, Group II, Group III, Group IV, and Group V, respectively. The dot represents the outlier. The cross represents the mean.
Figure 3. The box plot depicting the salt injury index of different ecotypes. The I, II, III, IV, and V indicates Group I, Group II, Group III, Group IV, and Group V, respectively. The dot represents the outlier. The cross represents the mean.
Ijms 23 08260 g003
Ijms 23 08260 g004 550
Figure 4. Genome-wide association analysis (GWAS) Manhattan plots of salt tolerance in two experiments: (a) Salt_2020 and (b) Salt_2021. The different colors represent the different chromosomes. The red dashed line represents the significance threshold. The red arrow indicates the position of the strong peaks that were repeatedly detected in the two experiments.
Figure 4. Genome-wide association analysis (GWAS) Manhattan plots of salt tolerance in two experiments: (a) Salt_2020 and (b) Salt_2021. The different colors represent the different chromosomes. The red dashed line represents the significance threshold. The red arrow indicates the position of the strong peaks that were repeatedly detected in the two experiments.
Ijms 23 08260 g004
Ijms 23 08260 g005 550
Figure 5. Identification of the candidate gene that was related to salt tolerance for the locus gST2.1. (a) Manhattan plot and LD heatmap of gST2.1. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of 11 functional annotated genes were predicted in the gST2.1 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) Gene structure of CsaV3_2G035120. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of C and G haplotype in highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_2G035120 in highly salt-sensitive accessions (R59, R132, CG13 and CG14) and highly salt-tolerant accessions (R45, R63, R78 and CG104) after salt treatment (0 h, 10 h, 48 h and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Figure 5. Identification of the candidate gene that was related to salt tolerance for the locus gST2.1. (a) Manhattan plot and LD heatmap of gST2.1. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of 11 functional annotated genes were predicted in the gST2.1 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) Gene structure of CsaV3_2G035120. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of C and G haplotype in highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_2G035120 in highly salt-sensitive accessions (R59, R132, CG13 and CG14) and highly salt-tolerant accessions (R45, R63, R78 and CG104) after salt treatment (0 h, 10 h, 48 h and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Ijms 23 08260 g005
Ijms 23 08260 g006 550
Figure 6. Identification of the candidate gene that was related to salt tolerance for the locus gST3.1. (a) Manhattan plot and LD heatmap of gST3.1. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of six functional annotated genes were predicted in the gST3.1 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) Gene structure of CsaV3_3G023710. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of GG and AA haplotype in the highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_3G023710 in the highly salt-sensitive accessions (R59, R132, CG13, and CG14) and highly salt-tolerant accessions (R45, R63, R78, and CG104) after salt treatment (0 h, 10 h, 48 h, and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Figure 6. Identification of the candidate gene that was related to salt tolerance for the locus gST3.1. (a) Manhattan plot and LD heatmap of gST3.1. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of six functional annotated genes were predicted in the gST3.1 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) Gene structure of CsaV3_3G023710. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of GG and AA haplotype in the highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_3G023710 in the highly salt-sensitive accessions (R59, R132, CG13, and CG14) and highly salt-tolerant accessions (R45, R63, R78, and CG104) after salt treatment (0 h, 10 h, 48 h, and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Ijms 23 08260 g006
Ijms 23 08260 g007 550
Figure 7. Identification of the candidate gene that was related to salt tolerance for the locus gST4.2. (a) Manhattan plot and LD heatmap of gST4.2. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of 13 functional annotated genes were predicted in the gST4.2 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) Gene structure of CsaV3_4G033150. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of G and C haplotypes in the highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_4G033150 in the highly salt-sensitive accessions (R59, R132, CG13, and CG14) and highly salt-tolerant accessions (R45, R63, R78, and CG104) after salt treatment (0 h, 10 h, 48 h, and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Figure 7. Identification of the candidate gene that was related to salt tolerance for the locus gST4.2. (a) Manhattan plot and LD heatmap of gST4.2. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of 13 functional annotated genes were predicted in the gST4.2 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) Gene structure of CsaV3_4G033150. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of G and C haplotypes in the highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_4G033150 in the highly salt-sensitive accessions (R59, R132, CG13, and CG14) and highly salt-tolerant accessions (R45, R63, R78, and CG104) after salt treatment (0 h, 10 h, 48 h, and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Ijms 23 08260 g007
Ijms 23 08260 g008 550
Figure 8. Identification of the candidate gene that was related to salt tolerance for the locus gST5.1. (a) Manhattan plot and LD heatmap of gST5.1. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of eight functional annotated genes were predicted in the gST5.1 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) The gene structure of CsaV3_5G023530. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of C and G haplotypes in the highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_5G023530 in the highly salt-sensitive accessions (R59, R132, CG13, and CG14) and highly salt-tolerant accessions (R45, R63, R78, and CG104) after salt treatment (0 h, 10 h, 48 h, and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Figure 8. Identification of the candidate gene that was related to salt tolerance for the locus gST5.1. (a) Manhattan plot and LD heatmap of gST5.1. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of eight functional annotated genes were predicted in the gST5.1 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) The gene structure of CsaV3_5G023530. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of C and G haplotypes in the highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_5G023530 in the highly salt-sensitive accessions (R59, R132, CG13, and CG14) and highly salt-tolerant accessions (R45, R63, R78, and CG104) after salt treatment (0 h, 10 h, 48 h, and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Ijms 23 08260 g008
Ijms 23 08260 g009 550
Figure 9. Identification of the candidate gene that as related to salt tolerance for the locus gST6.1. (a) Manhattan plot and LD heatmap of gST5.1. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of 12 functional annotated genes were predicted in the gST6.1 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) Gene structure of CsaV3_6G009810. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of AGG and GTT haplotype in the highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_6G009810 in the highly salt-sensitive accessions (R59, R132, CG13, and CG14) and highly salt-tolerant accessions (R45, R63, R78, and CG104) after salt treatment (0 h, 10 h, 48 h, and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Figure 9. Identification of the candidate gene that as related to salt tolerance for the locus gST6.1. (a) Manhattan plot and LD heatmap of gST5.1. The black dashed line represents the significance threshold and the red dashed lines represent the candidate region (100 kb). (b) A total of 12 functional annotated genes were predicted in the gST6.1 candidate region. The blue arrow indicates that the functional annotation of the gene which is related to salt tolerance. (c) Gene structure of CsaV3_6G009810. Blue rectangles and white rectangles indicate exons and introns, respectively. (d) Orange and blue blocks represent the proportion of AGG and GTT haplotype in the highly salt-tolerant or highly salt-sensitive accessions, respectively. (e) qRT-PCR was used to detect the expression changes of CsaV3_6G009810 in the highly salt-sensitive accessions (R59, R132, CG13, and CG14) and highly salt-tolerant accessions (R45, R63, R78, and CG104) after salt treatment (0 h, 10 h, 48 h, and 96 h). Significant difference (* p < 0.05, ** p < 0.01).
Ijms 23 08260 g009
Table
Table 1. The salt injury indexes of the germplasm in two years.
Table 1. The salt injury indexes of the germplasm in two years.
TreatmentMax. ScoreMin. ScoreMeanSTDSkewnessKurtosisC.VNormal Distribution Test (p Value)
Salt_202063.003.0031.1812.89−0.05−0.2541.35%0.83
Salt_202169.306.9528.8013.780.770.2247.84%0.06
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, D.; Dong, S.; Miao, H.; Liu, X.; Li, C.; Han, J.; Zhang, S.; Gu, X. A Large-Scale Genomic Association Analysis Identifies the Candidate Genes Regulating Salt Tolerance in Cucumber (Cucumis sativus L.) Seedlings. Int. J. Mol. Sci. 2022, 23, 8260. https://doi.org/10.3390/ijms23158260

AMA Style

Liu D, Dong S, Miao H, Liu X, Li C, Han J, Zhang S, Gu X. A Large-Scale Genomic Association Analysis Identifies the Candidate Genes Regulating Salt Tolerance in Cucumber (Cucumis sativus L.) Seedlings. International Journal of Molecular Sciences. 2022; 23(15):8260. https://doi.org/10.3390/ijms23158260

Chicago/Turabian Style

Liu, Dongrang, Shaoyun Dong, Han Miao, Xiaoping Liu, Caixia Li, Jianan Han, Shengping Zhang, and Xingfang Gu. 2022. "A Large-Scale Genomic Association Analysis Identifies the Candidate Genes Regulating Salt Tolerance in Cucumber (Cucumis sativus L.) Seedlings" International Journal of Molecular Sciences 23, no. 15: 8260. https://doi.org/10.3390/ijms23158260

APA Style

Liu, D., Dong, S., Miao, H., Liu, X., Li, C., Han, J., Zhang, S., & Gu, X. (2022). A Large-Scale Genomic Association Analysis Identifies the Candidate Genes Regulating Salt Tolerance in Cucumber (Cucumis sativus L.) Seedlings. International Journal of Molecular Sciences, 23(15), 8260. https://doi.org/10.3390/ijms23158260

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