Genome-wide association of drought-related and biomass traits with HapMap SNPs in Medicago truncatula

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.

The Genetic Dissection of Nitrogen Use-Related Traits in Flax (Linum usitatissimum L.) at the Seedling Stage through the Integration of Multi-Locus GWAS, RNA-seq and Genomic Selection

Nitrogen (N), the most important macro-nutrient for plant growth and development, is a key factor that determines crop yield. Yet its excessive applications pollute the environment and are expensive. Hence, studying nitrogen use efficiency (NUE) in crops is fundamental for sustainable agriculture. Here, an association panel consisting of 123 flax accessions was evaluated for 21 NUE-related traits at the seedling stage under optimum N (N+) and N deficiency (N-) treatments to dissect the genetic architecture of NUE-related traits using a multi-omics approach integrating genome-wide association studies (GWAS), transcriptome analysis and genomic selection (GS). Root traits exhibited significant and positive correlations with NUE under N- conditions (r = 0.33 to 0.43, p < 0.05). A total of 359 QTLs were identified, accounting for 0.11% to 23.1% of the phenotypic variation in NUE-related traits. Transcriptomic analysis identified 1034 differentially expressed genes (DEGs) under contrasting N conditions. DEGs involved in N metabolism, root development, amino acid transport and catabolism and others, were found near the QTLs. GS models to predict NUE stress tolerance index (NUE_STI) trait were tested using a random genome-wide SNP dataset and a GWAS-derived QTLs dataset. The latter produced superior prediction accuracy (r = 0.62 to 0.79) compared to the genome-wide SNP marker dataset (r = 0.11) for NUE_STI. Our results provide insights into the QTL architecture of NUE-related traits, identify candidate genes for further studies, and propose genomic breeding tools to achieve superior NUE in flax under low N input.

Genetic Architecture of Salt Tolerance in Cowpea (Vigna unguiculata (L.) Walp.) at Seedling Stage Using a Whole Genome Resequencing Approach

Cowpea (Vigna unguiculata (L.) Walp.) is a diploid legume crop used for human consumption, feed for livestock, and cover crops. Earlier reports have shown that salinity has been a growing threat to cowpea cultivation. The objectives of this study were to conduct a genome-wide association study (GWAS) to identify SNP markers and to investigate candidate genes for salt tolerance in cowpea. A total of 331 cowpea genotypes were evaluated for salt tolerance by supplying a solution of 200 mM NaCl in our previous work. The cowpea panel was genotyped using a whole genome resequencing approach, generating 14,465,516 SNPs. Moreover, 5,884,299 SNPs were used after SNP filtering. GWAS was conducted on a total of 296 cowpea genotypes that have high-quality SNPs. BLINK was used for conducting GWAS. Results showed (1) a strong GWAS peak on an 890-bk region of chromosome 2 for leaf SPAD chlorophyll under salt stress in cowpea and harboring a significant cluster of nicotinamide adenine dinucleotide (NAD) dependent epimerase/dehydratase genes such as Vigun02g128900.1, Vigun02g129000.1, Vigun02g129100.1, Vigun02g129200.1, and Vigun02g129500.1; (2) two GWAS peaks associated with relative tolerance index for chlorophyll were identified on chromosomes 1 and 2. The peak on chromosome 1 was defined by a cluster of 10 significant SNPs mapped on a 5 kb region and was located in the vicinity of Vigun01g086000.1, encoding for a GATA transcription factor. The GWAS peak on chromosome 2 was defined by a cluster of 53 significant SNPs and mapped on a 68 bk region of chromosome 2, and (3) the highest GWAS peak was identified on chromosome 3, and this locus was associated with leaf score injury. This peak was within the structure of a potassium channel gene (Vigun03g144700.1). To the best of our knowledge, this is one the earliest reports on the salt tolerance study of cowpea using whole genome resequencing data.

Genome-wide association study of drought tolerance traits in sugar beet germplasms at the seedling stage

Introduction: Sugar beets are an important crop for global sugar production. Intense drought and the increasing lack of water resources pose a great threat to sugar beet cultivation. It is a priority to investigate favourable germplasms and functional genes to improve the breeding of drought tolerant plants. Methods: Thus, in this study, 328 sugar beet germplasms were used in a genome-wide association study (GWAS) to identify single nucleotide polymorphism (SNP) markers and candidate genes associated with drought tolerance. Results: The results showed that under drought stress (9% PEG-6000), there were 11 significantly associated loci on chromosomes 2, 3, 5, 7, and 9 from the 108946 SNPs filtered using a mixed linear model (MLM). Genome-wide association analysis combined with qRT-PCR identified 13 genes that were significantly differentially expressed in drought-tolerant extreme materials. Discussion: These candidate genes mainly exhibited functions such as regulating sugar metabolism, maintaining internal environmental stability and participating in photosystem repair. This study provides valuable information for exploring the molecular mechanisms of drought tolerance and improvement in sugar beet.

Evolutionary dynamics and adaptive benefits of deleterious mutations in crop gene pools

Mutations with deleterious consequences in nature may be conditionally deleterious in crop plants. That is, while some genetic variants may reduce fitness under wild conditions and be subject to purifying selection, they can be under positive selection in domesticates. Such deleterious alleles can be plant breeding targets, particularly for complex traits. The difficulty of distinguishing favorable from unfavorable variants reduces the power of selection, while favorable trait variation and heterosis may be attributable to deleterious alleles. Here, we review the roles of deleterious mutations in crop breeding and discuss how they can be used as a new avenue for crop improvement with emerging genomic tools, including HapMaps and pangenome analysis, aiding the identification, removal, or exploitation of deleterious mutations.

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.

Outlier analyses and genome-wide association study identify glgC and ERD6-like 4 as candidate genes for foliar water-soluble carbohydrate accumulation in Trifolium repens

