Genome-wide association analysis of salinity responsive traits in Medicago truncatula

Genetic diversity, population structure, and genome-wide association study for the flowering trait in a diverse panel of 428 moth bean (Vigna aconitifolia) accessions using genotyping by sequencing

BACKGROUND

Moth bean (Vigna aconitifolia) is an underutilized, protein-rich legume that is grown in arid and semi-arid areas of south Asia and is highly resistant to abiotic stresses such as heat and drought. Despite its economic importance, the crop remains unexplored at the genomic level for genetic diversity and trait mapping studies. To date, there is no report of SNP marker discovery and association mapping of any trait in this crop. Therefore, this study aimed to dissect the genetic diversity, population structure and marker-trait association for the flowering trait in a diversity panel of 428 moth bean accessions using genotyping by sequencing (GBS) approach.

RESULTS

A total of 9078 high-quality single nucleotide polymorphisms (SNPs) were discovered by genotyping of 428 moth bean accessions. Model-based structure analysis and PCA grouped the moth bean accessions into two subpopulations. Cluster analysis revealed accessions belonging to the Northwestern region of India had higher variability than accessions from the other regions suggesting that this region represents its center of diversity. AMOVA revealed more variations within individuals (74%) and among the individuals (24%) than among the populations (2%). Marker-trait association analysis using seven multi-locus models including mrMLM, FASTmrEMMA FASTmrEMMA, ISIS EM-BLASSO, MLMM, BLINK and FarmCPU revealed 29 potential genomic regions for the trait days to 50% flowering, which were consistently detected in three or more models. Analysis of the allelic effect of the major genomic regions explaining phenotypic variance of more than 10% and those detected in at least 2 environments showed 4 genomic regions with significant phenotypic effect on this trait. Further, we also analyzed genetic relationships among the Vigna species using SNP markers. The genomic localization of moth bean SNPs on genomes of closely related Vigna species demonstrated that maximum numbers of SNPs were getting localized on Vigna mungo. This suggested that the moth bean is most closely related to V. mungo.

CONCLUSION

Our study shows that the north-western regions of India represent the center of diversity of the moth bean. Further, the study revealed flowering-related genomic regions/candidate genes which can be potentially exploited in breeding programs to develop early-maturity moth bean varieties.

Medicago truncatula resources to study legume biology and symbiotic nitrogen fixation

Medicago truncatula is a chosen model for legumes towards deciphering fundamental legume biology, especially symbiotic nitrogen fixation. Current genomic resources for M. truncatula include a completed whole genome sequence information for R108 and Jemalong A17 accessions along with the sparse draft genome sequences for other 226 M. truncatula accessions. These genomic resources are complemented by the availability of mutant resources such as retrotransposon (Tnt1) insertion mutants in R108 and fast neutron bombardment (FNB) mutants in A17. In addition, several M. truncatula databases such as small secreted peptides (SSPs) database, transporter protein database, gene expression atlas, proteomic atlas, and metabolite atlas are available to the research community. This review describes these resources and provide information regarding how to access these resources.

High-throughput phenotyping-based quantitative trait loci mapping reveals the genetic architecture of the salt stress tolerance of Brassica napus

Salt stress is a major limiting factor that severely affects the survival and growth of crops. It is important to understand the salt stress tolerance ability of Brassica napus and explore the underlying related genetic resources. We used a high-throughput phenotyping platform to quantify 2111 image-based traits (i-traits) of a natural population under three different salt stress conditions and an intervarietal substitution line (ISL) population under nine different stress conditions to monitor and evaluate the salt stress tolerance of B. napus over time. We finally identified 928 high-quality i-traits associated with the salt stress tolerance of B. napus. Moreover, we mapped the salt stress-related loci in the natural population via a genome-wide association study and performed a linkage analysis associated with the ISL population, respectively. These results revealed 234 candidate genes associated with salt stress response, and two novel candidate genes, BnCKX5 and BnERF3, were experimentally verified to regulate the salt stress tolerance of B. napus. This study demonstrates the feasibility of using high-throughput phenotyping-based quantitative trait loci mapping to accurately and comprehensively quantify i-traits associated with B. napus. The mapped loci could be used for genomics-assisted breeding to genetically improve the salt stress tolerance of B. napus.

Genetic Augmentation of Legume Crops Using Genomic Resources and Genotyping Platforms for Nutritional Food Security

