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.


Background
The need for a more sustainable and environmentally safe agriculture has reinforced interest in the cultivation of legumes which include some of the most important agricultural species such as alfalfa, clover, pea, soybean, bean and peanut. These species have the capacity to establish atmospheric dinitrogen fixing symbiose with soil bacteria collectively named rhizobia, and to form symbiotic root mycorrhizae with soil fungi, thus facilitating their uptake of phosphate, water and other soil nutrients [1,2]. However, genetic analysis of these processes remains difficult in the major crop legumes due to features such as tetraploidy, large genomes and/or the lack of efficient methods for transgenesis. Since the model plant Arabidopsis thaliana, as indeed for other Cruciferacae, is unable to establish either rhizobial or mycorhizal symbioses, the need to establish a model legume has been recognized for over a decade. Furthermore, studying a model legume offers the opportunity to compare in the same plant both symbiotic and plantpathogen interactions and also to analyse plant physiological processes which cannot be satisfactorily studied in A. thaliana [3].
Furthermore, the genus Medicago is part of the Galegoid phylum and is therefore related to a number of important crop legumes in addition to alfalfa such as pea, faba bean, chickpea, lentil and clover [22]. Members of this phylum are expected to show a high level of nucleotide sequence conservation and similar genetic organization, and hence the potential for retively easy transfer of genome information between member species.
The diploid Lotus japonicus has also been proposed as a model legume [5,23,24]. This plant has similar attributes to those of M. truncatula but it is phylogenetically distant from the Galegoid phylum and other legumes such as soybean or bean.
Up to now, five genetic linkage maps have been constructed for either diploid M. sativa [25][26][27][28][29] or tetraploid M. sativa [30]. In this paper, we describe the first genetic map of the diploid model legume M. truncatula, an essential tool for genetics and genomics in this species. We report on the selection of two polymorphic M. truncatula genotypes providing the basis for classical and molecular genetics. An improved protocol was used to cross these two lines, leading to the construction of a genetic map of M. truncatula. Based on an F2 population of 124 plants, this map comprises 289 molecular markers and shows a low level of distortion in segregation and no clustering of markers. We took advantage of the existing high density genetic map of the diploid allogamous species Medicago sativa [29] to map a number of orthologous genes and to identify homologous linkage groups. These have been similarly numbered for the two Medicago species. Finally, we illustrate the use of these new tools to precisely map two genes of M. truncatula: the Mtsym6 gene [31] involved in strain x cultivar specificity of nitrogen fixation and the SPC gene, determining the direction of pod coiling [32].

