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

Transcriptome Analysis for Salt-Responsive Genes in Two Different Alfalfa (Medicago sativa L.) Cultivars and Functional Analysis of MsHPCA1

Alfalfa (Medicago sativa L.) is an important forage legume and soil salinization seriously affects its growth and yield. In a previous study, we identified a salt-tolerant variety 'Gongnong NO.1' and a salt-sensitive variety 'Sibeide'. To unravel the molecular mechanism involved in salt stress, we conducted transcriptomic analysis on these two cultivars grown under 0 and 250 mM NaCl treatments for 0, 12, and 24 h. Totals of 336, and 548 differentially expressed genes (DEGs) in response to NaCl were, respectively, identified in the 'Gongnong NO.1' and 'Sibeide' varieties. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway enrichment analysis showed that the DEGs were classified in carbohydrate metabolism, energy production, transcription factor, and stress-associated pathway. Expression of MsHPCA1, encoding a putative H2O2 receptor, was responsive to both NaCl and H2O2 treatment. MsHPCA1 was localized in cell membrane and overexpression of MsHPCA1 in alfalfa increased salt tolerance and H2O2 content. This study will provide new gene resources for the improvement in salt tolerance in alfalfa and legume crops, which has important theoretical significance and potential application value.

MsSPL12 is a positive regulator in alfalfa (Medicago sativa L.) salt tolerance

KEY MESSAGE

Over expression of MsSPL12 improved alfalfa salt tolerance by reducing Na+ accumulation and increasing antioxidant enzyme activity and regulating down-stream gene expression. Improvement of salt tolerance is one of the major goals in alfalfa breeding. Here, we demonstrated that MsSPL12, an alfalfa transcription factor gene highly expressed in the stem cells, plays a positive role in alfalfa salt tolerance. MsSPL12 is localized in the nucleus and shows transcriptional activity in the presence of its C-terminus. To investigate MsSPL12 function in plant response to salt stress, we generated transgenic plants overexpressing either MsSPL12 or a chimeric MsSPL12-SRDX gene that represses the function of MsSPL12 by using the Chimeric REpressor gene-Silencing Technology (CRES-T), and observed that overexpression of MsSPL12 increased the salt tolerance of alfalfa transgenic plants associated with an increase in K+/Na+ ratio and relative water content (RWC) under salt stress treatment, but a reduction in electrolyte leakage (EL), reactive oxygen species (ROS), malondialdehyde (MDA), and proline (Pro) compared to wild type (WT) plants. However, transgenic plants overexpressing MsSPL12-SRDX showed an inhibited plant growth and a reduced salt tolerance. RNA-sequencing and quantitative real-time PCR analyses revealed that MsSPL12 affected the expression of plant abiotic resistance-related genes in multiple physiological pathways. The potential MsSPL12-mediated regulatory pathways based on the differentially expressed genes between the MsSPL12 overexpression transgenics and WT controls were predicted. In summary, our study proves that MsSPL12 is a positive regulator in alfalfa salt tolerance and can be used as a new candidate for manipulation to develop forage crops with enhanced salt tolerance.

Genome-wide identification and expression analysis of the glycosyl hydrolase family 1 genes in Medicago sativa revealed their potential roles in response to multiple abiotic stresses

Glycoside hydrolase family 1 (GH1) β-glucosidases (BGLUs), are encoded by a large number of genes, which participate in the development and stress response of plants, particularly under biotic and abiotic stresses through the activation of phytohormones. However, there are few studies systematically analyzing stress or hormone-responsive BGLU genes in alfalfa. In this study, a total of 179 BGLU genes of the glycoside hydrolase family 1 were identified in the genome of alfalfa, and then were classified into five distinct clusters. Sequence alignments revealed several conserved and unique motifs among these MsBGLU proteins. Many cis-acting elements related to abiotic stresses and phytohormones were identified in the promoter of some MsBGLUs. Moreover, RNA-seq and RT-qPCR analyses showed that these MsBGLU genes exhibited distinct expression patterns in response to different abiotic stress and hormonal treatments. In summary, this study suggests that MsBGLU genes play crucial roles in response to various abiotic stresses and hormonal responses, and provides candidate genes for stress tolerance breeding in alfalfa.

