Genome-Wide Association Study Identifies Loci for Salt Tolerance during Germination in Autotetraploid Alfalfa (Medicago sativa L.) Using Genotyping-by-Sequencing

Discovering new genes for alfalfa (Medicago sativa) growth and biomass resilience in combined salinity and Phoma medicaginis infection through GWAS

Salinity and Phoma medicaginis infection represent significant challenges for alfalfa cultivation in South Africa, Europe, Australia, and, particularly, Tunisia. These constraints have a severe impact on both yield and quality. The primary aim of this study was to establish the genetic basis of traits associated with biomass and growth of 129 Medicago sativa genotypes through genome-wide association studies (GWAS) under combined salt and P. medicaginis infection stresses. The results of the analysis of variance (ANOVA) indicated that the variation in these traits could be primarily attributed to genotype effects. Among the test genotypes, the length of the main stem, the number of ramifications, the number of chlorotic leaves, and the aerial fresh weight exhibited the most significant variation. The broad-sense heritability (H²) was relatively high for most of the assessed traits, primarily due to genetic factors. Cluster analysis, applied to morpho-physiological traits under the combined stresses, revealed three major groups of accessions. Subsequently, a GWAS analysis was conducted to validate significant associations between 54,866 SNP-filtered single-nucleotide polymorphisms (SNPs) and seven traits. The study identified 27 SNPs that were significantly associated with the following traits: number of healthy leaves (two SNPs), number of chlorotic leaves (five SNPs), number of infected necrotic leaves (three SNPs), aerial fresh weight (six SNPs), aerial dry weight (nine SNPs), number of ramifications (one SNP), and length of the main stem (one SNP). Some of these markers are related to the ionic transporters, cell membrane rigidity (related to salinity tolerance), and the NBS_LRR gene family (associated with disease resistance). These findings underscore the potential for selecting alfalfa genotypes with tolerance to the combined constraints of salinity and P. medicaginis infection.

Identification of Single Nucleotide Polymorphic Loci and Candidate Genes for Seed Germination Percentage in Okra under Salt and No-Salt Stresses by Genome-Wide Association Study

Excessive soil salinity is a major stressor inhibiting crops' growth, development, and yield. Seed germination is a critical stage of crop growth and development, as well as one of the most salt-sensitive stages. Salt stress has a significant inhibitory effect on seed germination. Okra is a nutritious vegetable, but its seed germination percentage (GP) is low, whether under salt stress conditions or suitable conditions. In this study, we used 180 okra accessions and conducted a genome-wide association study (GWAS) on the germination percentage using 20,133,859 single nucleotide polymorphic (SNP) markers under 0 (CK, diluted water), 70 (treatment 1, T1), and 140 mmol/L (treatment 2, T2) NaCl conditions. Using the mixed linear model (MLM) in Efficient Mixed-model Association eXpedated (EMMAX) and Genome-wide Efficient Mixed Model Association (GEMMA) software, 511 SNP loci were significantly associated during germination, of which 167 SNP loci were detected simultaneously by both programs. Among the 167 SNPs, SNP2619493 on chromosome 59 and SNP2692266 on chromosome 44 were detected simultaneously under the CK, T1, and T2 conditions, and were key SNP loci regulating the GP of okra seeds. Linkage disequilibrium block analysis revealed that nsSNP2626294 (C/T) in Ae59G004900 was near SNP2619493, and the amino acid changes caused by nsSNP2626294 led to an increase in the phenotypic values in some okra accessions. There was an nsSNP2688406 (A/G) in Ae44G005470 near SNP2692266, and the amino acid change caused by nsSNP2688406 led to a decrease in phenotypic values in some okra accessions. These results indicate that Ae59G004900 and Ae44G005470 regulate the GP of okra seeds under salt and no-salt stresses. The gene expression analysis further demonstrated these results. The SNP markers and genes that were identified in this study will provide reference for further research on the GP of okra, as well as new genetic markers and candidate genes for cultivating new okra varieties with high GPs under salt and no-salt stress conditions.

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.

Comparative time-course transcriptome analysis of two contrasting alfalfa (Medicago sativa L.) genotypes reveals tolerance mechanisms to salt stress

Salt stress is a major abiotic stress affecting plant growth and crop yield. For the successful cultivation of alfalfa (Medicago sativa L.), a key legume forage, in saline-affected areas, it's essential to explore genetic modifications to improve salt-tolerance.Transcriptome assay of two comparative alfalfa genotypes, Adina and Zhaodong, following a 4 h and 8 h's 300 mM NaCl treatment was conducted in this study in order to investigate the molecular mechanism in alfalfa under salt stress conditions. Results showed that we obtained 875,023,571 transcripts and 662,765,594 unigenes were abtained from the sequenced libraries, and 520,091 assembled unigenes were annotated in at least one database. Among them, we identified 1,636 differentially expression genes (DEGs) in Adina, of which 1,426 were up-regulated and 210 down-regulated, and 1,295 DEGs in Zhaodong, of which 565 were up-regulated and 730 down-regulated. GO annotations and KEGG pathway enrichments of the DEGs based on RNA-seq data indicated that DEGs were involved in (1) ion and membrane homeostasis, including ABC transporter, CLC, NCX, and NHX; (2) Ca²⁺ sensing and transduction, including BK channel, EF-hand domain, and calmodulin binding protein; (3) phytohormone signaling and regulation, including TPR, FBP, LRR, and PP2C; (4) transcription factors, including zinc finger proteins, YABBY, and SBP-box; (5) antioxidation process, including GST, PYROX, and ALDH; (6) post-translational modification, including UCH, ubiquitin family, GT, MT and SOT. The functional roles of DEGs could explain the variations in salt tolerance performance observed between the two alfalfa genotypes Adina and Zhaodong. Our study widens the understanding of the sophisticated molecular response and tolerance mechanism to salt stress, providing novel insights on candidate genes and pathways for genetic modification involved in salt stress adaptation in alfalfa.