Increasing water-soluble carbohydrate (WSC) content in white clover is important for improving nutritional quality and reducing environmental impacts from pastoral agriculture. Elucidation of genes responsible for foliar WSC variation would enhance genetic improvement by enabling molecular breeding approaches. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) associated with variation in foliar WSC in white clover. A set of 935 white clover individuals, randomly sampled from five breeding pools selectively bred for divergent (low or high) WSC content, were assessed with 14,743 genotyping-by-sequencing SNPs, using three outlier detection methods: PCAdapt, BayeScan and KGD-F_ST. These analyses identified 33 SNPs as discriminating between high and low WSC populations and putatively under selection. One SNP was located in the intron of ERD6-like 4, a gene coding for a sugar transporter located on the vacuole membrane. A genome-wide association study using a subset of 605 white clover individuals and 5,757 SNPs, identified a further 12 SNPs, one of which was associated with a starch biosynthesis gene, glucose-1-phosphate adenylyltransferase, glgC. Our results provide insight into genomic regions underlying WSC accumulation in white clover, identify candidate genomic regions for further functional validation studies, and reveal valuable information for marker-assisted or genomic selection in white clover.

Plastid phylogenomics uncovers multiple species in Medicago truncatula (Fabaceae) germplasm accessions

Medicago truncatula is a model legume that has been extensively investigated in diverse subdisciplines of plant science. Medicago littoralis can interbreed with M. truncatula and M. italica; these three closely related species form a clade, i.e. TLI clade. Genetic studies have indicated that M. truncatula accessions are heterogeneous but their taxonomic identities have not been verified. To elucidate the phylogenetic position of diverse M. truncatula accessions within the genus, we assembled 54 plastid genomes (plastomes) using publicly available next-generation sequencing data and conducted phylogenetic analyses using maximum likelihood. Five accessions showed high levels of plastid DNA polymorphism. Three of these highly polymorphic accessions contained sequences from both M. truncatula and M. littoralis. Phylogenetic analyses of sequences placed some accessions closer to distantly related species suggesting misidentification of source material. Most accessions were placed within the TLI clade and maximally supported the interrelationships of three subclades. Two Medicago accessions were placed within a M. italica subclade of the TLI clade. Plastomes with a 45-kb (rpl20-ycf1) inversion were placed within the M. littoralis subclade. Our results suggest that the M. truncatula accession genome pool represents more than one species due to possible mistaken identities and gene flow among closely related species.

Differential environmental and genomic architectures shape the natural diversity for trichome patterning and morphology in different Arabidopsis organs

Despite the adaptive and taxonomic relevance of the natural diversity for trichome patterning and morphology, the molecular and evolutionary mechanisms underlying these traits remain mostly unknown, particularly in organs other than leaves. In this study, we address the ecological, genetic and molecular bases of the natural variation for trichome patterning and branching in multiple organs of Arabidopsis (Arabidopsis thaliana). To this end, we characterized a collection of 191 accessions and carried out environmental and genome-wide association (GWA) analyses. Trichome amount in different organs correlated negatively with precipitation in distinct seasons, thus suggesting a precise fit between trichome patterning and climate throughout the Arabidopsis life cycle. In addition, GWA analyses showed small overlapping between the genes associated with different organs, indicating partly independent genetic bases for vegetative and reproductive phases. These analyses identified a complex locus on chromosome 2, where two adjacent MYB genes (ETC2 and TCL1) displayed differential effects on trichome patterning in several organs. Furthermore, analyses of transgenic lines carrying different natural alleles demonstrated that TCL1 accounts for the variation for trichome patterning in all organs, and for stem trichome branching. By contrast, two other MYB genes (TRY and GL1), mainly showed effects on trichome patterning or branching, respectively.

Natural variation identifies new effectors of water-use efficiency in Arabidopsis

Water-use efficiency (WUE) is the ratio of biomass produced per unit of water consumed; thus, it can be altered by genetic factors that affect either side of the ratio. In the present study, we exploited natural variation for WUE to discover loci affecting either biomass accumulation or water use as factors affecting WUE. Genome-wide association studies (GWAS) using integrated WUE measured through carbon isotope discrimination (δ¹³C) of Arabidopsis thaliana accessions identified genomic regions associated with WUE. Reverse genetic analysis of 70 candidate genes selected based on the GWAS results and transcriptome data identified 25 genes affecting WUE as measured by gravimetric and δ¹³C analyses. Mutants of four genes had higher WUE than wild type, while mutants of the other 21 genes had lower WUE. The differences in WUE were caused by either altered biomass or water consumption (or both). Stomatal density (SD) was not a primary cause of altered WUE in these mutants. Leaf surface temperatures indicated that transpiration differed for mutants of 16 genes, but generally biomass accumulation had a greater effect on WUE. The genes we identified are involved in diverse cellular processes, including hormone and calcium signaling, meristematic activity, photosynthesis, flowering time, leaf/vasculature development, and cell wall composition; however, none of them had been previously linked to WUE. Thus, our study successfully identified effectors of WUE that can be used to understand the genetic basis of WUE and improve crop productivity.

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.

Genome-wide identification of invertases in Fabaceae, focusing on transcriptional regulation of Pisum sativum invertases in seed subjected to drought

Invertases are key enzymes for carbon metabolism, cleaving sucrose into energy-rich and signaling metabolites, glucose and fructose. Invertases play pivotal roles in development and stress response, determining yield and quality of seed production. In this context, the repertoire of invertase gene families is critically scarce in legumes. Here, we performed a systematic search for invertase families in 16 Fabaceae genomes. For instance, we identified 19 invertase genes in the model plant Medicago and 17 accessions in the agronomic crop Pisum sativum. Our comprehensive phylogenetic analysis sets a milestone for the scientific community as we propose a new nomenclature to correctly name plant invertases. Thus, neutral invertases were classified into four clades of cytosolic invertase (CINV). Acid invertases were classified into two cell wall invertase clades (CWINV) and two vacuolar invertase clades (VINV). Then, we explored transcriptional regulation of the pea invertase family, focusing on seed development and water stress. Invertase expression decreased sharply from embryogenesis to seed-filling stages, consistent with higher sucrose and lower monosaccharide contents. The vacuolar invertase PsVINV1.1 clearly marked the transition between both developmental stages. We hypothesize that the predominantly expressed cell wall invertase, PsCWINV1.2, may drive sucrose unloading towards developing seeds. The same candidates, PsVINV1.1 and PsCWINV1.2, were also regulated by water deficit during embryonic stage. We suggest that PsVINV1.1 along with vacuolar sugar transporters maintain cellular osmotic pressure and PsCWINV1.2 control hexose provision, thereby ensuring embryo survival in drought conditions. Altogether, our findings provide novel insights into the regulation of plant carbon metabolism in a challenging environment.

