Comparative genomics to bridge Vicia faba with model and closely-related legume species: stability of QTLs for flowering and yield-related traits.

This study presents the development of an enhanced map in faba bean. The map contains 258 loci, mostly gene-based markers, organized in 16 linkage groups that expand 1,875 cM, with an average inter-marker distance of 7.26 cM. The combination of EST-derived markers with a number of markers physically located or previously ascribed to chromosomes by trisomic segregation, allowed the allocation of eight linkage groups (229 markers), to specific chromosomes. Moreover, this approach provided anchor points to establish a global homology among the faba bean chromosomes and those of closely-related legumes species. The map was used to identify and validate, for the first time, QTLs controlling five flowering and reproductive traits: days to flowering, flowering length, pod length, number of seeds per pod and number of ovules per pod. Twelve QTLs stable in the 2 years of evaluation were identified in chromosomes II, V and VI. Comparative mapping suggested the conservation of one of the faba bean genomic regions controlling the character days to flowering in other five legume species (Medicago, Lotus, pea, lupine, chickpea). Additional syntenic co-localizations of QTLs controlling pod length and number of seeds per pod between faba bean and Lotus japonicus are likely. The new genetic map opens the way for further translational studies between faba bean and related legume species, and provides an efficient tool for breeding applications such as QTL analysis and marker-assisted selection.

Validation of QTLs for Orobanche crenata resistance in faba bean (Vicia faba L.) across environments and generations.

Broomrape (Orobanche crenata Forsk.) is a major root-parasite of faba bean (Vicia faba L.), that seriously limits crop cultivation in the whole Mediterranean area. This parasitic weed is difficult to control, difficult to evaluate and the resistance identified so far is of polygenic nature. This study was conducted to identify genetic regions associated with broomrape resistance in recombinant inbred lines (RILs) and to validate their previous location in the original F(2) population derived from the cross between lines Vf6 and Vf136. A progeny consisting of 165 F(6) RILs was evaluated in three environments across two locations in 2003 and 2004. Two hundred seventy seven molecular markers were assigned to 21 linkage groups (9 of them assigned to specific chromosomes) that covered 2,856.7 cM of the V. faba genome. The composite interval mapping on the F(6) map detected more quantitative trait loci (QTL) than in the F(2) analysis. In this sense, four QTLs controlling O. crenata resistance (Oc2-Oc5) were identified in the RI segregant population in three different environments. Only Oc1, previously reported in the F(2) population, was not significant in the advanced lines. Oc2 and Oc3 were found to be associated with O. crenata resistance in at least two of the three environments, while the remaining two, Oc4 and Oc5, were only detected in Córdoba-04 and Mengíbar-04 and seemed to be environment dependent.

Detection of two quantitative trait loci for resistance to ascochyta blight in an intra-specific cross of chickpea (Cicer arietinum L.): development of SCAR markers associated with resistance.

Two quantitative trait loci (QTLs), (QTL(AR1) and QTL(AR2)) associated with resistance to ascochyta blight, caused by Ascochyta rabiei, have been identified in a recombinant inbred line population derived from a cross of kabulixdesi chickpea. The population was evaluated in two cropping seasons under field conditions and the QTLs were found to be located in two different linkage groups (LG4a and LG4b). LG4b was saturated with RAPD markers and four of them associated with resistance were sequenced to give sequence characterized amplified regions (SCARs) that segregated with QTL(AR2). This QTL explained 21% of the total phenotypic variation. However, QTL(AR1), located in LG4a, explained around 34% of the total phenotypic variation in reaction to ascochyta blight when scored in the second cropping season. This LG4a region only includes a few markers, the flower colour locus (B/b), STMS GAA47, a RAPD marker and an inter-simple-sequence-repeat and corresponds with a previously reported QTL. From the four SCARs tagging QTL(AR2), SCAR (SCY17(590)) was co-dominant, and the other three were dominant. All SCARs segregated in a 1:1 (presence:absence) ratio and the scoring co-segregated with their respective RAPD markers. QTL(AR2) on LG4b was mapped in a highly saturated genomic region covering a genetic distance of 0.8 cM with a cluster of nine markers (three SCARs, two sequence-tagged microsatellite sites (STMS) and four RAPDs). Two of the four SCARs showed significant alignment with genes or proteins related to disease resistance in other species and one of them (SCK13(603)) was sited in the highly saturated region linked to QTL(AR2). STMS TA72 and TA146 located in LG4b were described in previous maps where QTL for blight resistance were also localized in both inter and intraspecific crosses. These findings may improve the precision of molecular breeding for QTL(AR2) as they will allow the choice of as much polymorphism as possible in any population and could be the starting point for finding a candidate resistant gene for ascochyta blight resistance in chickpea.