Application of machine learning to explore the genomic prediction accuracy of fall dormancy in autotetraploid alfalfa

Fall dormancy (FD) is an essential trait to overcome winter damage and for alfalfa (Medicago sativa) cultivar selection. The plant regrowth height after autumn clipping is an indirect way to evaluate FD. Transcriptomics, proteomics, and quantitative trait locus mapping have revealed crucial genes correlated with FD; however, these genes cannot predict alfalfa FD very well. Here, we conducted genomic prediction of FD using whole-genome SNP markers based on machine learning-related methods, including support vector machine (SVM) regression, and regularization-related methods, such as Lasso and ridge regression. The results showed that using SVM regression with linear kernel and the top 3000 genome-wide association study (GWAS)-associated markers achieved the highest prediction accuracy for FD of 64.1%. For plant regrowth height, the prediction accuracy was 59.0% using the 3000 GWAS-associated markers and the SVM linear model. This was better than the results using whole-genome markers (25.0%). Therefore, the method we explored for alfalfa FD prediction outperformed the other models, such as Lasso and ElasticNet. The study suggests the feasibility of using machine learning to predict FD with GWAS-associated markers, and the GWAS-associated markers combined with machine learning would benefit FD-related traits as well. Application of the methodology may provide potential targets for FD selection, which would accelerate genetic research and molecular breeding of alfalfa with optimized FD.

A Genome-Wide Association Study Coupled With a Transcriptomic Analysis Reveals the Genetic Loci and Candidate Genes Governing the Flowering Time in Alfalfa (Medicago sativa L.)

The transition to flowering at the right time is very important for adapting to local conditions and maximizing alfalfa yield. However, the understanding of the genetic basis of the alfalfa flowering time remains limited. There are few reliable genes or markers for selection, which hinders progress in genetic research and molecular breeding of this trait in alfalfa. We sequenced 220 alfalfa cultivars and conducted a genome-wide association study (GWAS) involving 875,023 single-nucleotide polymorphisms (SNPs). The phenotypic analysis showed that the breeding status and geographical origin strongly influenced the alfalfa flowering time. Our GWAS revealed 63 loci significantly related to the flowering time. Ninety-five candidate genes were detected at these SNP loci within 40 kb (20 kb up- and downstream). Thirty-six percent of the candidate genes are involved in development and pollen tube growth, indicating that these genes are key genetic mechanisms of alfalfa growth and development. The transcriptomic analysis showed that 1,924, 2,405, and 3,779 differentially expressed genes (DEGs) were upregulated across the three growth stages, while 1,651, 2,613, and 4,730 DEGs were downregulated across the stages. Combining the results of our GWAS and transcriptome analysis, in total, 38 candidate genes (7 differentially expressed during the bud stage, 13 differentially expressed during the initial flowering stage, and 18 differentially expressed during the full flowering stage) were identified. Two SNPs located in the upstream region of the Msa0888690 gene (which is involved in isop renoids) were significantly related to flowering. The two significant SNPs within the upstream region of Msa0888690 existed as four different haplotypes in this panel. The genes identified in this study represent a series of candidate targets for further research investigating the alfalfa flowering time and could be used for alfalfa molecular breeding.

Whole-Genome Sequencing of KMR3 and Oryza rufipogon-Derived Introgression Line IL50-13 (Chinsurah Nona 2/Gosaba 6) Identifies Candidate Genes for High Yield and Salinity Tolerance in Rice

The genomes of an elite rice restorer line KMR3 (salinity-sensitive) and its salinity-tolerant introgression line IL50-13, a popular variety of coastal West Bengal, India, were sequenced. High-quality paired-end reads were obtained for KMR3 (147.6 million) and IL50-13 (131.4 million) with a sequencing coverage of 30X-39X. Scaffolds generated from the pre-assembled contigs of each sequenced genome were mapped separately onto the reference genome of Oryza sativa ssp. japonica cultivar Nipponbare to identify genomic variants in terms of SNPs and InDels. The SNPs and InDels identified for KMR3 and IL50-13 were then compared with each other to identify polymorphic SNPs and InDels unique and common to both the genomes. Functional enrichment analysis of the protein-coding genes with unique InDels identified GO terms involved in protein modification, ubiquitination, deubiquitination, peroxidase activity, and antioxidant activity in IL50-13. Linoleic acid metabolism, circadian rhythm, and alpha-linolenic acid metabolism pathways were enriched in IL50-13. These GO terms and pathways are involved in reducing oxidative damage, thus suggesting their role in stress responses. Sequence analysis of QTL markers or genes known to be associated with grain yield and salinity tolerance showed polymorphism in 20 genes, out of which nine were not previously reported. These candidate genes encoded Nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4 (NB-ARC) domain-containing protein, cyclase, receptor-like kinase, topoisomerase II-associated protein PAT1 domain-containing protein, ion channel regulatory protein, UNC-93 domain-containing protein, subunit A of the heteromeric ATP-citrate lyase, and three conserved hypothetical genes. Polymorphism was observed in the coding, intron, and untranslated regions of the genes on chromosomes 1, 2, 4, 7, 11, and 12. Genes showing polymorphism between the two genomes were considered as sequence-based new candidates derived from Oryza rufipogon for conferring high yield and salinity tolerance in IL50-13 for further functional studies.

High-LD SNP markers exhibiting pleiotropic effects on salt tolerance at germination and seedlings stages in spring wheat