Physiological and transcriptome analyses for assessing the effects of exogenous uniconazole on drought tolerance in hemp (Cannabis sativa L.)

Uniconazole (S-(+)-uniconazole), a plant growth retardant, exerts key roles in modulating growth and development and increasing abiotic stress tolerance in plants. However, the underlying mechanisms by which uniconazole regulates drought response remain largely unknown. Here, the effects of exogenous uniconazole on drought tolerance in hemp were studied via physiological and transcriptome analyses of the drought-sensitive industrial hemp cultivar Hanma No. 2 grown under drought stress. Exogenous uniconazole treatment increased hemp tolerance to drought-induced damage by enhancing chlorophyll content and photosynthesis capacity, regulating activities of enzymes involved in carbon and nitrogen metabolism, and altering endogenous hormone levels. Expression of genes associated with porphyrin and chlorophyll metabolism, photosynthesis-antenna proteins, photosynthesis, starch and sucrose metabolism, nitrogen metabolism, and plant hormone signal transduction were significantly regulated by uniconazole compared with that by control (distilled water) under drought stress. Numerous genes were differentially expressed to increase chlorophyll content, enhance photosynthesis, regulate carbon-nitrogen metabolism-related enzyme activities, and alter endogenous hormone levels. Thus, uniconazole regulated physiological and molecular characteristics of photosynthesis, carbon-nitrogen metabolism, and plant hormone signal transduction to enhance drought resistance in industrial hemp.

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

Distinct genetic basis for root responses to lipo-chitooligosaccharide signal molecules from different microbial origins

Lipo-chitooligosaccharides (LCOs) were originally found as symbiotic signals called Nod Factors (Nod-LCOs) controlling the nodulation of legumes by rhizobia. More recently, LCOs were also found in symbiotic fungi and, more surprisingly, very widely in the kingdom Fungi, including in saprophytic and pathogenic fungi. The LCO-V(C18:1, fucosylated/methyl fucosylated), hereafter called Fung-LCOs, are the LCO structures most commonly found in fungi. This raises the question of how legume plants such as Medicago truncatula can discriminate between Nod-LCOs and Fung-LCOs. To address this question, we performed a genome-wide association study on 173 natural accessions of M. truncatula, using a root branching phenotype and a newly developed local score approach. Both Nod-LCOs and Fung-LCOs stimulated root branching in most accessions, but the root responses to these two types of LCO molecules were not correlated. In addition, the heritability of the root response was higher for Nod-LCOs than for Fung-LCOs. We identified 123 loci for Nod-LCO and 71 for Fung-LCO responses, of which only one was common. This suggests that Nod-LCOs and Fung-LCOs both control root branching but use different molecular mechanisms. The tighter genetic constraint of the root response to Fung-LCOs possibly reflects the ancestral origin of the biological activity of these molecules.

Genomic prediction for growth using a low-density SNP panel in dromedary camels

For thousands of years, camels have produced meat, milk, and fiber in harsh desert conditions. For a sustainable development to provide protein resources from desert areas, it is necessary to pay attention to genetic improvement in camel breeding. By using genotyping-by-sequencing (GBS) method we produced over 14,500 genome wide markers to conduct a genome- wide association study (GWAS) for investigating the birth weight, daily gain, and body weight of 96 dromedaries in the Iranian central desert. A total of 99 SNPs were associated with birth weight, daily gain, and body weight (p-value < 0.002). Genomic breeding values (GEBVs) were estimated with the BGLR package using (i) all 14,522 SNPs and (ii) the 99 SNPs by GWAS. Twenty-eight SNPs were associated with birth weight, daily gain, and body weight (p-value < 0.001). Annotation of the genomic region (s) within ± 100 kb of the associated SNPs facilitated prediction of 36 candidate genes. The accuracy of GEBVs was more than 0.65 based on all 14,522 SNPs, but the regression coefficients for birth weight, daily gain, and body weight were 0.39, 0.20, and 0.23, respectively. Because of low sample size, the GEBVs were predicted using the associated SNPs from GWAS. The accuracy of GEBVs based on the 99 associated SNPs was 0.62, 0.82, and 0.57 for birth weight, daily gain, and body weight. This report is the first GWAS using GBS on dromedary camels and identifies markers associated with growth traits that could help to plan breeding program to genetic improvement. Further researches using larger sample size and collaboration of the camel farmers and more profound understanding will permit verification of the associated SNPs identified in this project. The preliminary results of study show that genomic selection could be the appropriate way to genetic improvement of body weight in dromedary camels, which is challenging due to a long generation interval, seasonal reproduction, and lack of records and pedigrees.

Genome-Wide Association Study of Maize Aboveground Dry Matter Accumulation at Seedling Stage

Dry matter accumulation and partitioning during the early phases of development could significantly affect crop growth and productivity. In this study, the aboveground dry matter (DM), the DM of different organs, and partition coefficients of a maize association mapping panel of 412 inbred lines were evaluated at the third and sixth leaf stages (V3 and V6). Further, the properties of these phenotypic traits were analyzed. Genome-wide association studies (GWAS) were conducted on the total aboveground biomass and the DM of different organs. Analysis of GWAS results identified a total of 1,103 unique candidate genes annotated by 678 significant SNPs (P value < 1.28e-6). A total of 224 genes annotated by SNPs at the top five of each GWAS method and detected by multiple GWAS methods were regarded as having high reliability. Pathway enrichment analysis was also performed to explore the biological significance and functions of these candidate genes. Several biological pathways related to the regulation of seed growth, gibberellin-mediated signaling pathway, and long-day photoperiodism were enriched. The results of our study could provide new perspectives on breeding high-yielding maize varieties.

A linkage disequilibrium-based statistical test for Genome-Wide Epistatic Selection Scans in structured populations

