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

Key message

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.

Abstract

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.

     

Genetic and physical localization of an anthracnose resistance gene in <Emphasis Type="Italic">Medicago truncatula</Emphasis>

Anthracnose of alfalfa, caused by the fungal pathogen Colletotrichum trifolii, is one of the most destructive diseases of alfalfa worldwide. An improved understanding of the genetic and molecular mechanisms underlying host resistance will facilitate the development of resistant alfalfa cultivars, thus providing the most efficient and environmentally sound strategy to control alfalfa diseases. Unfortunately, cultivated alfalfa has an intractable genetic system because of its tetrasomic inheritance and out-crossing nature. Nevertheless, the model legume Medicago truncatula, a close relative of alfalfa, has the potential to serve as a surrogate to map and clone the counterparts of agronomically important genes in alfalfa—particularly, disease resistance genes against economically important pathogens. Here we describe the high-resolution genetic and physical mapping of RCT1, a host resistance gene against C. trifolii race 1 in M. truncatula. We have delimited the RCT1 locus within a physical interval spanning ∼200 kb located on the top of M. truncatula linkage group 4. RCT1 is part of a complex locus containing numerous genes homologous to previously characterized TIR-NBS-LRR type resistance genes. The result presented in this paper will facilitate the positional cloning of RCT1 in Medicago.     

A high-density genetic map of the Medicago truncatula major freezing tolerance QTL on chromosome 6 reveals colinearity with a QTL related to freezing damage on Pisum sativum linkage group VI

Freezing is one of the most serious abiotic stress factors that affect cool-season legumes. It limits species geographic distribution and causes severe yield losses. Improving tolerance to freezing has long been a main concern for legume breeders. Medicago truncatula Gaertn. has been selected as a model species for legume biology. Various studies have shown significant macrosynteny between M. truncatula and agronomically important crop legumes. A major freezing tolerance quantitative trait locus (QTL), herein referred to as Mt-FTQTL6, was previously identified on M. truncatula chromosome 6. The physical location of this QTL was determined in this study and its corresponding chromosomal interval was enriched with additional markers. Markers were first developed using the draft sequence of M. truncatula euchromatin (release versions Mt3.0 and Mt3.5). Because Mt-FTQTL6 was found to coincide with an assembly gap, the Glycine max (L.) Merr. genome sequence was also used to generate markers. Five Mt-FTQTL6-linked markers were found to be common to a region on Pisum sativum L. linkage group VI harboring a QTL for freezing damage. A subset of markers was tested for transferability across 11 additional legume species. This study lays the groundwork for identifying the molecular basis of Mt-FTQTL6. Cross-legume markers will be useful in future efforts aiming to investigate the conservation of Mt-FTQTL6 in cool-season legumes and subsequently the existence of common mechanisms for response to freezing between M. truncatula and crop legumes.     

Expressed Sequence Tags from a Root-Hair-Enriched Medicago truncatula cDNA Library

The root hair is a specialized cell type involved in water and nutrient uptake in plants. In legumes the root hair is also the primary site of recognition and infection by symbiotic nitrogen-fixing Rhizobium bacteria. We have studied the root hairs of Medicago truncatula, which is emerging as an increasingly important model legume for studies of symbiotic nodulation. However, only 27 genes from M. truncatula were represented in GenBank/EMBL as of October, 1997. We report here the construction of a root-hair-enriched cDNA library and single-pass sequencing of randomly selected clones. Expressed sequence tags (899 total, 603 of which have homology to known genes) were generated and made available on the Internet. We believe that the database and the associated DNA materials will provide a useful resource to the community of scientists studying the biology of roots, root tips, root hairs, and nodulation.     

Patterns of divergence of a large family of nodule cysteine‐rich peptides in accessions of Medicago truncatula

The nodule cysteine‐rich (NCR) groups of defensin‐like (DEFL) genes are one of the largest gene families expressed in the nodules of some legume plants. They have only been observed in the inverted repeat loss clade (IRLC) of legumes, which includes the model legume Medicago truncatula. NCRs are reported to play an important role in plant–microbe interactions. To understand their diversity we analyzed their expression and sequence polymorphisms among four accessions of M. truncatula. A significant expression and nucleotide variation was observed among the genes. We then used 26 accessions to estimate the selection pressures shaping evolution among the accessions by calculating the nucleotide diversity at non‐synonymous and synonymous sites in the coding region. The mature peptides of the orthologous NCRs had signatures of both purifying and diversifying selection pressures, unlike the seed DEFLs, which predominantly exhibited purifying selection. The expression, sequence variation and apparent diversifying selection in NCRs within the Medicago species indicates rapid and recent evolution, and suggests that this family of genes is actively evolving to adapt to different environments and is acquiring new functions.     