Selection of parental genotypes
Since genetic analysis is based on the study of the segregation pattern of traits which differ between the two parents, two M. truncatula lines were selected on the basis of their phenotypic and molecular differences. M. truncatula is a self-fertilizing plant whose natural populations are composed of a variable number of different homozygous genotypes [33]. Cultivars can be genetically heterogeneous due either to a multiline selection process or to genetic pollution from other plant species during multiplication. In addition, in some natural populations, M. truncatula flowers occasionally cross-pollinate (<1%) [34,35]. For these reasons, single seeds from the cultivar Jemalong and from DZA315, an Algerian natural population (JMP, unpublished results) were selfed at least twice in order to ensure maximal homozygosity of the M. truncatula line. These two parental lines (Jemalong 6 and DZA315. 16) were named after their respective population or cultivar, followed by a number which identifies the particular line.
As stated above our choice of line DZA315.16 was based on the large number of morphological, developmental or symbiotic traits which are polymorphic with Jemalong (JMP and TH, unpuplished results). For instance, the two lines can be easily distinguished from each other by leaf pigmentation pattern: Jemalong 6 shows a typical pattern on the adaxial leaf surface whereas DZA315.16 leaves have only infrequent and randomly distributed spots (Figure 1A). Another striking difference lies in pod shape: Jemalong 6 pods have a spiny barrel shape characteristic of the M. truncatula ssp truncatula pods described by Lesins and Lesins [36], whereas DZA315.16 forms shorter and smaller pods and spines ( Figure 1B). Finally, the two lines show opposite coiling: anticlockwise for Jemalong 6 and clockwise for DZA315.16 ( Figure 1C).
The molecular diversity between the two lines was estimated by using anonymous AFLP banding patterns. We observed that 32% of AFLP bands are polymorphic between Jemalong 6 and DZA315.16 (not shown).
Another criteria for choosing these two lines was their identical DNA content since Blondon et al [7] have shown that M. truncatula genotypes can differ by up to 20% in their DNA content. The DNA content of the two parental lines was measured by flow cytometry and shown to be identical (1.16 pg/2C), corresponding approximately to a haploid genome size of 580 Mb. According to Arumuganathan and Earle [37], as well as Bennett et al [38], this value is similar to Lotus japonicus (1.0 pg/2C, [23]) and about 4 times greater than the 145 Mb/1 C of Arabidopsis thaliana (0.30 pg/2 C). It is worth noting that the diploid M. sativa (800 Mb/1 C, [7]) possesses a significantly higher DNA content than M. truncatula.
In studies carried out over the last decade, a number of other single-seed descent lines of M. truncatula have been described. Three of these are derived from the Jemalong cultivar. The Jemalong line A17 (TH, unpublished results) has been used for the construction of most cDNA and BAC libraries as well as for mutagenesis [39] and the Jemalong line J5 was used in a γ-ray mutagenesis program [12]. So far we have failed to observe any difference between the Jemalong lines A17, J5 and J6 despite the use of more than 4000 AFLP markers (TH and JMA, data not shown). Furthermore, we have never been able to identify any phenotypic (morphological or symbiotic) difference between these three Jemalong lines. Our conclusion is that lines, A17, J5 and J6 can be considered as having an identical genotype, even if very limited differences in nucleotide sequence cannot be ruled out. On the other hand, Rose et al [40] have isolated, after in vitro selection, a derivative, 2HA, from Jemalong with high regeneration capacity. In addition to Jemalong, Penmetsa and Cook [41] proposed A20, a line selected from natural populations (TH, unpublished results), for genetic studies. Similarly, the R108-1 (c3) line, isolated from cell culture selection and possessing a higher capacity for regeneration has also been proposed as a model genotype [17,20,42]. Since all these four M. truncatula genotypes (Jemalong, DZA315.16, A20 and R108-1 (c3) possess distinctive traits, it is clear that there is considerable potential in exploiting the natural genetic variability of M. truncatula for revealing plant developmental processes.

Construction of the genetic map
The 124 individuals of the mapping population were genotyped for 313 markers comprising, 292 dominant anonymous markers (72 RAPD and 220 AFLP), 19 genes with known functions (Table 3) and 2 codominant isoenzyme markers ( Table 3). The 72 polymorphic RAPD bands were generated using a set of 18 10-mer primers ( Table 1, 4 bands generated per primer on average). All the 13 AFLP primer combinations (Table 2) were shown to be polymorphic and generated an average of 16.9 markers/primer combination.
The M. truncatula genetic map was built using the MAP-MAKER software [43] with a minimum LOD score value of 5. Three hundred and four markers were organized in 8 linkage groups, corresponding to the haploid chromosome set of M. truncatula (2n = 16). Only 9 markers (2 RAPD and 7 AFLP markers) remained unlinked. Five pairs of RAPD and ten pairs of AFLP dominant markers were transformed into, respectively, 5 and 10 codominant markers when the following conditions were observed: the two bands should be in repulsing phase, have a simi-   (Table 4). Figure 3 shows the distribution diagram of interval distance between two adjacent markers with an average value of 4.4 cM and a standard deviation of 4.3 cM; 90% of the markers are located within an interval of less than 10 cM. According to the number, origin and nature of the markers their distribution is approximately similar between the linkage groups (Table 5).
A serious problem in creating F2 genetic maps comes from the presence of alternating markers in coupling and repulsion phase. To circumvent this problem, we separated the markers which are dominant into two sets, grouping separately either male (DZA315. 16) or female (Jemalong 6) dominant markers with codominant markers ( Table 6). A dominant marker is considered as a male (or female) marker if the recessive allelic form (absence) is male (or female). Male and female genetic maps were built that span 888 and 979 Kosambi cM respectively (Table 7). Since the order of codominant markers is not modified between the two maps (data not shown), we present here the easier-to-read global F2 representation ( Figure 2) with 289 markers. Table 6, about 27% of the markers do not show a typical Mendelian segregation ratio at α = 0.05. It is noteworthy that there is a high level of deviation from the theoretical ratio only in the case of 3 linkage groups. 59, 46 and 88% of the markers in linkage groups 1, 2 and 3, respectively, are distorted, together representing 86% of all the distorted markers with high χ 2 values for 1:3 segregation. On the other hand, each of the 5 other linkage groups has less than 3 moderately distorted markers (data not shown). Furthermore these distorted markers are not randomly distributed throughout chromosomes: for example, a nearly linear gradient of distortion can be observed on linkage group 3 ( Figure 4) ranging from Mendelian equality of male and female markers to a 5fold frequency in favor of male alleles. On a whole genome basis, the frequency of distorted male markers compared to female markers appears to be similar: 22% and 29% respectively ( Table 6).

Alignment of M. sativa and M. truncatula genetic maps
The agronomically important legume Medicago sativa (alfalfa) is taxonomically very close to M. truncatula. However, the colinearity of the two genomes has not been evaluated. For this reason, we mapped 18 gene or isoenzyme markers (Table 3, Figure 2) which tag the 8 linkage groups of the diploid M. sativa genetic map [29]. All these markers were similarly linked in the two species with the exception of the rDNA which was found on linkage group 5 of M. truncatula, compared to linkage group 6 of M. sativa. To facilitate comparison, we gave the same numbering to homologous linkage groups. The relative orientation of

Figure 2
Global F2 genetic map of Medicago truncatula. The number above the linkage groups refers to the homologous linkage group in M. sativa [29]. The code to the right of the linkage groups refers to the marker name. The numbers to the left of the linkage groups refers to the genetic distances (Kosambi cM) from the top and have been rounded up for clarity. The sign + indicates that two RAPD markers have been transformed into a codominant marker. In the case of codominant AFLP markers, the codes of the two bands are given together. Stars refer to known genes used for synteny studies. Circles refer to known genes not used for synteny studies.

Mapping of symbiotic (Mtsym6) and developmental (SPC) genes
Many polymorphisms exist between Jemalong 6 and DZA315. 16. One symbiotic and one developmental mo-nogenic trait were selected and genetically mapped to illustrate the use of our genetic map.
Tirichine et al [31] showed that a Fix -/Fix + polymorphism is present between the Jemalong 6 and DZA315.16 parents respectively following inoculation with the wild type Sinorhizobium meliloti A 145 strain. This polymorphism is controlled by a single recessive gene, named Mtsym6: the Fix + trait being dominant. In order to precisely map this gene, we determined the nitrogen-fixing genotype of 73 F2 plants of the mapping population by inoculating their  [31]. Symbiotic typing data together with molecular marker segregation allowed us to map the Mtsym6 gene at less than 1 cM (LOD > 10) from the cosegregating PI 12 and L12200 markers on linkage group 8 ( Figure 2).

Selection of parental lines
In order to construct a genetic map for M. truncatula, we selected two polymorphic lines, one from an Australian cultivar and the other from a natural Algerian population. These two lines are very polymorphic and have traits which allow them to be readily distinguished. For instance, Jemalong leaf markings are almost absent in the DZA315.16 line. Penmetsa and Cook [41] reported the existence of a single gene responsible for the leaf spots between Jemalong and A20 lines.
The cultivar Jemalong was the first M. truncatula cultivar to be released [44] and is widely used in agricultural leyfarming [45]. The cultivar Jemalong was initially selected for molecular studies because of its transformation/regeneration characteristics [4,19,46]. So far, most of the Medicago programs including mutagenesis, the construction of large scale cDNA and genomic libraries, DNA sequencing and genetic resources [39] have made use of Jemalong lines. A major interest of the Jemalong line is the fact that  LG = M. truncatula Linkage Group cM = Kosambi centimorgan SD = Standard Deviation  it can be nodulated by the Sinorhizobium meliloti strain 1021, whose genome has been recently sequenced [11] [http://sequence.toulouse.inra.fr/meliloti.html] .
The percentage of AFLP bands which are polymorphic between Jemalong and DZA315.16 (32%), is similar to the 30 to 34% observed between 3 ecotypes of Arabidopsis thaliana [47]. In addition to the limited level of marker distortion, this high % of polymorphism is convenient for generating genetic maps in the model legume M. truncatula.

The genetic map
We have built a genetic map using a combination of anonymous dominant markers and known genes of co-dominant evaluation. The reliability and value of this genetic map is demonstrated in the following ways: i) identical number of linkage groups and chromosomes as well as the small number of unlinked markers (2.8%) at a minimal LOD 5 ; ii) no discrepancies between the global F2 map and both male and female maps; iii) only a limited clustering of markers on a few linkage groups. The excellent distribution of markers is illustrated by the fact that 90% of the markers are located within an interval of less LG LG = M. truncatula Linkage Group LG = M. truncatula Linkage Group A dominant marker is considered as a male marker if the recessive allelic form is male (and the same for female). , it is now well known that the genetic size of genomes is not proportional to haploid DNA content. This leads to different values for the average ratio of physical to genetic distance (kb/cM): 263 kb/cM for Arabidopsis thaliana [47], 3000 kb/cM for Pisum sativum [48] 1900 kb/cM for Helianthus annuus [49] and 1000 kb/cM for diploid Medicago sativa [29].