The quest for signatures of selection using single nucleotide polymorphism (SNP) data has proven efficient to uncover genes involved in conserved and/or adaptive molecular functions, but none of the statistical methods were designed to identify interacting alleles as targets of selective processes. Here, we propose a statistical test aimed at detecting epistatic selection, based on a linkage disequilibrium (LD) measure accounting for population structure and heterogeneous relatedness between individuals. SNP-based ([Formula: see text]) and window-based ([Formula: see text]) statistics fit a Student distribution, allowing to test the significance of correlation coefficients. As a proof of concept, we use SNP data from the Medicago truncatula symbiotic legume plant and uncover a previously unknown gene coadaptation between the MtSUNN (Super Numeric Nodule) receptor and the MtCLE02 (CLAVATA3-Like) signaling peptide. We also provide experimental evidence supporting a MtSUNN-dependent negative role of MtCLE02 in symbiotic root nodulation. Using human HGDP-CEPH SNP data, our new statistical test uncovers strong LD between SLC24A5 (skin pigmentation) and EDAR (hairs, teeth, sweat glands development) world-wide, which persists after correction for population structure and relatedness in Central South Asian populations. This result suggests that epistatic selection or coselection could have contributed to the phenotypic make-up in some human populations. Applying this approach to genome-wide SNP data will facilitate the identification of coadapted gene networks in model or non-model organisms.

GWASpro: a high-performance genome-wide association analysis server

Summary

We present GWASpro, a high-performance web server for the analyses of large-scale genome-wide association studies (GWAS). GWASpro was developed to provide data analyses for large-scale molecular genetic data, coupled with complex replicated experimental designs such as found in plant science investigations and to overcome the steep learning curves of existing GWAS software tools. GWASpro supports building complex design matrices, by which complex experimental designs that may include replications, treatments, locations and times, can be accounted for in the linear mixed model. GWASpro is optimized to handle GWAS data that may consist of up to 10 million markers and 10 000 samples from replicable lines or hybrids. GWASpro provides an interface that significantly reduces the learning curve for new GWAS investigators.

Availability and implementation

GWASpro is freely available at https://bioinfo.noble.org/GWASPRO.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Physical Dormancy Release in Medicago truncatula Seeds Is Related to Environmental Variations

Seed dormancy and timing of its release is an important developmental transition determining the survival of individuals, populations, and species in variable environments. Medicago truncatula was used as a model to study physical seed dormancy at the ecological and genetics level. The effect of alternating temperatures, as one of the causes releasing physical seed dormancy, was tested in 178 M. truncatula accessions over three years. Several coefficients of dormancy release were related to environmental variables. Dormancy varied greatly (4-100%) across accessions as well as year of experiment. We observed overall higher physical dormancy release under more alternating temperatures (35/15 °C) in comparison with less alternating ones (25/15 °C). Accessions from more arid climates released dormancy under higher experimental temperature alternations more than accessions originating from less arid environments. The plasticity of physical dormancy can probably distribute the germination through the year and act as a bet-hedging strategy in arid environments. On the other hand, a slight increase in physical dormancy was observed in accessions from environments with higher among-season temperature variation. Genome-wide association analysis identified 136 candidate genes related to secondary metabolite synthesis, hormone regulation, and modification of the cell wall. The activity of these genes might mediate seed coat permeability and, ultimately, imbibition and germination.

A Genome-Wide Association Study Identifies Quantitative Trait Loci Affecting Hematological Traits in Camelus bactrianus

Simple Summary

Bactrian camels can adapt to harsh natural environments. This unique tolerance of camels is tightly linked to their hematological traits, which are related to their immune, metabolic, and disease status. Therefore, mapping genomic regions that affect blood cell traits can help identify genomic characteristics that can be used as biomarkers of immune, metabolic, and disease states. This knowledge will further our understanding of the camel’s tolerance mechanisms.

Abstract

Bactrian camels (Camelus bactrianus) are one of the few large livestock species that can survive in the Gobi Desert. Animal immunity and disease resistance are related to hematological traits, which are also associated with tolerance observed in Bactrian camels. However, no genome-wide association studies have examined the genetic mechanism of the immune capability of Bactrian camels. In the present study, we used genotyping-by-sequencing data generated from 366 Bactrian camel accessions to perform a genome-wide association study for 17 hematological traits. Of the 256,616 single-nucleotide polymorphisms (SNPs) obtained, 1,635 trait–SNP associations were among the top quantitative trait locus candidates. Lastly, 664 candidate genes associated with 13 blood traits were identified. The most significant were ZNF772, MTX2, ESRRG, MEI4, IL11, FRMPD4, GABPA, NTF4, CRYBG3, ENPP5, COL16A1, and CD207. The results of our genome-wide association study provide a list of significant SNPs and candidate genes, which offer valuable information for further dissection of the molecular mechanisms that regulate the camel’s hematological traits to ultimately reveal their tolerance mechanisms.

LeGOO: An Expertized Knowledge Database for the Model Legume Medicago truncatula