Genome-Wide Association Study Unravels Quantitative Trait Loci and Genes Associated with Yield-Related Traits in Sugarcane

Sugarcane, a major sugar and energy crop worldwide faces an increasing demand for higher yields. Identifying yield-related markers and candidate genes is valuable for breeding high-yield varieties using molecular techniques. In this work, seven yield-related traits were evaluated in a diversity panel of 159 genotypes, derived from Tripidium arundinaceum, Saccharum spontaneum, and modern sugarcane genotypes. All traits exhibited significant genetic variance with high heritability and high correlations. Genetic diversity analysis reveals a genomic decay of 23 kb and an average single nucleotide polymorphism (SNP) number of 25,429 per genotype. These 159 genotypes were divided into 4 subgroups. Genome-wide association analysis identified 47 SNPs associated with brix, spanning 36 quantitative trait loci (QTLs), and 138 SNPs for other traits across 104 QTLs, covering all 32 chromosomes. Interestingly, 12 stable QTLs associated with yield-related traits were identified, which contained 35 candidate genes. This work provides markers and candidate genes for marker-assisted breeding to improve sugarcane yields.

Association Analysis of Tiller-Related Traits with EST-SSR Markers in Psathyrostachys juncea

Psathyrostachys juncea is a long-lived perennial Gramineae grass with dense basal tillers and soft leaves. It is used widely in cold and dry areas of Eurasia and North America to establish grazing pasture and is even used as an ideal plant for revegetation and ecological restoration. Plant architecture, especially tillering traits, is critical for bunch grasses in breeding programs, and these traits in plants are mostly quantitative traits. In this study, the genetic diversity, population structure, and linkage disequilibrium of 480 individual lines were analyzed using 127 pairs of the EST-SSR marker, and a significant association between ten plant-architecture-related traits of P. juncea and molecular markers was found. The results of the genetic diversity analysis showed that the number of observed alleles was 1.957, the number of effective alleles was 1.682, Shannon's information index was 0.554, observed heterozygosity was 0.353, expected heterozygosity was 0.379, and the polymorphism information content was 0.300. A total of 480 individual lines were clustered into five groups based on population genetic structure, principal coordinate analysis, and unweighted pair group method with arithmetic mean analysis (UPGMA). The linkage disequilibrium coefficient (r2) was between 0.00 and 0.68, with an average of 0.04, which indicated a relatively low level of linkage disequilibrium among loci. The results of the association analysis revealed 55 significant marker-trait associations (MTA). Moreover, nine SSR markers were associated with multiple traits. This study provides tools with promising applications in the molecular selection and breeding of P. juncea germplasm.

Graphene nanoparticles improve alfalfa (Medicago sativa L.) growth through multiple metabolic pathways under salinity-stressed environment

Graphene, one of the emerging carbon nanomaterials, has many advantages and applications. Salinity stress seriously affects ecology and agroforestry worldwide. The effects of graphene on alfalfa under salinity stress were investigated. The results indicated that graphene promoted alfalfa growth under non-salinity stress but caused some degree of damage to root cells and leaf parameters. Graphene used in salinity stress had a positive effect on growth parameters, chlorophyll, photosynthetic gas parameters, stomatal opening, ion balance, osmotic homeostasis, cell membrane integrity and antioxidant system, while it decreased Na+, lipid peroxidation and reactive oxygen species levels. Correlation analysis revealed that most of the parameters were significantly correlated; and principal component analysis indicated that the first two dimensions (78.1% and 4.1%) explained 82.2% of the total variability, and the majority of them exceeded the average contribution. Additionally, Gene Ontology functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes signaling pathway enrichment analysis showed that there were numerous differentially expressed genes and pathways to regulate alfalfa responding to salinity stress. Taken together, the findings reveal that graphene does not enter the plant, but improves the properties and adsorption of soil to enhance salt tolerance and seedling growth of alfalfa through morphological, physiological, biochemical, and transcriptomic aspects. Furthermore, this study provides a reference for the application of graphene to improve soil environment and agricultural production.