Recent advances in next generation sequencing (NGS) technologies have led the surge of genomic resources for the improvement legume crops. Advances in high throughput genotyping (HTG) and high throughput phenotyping (HTP) enable legume breeders to improve legume crops more precisely and efficiently. Now, the legume breeder can reshuffle the natural gene combinations of their choice to enhance the genetic potential of crops. These genomic resources are efficiently deployed through molecular breeding approaches for genetic augmentation of important legume crops, such as chickpea, cowpea, pigeonpea, groundnut, common bean, lentil, pea, as well as other underutilized legume crops. In the future, advances in NGS, HTG, and HTP technologies will help in the identification and assembly of superior haplotypes to tailor the legume crop varieties through haplotype-based breeding. This review article focuses on the recent development of genomic resource databases and their deployment in legume molecular breeding programmes to secure global food security.

How salt stress-responsive proteins regulate plant adaptation to saline conditions

An overview is presented of recent advances in our knowledge of candidate proteins that regulate various physiological and biochemical processes underpinning plant adaptation to saline conditions. Salt stress is one of the environmental constraints that restrict plant distribution, growth and yield in many parts of the world. Increased world population surely elevates food demands all over the globe, which anticipates to add a great challenge to humanity. These concerns have necessitated the scientists to understand and unmask the puzzle of plant salt tolerance mechanisms in order to utilize various strategies to develop salt tolerant crop plants. Salt tolerance is a complex trait involving alterations in physiological, biochemical, and molecular processes. These alterations are a result of genomic and proteomic complement readjustments that lead to tolerance mechanisms. Proteomics is a crucial molecular tool that indicates proteins expressed by the genome, and also identifies the functions of proteins accumulated in response to salt stress. Recently, proteomic studies have shed more light on a range of promising candidate proteins that regulate various processes rendering salt tolerance to plants. These proteins have been shown to be involved in photosynthesis and energy metabolism, ion homeostasis, gene transcription and protein biosynthesis, compatible solute production, hormone modulation, cell wall structure modification, cellular detoxification, membrane stabilization, and signal transduction. These candidate salt responsive proteins can be therefore used in biotechnological approaches to improve tolerance of crop plants to salt conditions. In this review, we provided comprehensive updated information on the proteomic data of plants/genotypes contrasting in salt tolerance in response to salt stress. The roles of salt responsive proteins that are potential determinants for plant salt adaptation are discussed. The relationship between changes in proteome composition and abundance, and alterations observed in physiological and biochemical features associated with salt tolerance are also addressed.

High-Altitude Genetic Selection and Genome-Wide Association Analysis of Yield-Related Traits in Elymus sibiricus L. Using SLAF Sequencing

The genetic adaptations to harsh climatic conditions in high altitudes and genetic basis of important agronomic traits are poorly understood in Elymus sibiricus L. In this study, an association population of 210 genotypes was used for population structure, selective sweep analysis, and genome-wide association study (GWAS) based on 88,506 single nucleotide polymorphisms (SNPs). We found 965 alleles under the natural selection of high altitude, which included 7 hub genes involved in the response to UV, and flavonoid and anthocyanin biosynthetic process based on the protein-protein interaction (PPI) analysis. Using a mixed linear model (MLM), the GWAS test identified a total of 1,825 significant loci associated with 12 agronomic traits. Based on the gene expression data of two wheat cultivars and the PPI analysis, we finally identified 12 hub genes. Especially, in plant height traits, the top hub gene (TOPLESS protein) encoding auxins and jasmonic acid signaling pathway, shoot apical meristem specification, and xylem and phloem pattern formation was highly overexpressed. These genes might play essential roles in controlling the growth and development of E. sibiricus. Therefore, this study provides fundamental insights relevant to hub genes and will benefit molecular breeding and improvement in E. sibiricus and other Elymus species.

Genetic architecture of root and shoot ionomes in rice (Oryza sativa L.)

KEY MESSAGE