Number of markers Dominant markers
An F2 genetic map is difficult to construct using dominant markers because recombination between markers can be estimated with confidence only when the markers are in coupling phase. For this reason, we established two additional maps containing the codominant markers, as anchor markers, combined with either male or female markers in coupling phase. The comparison of these two maps to the global F2 map gives a reasonable estimate of the quality of the genetic map. Except for a shorter genetic length due to an underestimation of crossing overs, the estimated distances between codominant markers and their relative order is very similar in the two maps.
The 8 linkage groups of M. truncatula differ in genetic size by a factor of 2.4 from 92 to 219 cM, whereas M. truncatula chromosome length differs by a factor of 1.5 (3 to 4.5 µm) or 2.3 (29 to 68 µm) for metaphasic [50] or pachytene [51] chromosomes respectively. No obvious correlation could be found between genetic size of linkage groups and pachytene chromosome length measurements [51]. Finally, the range of genetic size of the M. truncatula linkage groups suggests that, on average, 1 to 2 crossingover occur per chromosome per meiosis.

Segregation distortion of markers and quality of the genetic map
About 27% of the markers do not have the expected Mendelian ratios (at α = 0.05). This can be compared with the 40% value observed by Jenczewski et al [52] in M. truncatula × M. tornata intraspecific crosses and to the 15 to 50% value for loci having distorted segregation in diploid alfalfas [26,27,29]. High frequencies of segregation distortion in the M. sativa genome can be explained by the exposure of deleterious recessive alleles during inbreeding which caused gametic and/or zygotic selection as hypothesized by Brouwer and Osborn [30]. Interestingly, the distorted markers are not scattered over the M. truncatula genome but concentrated on linkage groups I, II and III suggesting structural reasons for this distortion.
From the analysis carried out by Jenczewski et al [52], it is clear that variations from expected Mendelian ratios are common within both inter-and intraspecific crosses. Distorting factors can be deleterious recessive alleles [53], self-incompatibility alleles [54], structural rearrangements [55] or differences in DNA content [52]. In the case of Jemalong 6 and DZA315.16, the two genotypes have a similar overall DNA content. However, we do not know if DNA content is identical for homologous chromosomes of the two lines. However, the gradient in distortion in favor of male alleles along linkage group 3 suggests the presence of specific gene(s) interfering with meiosis. Whatever the cause, a consequence of segregation distortion is that chromosomes, or at least parts of chromosomes, are not transmitted equally to the progeny. It would be of interest to know if the segregation distortion frequency depends on the lines involved in the cross and if it is always limited to the same 3 linkage groups. Cytogenetics experiments are in progress to identify the linkage between meiotic pairing and segregation distortion. The average genetic distances between adjacent markers for male and female linkage groups are not significantly different (P = 0.95). cM = Kosambi Centimorgan SD = Standard Deviation

