Genetic determinants of seed protein plasticity in response to the environment in Medicago truncatula.

As the frequency of extreme environmental events is expected to increase with climate change, identifying candidate genes for stabilizing the protein composition of legume seeds or optimizing this in a given environment is increasingly important. To elucidate the genetic determinants of seed protein plasticity, major seed proteins from 200 ecotypes of Medicago truncatula grown in four contrasting environments were quantified after one-dimensional electrophoresis. The plasticity index of these proteins was recorded for each genotype as the slope of Finlay and Wilkinson's regression and then used for genome-wide association studies (GWASs), enabling the identification of candidate genes for determining this plasticity. This list was enriched in genes related to transcription, DNA repair and signal transduction, with many of them being stress responsive. Other over-represented genes were related to sulfur and aspartate family pathways leading to the synthesis of the nutritionally essential amino acids methionine and lysine. By placing these genes in metabolic pathways, and using a M. truncatula mutant impaired in regenerating methionine from S-methylmethionine, we discovered that methionine recycling pathways are major contributors to globulin composition establishment and plasticity. These data provide a unique resource of genes that can be targeted to mitigate negative impacts of environmental stresses on seed protein composition.

Exploring natural diversity of Medicago truncatula reveals physiotypes and loci associated with the response of seedling performance to nitrate supply.

Seedling pre-emergence is a critical phase of development for successful crop establishment because of its susceptibility to environmental conditions. In a context of reduced use of inorganic fertilizers, the genetic bases of the response of seedlings to nitrate supply received little attention. This issue is important even in legumes where nitrate absorption starts early after germination, before nodule development. Natural variation of traits characterizing seedling growth in the absence or presence of nitrate was investigated in a core collection of 192 accessions of Medicago truncatula. Plasticity indexes to the absence of nitrate were calculated. The genetic determinism of the traits was dissected by genome-wide association study (GWAS). The absence of nitrate affected seed biomass mobilization and root/shoot length ratio. However, the large range of genetic variability revealed different seedling performances within natural diversity. A principal component analysis (PCA) carried out with plasticity indexes highlighted four physiotypes of accessions differing in relationships between seedling elongation and seed biomass partitioning traits in response to the absence of nitrate. Finally, GWAS revealed 45 associations with single or combined traits corresponding to coordinates of accessions on PCA, as well as two clusters of genes encoding sugar transporters and glutathione transferases surrounding loci associated with seedling elongation traits.

Genome-wide association studies with proteomics data reveal genes important for synthesis, transport and packaging of globulins in legume seeds.

Improving nutritional seed quality is an important challenge in grain legume breeding. However, the genes controlling the differential accumulation of globulins, which are major contributors to seed nutritional value in legumes, remain largely unknown. We combined a search for protein quantity loci with genome-wide association studies on the abundance of 7S and 11S globulins in seeds of the model legume species Medicago truncatula. Identified genomic regions and genes carrying polymorphisms linked to globulin variations were then cross-compared with pea (Pisum sativum), leading to the identification of candidate genes for the regulation of globulin abundance in this crop. Key candidates identified include genes involved in transcription, chromatin remodeling, post-translational modifications, transport and targeting of proteins to storage vacuoles. Inference of a gene coexpression network of 12 candidate transcription factors and globulin genes revealed the transcription factor ABA-insensitive 5 (ABI5) as a highly connected hub. Characterization of loss-of-function abi5 mutants in pea uncovered a role for ABI5 in controlling the relative abundance of vicilin, a sulfur-poor 7S globulin, in pea seeds. This demonstrates the feasibility of using genome-wide association studies in M. truncatula to reveal genes that can be modulated to improve seed nutritional value.

Adaptation to climate through flowering phenology: a case study in Medicago truncatula.

