SSR analysis of the Medicago truncatula SARDI core collection reveals substantial diversity and unusual genotype dispersal throughout the Mediterranean basin

Medicago truncatula quantitative resistance to a new strain of Verticillium alfalfae from Iran revealed by a genome-wide association study

Verticillium wilt is a major threat to many crops, among them alfalfa (Medicago sativa). The model plant Medicago truncatula, a close relative of alfalfa was used to study the genetic control of resistance towards a new Verticillium alfalfae isolate. The accidental introduction of pathogen strains through global trade is a threat to crop production and such new strains might also be better adapted to global warming. Isolates of V. alfalfae were obtained from alfalfa fields in Iran and characterized. The Iranian isolate AF1 was used in a genome-wide association study (GWAS) involving 242 accessions from the Mediterranean region. Root inoculations were performed with conidia at 25°C and symptoms were scored regularly. Maximum Symptom Score and Area under Disease Progess Curve were computed as phenotypic traits to be used in GWAS and for comparison to a previous study with French isolate V31.2 at 20°C. This comparison showed high correlation with a shift to higher susceptibility, and similar geographical distribution of resistant and susceptible accessions to AF1 at 25°C, with resistant accessions mainly in the western part. GWAS revealed 30 significant SNPs linked to resistance towards isolate AF1. None of them were common to the previous study with isolate V31.2 at 20°C. To confirm these loci, the expression of nine underlying genes was studied. All genes were induced in roots following inoculation, in susceptible and resistant plants. However, in resistant plants induction was higher and lasted longer. Taken together, the use of a new pathogen strain and a shift in temperature revealed a completely different genetic control compared to a previous study that demonstrated the existence of two major QTLs. These results can be useful for Medicago breeding programs to obtain varieties better adapted to future conditions.

Plastid phylogenomics uncovers multiple species in Medicago truncatula (Fabaceae) germplasm accessions

Medicago truncatula is a model legume that has been extensively investigated in diverse subdisciplines of plant science. Medicago littoralis can interbreed with M. truncatula and M. italica; these three closely related species form a clade, i.e. TLI clade. Genetic studies have indicated that M. truncatula accessions are heterogeneous but their taxonomic identities have not been verified. To elucidate the phylogenetic position of diverse M. truncatula accessions within the genus, we assembled 54 plastid genomes (plastomes) using publicly available next-generation sequencing data and conducted phylogenetic analyses using maximum likelihood. Five accessions showed high levels of plastid DNA polymorphism. Three of these highly polymorphic accessions contained sequences from both M. truncatula and M. littoralis. Phylogenetic analyses of sequences placed some accessions closer to distantly related species suggesting misidentification of source material. Most accessions were placed within the TLI clade and maximally supported the interrelationships of three subclades. Two Medicago accessions were placed within a M. italica subclade of the TLI clade. Plastomes with a 45-kb (rpl20-ycf1) inversion were placed within the M. littoralis subclade. Our results suggest that the M. truncatula accession genome pool represents more than one species due to possible mistaken identities and gene flow among closely related species.

Genetic Diversity and Differentiation of Eleven Medicago Species from Campania Region Revealed by Nuclear and Chloroplast Microsatellites Markers

The species belonging to the genus Medicago are considered a very important genetic resource at global level both for planet’s food security and for sustainable rangelands management. The checklist of the Italian flora (2021) includes a total number of 40 Medicago species for Italy, and 27 for Campania region, with a number of doubtful records or related to species no more found in the wild. In this study, 10 Medicago species native to Campania region, and one archaeophyte (M. sativa), identified by means of morphological diagnostic characters, were analyzed in a blind test to assay the efficacy of nine microsatellite markers (five cp-SSRs and four n-SSRs). A total number of 33 individuals from 6 locations were sampled and genotyped. All markers were polymorphic, 40 alleles were obtained with n-SSRs ranging from 8–12 alleles per locus with an average of 10 alleles per marker, PIC values ranged from 0.672 to 0.847, and the most polymorphic SSR was MTIC 564. The cp-SSRs markers were highly polymorphic too; PIC values ranged from 0.644 to 0.891 with an average of 0.776, the most polymorphic cp-SSR was CCMP10. 56 alleles were obtained with cp-SSRs ranging from 7 to 17 alleles per locus with an average of 11. AMOVA analysis with n-SSR markers highlighted a great level of genetic differentiation among the 11 species, with a statistically significant fixation index (FST). UPGMA clustering and Bayesian-based population structure analysis assigned these 11 species to two main clusters, but the distribution of species within clusters was not the same for the two analyses. In conclusion, our results demonstrated that the combination of the used SSRs well distinguished the 11 Medicago species. Moreover, our results demonstrated that the use of a limited number of SSRs might be considered for further genetic studies on other Medicago species.

Phoma medicaginis Isolate Differences Determine Disease Severity and Phytoestrogen Production in Annual Medicago spp