Macrosynteny between Medicago species
Despite different genome sizes and evolutionary divergence, a remarkable degree of genome conservation has already been revealed by comparative genetic mapping experiments, as for instance in the Poaceae family [56]. We took advantage of the availability of a dense genetic map with more than 850 markers in diploid M. sativa [29], including many known genes, to address the question of macrosynteny between M. truncatula and alfalfa. At least two distant markers per chromosome were selected on the alfalfa map and mapped onto the M. truncatula segregation population. These results allowed identification of homologous linkage groups between the two species and consequently to name the M. truncatula linkage groups following the M. sativa nomenclature, as already stated in Kalo et al [29]. This nomenclature has been further extended to chromosome nomenclature [51]. The only observed discrepancy was the different localization of rDNA in the two species, on linkage group 5 for M. truncatula instead of group 6 for M. sativa [29]. The localisation of rDNA on linkage group 5 is ascertained from the LOD score for linkage between ENOD40 and rDNA and Lb1 and rDNA of, respectively, 7 and 10. Moreover, FISH experiments show that rDNA is localized on the same chromosome than ENOD40 [51]. The estimation of the degree of genome colinearity between these two closely related legume species (as well as with more distant legumes) is presently in progress. This should open the way to use M. truncatula as a nodal plant for a identifying legume genes of interest, characterized either after mutagenesis or within natural diversity. Unfortunately, due to its heterozygous character, alfalfa cannot be easily used either in mutagenesis or in genomic programs, thus reinforcing the interest in a closely related model legume.  [31]. Here, we precisely map this gene to linkage group 8. This is the first report of the genetic mapping of a symbiotic gene in the model legume M. truncatula. The map-based cloning of Mtsym6 is currently under progress in the laboratory.
Anticlockwise pod coiling is a characteristic of most Medicago species [36]. Bena et al [57] have shown that only a single monophyletic clade contains species (including M. truncatula), where pods can be coiled in either orientation, suggesting that a single change must have led to the clockwise coiling character. Lilienfeld and Kihara [32] have demonstrated that, in Medicago littoralis (a close relative of M. truncatula), the clockwise direction is dominant and under monogenic control. We took advantage of the fact that Jemalong 6 and DZA315.16 are polymorphic for this character, having anticlockwise and clockwise pods respectively, to make a genetic analysis of this trait within our mapping population. In agreement with Lilienfeld and Kihara [32], we reached the conclusion that pod coiling sense has a monogenic determinism with clockwise turning dominant. We named this gene SPC and mapped it to linkage group 7. The map-based cloning of SPC is currently under progress in the laboratory.

