Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation 1. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Mya). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species 2. Medicago truncatula (Mt) is a long-established model for the study of legume biology. Here we describe the draft sequence of the Mt euchromatin based on a recently completed BAC-assembly supplemented with Illumina-shotgun sequence, together capturing ~94% of all Mt genes. A whole-genome duplication (WGD) approximately 58 Mya played a major role in shaping the Mt genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the Mt genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max (Gm) and Lotus japonicus (Lj). Mt is a close relative of alfalfa (M. sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the Mt genome sequence provides significant opportunities to expand alfalfa’s genomic toolbox.

Legumes are simultaneously one of the largest families of crop plants and a cornerstone in the biological nitrogen cycle. We combined molecular and phylogenetic analyses to evaluate genome conservation both within and between the two major clades of crop legumes. Genetic mapping of orthologous genes identifies broad conservation of genome macrostructure, especially within the galegoid legumes, while also highlighting inferred chromosomal rearrangements that may underlie the variation in chromosome number between these species. As a complement to comparative genetic mapping, we compared sequenced regions of the model legume Medicago truncatula with those of the diploid Lotus japonicus and the polyploid Glycine max. High conservation was observed between the genomes of M. truncatula and L. japonicus, whereas lower levels of conservation were evident between M. truncatula and G. max. In all cases, conserved genome microstructure was punctuated by significant structural divergence, including frequent insertion͞deletion of individual genes or groups of genes and lineage-specific expansion͞contraction of gene families. These results suggest that comparative mapping may have considerable utility for basic and applied research in the legumes, although its predictive value is likely to be tempered by phylogenetic distance and genome duplication.

Microsatellites are tandemly repeated short DNA sequences that are favored as molecular-genetic markers due to their high polymorphism index. Plant genomes characterized to date exhibit taxonspecific differences in frequency, genomic location, and motif structure of microsatellites, indicating that extant microsatellites originated recently and turn over quickly. With the goal of using microsatellite markers to integrate the physical and genetic maps of Medicago truncatula, we surveyed the frequency and distribution of perfect microsatellites in 77 Mbp of gene-rich BAC sequences, 27 Mbp of nonredundant transcript sequences, 20 Mbp of random whole genome shotgun sequences, and 49 Mbp of BAC-end sequences. Microsatellites are predominantly located in gene-rich regions of the genome, with a density of one long (i.e., $20 nt) microsatellite every 12 kbp, while the frequency of individual motifs varied according to the genome fraction under analysis. A total of 1,236 microsatellites were analyzed for polymorphism between parents of our reference intraspecific mapping population, revealing that motifs (AT) n , (AG) n , (AC) n , and (AAT) n exhibit the highest allelic diversity. A total of 378 genetic markers could be integrated with sequenced BAC clones, anchoring 274 physical contigs that represent 174 Mbp of the genome and composing an estimated 70% of the euchromatic gene space.