Ideotype breeding and genome engineering for legume crop improvement

Ideotype breeding is a strategy whereby traits are modeled a priori and then introduced into a model or crop species to assess their impact on yield. Thus, knowledge about the connection between genotype and phenotype is required for ideotype breeding to be deployed successfully. The growing understanding of the genetic basis of yield-related traits, combined with increasingly efficient genome engineering tools, improved transformation efficiency, and high-throughput genotyping of regenerants paves the way for the widespread adoption of ideotype breeding as a complement to conventional breeding. We briefly discuss how ideotype breeding, coupled with such state-of-the-art biotechnological tools, could contribute to knowledge-based legume breeding and accelerate yield gains to ensure food security in the coming decades.

Comprehensive analysis of NAC transcription factors in Scutellaria baicalensis and their response to exogenous ABA and GA3

The NAC is a plant-specific family of transcription factor that plays important roles in various biological processes. Scutellaria baicalensis Georgi, belongs to the Lamiaceae family and has been widely used as a traditional herb with a wide range of pharmacological activities, including antitumor, heat-clearing, and detoxifying functions. However, no study on the NAC family in S. baicalensis has been conducted to date. In the present study, we identified 56 SbNAC genes using genomic and transcriptome analyses. These 56 SbNACs were unevenly distributed across nine chromosomes and were phylogenetically divided into six clusters. Cis-element analysis identified plant growth and development-, phytohormone-, light-, and stress-responsive elements were present in SbNAC genes promoter regions. Protein-protein interaction analysis was performed using Arabidopsis homologous proteins. Potential transcription factors, including bHLH, ERF, MYB, WRKY, and bZIP, were identified and constructed a regulatory network with SbNAC genes. The expression of 12 flavonoid biosynthetic genes was significantly upregulated with abscisic acid (ABA) and gibberellin (GA₃) treatments. Eight SbNAC genes (SbNAC9/32/33/40/42/43/48/50) also exhibited notable variation with two phytohormone treatments, among which SbNAC9 and SbNAC43 showed the most significant variation and deserved further study. Additionally, SbNAC44 displayed a positive correlation with C4H3, PAL5, OMT3, and OMT6, while SbNAC25 had negatively correlated with OMT2, CHI, F6H2, and FNSII-2. This study constitutes the first analysis of SbNAC genes and lays the basis foundation for further functional studies of SbNAC genes family members, while it may also facilitate the genetic improvement of plants and breeding of elite S. baicalensis varieties.

Role of bZIP Transcription Factors in Plant Salt Stress

Soil salinity has become an increasingly serious problem worldwide, greatly limiting crop development and yield, and posing a major challenge to plant breeding. Basic leucine zipper (bZIP) transcription factors are the most widely distributed and conserved transcription factors and are the main regulators controlling various plant response processes against external stimuli. The bZIP protein contains two domains: a highly conserved, DNA-binding alkaline region, and a diverse leucine zipper, which is one of the largest transcription factor families in plants. Plant bZIP is involved in many biological processes, such as flower development, seed maturation, dormancy, and senescence, and plays an important role in abiotic stresses such as salt damage, drought, cold damage, osmotic stress, mechanical damage, and ABA signal response. In addition, bZIP is involved in the regulation of plant response to biological stresses such as insect pests and pathogen infection through salicylic acid, jasmonic acid, and ABA signal transduction pathways. This review summarizes and discusses the structural characteristics and functional characterization of the bZIP transcription factor group, the bZIP transcription factor complex and its molecular regulation mechanisms related to salt stress resistance, and the regulation of transcription factors in plant salt stress resistance. This review provides a theoretical basis and research ideas for further exploration of the salt stress-related functions of bZIP transcription factors. It also provides a theoretical basis for crop genetic improvement and green production in agriculture.

Transcriptome and GWAS Analyses Reveal Candidate Gene for Root Traits of Alfalfa during Germination under Salt Stress