Conclusions
The identification of polymorphic Jemalong 6 and DZA315.16 lines together with the making of a F2 genetic map and the identification of synteny with alfalfa, described in this report are likely to provide a powerful tool for both fundamental and applied approaches, such as studying the conservation of synteny between genomes of legumes and other model species, and the positional cloning of agriculturally important genes.

Plant material
Individual seeds of the M. truncatula cultivar Jemalong and the Algerian natural population DZA315 (JMP, unpublished results) were selfed at least twice and then crossed manually with Jemalong being the female parent. A population of 124 F2 plants, derived from a single F1 plant, was used for genetic mapping and is designated thereafter as the "mapping population". The hybrid nature of the F1 plant was confirmed with molecular markers. More than 95% of F2 seeds germinated, out of which about 6% gave albino plantlets that were discarded. One plant from each parental line, one F1 plant and the F2 individuals of the mapping population were cultivated axenically in Magenta boxes on SHb10 medium [58] and maintained by cuttings every two months.

Crossing procedure
An efficient manual crossing procedure has been developped based on a microscopic study of pollination and fertilization kinetics in relation to the intact flower morphology of M. truncatula (E.-P. Journet, unpublished). This method is similar to method 3 described in Pathipanawat et al [59] and is outlined here. Parental plants are vernalized at 4°C for 15 days after germination and cultivated in a growth chamber (25°C, 16 h light photoperiod, 100 µE.m -2 .sec -1 ). For pollen sampling, flowers at the optimal stage (just past anthesis) are selected on the male parent and dissected under a stereomicroscope using curved extra fine forceps. At this stage, released pollen appears turgescent, moist and sticky and is packed around the stigma. On the female parent, flowers are selected at a developmental stage just prior to anther bursting and pollen release. The standard petal is longitudinally incised 0.5 mm below its central line using a scalpel. The 10 non-open anther bags are then carefully removed, and the tip of a freshly harvested sexual column applied to the stigma in order to saturate its sticky surface with exogenous pollen. The pollinated pistil is then gently placed back under the standard petal. The branch tip bearing the cross-pollinated flower is labelled, inserted into a ~25-ml clear plastic vial containing 1 ml water and gently secured with a cotton wool plug. Female plants are kept under indirect light until the protecting vial is removed 24-48 h later. The success of crosses is indicated by the development of small coiling pods that become visible 2-4 days after pollination. The efficiency of this method is close to 80 % on average with 5 viable seeds per pod. Significant variations in the optimal stage for crossing and in the rate of success have been observed for different M. truncatula genotypes and appear to depend on the characteristics of the maternal genotype and genetic distance between parental lines. A more detailed crossing protocol is available from EPJ upon request.