Local climatic conditions likely constitute an important selective pressure on genes underlying important fitness-related traits such as flowering time, and in many species, flowering phenology and climatic gradients strongly covary. To test whether climate shapes the genetic variation on flowering time genes and to identify candidate flowering genes involved in the adaptation to environmental heterogeneity, we used a large Medicago truncatula core collection to examine the association between nucleotide polymorphisms at 224 candidate genes and both climate variables and flowering phenotypes. Unlike genome-wide studies, candidate gene approaches are expected to enrich for the number of meaningful trait associations because they specifically target genes that are known to affect the trait of interest. We found that flowering time mediates adaptation to climatic conditions mainly by variation at genes located upstream in the flowering pathways, close to the environmental stimuli. Variables related to the annual precipitation regime reflected selective constraints on flowering time genes better than the other variables tested (temperature, altitude, latitude or longitude). By comparing phenotype and climate associations, we identified 12 flowering genes as the most promising candidates responsible for phenological adaptation to climate. Four of these genes were located in the known flowering time QTL region on chromosome 7. However, climate and flowering associations also highlighted largely distinct gene sets, suggesting different genetic architectures for adaptation to climate and flowering onset.

Cell length instead of cell number becomes the predominant factor contributing to hypocotyl length genotypic differences under abiotic stress in Medicago truncatula.

Hypocotyl elongation in the dark is a crucial process to ensure seedling emergence. It relies both on the cell number and cell length. The contribution of these two factors to the maximal hypocotyl length and the impact of environmental conditions on this contribution are not known. This is surprising considering the agronomic and economical importance of seedling emergence in crop establishment. Using 14 genotypes from a nested core collection representing Medicago truncatula (barrel medic) natural variation, we investigated how epidermal cell number and cell length contribute to hypocotyl length under optimal, low temperature (8°C) and water deficit (-0.50 MPa) conditions. Both cell number and length vary according to genotypes and contribute to maximal hypocotyl length differences between genotypes. This contribution, however, depends on growth conditions. Cell number is the major contributor under optimal conditions (60%) whereas cell length becomes the major determinant under stress. Maximal hypocotyl length is correlated with hypocotyl elongation rate under both stresses but not under optimal condition, revealing contrasted genotypes for cell elongation capacity under stress. To identify the genetic regulators determining cell number and cell length, quantitative trait loci (QTLs) were detected using a recombinant inbred lines population exhibiting segregation in maximal hypocotyl length. Two QTLs controlling cell number and three QTLs controlling cell length at low temperature were detected. One QTL for cell number and two for cell length were found to be associated with hypocotyl length under low temperature. This study provides new information to improve seedling emergence under abiotic stress.

High-density genome-wide association mapping implicates an F-box encoding gene in Medicago truncatula resistance to Aphanomyces euteiches.

• The use of quantitative disease resistance (QDR) is a promising strategy for promoting durable resistance to plant pathogens, but genes involved in QDR are largely unknown. To identify genetic components and accelerate improvement of QDR in legumes to the root pathogen Aphanomyces euteiches, we took advantage of both the recently generated massive genomic data for Medicago truncatula and natural variation of this model legume. • A high-density (≈5.1 million single nucleotide polymorphisms (SNPs)) genome-wide association study (GWAS) was performed with both in vitro and glasshouse phenotyping data collected for 179 lines. • GWAS identified several candidate genes and pinpointed two independent major loci on the top of chromosome 3 that were detected in both phenotyping methods. Candidate SNPs in the most significant locus (σ(A)²= 23%) were in the promoter and coding regions of an F-box protein coding gene. Subsequent qRT-PCR and bioinformatic analyses performed on 20 lines demonstrated that resistance is associated with mutations directly affecting the interaction domain of the F-box protein rather than gene expression. • These results refine the position of previously identified QTL to specific candidate genes, suggest potential molecular mechanisms, and identify new loci explaining QDR against A. euteiches.

Natural diversity in the model legume Medicago truncatula allows identifying distinct genetic mechanisms conferring partial resistance to Verticillium wilt.