Salt tolerance at germination and seedling growth stages was investigated. GWAS revealed nine genomic regions with pleiotropic effects on salt tolerance. Salt tolerant genotypes were identified for future breeding program. With 20% of the irrigated land worldwide affected by it, salinity is a serious threat to plant development and crop production. While wheat is the most stable food source worldwide, it has been classified as moderately tolerant to salinity. In several crop plants; such as barley, maize and rice, it has been shown that salinity tolerance at seed germination and seedling establishment is under polygenic control. As yield was the ultimate goal of breeders and geneticists, less attention has been paid to understanding the genetic architecture of salt tolerance at early stages. Thus, the genetic control of salt tolerance at these stages is poorly understood relative to the late stages. In the current study, 176 genotypes of spring wheat were tested for salinity tolerance at seed germination and seedling establishment. Genome-Wide Association Study (GWAS) has been used to identify the genomic regions/genes conferring salt tolerance at seed germination and seedling establishment. Salinity stress negatively impacted all germination and seedling development parameters. A set of 137 SNPs showed significant association with the traits of interest. Across the whole genome, 33 regions showed high linkage disequilibrium (LD). These high LD regions harbored 15 SNPs with pleiotropic effect (i.e. SNPs that control more than one trait). Nine genes belonging to different functional groups were found to be associated with the pleiotropic SNPs. Noteworthy, chromosome 2B harbored the gene TraesCS2B02G135900 that acts as a potassium transporter. Remarkably, one SNP marker, reported in an early study, associated with salt tolerance was validated in this study. Our findings represent potential targets of genetic manipulation to understand and improve salinity tolerance in wheat.

Genome assembly of alfalfa cultivar zhongmu-4 and identification of SNPs associated with agronomic traits

Alfalfa (Medicago sativa L.) is the most important legume forage crop worldwide with high nutritional value and yield. For a long time, the breeding of alfalfa was hampered by lacking reliable information on the autotetraploid genome and molecular markers linked to important agronomic traits. We herein reported the de novo assembly of the allele-aware chromosome-level genome of Zhongmu-4, a cultivar widely cultivated in China, and a comprehensive database of genomic variations based on resequencing of 220 germplasms. Approximate 2.74 Gb contigs (N50 of 2.06 Mb), accounting for 88.39% of the estimated genome, were assembled, and 2.56 Gb contigs were anchored to 32 pseudo-chromosomes. A total of 34,922 allelic genes were identified from the allele-aware genome. We observed the expansion of gene families, especially those related to the nitrogen metabolism, and the increase of repetitive elements including transposable elements, which probably resulted in the increase of Zhongmu-4 genome compared with Medicago truncatula. Population structure analysis revealed that the accessions from Asia and South America had relatively lower genetic diversity than those from Europe, suggesting that geography may influence alfalfa genetic divergence during local adaption. Genome-wide association studies identified 101 single nucleotide polymorphisms (SNPs) associated with 27 agronomic traits. Two candidate genes were predicted to be correlated with fall dormancy and salt response. We believe that the allele-aware chromosome-level genome sequence of Zhongmu-4 combined with the resequencing data of the diverse alfalfa germplasms will facilitate genetic research and genomics-assisted breeding in variety improvement of alfalfa.

Strategies to Increase Prediction Accuracy in Genomic Selection of Complex Traits in Alfalfa (Medicago sativa L.)

Agronomic traits such as biomass yield and abiotic stress tolerance are genetically complex and challenging to improve through conventional breeding approaches. Genomic selection (GS) is an alternative approach in which genome-wide markers are used to determine the genomic estimated breeding value (GEBV) of individuals in a population. In alfalfa (Medicago sativa L.), previous results indicated that low to moderate prediction accuracy values (<70%) were obtained in complex traits, such as yield and abiotic stress resistance. There is a need to increase the prediction value in order to employ GS in breeding programs. In this paper we reviewed different statistic models and their applications in polyploid crops, such as alfalfa and potato. Specifically, we used empirical data affiliated with alfalfa yield under salt stress to investigate approaches that use DNA marker importance values derived from machine learning models, and genome-wide association studies (GWAS) of marker-trait association scores based on different GWASpoly models, in weighted GBLUP analyses. This approach increased prediction accuracies from 50% to more than 80% for alfalfa yield under salt stress. Finally, we expended the weighted GBLUP approach to potato and analyzed 13 phenotypic traits and obtained similar results. This is the first report on alfalfa to use variable importance and GWAS-assisted approaches to increase the prediction accuracy of GS, thus helping to select superior alfalfa lines based on their GEBVs.

A global alfalfa diversity panel reveals genomic selection signatures in Chinese varieties and genomic associations with root development

Alfalfa (Medicago sativa L.) is an important forage crop worldwide. However, little is known about the effects of breeding status and different geographical populations on alfalfa improvement. Here, we sequenced 220 alfalfa core germplasms and determined that Chinese alfalfa cultivars form an independent group, as evidenced by comparisons of F_ST values between different subgroups, suggesting that geographical origin plays an important role in group differentiation. By tracing the influence of geographical regions on the genetic diversity of alfalfa varieties in China, we identified 350 common candidate genetic regions and 548 genes under selection. We also defined 165 loci associated with 24 important traits from genome-wide association studies. Of those, 17 genomic regions closely associated with a given phenotype were under selection, with the underlying haplotypes showing significant differences between subgroups of distinct geographical origins. Based on results from expression analysis and association mapping, we propose that 6-phosphogluconolactonase (MsPGL) and a gene encoding a protein with NHL domains (MsNHL) are critical candidate genes for root growth. In conclusion, our results provide valuable information for alfalfa improvement via molecular breeding.

Progressive Genomic Approaches to Explore Drought- and Salt-Induced Oxidative Stress Responses in Plants under Changing Climate

Drought and salinity are the major environmental abiotic stresses that negatively impact crop development and yield. To improve yields under abiotic stress conditions, drought- and salinity-tolerant crops are key to support world crop production and mitigate the demand of the growing world population. Nevertheless, plant responses to abiotic stresses are highly complex and controlled by networks of genetic and ecological factors that are the main targets of crop breeding programs. Several genomics strategies are employed to improve crop productivity under abiotic stress conditions, but traditional techniques are not sufficient to prevent stress-related losses in productivity. Within the last decade, modern genomics studies have advanced our capabilities of improving crop genetics, especially those traits relevant to abiotic stress management. This review provided updated and comprehensive knowledge concerning all possible combinations of advanced genomics tools and the gene regulatory network of reactive oxygen species homeostasis for the appropriate planning of future breeding programs, which will assist sustainable crop production under salinity and drought conditions.