DNA content measurements
DNA content was determined using an EPICS V flow cytometer at the Service de Cytometrie, Institut des Sciences Végétales, CNRS, Gif-sur-Yvette, France. Measurements were assessed as described by Blondon et al [7] using ethidium bromide fluorescence and Petunia hybrida PxPC6 (2C = 2.85 pg) as the routine internal standard. Each genome size estimation resulted from five independent measurements for each single plant.

DNA isolation and marker typing
DNA isolation and RAPD amplification were as described by Ghérardi et al [60]. RAPD markers were named as follows: the letter and the first two digits refer to the identification of the 10-mer primer (Table 1) from Operon Kits (Operon Technologies, Alameda, Calif, USA) and the last digits correspond to the molecular weight of the polymorphic bands.
AFLP markers were produced essentially as described by Vos et al [61] as modified by Moreau et al [62] with the EcoRI and MseI restriction enzymes and using, respectively, a +2/+3 combination. Polymorphic bands were identi-fied by the code of the Eco/Mse primer combinations ( Table 2) followed by a number corresponding to the band number in the parental profile. RAPD and AFLP markers were considered as dominant markers on a presence/absence alternative. The intensity of the bands was not taken into account.
Among the 19 genes used for mapping (Table 3) two genes were typed as RFLP markers. RFLP analysis was performed as follows: from each individual of the mapping population, DNA was isolated from a pool of more than 10 F3 plants, digested with restriction enzymes (Table 3) and Southern blotted. DNA probes (Table 3) were obtained from the respective authors, and labelled and hybridized as described in Kalo et al [29]. The polymorphism of the 17 other genes was identified after PCR amplification with primers anchored either on both sides of a microsatellite motif (SSR) or within exons and crossing introns (data not shown). When a length polymorphism could not be observed, the amplification products were digested with restriction enzymes (CAPS, [63]). RFLP, CAPS and PCR markers were scored as codominant.

Isozyme markers
Fresh pieces of actively growing young leaf tissue were ground in a cold 0.1 M Tris-HCl (pH = 7.2) extraction buffer. The supernatant was absorbed onto filter paper wicks, kept frozen at -80°C overnight, and then loaded onto 13% starch gels. The 6-phosphogluco dehydrogenase (PGD) and the phosphogluco-mutase (PGM) enzyme isoforms were resolved with a Tris-citrate (pH = 7.0) buffer system as described by Chaulet [34]. Banding intensity was used to identify heterozygotes. Isozymes were scored as codominant markers.

Map construction
The genetic map was constructed using the MapMaker/ Exp V3.0 software [43] with the Kosambi map function and was drawn with the Drawmap package [64]. The segregation data were tested for deviation from the expected Mendelian ratio using a Chi-square test. The markers showing significant deviation at α = 0.05 were nevertheless used to build the linkage map.

Gene mapping
For the Mtsym6 gene, each F2 plant of the mapping population was scored for its symbiotic phenotype by testing a minimum of 10 F3 plants inoculated with Sinorhizobium meliloti strain A145 as described by Tirichine et al [31]. The nitrogen fixation phenotype (Fix) of F3 individuals was scored 30 days after inoculation and treated as a codominant marker. Thus the heterozygous F2 were identified by the presence of both Fix+ and Fix-in their progeny. For the SPC gene, each F2 plant of the mapping population was scored for pod coiling (clockwise or anticlockwise) and the marker was treated as dominant.