Verticillium wilt is a major threat to alfalfa (Medicago sativa) and many other crops. The model legume Medicago truncatula was used as a host for studying resistance and susceptibility to Verticillium albo-atrum. In addition to presenting well-established genetic resources, this wild plant species enables to investigate biodiversity of the response to the pathogen and putative crosstalk between disease and symbiosis. Symptom scoring after root inoculation and modelling of disease curves allowed assessing susceptibility levels in recombinant lines of three crosses between susceptible and resistant lines, in a core collection of 32 lines, and in mutants affected in symbiosis with rhizobia. A GFP-expressing V. albo-atrum strain was used to study colonization of susceptible plants. Symptoms and colonization pattern in infected M. truncatula plants were typical of Verticillium wilt. Three distinct major quantitative trait loci were identified using a multicross, multisite design, suggesting that simple genetic mechanisms appear to control Verticillium wilt resistance in M. truncatula lines A17 and DZA45.5. The disease functional parameters varied largely in lines of the core collection. This biodiversity with regard to disease response encourages the development of association genetics and ecological approaches. Several mutants of the resistant line, impaired in different steps of rhizobial symbiosis, were affected in their response to V. albo-atrum, which suggests that mechanisms involved in the establishment of symbiosis or disease might have some common regulatory control points.

Genetic variability and QTL mapping of freezing tolerance and related traits in Medicago truncatula.

Freezing is a major environmental limitation to crop productivity for a number of species including legumes. We investigated the genetic determinism of freezing tolerance in the model legume Medicago truncatula Gaertn (M. truncatula). After having observed a large variation for freezing tolerance among 15 M. truncatula accessions, the progeny of a F6 recombinant inbred line population, derived from a cross between two accessions, was acclimated to low above-freezing temperatures and assessed for: (a) number of leaves (NOL), leaf area (LA), chlorophyll content index (CCI), shoot and root dry weights (SDW and RDW) at the end of the acclimation period and (b) visual freezing damage (FD) during the freezing treatment and 2 weeks after regrowth and foliar electrolyte leakage (EL) 2 weeks after regrowth. Consistent QTL positions with additive effects for FD were found on LG1, LG4 and LG6, the latter being the most explanatory (R (2) ≈ 40 %). QTL for NOL, QTL for EL, NOL and RDW, and QTL for EL and CCI colocalized with FD QTL on LG1, LG4 and LG6, respectively. Favorable alleles for these additive effects were brought by the same parent suggesting that this accession contributes to superior freezing tolerance by affecting plants' capacity to maintain growth at low above-freezing temperatures. No epistatic effects were found between FD QTL, but for each of the studied traits, 3-6 epistatic effects were detected between loci not detected directly as QTL. These results open the way to the assessment of syntenic relationships between QTL for frost tolerance in M. truncatula and cultivated legume species.

Using a physiological framework for improving the detection of quantitative trait loci related to nitrogen nutrition in Medicago truncatula.

Medicago truncatula is used as a model plant for exploring the genetic and molecular determinants of nitrogen (N) nutrition in legumes. In this study, our aim was to detect quantitative trait loci (QTL) controlling plant N nutrition using a simple framework of carbon/N plant functioning stemming from crop physiology. This framework was based on efficiency variables which delineated the plant's efficiency to take up and process carbon and N resources. A recombinant inbred line population (LR4) was grown in a glasshouse experiment under two contrasting nitrate concentrations. At low nitrate, symbiotic N(2) fixation was the main N source for plant growth and a QTL with a large effect located on linkage group (LG) 8 affected all the traits. Significantly, efficiency variables were necessary both to precisely localize a second QTL on LG5 and to detect a third QTL involved in epistatic interactions on LG2. At high nitrate, nitrate assimilation was the main N source and a larger number of QTL with weaker effects were identified compared to low nitrate. Only two QTL were common to both nitrate treatments: a QTL of belowground biomass located at the bottom of LG3 and another one on LG6 related to three different variables (leaf area, specific N uptake and aboveground:belowground biomass ratio). Possible functions of several candidate genes underlying QTL of efficiency variables could be proposed. Altogether, our results provided new insights into the genetic control of N nutrition in M. truncatula. For instance, a novel result for M. truncatula was identification of two epistatic interactions in controlling plant N(2) fixation. As such this study showed the value of a simple conceptual framework based on efficiency variables for studying genetic determinants of complex traits and particularly epistatic interactions.