Association analysis for ionomic concentrations of 20 elements identified independent genetic factors underlying the root and shoot ionomes of rice, providing a platform for selecting and dissecting causal genetic variants. Understanding the genetic basis of mineral nutrient acquisition is key to fully describing how terrestrial organisms interact with the non-living environment. Rice (Oryza sativa L.) serves both as a model organism for genetic studies and as an important component of the global food system. Studies in rice ionomics have primarily focused on above ground tissues evaluated from field-grown plants. Here, we describe a comprehensive study of the genetic basis of the rice ionome in both roots and shoots of 6-week-old rice plants for 20 elements using a controlled hydroponics growth system. Building on the wealth of publicly available rice genomic resources, including a panel of 373 diverse rice lines, 4.8 M genome-wide single-nucleotide polymorphisms, single- and multi-marker analysis pipelines, an extensive tome of 321 candidate genes and legacy QTLs from across 15 years of rice genetics literature, we used genome-wide association analysis and biparental QTL analysis to identify 114 genomic regions associated with ionomic variation. The genetic basis for root and shoot ionomes was highly distinct; 78 loci were associated with roots and 36 loci with shoots, with no overlapping genomic regions for the same element across tissues. We further describe the distribution of phenotypic variation across haplotypes and identify candidate genes within highly significant regions associated with sulfur, manganese, cadmium, and molybdenum. Our analysis provides critical insight into the genetic basis of natural phenotypic variation for both root and shoot ionomes in rice and provides a comprehensive resource for dissecting and testing causal genetic variants.

Genome-Wide Association Studies of Seed Performance Traits in Response to Heat Stress in Medicago truncatula Uncover MIEL1 as a Regulator of Seed Germination Plasticity

Legume seeds are an important source of proteins, minerals, and vitamins for human and animal diets and represent a keystone for food security. With climate change and global warming, the production of grain legumes faces new challenges concerning seed vigor traits that allow the fast and homogenous establishment of the crop in a wide range of environments. These seed performance traits are regulated during seed maturation and are under the strong influence of the maternal environment. In this study, we used 200 natural Medicago truncatula accessions, a model species of legumes grown in optimal conditions and under moderate heat stress (26°C) during seed development and maturation. This moderate stress applied at flowering onwards impacted seed weight and germination capacity. Genome-wide association studies (GWAS) were performed to identify putative loci or genes involved in regulating seed traits and their plasticity in response to heat stress. We identified numerous significant quantitative trait nucleotides and potential candidate genes involved in regulating these traits under heat stress by using post-GWAS analyses combined with transcriptomic data. Out of them, MtMIEL1, a RING-type zinc finger family gene, was shown to be highly associated with germination speed in heat-stressed seeds. In Medicago, we highlighted that MtMIEL1 was transcriptionally regulated in heat-stressed seed production and that its expression profile was associated with germination speed in different Medicago accessions. Finally, a loss-of-function analysis of the Arabidopsis MIEL1 ortholog revealed its role as a regulator of germination plasticity of seeds in response to heat stress.

Salinity Stress Alters the Secondary Metabolic Profile of M. sativa, M. arborea and Their Hybrid (Alborea)

Increased soil salinity, and therefore accumulation of ions, is one of the major abiotic stresses of cultivated plants that negatively affect their growth and yield. Among Medicago species, only Medicago truncatula, which is a model plant, has been extensively studied, while research regarding salinity responses of two important forage legumes of Medicago sativa (M. sativa) and Medicago arborea (M. arborea) has been limited. In the present work, differences between M. arborea, M. sativa and their hybrid Alborea were studied regarding growth parameters and metabolomic responses. The entries were subjected to three different treatments: (1) no NaCl application (control plants), (2) continuous application of 100 mM NaCl (acute stress) and (3) gradual application of NaCl at concentrations of 50-75-150 mM by increasing NaCl concentration every 10 days. According to the results, M. arborea maintained steady growth in all three treatments and appeared to be more resistant to salinity. Furthermore, results clearly demonstrated that M. arborea presented a different metabolic profile from that of M. sativa and their hybrid. In general, it was found that under acute and gradual stress, M. sativa overexpressed saponins in the shoots while M. arborea overexpressed saponins in the roots, which is the part of the plant where most of the saponins are produced and overexpressed. Alborea did not perform well, as more metabolites were downregulated than upregulated when subjected to salinity stress. Finally, saponins and hydroxycinnamic acids were key players of increased salinity tolerance.

Genome-Wide Association Studies of Salt Tolerance at Seed Germination and Seedling Stages in Brassica napus