Annual and perennial Medicago show signatures of parallel adaptation to climate and soil in highly conserved genes

Human induced environmental change may require rapid adaptation of plant populations and crops, but the genomic basis of environmental adaptation remain poorly understood. We analysed polymorphic loci from the perennial crop Medicago sativa (alfalfa or lucerne) and the annual legume model species M. truncatula to search for a common set of candidate genes that might contribute to adaptation to abiotic stress in both annual and perennial Medicago species. We identified a set of candidate genes of adaptation associated with environmental gradients along the distribution of the two Medicago species. Candidate genes for each species were detected in homologous genomic linkage blocks using genome-environment (GEA) and genome-phenotype association analyses. Hundreds of GEA candidate genes were species-specific, of these, 13.4% (M. sativa) and 24% (M. truncatula) were also significantly associated with phenotypic traits. A set of 168 GEA candidates were shared by both species, which was 25.4% more than expected by chance. When combined, they explained a high proportion of variance for certain phenotypic traits associated with adaptation. Genes with highly conserved functions dominated among the shared candidates and were enriched in gene ontology terms that have shown to play a central role in drought avoidance and tolerance mechanisms by means of cellular shape modifications and other functions associated with cell homeostasis. Our results point to the existence of a molecular basis of adaptation to abiotic stress in Medicago determined by highly conserved genes and gene functions. We discuss these results in light of the recently proposed omnigenic model of complex traits.

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

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.

The Chromosome-Level Genome Sequence of the Autotetraploid Alfalfa and Resequencing of Core Germplasms Provide Genomic Resources for Alfalfa Research

Alfalfa (Medicago sativa) is one of the most important forage crops in the world; however, its molecular genetics and breeding research are hindered due to the lack of a high-quality reference genome. Here, we report a de novo assembled 816-Mb high-quality, chromosome-level haploid genome sequence for 'Zhongmu No.1' alfalfa, a heterozygous autotetraploid. The contig N50 is 3.92 Mb, and 49 165 genes are annotated in the genome. The alfalfa genome is estimated to have diverged from M. truncatula approximately 8 million years ago. Genomic population analysis of 162 alfalfa accessions revealed high genetic diversity, weak population structure, and extensive gene flow from wild to cultivated alfalfa. Genome-wide association studies identified many candidate genes associated with important agronomic traits. Furthermore, we showed that MsFTa2, a Flowering Locus T homolog, whose expression is upregulated in salt-resistant germplasms, may be associated with fall dormancy and salt resistance. Taken together, these genomic resources will facilitate alfalfa genetic research and agronomic improvement.

Overexpression of MsSiR enhances alkali tolerance in alfalfa (Medicago sativa L.) by increasing the glutathione content

The sulfite reductase gene in Medicago sativa L. (MsSiR) encodes sulfite reductase (SiR) and catalyses the conversion of sulfite to sulfate in the sulfite assimilation pathway. In this study, we investigated the role of MsSiR in alfalfa by generating transgenic alfalfa that ectopically expressed MsSiR under the control of the CaMV35S promoter. The differences in alkali tolerance between the MsSiR-overexpressing and wild-type (WT) plants were analyzed, and the MsSiR-overexpressing plants exhibited an improved phenotype under alkali stress. Compared to WT plants, these plants demonstrated improved antioxidant activity as well as decreased H₂O₂ and O₂^- contents and increased glutathione reduced (GSH), Cysteine (Cys) and glutathione oxidized (GSSG) contents. MsSiR-overexpressing plants also exhibited high levels of adenosyl phosphosulfate reductases (APR), sulfite oxidase (SO) and MsSiR expression under alkali stress. It was speculated that MsSiR is involved in sulfur metabolism pathways, including the stabilization of sulfate and sulfite levels and the synthesis of GSH. These two processes achieve alkali tolerance by positively regulating the detoxification and antioxidant activities of alfalfa.

Identification of Genes Conferring Plant Salt Tolerance using GWAS: Current Success and Perspectives

An understanding of the molecular mechanisms that underlie plant salt tolerance is important for both economic and scientific interests. Genome-wide association study (GWAS) is a promising approach to pinpoint genes that confer plant salt tolerance. With the advancement of supporting technology and methodology, GWAS has enabled the discovery of genes that play central roles in regulating plant salt tolerance in the past decade. Here, I highlight recent successful GWAS work in unveiling the molecular factors underlying plant salt tolerance and discuss the concerns and opportunities in conducting such experiments. It is anticipated that GWAS will be increasingly successful in the identification of key genes that are useful for crop improvement.

Biotechnological Perspectives of Omics and Genetic Engineering Methods in Alfalfa

For several decades, researchers are working to develop improved major crops with better adaptability and tolerance to environmental stresses. Forage legumes have been widely spread in the world due to their great ecological and economic values. Abiotic and biotic stresses are main factors limiting legume production, however, alfalfa (Medicago sativa L.) shows relatively high level of tolerance to drought and salt stress. Efforts focused on alfalfa improvements have led to the release of cultivars with new traits of agronomic importance such as high yield, better stress tolerance or forage quality. Alfalfa has very high nutritional value due to its efficient symbiotic association with nitrogen-fixing bacteria, while deep root system can help to prevent soil water loss in dry lands. The use of modern biotechnology tools is challenging in alfalfa since full genome, unlike to its close relative barrel medic (Medicago truncatula Gaertn.), was not released yet. Identification, isolation, and improvement of genes involved in abiotic or biotic stress response significantly contributed to the progress of our understanding how crop plants cope with these environmental challenges. In this review, we provide an overview of the progress that has been made in high-throughput sequencing, characterization of genes for abiotic or biotic stress tolerance, gene editing, as well as proteomic and metabolomics techniques bearing biotechnological potential for alfalfa improvement.