Identification and characterization of NBS-LRR class resistance gene analogs in faba bean (Vicia faba L.) and chickpea (Cicer arietinum L.).

A PCR approach with degenerate primers designed from conserved NBS-LRR (nucleotide binding site-leucine-rich repeat) regions of known disease-resistance (R) genes was used to amplify and clone homologous sequences from 5 faba bean (Vicia faba) lines and 2 chickpea (Cicer arietinum) accessions. Sixty-nine sequenced clones showed homologies to various R genes deposited in the GenBank database. The presence of internal kinase-2 and kinase-3a motifs in all the sequences isolated confirm that these clones correspond to NBS-containing genes. Using an amino-acid sequence identity of 70% as a threshold value, the clones were grouped into 10 classes of resistance-gene analogs (RGA01 to RGA10). The number of clones per class varied from 1 to 30. RGA classes 1, 6, 8, and 9 were comprised solely of clones isolated from faba bean, whereas classes 2, 3, 4, 5, and 7 included only chickpea clones. RGA10, showing a within-class identity of 99%, was the only class consisting of both faba bean and chickpea clones. A phylogenetic tree, based on the deduced amino-acid sequences of 12 representative clones from the 10 RGA classes and the NBS domains of 6 known R genes (I2 and Prf from tomato, RPP13 from Arabidopsis, Gro1-4 from potato, N from tobacco, L6 from flax), clearly indicated the separation between TIR (Toll/interleukin-1 receptor homology: Gro1-4, L6, N, RGA05 to RGA10)- and non-TIR (I2, Prf, RPP13, RGA01 to RGA04)-type NBS-LRR sequences. The development of suitable polymorphic markers based on cloned RGA sequences to be used in genetic mapping will facilitate the assessment of their potential linkage relationships with disease-resistance genes in faba bean and chickpea. This work is the first to report on faba bean RGAs.

Genetic mapping of QTLs controlling horticultural traits in diploid roses.

A segregating progeny set of 96 F1 diploid hybrids (2n = 2x = 14) between "Blush Noisette" (D10), one of the first seedlings from the original "Champneys' Pink Cluster", and Rosa wichurana (E15), was used to construct a genetic linkage map of the rose genome following a "pseudo-testcross" mapping strategy. A total of 133 markers (130 RAPD, one morphological and two microsatellites) were located on the 14 linkage groups (LGs) of the D10 and E15 maps, covering total map lengths of 388 and 260 cM, respectively. Due to the presence of common biparental markers the homology of four LGs between parental maps (D10-1/E15-1 to D10-4/E15-4) could be inferred. Four horticulturally interesting quantitative traits, flower size (FS), days to flowering (DF), leaf size (LS), and resistance to powdery mildew (PM) were analysed in the progeny in order to map quantitative trait loci (QTLs) controlling these traits. A total of 13 putative QTLs (LOD > 3.0) were identified, four for FS, two for flowering time, five for LS, and two for resistance to PM. Possible homologies between QTLs detected in the D10 and E15 maps could be established between Fs1 and Fs3, Fs2 and Fs4, and Ls1 and Ls3. Screening for pairwise epistatic interactions between loci revealed additional, epistatic QTLs (EQTLs) for DF and LS that were not detected in the original QTL analysis. The genetic maps developed in this study will be useful to add new markers and locate genes for important traits in the genus providing a practical resource for marker-assisted selection programs in roses.

Cross-species amplification of Medicago truncatula microsatellites across three major pulse crops.