Alfalfa growth and production in China are negatively impacted by high salt concentrations in soils, especially in regions with limited water supplies. Few reliable genetic markers are currently available for salt tolerance selection. As a result, molecular breeding strategies targeting alfalfa are hindered. Therefore, with the continuous increase in soil salinity in agricultural lands, it is indispensable that a salt-tolerant variety of alfalfa is produced. We collected 220 alfalfa varieties around the world for resequencing and performed genome-wide association studies (GWASs). Alfalfa seeds were germinated in saline water with different concentrations of NaCl, and the phenotypic differences in several key root traits were recorded. In the phenotypic analysis, the breeding status and geographical origin strongly affected the salt tolerance of alfalfa. Forty-nine markers were significantly associated with salt tolerance, and 103 candidate genes were identified based on linkage disequilibrium. A total of 2712 differentially expressed genes were upregulated and 3570 were downregulated based on transcriptomic analyses. Some candidate genes that affected root development in the seed germination stage were identified through the combination of GWASs and transcriptome analyses. These genes could be used for molecular breeding strategies to increase alfalfa’s salt tolerance and for further research on salt tolerance in general.

Genome-Wide Identification and Phylogenetic and Expression Analyses of the PLATZ Gene Family in Medicago sativa L

The PLATZ family is a novel class of plant-specific zinc finger transcription factors with important roles in plant growth and development and abiotic stress responses. PLATZ members have been identified in many plants, including Oryza sativa, Zea mays, Triticum aestivum, Fagopyrum tataricum, and Arabidopsis thaliana; however, due to the complexity of the alfalfa reference genome, the members of the PLATZ gene family in alfalfa (Medicago sativa L.) have not been systematically identified and analyzed. In this study, 55 Medicago sativa PLATZ genes (MsPLATZs) were identified in the alfalfa "Xinjiangdaye" reference genome. Basic bioinformatic analysis was performed, including the characterization of sequence lengths, protein molecular weights, genomic positions, and conserved motifs. Expression analysis reveals that 7 MsPLATZs are tissue-specifically expressed, and 10 MsPLATZs are expressed in all examined tissues. The transcriptomic expression of these genes is obvious, indicating that these MsPLATZs have different functions in the growth and development of alfalfa. Based on transcriptome data analysis and real-time quantitative PCR (RT-qPCR), we identified 22, 22, and 21 MsPLATZ genes that responded to salt, cold, and drought stress, respectively, with 20 MsPLATZs responding to all three stresses. This study lays a foundation for further exploring the functions of MsPLATZs, and provides ideas for the improvement of alfalfa varieties and germplasm innovation.

Full-length transcriptome of in Medicago sativa L. roots in response to drought stress

Background: Alfalfa (Medicago sativa L.), serves as a legume with high drought tolerance, is a major forage crop with a high biomass of production. However, the molecular mechanism of Alfalfa in response to drought stress are still unclear. Results: We constructed the first full-length transcriptome for Alfalfa root. 21.53Gb clean data were obtained by further data filtering, in which incorporate 566,076 reads of Insert (ROI), and 409,291 full length reads non-Chimeric (FLNC) sequences. Combined with second-generation sequencing (SGS), there were 2615, 6011, and 4617 differentially expressed genes (DEGs) in three comparisons. KEGG pathway analysis showed enrichment of ribosome, glutathione metabolism, and biosynthesis of amino acids are among the DEGs. The majority of transcription factors (TFs) from DEGs were AP2/ERF-ERF (37), C2H2 (32), and bHLH (22) bZIP (22), followed by C3H (19), MYB (18), WRKY (18), GRAS (16), and NAC (15). 32 C2H2 genes were differentially expressed in three groups. In addition, TFs annotated as C3H (19), MYB (18), GRAS (16), and NAC (15) also changed significantly in expression in the three comparisons. We found 24 genes participate in the abscisic acid (ABA) and auxin hormone signaling pathway in response to drought stress, and monitored the expression patterns of these related genes. Conclusion: The present study enhanced our understanding of the genetic diversity and complexity, and provides greater insight into the fundamental transcriptome reprogramming of Alfalfa under drought.

Genome-wide association scan and transcriptome analysis reveal candidate genes for waterlogging tolerance in cultivated barley