Morphological, Physiological, and Genetic Responses to Salt Stress in Alfalfa: A Review

Alfalfa (Medicago sativa L.) is an important legume forage crop. However, its genetic improvement for salt tolerance is challenging, as alfalfa’s response to salt stress is genetically and physiologically complex. A review was made to update the knowledge of morphological, physiological, biochemical, and genetic responses of alfalfa plants to salt stress, and to discuss the potential of applying modern plant technologies to enhance alfalfa salt-resistant breeding, including genomic selection, RNA-Seq analysis, and cutting-edge Synchrotron beamlines. It is clear that alfalfa salt tolerance can be better characterized, genes conditioning salt tolerance be identified, and new marker-based tools be developed to accelerate alfalfa breeding for salt tolerance.

Genome-Wide Association Study of Salinity Tolerance During Germination in Barley (Hordeum vulgare L.)

Barley seeds need to be able to germinate and establish seedlings in saline soils in Mediterranean-type climates. Despite being a major cereal crop, barley has few reported quantitative trait loci (QTL) and candidate genes underlying salt tolerance at the germination stage. Breeding programs targeting salinity tolerance at germination require an understanding of genetic loci and alleles in the current germplasm. In this study, we investigated seed-germination-related traits under control and salt stress conditions in 350 diverse barley accessions. A genome-wide association study, using ~24,000 genetic markers, was undertaken to detect marker-trait associations (MTA) and the underlying candidate genes for salinity tolerance during germination. We detected 19 loci containing 52 significant salt-tolerance-associated markers across all chromosomes, and 4 genes belonging to 4 family functions underlying the predicted MTAs. Our results provide new genetic resources and information to improve salt tolerance at germination in future barley varieties via genomic and marker-assisted selection and to open up avenues for further functional characterization of the identified candidate genes.

A comparative morphological and transcriptomic study on autotetraploid Stevia rebaudiana (bertoni) and its diploid

Stevia rebaudiana is an important medical plant for producing steviol glycosides (SGs) or stevioside. Autotetraploids (4x = 44) show an increasing level of morphology, physiology and tolerances comparing to diploids (2x = 22). However, little information regarded on the comparative transcriptome analysis between diploid and autotetraploid S. rebaudiana was found. In this study, synthetic autotetraploid was induced and morphological features were confirmed. A comprehensive transcriptome of stevia leaf, stem and root from the diploids and autotetraploids was constructed based on RNA-seq, yielded 1,000,892,422 raw reads and subsequently assembled into 251,455 transcripts, corresponded to 146,130 genes. Pairwise comparisons of the six leaf libraries between the diploids and autotetraploids revealed 4114 differentially expression genes (DEGs), in which 2105 (51.17%) were up-regulated in autotetraploids and associated with SGs biosynthesis, plant growth and secondary metabolism. Moreover, weighted gene co-expression network analysis showed co-expressed genes of fifteen genes of SG biosynthesis pathway were enriched in photosynthesis, flavonoid and secondary metabolic process, plant growth and morphogenesis. A hundred of DEGs related to plant resistance were identified by interviewing PlantPReS database. This study has highlighted molecular changes related to SGs metabolism of polyploidy, and advanced our understanding in plant resistance responsible for phenotypic change of autotetraploids.

Genetic Loci Associated with Salt Tolerance in Advanced Breeding Populations of Tetraploid Alfalfa Using Genome-Wide Association Studies

Many agricultural lands in the western United States consist of soil with high concentrations of salt, which is detrimental to alfalfa ( L.) growth and production, especially in the region where water resource is limited. Developing alfalfa varieties with salt tolerance is imperative for sustainable production under increasing soil salinity. In the present study, we used advanced alfalfa breeding populations and evaluated five traits related to salt tolerance including biomass dry weight (DW) and fresh weight (FW), plant height (PH), leaf relative water content (RWC), and stomatal conductance (SC) under control and salt stress. Stress susceptibility index (SSI) of each trait and single-nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing (GBS) were used for genome-wide association studies (GWAS) to identify loci associated with salt tolerance. A total of 53 significant SNPs associated with salt tolerance were identified and they were located at 49 loci through eight chromosomes. A Basic Local Alignment Search Tool (BLAST) search of the regions surrounding the SNPs revealed 21 putative candidate genes associated with salt tolerance. The genetic architecture for traits related to salt tolerance characterized in this report could help in understanding the genetic mechanism by which salt stress affects plant growth and production in alfalfa. The markers and candidate genes identified in the present study would be useful for marker-assisted selection (MAS) in breeding salt-tolerant alfalfa after validation of the markers.

Understanding the Complexity of Cold Tolerance in White Clover using Temperature Gradient Locations and a GWAS Approach

White clover ( L.) is the most important grazing perennial forage legume in temperate climates. However, its limited capacity to survive and restore growth after low temperatures during winter constrains the productivity and wide adoption of the crop. Despite the importance of cold tolerance for white clover cultivar development, the genetic basis of this trait remains largely unknown. Hence, in this study, we performed the first genome-wide association study (GWAS) analyses in white clover to identify quantitative trait loci (QTL) for cold-tolerance-related traits. Seeds from 192 divergent genotypes from six populations in the Patagonia region of South America were collected and seed-derived plants were further clonally propagated. Clonal trials were established in three locations representing temperature gradient associated with elevation. Given the allotetraploid nature of the white clover genome, distinct genetic models (diploid and tetraploid) were tested. Only the tetraploid parameterization was able to detect the 53 loci associated with cold-tolerance traits. Out of the 53 single nucleotide polymorphism (SNP) trait associations, 17 controlled more than one trait or were stable across multiple sites. This work represents the first report of QTL for cold-tolerance-related traits, providing insights into its genetic basis and candidate genomic regions for further functional validation studies.