A genome-wide analysis suggests pleiotropic effects of Green Revolution genes on shade avoidance in wheat.

A classic example of phenotypic plasticity in plants is the suit of phenotypic responses induced by a change in the ratio of red to far-red light (R∶FR) as a result of shading, also known as the shade avoidance syndrome (SAS). While the adaptive consequences of this syndrome have been extensively discussed in natural ecosystems, how SAS varies within crop populations and how SAS evolved during crop domestication and breeding remain poorly known. In this study, we grew a panel of 180 durum wheat (Triticum turgidum ssp. durum) genotypes spanning diversity from wild, early domesticated, and elite genetic compartments under two light treatments: low R:FR light (shaded treatment) and high R:FR light (unshaded treatment). We first quantified the genetic variability of SAS, here measured as a change in plant height at the seedling stage. We then dissected the genetic basis of this variation through genome-wide association mapping. Genotypes grown in shaded conditions were taller than those grown under unshaded conditions. Interaction between light quality and genotype did not affect plant height. We found six QTLs affecting plant height. Three significantly interacted with light quality among which the well-known Rht1 gene introgressed in elite germplasm during the Green Revolution. Interestingly at three loci, short genotypes systematically expressed reduced SAS, suggesting a positive genetic correlation between plant height and plant height plasticity. Overall, our study sheds light on the evolutionary history of crops and illustrates the relevance of genetic approaches to tackle agricultural challenges.

Quantitative trait loci analysis reveals a correlation between the ratio of sucrose/raffinose family oligosaccharides and seed vigour in Medicago truncatula.

Seed vigour is important for successful establishment and high yield, especially under suboptimal environmental conditions. In legumes, raffinose oligosaccharide family (RFO) sugars have been proposed as an easily available energy reserve for seedling establishment. In this study, we investigated whether the composition or amount of soluble sugars (sucrose and RFO) is part of the genetic determinants of seed vigour of Medicago truncatula using two recombinant inbred line (RIL) populations. Quantitative trait loci (QTL) mapping for germination rate, hypocotyl and radicle growth under water deficit and nutritional stress, seed weight and soluble sugar content was performed using RIL populations LR1 and LR4. Seven of the 12 chromosomal regions containing QTL for germination rate or post-germinative radicle growth under optimal or stress conditions co-located with Suc/RFO QTL. A significant negative correlation was also found between seed vigour traits and Suc/RFO. In addition, one QTL that explained 80% of the variation in the ratio stachyose/verbascose co-located with a stachyose synthase gene whose expression profile in the parental lines could explain the variation in oligosaccharide composition. The correlation and co-location of Suc/RFO ratio with germination and radicle growth QTL suggest that an increased Suc/RFO ratio in seeds of M. truncatula might negatively affect seed vigour.

A complex genetic network involving a broad-spectrum locus and strain-specific loci controls resistance to different pathotypes of Aphanomyces euteiches in Medicago truncatula.