Waterlogging is the primary abiotic factor that destabilizes the yield and quality of barley (Hordeum vulgare L.). However, the genetic basis of waterlogging tolerance remains poorly understood. In this study, we conducted a genome-wide association study (GWAS) by involving 106,131 single-nucleotide polymorphisms (SNPs) with a waterlogging score (WLS) of 250 barley accessions in two years. Out of 72 SNPs that were found to be associated with WLS, 34 were detected in at least two environments. We further performed the transcriptome analysis in root samples from TX9425 (waterlogging tolerant) and Franklin (waterlogging sensitive), resulting in the identification of 5,693 and 8,462 differentially expressed genes (DEGs) in these genotypes, respectively. The identified DEGs included various transcription factor (TF) genes, primarily including AP2/ERF, bZIP and MYB. By combining GWAS and RNA-seq, we identified 27 candidate genes associated with waterlogging, of which three TFs (HvDnaJ, HvMADS and HvERF1) were detected in multiple treatments. Moreover, by overexpressing barley HvERF1 in Arabidopsis, the transgenic lines were detected with enhanced waterlogging tolerance. Altogether, our results provide new insights into the genetic mechanisms of waterlogging, which have implications in the molecular breeding of waterlogging-tolerant barley varieties.

Genome-Wide Identification and Expression Analysis of HSF Transcription Factors in Alfalfa (Medicago sativa) under Abiotic Stress

Alfalfa (Medicago sativa) is one of the most important legume forage species in the world. It is often affected by several abiotic stressors that result in reduced yields and poor growth. Therefore, it is crucial to study the resistance of M. sativa to abiotic stresses. Heat shock transcription factors (HSF) are key players in a number of transcriptional regulatory pathways. These pathways play an essential role in controlling how plants react to different abiotic stressors. Studies on the HSF gene family have been reported in many species but have not yet undergone a thorough analysis in M. sativa. Therefore, in order to identify a more comprehensive set of HSF genes, from the genomic data, we identified 16 members of the MsHSF gene, which were unevenly distributed over six chromosomes. We also looked at their gene architectures and protein motifs, and phylogenetic analysis allowed us to divide them into 3 groups with a total of 15 subgroups. Along with these aspects, we then examined the physicochemical properties, subcellular localization, synteny analysis, GO annotation and enrichment, and protein interaction networks of amino acids. Finally, the analysis of 16 MsHSF genes' expression levels across all tissues and under four abiotic stresses using publicly available RNA-Seq data revealed that these genes had significant tissue-specific expression. Moreover, the expression of most MsHSF genes increased dramatically under abiotic stress, further validating the critical function played by the MsHSF gene family in abiotic stress. These results provided basic information about MsHSF gene family and laid a foundation for further study on the biological role of MsHSF gene in response to stress in M. sativa.

Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Alfalfa Revealed Its Potential Roles in Response to Multiple Abiotic Stresses

NAC (NAM, ATAF1/2, and CUC2) transcription factors compose one of the largest families of plant-specific transcription factors; they are widely involved in plant growth and development and have especially important roles in improving stress resistance in plants. However, NAC gene family members in alfalfa (Medicago sativa L.) have not been systematically identified and analyzed genome-wide due to the complexity of the alfalfa reference genome. In this study, a total of 421 M. sativa NAC genes (MsNACs) were identified from the alfalfa “Xinjiangdaye” reference genome. Basic bioinformatics analysis, including characterization of sequence length, protein molecular weight and genome position and conserved motif analysis, was conducted. Expression analysis showed that 47 MsNACs had tissue-specific expression, and 64 MsNACs were expressed in all tissues. The transcriptomic profiles of the genes were very different, indicating that these MsNACs have various functions in alfalfa growth and development. We identified 25, 42 and 47 MsNACs that respond to cold, drought and salt stress based on transcriptome data analysis and real-time quantitative PCR (RT−qPCR). Furthermore, 22 MsNACs were found to respond to both salt and drought stress, and 15 MsNACs were found to respond to cold, salt and drought stress. The results of this study could provide valuable information for further functional analysis of MsNACs and for the improvement of stress resistance in alfalfa.