Phoma black stem and leaf spot disease of annual Medicago spp., caused by Phoma medicaginis, not only can devastate forage and seed yield but can reduce herbage quality by inducing production of phytoestrogens (particularly coumestrol and 4'-O-methylcoumestrol), which can also reduce the ovulation rates of animals grazing infected forage. We determined the consequent phytoestrogen levels on three different annual Medicago species/cultivars (Medicago truncatula cultivar Cyprus, Medicago polymorpha var. brevispina cultivar Serena, and Medicago murex cultivar Zodiac) after inoculation with 35 isolates of P. medicaginis. Across the isolate × cultivar combinations, leaf disease incidence, petiole/stem disease incidence, leaf disease severity, petiole disease severity, and leaf yellowing severity ranged up to 100, 89.4, 100, 58.1, and 61.2%, respectively. Cultivars Cyprus and Serena were the most susceptible and cultivar Zodiac was the most resistant to P. medicaginis. Isolates WAC3653, WAC3658, and WAC4252 produced the most severe disease. Levels of phytoestrogens in stems ranged from 25 to 1,995 mg/kg for coumestrol and from 0 to 418 mg/kg for 4'-O-methylcoumestrol. There was a significant positive relationship of disease incidence and severity parameters with both coumestrol and 4'-O-methylcoumestrol contents, as noted across individual cultivars and across the three cultivars overall, where r = 0.39 and 0.37 for coumestrol and 4'-O-methylcoumestrol, respectively (P < 0.05). Although cultivar Serena was most susceptible to P. medicaginis and produced the highest levels of phytoestrogens in the presence of P. medicaginis, cultivar Zodiac contained the highest levels of phytoestrogens in comparison with other cultivars in the absence of P. medicaginis. There was a 15-fold increase in coumestrol in cultivar Serena but only a 7-fold increase in cultivar Zodiac from infection of P. medicaginis. The study highlights that the intrinsic ability of a particular cultivar to produce phytoestrogens in the absence of the pathogen, and its comparative ability to produce phytoestrogens in the presence of the P. medicaginis, are both important and highly relevant to developing new annual Medicago spp. cultivars that offer improved disease resistance and better animal reproductive outcomes.

Candidate genes and SNPs associated with stomatal conductance under drought stress in Vitis

BACKGROUND

Understanding the complexity of the vine plant's response to water deficit represents a major challenge for sustainable winegrowing. Regulation of water use requires a coordinated action between scions and rootstocks on which cultivars are generally grafted to cope with phylloxera infestations. In this regard, a genome-wide association study (GWAS) approach was applied on an 'ad hoc' association mapping panel including different Vitis species, in order to dissect the genetic basis of transpiration-related traits and to identify genomic regions of grape rootstocks associated with drought tolerance mechanisms. The panel was genotyped with the GrapeReSeq Illumina 20 K SNP array and SSR markers, and infrared thermography was applied to estimate stomatal conductance values during progressive water deficit.

RESULTS

In the association panel the level of genetic diversity was substantially lower for SNPs loci (0.32) than for SSR (0.87). GWAS detected 24 significant marker-trait associations along the various stages of drought-stress experiment and 13 candidate genes with a feasible role in drought response were identified. Gene expression analysis proved that three of these genes (VIT_13s0019g03040, VIT_17s0000g08960, VIT_18s0001g15390) were actually induced by drought stress. Genetic variation of VIT_17s0000g08960 coding for a raffinose synthase was further investigated by resequencing the gene of 85 individuals since a SNP located in the region (chr17_10,497,222_C_T) was significantly associated with stomatal conductance.

CONCLUSIONS

Our results represent a step forward towards the dissection of genetic basis that modulate the response to water deprivation in grape rootstocks. The knowledge derived from this study may be useful to exploit genotypic and phenotypic diversity in practical applications and to assist further investigations.

Ethylene Is Not Essential for R-Gene Mediated Resistance but Negatively Regulates Moderate Resistance to Some Aphids in Medicago truncatula

Ethylene is important for plant responses to environmental factors. However, little is known about its role in aphid resistance. Several types of genetic resistance against multiple aphid species, including both moderate and strong resistance mediated by R genes, have been identified in Medicago truncatula. To investigate the potential role of ethylene, a M. truncatula ethylene- insensitive mutant, sickle, was analysed. The sickle mutant occurs in the accession A17 that has moderate resistance to Acyrthosiphon kondoi, A. pisum and Therioaphis trifolii. The sickle mutant resulted in increased antibiosis-mediated resistance against A. kondoi and T. trifolii but had no effect on A. pisum. When sickle was introduced into a genetic background carrying resistance genes, AKR (A. kondoi resistance), APR (A. pisum resistance) and TTR (T. trifolii resistance), it had no effect on the strong aphid resistance mediated by these genes, suggesting that ethylene signaling is not essential for their function. Interestingly, for the moderate aphid resistant accession, the sickle mutant delayed leaf senescence following aphid infestation and reduced the plant biomass losses caused by both A. kondoi and T. trifolii. These results suggest manipulation of the ethylene signaling pathway could provide aphid resistance and enhance plant tolerance against aphid feeding.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

Medicago truncatula: Genetic and Genomic Resources

Medicago truncatula was chosen by the legume community, along with Lotus japonicus, as a model plant to study legume biology. Since then, numerous resources and tools have been developed for M. truncatula. These include, for example, its genome sequence, core ecotype collections, transformation/regeneration methods, extensive mutant collections, and a gene expression atlas. This review aims to describe the different genetic and genomic tools and resources currently available for M. truncatula. We also describe how these resources were generated and provide all the information necessary to access these resources and use them from a practical point of view. © 2017 by John Wiley & Sons, Inc.

Is there genetic variation in mycorrhization of Medicago truncatula?

Differences in the plant's response among ecotypes or accessions are often used to identify molecular markers for the respective process. In order to analyze genetic diversity of Medicago truncatula in respect to interaction with the arbuscular mycorrhizal (AM) fungus Rhizophagus irregularis, mycorrhizal colonization was evaluated in 32 lines of the nested core collection representing the genetic diversity of the SARDI collection. All studied lines and the reference line Jemalong A17 were inoculated with R. irregularis and the mycorrhization rate was determined at three time points after inoculation. There were, however, no reliable and consistent differences in mycorrhization rates among all lines. To circumvent possible overlay of potential differences by use of the highly effective inoculum, native sandy soil was used in an independent experiment. Here, significant differences in mycorrhization rates among few of the lines were detectable, but the overall high variability in the mycorrhization rate hindered clear conclusions. To narrow down the number of lines to be tested in more detail, root system architecture (RSA) of in vitro-grown seedlings of all lines under two different phosphate (Pi) supply condition was determined in terms of primary root length and number of lateral roots. Under high Pi supply (100 µM), only minor differences were observed, whereas in response to Pi-limitation (3 µM) several lines exhibited a drastically changed number of lateral roots. Five lines showing the highest alterations or deviations in RSA were selected and inoculated with R. irregularis using two different Pi-fertilization regimes with either 13 mM or 3 mM Pi. Mycorrhization rate of these lines was checked in detail by molecular markers, such as transcript levels of RiTubulin and MtPT4. Under high phosphate supply, the ecotypes L000368 and L000555 exhibited slightly increased fungal colonization and more functional arbuscules, respectively. To address the question, whether capability for mycorrhizal colonization might be correlated to general invasion by microorganisms, selected lines were checked for infection by the root rot causing pathogen, Aphanoymces euteiches. The mycorrhizal colonization phenotype, however, did not correlate with the resistance phenotype upon infection with two strains of A. euteiches as L000368 showed partial resistance and L000555 exhibited high susceptibility as determined by quantification of A. euteiches rRNA within infected roots. Although there is genetic diversity in respect to pathogen infection, genetic diversity in mycorrhizal colonization of M. truncatula is rather low and it will be rather difficult to use it as a trait to access genetic markers.

Genetic Mapping of a Major Resistance Gene to Pea Aphid (Acyrthosipon pisum) in the Model Legume Medicago truncatula

Resistance to the Australian pea aphid (PA; Acyrthosiphon pisum) biotype in cultivar Jester of the model legume Medicago truncatula is mediated by a single dominant gene and is phloem-mediated. The genetic map position for this resistance gene, APR (Acyrthosiphon pisum resistance), is provided and shows that APR maps 39 centiMorgans (cM) distal of the A. kondoi resistance (AKR) locus, which mediates resistance to a closely related species of the same genus bluegreen aphid (A. kondoi). The APR region on chromosome 3 is dense in classical nucleotide binding site leucine-rich repeats (NLRs) and overlaps with the region harbouring the RAP1 gene which confers resistance to a European PA biotype in the accession Jemalong A17. Further screening of a core collection of M. truncatula accessions identified seven lines with strong resistance to PA. Allelism experiments showed that the single dominant resistance to PA in M. truncatula accessions SA10481 and SA1516 are allelic to SA10733, the donor of the APR locus in cultivar Jester. While it remains unclear whether there are multiple PA resistance genes in an R-gene cluster or the resistance loci identified in the other M. truncatula accessions are allelic to APR, the introgression of APR into current M. truncatula cultivars will provide more durable resistance to PA.

Naturally occurring diversity helps to reveal genes of adaptive importance in legumes

Environmental changes challenge plants and drive adaptation to new conditions, suggesting that natural biodiversity may be a source of adaptive alleles acting through phenotypic plasticity and/or micro-evolution. Crosses between accessions differing for a given trait have been the most common way to disentangle genetic and environmental components. Interestingly, such man-made crosses may combine alleles that never meet in nature. Another way to discover adaptive alleles, inspired by evolution, is to survey large ecotype collections and to use association genetics to identify loci of interest. Both of these two genetic approaches are based on the use of biodiversity and may eventually help us in identifying the genes that plants use to respond to challenges such as short-term stresses or those due to global climate change. In legumes, two wild species, Medicago truncatula and Lotus japonicus, plus the cultivated soybean (Glycine max) have been adopted as models for genomic studies. In this review, we will discuss the resources, limitations and future plans for a systematic use of biodiversity resources in model legumes to pinpoint genes of adaptive importance in legumes, and their application in breeding.

Genome variations account for different response to three mineral elements between Medicago truncatula ecotypes Jemalong A17 and R108

BACKGROUND

Resequencing can be used to identify genome variations underpinning many morphological and physiological phenotypes. Legume model plant Medicago truncatula ecotypes Jemalong A17 (J. A17) and R108 differ in their responses to mineral toxicity of aluminum and sodium, and mineral deficiency of iron in growth medium. The difference may result from their genome variations, but no experimental evidence supports this hypothesis.

RESULTS

A total of 12,750 structure variations, 135,045 short insertions/deletions and 764,154 single nucleotide polymorphisms were identified by resequencing the genome of R108. The suppressed expression of MtAACT that encodes a putative aluminum-induced citrate efflux transporter by deletion of partial sequence of the second intron may account for the less aluminum-induced citrate exudation and greater accumulation of aluminum in roots of R108 than in roots of J. A17, thus rendering R108 more sensitive to aluminum toxicity. The higher expression-level of MtZpt2-1 encoding a TFIIIA-related transcription factor in J. A17 than R108 under conditions of salt stress can be explained by the greater number of stress-responsive elements in its promoter sequence, thus conferring J. A17 more tolerant to salt stress than R108 plants by activating the expression of downstream stress-responsive genes. YSLs (Yellow Stripe-Likes) are involved in long-distance transport of iron in plants. We found that an YSL gene was deleted in the genome of R108 plants, thus rendering R108 less tolerance to iron deficiency than J. A17 plants.

CONCLUSIONS

The deletion or change in several genes may account for the different responses of M. truncatula ecotypes J. A17 and R108 to mineral toxicity of aluminum and sodium as well as iron deficiency. Uncovering genome variations by resequencing is an effective method to identify different traits between species/ecotypes that are genetically related. These findings demonstrate that analyses of genome variations by resequencing can shed important light on differences in responses of M. truncatula ecotypes to abiotic stress in general and mineral stress in particular.

Medicago truncatula ecotypes A17 and R108 differed in their response to iron deficiency

Medicago truncatula Gaertn is a model legume species with a wide genetic diversity. To evaluate the responses of the two M. truncatula ecotypes, the effect of Fe deficiency on ecotype A17 and ecotype R108, which have been widely used in physiological and molecular studies, was investigated. A greater reduction in shoot Fe concentration of R108 plants than that of A17 plants was observed under Fe-deficient conditions. Exposure to Fe-deficient medium led to a greater increase in ferric chelate reductase (FCR) activity in roots of A17 than those of R108 plants, while expression of genes encoding FCR in roots of A17 and R108 plants was similarly up-regulated by Fe deficiency. Exposure of A17 plants to Fe-deficient medium evoked an ethylene evolution from roots, while the same treatment had no effect on ethylene evolution from R108 roots. There was a significant increase in expression of MtIRT encoding a Fe transporter in A17, but not in R108 plants, upon exposure to Fe-deficient medium. Transcripts of MtFRD3 that is responsible for loading of iron chelator citrate into xylem were up-regulated by Fe deficiency in A17, but not in R108 plants. These results suggest that M. truncatula ecotypes A17 and R108 differed in their response and adaptation to Fe deficiency, and that ethylene may play an important role in regulation of greater tolerance of A17 plant to Fe deficiency. These findings provide important clues for further elucidation of molecular mechanism by which legume plants respond and adapt to low soil Fe availability.

Plant-aphid interactions with a focus on legumes

Sap-sucking insects such as aphids cause substantial yield losses in agriculture by draining plant nutrients as well as vectoring viruses. The main method of control in agriculture is through the application of insecticides. However, aphids rapidly evolve mechanisms to detoxify these, so there is a need to develop durable plant resistance to these damaging insect pests. The focus of this review is on aphid interactions with legumes, but work on aphid interactions with other plants, particularly Arabidopsis and tomato is also discussed. This review covers advances on the plant side of the interaction, including the identification of major resistance genes and quantitative trait loci conferring aphid resistance in legumes, basal and resistance gene mediated defence signalling following aphid infestation and the role of specialised metabolites. On the aphid side of the interaction, this review covers what is known about aphid effector proteins and aphid detoxification enzymes. Recent advances in these areas have provided insight into mechanisms underlying resistance to aphids and the strategies used by aphids for successful infestations and have significant impacts for the delivery of durable resistance to aphids in legume crops.

The population structure and diversity of eggplant from Asia and the Mediterranean Basin

A collection of 238 eggplant breeding lines, heritage varieties and selections within local landraces provenanced from Asia and the Mediterranean Basin was phenotyped with respect to key plant and fruit traits, and genotyped using 24 microsatellite loci distributed uniformly throughout the genome. STRUCTURE analysis based on the genotypic data identified two major sub-groups, which to a large extent mirrored the provenance of the entries. With the goal to identify true-breeding types, 38 of the entries were discarded on the basis of microsatellite-based residual heterozygosity, along with a further nine which were not phenotypically uniform. The remaining 191 entries were scored for a set of 19 fruit and plant traits in a replicated experimental field trial. The phenotypic data were subjected to principal component and hierarchical principal component analyses, allowing three major morphological groups to be identified. All three morphological groups were represented in both the “Occidental” and the “Oriental” germplasm, so the correlation between the phenotypic and the genotypic data sets was quite weak. The relevance of these results for evolutionary studies and the further improvement of eggplant are discussed. The population structure of the core set of germplasm shows that it can be used as a basis for an association mapping approach.

MtQRRS1, an R-locus required for Medicago truncatula quantitative resistance to Ralstonia solanacearum

Ralstonia solanacearum is a major soilborne pathogen that attacks > 200 plant species, including major crops. To characterize MtQRRS1, a major quantitative trait locus (QTL) for resistance towards this bacterium in the model legume Medicago truncatula, genetic and functional approaches were combined. QTL analyses together with disease scoring of heterogeneous inbred families were used to define the locus. The candidate region was studied by physical mapping using a bacterial artificial chromosome (BAC) library of the resistant line, and sequencing. In planta bacterial growth measurements, grafting experiments and gene expression analysis were performed to investigate the mechanisms by which this locus confers resistance to R. solanacearum. The MtQRRS1 locus was localized to the same position in two recombinant inbred line populations and was narrowed down to a 64 kb region. Comparison of parental line sequences revealed 15 candidate genes with sequence polymorphisms, but no evidence of differential gene expression upon infection. A role for the hypocotyl in resistance establishment was shown. These data indicate that the quantitative resistance to bacterial wilt conferred by MtQRRS1, which contains a cluster of seven R genes, is shared by different accessions and may act through intralocus interactions to promote resistance.

Genotype delimitation in the Nod-independent model legume Aeschynomene evenia

Research on the nitrogen-fixing symbiosis has been so far focused on two model legumes, Medicago truncatula and Lotus japonicus, which use a sophisticated infection process involving infection thread formation. However, in 25% of the legumes, the bacterial entry occurs more simply in an intercellular fashion. Among them, some semi-aquatic Aeschynomene species present the distinctive feature to form nitrogen-fixing nodules on both roots and stems following elicitation by photosynthetic bradyrhizobia that do not produce Nod factors. This interaction is believed to represent a living testimony of the ancestral state of the rhizobium-legume symbiosis. To decipher the molecular mechanisms of this unique Nod-independent nitrogen-fixing symbiosis, we previously identified A. evenia C. Wright as an appropriate model legume, because it displays all the requisites for molecular and genetic approaches. To advance the use of this new model legume species, here we characterized the intraspecific diversity found in A. evenia. For this, the accessions available in germplasm banks were collected and subjected to morphological investigations, genotyping with RAPD and SSR markers, molecular phylogenies using ITS and single nuclear gene sequences, and cross-compatibility tests. These combined analyses revealed an important intraspecific differentiation that led us to propose a new taxonomic classification for A. evenia comprising two subspecies and four varieties. The A. evenia ssp. evenia contains var. evenia and var. pauciciliata whereas A. evenia ssp. serrulata comprises var. serrulata and var. major. This study provides information to exploit efficiently the diversity encountered in A. evenia and proposes subsp. evenia as the most appropriate subspecies for future projects aimed at identifying plant determinants of the Nod-independent symbiotic process.

Genetic diversity and population structure assessed by SSR and SNP markers in a large germplasm collection of grape

BACKGROUND

The economic importance of grapevine has driven significant efforts in genomics to accelerate the exploitation of Vitis resources for development of new cultivars. However, although a large number of clonally propagated accessions are maintained in grape germplasm collections worldwide, their use for crop improvement is limited by the scarcity of information on genetic diversity, population structure and proper phenotypic assessment. The identification of representative and manageable subset of accessions would facilitate access to the diversity available in large collections. A genome-wide germplasm characterization using molecular markers can offer reliable tools for adjusting the quality and representativeness of such core samples.

RESULTS

We investigated patterns of molecular diversity at 22 common microsatellite loci and 384 single nucleotide polymorphisms (SNPs) in 2273 accessions of domesticated grapevine V. vinifera ssp. sativa, its wild relative V. vinifera ssp. sylvestris, interspecific hybrid cultivars and rootstocks. Despite the large number of putative duplicates and extensive clonal relationships among the accessions, we observed high level of genetic variation. In the total germplasm collection the average genetic diversity, as quantified by the expected heterozygosity, was higher for SSR loci (0.81) than for SNPs (0.34). The analysis of the genetic structure in the grape germplasm collection revealed several levels of stratification. The primary division was between accessions of V. vinifera and non-vinifera, followed by the distinction between wild and domesticated grapevine. Intra-specific subgroups were detected within cultivated grapevine representing different eco-geographic groups. The comparison of a phenological core collection and genetic core collections showed that the latter retained more genetic diversity, while maintaining a similar phenotypic variability.

CONCLUSIONS

The comprehensive molecular characterization of our grape germplasm collection contributes to the knowledge about levels and distribution of genetic diversity in the existing resources of Vitis and provides insights into genetic subdivision within the European germplasm. Genotypic and phenotypic information compared in this study may efficiently guide further exploration of this diversity for facilitating its practical use.

Phoma medicaginis stimulates the induction of the octadecanoid and phenylpropanoid pathways in Medicago truncatula

Gene expression changes and metabolite abundances were measured during the interaction of Medicago truncatula with the fungal necrotrophic pathogen Phoma medicaginis in leaf tissue of susceptible and resistant accessions. Over 330 genes were differentially expressed in plants infected with P. medicaginis relative to mock-inoculated plants at 12 h post-inoculation. Of these, 191 were induced in either the resistant or the susceptible accession, with 143 genes repressed. Expression changes were observed in genes involved in the oxidative burst, cell wall strengthening and lipid metabolism, as well as several transcription factors. Genes related to salicylic acid, jasmonate and ethylene responses were up-regulated, as well as genes leading to the production of jasmonic acid. Significant induction of genes in the phenylpropanoid pathway leading to lignin and isoflavonoid biosynthesis occurred. High-pressure liquid chromatography with UV detection (HPLC-UV) identified several phenolic compounds induced by P. medicaginis, as well as constitutively higher levels of phenolic compounds, in the resistant M. truncatula accession. Differentially regulated genes induced in both the resistant and susceptible accessions, but with different kinetics, and constitutively more highly expressed and induced phenolic compounds provide candidates for functional analysis. Taken together, these results highlight the importance of the octadecanoid and phenylpropanoid pathways in defence against this necrotrophic pathogen.

Identification and characterization of resistance to cowpea aphid (Aphis craccivora Koch) in Medicago truncatula

BACKGROUND

Cowpea aphid (CPA; Aphis craccivora) is the most important insect pest of cowpea and also causes significant yield losses in other legume crops including alfalfa, beans, chickpea, lentils, lupins and peanuts. In many of these crops there is no natural genetic resistance to this sap-sucking insect or resistance genes have been overcome by newly emerged CPA biotypes.

RESULTS

In this study, we screened a subset of the Medicago truncatula core collection of the South Australian Research and Development Institute (SARDI) and identified strong resistance to CPA in a M. truncatula accession SA30199, compared to all other M. truncatula accessions tested. The biology of resistance to CPA in SA30199 plants was characterised compared to the highly susceptible accession Borung and showed that resistance occurred at the level of the phloem, required an intact plant and involved a combination of antixenosis and antibiosis. Quantitative trait loci (QTL) analysis using a F2 population (n = 150) from a cross between SA30199 and Borung revealed that resistance to CPA is controlled in part by a major quantitative trait locus (QTL) on chromosome 2, explaining 39% of the antibiosis resistance.

CONCLUSIONS

The identification of strong CPA resistance in M. truncatula allows for the identification of key regulators and genes important in this model legume to give effective CPA resistance that may have relevance for other legume crops. The identified locus will also facilitate marker assisted breeding of M. truncatula for increased resistance to CPA and potentially other closely related Medicago species such as alfalfa.

Simple sequence repeat polymorphisms (SSRPs) for evaluation of molecular diversity and germplasm classification of minor crops

Evaluation of the genetic diversity among populations is an essential prerequisite for the preservation of endangered species. Thousands of new accessions are introduced into germplasm institutes each year, thereby necessitating assessment of their molecular diversity before elimination of the redundant genotypes. Of the protocols that facilitate the assessment of molecular diversity, SSRPs (simple sequence repeat polymorphisms) or microsatellite variation is the preferred system since it detects a large number of DNA polymorphisms with relatively simple technical complexity. The paucity of information on DNA sequences has limited their widespread utilization in the assessment of genetic diversity of minor or neglected crop species. However, recent advancements in DNA sequencing and PCR technologies in conjunction with sophisticated computer software have facilitated the development of SSRP markers in minor crops. This review examines the development and molecular nature of SSR markers, and their utilization in many aspects of plant genetics and ecology.

Population genetic structure of two Medicago species shaped by distinct life form, mating system and seed dispersal

BACKGROUND AND AIMS

Life form, mating system and seed dispersal are important adaptive traits of plants. In the first effort to characterize in detail the population genetic structure and dynamics of wild Medicago species in China, a population genetic study of two closely related Medicago species, M. lupulina and M. ruthenica, that are distinct in these traits, are reported. These species are valuable germplasm resources for the improvement of Medicago forage crops but are under threat of habitat destruction.

METHODS

Three hundred and twenty-eight individuals from 16 populations of the annual species, M. lupulina, and 447 individuals from 15 populations of the perennial species, M. ruthenica, were studied using 15 and 17 microsatellite loci, respectively. Conventional and Bayesian-clustering analyses were utilized to estimate population genetic structure, mating system and gene flow.

KEY RESULTS

Genetic diversity of M. lupulina (mean H(E)=0.246) was lower than that of M. ruthenica (mean H(E)=0.677). Populations of M. lupulina were more highly differentiated (F(ST)=0.535) than those of M. ruthenica (F(ST)=0.130). For M. lupulina, 55.5% of the genetic variation was partitioned among populations, whereas 76.6% of the variation existed within populations of M. ruthenica. Based on the genetic data, the selfing rates of M. lupulina and M. ruthenica were estimated at 95.8% and 29.5%, respectively. The genetic differentiation among populations of both species was positively correlated with geographical distance.

CONCLUSIONS

The mating system differentiation estimated from the genetic data is consistent with floral morphology and observed pollinator visitation. There was a much higher historical gene flow in M. ruthenica than in M. lupulina, despite more effective seed dispersal mechanisms in M. lupulina. The population genetic structure and geographical distribution of the two Medicago species have been shaped by life form, mating systems and seed dispersal mechanisms.

Morphological and microsatellite diversity associated with ecological factors in natural populations of Medicago laciniata Mill. (Fabaceae)

Genetic variability in 10 natural Tunisian populations of Medicago laciniata were analysed using 19 quantitative traits and 12 polymorphic microsatellite loci. A large degree of genetic variability within-populations and among-populations was detected for both quantitative characters and molecular markers. High genetic differentiation among populations for quantitative traits was seen, with Q(ST) = 0.47, and F(ST) = 0.47 for microsatellite markers. Several quantitative traits displayed no statistical difference in the levels of Q(ST) and F(ST). Further, significant correlations between quantitative traits and eco-geographical factors suggest that divergence in the traits among populations may track environmental differences. There was no significant correlation between genetic variability at quantitative traits and microsatellite markers within populations. The site-of-origin of eco-geographical factors explain between 18.13% and 23.40% of genetic variance among populations at quantitative traits and microsatellite markers, respectively. The environmental factors that most influence variation in measured traits among populations are assimilated phosphorus (P(2)0(5)) and mean annual rainfall, followed by climate and soil texture, altitude and organic matter. Significant associations between eco-geographical factors and gene diversity, He, were established in five microsatellite loci suggesting that these simple sequence repeats (SSRs) are not necessarily biologically neutral.

The use of neutral and non-neutral SSRs to analyse the genetic structure of a Tunisian collection of Medicago truncatula lines and to reveal associations with eco-environmental variables

In this study, we investigated the genetic diversity of a collection of 136 Medicago truncatula lines from 10 Tunisian natural populations collected in well-defined locations and in various ecological conditions of soil, salinity and water availability. The genetic diversity was evaluated using a set of 18 microsatellites (SSRs), representing the 8 chromosomes of M. truncatula. A neutrality test showed that 7 SSRs were non-neutral with evidence of balancing selection. The 11 neutral SSRs revealed a geographical pooling with the Tunisian Dorsale axis restricting migration of alleles. The 7 non-neutral alleles demonstrate a correlation with rainfall, altitude and salinity environmental variables suggesting that these SSRs are linked to genes involved in water use efficiency, resistance to salinity or adaptation to altitude, and that there is local adaptation of M. truncatula to these variables. This demonstrates that the choice of so-called neutral markers should be carefully evaluated in population genetic studies. This study illustrates the genetic diversity occurring in natural Tunisian populations of M. truncatula and describes the first collection of this species dedicated to natural variation involved in adaptation to the environment.

Construction of a comparative genetic map in faba bean (Vicia faba L.); conservation of genome structure with Lens culinaris

Background

The development of genetic markers is complex and costly in species with little pre-existing genomic information. Faba bean possesses one of the largest and least studied genomes among cultivated crop plants and no gene-based genetic maps exist. Gene-based orthologous markers allow chromosomal regions and levels of synteny to be characterised between species, reveal phylogenetic relationships and chromosomal evolution, and enable targeted identification of markers for crop breeding. In this study orthologous codominant cross-species markers have been deployed to produce the first exclusively gene-based genetic linkage map of faba bean (Vicia faba), using an F₆ population developed from a cross between the lines Vf6 (equina type) and Vf27 (paucijuga type).

Results

Of 796 intron-targeted amplified polymorphic (ITAP) markers screened, 151 markers could be used to construct a comparative genetic map. Linkage analysis revealed seven major and five small linkage groups (LGs), one pair and 12 unlinked markers. Each LG was comprised of three to 30 markers and varied in length from 23.6 cM to 324.8 cM. The map spanned a total length of 1685.8 cM. A simple and direct macrosyntenic relationship between faba bean and Medicago truncatula was evident, while faba bean and lentil shared a common rearrangement relative to M. truncatula. One hundred and four of the 127 mapped markers in the 12 LGs, which were previously assigned to M. truncatula genetic and physical maps, were found in regions syntenic between the faba bean and M. truncatula genomes. However chromosomal rearrangements were observed that could explain the difference in chromosome numbers between these three legume species. These rearrangements suggested high conservation of M. truncatula chromosomes 1, 5 and 8; moderate conservation of chromosomes 2, 3, 4 and 7 and no conservation with M. truncatula chromosome 6. Multiple PCR amplicons and comparative mapping were suggestive of small-scale duplication events in faba bean. This study also provides a preliminary indication for finer scale macrosynteny between M. truncatula, lentil and faba bean. Markers originally designed from genes on the same M. truncatula BACs were found to be grouped together in corresponding syntenic areas in lentil and faba bean.

Conclusion

Despite the large size of the faba bean genome, comparative mapping did not reveal evidence for polyploidisation, segmental duplication, or significant rearrangements compared to M. truncatula, although a bias in the use of single locus markers may have limited the detection of duplications. Non-coding repetitive DNA or transposable element content provides a possible explanation for the difference in genome sizes. Similar patterns of rearrangements in faba bean and lentil compared to M. truncatula support phylogenetic studies dividing these species into the tribes Viceae and Trifoliae. However, substantial macrosynteny was apparent between faba bean and M. truncatula, with the exception of chromosome 6 where no orthologous markers were found, confirming previous investigations suggesting chromosome 6 is atypical. The composite map, anchored with orthologous markers mapped in M. truncatula, provides a central reference map for future use of genomic and genetic information in faba bean genetic analysis and breeding.

Detection of QTLs for flowering date in three mapping populations of the model legume species Medicago truncatula

Adaptation to the environment and reproduction are dependent on the date of flowering in the season. The objectives of this paper were to evaluate the effect of photoperiod on flowering date of the model species for legume crops, Medicago truncatula and to describe genetic architecture of this trait in multiple mapping populations. The effect of photoperiod (12 and 18 h) was analysed on eight lines. Quantitative variation in three recombinant inbred lines (RILs) populations involving four parental lines was evaluated, and QTL detection was carried out. Flowering occurred earlier in long than in short photoperiods. Modelling the rate of progression to flowering with temperature and photoperiod gave high R2, with line-specific parameters that indicated differential responses of the lines to both photoperiod and temperature. QTL detection showed a QTL on chromosome 7 that was common to all populations and seasons. Taking advantage of the multiple mapping populations, it was condensed into a single QTL with a support interval of only 0.9 cM. In a bioanalysis, six candidate genes were identified in this interval. This design also indicated other genomic regions that were involved in flowering date variation more specifically in one population or one season. The analysis on three different mapping populations detected more QTLs than on a single population, revealed more alleles and gave a more precise position of the QTLs that were common to several populations and/or seasons. Identification of candidate genes was a result of integration of QTL analysis and genomics in M. truncatula.

Medicago truncatula as a model for understanding plant interactions with other organisms, plant development and stress biology: past, present and future

Medicago truncatula Gaertn. cv. Jemalong, a pasture species used in Australian agriculture, was first proposed as a model legume in 1990. Since that time M. truncatula, along with Lotus japonicus (Regal) Larsen, has contributed to major advances in understanding rhizobia Nod factor perception and the signalling pathway involved in nodule formation. Research using M. truncatula as a model has expanded beyond nodulation and the allied mycorrhizal research to investigate interactions with insect pests, plant pathogens and nematodes. In addition to biotic stresses the genetic mechanisms to ameliorate abiotic stresses such as salinity and drought are being investigated. Furthermore, M. truncatula is being used to increase understanding of plant development and cellular differentiation, with nodule differentiation providing a different perspective to organogenesis and meristem biology. This legume plant represents one of the major evolutionary success stories of plant adaptation to its environment, and it is particularly in understanding the capacity to integrate biotic and abiotic plant responses with plant growth and development that M. truncatula has an important role to play. The expanding genomic and genetic toolkit available with M. truncatula provides many opportunities for integrative biological research with a plant which is both a model for functional genomics and important in agricultural sustainability.

Two alternative recessive quantitative trait loci influence resistance to spring black stem and leaf spot in Medicago truncatula

BACKGROUND

Knowledge of the genetic basis of plant resistance to necrotrophic pathogens is incomplete and has been characterised in relatively few pathosystems. In this study, the cytology and genetics of resistance to spring black stem and leaf spot caused by Phoma medicaginis, an economically important necrotrophic pathogen of Medicago spp., was examined in the model legume M. truncatula.

RESULTS

Macroscopically, the resistant response of accession SA27063 was characterised by small, hypersensitive-like spots following inoculation while the susceptible interaction with accessions A17 and SA3054 showed necrotic lesions and spreading chlorosis. No unique cytological differences were observed during early infection (<48 h) between the resistant and susceptible genotypes, except pathogen growth was restricted to one or a few host cells in SA27063. In both interactions reactive oxygen intermediates and phenolic compounds were produced, and cell death occurred. Two F2 populations segregating for resistance to spring black stem and leaf spot were established between SA27063 and the two susceptible accessions, A17 and SA3054. The cross between SA27063 and A17 represented a wider cross than between SA27063 and SA3054, as evidenced by higher genetic polymorphism, reduced fertility and aberrant phenotypes of F2 progeny. In the SA27063 x A17 F2 population a highly significant quantitative trait locus (QTL, LOD = 7.37; P < 0.00001) named resistance to the necrotroph Phoma medicaginis one (rnpm1) genetically mapped to the top arm of linkage group 4 (LG4). rnpm1 explained 33.6% of the phenotypic variance in the population's response to infection depicted on a 1-5 scale and was tightly linked to marker AW256637. A second highly significant QTL (LOD = 6.77; P < 0.00001), rnpm2, was located on the lower arm of LG8 in the SA27063 x SA3054 map. rnpm2 explained 29.6% of the phenotypic variance and was fine mapped to a 0.8 cM interval between markers h2_16a6a and h2_21h11d. rnpm1 is tightly linked to a cluster of Toll/Interleukin1 receptor-nucleotide binding site-leucine-rich repeat (TIR-NBS-LRR) genes and disease resistance protein-like genes, while no resistance gene analogues (RGAs) are apparent in the genomic sequence of the reference accession A17 at the rnpm2 locus.

CONCLUSION

The induction of defence responses and cell death in the susceptible interaction following infection by P. medicaginis suggested this pathogen is not negatively affected by these responses and may promote them. A QTL for resistance was revealed in each of two populations derived from crosses between a resistant accession and two different susceptible accessions. Both loci are recessive in nature, and the simplest explanation for the existence of two separate QTLs is the occurrence of host genotype-specific susceptibility loci that may interact with undetermined P. medicaginis virulence factors.

Recent Advances in Medicago truncatula Genomics

Legume rotation has allowed a consistent increase in crop yield and consequently in human population since the antiquity. Legumes will also be instrumental in our ability to maintain the sustainability of our agriculture while facing the challenges of increasing food and biofuel demand. Medicago truncatula and Lotus japonicus have emerged during the last decade as two major model systems for legume biology. Initially developed to dissect plant-microbe symbiotic interactions and especially legume nodulation, these two models are now widely used in a variety of biological fields from plant physiology and development to population genetics and structural genomics. This review highlights the genetic and genomic tools available to the M. truncatula community. Comparative genomic approaches to transfer biological information between model systems and legume crops are also discussed.

The first genetic and comparative map of white lupin (Lupinus albus L.): identification of QTLs for anthracnose resistance and flowering time, and a locus for alkaloid content

We report the first genetic linkage map of white lupin (Lupinus albus L.). An F8 recombinant inbred line population developed from Kiev mutant x P27174 was mapped with 220 amplified fragment length polymorphism and 105 gene-based markers. The genetic map consists of 28 main linkage groups (LGs) that varied in length from 22.7 cM to 246.5 cM and spanned a total length of 2951 cM. There were seven additional pairs and 15 unlinked markers, and 12.8% of markers showed segregation distortion at P < 0.05. Syntenic relationships between Medicago truncatula and L. albus were complex. Forty-five orthologous markers that mapped between M. truncatula and L. albus identified 17 small syntenic blocks, and each M. truncatula chromosome aligned to between one and six syntenic blocks in L. albus. Genetic mapping of three important traits: anthracnose resistance, flowering time, and alkaloid content allowed loci governing these traits to be defined. Two quantitative trait loci (QTLs) with significant effects were identified for anthracnose resistance on LG4 and LG17, and two QTLs were detected for flowering time on the top of LG1 and LG3. Alkaloid content was mapped as a Mendelian trait to LG11.

Extensive macrosynteny between Medicago truncatula and Lens culinaris ssp. culinaris

The first predominantly gene-based genetic linkage map of lentil (Lens culinaris ssp. culinaris) was constructed using an F5 population developed from a cross between the cultivars Digger (ILL5722) and Northfield (ILL5588) using 79 intron-targeted amplified polymorphic (ITAP) and 18 genomic simple sequence repeat (SSR) markers. Linkage analysis revealed seven linkage groups (LGs) comprised of 5-25 markers that varied in length from 80.2 to 274.6 cM. The genome map spanned a total length of 928.4 cM. Clear evidence of a simple and direct macrosyntenic relationship between lentil and Medicago truncatula was observed. Sixty-six out of the 71 gene-based markers, which were previously assigned to M. truncatula genetic and physical maps, were found in regions syntenic between the Lens c. ssp. culinaris and M. truncatula genomes. However, there was evidence of moderate chromosomal rearrangements which may account for the difference in chromosome numbers between these two legume species. Eighteen common SSR markers were used to connect the current map with the most comprehensive and recent map that exists for lentil, providing the syntenic context of four important domestication traits. The composite map presented, anchored with orthologous markers mapped in M. truncatula, provides a strong foundation for the future use of genomic and genetic information in lentil genetic analysis and breeding.

The Medicago truncatula reference accession A17 has an aberrant chromosomal configuration

Medicago truncatula (barrel medic) has emerged as a model legume and accession A17 is the reference genotype selected for the sequencing of the genome. In the present study we compare the A17 chromosomal configuration with that of other accessions by examining pollen viability and genetic maps of intraspecific hybrids. Hybrids derived from crosses between M. truncatula accessions, representative of the large genetic variation within the germplasm collection, were evaluated for pollen viability using Alexander's stain. Genetic maps were generated for the following crosses: SA27063 x SA3054 (n = 94), SA27063 x A17 (n = 92), A17 x Borung (n = 99) and A17 x A20 (n = 69). All F(1) individuals derived from crosses involving A17 showed 50% pollen viability or less. Examination of the recombination frequencies between markers of chromosomes 4 and 8 revealed an apparent genetic linkage between the lower arms of these chromosomes in genetic maps derived from A17. Semisterility and unexpected linkage relationship are both good indicators of a reciprocal translocation. The implications of the A17 distinctive chromosomal rearrangement on studies of genetic mapping, genome sequencing and synteny between species are discussed.

Assessment of molecular diversity in landraces of bread wheat (Triticum aestivum L.) held in an ex situ collection with Diversity Arrays Technology (DArT™)

Diversity Arrays Technology (DArT™) was evaluated as a tool for determining molecular diversity of wheat landraces held within the Australian Winter Cereals Collection (AWCC). Initially, a set of 44 wheat landraces was evaluated with 256 DArT markers. The dataset was compared with the results obtained using 16 amplified fragment length polymorphism (AFLP) primer combinations and 63 simple sequence repeat (SSR) markers, mapped on the 21 chromosomes, from a previous study. The DArT markers exhibited a strong positive correlation with AFLP and SSR, with each marker type distinguishing similar relationships among the 44 landrace accessions. The DArT markers exhibited a higher polymorphic information content than AFLP, and were comparable with that obtained with SSR. Three hundred and fifty-five DArT markers were then used to evaluate genetic diversity among 705 wheat landrace accessions from within the AWCC, chosen to represent 5 world regions. DArT analysis was capable of distinguishing accessions from different geographic regions, and suggested that accessions originating from Nepal represent a unique gene pool within the collection. A statistical resampling of DArT loci indicated that 10–20 loci were enough to distinguish the maximum molecular diversity present within the collection. This research demonstrates the efficacy of the DArT platform as a tool for efficient examination of wheat diversity. As an ex situ germplasm repository, the AWCC contains wheat accessions of high genetic diversity, from genetically differentiated collection sites, even though diversity was under-represented in some countries represented in the repository.

Identification of Sources of Resistance to Phoma medicaginis Isolates in Medicago truncatula SARDI Core Collection Accessions, and Multigene Differentiation of Isolates

ABSTRACT Phoma medicaginis is a necrotrophic fungal pathogen, commonly found infecting the annual medic Medicago truncatula. To differentiate eight P. medicaginis isolates, five gene regions were examined: actin, beta-tubulin, calmodulin, translation elongation factor 1-alpha (EF-1alpha), and the internal transcribed spacer ribosomal DNA. Sequence comparisons showed that specimens isolated from M. truncatula in Western Australia formed a group that was consistently different from, but allied to, a P. medicaginis var. medicaginis type specimen. EF-1alpha contained a hyper-variable 55-bp repeat unit, which forms the basis of a rapid polymerase chain reaction-based method of reliably distinguishing isolates. Characterization of three isolates showed that all exhibited a narrow host range, causing disease only in M. sativa and M. truncatula among eight commonly cultivated legume species sampled. Infection of 86 M. truncatula single-seeded accessions showed a continuous distribution in disease phenotypes, with the majority of accessions susceptible. On a 1-to-5 disease reaction scale increasing in severity, individual fungal isolates showed means of 2.6 to 3.2, and scores ranged from 1 to 4.8 among accessions. The results presented here suggest that M. truncatula harbors specific and diverse sources of resistance to individual P. medicaginis genotypes.