A higher understanding of genetic and genomic bases of partial resistance in plants and their diversity regarding pathogen variability is required for a more durable management of resistance genetic factors in sustainable cropping systems. In this study, we investigated the diversity of genetic factors involved in partial resistance to Aphanomyces euteiches, a very damaging pathogen on pea and alfalfa, in Medicago truncatula. A mapping population of 178 recombinant inbred lines, from the cross F83005.5 (susceptible) and DZA045.5 (resistant), was used to identify quantitative trait loci for resistance to four A. euteiches reference strains belonging to the four main pathotypes currently known on pea and alfalfa. A major broad-spectrum genomic region, previously named AER1, was localized to a reduced 440 kb interval on chromosome 3 and was involved in complete or partial resistance, depending on the A. euteiches strain. We also identified 21 additive and/or epistatic genomic regions specific to one or two strains, several of them being anchored to the M. truncatula physical map. These results show that, in M. truncatula, a complex network of genetic loci controls partial resistance to different pea and alfalfa pathotypes of A. euteiches, suggesting a diversity of molecular mechanisms underlying partial resistance.

The RAP1 gene confers effective, race-specific resistance to the pea aphid in Medicago truncatula independent of the hypersensitive reaction.

Plant resistance to pathogens is commonly associated with a hypersensitive response (HR), but the degree to which the HR is responsible for incompatibility is subject to debate. Resistance to aphids is likely to share features with resistance to pathogens but is less well understood. Here, we report effective resistance to the pea aphid Acyrthosiphon pisum in Medicago truncatula. Aphids lost weight and died rapidly (within two days) on the resistant genotype Jemalong, which developed necrotic lesions following infestation. Lesions were induced by nonvascular intracellular stylet punctures by aphids, remained localized to the site of stylet entry, stained for the presence of reactive oxygen species, and were similar to the HR induced by the bacterial pathogen Pseudomonas syringae pv. phaseolicola. The implication that aphid-induced lesions confer resistance was tested by quantitative trait loci analysis using recombinant inbred lines derived from a cross between Jemalong and the susceptible genotype DZA315.16. One major locus, RAP1, was identified that was sufficient to confer race-specific resistance against the pea aphid and was mapped to the middle of chromosome 3. Surprisingly, a separate locus, mapping to the top of chromosome 3, governed aphid-induced HR, indicating that the HR-like lesions are not required for RAP1-mediated aphid resistance.

Characterization of Root-Knot Nematode Resistance in Medicago truncatula.

Root knot (Meloidogyne spp.) and cyst (Heterodera and Globodera spp.) nematodes infect all important crop species, and the annual economic loss due to these pathogens exceeds $90 billion. We screened the worldwide accession collection with the root-knot nematodes Meloidogyne incognita, M. arenaria and M. hapla, soybean cyst nematode (SCN-Heterodera glycines), sugar beet cyst nematode (SBCN-Heterodera schachtii) and clover cyst nematode (CLCN-Heterodera trifolii), revealing resistant and susceptible accessions. In the over 100 accessions evaluated, we observed a range of responses to the root-knot nematode species, and a non-host response was observed for SCN and SBCN infection. However, variation was observed with respect to infection by CLCN. While many cultivars including Jemalong A17 were resistant to H. trifolii, cultivar Paraggio was highly susceptible. Identification of M. truncatula as a host for root-knot nematodes and H. trifolii and the differential host response to both RKN and CLCN provide the opportunity to genetically and molecularly characterize genes involved in plant-nematode interaction. Accession DZA045, obtained from an Algerian population, was resistant to all three root-knot nematode species and was used for further studies. The mechanism of resistance in DZA045 appears different from Mi-mediated root-knot nematode resistance in tomato. Temporal analysis of nematode infection showed that there is no difference in nematode penetration between the resistant and susceptible accessions, and no hypersensitive response was observed in the resistant accession even several days after infection. However, less than 5% of the nematode population completed the life cycle as females in the resistant accession. The remainder emigrated from the roots, developed as males, or died inside the roots as undeveloped larvae. Genetic analyses carried out by crossing DZA045 with a susceptible French accession, F83005, suggest that one gene controls resistance in DZA045.

Microsatellite diversity and broad scale geographic structure in a model legume: building a set of nested core collection for studying naturally occurring variation in Medicago truncatula.