Model plants are facilitating the genetic characterization and comparative mapping of a number of traditional crops. Medicago truncatula has been widely accepted as a model plant to this end as it provides the essential tools for multiple aspects of legume genetics and genomics. A large set of markers from highly conserved M. truncatula gene regions is being created and used to establish a worldwide framework for comparative genomic studies in legumes. We have investigated the potential for cross-species amplification of 209 expressed sequence tag (EST)-based and 33 bacterial artificial chromosome (BAC)-based microsatellites from M. truncatula in the three most important European legume pulses-pea, faba bean and chickpea-that might facilitate future comparative mapping. Our results revealed significant transferability of M. truncatula microsatellites to the three pulses (40% in faba bean, 36.3% in chickpea and 37.6% in pea). The percentage of M. truncatula EST-SSRs (simple sequence repeats) amplified in the three crops (39-43%) was twofold higher than that of the genomic SSRs (21-24%). Sequence analysis determined that the level of conservation in the microsatellite motif was very low, while the flanking regions were generally well conserved. The variations in the sequences were mainly due to changes in the number of repeat motifs in the microsatellite region combined with indel and base substitutions. None of the functional microsatellites showed direct polymorphism among the parental genotypes tested, consequently preventing their immediate use for mapping purposes.

Interaction between Orobanche crenata and its host legumes: unsuccessful haustorial penetration and necrosis of the developing parasite.

BACKGROUND AND AIMS: Orobanche species represent major constraints to crop production in many parts of the world as they reduce yield and alter root/shoot allometry. Although much is known about the histology and effect of Orobanche spp. on susceptible hosts, less is known about the basis of host resistance to these parasites. In this work, histological aspects related to the resistance of some legumes to Orobanche crenata have been investigated in order to determine which types of resistance responses are involved in the unsuccessful penetration of O. crenata. METHODS: Samples of resistance reactions against O. crenata on different genotypes of resistant legumes were collected. The samples were fixed, sectioned and stained using different procedures. Sections were observed using a transmission light microscope and by epi-fluorescence. KEY RESULTS: Lignification of endodermal and pericycle host cells seems to prevent parasite intrusion into the root vascular cylinder at early infection stages. But in other cases, established tubercles became necrotic and died. Contrary to some previous studies, it was found that darkening at the infection site in these latter cases does not correspond to death of host tissues, but to the secretion of substances that fill the apoplast in the host-parasite interface and in much of the infected host tissues. The secretions block neighbouring host vessels. This may interfere with the nutrient flux between host and parasite, and may lead to necrosis and death of the developing parasite. CONCLUSIONS: The unsuccessful penetration of O. crenata seedlings into legume roots cannot be attributed to cell death in the host. It seems to be associated with lignification of host endodermis and pericycle cells at the penetration site. The accumulation of secretions at the infection site, may lead to the activation of xylem occlusion, another defence mechanism, which may cause further necrosis of established tubercles.

Development of a composite map in Vicia faba, breeding applications and future prospects.

A composite map of the Vicia faba genome based on morphological markers, isozymes, RAPDs, seed protein genes and microsatellites was constructed. The map incorporates data from 11 F(2) families for a total of 654 individuals all sharing the common female parent Vf 6. The integrated map is arranged in 14 major linkage groups (five of which were located in specific chromosomes). These linkage groups include 192 loci and cover 1559 cM with an overall average marker interval of 8 cM. By joining data of a new F(2) population segregating for resistance to ascochyta, broomrape and others traits of agronomic interest, have been saturated new areas of the genome. The combination of trisomic segregation, linkage analysis among loci from different families with a recurrent parent, and the analysis of new physically located markers, has allowed the establishment of the present status of the V. faba map with a wide coverage. This map provides an efficient tool in breeding applications such as disease-resistance mapping, QTL analyses and marker-assisted selection.

Phylogenetic analysis in the genus Cicer and cultivated chickpea using RAPD and ISSR markers.

Seventy five accessions belonging to 14 species of the genus Cicer were analysed with PCR-based molecular markers to determine their phylogenetic relationships. Eight of the species were annuals and included the Section Monocicer which contains cultivated chickpea ( Cicer arietinum L.). The remaining six species were perennials (five from Section Polycicer and one from Section Acanthocicer). More than one accession per species was analysed in most of the wild species; within C. arietinum, 26 accessions including Kabuli and Desi types, were studied. RAPD analyses using 12 primers gave 234 polymorphic fragments. Variability within species was detected. A dendrogram based on the Jaccard similarity index showed that the distribution pattern of variability between species was related to both growth habit and geographical origin. An accession of Cicer reticulatum was closer to accessions of Cicer echinospermum than to the four remaining of C. reticulatum, suggesting the possibility of gene flow between species. Cluster analysis for cultivated chickpea differentiated Kabuli and Desi types but we did not detect a clear relationship between groups and the geographical origin of the accessions.

