A correction has been applied to this article in:
Corrigendum: Integrated GWAS and transcriptomic analysis reveal the candidate salt-responding genes regulating Na+/K+ balance in barley (Hordeum vulgare L.)
Tingting Xu
Shan Meng
Xiaopin Zhu
Jiachun Di
Yin Zhu
Xin Yang
Wei Yan- Institute of Germplasm Resources and Biotechnology, Jiangsu Key Laboratory for Agrobiolog, Jiangsu Provincial Platform for Conservation and Utilization of Agricultural Germplasm, Jiangsu Academy of Agricultural Sciences, Nanjing, China
Salt stress is one of the main abiotic stresses affecting crop yield and quality. Barley has strong salt tolerance, however, the underlying genetic basis is not fully clear, especially in the seedling stage. This study examined the ionic changes in barley core germplasms under the control and salt conditions. Genome-wide association study (GWAS) analysis revealed 54 significant SNPs from a pool of 25,342 SNPs distributed in 7 chromosomes (Chr) of the Illumina Barley 50K SNP array. These SNPs are associated with ion homeostasis traits, sodium (Na+) and potassium (K+) content, and Na+/K+ ratio representing five genomic regions on Chr 2, 4, 5, 6, and 7 in the leaves of worldwide barley accessions. And there are 3 SNP peaks located on the Chr 4, 6, and 7, which could be the “hot spots” regions for mining and identifying candidate genes for salt tolerance. Furthermore, 616 unique candidate genes were screened surrounding the significant SNPs, which are associated with transport proteins, protein kinases, binding proteins, and other proteins of unknown function. Meanwhile, transcriptomic analysis (RNA-Seq) was carried out to compare the salt-tolerant (CM72) and salt-sensitive (Gairdner) genotypes subjected to salt stress. And there was a greater accumulation of differentially expressed genes(DEGs) in Gairdner compared to CM72, mainly enriched in metabolic pathway, biosynthesis of secondary metabolites, photosynthesis, signal transduction,emphasizing the different transcriptional response in both genotypes following salt exposure. Combined GWAS and RNA-Seq analysis revealed 5 promising salt-responding genes (PGK2, BASS3, SINAT2, AQP, and SYT3) from the hot spot regions, which were verified between the salt-tolerant and salt-sensitive varieties by qRT-PCR. In all, these results provide candidate SNPs and genes responsible for salinity responding in barley, and a new idea for studying such genetic basis in similar crops.
1 Introduction
With increasing soil salinization, the area of arable land is decreasing worldwide threatening crop production and food safety (Han et al., 2015; Ramamoorthy et al., 2017). Though the use of modern farming and irrigation techniques can alleviate soil salinization to some extent, these methods require high input and have a slow response. Therefore, screening or cultivating new varieties of salt-tolerant crops can fundamentally solve the problem of saline-alkali land improving crop yield and quality.
Excessive accumulation of Na+ in plant tissues, on the one hand, causes osmotic stress, making it more difficult for the root system to absorb water, resulting in accelerated leaf senescence and death. On the other hand, it leads to a serious imbalance of intracellular ions, which affects the absorption of nutrients such as K, N, P, Ca, Mg, and Fe, resulting in plant physiology and metabolic disorders (Ruiz et al., 1997). Especially, since Na+ and K+ have similar hydration radii, non-selective cation channels or transporters cannot distinguish them, resulting in a large amount of Na+ competitively inhibiting K+ uptake (Assaha et al., 2017), disturbing the intracellular Na+/K+ ratio, leading to programmed cell death (PCD) (Shabala et al., 2007), indicating that the main toxic ions in saline-alkali soil cause osmotic stress and ion toxicity in plant cells, trigger oxidative stress, and affect mineral nutrient metabolism, plant seed germination, crop yield, and quality (Munns and Tester, 2008; Janz and Polle, 2012; Zhu et al., 2017). To counter salt stress, plants have evolved several physiological and biochemical strategies (Slama et al., 2015), such as the regionalization of Na+ in vacuoles by Na+/H+ antiporter (NHXs) (Jiang et al., 2010; Zhang et al., 2018; Sun et al., 2019) and efflux of Na+ from the cytosol to the apoplast via the Salt Overly Sensitive (SOS) signal pathway (Yang et al., 2009; Nutan et al., 2018; Zhang et al., 2018). High-affinity potassium transporter (HKT), which regulates the Na+/K+ transport, reduces Na+ accumulation in leaves by regulating phloem loading of Na+ (Hamamoto et al., 2015; Suzuki et al., 2016; Kobayashi et al., 2017; Huang et al., 2020). In addition, plants can reduce salt stress-induced osmotic stress by accumulating osmotic protective substances such as proline, betaine, late embryogenesis-enriched protein (LEA), and carbohydrates (Cuin and Shabala, 2005; Cuin and Shabala, 2007; Jia et al., 2020). In addition, phytohormones (Verma et al., 2016), NO (Nitric oxide) (Yemets et al., 2011), and polyamines (Shi and Chan, 2014) play important roles in maintaining ion homeostasis and signaling in plant salt tolerance networks.
Barley (Hordeum vulgare L.) is an ancient cereal crop considered to be the world’s most tolerant crop to salt, drought, low temperature, and cold stresses, which can be used as a model crop for studying the molecular regulation mechanism of salt tolerance (Nevo and Chen, 2010). At present, a lot of quantitative trait loci (QTLs) of physiological and biochemical indicators under salt stress have been identified using bi-parental populations (Mano and Takeda, 1997; Huang et al., 2009; Shavrukov et al., 2010; Xue et al., 2017). However, QTL mapping is a time-consuming and laborious process, resulting in the limited elucidation of the genetic mechanisms of plant salt tolerance.
In contrast, genome-wide association studies (GWAS) can find SNPs positioned across the genome as molecular markers and, therefore, is an efficient and convenient tool to dissect the genetic basis of complex traits in the plant, such as Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), maize (Zea mays), wheat (Triticum aestivum), and cotton (Gossypium hirsutum L.) (Atwell et al., 2010; Chen et al., 2015; Fang et al., 2017; Sun Z. et al., 2018). In barley, the major salt tolerance controlling QTL is on Chr 6, and a strong QTL for ion content is on Chr 4, which were identified using 192 spring barley genotypes with 1,000 SNPs (Long et al., 2013). Sbei et al. (2014) identified 7 significant QTLs related to salt tolerance using 296 Asian barley and 384 SNPs. Recently, using the Illumina 9K SNP array (Comadran et al., 2012), HKT1;5 was mapped on the distal part of Chr 4, which negatively regulates the salt tolerance in barley (Hazzouri et al., 2018; van Bezouw et al., 2019; Huang et al., 2020). Mwando et al. (2020) identified 19 loci containing 52 significant salt-tolerance-associated markers for salinity tolerance during germination. In recent years, with the publication of the whole genome and pan-genome of barley (Mayer et al., 2012; Mascher et al., 2017; Jayakodi et al., 2020), the first commercial high-density Barley SNP50K array was developed recently (Bayer et al., 2017), which can be used to accelerate the mining of salt tolerance loci or genes that were used in crop breeding for breeders.
Accordingly, this study identified the natural allelic variation in the barley genome and candidate genes that are significantly associated with salt tolerance. Here, we performed a GWAS using a natural population of 288 barley core collections and the Illumina Barley 50K SNP array based on ion traits for salt tolerance. In addition, we performed transcriptome analysis and qRT-PCR to verify the candidate genes identified by GWAS. In total, we identified 54 SNPs and 5 putative candidate genes associated with Na+/K+ balance in barley leaves. The results might provide useful candidate SNPs and genes for high-quality barley salt-tolerance breeding and for understanding the mechanisms for salt tolerance in barley.
2 Materials and methods
2.1 Barley germplasm
To evaluate salinity tolerance, a natural population consisting of a global barley core collection of 288 barley accessions (149 two-rowed and 139 six-rowed) collected from different countries or regions was used in this study. Among the accessions, 129 were from different provinces in China, 156 were from Australia, Germany, Japan, and other countries, and 3 were from unknown geographical sources (Supplementary Table 1) stored in the Jiangsu Provincial Platform for Conservation and Utilization of Agricultural Germplasm.
2.2 Seedlings hydroponic culture
The experiment was carried out in the greenhouse at Jiangsu Academy of Agricultural Science, Nanjing, China. Seeds with full and non-destructive granules were selected, and germinated in a plug with one seed per hole at 22/18 °C day/night conditions. Later, the seedlings (at the 1.5-2.0 leaf stage; 7-10 days old) were transferred to a turnover box containing modified Hoagland’s nutrient solution as described by Qiu et al. (2011). The box was covered with a plastic lid with evenly spaced holes; each hole contained 2 seedlings bundled with foam. After 10 days of transplanting, NaCl was added to the nutrient solution by 100 mM NaCl per day to reach a final concentration of 300 mM NaCl. Finally, two treatments were produced: (1) control treatment: 0 mM NaCl. (2) salt treatment: 300 mM NaCl. The nutrient solution was continuously oxygenated and changed weekly. After 3 weeks of salt stress and control treatment, the barley samples were taken for measurement. The whole experiment used a completely random design (CRD) under factorial arrangement with 3 replicates and four plants per replication.
2.3 Evaluation of salt tolerance
To estimate salt tolerance in barley, data of seedlings phenotypes were collected and statistically analyzed after 3 weeks of different treatment conditions. After 3 weeks of salt treatment, the sampled seedlings were rinsed with deionized water and dried with absorbent papers. Immediately after sample collection, plant height, fresh and dry weights, and other phenotypes were measured. Na+ and K+ contents in dried leaves were determined using an ICP-OES/MS spectrometer as described by Qiu et al. (2011). Salt tolerance (ST) was estimated as treatment/control×100% based on Na+ content, K+content, and Na+/K+ ratio. Statistical analysis, Pearson correlation between traits, ANOVA, and principal component analysis were performed by SPSS 20 software.
2.4 Genotyping and screening of SNP markers
The genomic DNA was extracted from young leaf tissues for genotyping using the CTAB method (Attitalla, 2011). The barley germplasms were genotyped using the Illumina BarleySNP50K array containing 44,040 SNPs (Bayer et al., 2017). The SNP data was filtered and processed by plink 1.90, with the criteria of calling rate < 0.9 and MAF < 0.05 (Chang et al., 2015).
2.5 Population structure analysis
The genetic structure of the population was estimated using the STRUCTURE 2.3.4 software; the K value was set from 2 to 10 and 15 iterations were used in an admixture model (Evanno et al., 2005). The obtained results were uploaded to the STRUCTURE Harvester website (http://taylor0.biology.ucla.Edu/structureHarvester/), when △K achieved maximum value, the optimal K was the number of clusters. Principal component analysis (PCA) was used to assess the population structure by GAPIT software (Tang et al., 2016). A neighbor-joining (NJ) phylogenetic tree was constructed by TASSEL software using individual sequences.
2.6 Genome-wide association analysis
The LD parameter (r2) between pairwise SNPs (MAF > 0.05) was estimated using the PopLDdecay software (Zhang et al., 2019). The association analysis was performed using the genome association and prediction integrated tool package (GAPIT). The significantly associated SNPs for the target traits were identified by the criteria -log10 (P) ≥-log10 (1/n), where n denotes the total filtered SNPs.
2.7 Identification of candidate genes
Candidate genes located within and/or adjacent to the associated SNPs were identified through the Barleymap (http://floresta.eead.csic.es/barleymap/) and Gramene (http://ensembl.gramene.org/Hordeum_vulgare/Info/Index) using the Morex genome. GO (Gene Ontology) and KEGG PATHWAY enrichment analyses were performed through the GENE ONTOLOGY (http://geneontology.org/).
2.8 Transcriptome sequencing
The seeds of salt-tolerant barley CM72 (C) and salt-sensitive barley Gairdner (G) were grown in a greenhouse, in Nanjing, China. Seedlings were collected after 0h, 3h, 12h, and 48 h with three biological replicates under 300 mM salt concentration and frozen in liquid nitrogen. Total RNA was extracted from the leaf samples using the EASY spin Plant RNA kit (Aidlab, Beijing, China), and sequencing libraries were generated using NEBNextR UltraTM Directional RNA Library Prep Kit for IlluminaR (NEB, USA) following manufacturer’s recommendations. The libraries were sequenced on an Illumina Hiseq Xten platform and paired-end reads were generated. Raw data (raw reads) of fastq format were firstly processed through in-house perl scripts by removing reads containing adapter, reads containing ploy-N and low quality reads from raw data. Then the clean reads were mapped to the reference genome sequence (https://webblast.ipk-gatersleben.de/barley_ibsc/downloads/) by HISAT2 tools software. Only reads with a perfect match or one mismatch were further analyzed and annotated based on the reference genome. All the downstream analyses were based on clean data with high quality. The original datas of the transcription sequencing were submitted to SRA database at NCBI (accession number: PRJNA866193 and PRJNA892763).
Gene expression levels were estimated by fragments per kilobase of transcript per million fragments mapped (FPKM). The false discovery rate (FDR) < 0.05 and |log2 (foldchange (FC))| ≥1 were set as the threshold by DESeq2 softwore (Tu et al., 2021) for significantly differential expression to identify differentially expressed gene (DEG) in response to salt stress. GO enrichment analysis of the DEGs was implemented by the clusterProfiler R package. We used KOBAS (Mao et al., 2005) software to test the statistical enrichment of differential expression genes in KEGG pathways and the clusterProfiler R packages to find KEGG pathway that are significantly enriched compared to the entire genome background.
2.9 qRT-PCR of the putative candidate genes
Six salt-tolerant and six salt-sensitive varieties were selected and grown in the greenhouse for 48 h under 300 mM salt stress, among them, CM72 and Gairdner were used as reference varieties (Supplementary Table 8). Three biological replicates of seedlings were collected and frozen in liquid nitrogen for RNA extraction. Next, total cDNA was synthesized using the PrimeScriptTM RT reagent Kit (TaKaRa).
qRT-PCRs were performed in a 20 ml reaction system for 1 min at 95 °C, followed by 40 cycles at 95 °C for 15 s, 60 °C for 15 s, and 72 °C for 45 s: the reaction mixture included 2.0 μL of cDNA, 0.6 μL of primer, 10 μL of SYBR Premix DimerEraser (TaKaRa) and 6.8 μL ddH2O. The used primers are listed in Supplementary Table 9. HvUBQ was used as the reference gene. All reactions were performed using the ABI7500 Real-Time PCR System with three independent biological replicates.
3 Results
3.1 Phenotypic evaluation of salt-tolerance traits of barley at the seedling stage under contrasting growth conditions
To find the new salt-tolerant variety of barley, a set of 288 barley germplasms from all over the world were examined for salt tolerance at the seedling stage. After salt treatment (300 mM NaCl) for 3 weeks, the barley genotypes showed typical symptoms of salt damage, such as stunting, leaf wilting, and yellowing (Supplementary Figure 1; Figures 1A, B).
FIGURE 1
Figure 1 Phenotypic analysis of barley accessions under control (CK: 0 mM) and salt treatment (300 mM NaCl) conditions after 3 weeks. Seedlings growth performance of (A) salt-tolerant and (B) salt-sensitive varieties of barley. Bars, 10 cm. (C) Raincloud plots for Na+ content, K+ content, and Na+/K+ ratio traits under contrasting growth conditions. **p < 0.01 (D) Frequency distribution of Na+, K+content, and Na+/K+ ratio under contrasting growth conditions and relative value of Na+, K+, Na+/K+ ratio. CK: 0 mM NaCl, Salt: 300 mM NaCl.
At the seedling stage, salt stress can affect the content of Na+ and K+ disturbing the ion balance, which affects the normal growth and development of plants. Compared with the control conditions, salt stress significantly increased the Na+ content (Salt-Na+) and Na+/K+ ratio (Salt-Na+/K+) from 5.46 to 48.66 mg/g (7.92-fold) and 0.09 to 2.25 (23.52-fold), respectively. Meanwhile, K+ content (Salt-K+) decreased significantly from 61.36 to 22.75 mg/g showing about a 62.92% reduction (Figure 1C and Table 1). Moreover, the relative value of Na+ (R-Na+), K+ (R-K+), and Na+/K+ (R-Na+/K+), which can be used to evaluate salt tolerance, ranged from 2.95 to 33.13, 0.21 to 0.66, and 6.42 to 115.75, respectively (Table 1), with a continuous distribution (Figure 1D), suggesting that the ion balance were controlled by multi-QTLs. The variance analysis based on ion-related traits showed significant differences among genotypes (Supplementary Table 2). Correlation analysis showed that Na+ and Na+/K+ were significantly negatively correlated with morphological traits (except YLN) and a positive correlation was also found between K+ and morphological traits (except YLN) (Supplementary Table 3). These results showed that salt stress conditions significantly increased the content of Na+ and reduced the content of K+, which disturbed the ion balance in barley and impaired the growth and development of plants.
TABLE 1
Table 1 Phenotypic variation for the salt-tolerance traits under salt treatment in the panel of barley accessions.
We further identified salt-tolerant (Figure 1B) and salt-sensitive barley germplasm (Figure 1B) based on the value of Salt-Na+, Salt-Na+/K+, R-Na+, and R-Na+/K+ as indicators of salt tolerance, which was obtained by PCA analysis (Table 2), indicating that these could be valuable germplasm for the development of salt-tolerant cultivars.
TABLE 2
Table 2 Principal component analysis of ionic traits.
3.2 Genotyping and screening of SNPs
The barley accessions were genotyped using the Barley 50K Illumina SNP array. A set of 25,342 mapped and high-quality SNPs were obtained using the criteria of MAF >0.05 and call rate >0.9 for genetic variation and GWAS analysis. These SNPs are evenly distributed on 7 chromosomes of barley (Figure 2A); Chr 5 has the maximum SNPs (4,707) while Chr 1 has the minimum (2,780). The polymorphism information content (PIC) values ranged from 0.354 to 0.369 among the 7 chromosomes (Table 3). At the same time, the T/C allele SNPs are the maximum (8,422), and the T/A allele SNPs are the minimum (1,656). This is consistent with the fact that each chromosome has the maximum number of T/C, while T/A is the least (Supplementary Table 4).
FIGURE 2
Figure 2 Population structure analysis of barley accessions. (A) Genome-wide SNP density in the entire association mapping panel. (B) Geographical distribution of barley accessions. (C) Population structure analysis revealed that barley accessions belong to four subgroups (K = 4). (D) Principal component analysis showing the four subpopulations. (E) A neighboring joining tree shows the four different subpopulations. (F) Genome-wide average LD decay distance in the barley genome.
TABLE 3
Table 3 The summary of the number of polymorphic SNPs mapped in barley genome.
3.3 Population structure and linkage disequilibrium analysis
We investigated the population structure of the barley germplasms based on 25,342 SNPs by STRUCTURE 2.3.4 software. At K = 4, the barley varieties were divided into 4 sub-populations (Figures 2B, C; Supplementary Table 1). Subgroup1 contained 74 barley accessions, about 70% of which were from abroad. Subgroup2 contained 77 barley accessions, of which, most were from Qinghai and Tibet. Subgroup3 contained 59 barley accessions mainly from China (44) and Japan (11), including 15 from Jiangsu. 78 barley accessions of subgroup4 were from Australia, North America, and Europe. Furthermore, the six-rowed barley belonged to subgroup2 and 4, while the two-rowed barley belonged to subgroup1 and 3. Meanwhile, PCA analysis (Figure 2D) and NJ tree also showed similar clustering results (Figure 2E).
We then analyzed LD in barley germplasm using filtered SNPs, and LD decay rates provided a moderate resolution for the identification of candidate gene regions (Bulik-Sullivan et al., 2015; Wijmenga, 2021). In this study, when the r2 value drops to half of the maximum value, the LD decay distance was approximately 5Kb, showing that the natural population had high genetic diversity. However, considering the large size of the barley genome, and the previously reported LD decay distance of barley was 2-8M (Long et al., 2013), we fixed 500 Kb as an interval distance for screening candidate genes around significantly associated SNPs (Figure 2F).
3.4 Identification of significant SNP loci by GWAS
To explore the genetic basis of salt tolerance in barley, 288 barley accessions, and 25,342 SNP markers were used for GWAS. In total, 54 significantly associated SNPs were found (criteria -Log10(P) ≥ 4.40) located mainly on Chr 2, 4, 5, 6, and 7. Some of these were detected repeatedly, indicating their true association with traits (Figure 3 and Table 4).
FIGURE 3
Figure 3 GWAS analysis for the ion-related traits under different conditions. (A) Manhattan plot for CK-Na+, CK-K+, and CK-Na+/K+, respectively. (B) Manhattan plot for Salt-Na+, Salt-K+, and Salt-Na+/K+, respectively. (C) Manhattan plot for R-Na+, R-K+, and R-Na+/K+, respectively.
TABLE 4
Table 4 The summary of SNPs significantly associated with ion-related traits.
Under control treatment, we detected 8 significant SNPs that correlated with ion-related traits (Figure 3A). 2 SNPs for K+ content (CK-K+) were detected on Chr 2 located at 6.193 and 647.542 Mb, respectively. There were 6 SNPs associated with Na+ content (CK-Na+), including 2 on Chr 5 and 4 on Chr 6, respectively. Under salt treatment, 20 significant SNPs were identified including 7 on Chr 4 and 9 on Chr 7 for Salt-Na+, which were nearby on their respective chromosomes (Figure 3B). A single SNP explained 14.85 to 17.81% of phenotypic variation. The strongest associated SNP, JHI-Hv50k-2016-272653 (-Log10 (P)=6.55) was found at 639.94 Mb on Chr 4 (Table 4). For Salt-Na+/K+, 30 significant SNPs was found on Chr 4 and 7 (Figure 3B), and a single SNP explained 17.85 to 27.73% of the phenotypic variation. Among these SNPs, 20 have been previously detected for Salt-Na+, including the strongest associated SNP JHI-Hv50k-2016-272653 (-Log10 (P)=7.35).
For the other two traits, R-Na+ and R-Na+/K+, 19 and 21 loci were detected, respectively. 10 SNPs for R-Na+ were found on Chr 6 (Figure 3C). 21 SNPs for R-Na+/K+ were found on Chr 4. Of these, JHI-Hv50k-2016-272613 and BOPA2_12_31219 were the most strongly associated SNP (-Log10 (P) = 6.77) located at 639.75 Mb. However, we did not find any SNP sites significantly associated with Na+/K+ ratio under control treatment (CK-Na+/K+), Salt-K+, and R-K+ (Figure 3). Overall, the 3 SNP peaks located on Chr 4, 6, and 7 were detected repeatedly, indicating that these regions may be the genomic hot spot related to salt tolerance in barley (Table 4).
3.5 Analysis of major SNP loci in SNP peak relevant to the ion-related traits
Next, we analyzed major SNP loci that are relevant to the Salt-Na+, Salt-Na+/K+, R-Na+, and R-Na+/K+. We found 25 significant SNPs within a significant peak located at the end of chromosome 4 (Table 4 and Figure 4A), which was detected in the previous studies (Hazzouri et al., 2018). There were eight SNPs detected repeatedly to further investigate the allelic variation. Haplotype analysis showed a high level of LD (r2 = 0.77–1.0) between the eight associated SNPs (Figure 4B), which resulted in two haplotypes (Figure 4C). The average Salt-Na+ of Hap1 (45.58 mg/g) was less than that of Hap2 (Figure 4D). Similar to Salt-Na+ data, the average values of Salt-Na+/K+, R-Na+, and R-Na+/K+ for Hap1 were 1.99, 9.58, and 25.51, respectively, which were all lower than that of Hap2 (Figure 4D).
FIGURE 4
Figure 4 Analysis of the SNP peak and the candidate genes on chromosome 4. (A) Manhattan plots for Chr 4. The black line represents the significance threshold (P < 10-4.40) and a red line indicates the position of the strong SNP peak. (B) LD based on pairwise r2 values between the SNPs estimated on Chr 4. The black inverted triangles indicate 8 significantly associated SNPs that were repeatedly detected. The green five-pointed star indicates the strongest SNP with the highest threshold. (C) Haplotypes were found among the barley accessions using the 8 SNPs. (D) Phenotypic differences of Salt-Na+, Salt-Na+/K+, R-Na+, and R-Na+/K+ between the two haplotypes.
Ten significant SNPs with a high level of LD located between 432 and 435 Mb at the end of Chr 6 were analyzed for R-Na+ (Table 3; Figures 5A, B). Among the 288 accessions, there were two haplotypes with the ten SNPs (Figure 5C), and the average R-Na+ of Hap1 (9.32) was 26.21% lower than that of Hap2 (Figure 5D). There was an SNP peak located at 608.8 Mb on Chr 7 (Figure 6A) containing 8 SNPs that were detected repeatedly and were in proximity related to Salt-Na+ and Salt-Na+/K+ traits. Haplotype analysis showed two haplotypes with a high level of LD (Figures 6B, C). The average Salt-Na+ and Salt-Na+/K+ of Hap1 were 46.46 mg/g and 2.07, respectively, which was significantly lower than the other one (Figure 6D).
FIGURE 5
Figure 5 Analysis of the SNP peak on chromosome 6. (A) Manhattan plots for Chr 6. The black line represents the significance threshold (P < 10-4.40) and a red line indicates the position of the strongest SNP peak. (B) LD based on pairwise r2 values between the SNPs estimated on Chr 6. The black inverted triangles show 10 significantly associated SNPs. The green five-pointed star indicates the strongest SNP with the highest threshold. (C) Haplotypes were observed in the barley accessions using the 10 SNPs. (D) Phenotypic differences of R-Na+ between the two haplotypes.
FIGURE 6
Figure 6 Analysis of the SNP peak on chromosome 7. (A) Manhattan plots for Chr 7. The black line represents the significance threshold (P < 10-4.40) and a red line indicates the position of the strongest SNP peak. (B) LD based on pairwise r2 values between the SNPs estimated on Chr 7. The black inverted triangles show 8 significantly associated SNPs that were detected repeatedly. The green five-pointed star indicates the strongest SNP with the highest threshold. (C) Haplotypes were observed in the barley accessions using the 8 SNPs. (D) Phenotypic differences of Salt-Na+ and Salt-Na+/K+ between the two haplotypes.
For the genetic effect of each SNP in Hap1 of Chr4, for example, JHI-Hv50k-2016-272613 could be classified into three genotypes. And the TT genotype displayed a lower average Salt-Na+, Salt-Na+/K+, R-Na+, and R-Na+/K+, respectively, than the other nonincluded genotypes (Supplementary Figure 2). The other SNPs in Hap1 of Chr6 and Chr7 had the same effect on the phenotype as the locus JHI-Hv50k-2016-272613 with the lower Salt-Na+ and Salt-Na+/K+, respectively (Supplementary Figure 3) and R-Na+ (Supplementary Figure 4).
3.6 Exploration of candidate genes associated with SNPs
Based on the criteria of 500 Kb distance from significant SNPs, we found 616 candidate genes. More than half of these genes have no functional annotation. The rest of the annotated genes are involved in binding proteins, transport proteins, protein kinases, and other unknown function proteins (Supplementary Table 5). For example, HORVU4Hr1G087960 code for a high-affinity K+ transporter (HvHKT1;5), which is involved in the ionic stress signaling pathway (Qiu et al., 2011; Hazzouri et al., 2018; van Bezouw et al., 2019; Huang et al., 2020; Houston et al., 2020; Wege et al., 2021). HORVU5Hr1G124800 is annotated for serine/threonine-protein kinase and may be involved in MAPK cascades responding to salt or osmotic stress. Furthermore, GO enrichment and KEGG pathway analysis were carried out to explore the role of these genes.
GO analysis indicated that in terms of molecular function, the rRNA N-glycosylase activity (GO:0030598) and RNA glycosylase activity (GO:0030597) with the common four genes were the most significantly enriched categories (Supplementary Table 6). GO cellular component analysis suggested that the candidate genes are mainly located in the NADH dehydrogenase complex (GO:0030964), autophagosome membrane (GO:0000421), and mitochondrial membrane (GO:0031966). Meanwhile, 5 genes showed significant enrichment for biological processes such as negative regulation of translation (GO:0017148) and cellular amide metabolic process (GO:0034249). Likewise, KEGG pathway analysis revealed that only 5 genes were enriched in 9 pathways. HORVU4Hr1G087100 encodes branched-chain-amino-acid aminotransferase 6, which participates in the synthesis of alanine, valine, leucine, and isoleucine. The other pathways include cytoskeletal regulation by Rho GTPase (HORVU4Hr1G087740), Huntington’s disease (HORVU4Hr1G087740, HORVU4Hr1G088610, and HORVU6Hr1G016100) and so on (Supplementary Table 6).
3.7 RNA-Seq analysis in response to salt stress in barley cultivars
To understand the regulatory mechanisms of salt stress and explore promising candidate genes, we determined the transcriptome changes in the leaves of CM72 and Gairdner under salt stress using RNA-seq. There were 21.54 million raw reads on average and more than 91.98% of the reads of each sample (0, 3, 12, and 48 h) aligned with the specified reference genome and led to the identification of 90,591 genes (Supplementary Table 7). With the increase of salt treatment time, the number of DEGs also gradually increased (Supplementary Figure 5A). Venn diagrams showing the number of distinct and common salt-responsive genes (DEGs) in the two contrasting genotypes in different point (Supplementary Figure 5B). Principal Component Analysis (PCA) showed that the samples separated into four major clusters, indicating that the salt treatment time was likely the largest source of variation in the data (Supplementary Figure 5C). Therefore, we chose the data of controls (0 h) and 48 h salt-treated samples of both genotypes to further bioinformatic analyses.
After 48h of salt treatment, 17,437 genes were significantly differentially expressed under salt stress, and of these, 5368 genes were commonly deferentially expressed in both varieties, while 4,125 and 8,004 genes were specifically expressed in CM72 and Gairdner, respectively. Moreover, there were 90 up-regulated and 28 down-regulated genes in CM72 while which were down-regulated and up-regulated in Gairdner (Figures 7A-C).
FIGURE 7
Figure 7 Transcriptome analysis in CM72 and Gairdner under salt conditions. (A) Venn diagram showing the number of transcripts detected as up and down-regulated (P < 0.01) in CM72 and Gairdner under salt conditions. (B) Heat maps showing differentially expressed transcripts in CM72 under control and salt conditions. (C) Heat maps showing differentially expressed transcripts in Gairdner under control and salt conditions. (D) The Gene Ontology (GO) enrichment analysis in CM72 under salt conditions. (E) The Gene Ontology (GO) enrichment analysis in Gairdner under salt conditions. (F) The top 20 significantly KEGG enrichment pathways in CM72 under salt conditions. (G) The top 20 significantly KEGG enrichment pathways in Gairdner under salt conditions.
Gene ontology (GO) enrichment analysis is performed to determine the potential function of DEGs in response to salt stress. The results of GO enrichment analysis indicated that in terms of cellular component the enriched same GO terms were cell, cell part, membrane, organelle in both CM72 and Gairdner. GO molecular function analysis suggested that the DEGs mainly richer in binding, catalytic activity, structural molecule activity, transporter activity in both CM72 and Gairdner. Within the biological process, the mainly enriched GO terms were biological regulation, cellular process, metabolic process, response to stimulus) in both CM72 and Gairdner (Figures 7D, E).
KEGG pathway enrichment analysis showed that the DEGs were mainly enriched into the pathway of “ribosome”, “phenylpropanoid biosynthesis” and “carbon fixation in photosynthetic organisms” in both varieties. Otherwise, the specific term KEGG pathway were “glyoxylate and dicarboxylate metabolism” and “photosynthesis - antenna proteins” in CM72, and “cutin, suberine and wax biosynthesis” and “plant hormone signal transduction” in Gairdner, respectively (Figures 7F, G).
3.8 Putative candidate genes identified by GWAS and RNA-seq analysis
Further, we investigated the expression level of the 616 genes and found that most of these showed no expression at the four-time points (Supplementary Table 5). A total of 59, 59, and 17 candidate genes were within the LD decay region (hot spots) on Chr 4, 6, and 7, respectively (Figure 8; Supplementary Table 5). Based on GWAS and RNA-seq analysis, we selected 5 candidate genes that are associated with the regulation of Na+/K+ balance (Figure 9 and Table 5). Among them, HORVU4Hr1G088190 is located in the SNPs peak with three significant SNPs JHI-Hv50k-2016-272632, JHI-Hv50k-2016-272635, and JHI-Hv50k-2016-272653 are at 639.7-639.95 Mb (Figure 8A). HORVU6Hr1G016070 and HORVU6Hr1G016120 are in proximity of peak SNP JHI-Hv50k-2016-383164 at 35.79-36.06 Mb. HORVU6Hr1G064140 located near the peak SNP JHI-Hv50k-2016-405919 at 433.30-433.71 Mb (Figure 8B). These three genes were annotated as sodium bile acid symporter family protein, SINA (seven in absentia) E3 ubiquitin-protein ligase SINAT2, and aquaporin-like superfamily protein, respectively. HORVU7Hr1G101310, located within the SNPs peak JHI-Hv50k-2016-497584, JHI-Hv50k-2016-497604, and JHI-Hv50k-2016-497611 at 608.81 Mb (Figure 8C), encodes a calcium-dependent lipid-binding (CaLB domain) family protein. The other three selected genes, HORVU2Hr1G002690, HORVU5Hr1G124800, and HORVU5Hr1G124880, are located near the association peak SNP loci of other chromosomes for CK-Na+ and CK-K+. In addition, analysis of the eight candidate genes by qRT-PCR appeared to follow a similar pattern as RNA-Seq (Figure 9; Supplementary Figure 6A).
FIGURE 8
Figure 8 Expression pattern of candidate genes in the peak regions of Chr4, Chr6 and Chr7. (A) Heat map showing the expression levels of 59 candidate genes in Chr4 among two barley cultivars. (B) Heat map showing the expression levels of 59 candidate genes in Chr6 among two barley cultivars. (C) Heat map showing the expression levels of 17 candidate genes in Chr6 among two barley cultivars. The barley leaf samples were collected at 0, 3, 12, and 48 h under salt stress.
FIGURE 9
Figure 9 Expression pattern validation of five selected genes (HORVU4Hr1G088190 (PGK2), HORVU6Hr1G016070 (BASS3), HORVU6Hr1G016120 (SINAT2), HORVU6Hr1G064140 (AQP), and HORVU7Hr1G101310 (SYT3)), which were also detected by RNA-Seq, in CM72 and Gairdner leaves by qRT-PCR. The samples were collected at 0, 3, 12, and 48 h under salt stress.
TABLE 5
Table 5 The putative candidate genes in SNP peaks associated with salt-tolerance related traits identified by GWAS and RNA-seq.
The expression levels of the above 5 genes in the SNP peak region and other 3 genes were measured between the six salt-tolerant and six salt-sensitive varieties by qRT-PCR (Figure 10; Supplementary Figure 6B). HORVU4Hr1G088190, HORVU7Hr1G101310, and HORVU5Hr1G124880 were up-regulated in the salt-tolerant varieties, while HORVU6Hr1G016070, HORVU6Hr1G064140, and HORVU2Hr1G002690 were higher in the salt-sensitive varieties. This is consistent with the RNA-seq results showing that HORVU7Hr1G101310 and HORVU5Hr1G124880 were up-regulated and HORVU6Hr1G016070, HORVU6Hr1G064140 were down-regulated (Supplementary Table 5). However, the DEGs HORVU6Hr1G016120 and HORVU5Hr1G124800 showed no difference in their expression level between the salt-sensitive and salt-tolerant varieties. In all, our data suggest that the above genes could be promising candidates genes as their orthologs also play a role in salt tolerance in other species (Hu et al., 2012; Liu et al., 2013; Joshi et al., 2016; Lee et al., 2019, Krausko et al., 2022, Myo et al., 2020, Wang et al., 2020).
FIGURE 10
Figure 10 Relative expression of 5 candidate genes by RT-PCR in six salt-sensitive and six salt-tolerant varieties. HORVU4Hr1G088190 (PGK2). HORVU6Hr1G016070 (BASS3). HORVU6Hr1G016120 (SINAT2). HORVU6Hr1G064140 (AQP). HORVU7Hr1G101310 (SYT3). Six salt-sensitive varieties, S1, Shan 3; S2, AKHELOL 1; S3, Boyer; S4, Rainbow; S5, Birka; S6, Gairdner; six salt-tolerant varieties; T1, PALLIDUM 043; T2, KIMALUNG; T3, Zhu Damai 5; T4, GURZAN; T5, Kunlun 14; T6, CM72.
4 Discussion
Increasing soil salinization has become a serious environmental challenge affecting crop yield and quality. Therefore, screening salt-tolerant crop varieties, identifying salt-tolerant genes, and breeding salt-tolerant varieties can be the most economical and effective strategies to alleviate the problem of soil salinization. The seedling stage which is most sensitive to stress is very important in the barley growth period. Therefore, the salt tolerance of barley measured by morphological and ion traits at the seedling stage reflects the salt tolerance of the variety to a certain extent. In our study, by phenotyping the 288 barley accessions under salt treatment, relatively abundant variation and significant differences among genotypes were observed for the ion traits that were selected to evaluate the salt tolerance of barley seedlings (Supplementary Table 2). Here, we identified some salt-tolerant varieties of barley with lower Na+ accumulation, R-Na+, Salt-Na+/K+, and Salt-Na+/K+ in leaves compared with salt-sensitive varieties (Figure 1A; Supplementary Table 9). Notably, previous studies have shown that salt-tolerant lines had higher Na+ content in roots (Tu et al., 2021), indicating the difference in translocation of Na+ from root to shoot among the two genotypes. Interestingly, the salt-tolerant varieties contain excellent Hap1 SNPs, while salt-sensitive varieties contain Hap2 SNPs. These findings can be further used to study the mechanism of salt tolerance in barley and prepare parents for salt tolerance breeding.
Some studies used QTL mapping to identify candidate intervals and molecular markers of salt tolerance in barley and the population constructed by bio-parents populations. Many QTL loci were identified (Mano and Takeda, 1997; Huang et al., 2009; Shavrukov et al., 2010; Angessa et al., 2017; Liu et al., 2017; Xue et al., 2017), but specific candidate genes related to salt tolerance could not be identified that slow the progress of salt-tolerant breeding in barley. GWAS, a novel and effective method, has been widely used to determine the genetic basis of crop plants. In recent years, several efforts have been made to detect QTLs related to salt tolerance using GWAS with simple sequence repeat (SSR) or SNP markers in barley (Long et al., 2013; Sbei et al., 2014; Elakhdar et al., 2016; Fan et al., 2016; Tu et al., 2021). However, the number of SNP markers was small which limited the mapping resolution and genome coverage, and also the SNPs were not found relevant to salt-tolerant loci on other chromosomes.
Here, we conducted a GWAS of ion-related salt-tolerant traits at the seedling stage of barley accessions having rich genetic diversity and 25,342 high-quality SNPs from the BarleySNP50K array for the first time. In total, 54 significant SNPs (P < 10-4.40) were identified for ion traits, and there were three SNP peaks located on Chr 4, 6, and 7 (Figure 3 and Table 3). Furthermore, the SNP peak located on Chr 4 was consistent with the genetic interval noticed by Hazzouri et al. (2018) (Figure 4). Also, we found that the SNPs in the three SNP peaks can produce two haplotypes. Most of the barley germplasms belonged to the Hap1 group having lower Salt-Na+ content, R-Na+, Salt-Na+/K+ ratio, and/or R-Na+/K+ ratio (Figures 4C, D, 5C, D and 6C, D). This indicated that the Hap1 group has superior salt tolerance and may provide new insights into the genetic basis of salt tolerance in barley. In addition, the genetic effect of each SNP located in candidate genes, such as i.e.PGK2 with 3 significant SNPs (JHI-Hv50k-2016-272632, JHI-Hv50k-2016-272635, and JHI-Hv50k-2016-272653) (Supplementary Figure 2) on salt tolerance was manifested in that the genotypes in haplotype 1 reduced value of ion traits (Supplementary Figures 3, 4), and these SNPs could be used to develop functional molecular markers for further screening and identification of barley salt-tolerant germplasm from ion traits.
Meanwhile, we screened 616 candidate genes surrounding the significant SNPs that are involved in kinases, transcription factors, transporters, and cellular protein components (Supplementary Table 5). To further select possible salt tolerance genes, we performed transcriptome analysis in CM72 and Gairdner leaves. We found that the number of DEGs increased with salt treatment time, mainly enriched into the metabolic pathway, biosynthesis of secondary metabolites, photosynthesis, signal transduction and other pathways (Figures 7F, G), indicating that NaCl stress can promote the up-regulation or down-regulation of related genes in vivo, and adjust related pathways by changing the expression of genes to adapt to salt stress. Taken together, 8 promising salt-responding genes (PGK2, BASS3, SINAT2, AQP, SYT3, receptor kinase 2, GsSRK, and MATE) from the hot spot regions as protein kinases and transporters were identified by GWAS and RNA-Seq integrating analysis (Figures 9, 10; Supplementary Figure 5), which might play roles in controlling ion content under salt stress. In addition, a total of 119 SNPs were found in the 8 putative candidate genes for alignment with the Morex genome in CM72 and Gairdner by RNA-seq data (Supplementary Table 10), showed the allelic variation among contrasting genotypes. However, whether these allelic variants could explain the salt tolerance of barley needs further study.
In previous studies, many ion transporter genes have been reported, involving Na+ absorption in roots, root-to-shoot Na+ transportation, Na+ translocation to shoots, and the localization of Na+ in the vacuole regionalization (Bahieldin et al., 2015; Gharaghanipor et al., 2022; Ouertani et al., 2022). In the current study, we have identified sodium transporter (HKT1), which was a down-regulated gene by RNA-seq, consistent with the previous findings that HKT negatively regulates salt tolerance when HKT;5 was knockdown, Na+ translocation from roots to shoots decrease, and increases in K+/Na+ in barley (Huang et al., 2020). However, some other well-known Na+ transporters (i.e. NHX, and SOS family genes) were not identified as common candidates (Supplementary Table 5), but we found that some novel genes in barley salt tolerance have not been reported and remain to be further verified and discussed in future research.
More importantly, some novel genes, including ion transporters and signal regulatory genes, were widely identified in the leaves of barley collections in salt response (Figure 8). HORVU4Hr1G088190 encodes 2-phosphoglycerate kinase (PGK2), which is an important enzyme localized in cytoplasm as well as stroma of chloroplast in photosynthetic carbon metabolism, glycolysis, gluconeogenesis, and reductive pentose phosphate cycle (McMorrow and Bradbeer, 1990). It was demonstrated that ectopic expression of OsPGK2-P improved salt tolerance in tobacco by maintaining better ion homeostasis, and higher chlorophyll retention and proline accumulation (Kumari et al., 2009; Joshi et al., 2016). Earlier studies reported that the PGK2 is an isozyme of the glycolytic pathway that provides ATP in plants, also known as ATPase with ATP binding molecular function (GO:0005524) (Supplementary Table 6), including vacuolar H+-ATPase (V-ATPase) to transport H+ in the cytoplasm to the vacuole, and plasma membrane ATPase (plasma membrane H+-ATPase, PM-ATPase) to transport the H+ in the cytoplasm to the outside of the cell to maintain cellular pH homeostasis (Sawyer et al., 2008 Danshina et al., 2010; Krebs et al., 2010). PGK2 was supposed to regulate ion homeostasis and maintain the intracellular pH balance with a high expression level in salt-tolerant varieties, although which was not the DEG in RNA-Seq analysis (Table 5). Additionally, the three genes HORVU4Hr1G087760, HORVU4Hr1G087960, and HORVU4Hr1G088140 have been identified previously (Hazzouri et al., 2018),encoding bifunctional inhibitor/lipid-transfer protein/seed storage 2S albumin superfamily protein, sodium transporter HKT1 and expansin B2, respectively. In all, our results indicate that there could be a hot spot of salt tolerance candidate genes at the end of Chr 4 located at 639.94 Mb.
Additionally, the three genes HORVU6Hr1G016070, HORVU6Hr1G016120, and HORVU6Hr1G064140 are homologous genes of BASS3, E3 ubiquitin-protein ligase SINAT2, and aquaporin PIP1 in rice, respectively. BASS family proteins encode a class of sodium/solute symporters with Na+-coupled metabolite transport properties which is predicted a Na+-dependent transporter (Zhou et al., 2014), so it may have a potential role in the Na+ physiology in plant salt tolerance. In wheat, TaBASS2, a putative pyruvate transporter, acts as a positive regulator of salinity tolerance via modulation of ABI4 expression (Zhao et al., 2016). The over-expression of GhBASS5 in Arabidopsis impairs salt tolerance by increasing Na+ loading and accumulation (Myo et al., 2020). SINA E3 ubiquitin ligase is a RING (really interesting new gene) type E3 ubiquitin ligase, which is involved in plant growth and development, response to stress, interactions between plants and other organisms, and autophagy in plant cells. In rice, OsDIS1 (Oryza sativa drought-induced SINA protein 1) inhibited the drought resistance of rice by interacting with OsSKIPa (a positive regulator of drought and salt stress), affecting the expression of downstream stress response genes to regulate the drought resistance (Ning et al., 2011a; Ning et al., 2011b). In Arabidopsis, AtSINA2 can interact with CDKG1 (cyclin-dependent protein kinase G1), which phosphorylates AtSINA2 and affects its stability, thereby regulating the plant response to ABA and osmotic stress in Arabidopsis (Bao et al., 2014; Chen et al., 2018). Aquaporin (AQPs) plays a major role in maintaining the water status and transport across the plant which has water channel activity (GO:0015250) and involves ion transmembrane transport (GO:0034220) (Supplementary Table 6), that is crucial for plants to combat salt stress (Kumari and Bhatla, 2021). TaAQP8 confers salt stress tolerance by retaining a high K+/Na+ ratio and Ca2+ content and enhancing the antioxidant system (Hu et al., 2012). Aquaporin OsPIP1;1 functions as an active water channel that promotes rice salt resistance and seed germination (Liu et al., 2013). In tobacco, ectopic expression of SpAQP1 from halophyte Sesuvium portulacastrum in transgenic tobacco could increase salt tolerance by enhancing the anti-oxidative activity of plants (Chang et al., 2016).
On Chr 7, the gene HORVU7Hr1G101310 encodes a calcium-dependent lipid-binding (CaLB domain) family protein with molecular function of metal ion binding (GO:0046872) that homologous to Arabidopsis synaptotagmin-3. In Arabidopsis, the stress-induced SYT3 is an endoplasmic reticulum-plasma membrane (ER-PM) tether that also participates in maintaining PM integrity along with SYT1 (Ruiz-Lopez et al., 2021). Under salt stress, loss of AtSYT1 function resulted in a decrease in photosynthetic efficiency accompanied by a decrease in chlorophyll and carotenoids and an increase in flavonol content in leaves, aggravating the adverse effects of salt stress (Krausko et al., 2022). However, analysis of the above genes in barley salt tolerance has not been reported.
In addition, the gene HORVU2Hr1G002690 encodes receptor kinase 2 and plays a crucial role in plant responses to abscisic acid (ABA), salt, and drought stresses (Sun et al., 2013a; Shi et al., 2014). The gene HORVU5Hr1G124800 encodes a serine/threonine-protein kinase, also called G-type lectin S-receptor-like kinase (GsSRK) with ATP binding (GO:0005524), is involved in ion homeostasis, ROS scavenging, and osmotic regulation, and have a dual role in regulating both plant architecture and salt stress responses (Sun et al., 2013b; Sun M. Z. et al., 2018). The gene HORVU5Hr1G124880 encodes the MATE efflux family protein with antiporter activity (GO:0015297) which modulates ABA sensitivity and increases tolerance to drought, salt, and cold stress with high production of antioxidant enzymes and significantly reduced levels of oxidants in plants (Zhang et al., 2014; Lu et al., 2018; Lu et al., 2019). On the other hand, the two other genes, HORVU7Hr1G101590 and HORVU7Hr1G103510 encode Cytochrome P450 superfamily protein and transducin/WD40 repeat-like superfamily protein (WD40), respectively. Both of these gene products participate in response to salt stress (Mao et al., 2013, Sun Z et al., 2018, Sun et al., 2021; Xia et al., 2021; Xin et al., 2021). However, the detailed function of these genes concerning salt tolerance in barley remains to be investigated by gene-editing and gain-of-function-based genetic engineering.
5 Conclusions
In this study, we first genotyped 288 diverse accessions having rich genetic diversity using the BarleySNP50K array. In total, 54 SNPs were significantly associated with salt tolerance traits and 616 candidate genes were identified by GWAS. There were 3 SNP peaks located on Chr 4, 6, and 7 that are associated with two and more traits. Furthermore, using integrated GWAS, RNA-seq, and qPCR analysis, we screened 5 putative genes (PGK2, BASS3, SINAT2, AQP, and SYT3) that may be related to ion balance in leaves from phenotypic and haplotype screening. This study provides novel SNPs and candidate genes for the study of salt tolerance in barley. Our results provide a basis for breeding new salt-tolerant barley varieties with marker-assisted selection.
Data availability statement
The data presented in the study are deposited in the SRA database at NCBI repository, accession number PRJNA866193 and PRJNA892763 for RNA-seq data; and European Variation Archive repository for SNP data, accession number PRJEB59041.
Author contributions
TX and WY conceived and designed the idea. TX, SM, XZ, JD, YZ, and XY performed the experiments. TX and WY analyzed the data. TX wrote and revised the manuscript. SM, XZ, JD, YZ, XY, and WY critically revised the article. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (32001543), Jiangsu Agricultural Science and Technology Innovation Fund (CX(22)3139), the Jiangsu Agricultural Science and technology Innovation and Popularization Special project (2021-SJ-004), and the Jiangsu Science and Technology development subsidy special project (26012101).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1004477/full#supplementary-material
References
Angessa, T. T., Zhang, X. Q., Zhou, G., Broughton, S., Zhang, W., Li, C. (2017). Early growth stages salinity stress tolerance in CM72 x gairdner doubled haploid barley population. PloS One 12, e179715. doi: 10.1371/journal.pone.0179715
Assaha, D. V.M., Ueda, A., Saneoka, H., Al-Yahyai, R., Yaish, M. W. (2017). The role of na+ and k+ transporters in salt stress adaptation in glycophytes. Front. Physiol. 18. doi: 10.3389/fphys.2017.00509
Attitalla, I. H. (2011). Modified CTAB method for high quality genomic DNA extraction from medicinal plants. Pak. J. Biol. Sci. 14, 998–999. doi: 10.3923/pjbs.2011.998.999
Atwell, S., Huang, Y. S., Vilhjalmsson, B. J., Willems, G., Horton, M. Y. L., et al. (2010). Genome-wide association study of 107 phenotypes in a common set of arabidopsis thaliana inbred lines. Nature 465, 627–631. doi: 10.1038/nature08800
Bahieldin, A., Atef, A., Sabir, J. S., Gadalla, N. O., Edris, S., Alzohairy, A. M., et al. (2015). RNA-Seq analysis of the wild barley (H. spontaneum) leaf transcriptome under salt stress. Comptes Rendus. Biol. 338, 285–297. doi: 10.1016/j.crvi.2015.03.010
Bao, Y., Wang, C. T., Jiang, C. M., Pan, J., Zhang, G. B., Liu, H., et al. (2014). The tumor necrosis factor receptor-associated factor (TRAF)-like family protein SEVEN IN ABSENTIA 2 (SINA2) promotes drought tolerance in an ABA-dependent manner in arabidopsis. New Phytol. 202, 174–187. doi: 10.1111/nph.12644
Bayer, M. M., Rapazote-Flores, P., Ganal, M., Hedley, P. E., Macaulay, M., Plieske, J., et al. (2017). Development and evaluation of a barley 50k iSelect SNP array. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.01792
Bulik-Sullivan, B. K., Loh, P. R., Finucane, H. K., Ripke, S., Yang, J., Patterson, N., et al. (2015). LD score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 47, 291–295. doi: 10.1038/ng.3211
Chang, C. C., Chow, C. C., Tellier, L. C., Vattikuti, S., Purcell, S. M., Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4, 7. doi: 10.1186/s13742-015-0047-8
Chang, W. J., Liu, X. W., Zhu, J. H., Fan, W., Zhang, Z. L. (2016). An aquaporin gene from halophyte sesuvium portulacastrum, SpAQP1, increases salt tolerance in transgenic tobacco. Plant Cell Rep. 35, 385–395. doi: 10.1007/s00299-015-1891-9
Chen, Y., Fokar, M., Kang, M. Y., Chen, N. C., Allen, D. R., Chen, Y. F. (2018). Phosphorylation of arabidopsis SINA2 by CDKG1 affects its ubiquitin ligase activity. BMC Plant Biol. 18, 147. doi: 10.1186/s12870-018-1364-8
Chen, G., Wang, X., Hao, J., Yan, J., Ding, J. (2015). Genome-wide association implicates candidate genes conferring resistance to maize rough dwarf disease in maize. PloS One 10, e142001. doi: 10.1371/journal.pone.0142001
Comadran, J., Kilian, B., Russell, J., Ramsay, L., Stein, N., Ganal, M., et al. (2012). Natural variation in a homolog of antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley. Nat. Genet. 44, 1388–1392. doi: 10.1038/ng.2447
Cuin, T. A., Shabala, S. (2005). Exogenously supplied compatible solutes rapidly ameliorate NaCl-induced potassium efflux from barley roots. Plant Cell Physiol. 46, 1924–1933. doi: 10.1093/pcp/pci205
Cuin, T. A., Shabala, S. (2007). Amino acids regulate salinity-induced potassium efflux in barley root epidermis. Planta 225, 753–761. doi: 10.1007/s00425-006-0386-x
Danshina, P. V., Geyer, C. B., Dai, Q., Goulding, E. H., Willis, W. D., Kitto, G. B., et al. (2010). Phosphoglycerate kinase 2 (PGK2) is essential for sperm function and male fertility in mice. Biol. Reprod. 82, 136–145. doi: 10.1095/biolreprod.109.079699
Elakhdar, A., EL-Sattar, M. A., Amer, K., Rady, A., Kumamaru, T. (2016). Population structure and marker–trait association of salt tolerance in barley (Hordeum vulgare l.). Comptes Rendus Biologies 339, 454–461. doi: 10.1016/j.crvi.2016.06.006
Evanno, G., Regnaut, S., Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x
Fang, L., Wang, Q., Hu, Y., Jia, Y., Chen, J., Liu, B., et al. (2017). Genomic analyses in cotton identify signatures of selection and loci associated with fiber quality and yield traits. Nat. Genet. 49, 1089–1098. doi: 10.1038/ng.3887
Fan, Y., Zhou, G., Shabala, S., Chen, Z. H., Cai, S., Li, C., et al. (2016). Genome-wide association study reveals a new QTL for salinity tolerance in barley (Hordeum vulgare l.). Front. Plant Sci. 7. doi: 10.3389/fpls.2016.00946
Gharaghanipor, N., Arzani, A., Rahimmalek, M., Ravash, R. (2022). Physiological and transcriptome indicators of salt tolerance in wild and cultivated barley. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.819282
Hamamoto, S., Horie, T., Hauser, F., Deinlein, U., Schroeder, J. I., Uozumi, N. (2015). HKT transporters mediate salt stress resistance in plants: from structure and function to the field. Curr. Opin. Biotechnol. 32, 113–120. doi: 10.1016/j.copbio.2014.11.025
Han, J., Shi, J., Zeng, L., Xu, J., Wu, L. (2015). Effects of nitrogen fertilization on the acidity and salinity of greenhouse soils. Environ. Sci. pollut. Res. Int. 22, 2976–2986. doi: 10.1007/s11356-014-3542-z
Hazzouri, K. M., Khraiwesh, B., Amiri, K., Pauli, D., Blake, T., Shahid, M., et al. (2018). Mapping of HKT1;5 gene in barley using GWAS approach and its implication in salt tolerance mechanism. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00156
Houston, K., Qiu, J., Wege, S., Hrmova, M., Oakey, H., Qu, Y., et al. (2020). Barley sodium content is regulated by natural variants of the na+ transporter HvHKT1;5. Commun. Biol. 3, 258. doi: 10.1038/s42003-020-0990-5
Huang, L., Kuang, L., Wu, L., Shen, Q., Han, Y., Jiang, L., et al. (2020). The HKT transporter HvHKT1;5 negatively regulates salt tolerance. Plant Physiol. 182, 584–596. doi: 10.1104/pp.19.00882
Huang, D. X. Y., Li, K. W. S. W., Lance, M. C. G. Z. (2009). Identification of QTLs associated with salinity tolerance at late growth stage in barley. Euphytica 169, 187–196. doi: 10.1007/s10681-009-9919-2
Hu, W., Yuan, Q., Wang, Y., Cai, R., Deng, X., Wang, J., et al. (2012). Overexpression of a wheat aquaporin gene, TaAQP8, enhances salt stress tolerance in transgenic tobacco. Plant Cell Physiol. 53, 2127–2141. doi: 10.1093/pcp/pcs154
Janz, D., Polle, A. (2012). Harnessing salt for woody biomass production. Tree Physiol. 32, 1–3. doi: 10.1093/treephys/tpr127
Jayakodi, M., Padmarasu, S., Haberer, G., Bonthala, V. S., Gundlach, H., Monat, C., et al. (2020). The barley pan-genome reveals the hidden legacy of mutation breeding. Nature 588, 284–289. doi: 10.1038/s41586-020-2947-8
Jiang, X., Leidi, E. O., Pardo, J. M. (2010). How do vacuolar NHX exchangers function in plant salt tolerance? Plant Signal Behav. 5, 792–795. doi: 10.4161/psb.5.7.11767
Jia, H., Wang, X., Shi, Y., Wu, X., Wang, Y., Liu, J., et al. (2020). Overexpression of medicago sativa LEA4-4 can improve the salt, drought, and oxidation resistance of transgenic arabidopsis. PloS One 15, e234085. doi: 10.1371/journal.pone.0234085
Joshi, R., Karan, R., Singla-Pareek, S. L., Pareek, A. (2016). Ectopic expression of pokkali phosphoglycerate kinase-2 (OsPGK2-p) improves yield in tobacco plants under salinity stress. Plant Cell Rep. 35, 27–41. doi: 10.1007/s00299-015-1864-z
Kobayashi, N. I., Yamaji, N., Yamamoto, H., Okubo, K., Ueno, H., Costa, A., et al. (2017). OsHKT1;5 mediates na+ exclusion in the vasculature to protect leaf blades and reproductive tissues from salt toxicity in rice. Plant J. 91, 657–670. doi: 10.1111/tpj.13595
Krausko, M., Kusa, Z., Peterkova, D., Labajova, M., Kumar, A., Pavlovic, A., et al. (2022). The absence of the AtSYT1 function elevates the adverse effect of salt stress on photosynthesis in arabidopsis. Int. J. Mol. Sci. 23, 1751. doi: 10.3390/ijms23031751
Krebs, M., Beyhl, D., Gorlich, E., Al-Rasheid, K. A., Marten, I., Stierhof, Y. D., et al. (2010). Arabidopsis V-ATPase activity at the tonoplast is required for efficient nutrient storage but not for sodium accumulation. Proc. Natl. Acad. Sci. U.S.A. 107, 3251–3256. doi: 10.1073/pnas.0913035107
Kumari, A., Bhatla, S. C. (2021). Regulation of salt-stressed sunflower (Helianthus annuus) seedling’s water status by the coordinated action of Na+/K+ accumulation, nitric oxide, and aquaporin expression. Funct. Plant Biol. 48, 573–587. doi: 10.1071/FP20334
Kumari, S., Sabharwal, V. P., Kushwaha, H. R., Sopory, S. K., Singla-Pareek, S. L., Pareek, A. (2009). Transcriptome map for seedling stage specific salinity stress response indicates a specific set of genes as candidate for saline tolerance in Oryza sativa l. Funct. Integr. Genomics 9, 109–123. doi: 10.1007/s10142-008-0088-5
Lee, E., Vanneste, S., Perez-Sancho, J., Benitez-Fuente, F., Strelau, M., Macho, A. P., et al. (2019). Ionic stress enhances ER-PM connectivity via phosphoinositide-associated SYT1 contact site expansion in arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 116, 1420–1429. doi: 10.1073/pnas.1818099116
Liu, X., Fan, Y., Mak, M., Babla, M., Holford, P., Wang, F., et al. (2017). QTLs for stomatal and photosynthetic traits related to salinity tolerance in barley. BMC Genomics 18, 9. doi: 10.1186/s12864-016-3380-0
Liu, C., Fukumoto, T., Matsumoto, T., Gena, P., Frascaria, D., Kaneko, T., et al. (2013). Aquaporin OsPIP1;1 promotes rice salt resistance and seed germination. Plant Physiol. Biochem. 63, 151–158. doi: 10.1016/j.plaphy.2012.11.018
Long, N. V., Dolstra, O., Malosetti, M., Kilian, B., Graner, A., Visser, R. G., et al. (2013). Association mapping of salt tolerance in barley (Hordeum vulgare l.). Theor. Appl. Genet. 126, 2335–2351. doi: 10.1007/s00122-013-2139-0
Lu, P., Magwanga, R. O., Guo, X., Kirungu, J. N., Lu, H., Cai, X., et al. (2018). Genome-wide analysis of multidrug and toxic compound extrusion (MATE) family in gossypium raimondii and gossypium arboreum and its expression analysis under salt, cadmium, and drought stress. G3 (Bethesda) 8, 2483–2500. doi: 10.1534/g3.118.200232
Lu, P., Magwanga, R. O., Kirungu, J. N., Hu, Y., Dong, Q., Cai, X., et al. (2019). Overexpression of cotton a DTX/MATE gene enhances drought, salt, and cold stress tolerance in transgenic arabidopsis. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00299
Mano, Y. O. U. K., Takeda, K. (1997). Mapping quantitative trait loci for salt tolerance at germination and the seedling stage in barley (Hordeum vulgare l.). Euphytica 94, 263–272. doi: 10.1023/A:1002968207362
Mao, X. Z., Cai, T., Olyarchuk, J. G., Wei, L. P. (2005). Automated genome annotation and pathway identification using the KEGG orthology (KO) as a controlled vocabulary. Bioinformatics 21, 3787–3793. doi: 10.1093/bioinformatics/bti430
Mao, G., Seebeck, T., Schrenker, D., Yu, O. (2013). CYP709B3, a cytochrome P450 monooxygenase gene involved in salt tolerance in Arabidopsis thaliana. BMC Plant Biol. 13, 169. doi: 10.1186/1471-2229-13-169
Mascher, M., Gundlach, H., Himmelbach, A., Beier, S., Twardziok, S. O., Wicker, T., et al. (2017). A chromosome conformation capture ordered sequence of the barley genome. Nature 544, 427–433. doi: 10.1038/nature22043
Mayer, K. F., Waugh, R., Brown, J. W., Schulman, A., Langridge, P., Platzer, M., et al. (2012). A physical, genetic and functional sequence assembly of the barley genome. Nature 491, 711–716. doi: 10.1038/nature11543
McMorrow, E. M., Bradbeer, J. W. (1990). Separation, purification, and comparative properties of chloroplast and cytoplasmic phosphoglycerate kinase from barley leaves. Plant Physiol. 93, 374–383. doi: 10.1104/pp.93.2.374
Munns, R., Tester, M. (2008). Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 59, 651–681. doi: 10.1146/annurev.arplant.59.032607.092911
Mwando, E., Han, Y., Angessa, T. T., Zhou, G. F., Hill, C. B., Zhang, X. Q., et al. (2020). Genome-wide association study of salinity tolerance during germination in barley (Hordeum vulgare l.). Front. Plant Sci. 11. doi: 10.3389/fpls.2020.00118
Myo, T., Tian, B., Zhang, Q., Niu, S., Liu, Z., Shi, Y., et al. (2020). Ectopic overexpression of a cotton plastidial na+ transporter GhBASS5 impairs salt tolerance in arabidopsis via increasing na+ loading and accumulation. Planta 252, 41. doi: 10.1007/s00425-020-03445-8
Nevo, E., Chen, G. (2010). Drought and salt tolerances in wild relatives for wheat and barley improvement. Plant Cell Environ. 33, 670–685. doi: 10.1111/j.1365-3040.2009.02107.x
Ning, Y. S., Jantasuriyarat, C., Zhao, Q. Z., Zhang, H. W., Chen, S. B., Liu, J. L., et al. (2011a). The SINA E3 ligase OsDIS1 negatively regulates drought response in rice. Plant Physiol. 157, 242–255. doi: 10.1104/pp.111.180893
Ning, Y. S., Xie, Q., Wang, G. L. (2011b). OsDIS1-mediated stress response pathway in rice. Plant Signal Behav. 6, 1684–1686. doi: 10.4161/psb.6.11.17916
Nutan, K. K., Kumar, G., Singla-Pareek, S. L., Pareek, A. (2018). A salt overly sensitive pathway member from brassica juncea BjSOS3 can functionally complement ΔAtsos3 in arabidopsis. Curr. Genomics 19, 60–69. doi: 10.2174/1389202918666170228133621
Ouertani, R.N., Arasappan, D., Ruhlman, T. A., Chikha., M. B., Abid, G., Mejri, S., et al. (2022). Effects of salt stress on transcriptional and physiological responses in barley leaves with contrasting salt tolerance. Int. J. Mol. Sci. 23, 5006. doi: 10.3390/ijms23095006
Qiu, L., Wu, D., Ali, S., Cai, S., Dai, F., Jin, X., et al. (2011). Evaluation of salinity tolerance and analysis of allelic function of HvHKT1 and HvHKT2 in Tibetan wild barley. Theor. Appl. Genet. 122, 695–703. doi: 10.1007/s00122-010-1479-2
Ramamoorthy, P., Lakshmanan, K., Upadhyaya, H. D., Vadez, V., Varshney, R. K. (2017). Root traits confer grain yield advantages under terminal drought in chickpea (Cicer arietinum l.). Field Crops Res. 201, 146–161. doi: 10.1016/j.fcr.2016.11.004
Ruiz, D., Martínez, V., Cerdá, V. (1997). Citrus response to salinity: Growth and nutrient uptake. Tree Physiol 17, 141–150. doi: 10.1093/treephys/17.3.141
Ruiz-Lopez, N., Perez-Sancho, J., Del, V. A., Haslam, R. P., Vanneste, S., Catala, R., et al. (2021). Synaptotagmins at the endoplasmic reticulum-plasma membrane contact sites maintain diacylglycerol homeostasis during abiotic stress. Plant Cell 33, 2431–2453. doi: 10.1093/plcell/koab122
Sawyer, G. M., Monzingo, A. F., Poteet, E. C., O'Brien, D. A., Robertus, J. D. (2008). X-Ray analysis of phosphoglycerate kinase 2, a sperm-specific isoform from mus musculus. Proteins 71, 1134–4410. doi: 1002/prot.21801
Sbei, H., Sato, K., Shehzad, T., Harrabi, M., Okuno, K. (2014). Detection of QTLs for salt tolerance in Asian barley (Hordeum vulgare l.) by association analysis with SNP markers. Breed Sci. 64, 378–388. doi: 10.1270/jsbbs.64.378
Shabala, S., Cuin, T. A., Prismall, L., Nemchinov, L. G. (2007). Expression of animal CED-9 anti-apoptotic gene in tobacco modifies plasma membrane ion fluxes in response to salinity and oxidative stress. Planta 27, 189–197. doi: 10.1007/s00425-007-0606-z
Shavrukov, Y., Gupta, N. K., Miyazaki, J., Baho, M. N., Chalmers, K. J., Tester, M., et al. (2010). HvNax3–a locus controlling shoot sodium exclusion derived from wild barley (Hordeum vulgare ssp. spontaneum). Funct. Integr. Genomics 10, 277–291. doi: 10.1007/s10142-009-0153-8
Shi, H. T., Chan, Z. L. (2014). Improvement of plant abiotic stress tolerance through modulation of the polyamine pathway. J. Integr. Plant Biol. 56, 114–121. doi: 10.1111/jipb.12128
Shi, C. C., Feng, C. C., Yang, M. M., Li, J. L., Li, X. X., Zhao, B. C., et al. (2014). Overexpression of the receptor-like protein kinase genes AtRPK1 and OsRPK1 reduces the salt tolerance of arabidopsis thaliana. Plant Sci. 63-70, 217–218, 63-70. doi: 10.1016/j.plantsci.2013.12.002
Slama, I., Abdelly, C., Bouchereau, A., Flowers, T., Savoure, A. (2015). Diversity, distribution and roles of osmoprotective compounds accumulated in halophytes under abiotic stress. Ann. Bot. 115, 433–447. doi: 10.1093/aob/mcu239
Sun, K., Fang, H., Chen, Y., Zhuang, Z., Chen, Q., Shan, T., et al. (2021). Genome-wide analysis of the Cytochrome P450 gene family involved in salt tolerance in Gossypium hirsutum. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.685054
Sun, T. J., Fan, L., Yang, J., Cao, R. Z., Yang, C. Y., Zhang, J., et al. (2019). A glycine max sodium/hydrogen exchanger enhances salt tolerance through maintaining higher na+ efflux rate and K+/Na+ ratio in arabidopsis. BMC Plant Biol. 19, 469. doi: 10.1186/s12870-019-2084-4
Sun, Z., Li, H., Zhang, Y., Li, Z., Ke, H., Wu, L., et al. (2018). Identification of SNPs and candidate genes associated with salt tolerance at the seedling stage in cotton (Gossypium hirsutum l.). Front. Plant Sci. 9. doi: 10.3389/fpls.2018.01011
Sun, M. Z., Qian, X., Chen, C., Cheng, S. F., Jia, B. W., Zhu, Y. M., et al. (2018). Ectopic expression of GsSRK in medicago sativa reveals its involvement in plant architecture and salt stress responses. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00226
Sun, X. L., Sun, M. Z., Luo, X., Ding, X. D., Cai, H., Bai, X., et al. (2013a). A glycine soja ABA-responsive receptor-like cytoplasmic kinase, GsRLCK, positively controls plant tolerance to salt and drought stresses. Planta 237, 1527–1545. doi: 10.1007/s00425-013-1864-6
Sun, X. L., Yu, Q. Y., Tang, L. L., Ji, W., Bai, X., Cai, H., et al. (2013b). GsSRK, a G-type lectin s-receptor-like serine/threonine protein kinase, is a positive regulator of plant tolerance to salt stress. J. Plant Physiol. 170, 505–515. doi: 10.1016/j.jplph.2012.11.017
Suzuki, K., Yamaji, N., Costa, A., Okuma, E., Kobayashi, N. I., Kashiwagi, T., et al. (2016). OsHKT1;4-mediated na+ transport in stems contributes to na+exclusion from leaf blades of rice at the reproductive growth stage upon salt stress. BMC Plant Biol. 16, 22. doi: 10.1186/s12870-016-0709-4
Tang, Y., Liu, X., Wang, J., Li, M., Wang, Q., Tian, F., et al. (2016). GAPIT version 2: An enhanced integrated tool for genomic association and prediction. Plant Genome 9, 1–9. doi: 10.3835/plantgenome2015.11.0120
Tu, Y., Fu, L., Wang, F., Wu, D., Shen, Q., Zhang, G. (2021). GWAS and transcriptomic integrating analysis reveals key salt-responding genes controlling na+ content in barley roots. Plant Physiol. Biochem. 167, 596–606. doi: 10.1016/j.plaphy.2021.08.038
van Bezouw, R., Janssen, E. M., Ashrafuzzaman, M., Ghahramanzadeh, R., Kilian, B., Graner, A., et al. (2019). Shoot sodium exclusion in salt stressed barley (Hordeum vulgare l.) is determined by allele specific increased expression of HKT1;5. J. Plant Physiol. 241, 153029. doi: 10.1016/j.jplph.2019.153029
Verma, V., Ravindran, P., Kumar, P. P. (2016). Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 16, 86. doi: 10.1186/s12870-016-0771-y
Wang, W., Wang, W., Wu, Y., Li, Q., Zhang, G., Shi, R., et al. (2020). The involvement of wheat U-box E3 ubiquitin ligase TaPUB1 in salt stress tolerance. J. Integr. Plant Biol. 62, 631–651. doi: 10.1111/jipb.12842
Wege, S., Qiu, J., Byrt, C., Houston, K., Waugh, R., Gilliham, M., et al. (2021). A single residue deletion in the barley HKT1;5 P189 variant restores plasma membrane localisation but not na+ conductance. Biochim. Biophys. Acta Biomembr. 1863, 183669. doi: 10.1016/j.bbamem.2021.183669
Wijmenga, C. (2021). From LD-based mapping to GWAS. Nat. Rev. Genet. 22, 480–481. doi: 10.1038/s41576-021-00366-4
Xia, Y., Yang, J., Ma, L., Yan, S., Pang, Y. (2021). Genome-wide identification and analyses of Drought/Salt-responsive Cytochrome P450 genes in medicago truncatula. Int. J. Mol. Sci. 22, 9957. doi: 10.3390/ijms22189957
Xin, Y., Wu, Y., Han, X., Xu, L. A. (2021). Overexpression of the ginkgo biloba WD40 gene GbLWD1-like improves salt tolerance in transgenic populus. Plant Sci. 313, 111092. doi: 10.1016/j.plantsci.2021.111092
Xue, W., Yan, J., Zhao, G., Jiang, Y., Cheng, J., Cattivelli, L., et al. (2017). A major QTL on chromosome 7HS controls the response of barley seedling to salt stress in the nure x tremois population. BMC Genet. 18, 79. doi: 10.1186/s12863-017-0545-z
Yang, Q., Chen, Z. Z., Zhou, X. F., Yin, H. B., Li, X., Xin, X. F., et al. (2009). Overexpression of SOS (Salt overly sensitive) genes increases salt tolerance in transgenic arabidopsis. Mol. Plant 2, 22–31. doi: 10.1093/mp/ssn058
Yemets, A. I., Krasylenko, Y. A., Lytvyn, D. I., Sheremet, Y. A., Blume, Y. B. (2011). Nitric oxide signalling via cytoskeleton in plants. Plant Sci. 181, 545–554. doi: 10.1016/j.plantsci.2011.04.017
Zhang, C., Dong, S. S., Xu, J. Y., He, W. M., Yang, T. L. (2019). PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 35, 1786–1788. doi: 10.1093/bioinformatics/bty875
Zhang, Y., Fang, J., Wu, X., Dong, L. (2018). Na+/K+ balance and transport regulatory mechanisms in weedy and cultivated rice (Oryza sativa l.) under salt stress. BMC Plant Biol. 18, 375. doi: 10.1186/s12870-018-1586-9
Zhang, H., Zhu, H., Pan, Y., Yu, Y., Luan, S., Li, L. (2014). A DTX/MATE-type transporter facilitates abscisic acid efflux and modulates ABA sensitivity and drought tolerance in arabidopsis. Mol. Plant 7, 1522–1532. doi: 10.1093/mp/ssu063
Zhao, Y., Ai, X., Wang, M., Xiao, L., Xia, G. (2016). A putative pyruvate transporter TaBASS2 positively regulates salinity tolerance in wheat via modulation of ABI4 expression. BMC Plant Biol. 16, 109. doi: 10.1186/s12870-016-0795-3
Zhou, X., Levin, E. J., Pan, Y., McCoy, J. G., Sharma, R., Kloss, B., et al. (2014). Structural basis of the alternating-access mechanism in a bile acid transporter. Nature 505, 569–573. doi: 10.1038/nature12811
Keywords: barley, GWAS, RNA-Seq, salt tolerance, Na +/K + balance, candidate genes
Citation: Xu T, Meng S, Zhu X, Di J, Zhu Y, Yang X and Yan W (2023) Integrated GWAS and transcriptomic analysis reveal the candidate salt-responding genes regulating Na+/K+ balance in barley (Hordeum vulgare L.). Front. Plant Sci. 13:1004477. doi: 10.3389/fpls.2022.1004477
Received: 27 July 2022; Accepted: 29 November 2022;
Published: 20 January 2023.
Edited by:
Zulfiqar Ali, University of Agriculture, Faisalabad, PakistanReviewed by:
Rana Muhammad Atif, University of Agriculture, Faisalabad, PakistanKlára Kosová, Crop Research Institute (CRI), Czechia
Copyright © 2023 Xu, Meng, Zhu, Di, Zhu, Yang and Yan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Yan, eWFud2VpQGphYXMuYWMuY24=