Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 03 February 2022
Sec. Plant Breeding
This article is part of the Research Topic Marker-Assisted Selection (MAS) in Crop Plants View all 45 articles

Genome-Wide Association Analysis Coupled With Transcriptome Analysis Reveals Candidate Genes Related to Salt Stress in Alfalfa (Medicago sativa L.)

\r\nFei HeFei HeChunxue WeiChunxue WeiYunxiu ZhangYunxiu ZhangRuicai LongRuicai LongMingna LiMingna LiZhen WangZhen WangQingchuan YangQingchuan YangJunmei Kang*Junmei Kang*Lin Chen*\r\nLin Chen*
  • Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China

Salt stress is the main abiotic factor affecting alfalfa yield and quality. However, knowledge of the genetic basis of the salt stress response in alfalfa is still limited. Here, a genome-wide association study (GWAS) involving 875,023 single-nucleotide polymorphisms (SNPs) was conducted on 220 alfalfa varieties under both normal and salt-stress conditions. Phenotypic analysis showed that breeding status and geographical origin play important roles in the alfalfa salt stress response. For germination ability under salt stress, a total of 15 significant SNPs explaining 9%–14% of the phenotypic variation were identified. For tolerance to salt stress in the seedling stage, a total of 18 significant SNPs explaining 12%–23% of the phenotypic variation were identified. Transcriptome analysis revealed 2,097 and 812 differentially expressed genes (DEGs) that were upregulated and 2,445 and 928 DEGs that were downregulated in the leaves and roots, respectively, under salt stress. Among these DEGs, many encoding transcription factors (TFs) were found, including MYB-, CBF-, NAC-, and bZIP-encoding genes. Combining the results of our GWAS analysis and transcriptome analysis, we identified a total of eight candidate genes (five candidate genes for tolerance to salt stress and three candidate genes for germination ability under salt stress). Two SNPs located within the upstream region of MsAUX28, which encodes an auxin response protein, were significantly associated with tolerance to salt stress. The two significant SNPs within the upstream region of MsAUX28 existed as three different haplotypes in this panel. Hap 1 (G/G, A/A) was under selection in the alfalfa domestication and improvement process.

Introduction

Soil salinization is one of the major abiotic stresses limiting agricultural production in many areas of the world (Machado and Serralheiro, 2017). As a global environmental problem, soil salinization is mainly caused by irrigation and natural mineralization (Jin et al., 2019). More than 800 million hectares of land worldwide are affected by soil salinization, which accounts for more than 6% of global land area (Shahid et al., 2018). According to the latest prediction, more than 50% of land will be affected by soil salinization by 2050 (Wang et al., 2003). Various types of ions, such as sodium, potassium, calcium, magnesium and chloride ions, can cause soil salinization (Haq et al., 2018). With respect to these ions, sodium chloride is the most abundant salt in the soil. Soil salinization can severely affect the growth and development of plants; for example, seeds cannot germinate, plant height decreases, the root system becomes shorter, and normal flowering and fruiting cannot be achieved, which eventually leads to a decline in yield and quality (Deinlein et al., 2014). Moreover, some salt-sensitive plants or crops can die even under slight soil salt stress. Modern molecular genetics has proven that the main reason for this phenomenon is that under salt stress, the normal physiological and metabolic functions of plants are damaged by ionic, osmotic and oxidative stress (Zhang et al., 2010).

Alfalfa (Medicago sativa L.), an important legume crop species, plays an important role in improving saline-alkali soil and agricultural production because this species can fix atmospheric nitrogen (Shen et al., 2020). However, severe salt stress greatly reduces the growth and productivity of alfalfa. Therefore, the development of salt-tolerant alfalfa varieties is urgently needed for the sustainable production of alfalfa worldwide.

In the long-term evolutionary process, plants, especially some salt-tolerant plants, have developed adaptive mechanisms to address damage caused by salt stress (Zhao et al., 2020). Previous studies on model plant species such as Arabidopsis, rice and maize have revealed that transcription factors (TFs) and microRNAs play a large role in salt tolerance (Kreps et al., 2002 and Zhao et al., 2020). These molecular modules form a large gene regulatory network in plants to cope with salt stress (Zelm et al., 2020). Within this gene regulatory network, the plant response to high salt stress involves a variety of signal transduction pathways, including two main ones: an abscisic acid (ABA)-dependent pathway and an ABA-independent pathway. In the ABA-dependent pathway, many TFs, such as MYB (Dubos et al., 2010; Li J. et al., 2019 and Zhao et al., 2019), bZIP (Zhang M. et al., 2020; Bi et al., 2021), NAC (Joshi et al., 2016; Khedia et al., 2018; Li et al., 2021; Zhang X. et al., 2021) TFs, which have ABA-responsive element-binding factors in their promoter regions, play an important role in the plant response to salt stress. For example, Arabidopsis NAC069 mutants are insensitive to salt and cold stresses (He et al., 2017). Dehydration responsive element binding protein (DREB) TFs are involved in the salt stress response of plants mainly through the ABA-independent pathway (Zhao et al., 2016). Many studies have revealed that overexpression of C-repeat binding factor (CBF) enhances plant tolerance to salt stress, which indicates that CBFs may play a positive role in salt stress tolerance (Li et al., 2011; Yin et al., 2018; Verma et al., 2020).

Little is known about the genetic architecture, molecular mechanism and physiology underlying the salt tolerance of alfalfa compared with other crop species. miR156, which is an important regulator controlling flowering time and plant architecture, has been studied in depth in several plant species (Wang et al., 2009; Wei et al., 2018; Jerome Jeyakumar et al., 2020). Recent studies in alfalfa have shown that miR156 plays a key role in the salt stress response by altering the expression of its target genes, such as SPL gene family members (Aung et al., 2015; Arshad et al., 2017a,b, 2018; Feyissa et al., 2021). Other genes such as MsLEA4-4 (Jia et al., 2020) and MsCBL4 (An et al., 2020) involved in the salt response have been identified and function in alfalfa.

The use of genome-wide association studies (GWASs) constitutes an effective way to study genetic variations and important complex traits, including those related to salt stress and drought stress tolerance, flowering time and yield in maize, rice and soybean (Qi et al., 2014; Shi et al., 2017; Yu et al., 2017; Zeng et al., 2017; Patishtan et al., 2018; Luo et al., 2019, 2021). In alfalfa, many functional loci associated with salt stress tolerance, drought stress tolerance, forage quality, and disease resistance have been identified via GWASs and quantitative trait loci (QTLs) (Yu et al., 2016, 2020; Biazzi et al., 2017; Jia et al., 2017; Liu and Yu, 2017). However, due to the lack of a reference genome, these significantly associated loci can be aligned to the reference genome of only Medicago truncatula (Wang et al., 2020). Therefore, some key genes for these traits cannot be found or identified. RNA sequencing (RNA-seq) is a powerful tool for discovering candidate genes and pathways to target traits and has been used widely in many plant (including crop plant) species (Cui et al., 2016; Fang et al., 2017; Yao et al., 2020). The alfalfa reference genomes that have been published in 2020 can greatly improve the accuracy and precision of GWASs and RNA-seq analysis (Chen et al., 2020; Shen et al., 2020; Long et al., 2021).