Estimation of linkage in trisomic inheritance.

Based on F2 families derived from selfed F1 trisomic plants we have developed a genetic model to estimate linkage relationships between pairs of loci located on the extra chromosome. Genotypic frequencies of each class expected in a trisomic F2 family have been calculated and the maximum-likelihood equations for recombination-fraction estimation have been derived for a variety of genetic situations. Morton's test of homogeneity was used to compare recombination fractions estimated between loci exhibiting trisomic segregation to those obtained in families where the same loci showed Mendelian segregation. This method has been applied to an analysis of morphological, isozyme and RAPD data from faba bean (Vicia faba L.).

Genetic mapping of new morphological, isozyme and RAPD markers in Vicia faba L. using trisomics.

Thirteen F2 families of faba bean (Vicia faba L.), descended from plants trisomic for chromosomes 3, 4, 5 and 6, have been analyzed for morphological, isozyme and RAPD markers. This allowed the establishment of linkage relationships among these markers as well as the assignment of some markers and/or linkage groups to their respective chromosomes. The linkage analysis of partially overlapping sets of informative genetic markers for the data pooled from 13 F2 families has revealed 48 linkage groups, six of which have been precisely assigned to specific chromosomes. A statistical procedure to analyze the data of joint segregation analysis in families derived from trisomic plants has been developed.

Using RAPDs to study phylogenetic relationships in Rosa.

Nineteen species of rose (Rosa sp.) were analysed using Random Amplified Polymorphic DNA markers (RAPD). Each 10-base-long arbitrary primer produced a specific DNA banding pattern that grouped plants belonging to the same species and botanical sections as predicted from their genetic background. One hundred and seventy-five amplification products were examined by cluster analysis to assess the genetic relationships among species and their genetic distances. All of the accessions belonging to 1 species grouped together before branching to other species. Dendrograms constructed for intra- and inter-specific studies showed a good correlation with previous classifications by different authors based on morphological and cariological studies. Our results show that the RAPD technique is a sensitive and precise tool for genomic analysis in rose, being useful in assigning unclassified accessions to specific taxonomic groups or else allowing accessions classified by traditional criteria to be re-classified.

Variation of favism-inducing factors (vicine, convicine and L-DOPA) during pod development in Vicia faba L.

Changes in the concentrations of vicine, convicine and L-DOPA in two cultivars of Vicia faba L. seeds in different stages of pod development were determined by high performance liquid chromatography (HPLC). The vicine and convicine content was highest in fresh green cotyledons (moisture content about 80%) and gradually declined until a constant level was reached when seed dry matter percentage was around 40%. A similar pattern of variation in glucoside concentration was observed for the seed coat. The pods contained neither vicine nor convicine but they were particularly rich in L-DOPA. These compounds were not homogeneously distributed in the seeds.

Genetics and mapping of new isozyme loci in Vicia faba L using trisomics.

Polymorphism in ten enzyme systems (ACO, ACP, AAT, EST, FK, ME, NAG, PRX, 6PGD, and SOD) in Vicia faba L. was analyzed, revealing 13 loci, six of which have not been reported before. Inheritance, genetics, possible location, and linkage analysis were studied in 13 different F(2) populations trisomic for four of the six chromosomes (nos. 3, 4, 5, and 6) of the species. Each of these loci exhibited typical Mendelian inheritance except for those involved in the trisomic chromosome. Five loci have been assigned to a specific chromosome: Est-2 to chromosome 3, Fk-2 to chromosome 4, Prx-1 to chromosome 5, and Sod-1 and Pgd-p to chromosome 6. Nag-1 and Pgd-c displayed a linkage of 22.8 cM indicating a clear homology with chromosome 5 of garden pea on which both markers are syntenic.

Genetic changes under domestication in Vicia faba.

The components of variation within each one of two sets of landraces and/or cultivars of Vicia faba, respectively constituted of primitive and advanced morphological types, were studied by means of two sets of 8 × 8 diallel crosses with two repetitions. The results show that primitive and modern forms differ from each other in both the intensity and the kind of selective pressures acting on them, mainly on those characters more modified through the domestication process: i.e., seed morphology and the number of flowers per node. Because of the paramount importance of the additive component in the primitive forms, it is suggested that the most important type of selection on them is the stabilizing one. On the contrary, in the most advanced forms the selection is directional and disruptive : directional towards greater yields, and disruptive separating two morphological types, 'major' and 'equina'. The plant response to these different selective pressures has been to modify the genetic control of different characters: thus the primitive forms generally show only additivity while the most advanced forms show additivity as well as directional and asymmetrical dominance.