A fully automated pipeline for quantitative genotype calling from next generation sequencing data in autopolyploids

BACKGROUND

Genotyping-by-sequencing (GBS) has been used broadly in genetic studies for several species, especially those with agricultural importance. However, its use is still limited in autopolyploid species because genotype calling software generally fails to properly distinguish heterozygous classes based on allele dosage.

RESULTS

VCF2SM is a Python script that integrates sequencing depth information of polymorphisms in variant call format (VCF) files and SUPERMASSA software for quantitative genotype calling. VCFs can be obtained from any variant discovery software that outputs exact allele sequencing depth, such as a modified version of the TASSEL-GBS pipeline provided here. VCF2SM was successfully applied in analyzing GBS data from diverse panels (alfalfa and potato) and full-sib mapping populations (alfalfa and switchgrass) of polyploid species.

CONCLUSIONS

We demonstrate that our approach can help plant geneticists working with autopolyploid species to advance their studies by distinguishing allele dosage from GBS data.

Pedigree-Based Deciphering of Genome-Wide Conserved Patterns in an Elite Potato Parental Line

Elite parental lines are more likely to breed fine varieties, but knowledge about elite parents and their genetic backgrounds is limited. In this paper, we investigated the pedigree relationships of potato varieties bred worldwide and in China. Several elite parents were identified, and these parents were more frequently used as parents in breeding programs across different time periods and countries. We next used 2b-RAD, a reduced-representation sequencing method, to genotype the elite parent Mira and 24 of its offspring. These cultivars span 5 generations, making this lineage the longest continuous pedigree among Chinese bred potatoes. A total of 47,314 tetraploid single nucleotide polymorphisms (SNPs) identified by FreeBayes were used to trace the conserved segments of the Mira genome. The conserved segments had identical or similar allele-specific SNPs across the analyzed genotypes. In Mira, 3,788 segments comprising over 10,000 bp, or 20.8% of the genome, were defined as conserved segments. These segments contain genes involved in crucial biological processes that are of special interest to breeders. These regions, which have been conserved across generations of highly selective breeding, may be helpful for further breeding and performing genome-wide breeding by design.

Genetic-based dissection of arsenic accumulation in maize using a genome-wide association analysis method

Understanding the mechanism of arsenic (As) accumulation in plants is important in reducing As's toxicity to plants and its potential risks to human health. Here, we performed a genome-wide association study to dissect the genetic basis of the As contents of different maize tissues in Xixian, which was irrigated with As-rich surface water, and Changge using an association population consisting of 230 representative maize inbred lines. Phenotypic data revealed a wide normal distribution and high repeatability for the As contents in maize tissues. The As concentrations in maize tissues followed the same trend in the two locations: kernels < axes < stems < bracts < leaves. In total, 15, 16 and 15 non-redundant quantitative trait loci (QTLs) associated with As concentrations were identified (P ≤ 2.04 × 10-6 ) in five tissues from Xixian, Changge, and the combination of the locations, respectively, explaining 9.70%-24.65% of the phenotypic variation for each QTL, on average. Additionally, four QTLs [involving 15 single nucleotide polymorphisms (SNPs)] were detected in the single and the combined locations, indicating that these loci/SNPs might be stable across different environments. The candidate genes associated with these four loci were predicted. In addition, four non-redundant QTLs (6 SNPs), including a QTL that was detected in multiple locations according to the genome-wide association study, were found to co-localize with four previously reported QTL intervals. These results are valuable to understand the genetic architecture of As mechanism in maize and facilitate the genetic improvement of varieties without As toxicity.

Ectopic Expression of GsSRK in Medicago sativa Reveals Its Involvement in Plant Architecture and Salt Stress Responses

Receptor-like kinases (RLK) play fundamental roles in plant growth and stress responses. Compared with other RLKs, little information is provided concerning the S-locus LecRLK subfamily, which is characterized by an extracellular G-type lectin domain and an S-locus-glycop domain. Until now, the function of the G-type lectin domain is still unknown. In a previous research, we identified a Glycine soja S-locus LecRLK gene GsSRK, which conferred increased salt stress tolerance in transgenic Arabidopsis. In this study, to investigate the role of the G-type lectin domain and to breed transgenic alfalfa with superior salt stress tolerance, we transformed the full-length GsSRK (GsSRK-f) and a truncated version of GsSRK (GsSRK-t) deleting the G-type lectin domain into alfalfa. Our results showed that overexpression of GsSRK-t, but not GsSRK-f, resulted in changes of plant architecture, as evidenced by more branches but shorter shoots of GsSRK-t transgenic alfalfa, indicating a potential role of the extracellular G-type lectin domain in regulating plant architecture. Furthermore, we also found that transgenic alfalfa overexpressing either GsSRK-f or GsSRK-t showed increased salt stress tolerance, and GsSRK-t transgenic alfalfa displayed better growth (more branches and higher fresh weight) than GsSRK-f lines under salt stress. In addition, our results suggested that both GsSRK-f and GsSRK-t were involved in ion homeostasis, ROS scavenging, and osmotic regulation. Under salt stress, the Na+ content in the transgenic lines was significantly lower, while the K+ content was slightly higher than that in WT. Moreover, the transgenic lines displayed reduced ion leakage and MDA content, but increased SOD activity and proline content than WT. Notably, no obvious difference in these physiological indices was observed between GsSRK-f and GsSRK-t transgenic lines, implying that deletion of the GsSRK G-type lectin domain does not affect its physiological function in salt stress responses. In conclusion, results in this research reveal the dual role of GsSRK in regulating both plant architecture and salt stress responses.

Comparative analysis of alfalfa (Medicago sativa L.) leaf transcriptomes reveals genotype-specific salt tolerance mechanisms

BACKGROUND