In this study, we first conducted a GWAS analysis of salt-related traits among a panel that consisted of 220 alfalfa varieties. Then, RNA-seq analysis was used to identify the genes involved in the salt response in alfalfa. Finally, we identified candidate genes associated with salt stress tolerance through a combination of a GWAS and RNA-seq analysis. The results provide valuable information to understand the genetic basis and molecular mechanism of the salt stress response in alfalfa.

Materials and Methods

Plant Materials and Phenotypic Analysis

The GWAS panel used in this study consisted of 220 alfalfa varieties from fifty countries throughout the world (Supplementary Table 1). The seeds of this panel were collected from the national grass seed resources of China and the national plant germplasm of the U.S. The plants used for phenotyping in this study were collected from the U.S. National Plant Germplasm System. For tolerance to salt stress, the seedlings of this panel were grown in 20 × 5 cm containers in a greenhouse. The fresh weight of every variety (10 individual plants) under normal conditions and salt-stress conditions was recorded. The tolerance to salt stress of every variety was evaluated by comparing the fresh weight under salt stress with that of an unstressed control; these evaluations were repeated three times for each variety. A relatively high fresh weight suggests that a variety has a relatively high salt tolerance. Germination ability under salt stress was conducted as described in a previous study (Rumbaugh and Pendery, 1990). A relatively high absolute value of germination ability under salt stress suggests that a variety has a relatively high salt tolerance in the germination stage.

For RNA-seq, the Zhongmu No. 1 alfalfa variety was used in this study. Seeds of this variety were placed in a petri dish and subjected to 4°C for 3 days. Then, the seeds were grown in a greenhouse at 24°C (day)/20°C (night) with a 16 h (day)/8 h (night) photoperiod for one week. Afterward, the plants were transplanted into pots and allowed to grow for one week. The alfalfa seedlings were subsequently treated with NaCl solution (300 mM) for 14 days. The roots and leaves of 10 individual alfalfa seedlings under normal conditions and salt-stress conditions were collected. Two replicates were included in this study.

Genome-Wide Association Study

The genotype of this panel was published in our previous study (Chen et al., 2021). A total of 875,023 high-quality single-nucleotide polymorphism (SNPs) based on the Zhongmu 4 reference haploid genome were used in this study. A GWAS analysis was performed using TASSEL software with the general linear model (GLM) and the mixed linear model (MLM) in conjunction with population genetic structure (Q) and kinship (K), which were analyzed in a previous study (Bradbury et al., 2007; Chen et al., 2021). The thresholds for significantly associated loci were a logarithm of odds (LOD) score ≥ 4.5 (GLM) and 5 (MLM).

RNA-seq and Transcription Analysis

Total RNA was extracted by using TRIzol reagent (Invitrogen, CA, United States). cDNA library construction and sequencing via an Illumina HiSeq 2000 platform were conducted by Novogene (Beijing, China). The clean reads obtained by the standard process were mapped to the Zhongmu No. 1 genome by using TopHat2 software (Trapnell et al., 2012). Only the reads that were unique hits were retained for further analysis. The expected number of fragments per kilobase of transcript sequence per million base pairs sequenced (FPKM) was used to estimate gene expression levels. The Pearson correlation coefficient between replications of different samples was estimated according to the log10(FPKM + 1) value. Genes with different expression levels were identified by using DESeq, and Benjamini and Hochberg’s approach was used to control the false discovery rate (FDR) (Benjamini and Hochberg, 1995; Anders and Huber, 2010). Significant DEGs were obtained according to Padj < 0.05 and | log2(fold change [FC])| ≥ 1.

Candidate Gene Analysis

The reported alfalfa Zhongmu No. 1 and Zhongmu 4 reference haploid genomes were used to identify genes with significantly associated loci. The UniProtKB1 and The Arabidopsis Information Resource (TAIR2) databases were used for gene annotation. All genes within 200 kb (100 kb up- and downstream) of the significant loci were identified according to the linkage disequilibrium (LD) of the association panel. The genes near significantly associated loci received increased attention. The results of the GWAS and RNA-seq analysis were combined to identify candidate genes for salt stress tolerance in alfalfa.

Quantitative Real-Time RT-PCR

Quantitative reverse transcription-PCR (RT-PCR) was carried out in triplicate for each sample using SYBR Premix Ex Taq (Takara, Japan) on a 7500 Real-Time PCR System (Applied Biosystems, CA, United States). The gene expression levels in these samples were normalized to that of alfalfa GAPDH, and the 2–ΔΔCt method was used to calculate relative gene expression levels. The primers used in this study are listed in the Supplementary Tables.

Results

Analysis of Phenotypic Variation

Phenotypic data of plant tolerance to salt stress and germination ability under salt stress were obtained from the USDA national database3. The data showed that the two traits were normally distributed, ranging from 0.599 to 2.616 for tolerance to salt stress and from −0.99 to −0.27 for germination ability under salt stress (Supplementary Figures 1, 2). Correlation analysis showed there was no correlation between the two traits (r = −0.01), which means that the genetic mechanisms underlying salt tolerance between the germination and seedling stages are different.