Nucleotide polymorphism of the Srlk gene that determines salt stress tolerance in alfalfa (Medicago sativa L.)

Based on legume genome syntheny, the nucleotide sequence of Srlk gene, key role of which in response to salt stress was demonstrated for the model species Medicago truncatula, was identified in the major forage and siderate crop alfalfa (Medicago sativa). In twelve alfalfa samples originating from regions with contrasting growing conditions, 19 SNPs were revealed in the Srlk gene. For two nonsynonymous SNPs, molecular markers were designed that could be further used to analyze the association between Srlk gene nucleotide polymorphism and the variability in salt stress tolerance among alfalfa cultivars.     

New evidence of ancestral polyploidy in the Genistoid legume Lupinus angustifolius L. (narrow-leafed lupin)

Key message

This is the first clear evidence of duplication and/or triplication of large chromosomal regions in a genome of a Genistoid legume, the most basal clade of Papilionoid legumes.

Abstract

Lupinus angustifolius L. (narrow-leafed lupin) is the most widely cultivated species of Genistoid legume, grown for its high-protein grain. As a member of this most basal clade of Papilionoid legumes, L. angustifolius serves as a useful model for exploring legume genome evolution. Here, we report an improved reference genetic map of L. angustifolius comprising 1207 loci, including 299 newly developed Diversity Arrays Technology markers and 54 new gene-based PCR markers. A comparison between the L. angustifolius and Medicago truncatula genomes was performed using 394 sequence-tagged site markers acting as bridging points between the two genomes. The improved L. angustifolius genetic map, the updated M. truncatula genome assembly and the increased number of bridging points between the genomes together substantially enhanced the resolution of synteny and chromosomal colinearity between these genomes compared to previous reports. While a high degree of syntenic fragmentation was observed that was consistent with the large evolutionary distance between the L. angustifolius and M. truncatula genomes, there were striking examples of conserved colinearity of loci between these genomes. Compelling evidence was found of large-scale duplication and/or triplication in the L. angustifolius genome, consistent with one or more ancestral polyploidy events.     

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

Symbiotic diversity of Ensifer meliloti strains recovered from various legume species in Tunisia

Ensifer meliloti (formerly Sinorhizobium meliloti) was first considered as a specific microsymbiont of Medicago, Melilotus and Trigonella. However, strains of E. meliloti were recovered from root nodules of various legume species and their symbiotic status still remains unclear. Here, we further investigate the specificity of these strains. A collection of 47 E. meliloti strains isolated in Tunisia from root nodules of Medicago truncatula, Medicago sativa, Medicago ciliaris, Medicago laciniata, Medicago marina, Medicago scutellata, Phaseolus vulgaris, Cicer arietinum, Argyrolobium uniflorum, Lotus creticus, Lotus roudairei, Ononis natrix, Retama raetam, Genista saharae, Acacia tortilis, Hedysarum carnosum and Hippocrepis bicontorta were examined by REP-PCR fingerprinting, PCR-RFLPs of the 16S-23S rDNA IGS, the nifH gene and nifD-K intergenic spacer, and sequencing of 16S rRNA and nodA genes. Their nodulation range was also assessed by cross-inoculation experiments. No clear correlation was found between chromosomal backgrounds and host plants of origin. The nodulation polyvalence of the species E. meliloti was associated with a high symbiotic heterogeneity. On the basis of PCR-RFLP data from the nifH gene and nifD-K intergenic spacer, E. meliloti strains isolated from non-Medicago legumes harboured distinct genes and possessed wider host ranges. Some strains did not nodulate Medicago species. On the basis of nodA phylogeny, the majority of the Tunisian strains, including strains from Medicago, harboured distinct nodA alleles more related to those found in E. medicae than those found in E. meliloti. However, more work is still needed to characterize this group further. The diversity observed among M. laciniata isolates, which was supported by nodA phylogeny, nifH typing and the efficiency profile on M. ciliaris, indicated that what was thought to be bv. medicaginis is certainly heterogeneous.     

An Italian functional genomic resource for <Emphasis Type="Italic">Medicago truncatula</Emphasis>

Background

Medicago truncatula is a model species for legumes. Its functional genomics have been considerably boosted in recent years due to initiatives based both in Europe and US. Collections of mutants are becoming increasingly available and this will help unravel the genetic control of important traits for many species of legumes.