Soil salinity is an important factor affecting growth, development, and productivity of almost all land plants, including the forage crop alfalfa (Medicago sativa). However, little is known about how alfalfa responds and adapts to salt stress, particularly among different salt-tolerant cultivars.

RESULTS

Among seven alfalfa cultivars, we found that Zhongmu-1 (ZM) is relatively salt-tolerant and Xingjiang Daye (XJ) is salt-sensitive. Compared to XJ, ZM showed slower growth under low-salt conditions, but exhibited stronger tolerance to salt stress. RNA-seq analysis revealed 2237 and 1125 differentially expressed genes (DEGs) between ZM and XJ in the presence and absence of salt stress, among which many genes are involved in stress-related pathways. After salt treatment, compared with the controls, the number of DEGs in XJ (19373) was about four times of that in ZM (4833). We also detected specific differential gene expression patterns: In response to salt stress, compared with XJ, ZM maintained relatively more stable expression levels of genes related to the ROS and Ca²⁺ pathways, phytohormone biosynthesis, and Na⁺/K⁺ transport. Notably, several salt resistance-associated genes always showed greater levels of expression in ZM than in XJ, including a transcription factor. Consistent with the suppression of plant growth resulting from salt stress, the expression of numerous photosynthesis- and growth hormone-related genes decreased more dramatically in XJ than in ZM. By contrast, the expression levels of photosynthetic genes were lower in ZM under low-salt conditions.

CONCLUSIONS

Compared with XJ, ZM is a salt-tolerant alfalfa cultivar possessing specific regulatory mechanisms conferring exceptional salt tolerance, likely by maintaining high transcript levels of abiotic and biotic stress resistance-related genes. Our results suggest that maintaining this specific physiological status and/or plant adaptation to salt stress most likely arises by inhibition of plant growth in ZM through plant hormone interactions. This study identifies new candidate genes that may regulate alfalfa tolerance to salt stress and increases the understanding of the genetic basis for salt tolerance.

Salt-tolerant and -sensitive alfalfa (Medicago sativa) cultivars have large variations in defense responses to the lepidopteran insect Spodoptera litura under normal and salt stress condition

In nature, plants are often exposed to multiple stress factors at the same time. Yet, little is known about how plants modulate their physiology to counteract simultaneous abiotic and biotic stresses, such as soil salinity and insect herbivory. In this study, insect performance bioassays, phytohormone measurements, quantification of transcripts, and protein determination were employed to study the phenotypic variations of two alfalfa (Medicago sativa) cultivars in response to insect Spodoptera litura feeding under normal and salt stress condition. When being cultivated in normal soil, the salt-tolerant alfalfa cultivar Zhongmu-1 exhibited lower insect resistance than did the salt-sensitive cultivar Xinjiang Daye. Under salinity stress, the defense responses of Xinjiang Daye were repressed, whereas Zhongmu-1 did not show changes in resistance levels. It is likely that salinity influenced the resistance of Xinjiang Daye through suppressing the accumulation of jasmonic acid-isoleucine (JA-Ile), which is the bioactive hormone inducing herbivore defense responses, leading to attenuated trypsin proteinase inhibitor (TPI) activity. Furthermore, exogenous ABA supplementation suppressed the insect herbivory-induced JA/JA-Ile accumulation and levels of JAR1 (jasmonate resistant 1) and TPI, and further decreased the resistance of Xinjiang Daye, whereas Zhongmu-1 showed very little response to the increased ABA level. We propose a mechanism, in which high levels of abscisic acid induced by salt treatment may affect the expression levels of JAR1 and consequently decrease JA-Ile accumulation and thus partly suppress the defense of Xinjiang Daye against insects under salt stress. This study provides new insight into the mechanism by which alfalfa responds to concurrent abiotic and biotic stresses.

Identification of genetic loci associated with crude protein and mineral concentrations in alfalfa (Medicago sativa) using association mapping

BACKGROUND

Alfalfa (Medicago sativa) is one of the most important legume forage species in China and many other countries of the world. It provides a quality source of proteins and minerals to animals. Genetic underpinnings for these important traits, however, are elusive. An alfalfa (M. sativa) association mapping study for six traits, namely crude protein (CP), rumen undegraded protein (RUP), and four mineral elements (Ca, K, Mg and P), was conducted in three consecutive years using a large collection encompassing 336 genotypes genotyped with 85 simple sequence repeat (SSR) markers.

RESULTS

All the traits were significantly influenced by genotype, environment, and genotype × environment interaction. Eight-five significant associations (P < 0.005) were identified. Among these, five associations with Ca were repeatedly observed and six co-localized associations were identified.

CONCLUSIONS

The identified marker alleles significantly associated with the traits provided important information for understanding genetic controls of alfalfa quality. The markers could be used in assisting selection for the individual traits in breeding populations for developing new alfalfa cultivars.

Genome-Wide Association Mapping of Loci Associated with Plant Growth and Forage Production under Salt Stress in Alfalfa (Medicago sativa L.)

Salinity tolerance is highly desirable to sustain alfalfa production in marginal lands that have been rendered saline. In this study, we used a diverse panel of 198 alfalfa accessions for mapping loci associated with plant growth and forage production under salt stress using genome-wide association studies (GWAS). The plants were genotyped using genotyping-by-sequencing (GBS). A greenhouse procedure was used for phenotyping four agronomic and physiological traits affected by salt stress, including dry weight (DW), plant height (PH), leaf chlorophyll content (LCC), and stomatal conductance (SC). For each trait, a stress susceptibility index (SSI) was used to evaluate plant performance under stressed and non-stressed conditions. Marker-trait association identified a total of 42 markers significantly associated with salt tolerance. They were located on all chromosomes except chromosome 2 based on the alignment of their flanking sequences to the reference genome (Medicago truncatula). Of those identified, 13 were associated with multiple traits. Several loci identified in the present study were also identified in previous reports. BLAST search revealed that 19 putative candidate genes linked to 24 significant markers. Among them, B3 DNA-binding protein, Thiaminepyrophosphokinase and IQ calmodulin-binding motif protein were identified among multiple traits in the present and previous studies. With further investigation, these markers and candidates would be useful for developing markers for marker-assisted selection in breeding programs to improve alfalfa cultivars with enhanced tolerance to salt stress.