BACKGROUND: Exploiting genetic diversity requires previous knowledge of the extent and structure of the variation occurring in a species. Such knowledge can in turn be used to build a core-collection, i.e. a subset of accessions that aim at representing the genetic diversity of this species with a minimum of repetitiveness. We investigate the patterns of genetic diversity and population structure in a collection of 346 inbred lines representing the breadth of naturally occurring diversity in the Legume plant model Medicago truncatula using 13 microsatellite loci distributed throughout the genome. RESULTS: We confirm the uniqueness of all these genotypes and reveal a large amount of genetic diversity and allelic variation within this autogamous species. Spatial genetic correlation was found only for individuals originating from the same population and between neighbouring populations. Using a model-based clustering algorithm, we identified four main genetic clusters in the set of individuals analyzed. This stratification matches broad geographic regions. We also identified a set of "admixed" individuals that do not fit with this population structure scheme. CONCLUSION: The stratification inferred is discussed considering potential historical events like expansion, refuge history and admixture between neighbouring groups. Information on the allelic richness and the inferred population structure are used to build a nested core-collection. The set of inbred lines and the core collections are publicly available and will help coordinating efforts for the study of naturally occurring variation in the growing Medicago truncatula community.

Potential Implications of Climate Change on Aegilops Species Distribution: Sympatry of These Crop Wild Relatives with the Major European Crop Triticum aestivum and Conservation Issues.

Gene flow from crop to wild relatives is a common phenomenon which can lead to reduced adaptation of the wild relatives to natural ecosystems and/or increased adaptation to agrosystems (weediness). With global warming, wild relative distributions will likely change, thus modifying the width and/or location of co-occurrence zones where crop-wild hybridization events could occur (sympatry). This study investigates current and 2050 projected changes in sympatry levels between cultivated wheat and six of the most common Aegilops species in Europe. Projections were generated using MaxEnt on presence-only data, bioclimatic variables, and considering two migration hypotheses and two 2050 climate scenarios (RCP4.5 and RCP8.5). Overall, a general decline in suitable climatic conditions for Aegilops species outside the European zone and a parallel increase in Europe were predicted. If no migration could occur, the decline was predicted to be more acute outside than within the European zone. The potential sympatry level in Europe by 2050 was predicted to increase at a higher rate than species richness, and most expansions were predicted to occur in three countries, which are currently among the top four wheat producers in Europe: Russia, France and Ukraine. The results are also discussed with regard to conservation issues of these crop wild relatives.

Development and composition of the seeds of nine genotypes of the Medicago truncatula species complex.

The seed development and composition of Medicago truncatula Gaertn., the new model plant for grain legumes, was studied using nine genotypes of the species complex: M. truncatula-Medicago littoralis (M. truncatula). The seed development of M. truncatula was very similar to that of other legumes, the only notable exception being the presence, in the mature seed, of an endosperm layer that is absent in grain legumes. During early embryogenesis and until mid-maturation, transient storage of starch occurred in the seed coat and embryo. This temporary storage probably contributed to the early development of the embryo and reserve synthesis. During maturation the synthesis and accumulation of proteins and oil took place at quasi-constant rates. Conversely oligosaccharides, mainly stachyose, were synthesised only during late maturation and at the beginning of desiccation. Proteins represented the major class of storage compounds and their average amino acid composition was found to be very close to that of pea and robust in various environmental conditions. Similar compositions between the two species and other grain legumes were also found for the fatty acids and the soluble sugars; most of these characters varied depending on the various environmental conditions used for seed production. All these similarities fully justify the use of M. truncatula as a model plant for genomic approaches to grain legume improvement. The major difference between M. truncatula seeds and European grain legume seeds resides in the nature of their carbon storage namely triacylglycerides for M. truncatula and starch for pea and faba bean.

Cytological, genetic, and molecular analysis to characterize compatible and incompatible interactions between Medicago truncatula and Colletotrichum trifolii.