Genotype-environment interaction in tomato.

Twelve varieties of tomato of economic importance and their hybrids (including reciprocals) were studied in four environments: inside and outside of greenhouses and with and without plastic mulching. Seven characters were recorded per plant per environment: (1) total yield, (2) fruit weight, (3) locules/fruit, (4) fruits/cluster, (5) earliness in maturity, (6) earliness in harvesting and (7) leaves between clusters. There was an almost general tendency for hybrids to show higher values than the parentals for characters (1), (4) and (7); the opposite was true for (2) and (3), even when the differences were not statistically significant at the 5% level. Environments were always highly significant; the effect of the greenhouse explained most of the variation. Genotype-environment interaction by regression analysis showed that the performance of the hybrids was generally higher than that of the parents for characters (1), (4) and (7). Total yield was higher, in general, in the most protected environments. Locules per fruit was very constant but when interaction did exist, the number of locules was higher in the less protected environments. Hybrids interacted with environments more strongly than parent lines. Earliness was the most environmental dependent characteristic the choice of early harvesting being irrelevant. Some of the hybrids obtained seem promising from a commercial point of view.

Quantitative inheritance and barriers to crossability in Cicer arietinum L.

Three different diallel crosses were studied in Cicer arietinum; two of size 6×6, one within each of the two botanical groups 'macrosperma' and 'microsperma' of the cultivated subspecies, and one of 9×9 involving lines covering most of the morphological variation of chickpea. Barriers to crossability present neither a botanical nor a geographical pattern, being probably a direct consequence of interactions between genotypes. The genetic systems of twelve quantitative characters were analysed. Full dominance in a negative sense (small values dominant) is shown by leaflet length, width and shape index, rachis length, leaflet density on the rachis and pod length. Full dominance in a positive sense is shown by seeds per pod. Overdominance (in a positive sense) is evident for pods, seeds and yield per plant. Weak reciprocal differences were manifested by pod length, and pods, seeds and yield per plant. The system controlling number of leaflets per leaf is not clear. Dominance of 'primitive' over 'selected' characters seems to be the rule. As far as the environmental effects have permitted the analysis, no differences in genetic systems were observed between botanical groups.

On the evolution of Vicia faba L.

Diallel crosses between seven inbred lines obtained from different botanical groups including the two known subspecies of V. faba permitted study of the genetics of seventeen quantitative characteristics; there were very few significant reciprocal differences, and in all cases but one they were only weakly dependent on genotype. These results accord very well with those obtained from a study on population distances among thirty unselected populations of V. faba; many of these lie very close, forming a strong nucleus, which carries the maximum of potentialities of the species and from which different populations branch. Nevertheless differentiation occurs: the genetic system which regulates seed length shows overdominance in the positive sense in some lines, but partial or complete dominance in the negative sense in the rest of the studied lines. Interpretations suggesting that V. faba has suffered very little intraspecific differentiation are substantiated by the studies showing the presence of a partial incompatibility system; this is stronger in the Central European populations studied, weak (to various degrees) in the Spanish ones and absent in at least one population of the paucijuga group.

Evolutionary trends in Vicia faba L.

1. Distances have been calculated between thirty natural populations of Vicia faba. 2. The results strongly suggest the existence of a Central Nucleus from which different trends of evolution have arisen. 3. The Central Nucleus represents the equilibrium between different genetic tendencies; this equilibrium is postulated for the primitive populations of this species. 4. Disturbance of this balanced system can be produced by cycles of cross and self-fertilization; recombinations following these cycles are the origin of new genetic combinations so that a new equilibrium is reached with a new cycle of disturbance. 5. The Muratova classification is discussed. It is not possible to give a new classification based on these concepts, because of the great number of possible lines of evolution. The species must be considered as a whole, because there are no genetic barriers of isolation, only geographical ones and no sterility barriers have been produced between the subspecies. 6. A description of the Central Nucleus of Evolution is given.