Medicago truncatula was proposed, about three decades ago, as a model legume to study the Rhizobium-legume symbiosis. It has now been adopted to study a wide range of biological questions, including various developmental processes (in particular root, symbiotic nodule and seed development), symbiotic (nitrogen-fixing and arbuscular mycorrhizal endosymbioses) and pathogenic interactions, as well as responses to abiotic stress. With a number of tools and resources set up in M. truncatula for omics, genetics and reverse genetics approaches, massive amounts of data have been produced, as well as four genome sequence releases. Many of these data were generated with heterogeneous tools, notably for transcriptomics studies, and are consequently difficult to integrate. This issue is addressed by the LeGOO (for Legume Graph-Oriented Organizer) knowledge base (https://www.legoo.org), which finds the correspondence between the multiple identifiers of the same gene. Furthermore, an important goal of LeGOO is to collect and represent biological information from peer-reviewed publications, whatever the technical approaches used to obtain this information. The information is modeled in a graph-oriented database, which enables flexible representation, with currently over 200,000 relations retrieved from 298 publications. LeGOO also provides the user with mining tools, including links to the Mt5.0 genome browser and associated information (on gene functional annotation, expression, methylome, natural diversity and available insertion mutants), as well as tools to navigate through different model species. LeGOO is, therefore, an innovative database that will be useful to the Medicago and legume community to better exploit the wealth of data produced on this model species.

Dissection of Phenotypic and Genetic Variation of Drought-Related Traits in Diverse Chinese Wheat Landraces

Variations in 16 seedling traits under normal and drought conditions were investigated. Extremely resistant and sensitive accessions were identified for future analyses. Under normal and drought conditions, 57 and 29 QTL were identified, respectively. A total of 77 candidate genes were identified, and four were validated by qRT-PCR. Drought is one of the most important abiotic stressors affecting wheat (Triticum aestivum L.) production. To improve wheat yield, a better understanding of the genetic control of traits governing drought resistance is paramount. Here, using 645 wheat landraces, we evaluated 16 seedling traits related to root and shoot growth and water content under normal and drought (induced by polyethylene glycol) conditions. Extremely resistant and sensitive accessions were identified for future drought-resistance breeding and further genetic analyses. A genome-wide association study was performed for the 16 traits using 52,118 diversity arrays technology sequencing (DArT-seq) markers. A total of 57 quantitative trait loci (QTL) were detected for seven traits under normal conditions, whereas 29 QTL were detected for eight traits under drought conditions. On the basis of these markers, we identified 56 candidate genes for six seedling traits under normal conditions, and 21 candidate genes for seven seedling traits under drought conditions. Four candidate genes were validated under normal and drought conditions using quantitative reverse transcription polymerase chain reaction (qRT-PCR) data. The co-localization of the flowering date and drought-related traits indicates that the regulatory networks of flowering may also respond to drought stress or are associated with the correlated responses of these traits. The phenotypic and genetic elucidation of drought-related traits will assist future gene discovery efforts and provide a basis for breeding drought-resistant wheat cultivars.

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.

The genetic basis of drought tolerance in the high oil crop Sesamum indicum

Unlike most of the important food crops, sesame can survive drought but severe and repeated drought episodes, especially occurring during the reproductive stage, significantly curtail the productivity of this high oil crop. Genome-wide association study was conducted for traits related to drought tolerance using 400 diverse sesame accessions, including landraces and modern cultivars. Ten stable QTLs explaining more than 40% of the phenotypic variation and located on four linkage groups were significantly associated with drought tolerance related traits. Accessions from the tropical area harboured higher numbers of drought tolerance alleles at the peak loci and were found to be more tolerant than those from the northern area, indicating a long-term genetic adaptation to drought-prone environments. We found that sesame has already fixed important alleles conferring survival to drought which may explain its relative high drought tolerance. However, most of the alleles crucial for productivity and yield maintenance under drought conditions are far from been fixed. This study also revealed that pyramiding the favourable alleles observed at the peak loci is of high potential for enhancing drought tolerance in sesame. In addition, our results highlighted two important pleiotropic QTLs harbouring known and unreported drought tolerance genes such as SiABI4, SiTTM3, SiGOLS1, SiNIMIN1 and SiSAM. By integrating candidate gene association study, gene expression and transgenic experiments, we demonstrated that SiSAM confers drought tolerance by modulating polyamine levels and ROS homeostasis, and a missense mutation in the coding region partly contributes to the natural variation of drought tolerance in sesame.

Analysis of Drought Tolerance and Associated Traits in Upland Cotton at the Seedling Stage

(1) Background: Upland cotton (Gossypium hirsutum L.) is the most important natural fiber worldwide, and it is extensively planted and plentifully used in the textile industry. Major cotton planting regions are frequently affected by abiotic stress, especially drought stress. Drought resistance is a complex, quantitative trait. A genome-wide association study (GWAS) constitutes an efficient method for dissecting the genetic architecture of complex traits. In this study, the drought resistance of a population of 316 upland cotton accessions was studied via GWAS. (2) Methods: GWAS methodology was employed to identify relationships between molecular markers or candidate genes and phenotypes of interest. (3) Results: A total of 8, 3, and 6 SNPs were associated with the euphylla wilting score (EWS), cotyledon wilting score (CWS), and leaf temperature (LT), respectively, based on a general linear model and a factored spectrally transformed linear mixed model. For these traits, 7 QTLs were found, of which 2 each were located on chromosomes A05, A11, and D03, and of which 1 was located on chromosome A01. Importantly, in the candidate regions WRKY70, GhCIPK6, SnRK2.6, and NET1A, which are involved in the response to abscisic acid (ABA), the mitogen-activated protein kinase (MAPK) signaling pathway and the calcium transduction pathway were identified in upland cotton at the seedling stage under drought stress according to annotation information and linkage disequilibrium (LD) block analysis. Moreover, RNA sequencing analysis showed that WRKY70, GhCIPK6, SnRK2.6, and NET1A were induced by drought stress, and the expression of these genes was significantly different between normal and drought stress conditions. (4) Conclusions: The present study should provide some genomic resources for drought resistance in upland cotton. Moreover, the germplasm of the different phenotypes, the detected SNPs and, the potential candidate genes will be helpful for molecular marker-assisted breeding studies about increased drought resistance in upland cotton.

Drought stress stimulates endocytosis and modifies membrane lipid order of rhizodermal cells of Medicago truncatula in a genotype-dependent manner

BACKGROUND

Drought stress negatively affects plant growth and productivity. Plants sense soil drought at the root level but the underlying mechanisms remain unclear. At the cell level, we aim to reveal the short-term root perception of drought stress through membrane dynamics.

RESULTS

In our study, 15 Medicago truncatula accessions were exposed to a polyethylene glycol (PEG)-induced drought stress, leading to contrasted ecophysiological responses, in particular related to root architecture plasticity. In the reference accession Jemalong A17, identified as drought susceptible, we analyzed lateral roots by imaging of membrane-localized fluorescent probes using confocal microscopy. We found that PEG stimulated endocytosis especially in cells belonging to the growth differentiation zone (GDZ). The mapping of membrane lipid order in cells along the root apex showed that membranes of root cap cells were more ordered than those of more differentiated cells. Moreover, PEG triggered a significant increase in membrane lipid order of rhizodermal cells from the GDZ. We initiated the membrane analysis in the drought resistant accession HM298, which did not reveal such membrane modifications in response to PEG.

CONCLUSIONS

Our data demonstrated that the plasma membranes of root cells from a susceptible genotype perceived drought stress by modulating their physical state both via a stimulation of endocytosis and a modification of the degree of lipid order, which could be proposed as mechanisms required for signal transduction.

A local score approach improves GWAS resolution and detects minor QTL: application to Medicago truncatula quantitative disease resistance to multiple Aphanomyces euteiches isolates

Quantitative trait loci (QTL) with small effects, which are pervasive in quantitative phenotypic variation, are difficult to detect in genome-wide association studies (GWAS). To improve their detection, we propose to use a local score approach that accounts for the surrounding signal due to linkage disequilibrium, by accumulating association signals from contiguous single markers. Simulations revealed that, in a GWAS context with high marker density, the local score approach outperforms single SNP p-value-based tests for detecting minor QTL (heritability of 5-10%) and is competitive with regard to alternative methods, which also aggregate p-values. Using more than five million SNPs, this approach was applied to identify loci involved in Quantitative Disease Resistance (QDR) to different isolates of the plant root rot pathogen Aphanomyces euteiches, from a GWAS performed on a collection of 174 accessions of the model legume Medicago truncatula. We refined the position of a previously reported major locus, underlying MYB/NB-ARC/tyrosine kinase candidate genes conferring resistance to two closely related A. euteiches isolates belonging to pea pathotype I. We also discovered a diversity of minor resistance QTL, not detected using p-value-based tests, some of which being putatively shared in response to pea (pathotype I and III) and/or alfalfa (race 1 and 2) isolates. Candidate genes underlying these QTL suggest pathogen effector recognition and plant proteasome as key functions associated with M. truncatula resistance to A. euteiches. GWAS on any organism can benefit from the local score approach to uncover many weak-effect QTL.

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

Salinity stress is an important cause of crop yield loss in many parts of the world. Here, we performed genome-wide association studies of salinity-stress responsive traits in 132 HapMap genotypes of the model legume Medicago truncatula. Plants grown in soil were subjected to a step-wise increase in NaCl concentration, from 0 through 0.5% and 1.0% to 1.5%, and the following traits were measured: vigor, shoot biomass, shoot water content, leaf chlorophyll content, leaf size, and leaf and root concentrations of proline and major ions (Na+ , Cl- , K+ , Ca2+ , etc.). Genome-wide association studies were carried out using 2.5 million single nucleotide polymorphisms, and 12 genomic regions associated with at least four traits each were identified. Transcript-level analysis of the top eight candidate genes in five extreme genotypes revealed association between salinity tolerance and transcript-level changes for seven of the genes, encoding a vacuolar H+ -ATPase, two transcription factors, two proteins involved in vesicle trafficking, one peroxidase, and a protein of unknown function. Earlier functional studies on putative orthologues of two of the top eight genes (a vacuolar H+ -ATPase and a peroxidase) demonstrated their involvement in plant salinity tolerance.

Functional Genomics and Flowering Time in Medicago truncatula: An Overview

Flowering time is an important trait that influences adaptation and yield in many crop legumes. Both the inherent earliness of flowering and the degree to which it is responsive to environmental factors determine both the eco-geographic range across which crops can be successfully grown and the seasonal cycles most suitable for production. This chapter will provide a brief review of studies investigating the genetic control of flowering time in Medicago truncatula.

Soil environment is a key driver of adaptation in Medicago truncatula: new insights from landscape genomics

Spatial differences in environmental selective pressures interact with the genomes of organisms, ultimately leading to local adaptation. Landscape genomics is an emergent research area that uncovers genome-environment associations, thus allowing researchers to identify candidate loci for adaptation to specific environmental variables. In the present study, we used latent factor mixed models (LFMMs) and Moran spectral outlier detection/randomization (MSOD-MSR) to identify candidate loci for adaptation to 10 environmental variables (climatic, soil and atmospheric) among 43 515 single nucleotide polymorphisms (SNPs) from 202 accessions of the model legume Medicago truncatula. Soil variables were associated with a large number of candidate loci identified through both LFMMs and MSOD-MSR. Genes tagged by candidate loci associated with drought and salinity are involved in the response to biotic and abiotic stresses, while those tagged by candidates associated with soil nitrogen and atmospheric nitrogen, participate in the legume-rhizobia symbiosis. Candidate SNPs identified through both LFMMs and MSOD-MSR explained up to 56% of variance in flowering traits. Our findings highlight the importance of soil in driving adaptation in the system and elucidate the basis of evolutionary potential of M. truncatula to respond to global climate change and anthropogenic disruption of the nitrogen cycle.

Functional Genomics and Seed Development in Medicago truncatula: An Overview

The study of seed development in the model species Medicago truncatula has made a significant contribution to our understanding of this process in crop legumes. Thanks to the availability of comprehensive proteomics and transcriptomics databases, coupled with exhaustive mutant collections, the roles of several regulatory genes in development and maturation are beginning to be deciphered and functionally validated. Advances in next-generation sequencing and the availability of a genomic sequence have made feasible high-density SNP genotyping, allowing the identification of markers tightly linked to traits of agronomic interest. A further major advance is to be expected from the integration of omics resources in functional network construction, which has been used recently to identify "hub" genes central to important traits.

The Medicago truncatula GRAS protein RAD1 supports arbuscular mycorrhiza symbiosis and Phytophthora palmivora susceptibility

The roots of most land plants are colonized by symbiotic arbuscular mycorrhiza (AM) fungi. To facilitate this symbiosis, plant genomes encode a set of genes required for microbial perception and accommodation. However, the extent to which infection by filamentous root pathogens also relies on some of these genes remains an open question. Here, we used genome-wide association mapping to identify genes contributing to colonization of Medicago truncatula roots by the pathogenic oomycete Phytophthora palmivora. Single-nucleotide polymorphism (SNP) markers most significantly associated with plant colonization response were identified upstream of RAD1, which encodes a GRAS transcription regulator first negatively implicated in root nodule symbiosis and recently identified as a positive regulator of AM symbiosis. RAD1 transcript levels are up-regulated both in response to AM fungus and, to a lower extent, in infected tissues by P. palmivora where its expression is restricted to root cortex cells proximal to pathogen hyphae. Reverse genetics showed that reduction of RAD1 transcript levels as well as a rad1 mutant are impaired in their full colonization by AM fungi as well as by P. palmivora. Thus, the importance of RAD1 extends beyond symbiotic interactions, suggesting a general involvement in M. truncatula microbe-induced root development and interactions with unrelated beneficial and detrimental filamentous microbes.

Medicago truncatula: Genetic and Genomic Resources

Medicago truncatula was chosen by the legume community, along with Lotus japonicus, as a model plant to study legume biology. Since then, numerous resources and tools have been developed for M. truncatula. These include, for example, its genome sequence, core ecotype collections, transformation/regeneration methods, extensive mutant collections, and a gene expression atlas. This review aims to describe the different genetic and genomic tools and resources currently available for M. truncatula. We also describe how these resources were generated and provide all the information necessary to access these resources and use them from a practical point of view. © 2017 by John Wiley & Sons, Inc.

Is there genetic variation in mycorrhization of Medicago truncatula?

Differences in the plant's response among ecotypes or accessions are often used to identify molecular markers for the respective process. In order to analyze genetic diversity of Medicago truncatula in respect to interaction with the arbuscular mycorrhizal (AM) fungus Rhizophagus irregularis, mycorrhizal colonization was evaluated in 32 lines of the nested core collection representing the genetic diversity of the SARDI collection. All studied lines and the reference line Jemalong A17 were inoculated with R. irregularis and the mycorrhization rate was determined at three time points after inoculation. There were, however, no reliable and consistent differences in mycorrhization rates among all lines. To circumvent possible overlay of potential differences by use of the highly effective inoculum, native sandy soil was used in an independent experiment. Here, significant differences in mycorrhization rates among few of the lines were detectable, but the overall high variability in the mycorrhization rate hindered clear conclusions. To narrow down the number of lines to be tested in more detail, root system architecture (RSA) of in vitro-grown seedlings of all lines under two different phosphate (Pi) supply condition was determined in terms of primary root length and number of lateral roots. Under high Pi supply (100 µM), only minor differences were observed, whereas in response to Pi-limitation (3 µM) several lines exhibited a drastically changed number of lateral roots. Five lines showing the highest alterations or deviations in RSA were selected and inoculated with R. irregularis using two different Pi-fertilization regimes with either 13 mM or 3 mM Pi. Mycorrhization rate of these lines was checked in detail by molecular markers, such as transcript levels of RiTubulin and MtPT4. Under high phosphate supply, the ecotypes L000368 and L000555 exhibited slightly increased fungal colonization and more functional arbuscules, respectively. To address the question, whether capability for mycorrhizal colonization might be correlated to general invasion by microorganisms, selected lines were checked for infection by the root rot causing pathogen, Aphanoymces euteiches. The mycorrhizal colonization phenotype, however, did not correlate with the resistance phenotype upon infection with two strains of A. euteiches as L000368 showed partial resistance and L000555 exhibited high susceptibility as determined by quantification of A. euteiches rRNA within infected roots. Although there is genetic diversity in respect to pathogen infection, genetic diversity in mycorrhizal colonization of M. truncatula is rather low and it will be rather difficult to use it as a trait to access genetic markers.

Genome-Wide Association Study Reveals Natural Variations Contributing to Drought Resistance in Crops

Crops are often cultivated in regions where they will face environmental adversities; resulting in substantial yield loss which can ultimately lead to food and societal problems. Thus, significant efforts have been made to breed stress tolerant cultivars in an attempt to minimize these problems and to produce more stability with respect to crop yields across broad geographies. Since stress tolerance is a complex and multi-genic trait, advancements with classical breeding approaches have been challenging. On the other hand, molecular breeding, which is based on transgenics, marker-assisted selection and genome editing technologies; holds great promise to enable farmers to better cope with these challenges. However, identification of the key genetic components underlying the trait is critical and will serve as the foundation for future crop genetic improvement. Recently, genome-wide association studies have made significant contributions to facilitate the discovery of natural variation contributing to stress tolerance in crops. From these studies, the identified loci can serve as targets for genomic selection or editing to enable the molecular design of new cultivars. Here, we summarize research progress on this issue and focus on the genetic basis of drought tolerance as revealed by genome-wide association studies and quantitative trait loci mapping. Although many favorable loci have been identified, elucidation of their molecular mechanisms contributing to increased stress tolerance still remains a challenge. Thus, continuous efforts are still required to functionally dissect this complex trait through comprehensive approaches, such as system biological studies. It is expected that proper application of the acquired knowledge will enable the development of stress tolerant cultivars; allowing agricultural production to become more sustainable under dynamic environmental conditions.

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.

Identification of loci controlling forage yield and nutritive value in diploid alfalfa using GBS-GWAS

We attempted to identify genomic regions controlling forage yield and nutritive value in alfalfa. Several candidate genes and associated genetic markers were identified that could potentially be useful for alfalfa breeding to more efficiently develop improved cultivars. Alfalfa is one of the most widely cultivated forage legumes worldwide and improving alfalfa forage yield and nutritive value is a major global breeding goal. Genotyping-by-sequencing (GBS) provides cost-effective molecular marker genotyping for genome-wide association studies (GWAS). Using more than 15,000 genome-wide single nucleotide polymorphisms (SNP) identified from GBS, we conducted a GWAS to investigate forage yield and nutritive value-related traits. We have detected a number of associations for all the traits evaluated and a number of associations detected were located on the Medicago truncatula genome. The SNP in a coding region of a cell wall biosynthesis gene was associated with several cell wall-related traits, and we suggest that it may be the causative polymorphism. Two other SNPs residing in meristematic development and early growth genes were found to associate with the total biomass yield. None of the SNPs associated with regrowth after harvest or with spring regrowth were mapped to the M. truncatula genome, possibly reflecting the fact that M. truncatula is an annual species related to alfalfa that typically has limited ability to regrow. The alleles we identify with the major impact on forage yield and nutritive value can be rapidly incorporated into our breeding program.

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.

Disulfide isomerase-like protein AtPDIL1-2 is a good candidate for trichlorophenol phytodetoxification

Trichlorophenol (TCP) is a widely used and persistent environmentally toxic compound that poses a carcinogenic risk to humans. Phytoremediation is a proficient cleanup technology for organic pollutants. In this study, we found that the disulfide isomerase-like protein AtPDIL1-2 in plants is a good candidate for enhancing 2,4,6-TCP phytoremediation. The expression of AtPDIL1-2 in Arabidopsis was induced by 2,4,6-TCP. The heterologously expressed AtPDIL1-2 in Escherichia coli exhibited both oxidase and isomerase activities as protein disulfide isomerase and improved bacteria tolerance to 2,4,6-TCP. Further research revealed that transgenic tobacco overexpressing AtPDIL1-2 was more tolerant to high concentrations of 2,4,6-TCP and removed the toxic compound at far greater rates than the control plants. To elucidate the mechanism of action of AtPDIL1-2, we investigated the chemical interaction of AtPDIL1-2 with 2,4,6-TCP for the first time. HPLC analysis implied that AtPDIL1-2 exerts a TCP-binding activity. A suitable configuration of AtPDIL1-2-TCP binding was obtained by molecular docking studies using the AutoDock program. It predicted that the TCP binding site is located in the b-b' domain of AtPDIL1-2 and that His254 of the protein is critical for the binding interaction. These findings imply that AtPDIL1-2 can be used for TCP detoxification by the way of overexpression in plants.

Genome wide association study (GWAS) for grain yield in rice cultivated under water deficit

The identification of rice drought tolerant materials is crucial for the development of best performing cultivars for the upland cultivation system. This study aimed to identify markers and candidate genes associated with drought tolerance by Genome Wide Association Study analysis, in order to develop tools for use in rice breeding programs. This analysis was made with 175 upland rice accessions (Oryza sativa), evaluated in experiments with and without water restriction, and 150,325 SNPs. Thirteen SNP markers associated with yield under drought conditions were identified. Through stepwise regression analysis, eight SNP markers were selected and validated in silico, and when tested by PCR, two out of the eight SNP markers were able to identify a group of rice genotypes with higher productivity under drought. These results are encouraging for deriving markers for the routine analysis of marker assisted selection. From the drought experiment, including the genes inherited in linkage blocks, 50 genes were identified, from which 30 were annotated, and 10 were previously related to drought and/or abiotic stress tolerance, such as the transcription factors WRKY and Apetala2, and protein kinases.

Evaluating Agronomic Performance and Investigating Molecular Structure of Drought and Heat Tolerant Wild Alfalfa (Medicago sativa L.) Collection from the Southeastern Turkey

Drought is a major stress factor for agricultural production including alfalfa production. One way to counterbalance the yield losses is the introgression of drought tolerant germplasm into breeding programs. As an effort to exploit such germplasm, 16 individual plants were selected from the Southeastern Turkey from their natural habitat and clonally propagated in field trials with an ultimate goal to use the germplasm as parents for releasing a synthetic cultivar. Forage yield and forage quality traits were evaluated and molecular genetic diversity among genotypes were determined using inter simple sequence repeat markers. Genotypes showed a variation from growth habit to yield and quality traits indicating sufficient phenotypic variation for diverse breeding efforts (for grazing or harvesting) and long term selection schemes. A large amount of genetic variation was observed even with a limited number of marker and genotypes. However, no pattern of spatial genetic structure was observed for the scale of the study when genetic variation is linked to the geographic origin. We conclude that ex situ natural variation provides a wealth of germplasm that could be incorporated into breeding programs aiming to improve drought tolerance. We also suggest an extensive collection of seeds/plant tissue from unique plants with desirable traits rather than putting more efforts to create a spatial germplasm sampling efforts in narrow regions.

Genetic Mapping of a Major Resistance Gene to Pea Aphid (Acyrthosipon pisum) in the Model Legume Medicago truncatula

Resistance to the Australian pea aphid (PA; Acyrthosiphon pisum) biotype in cultivar Jester of the model legume Medicago truncatula is mediated by a single dominant gene and is phloem-mediated. The genetic map position for this resistance gene, APR (Acyrthosiphon pisum resistance), is provided and shows that APR maps 39 centiMorgans (cM) distal of the A. kondoi resistance (AKR) locus, which mediates resistance to a closely related species of the same genus bluegreen aphid (A. kondoi). The APR region on chromosome 3 is dense in classical nucleotide binding site leucine-rich repeats (NLRs) and overlaps with the region harbouring the RAP1 gene which confers resistance to a European PA biotype in the accession Jemalong A17. Further screening of a core collection of M. truncatula accessions identified seven lines with strong resistance to PA. Allelism experiments showed that the single dominant resistance to PA in M. truncatula accessions SA10481 and SA1516 are allelic to SA10733, the donor of the APR locus in cultivar Jester. While it remains unclear whether there are multiple PA resistance genes in an R-gene cluster or the resistance loci identified in the other M. truncatula accessions are allelic to APR, the introgression of APR into current M. truncatula cultivars will provide more durable resistance to PA.

A Snapshot of Functional Genetic Studies in Medicago truncatula

In the current context of food security, increase of plant protein production in a sustainable manner represents one of the major challenges of agronomic research, which could be partially resolved by increased cultivation of legume crops. Medicago truncatula is now a well-established model for legume genomic and genetic studies. With the establishment of genomics tools and mutant populations in M. truncatula, it has become an important resource to answer some of the basic biological questions related to plant development and stress tolerance. This review has an objective to overview a decade of genetic studies in this model plant from generation of mutant populations to nowadays. To date, the three biological fields, which have been extensively studied in M. truncatula, are the symbiotic nitrogen fixation, the seed development, and the abiotic stress tolerance, due to their significant agronomic impacts. In this review, we summarize functional genetic studies related to these three major biological fields. We integrated analyses of a nearly exhaustive list of genes into their biological contexts in order to provide an overview of the forefront research advances in this important legume model plant.