Most crops are sensitive to salt stress, but their degree of susceptibility varies among species and cultivars. In order to understand the salt stress adaptability of Brassica napus to salt stress, we collected the phenotypic data of 505 B. napus accessions at the germination stage under 150 or 215 mM sodium chloride (NaCl) and at the seedling stage under 215 mM NaCl. Genome-wide association studies (GWAS) of 16 salt tolerance coefficients (STCs) were applied to investigate the genetic basis of salt stress tolerance of B. napus. In this study, we mapped 31 salts stress-related QTLs and identified 177 and 228 candidate genes related to salt stress tolerance were detected at germination and seedling stages, respectively. Overexpression of two candidate genes, BnCKX5 and BnERF3 overexpression, were found to increase the sensitivity to salt and mannitol stresses at the germination stage. This study demonstrated that it is a feasible method to dissect the genetic basis of salt stress tolerance at germination and seedling stages in B. napus by GWAS, which provides valuable loci for improving the salt stress tolerance of B. napus. Moreover, these candidate genes are rich genetic resources for the following exploration of molecular mechanisms in adaptation to salt stress in B. napus.

Genome-Wide Association Study Reveals Complex Genetic Architecture of Cadmium and Mercury Accumulation and Tolerance Traits in Medicago truncatula

Heavy metals are an increasing problem due to contamination from human sources that and can enter the food chain by being taken up by plants. Understanding the genetic basis of accumulation and tolerance in plants is important for reducing the uptake of toxic metals in crops and crop relatives, as well as for removing heavy metals from soils by means of phytoremediation. Following exposure of Medicago truncatula seedlings to cadmium (Cd) and mercury (Hg), we conducted a genome-wide association study using relative root growth (RRG) and leaf accumulation measurements. Cd and Hg accumulation and RRG had heritability ranging 0.44 - 0.72 indicating high genetic diversity for these traits. The Cd and Hg trait associations were broadly distributed throughout the genome, indicated the traits are polygenic and involve several quantitative loci. For all traits, candidate genes included several membrane associated ATP-binding cassette transporters, P-type ATPase transporters, oxidative stress response genes, and stress related UDP-glycosyltransferases. The P-type ATPase transporters and ATP-binding cassette protein-families have roles in vacuole transport of heavy metals, and our findings support their wide use in physiological plant responses to heavy metals and abiotic stresses. We also found associations between Cd RRG with the genes CAX3 and PDR3, two linked adjacent genes, and leaf accumulation of Hg associated with the genes NRAMP6 and CAX9. When plant genotypes with the most extreme phenotypes were compared, we found significant divergence in genomic regions using population genomics methods that contained metal transport and stress response gene ontologies. Several of these genomic regions show high linkage disequilibrium (LD) among candidate genes suggesting they have evolved together. Minor allele frequency (MAF) and effect size of the most significant SNPs was negatively correlated with large effect alleles being most rare. This is consistent with purifying selection against alleles that increase toxicity and abiotic stress. Conversely, the alleles with large affect that had higher frequencies that were associated with the exclusion of Cd and Hg. Overall, macroevolutionary conservation of heavy metal and stress response genes is important for improvement of forage crops by harnessing wild genetic variants in gene banks such as the Medicago HapMap collection.

Genome-wide association study in accessions of the mini-core collection of mungbean (Vigna radiata) from the World Vegetable Gene Bank (Taiwan)

BACKGROUND

Mungbean (Vigna radiata (L.) R. Wilczek, or green gram) is important tropical and sub-tropical legume and a rich source of dietary protein and micronutrients. In this study we employ GWAS to examine the genetic basis of variation in several important traits in mungbean, using the mini-core collection established by the World Vegetable Center, which includes 296 accessions that represent the major market classes. This collection has been grown in a common field plot in southern European part of Russia in 2018.

RESULTS

We used 5041 SNPs in 293 accessions that passed strict filtering for genetic diversity, linkage disequilibrium, population structure and GWAS analysis. Polymorphisms were distributed among all chromosomes, but with variable density. Linkage disequilibrium decayed in approximately 105 kb. Four distinct subgroups were identified within 293 accessions with 70% of accessions attributed to one of the four populations. By performing GWAS on the mini-core collection we have found several loci significantly associated with two important agronomical traits. Four SNPs associated with possibility of maturation in Kuban territory of Southern Russia in 2018 were identified within a region of strong linkage which contains genes encoding zinc finger A20 and an AN1 domain stress-associated protein.

CONCLUSIONS