Our association panel could be divided into three subgroups according to breeding status: Wild, Landrace and Cultivar. Statistical analysis of the tolerance to salt stress showed that there were significant differences among the three subgroups (R2 = 6.20%; P = 0.01) and that the wild materials had the highest tolerance to salt stress in the seedling stage (Figure 1A). For germination ability under salt stress, there was no significant difference among the three subgroups (R2 = 1.99%; P = 0.21), and compared with the landrace and wild materials, the cultivated materials had a higher germination rate under salt stress (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1. Box plot showing the variation in tolerance to salt stress and germination ability under salt stress evaluated for the association panel. (A) Box plot of the tolerance to salt stress among subgroups (Wild, Landrace, Cultivar) according to breeding status. (B) Box plot of germination ability under salt stress among subgroups (Wild, Landrace, Cultivar). (C) Box plot of the tolerance to salt stress among subgroups (Europe, Asia, America and Africa) according to geographical origin. (D) Box plot of germination ability under salt stress among subgroups (Europe, Asia, America and Africa). The different letters above the boxplots in (C,D) indicate significant differences at the P < 0.05 level according to Duncan’s multiple comparisons test.

The association panel could also be divided into four subgroups according to geographical origin: Asia, Africa, Europe and America. For tolerance to salt stress, there was no significant difference among the four subgroups (R2 = 3.51%; P = 0.16). The alfalfa varieties from the African region had the lowest tolerance to salt stress among the four subgroups (P < 0.01, Figure 1C); the European and Asian subgroups showed a higher tolerance to salt stress. For germination under salt stress, there were significant differences among the four subgroups (R2 = 7.55%; P = 3.60 × 10–3), and the alfalfa varieties from Asian and African regions had a higher germination rate under salt stress than did the varieties from European and American regions (P < 0.01, Figure 1D). Integrating these results, we found that breeding status had a significant effect on tolerance to salt stress at the seedling stage, while geographical origin mainly affected germination under salt stress. In addition, compared with those from the other regions, the alfalfa varieties from the Asian region had a higher tolerance to salt stress at the seedling stage and a higher germination under salt stress, which may be related to the growth environment of alfalfa in Asia.

Genome-Wide Association Study of Salt Stress Tolerance

Based on the resequencing results, we obtained 875,023 high-quality SNPs distributed across 8 alfalfa chromosomes. The significant SNPs were identified by using the GLM and MLM for tolerance to salt stress and germination ability under salt stress. When the GLM was used, a total of 20 and 24 significant loci were found according to the LOD > 5.00 threshold for tolerance to salt stress (Supplementary Figures 3A,B and Supplementary Table 2) and germination ability under salt stress (Supplementary Figures 3C,D and Supplementary Table 2).

To avoid false-positive results, the MLM was conducted for the two traits. Under the LOD > 4.50 threshold, in total, 18 and 15 significant SNPs were found for tolerance to salt stress and germination ability under salt stress, respectively (Figures 2A,B and Table 1). For tolerance to salt stress, the proportion of phenotypic variance explained (PVE) by the 18 significant SNPs ranged from 12 to 23% (Table 1). The SNP located at 37.34 Mb on chromosome 2 had the largest PVE (R2 = 23%) for tolerance to salt stress. For germination ability under salt stress, the PVE by the 15 significant SNPs ranged from 9 to 14%, and the SNP located at 29.42 Mb at chromosome 5 had the largest PVE (R2 = 14%). Taken together, these results provide valuable information for salt tolerance-related gene mining and molecular breeding for alfalfa in the future.

FIGURE 2
www.frontiersin.org

Figure 2. GWAS of tolerance to salt stress and germination ability under salt stress. (A,B), Manhattan plots of the GWAS results of the tolerance to salt stress and germination ability under salt stress among the association panel. (C,D), Correlations between phenotypic traits (C) for tolerance to salt stress and (D) for germination ability under salt stress and number of excellent alleles identified by the GWAS analysis.

TABLE 1
www.frontiersin.org

Table 1. Significant SNPs identified via GWAS.

To further verify the accuracy of the GWAS results, we analyzed the correlations between the number of excellent alleles carried by each material and the two-salt stress-related traits within our association panel. The results showed that the correlation between the number of excellent alleles and tolerance to salt stress was significantly positive, with a coefficient of 0.82 (Figure 2C). For example, Irmsskaja alfalfa from the European region has 16 excellent alleles, and its tolerance to salt stress is 1.97. For germination ability under salt stress, a similar result was found (Figure 2D). The No. 12835 alfalfa from the Asian region has 14 excellent alleles, and its germination ability under salt stress is −0.99, which means that this alfalfa variety has a higher germination under salt stress. These results show that our GWAS results are accurate and can be used for further analysis.

RNA-seq for Salt Stress Tolerance

We sampled the leaves and roots of 3 Zhongmu No. 1 plants growing for 21 days to construct cDNA libraries and conduct transcriptome sequencing. In total, 2,097 and 812 differentially expressed genes (DEGs) were identified as being upregulated in the leaves and roots, respectively, and 2445 and 928 DEGs were downregulated, respectively (Figure 3A and Supplementary Tables 36). These DEGs were distributed across the alfalfa genome, as shown in Figure 3B. Among the 2,909 (2,097 + 812) upregulated DEGs, 242 were common to the leaves and roots under salt stress (Figure 3C and Supplementary Table 7). Among the 3,373 (2,445 + 928) downregulated DEGs, a total of 154 were common to the leaves and roots under salt stress (Figure 3D and Supplementary Table 8). An overview of the expression of these DEGs is shown in a heatmap, and their expression could be divided into six patterns (Figure 3E). Gene Ontology (GO) analysis showed that these DEGs could be divided into three GO categories: biological processes, cellular components and molecular functions (Supplementary Figure 4). In the biological process category, the three largest subcategories were metabolic processes, cellular processes and single-organism processes. For molecular functions, the two most common subcategories were binding and catalytic activity.

FIGURE 3
www.frontiersin.org

Figure 3. Transcriptome analysis. (A) Number of upregulated genes (red) and downregulated genes (blue) between alfalfa seedlings under salt stress and those under normal conditions. (B) Distribution of DEGs in the leaves and roots of alfalfa. The inner circle indicates the DEGs in the roots, and the outer circle indicates the DEGs in the leaves. (C) Upregulated DEGs common to the leaves and roots. (D) Downregulated DEGs common to the leaves and roots. (E) Heatmap clustering of all the DEGs according to their expression abundance. The different colors indicate different levels of expression abundance.

Among these DEGs, many encoding TFs were identified, including the AP2, CBF, NAC and WRKY genes (Figure 4A). A total of eight genes encoding MYB TFs were detected; two of the MYB TF-encoding genes (MsG0180000292.01 and MsG0880045184.01) were downregulated in the leaves under salt stress, and MsG0280010554.01 was upregulated. The expression of the other two genes encoding MYB TFs changed: MsG0180003989.01 was upregulated in the roots under salt stress, and MsG0580025527.01 was downregulated (Figure 4B). Three genes encoding CBF TFs were found, two of which (MsG0380016957.01 and MsG0580029966.01) were downregulated in leaves under salt stress, and another CBF TF-encoding gene (MsG0280009550.01) was upregulated in the roots (Figure 4C). Two genes encoding bZIP TFs were found, and their expression changed only in the leaves. MsG0480020539.01, one of the bZIP TF-encoding genes, was significantly upregulated in alfalfa leaves under salt stress, and MsG0580027503.01, another bZIP TF-encoding gene, was significantly downregulated (Figure 4D). In addition, we found that the expression of 20 NAC genes changed under salt stress. Fourteen NAC genes were upregulated in the leaves, 6 NAC genes were downregulated in the leaves, and 5 NAC genes were upregulated in the roots of alfalfa (Figure 4E). Taken together, these results showed that these TFs may play an important role in alfalfa under salt stress.

FIGURE 4
www.frontiersin.org

Figure 4. Transcriptome analysis of differentially expressed TFs. (A) Number of differentially expressed TFs. (B) Differentially expressed MYB-encoding genes. (C) Differentially expressed CBF-encoding genes. (D) Differentially expressed bZIP-encoding genes. (E) Differentially expressed NAC-encoding genes. The color in each cell in panel (B–D) indicates the value of expression abundance.

Candidate Gene Analysis for Salt Stress in Alfalfa

We combined the results of our GWAS and DEG analysis to further analyze the genes associated with salt stress tolerance in alfalfa. According to the LD, a total of 366 candidate genes (211 genes for tolerance to salt stress and 155 genes for germination ability under salt stress) were found for the 33 (18 for tolerance to salt stress and 15 for germination ability under salt stress) significant SNPs (Supplementary Table 9). For the 155 genes associated with germination ability under salt stress, 11 were differentially expressed. Among the 211 candidate genes for tolerance to salt stress, 27 genes were differentially expressed under salt stress. Most of the significant SNPs were located in intergenic spacer regions. To further identify the candidate genes for salt stress in alfalfa, the DEGs closest to the significantly associated SNPs were analyzed by candidate gene association mapping. According to this criterion, there were 5 DEGs associated with tolerance to salt stress and 3 DEGs associated with germination ability under salt stress for analysis (Table 2). To validate the RNA-seq results, we designed primers and performed qRT-PCR on the eight candidate genes (Figure 5 and Supplementary Table 10). The results showed that the RNA-seq results and qRT-PCR results were highly consistent. Among the 8 candidate genes, Msa0033600 (referred to as MsG0180003906.01 in the Zhongmu No. 1 genome annotation), which encodes auxin response protein 28 and is named MsAUX28, was selected for analysis because there were enough SNP sites around this gene.

TABLE 2
www.frontiersin.org

Table 2. Candidate genes associated with tolerance to salt stress and germination ability under salt stress.

FIGURE 5
www.frontiersin.org

Figure 5. qRT-PCR analysis of eight candidate genes associated with the salt stress response in alfalfa. “**” means a significant difference at the P < 0.01 level.

We found that two SNPs located in the upstream region (∼10 kb) of the MsAUX28 gene were significantly associated with tolerance to salt stress (Figure 6A). The two significant SNPs in the upstream region of MsAUX28 existed as three different haplotypes in the association panel (Figure 6B). Hap1 (G/G, A/A) had the lowest frequency in the panel but had the highest tolerance to salt stress (Figure 6C). In addition, we found that Hap1 did not exist in the wild alfalfa materials and that the frequency of Hap1 increased in the cultivated materials, which suggests that Hap1 was under selection in the alfalfa breeding process (Figure 6D).

FIGURE 6
www.frontiersin.org

Figure 6. MsAux28 involvement in the salt stress response in alfalfa. (A) Association analysis of the SNPs located within the MsAUX28 gene region and tolerance to salt stress and germination ability under salt stress. R2 means the phenotypic variation that can be explained by this locus. (B) Haplotype analysis of the two significant SNPs. (C) Comparison of salt stress tolerance the three haplotypes. And, there are three haplotypes (Hap1: G/G:A/A; Hap2: G/G:G/A; Hap3: G/A:G/A). (D) Allele frequency of the three haplotypes in the Wild, Landrace and Cultivar subgroups.

Discussion

Alfalfa is often planted on saline-alkali land in Asian countries such as China. Therefore, it is necessary to clarify the genetic basis of alfalfa salt tolerance, which could help improve the salt tolerance of this species (Yu et al., 2016). In this study, GWAS analysis based on resequencing was conducted to dissect the genetic architecture of salt stress-related traits, including germination ability under salt stress and tolerance to salt stress. Multiple significant SNPs were identified for these two salt-related traits. Furthermore, several candidate genes involved in salt tolerance were analyzed by a combination of transcriptome analysis and a GWAS. These candidate genes can serve as potential targets of the molecular mechanism of salt tolerance in alfalfa.

Salt stress severely affects the growth and development of alfalfa (Yu et al., 2016; Liu and Yu, 2017). Salt stress causes a decrease in seed germination, plant height and root development in alfalfa (El-Sharkawy et al., 2017; Lei et al., 2018). However, few studies have focused on the effects of salt stress on alfalfa varieties from different geographical regions or on the improvement of breeding status. Based on our results, we found a significant difference in germination ability under salt stress among the alfalfa varieties from different geographical regions. However, there was no difference in this trait among the alfalfa varieties with different breeding status improvements. Compared with those from other geographical regions, the alfalfa varieties from the Asian region had a higher germination ability under salt stress. The reason for this phenomenon may be related to most alfalfa in Asia, especially in China, being planted in saline-alkali soil because most good-quality farmland is used to grow staple crop species such as maize, rice and wheat for human consumption. However, this phenomenon is rare in European and American regions. We also found that, compared with the landrace and cultivated alfalfa, the wild alfalfa varieties were more tolerant to salt stress. These results indicate that many salt tolerance-related genes exist in wild populations. In the future, more salt tolerance-related genes should be mined and transferred to cultivated varieties to improve their salt tolerance. This phenomenon is widespread in other plant species such as maize and rice (Xu et al., 2019; Xu and Sun, 2021).

Many GWASs have been performed on salt stress-related traits in alfalfa. A total of 36 significant markers were found to be associated with salt tolerance during germination under different levels of salt treatments when 198 alfalfa cultivars and landraces were subjected to genotyping by sequencing (GBS) (Yu et al., 2016). Using the same population, Liu and Yu (2017) identified 42 markers significantly associated with salt tolerance-related traits, such as dry weight, plant height, stomatal conductance and leaf chlorophyll content. Due to the lack of a reference genome for alfalfa, these significant markers were aligned only to the reference genome (Medicago truncatula) (Wang et al., 2020). Therefore, much important information or presence of significant markers may be lost. In 2020, several reference genomes of alfalfa were published, which allowed the ability to mine candidate genes for salt stress based on GWAS results (Chen et al., 2020; Shen et al., 2020; Long et al., 2021). Moreover, based on the resequencing results, the number of SNPs in this study is more than that in previous studies in which SNPs were identified based on GBS technology (Yu et al., 2016; Liu et al., 2017). In this study, 18 and 15 significant SNPs were identified as being associated with tolerance to salt stress and germination ability under salt stress, respectively. No consistent loci were detected, however, which suggests that there are different genetic structures and molecular mechanisms underlying germination ability under salt stress and salt tolerance in the seedling stage for alfalfa.

Transcriptome analysis has been proven to be a powerful tool for identifying genes associated with plant growth and development. Recently, RNA-seq technology was used to identify key genes related to many important traits in alfalfa, such as fall dormancy (via transcriptome profiling of gene expression in fall dormant and non-dormant alfalfa) (Liu et al., 2019), freezing stress (via deep-sequencing transcriptome analysis of field-grown crown buds acclimated to freezing stress) (Song et al., 2016) and aluminum stress (via a transcriptome analysis of candidate genes potentially involved in the Al stress response) (Liu et al., 2017). In this study, 2,097 and 812 DEGs were upregulated and 2,445 and 928 DEGs were downregulated in the leaves and roots under salt stress, respectively. These genes can help us understand the molecular mechanism of salt tolerance in alfalfa. Previous studies have revealed many genes and pathways controlling salt stress in plants, such as some TF-encoding genes and plant hormone signal pathways (Zhao et al., 2020). We also detected some important TF-encoding genes that were significantly differently expressed under salt stress in alfalfa leaves and roots. In Arabidopsis, the AtMYB49 gene is induced in response to salinity and acts as a positive regulator of salt stress (Zhang P. et al., 2020). MYB TFs are involved in the salt tolerance pathway in many plant species, such as maize, rice and apple (Wang et al., 2014; Tang et al., 2019; Wu et al., 2019). In this study, we also found that MsG0280010554.01, which encodes a MYB TF in the alfalfa genome, was upregulated in the leaves under salt stress. In addition to MYB TFs, NAC TFs also play an important role in plant drought and salt stress responses. For example, GmSIN1, which encodes a NAC TF in soybean, promotes root growth and salt tolerance and improves grain yield under salt stress (Li S. et al., 2019). Based on our results, we also found that many NAC TFs were upregulated in alfalfa leaves under salt stress, which suggests that these NAC TFs play an important role in alfalfa salt tolerance.

Many previous studies in other plant species have proven that the combination of GWAS and RNA-seq results can be used as an effective way to identify candidate loci or genes. For example, five important genes were identified for kernel row number in maize by a combination of transcriptome analysis and regional association mapping (An et al., 2019). Seven candidate genes for root-related traits in maize were also identified via this strategy (Guo et al., 2019). To our knowledge, this study is the first to use this strategy to identify candidate genes related to salt tolerance in alfalfa. Here, we identified 8 DEGs close to the six significant SNPs (four SNPs located within four DEG sequences and the other two SNPs located in intergenic regions) for germination ability under salt stress and salt tolerance to salt stress.

Plant hormones such as auxin, gibberellic acid (GA), and ABA also play important roles in the salt stress response (Zhao et al., 2020). In rice, OsIAA20, which encodes an Aux/IAA protein, improves rice salt tolerance through an ABA-related pathway (Zhang A. et al., 2021). Overexpression of the auxin signaling F-Box 3 receptor (AFB3) can increase salt tolerance in Arabidopsis through the auxin pathway (Garrido-vargas et al., 2020). In tomato, a mutation in SlARF4, which encodes an auxin response factor, results in improved tolerance to salt stress and osmotic stress (Bouzroud et al., 2020). Here, we also found that MsG0180003906.01 (MsAux28), which encodes an Aux/IAA protein, was significantly associated with salt tolerance in alfalfa. The expression of MsAux28 was downregulated in alfalfa leaves under salt stress. Therefore, both the overexpression and the knockdown of this gene to validate the function of this gene in alfalfa under salt stress are needed.

In our study, we used a combination of a GWAS and RNA-seq analysis to mine important genes involved in salt tolerance in alfalfa. The results could provide valuable information on the critical regulators of salt tolerance in alfalfa molecular breeding.

Data Availability Statement

RNA sequence data from roots and leaves after 14 days of salt treatment in a greenhouse has been submitted to The NCBI Sequence Read Archive (BioProject: PRJNA777963). The raw sequencing data can be obtained from the National Genomics Data Center (PRJCA004024, https://bigd.big.ac.cn/).

Author Contributions

JK and LC designed the study and revised the manuscript. FH, RL, ML, ZW, CW, and YZ analyzed these data and generated the figures and tables. FH and LC analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.

Funding

The authors thank the reviewers for their valuable comments on this manuscript and gratefully acknowledge the financial support for this study provided by the Ministry of Science and Technology of China (2017YFE0111000/EUCLEG 727312), the Key Research Project of Ningxia Province for the Alfalfa Breeding Program (2019NYY203), the earmarked fund for the China Agriculture Research System (CARS-35-04), and the Agricultural Science and Technology Innovation Program (ASTIP-IAS14), and the Key Projects in Science and Technology of Inner Mongolia (2021ZD0031).

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.2021.826584/full#supplementary-material

Supplementary Figure 1 | Distribution of germination ability under salt stress among the association panel.

Supplementary Figure 2 | Distribution of tolerance to salt stress among the association panel.

Supplementary Figure 3 | GWAS analysis of tolerance to salt stress and germination ability under salt stress via a GLM.

Supplementary Figure 4 | GO analysis of the total DEGs.

Supplementary Table 1 | The association panel related information used in this study.

Supplementary Table 2 | Significant SNPs identified via a GLM.

Supplementary Table 3 | Upregulated DEGs in leaves under salt stress.

Supplementary Table 4 | Downregulated DEGs in leaves under salt stress.

Supplementary Table 5 | Upregulated DEGs in roots under salt stress.

Supplementary Table 6 | Downregulated DEGs in roots under salt stress.

Supplementary Table 7 | Upregulated genes common to both leaves and roots under salt stress.

Supplementary Table 8 | Downregulated genes common to both leaves and roots under salt stress.

Supplementary Table 9 | Candidate genes identified via GWAS and transcriptome analysis.

Supplementary Table 10 | Primers used in this study.

Footnotes

  1. ^ https://www.UniProt.org/
  2. ^ https://www.arabidopsis.org/
  3. ^ https://www.usda.gov/

References

An, Y., Chen, L., Li, Y. X., Li, C., Shi, Y., Song, Y., et al. (2019). Candidate loci for the kernel row number in maize revealed by a combination of transcriptome analysis and regional association mapping. BMC Plant Biol. 19:201. doi: 10.1186/s12870-019-1811-1

PubMed Abstract | CrossRef Full Text | Google Scholar

An, Y., Yang, X. X., Zhang, L., Zhang, J., Du, B., Yao, L., et al. (2020). Alfalfa MsCBL4 enhances calcium metabolism but not sodium transport in transgenic tobacco under salt and saline-alkali stress. Plant Cell Rep. 39, 997–1011. doi: 10.1007/s00299-020-02543-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1038/npre.2010.4282.1

CrossRef Full Text | Google Scholar

Arshad, M., Feyissa, B. A., Amyot, L., Aung, B., and Hannoufa, A. (2017a). MicroRNA156 improves drought stress tolerance in alfalfa (Medicago sativa) by silencing SPL13. Plant Sci. 258, 122–136. doi: 10.1016/j.plantsci.2017.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Arshad, M., Gruber, M. Y., Wall, K., and Hannoufa, A. (2017b). An Insight into microRNA156 role in salinity stress responses of Alfalfa. Front. Plant Sci. 8:356. doi: 10.3389/fpls.2017.00356

PubMed Abstract | CrossRef Full Text | Google Scholar

Arshad, M., Gruber, M. Y., and Hannoufa, A. (2018). Transcriptome analysis of microRNA156 overexpression alfalfa roots under drought stress. Sci Rep. 8:9363. doi: 10.1038/s41598-018-27088-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Aung, B., Gruber, M. Y., Amyot, L., Omari, K., Bertrand, A., and Hannoufa, A. (2015). MicroRNA156 as a promising tool for alfalfa improvement. Plant Biotechnol. J. 13, 779–790. doi: 10.1111/pbi.12308

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Stat. Methodol. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Bi, C., Yu, Y., Dong, C., Yang, Y., Zhai, Y., Du, F., et al. (2021). The bZIP transcription factor TabZIP15 improves salt stress tolerance in wheat. Plant Biotechnol. J. 19, 209–211. doi: 10.1111/pbi.13453

PubMed Abstract | CrossRef Full Text | Google Scholar

Biazzi, E., Nazzicari, N., Pecetti, L., Brummer, E. C., Palmonari, A., Tava, A., et al. (2017). Genome-wide association mapping and genomic selection for alfalfa (Medicago sativa) forage quality traits. PLoS One 12:e0169234. doi: 10.1371/journal.pone.0169234

PubMed Abstract | CrossRef Full Text | Google Scholar

Bouzroud, S., Gasparini, K., Hu, G., Barbosa, M. A. M., Rosa, B. L., Fahr, M., et al. (2020). Down regulation and loss of auxin response factor 4 function using CRISPR/Cas9 alters plant growth, stomatal function and improves tomato tolerance to salinity and osmotic stress. Genes 11:272. doi: 10.3390/genes11030272

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). Tassel: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–5. doi: 10.1093/bioinformatics/btm308

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Zeng, Y., Yang, Y., Huang, L., Tang, B., Zhang, H., et al. (2020). Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat. Commun. 11:2494. doi: 10.1038/s41467-020-16338-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., He, F., Long, R. C., Zhang, F., Li, M. N., Wang, Z., et al. (2021). A global alfalfa diversity panel reveals genomic selection signatures in Chinese varieties and genomic associations with root development. J. Integr. Plant Biol. 63, 1937–1951. doi: 10.1111/jipb.13172

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, Z., Luo, J., Qi, C., Ruan, Y., Li, J., Zhang, A., et al. (2016). Genome-wide association study (GWAS) reveals the genetic architecture of four husk traits in maize. BMC Genomics 17:946. doi: 10.1186/s12864-016-3229-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Deinlein, U., Stephan, A. B., Horie, T., Luo, W., Xu, G., and Schroeder, J. I. (2014). Plant salt-tolerance mechanisms. Trends Plant Sci. 19, 371–379. doi: 10.1016/j.tplants.2014.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubos, C., Stracke, R., Grotewold, E., Weisshaar, B., Martin, C., and Lepiniec, L. (2010). MYB transcription factors in Arabidopsis. Trends Plant Sci. 15, 573–581. doi: 10.1016/j.tplants.2010.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Sharkawy, M., El-Beshsbeshy, T., Al-Shal, R., and Missaoui, A. (2017). Effect of plant growth stimulants on alfalfa response to salt stress. Agric. Sci. China 8, 267–291. doi: 10.4236/as.2017.84020

CrossRef Full Text | Google Scholar

Fang, L., Sahana, G., Su, G., Yu, Y., Zhang, S., Lund, M. S., et al. (2017). Integrating sequence-based GWAS and RNA-Seq provides novel insights into the genetic basis of mastitis and milk production in dairy cattle. Sci. Rep. 7:45560. doi: 10.1038/srep45560

PubMed Abstract | CrossRef Full Text | Google Scholar

Feyissa, B. A., Amyot, L., Nasrollahi, V., Papadopoulos, Y., Kohalmi, S. E., and Hannoufa, A. (2021). Involvement of the miR156/SPL module in flooding response in Medicago sativa. Sci Rep. 11:3243. doi: 10.1038/s41598-021-82450-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Garrido-vargas, F., Godoy, T., Tejos, R., and O’Brien, J. (2020). Overexpression of the auxin receptor AFB3 in Arabidopsis results in salt stress resistance and the modulation of NAC4 and SZF1. Int. J. Mol. Sci. 21:9528. doi: 10.3390/ijms21249528

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, J., Li, C., Zhang, X., Li, Y. X., Zhang, D., Shi, Y., et al. (2019). Transcriptome and GWAS analyses reveal candidate gene for seminal root length of maize seedlings under drought stress. Plant Sci. 292:110380. doi: 10.1016/j.plantsci.2019.110380

PubMed Abstract | CrossRef Full Text | Google Scholar

Haq, S., Kaushal, S. S., and Duan, S. (2018). Episodic salinization and freshwater salinization syndrome mobilize base cations, carbon, and nutrients to streams across urban regions. Biogeochemistry 141, 463–486. doi: 10.1007/s10533-018-0514-2

CrossRef Full Text | Google Scholar

He, L., Shi, X., Wang, Y., Guo, Y., Yang, K., and Wang, Y. (2017). Arabidopsis ANAC069 binds to C[A/G]CG[T/G] sequences to negatively regulate salt and osmotic stress tolerance. Plant Mol. Biol. 93, 369–387. doi: 10.1007/s11103-016-0567-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Jerome Jeyakumar, J. M., Ali, A., Wang, W. M., and Thiruvengadam, M. (2020). Characterizing the role of the miR156-SPL network in plant development and stress response. Plants 9:1206. doi: 10.3390/plants9091206

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, C., Wu, X., Chen, M., Wang, Y., Liu, X., Gong, P., et al. (2017). Identification of genetic loci associated with crude protein and mineral concentrations in alfalfa (Medicago sativa) using association mapping. BMC Plant Biol. 17:97. doi: 10.1186/s12870-017-1047-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, H., Wang, X., Shi, Y., Wu, X., and Dong, K. (2020). Overexpression of Medicago sativa LEA4-4 can improve the salt, drought, and oxidation resistance of transgenic Arabidopsis. PLoS One 15:e0234085. doi: 10.1371/journal.pone.0234085

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, X., Liu, T., Xu, J., Gao, Z., and Hu, X. (2019). Exogenous GABA enhances muskmelon tolerance to salinity-alkalinity stress by regulating redox balance and chlorophyll biosynthesis. BMC Plant Biol. 19:48. doi: 10.1186/s12870-019-1660-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Joshi, R., Wani, S. H., Singh, B., Bohra, A., Dar, Z. A., Lone, A. A., et al. (2016). Transcription factors and plants response to drought stress: current understanding and future directions. Front. Plant Sci. 7:1029. doi: 10.3389/fpls.2016.01029

PubMed Abstract | CrossRef Full Text | Google Scholar

Khedia, J., Agarwal, P., and Agarwal, P. K. (2018). AlNAC4 transcription factor from halophyte Aeluropus Lagopoides mitigates oxidative stress by maintaining ROS homeostasis in transgenic tobacco. Front. Plant Sci. 9:1522. doi: 10.3389/fpls.2018.01522

PubMed Abstract | CrossRef Full Text | Google Scholar

Kreps, J. A., Wu, Y., Chang, H. S., Zhu, T., and Harper, W. J. F. (2002). Transcriptome changes for Arabidopsis in response to salt, osmotic, and cold stress. Plant Physiol. 130, 2129–2141. doi: 10.2307/4280849

CrossRef Full Text | Google Scholar

Lei, Y., Xu, Y., Hettenhausen, C., Lu, C., Shen, G., Zhang, C., et al. (2018). Comparative analysis of alfalfa (Medicago sativa L.) leaf transcriptomes reveals genotype-specific salt tolerance mechanisms. BMC Plant Biol. 18:35. doi: 10.1186/s12870-018-1250-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Zhang, Y., Hu, X., Shen, X., Ma, L., Su, Z., et al. (2011). Transcriptional profiling of Medicago truncatula under salt stress identified a novel CBF transcription factor MtCBF4 that plays an important role in abiotic stress responses. BMC Plant Biol. 11:109. doi: 10.1186/1471-2229-11-109

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Han, G., Sun, C., and Sui, N. (2019). Research advances of MYB transcription factors in plant stress resistance and breeding. Plant Signal. Behav. 14:1613131. doi: 10.1080/15592324.2019.1613131

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Wang, N., Ji, D., Zhang, W., Wang, Y., Yu, Y., et al. (2019). A GmSIN1/GmNCED3s/GmRbohBs feed-forward loop acts as a signal amplifier that regulates root growth in soybean exposed to salt stress. Plant Cell 31, 2107–2130. doi: 10.1105/tpc.18.00662

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Wu, Z., Gu, H., Cheng, D., Guo, X., Li, L., et al. (2021). AvNAC030, a NAC domain transcription factor, enhances salt stress tolerance in Kiwifruit. Int. J. Mol. Sci. 22:11897. doi: 10.3390/ijms222111897

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., Xiong, C., Yan, L., Zhang, Z., Ma, L., Wang, Y., et al. (2017). Transcriptome analyses reveal candidate genes potentially involved in Al stress response in alfalfa. Front. Plant Sci. 8:26. doi: 10.3389/fpls.2017.00026

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X. P., and Yu, L. X. (2017). Genome-wide association mapping of loci associated with plant growth and forage production under salt stress in alfalfa (Medicago sativa L.). Front. Plant Sci. 8:853. doi: 10.3389/fpls.2017.00853

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z. Y., Baoyin, T., Li, X. L., and Wang, Z. L. (2019). How fall dormancy benefits alfalfa winter-survival? Physiologic and transcriptomic analyses of dormancy process. BMC Plant Biol. 19:205. doi: 10.1186/s12870-019-1773-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, R. C., Zhang, F., Zhang, Z., Li, M., Chen, L., Wang, X., et al. (2021). Assembly of chromosome-scale and allele-aware autotetraploid genome of the Chinese alfalfa cultivar Zhongmu-4 and identification of SNP loci associated with 27 agronomic traits. bioRxiv [Preprint]. doi: 10.1101/2021.02.21.428692

CrossRef Full Text | Google Scholar

Luo, M., Zhang, Y., Li, J., Zhang, P., Chen, K., Song, W., et al. (2021). Molecular dissection of maize seedling salt tolerance using a genome-wide association analysis method. Plant Biotechnol. J. 19, 1937–1951. doi: 10.1111/pbi.13607

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, X., Wang, B., Gao, S., Zhang, F., Terzaghi, W., and Dai, M. (2019). Genome-wide association study dissects the genetic bases of salt tolerance in maize seedlings. J. Integr. Plant Biol. 61, 658–674. doi: 10.1111/jipb.12797

PubMed Abstract | CrossRef Full Text | Google Scholar

Machado, R. M. A., and Serralheiro, R. P. (2017). Soil salinity: effect on vegetable crop growth. Management practices to prevent and mitigate soil salinization. Horticulturae 3:30. doi: 10.3390/horticulturae3020030

CrossRef Full Text | Google Scholar

Patishtan, J., Hartley, T. N., Fonseca de Carvalho, R., and Maathuis, F. J. M. (2018). Genome-wide association studies to identify rice salt-tolerance markers. Plant Cell Environ. 41, 970–982. doi: 10.1111/pce.12975

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, X., Li, M. W., Xie, M., Liu, X., Ni, M., Shao, G., et al. (2014). Identification of a novel salt tolerance gene in wild soybean by whole-genome sequencing. Nat. Commun. 5:4340. doi: 10.1038/ncomms5340

PubMed Abstract | CrossRef Full Text | Google Scholar

Rumbaugh, M. D., and Pendery, B. M. (1990). Germination salt resistance of alfalfa (Medicago sativa L.) germplasm in relation to subspecies and centers of diversity. Plant Soil 124, 47–51. doi: 10.1007/BF00010930

CrossRef Full Text | Google Scholar

Shahid, S. A., Zaman, M., and Heng, L. (2018). “Soil salinity: historical perspectives and a world overview of the problem,” in Guideline For Salinity Assessment, Mitigation and Adaptation Using Nuclear and Related Techniques, eds M. Zaman, S. A. Shahid, and L. Heng (Cham: Springer), 43–53. doi: 10.1007/978-3-319-96190-3_2

CrossRef Full Text | Google Scholar

Shen, C., Du, H., Chen, Z., Lu, H., Zhu, F., Chen, H., et al. (2020). The chromosome-level genome sequence of the autotetraploid alfalfa and resequencing of core germplasms provide genomic resources for alfalfa research. Mol. Plant 13, 1250–1261. doi: 10.1016/j.molp.2020.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Y., Gao, L., Wu, Z., Zhang, X., Wang, M., Zhang, C., et al. (2017). Genome-wide association study of salt tolerance at the seed germination stage in rice. BMC Plant Biol. 17:92. doi: 10.1186/s12870-017-1044-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, L., Jiang, L., Chen, Y., Shu, Y., Bai, Y., and Guo, C. (2016). Deep-sequencing transcriptome analysis of field-grown Medicago sativa L. crown buds acclimated to freezing stress. Funct. Integr. Genomics 16, 495–511. doi: 10.1007/s10142-016-0500-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Y., Bao, X., Zhi, Y., Wu, Q., Guo, Y., Yin, X., et al. (2019). Overexpression of a MYB family gene, OsMYB6, Increases drought and salinity stress tolerance in transgenic rice. Front. Plant Sci. 10:168. doi: 10.3389/fpls.2019.00168

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Verma, R. K., Venkata, V., Kumar, V. V. S., Yadav, S. K., and Chinnusamy, V. (2020). Overexpression of Arabidopsis ICE1 enhances yield and multiple abiotic stress tolerance in indica rice. Plant Signal Behav. 15:1814547. doi: 10.1080/15592324.2020.1814547

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J. W., Czech, B., and Weigel, D. (2009). miR156-regulated SPL transcription factors define an endogenous flowering pathway in Arabidopsis thaliana. Cell 138, 738–749. doi: 10.1016/j.cell.2009.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R. K., Cao, Z. H., and Hao, Y. J. (2014). Overexpression of a R2R3 MYB gene MdSIMYB1 increases tolerance to multiple stresses in transgenic tobacco and apples. Physiol. Plant 150, 76–87. doi: 10.1111/ppl.12069

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Vinocur, B., and Altman, A. (2003). Plant responses to drought, salinity and extreme temperatures: towards genetic engineering for stress tolerance. Planta 218, 1–14. doi: 10.1007/s00425-003-1105-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Wang, X., Zhang, H., Ma, L., Zhao, H., Jones, C. S., et al. (2020). A genome-wide association study approach to the identification of candidate genes underlying agronomic traits in alfalfa (Medicago sativa L.). Plant Biotechnol. J. 18, 611–613. doi: 10.1111/pbi.13251

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, H., Zhao, Y., Xie, Y., and Wang, H. (2018). Exploiting SPL genes to improve maize plant architecture tailored for high-density planting. J. Exp. Bot. 69, 4675–4688. doi: 10.1093/jxb/ery258

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J., Jiang, Y., Liang, Y., Chen, L., Chen, W., and Cheng, B. (2019). Expression of the maize MYB transcription factor ZmMYB3R enhances drought and salt stress tolerance in transgenic plants. Plant Physiol. Biochem. 137, 179–188. doi: 10.1016/j.plaphy.2019.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, G., Cao, J., Wang, X., Chen, Q., and Tian, F. (2019). Evolutionary metabolomics identifies substantial metabolic divergence between maize and its wild ancestor, teosinte. Plant Cell 31, 1990–2019. doi: 10.1105/tpc.19.00111

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, R., and Sun, C. (2021). What happened during domestication of wild to cultivated rice. Crop J. 9, 564–576. doi: 10.1016/j.cj.2021.02.005

CrossRef Full Text | Google Scholar

Yao, L., Li, Y., Ma, C., Tong, L., Du, F., and Xu, M. (2020). Combined genome-wide association study and transcriptome analysis reveal candidate genes for resistance to Fusarium ear rot in maize. J. Integr. Plant Biol. 62, 1535–1551. doi: 10.1111/jipb.12911

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, S., Han, Y., Huang, L., Hong, Y., and Zhang, G. (2018). Overexpression of HvCBF7 and HvCBF9 changes salt and drought tolerance in Arabidopsis. Plant Growth Regul. 85, 281–292. doi: 10.1007/s10725-018-0394-4

CrossRef Full Text | Google Scholar

Yu, J., Zao, W., He, Q., Kim, T. S., and Park, Y. J. (2017). Genome-wide association study and gene set analysis for understanding candidate genes involved in salt tolerance at the rice seedling stage. Mol. Genet. Genomics 292, 1391–1403. doi: 10.1007/s00438-017-1354-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, L. X., Liu, X., Boge, W., and Liu, X. P. (2016). Genome-wide association study identifies loci for salt tolerance during germination in autotetraploid Alfalfa (Medicago sativa L.) using genotyping-by-sequencing. Front. Plant Sci. 7:956. doi: 10.3389/fpls.2016.00956

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, L. X., Zhang, F., Culma, C. M., Lin, S., Niu, Y., Zhang, T., et al. (2020). Construction of high-density linkage maps and identification of quantitative trait loci associated with verticillium wilt resistance in autotetraploid alfalfa (Medicago sativa L.). Plant Dis. 104, 1439–1444. doi: 10.1094/PDIS-08-19-1718-RE

PubMed Abstract | CrossRef Full Text | Google Scholar

Zelm, E. V., Zhang, Y., and Testerink, C. (2020). Salt tolerance mechanisms of plants. Annu. Rev. Plant Biol. 71, 403–433. doi: 10.1146/annurev-arplant-050718-100005

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, A., Chen, P., Korth, K., Hancock, F., Pereira, A., Brye, K., et al. (2017). Genome-wide association study (GWAS) of salt tolerance in worldwide soybean germplasm lines. Mol. Breed. 37:30. doi: 10.1007/s11032-017-0634-8

CrossRef Full Text | Google Scholar

Zhang, A., Yang, X., Lu, J., Song, F., Sun, J., Wang, C., et al. (2021). OsIAA20, an Aux/IAA protein, mediates abiotic stress tolerance in rice through an ABA pathway. Plant Sci. 308:110903. doi: 10.1016/j.plantsci.2021.110903

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Long, Y., and Chen, X. (2021). A NAC transcription factor OsNAC3 positively regulates ABA response and salt tolerance in rice. BMC Plant Biol. 21:546. doi: 10.1186/s12870-021-03333-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Irving, L. J., Mcgill, C., Matthew, C., and Zhou, D. (2010). The effects of salinity and osmotic stress on barley germination rate: sodium as an osmotic regulator. Ann. Bot. 106, 1027–1035. doi: 10.1093/aob/mcq204

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, M., Liu, Y., Cai, H., Guo, M., Chai, M., She, Z., et al. (2020). The bZIP transcription factor GmbZIP15 negatively regulates salt- and drought-stress responses in soybean. Int. J. Mol. Sci. 21:7778. doi: 10.3390/ijms21207778

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, P., Wang, R., Yang, X., Ju, Q., Li, W., Lü, S., et al. (2020). The R2R3-MYB transcription factor AtMYB49 modulates salt tolerance in Arabidopsis by modulating the cuticle formation and antioxidant defence. Plant Cell Environ. 43, 1925–1943. doi: 10.1111/pce.13784

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Zhang, H., Song, C., Zhu, J. K., and Shabala, S. (2020). Mechanisms of plant responses and adaptation to soil salinity. Innovation 1:41. doi: 10.1016/j.xinn.2020.100017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Zhang, Z., Xie, S., Si, T., Li, Y., and Zhu, J. K. (2016). Mutational evidence for the critical role of CBF transcription factors in cold acclimation in Arabidopsis. Plant Physiol. 171, 2744–2759. doi: 10.1104/pp.16.00533

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Yang, Z., Ding, Y., Liu, L., and Ge, X. (2019). Over-expression of an R2R3 MYB gene, GhMYB73, increases tolerance to salt stress in transgenic Arabidopsis. Plant Sci. 286, 28–36. doi: 10.1016/j.plantsci.2019.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: alfalfa, GWAS, salt stress, RNA-seq, auxin

Citation: He F, Wei C, Zhang Y, Long R, Li M, Wang Z, Yang Q, Kang J and Chen L (2022) Genome-Wide Association Analysis Coupled With Transcriptome Analysis Reveals Candidate Genes Related to Salt Stress in Alfalfa (Medicago sativa L.). Front. Plant Sci. 12:826584. doi: 10.3389/fpls.2021.826584

Received: 01 December 2021; Accepted: 28 December 2021;
Published: 03 February 2022.

Edited by:

Baohua Wang, Nantong University, China

Reviewed by:

Sibu Simon, Mar Athanasios College for Advanced Studies Tiruvalla (MACFAST), India
Cesar Augusto Medina, Agricultural Research Service, United States Department of Agriculture (USDA), United States

Copyright © 2022 He, Wei, Zhang, Long, Li, Wang, Yang, Kang and Chen. 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: Lin Chen, chenlin@caas.cn; Junmei Kang, kangjmei@126.com

Disclaimer: 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.