Systems biology and genome-wide approaches to unveil the molecular players involved in the pre-germinative metabolism: implications on seed technology traits

The pre-germinative metabolism is among the most fascinating aspects of seed biology. The early seed germination phase, or pre-germination, is characterized by rapid water uptake (imbibition), which directs a series of dynamic biochemical events. Among those are enzyme activation, DNA damage and repair, and use of reserve storage compounds, such as lipids, carbohydrates and proteins. Industrial seedling production and intensive agricultural production systems require seed stocks with high rate of synchronized germination and low dormancy. Consequently, seed dormancy, a quantitative trait related to the activation of the pre-germinative metabolism, is probably the most studied seed trait in model species and crops. Single omics, systems biology, QTLs and GWAS mapping approaches have unveiled a list of molecules and regulatory mechanisms acting at transcriptional, post-transcriptional and post-translational levels. Most of the identified candidate genes encode for regulatory proteins targeting ROS, phytohormone and primary metabolisms, corroborating the data obtained from simple molecular biology approaches. Emerging evidences show that epigenetic regulation plays a crucial role in the regulation of these mentioned processes, constituting a still unexploited strategy to modulate seed traits. The present review will provide an up-date of the current knowledge on seed pre-germinative metabolism, gathering the most relevant results from physiological, genetics, and omics studies conducted in model and crop plants. The effects exerted by the biotic and abiotic stresses and priming are also addressed. The possible implications derived from the modulation of pre-germinative metabolism will be discussed from the point of view of seed quality and technology.

Genome-Wide Single-Nucleotide Polymorphisms in CMS and Restorer Lines Discovered by Genotyping Using Sequencing and Association with Marker-Combining Ability for 12 Yield-Related Traits in Oryza sativa L. subsp. Japonica

Heterosis or hybrid vigor is closely related with general combing ability (GCA) of parents and special combining ability (SCA) of combinations. The evaluation of GCA and SCA facilitate selection of parents and combinations in heterosis breeding. In order to improve combining ability (CA) by molecular marker assist selection, it is necessary to identify marker loci associated with the CA. To identify the single nucleotide polymorphisms (SNP) loci associated with CA in the parental genomes of japonica rice, genome-wide discovered SNP loci were tested for association with the CA of 18 parents for 12 yield-related traits. In this study, 81 hybrids were created and evaluated to calculate the CA of 18 parents. The parents were sequenced by genotyping by sequencing (GBS) method for identification of genome-wide SNPs. The analysis of GBS indicated that the successful mapping of 9.86 × 10⁶ short reads in the Nipponbare reference genome consists of 39,001 SNPs in parental genomes at 11,085 chromosomal positions. The discovered SNPs were non-randomly distributed within and among the 12 chromosomes of rice. Overall, 20.4% (8026) of the discovered SNPs were coding types, and 8.6% (3344) and 9.9% (3951) of the SNPs revealed synonymous and non-synonymous changes, which provide valuable knowledge about the underlying performance of the parents. Furthermore, the associations between SNPs and CA indicated that 362 SNP loci were significantly related to the CA of 12 parental traits. The identified SNP loci of CA in our study were distributed genome wide and caused a positive or negative effect on the CA of traits. For the yield-related traits, such as grain thickness, days to heading, panicle length, grain length and 1000-grain weight, a maximum number of positive SNP loci of CA were found in CMS A171 and in the restorers LC64 and LR27. On an individual basis, some of associated loci that resided on chromosomes 2, 5, 7, 9, and 11 recorded maximum positive values for the CA of traits. From our results, we suggest that heterosis in japonica rice would be improved by pyramiding the favorable SNP loci of CA and eliminating the unfavorable loci from parental genomes.

Genome-Wide Association Mapping and Genomic Selection for Alfalfa (Medicago sativa) Forage Quality Traits

Genetic progress for forage quality has been poor in alfalfa (Medicago sativa L.), the most-grown forage legume worldwide. This study aimed at exploring opportunities for marker-assisted selection (MAS) and genomic selection of forage quality traits based on breeding values of parent plants. Some 154 genotypes from a broadly-based reference population were genotyped by genotyping-by-sequencing (GBS), and phenotyped for leaf-to-stem ratio, leaf and stem contents of protein, neutral detergent fiber (NDF) and acid detergent lignin (ADL), and leaf and stem NDF digestibility after 24 hours (NDFD), of their dense-planted half-sib progenies in three growing conditions (summer harvest, full irrigation; summer harvest, suspended irrigation; autumn harvest). Trait-marker analyses were performed on progeny values averaged over conditions, owing to modest germplasm × condition interaction. Genomic selection exploited 11,450 polymorphic SNP markers, whereas a subset of 8,494 M. truncatula-aligned markers were used for a genome-wide association study (GWAS). GWAS confirmed the polygenic control of quality traits and, in agreement with phenotypic correlations, indicated substantially different genetic control of a given trait in stems and leaves. It detected several SNPs in different annotated genes that were highly linked to stem protein content. Also, it identified a small genomic region on chromosome 8 with high concentration of annotated genes associated with leaf ADL, including one gene probably involved in the lignin pathway. Three genomic selection models, i.e., Ridge-regression BLUP, Bayes B and Bayesian Lasso, displayed similar prediction accuracy, whereas SVR-lin was less accurate. Accuracy values were moderate (0.3-0.4) for stem NDFD and leaf protein content, modest for leaf ADL and NDFD, and low to very low for the other traits. Along with previous results for the same germplasm set, this study indicates that GBS data can be exploited to improve both quality traits (by genomic selection or MAS) and forage yield.