In this study, a new pathosystem was established using the model plant Medicago truncatula and Colletotrichum trifolii, the causal agent of anthracnose on Medicago sativa. Screening of a few M. truncatula lines identified Jemalong and F83005.5 as resistant and susceptible to Colletotrichum trifolii race 1, respectively. Symptom analysis and cytological studies indicated that resistance of Jemalong was associated with a hypersensitive response of the plant. The two selected lines were crossed, and inoculations with C. trifolii were performed on the resulting F1 and F2 progenies. Examination of the disease phenotypes indicated that resistance was dominant and was probably due to a major resistance gene. Molecular components of the resistance were analyzed through macroarray experiments. Expression profiling of 126 expressed sequence tags corresponding to 92 genes, which were selected for their putative functions in plant defense or signal transduction, were compared in Jemalong and F83005.5 lines. A strong correlation was observed between the number of up-regulated genes and the resistance phenotype. Large differences appeared at 48 h postinoculation; more than 40% of the tested genes were up-regulated in the Jemalong line compared with only 10% in the susceptible line. Interestingly, some nodulin genes were also induced in the resistant line upon inoculation with C. trifolii.

The molecular genetic linkage map of the model legume Medicago truncatula: an essential tool for comparative legume genomics and the isolation of agronomically important genes.

BACKGROUND: The legume Medicago truncatula has emerged as a model plant for the molecular and genetic dissection of various plant processes involved in rhizobial, mycorrhizal and pathogenic plant-microbe interactions. Aiming to develop essential tools for such genetic approaches, we have established the first genetic map of this species. Two parental homozygous lines were selected from the cultivar Jemalong and from the Algerian natural population (DZA315) on the basis of their molecular and phenotypic polymorphism. RESULTS: An F2 segregating population of 124 individuals between these two lines was obtained using an efficient manual crossing technique established for M. truncatula and was used to construct a genetic map. This map spans 1225 cM (average 470 kb/cM) and comprises 289 markers including RAPD, AFLP, known genes and isoenzymes arranged in 8 linkage groups (2n = 16). Markers are uniformly distributed throughout the map and segregation distortion is limited to only 3 linkage groups. By mapping a number of common markers, the eight linkage groups are shown to be homologous to those of diploid alfalfa (M. sativa), implying a good level of macrosynteny between the two genomes. Using this M. truncatula map and the derived F3 populations, we were able to map the Mtsym6 symbiotic gene on linkage group 8 and the SPC gene, responsible for the direction of pod coiling, on linkage group 7. CONCLUSIONS: These results demonstrate that Medicago truncatula is amenable to diploid genetic analysis and they open the way to map-based cloning of symbiotic or other agronomically-important genes using this model plant.

Single-nucleotide polymorphism discovery and diversity in the model legume Medicago truncatula.

Extensive genomic resources are available in the model legume Medicago truncatula. Here, we present the discovery and design of the first array of single-nucleotide polymorphism (SNP) markers in M. truncatula through large-scale Sanger resequencing of genomic fragments spanning the genome, in a diverse panel of 16 M. truncatula accessions. Both anonymous fragments and fragments targeting candidate genes for flowering phenology and symbiosis were surveyed for nucleotide variation in almost 230 kb of unique genomic regions. A set of 384 SNP markers was designed for an Illumina's GoldenGate assay, genotyped on a collection of 192 inbred lines (CC192) representing the geographical range of the species and used to survey the diversity of two natural populations. Finally, 86% of the tested SNPs were of high quality and exhibited polymorphism in the CC192 collection. Even at the population level, we detected polymorphism for more than 50% of the selected SNPs. Analysis of the allele frequency spectrum in the CC192 showed a reduced ascertainment bias, mostly limited to very rare alleles (frequency <0.01). The substantial polymorphism detected at the species and population levels, the high marker quality and the potential to survey large samples of individuals make this set of SNP markers a valuable tool to improve our understanding of the effect of demographic and selective factors that shape the natural genetic diversity within the selfing species Medicago truncatula.