Genome-Wide Identification and Analysis of the NF-Y Transcription Factor Family Reveal Its Potential Roles in Salt Stress in Alfalfa (Medicago sativa L.)

Nuclear factor Y (NF-Y) is a heterotrimeric transcription factor that plays an important role in various biological processes in plants, such as flowering regulation, drought resistance, and salt stress. However, few in-depth studies investigated the alfalfa NF-Y gene family. In this study, in total, 60 MsNF-Y genes, including 9 MsNF-YAs, 26 MsNF-YBs, and 25 MsNF-YCs, were identified in the alfalfa genome. The genomic locations, gene structures, protein molecular weights, conserved domains, phylogenetic relationships, and gene expression patterns in different tissues and under different stresses (cold stress, drought stress, and salt stress) of these NF-Y genes were analyzed. The illustration of the conserved domains and specific domains of the different subfamilies of the MsNF-Y genes implicates the conservation and diversity of their functions in alfalfa growth, development, and stress resistance. The gene expression analysis showed that 48 MsNF-Y genes (7 MsNF-YAs, 22 MsNF-YBs, and 19 MsNF-YCs) were expressed in all tissues at different expression levels, indicating that these genes have tissue expression specificity and different biological functions. In total, seven, seven, six, and eight MsNF-Y genes responded to cold stress, the ABA treatment, drought stress, and salt stress in alfalfa, respectively. According to the WGCNA, molecular regulatory networks related to salt stress were constructed for MsNF-YB2, MsNF-YB5, MsNF-YB7, MsNF-YB15, MsNF-YC5, and MsNF-YC6. This study could provide valuable information for further elucidating the biological functions of MsNF-Ys and improving salt tolerance and other abiotic stress resistance in alfalfa.

Combined GWAS and Transcriptome Analyses Provide New Insights Into the Response Mechanisms of Sunflower Against Drought Stress

Sunflower is one of the most important oil crops in the world, and drought stress can severely limit its production and quality. To understand the underlying mechanism of drought tolerance, and identify candidate genes for drought tolerance breeding, we conducted a combined genome-wide association studies (GWAS) and RNA-seq analysis. A total of 226 sunflower inbred lines were collected from different regions of China and other countries. Eight phenotypic traits were evaluated under control and drought stress conditions. Genotyping was performed using a Specific-Locus Amplified Fragment Sequencing (SLAF-seq) approach. A total of 934.08 M paired-end reads were generated, with an average Q30 of 91.97%. Based on the 243,291 polymorphic SLAF tags, a total of 94,162 high-quality SNPs were identified. Subsequent analysis of linkage disequilibrium (LD) and population structure in the 226 accessions was carried out based on the 94,162 high-quality SNPs. The average LD decay across the genome was 20 kb. Admixture analysis indicated that the entire population most likely originated from 11 ancestors. GWAS was performed using three methods (MLM, FarmCPU, and BLINK) simultaneously. A total of 80 SNPs showed significant associations with the 8 traits (p < 1.062 × 10-6). Next, a total of 118 candidate genes were found. To obtain more reliable candidate genes, RNA-seq analysis was subsequently performed. An inbred line with the highest drought tolerance was selected according to phenotypic traits. RNA was extracted from leaves at 0, 7, and 14 days of drought treatment. A total of 18,922 differentially expressed genes were obtained. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis showed up-regulated genes were mainly enriched in the branched-chain amino acid catabolic process, while the down-regulated genes were mainly enriched in the photosynthesis-related process. Six DEGs were randomly selected from all DEGs for validation; these genes showed similar patterns in RNA-seq and RT-qPCR analysis, with a correlation coefficient of 0.8167. Through the integration of the genome-wide association study and the RNA-sequencing, 14 candidate genes were identified. Four of them (LOC110885273, LOC110872899, LOC110891369, LOC110920644) were abscisic acid related protein kinases and transcription factors. These genes may play an important role in sunflower drought response and will be used for further study. Our findings provide new insights into the response mechanisms of sunflowers against drought stress and contribute to further genetic breeding.