The core collection of mungbean established by the World Vegetable Center is a valuable resource for mungbean breeding. The collection has been grown in southern European part of Russia in 2018 under incidental stresses caused by abnormally hot weather and different photoperiod. We have found several loci significantly associated with color of hypocotyl and possibility of maturation under these stressful conditions. SNPs associated with possibility of maturation localize to a region on chromosome 2 with strong linkage, in which genes encoding zinc finger A20 and AN1 domain stress associated protein (SAP) are located. Phenotyping of WorldVeg collection for maturation traits in temperate climatic locations is important as phenology remains a critical breeding target for mungbean. As demand rises for mungbean, production in temperate regions with shorter growing seasons becomes crucial to keep up with needs. Uncovering SNPs for phenology traits will speed breeding efforts.

Review: Climate change impacts on food security- focus on perennial cropping systems and nutritional value

Anthropogenic increases in fossil fuel emissions have been a primary driver of increased concentrations of atmospheric carbon dioxide ([CO₂]) and other greenhouse gases resulting in warmer temperatures, alterations in precipitation patterns, and increased occurrence of extreme weather events in terrestrial areas across the globe. In agricultural growing regions, alterations in climate can challenge plant productivity in ways that impact the ability of the world to sustain adequate food production for a growing and increasingly affluent population with shifting access to affordable and nutritious food. While the knowledge gap that exists regarding potential climate change impacts is large across agriculture, it is especially large in specialty cropping systems. This includes fruit and vegetable crops, and perennial cropping systems which also contribute (along with row crops) to our global diet. In order to obtain a comprehensive view of the true impact of climate change on our global food supply, we must expand our narrow focus from improving yield and plant productivity to include the impact of climate change on the nutritional value of these crops. In order to address these questions, we need a multi-faceted approach that integrates physiology and genomics tools and conducts comprehensive experiments under realistic depictions of future projected climate. This review describes gaps in our knowledge in relation to these responses, and future questions and actions that are needed to develop a sustainable future food supply in light of global climate change.

Loss of salt tolerance during tomato domestication conferred by variation in a Na+ /K+ transporter

Domestication has resulted in reduced salt tolerance in tomato. To identify the genetic components causing this deficiency, we performed a genome-wide association study (GWAS) for root Na⁺ /K⁺ ratio in a population consisting of 369 tomato accessions with large natural variations. The most significant variations associated with root Na⁺ /K⁺ ratio were identified within the gene SlHAK20 encoding a member of the clade IV HAK/KUP/KT transporters. We further found that SlHAK20 transports Na⁺ and K⁺ and regulates Na⁺ and K⁺ homeostasis under salt stress conditions. A variation in the coding sequence of SlHAK20 was found to be the causative variant associated with Na⁺ /K⁺ ratio and confer salt tolerance in tomato. Knockout mutations in tomato SlHAK20 and the rice homologous genes resulted in hypersensitivity to salt stress. Together, our study uncovered a previously unknown molecular mechanism of salt tolerance responsible for the deficiency in salt tolerance in cultivated tomato varieties. Our findings provide critical information for molecular breeding to improve salt tolerance in tomato and other crops.

Method for Genome-Wide Association Study: A Soybean Example

Genome-wide association studies (GWAS) are a valuable approach to identify single nucleotide polymorphisms (SNPs) associated with a phenotype of interest. There are now a variety of R-packages and command line tools available to perform GWAS. Here, we provide an example downloading and filtering SNP data, followed by GWAS analysis using the R-package rMVP.

Genomics of Plant Disease Resistance in Legumes

The constant interactions between plants and pathogens in the environment and the resulting outcomes are of significant importance for agriculture and agricultural scientists. Disease resistance genes in plant cultivars can break down in the field due to the evolution of pathogens under high selection pressure. Thus, the protection of crop plants against pathogens is a continuous arms race. Like any other type of crop plant, legumes are susceptible to many pathogens. The dawn of the genomic era, in which high-throughput and cost-effective genomic tools have become available, has revolutionized our understanding of the complex interactions between legumes and pathogens. Genomic tools have enabled a global view of transcriptome changes during these interactions, from which several key players in both the resistant and susceptible interactions have been identified. This review summarizes some of the large-scale genomic studies that have clarified the host transcriptional changes during interactions between legumes and their plant pathogens while highlighting some of the molecular breeding tools that are available to introgress the traits into breeding programs. These studies provide valuable insights into the molecular basis of different levels of host defenses in resistant and susceptible interactions.