Findings

Our report is on the production of three complementary mutant collections of the model species Medicago truncatula produced in Italy in the frame of a national genomic initiative. Well established strategies were used: Tnt1 mutagenesis, TILLING and activation tagging. Both forward and reverse genetics screenings proved the efficiency of the mutagenesis approaches adopted, enabling the isolation of interesting mutants which are in course of characterization. We anticipate that the reported collections will be complementary to the recently established functional genomics tools developed for Medicago truncatula both in Europe and in the United States.

     

Recent Progress in Development of Tnt1 Functional Genomics Platform for Medicago truncatula and Lotus japonicus in Bulgaria

Legumes, as protein-rich crops, are widely used for human food, animal feed and vegetable oil production. Over the past decade, two legume species, Medicago truncatula and Lotus japonicus, have been adopted as model legumes for genomics and physiological studies. The tobacco transposable element, Tnt1, is a powerful tool for insertional mutagenesis and gene inactivation in plants. A large collection of Tnt1-tagged lines of M. truncatula cv. Jemalong was generated during the course of the project 'GLIP': Grain Legumes Integrated Project, funded by the European Union (www.eugrainlegumes.org). In the project 'IFCOSMO': Integrated Functional and COmparative genomics Studies on the MOdel Legumes Medicago truncatula and Lotus japonicus, supported by a grant from the Ministry of Education, Youth and Science, Bulgaria, these lines are used for development of functional genomics platform of legumes in Bulgaria. This review presents recent advances in the evaluation of the M. truncatula Tnt1 mutant collection and outlines the steps that are taken in using the Tnt1-tagging for generation of a mutant collection of the second model legume L. japonicus. Both collections will provide a number of legume-specific mutants and serve as a resource for functional and comparative genomics research on legumes. Genomics technologies are expected to advance genetics and breeding of important legume crops (pea, faba bean, alfalfa and clover) in Bulgaria and worldwide.     

Genome sequence of Ensifer medicae strain WSM1369; an effective microsymbiont of the annual legume Medicago sphaerocarpos

Ensifer medicae WSM1369 is an aerobic, motile, Gram-negative, non-spore-forming rod that can exist as a soil saprophyte or as a legume microsymbiont of Medicago. WSM1369 was isolated in 1993 from a nodule recovered from the roots of Medicago sphaerocarpos growing at San Pietro di Rudas, near Aggius in Sardinia (Italy). WSM1369 is an effective microsymbiont of the annual forage legumes M. polymorpha and M. sphaerocarpos. Here we describe the features of E. medicae WSM1369, together with genome sequence information and its annotation. The 6,402,557 bp standard draft genome is arranged into 307 scaffolds of 307 contigs containing 6,656 protein-coding genes and 79 RNA-only encoding genes. This rhizobial genome is one of 100 sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.     

Phosphorus deprivation affects composition and spatial distribution of membrane lipids in legume nodules

In legumes, symbiotic nitrogen (N) fixation (SNF) occurs in specialized organs called nodules after successful interactions between legume hosts and rhizobia. In a nodule, N-fixing rhizobia are surrounded by symbiosome membranes, through which the exchange of nutrients and ammonium occurs between bacteria and the host legume. Phosphorus (P) is an essential macronutrient, and N₂-fixing legumes have a higher requirement for P than legumes grown on mineral N. As in the previous studies, in P deficiency, barrel medic (Medicago truncatula) plants had impaired SNF activity, reduced growth, and accumulated less phosphate in leaves, roots, and nodules compared with the plants grown in P sufficient conditions. Membrane lipids in M. truncatula tissues were assessed using electrospray ionization–mass spectrometry. Galactolipids were found to increase in P deficiency, with declines in phospholipids (PL), especially in leaves. Lower PL losses were found in roots and nodules. Subsequently, matrix-assisted laser desorption/ionization–mass spectrometry imaging was used to spatially map the distribution of the positively charged phosphatidylcholine (PC) species in nodules in both P-replete and P-deficient conditions. Our results reveal heterogeneous distribution of several PC species in nodules, with homogeneous distribution of other PC classes. In P poor conditions, some PC species distributions were observed to change. The results suggest that specific PC species may be differentially important in diverse nodule zones and cell types, and that membrane lipid remodeling during P stress is not uniform across the nodule.

ESI–MS and matrix-assisted laser desorption ionization–mass spectrometry imaging reveal alterations in Medicago truncatula nodules membrane lipid composition and spatial distribution in phosphorus deficiency.