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

Genome assembly of Medicago truncatula accession SA27063 provides insight into spring black stem and leaf spot disease resistance

Medicago truncatula, model legume and alfalfa relative, has served as an essential resource for advancing our understanding of legume physiology, functional genetics, and crop improvement traits. Necrotrophic fungus, Ascochyta medicaginicola, the causal agent of spring black stem (SBS) and leaf spot is a devasting foliar disease of alfalfa affecting stand survival, yield, and forage quality. Host resistance to SBS disease is poorly understood, and control methods rely on cultural practices. Resistance has been observed in M. truncatula accession SA27063 (HM078) with two recessively inherited quantitative-trait loci (QTL), rnpm1 and rnpm2, previously reported. To shed light on host resistance, we carried out a de novo genome assembly of HM078. The genome, referred to as MtHM078 v1.0, is comprised of 23 contigs totaling 481.19 Mbp. Notably, this assembly contains a substantial amount of novel centromere-related repeat sequences due to deep long-read sequencing. Genome annotation resulted in 98.4% of BUSCO fabales proteins being complete. The assembly enabled sequence-level analysis of rnpm1 and rnpm2 for gene content, synteny, and structural variation between SBS-resistant accession SA27063 (HM078) and SBS-susceptible accession A17 (HM101). Fourteen candidate genes were identified, and some have been implicated in resistance to necrotrophic fungi. Especially interesting candidates include loss-of-function events in HM078 because they fit the inverse gene-for-gene model, where resistance is recessively inherited. In rnpm1, these include a loss-of-function in a disease resistance gene due to a premature stop codon, and a 10.85 kbp retrotransposon-like insertion disrupting a ubiquitin conjugating E2. In rnpm2, we identified a frameshift mutation causing a loss-of-function in a glycosidase, as well as a missense and frameshift mutation altering an F-box family protein. This study generated a high-quality genome of HM078 and has identified promising candidates, that once validated, could be further studied in alfalfa to enhance disease resistance.

A pan-genome and chromosome-length reference genome of narrow-leafed lupin (Lupinus angustifolius) reveals genomic diversity and insights into key industry and biological traits

Narrow-leafed lupin (NLL; Lupinus angustifolius) is a key rotational crop for sustainable farming systems, whose grain is high in protein content. It is a gluten-free, non-genetically modified, alternative protein source to soybean (Glycine max) and as such has gained interest as a human food ingredient. Here, we present a chromosome-length reference genome for the species and a pan-genome assembly comprising 55 NLL lines, including Australian and European cultivars, breeding lines and wild accessions. We present the core and variable genes for the species and report on the absence of essential mycorrhizal associated genes. The genome and pan-genomes of NLL and its close relative white lupin (Lupinus albus) are compared. Furthermore, we provide additional evidence supporting LaRAP2-7 as the key alkaloid regulatory gene for NLL and demonstrate the NLL genome is underrepresented in classical NLR disease resistance genes compared to other sequenced legume species. The NLL genomic resources generated here coupled with previously generated RNA sequencing datasets provide new opportunities to fast-track lupin crop improvement.

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.

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.

Genotyping-by-sequencing markers facilitate the identification of quantitative trait loci controlling resistance to Penicillium expansum in Malus sieversii

Blue mold caused by Penicillium expansum is the most important postharvest disease of apple worldwide and results in significant financial losses. There are no defined sources of resistance to blue mold in domesticated apple. However, resistance has been described in wild Malus sieversii accessions, including plant introduction (PI)613981. The objective of the present study was to identify the genetic loci controlling resistance to blue mold in this accession. We describe the first quantitative trait loci (QTL) reported in the Rosaceae tribe Maleae conditioning resistance to P. expansum on genetic linkage group 3 (qM-Pe3.1) and linkage group 10 (qM-Pe10.1). These loci were identified in a M.× domestica 'Royal Gala' X M. sieversii PI613981 family (GMAL4593) based on blue mold lesion diameter seven days post-inoculation in mature, wounded apple fruit inoculated with P. expansum. Phenotypic analyses were conducted in 169 progeny over a four year period. PI613981 was the source of the resistance allele for qM-Pe3.1, a QTL with a major effect on blue mold resistance, accounting for 27.5% of the experimental variability. The QTL mapped from 67.3 to 74 cM on linkage group 3 of the GMAL4593 genetic linkage map. qM-Pe10.1 mapped from 73.6 to 81.8 cM on linkage group 10. It had less of an effect on resistance, accounting for 14% of the experimental variation. 'Royal Gala' was the primary contributor to the resistance effect of this QTL. However, resistance-associated alleles in both parents appeared to contribute to the least square mean blue mold lesion diameter in an additive manner at qM-Pe10.1. A GMAL4593 genetic linkage map composed of simple sequence repeats and 'Golden Delicious' single nucleotide polymorphism markers was able to detect qM-Pe10.1, but failed to detect qM-Pe3.1. The subsequent addition of genotyping-by-sequencing markers to the linkage map provided better coverage of the PI613981 genome on linkage group 3 and facilitated discovery of qM-Pe3.1. A DNA test for qM-Pe3.1 has been developed and is currently being evaluated for its ability to predict blue mold resistance in progeny segregating for qM-Pe3.1. Due to the long juvenility of apple, the availability of a DNA test to screen for the presence of qM-Pe3.1 at the seedling stage will greatly improve efficiency of breeding apple for blue mold resistance.

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.

A rapid method for profiling of volatile and semi-volatile phytohormones using methyl chloroformate derivatisation and GC-MS

Phytohormones are central components of complex signalling networks in plants. The interplay between these metabolites, which include abscisic acid (ABA), auxin (IAA), ethylene, jasmonic acid (JA) and salicylic acid (SA), regulate plant growth and development and modulate responses to biotic and abiotic stress. Few methods of phytohormone profiling can adequately quantify a large range of plant hormones simultaneously and without the requirement for laborious or highly specialised extraction protocols. Here we describe the development and validation of a phytohormone profiling protocol, based on methyl-chloroformate derivatisation of the plant metabolites and analysis by gas chromatography/mass spectrometry (GC-MS). We describe the analysis of 11 metabolites, either plant phytohormones or intermediates of phytohormone metabolism; ABA, azelaic acid, IAA, JA and SA, and the phytohormone precursors 1-aminocyclopropane 1-carboxylic acid, benzoic acid, cinnamic acid, 13-epi-12-oxophytodienoic acid (13-epi-OPDA), linoleic acid and linolenic acid, and validate the isolation from foliar tissue of the model legume Medicago truncatula. The preparation is insensitive to the presence of water, facilitating measurement of the volatile metabolites. Quantitation was linear over four orders of magnitude, and the limits of detection between two and 10 ng/mL for all measured metabolites using a single quadrupole GC-MS.

Transcriptome sequencing of different narrow-leafed lupin tissue types provides a comprehensive uni-gene assembly and extensive gene-based molecular markers

Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is an important grain legume crop that is valuable for sustainable farming and is becoming recognized as a human health food. NLL breeding is directed at improving grain production, disease resistance, drought tolerance and health benefits. However, genetic and genomic studies have been hindered by a lack of extensive genomic resources for the species. Here, the generation, de novo assembly and annotation of transcriptome datasets derived from five different NLL tissue types of the reference accession cv. Tanjil are described. The Tanjil transcriptome was compared to transcriptomes of an early domesticated cv. Unicrop, a wild accession P27255, as well as accession 83A:476, together being the founding parents of two recombinant inbred line (RIL) populations. In silico predictions for transcriptome-derived gene-based length and SNP polymorphic markers were conducted and corroborated using a survey assembly sequence for NLL cv. Tanjil. This yielded extensive indel and SNP polymorphic markers for the two RIL populations. A total of 335 transcriptome-derived markers and 66 BAC-end sequence-derived markers were evaluated, and 275 polymorphic markers were selected to genotype the reference NLL 83A:476 × P27255 RIL population. This significantly improved the completeness, marker density and quality of the reference NLL genetic map.

Necrosis-inducing proteins of Rhynchosporium commune, effectors in quantitative disease resistance

The barley pathogen Rhynchosporium commune secretes necrosis-inducing proteins NIP1, NIP2, and NIP3. Expression analysis revealed that NIP1 transcripts appear to be present in fungal spores already, whereas NIP2 and NIP3 are synthesized after inoculation of host plants. To assess the contribution of the three effector proteins to disease development, deletion mutants were generated. The development of these fungal mutants on four barley cultivars was quantified in comparison with that of the parent wild-type strain and with two fungal strains failing to secrete an "active" NIP1 avirulence protein, using quantitative polymerase chain reaction as well as microscopic imaging after fungal green fluorescent protein tagging. The impact of the three deletions varied quantitatively depending on the host genotype, suggesting that the activities of the fungal effectors add up to produce stronger growth patterns and symptom development. Alternatively, recognition events of differing intensities may be converted into defense gene expression in a quantitative manner.

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.

Identification of distinct quantitative trait loci associated with defence against the closely related aphids Acyrthosiphon pisum and A. kondoi in Medicago truncatula

Aphids are a major family of plant insect pests. Medicago truncatula and Acyrthosiphon pisum (pea aphid, PA) are model species with a suite of resources available to help dissect the mechanism underlying plant-aphid interactions. A previous study focused on monogenic and relatively strong resistance in M. truncatula to PA and other aphid species. In this study a moderate resistance to PA was characterized in detail in the M. truncatula line A17 and compared with the highly susceptible line A20 and the more resistant line Jester. The results show that PA resistance in A17 involves both antibiosis and tolerance, and that resistance is phloem based. Quantitative trait locus (QTL) analysis using a recombinant inbred line (RIL) population (n=114) from a cross between A17 and A20 revealed that one locus, which co-segregated with AIN (Acyrthosiphon-induced necrosis) on chromosome 3, is responsible for the reduction of aphid biomass (indicator of antibiosis) for both PA and bluegreen aphid (BGA, A. kondoi), albeit to a lesser degree for PA than BGA. Interestingly, two independent loci on chromosomes 5 and 3 were identified for the plant biomass reduction (indicator of plant tolerance) by PA and BGA, respectively, demonstrating that the plant's tolerance response to these two closely related aphid species is distinct. Together with previously identified major resistant (R) genes, the QTLs identified in this study are powerful tools to understand fully the spectrum of plant defence against sap-sucking insects and provide opportunities for breeders to generate effective and sustainable strategies for aphid control.

Using a physiological framework for improving the detection of quantitative trait loci related to nitrogen nutrition in Medicago truncatula

Medicago truncatula is used as a model plant for exploring the genetic and molecular determinants of nitrogen (N) nutrition in legumes. In this study, our aim was to detect quantitative trait loci (QTL) controlling plant N nutrition using a simple framework of carbon/N plant functioning stemming from crop physiology. This framework was based on efficiency variables which delineated the plant's efficiency to take up and process carbon and N resources. A recombinant inbred line population (LR4) was grown in a glasshouse experiment under two contrasting nitrate concentrations. At low nitrate, symbiotic N(2) fixation was the main N source for plant growth and a QTL with a large effect located on linkage group (LG) 8 affected all the traits. Significantly, efficiency variables were necessary both to precisely localize a second QTL on LG5 and to detect a third QTL involved in epistatic interactions on LG2. At high nitrate, nitrate assimilation was the main N source and a larger number of QTL with weaker effects were identified compared to low nitrate. Only two QTL were common to both nitrate treatments: a QTL of belowground biomass located at the bottom of LG3 and another one on LG6 related to three different variables (leaf area, specific N uptake and aboveground:belowground biomass ratio). Possible functions of several candidate genes underlying QTL of efficiency variables could be proposed. Altogether, our results provided new insights into the genetic control of N nutrition in M. truncatula. For instance, a novel result for M. truncatula was identification of two epistatic interactions in controlling plant N(2) fixation. As such this study showed the value of a simple conceptual framework based on efficiency variables for studying genetic determinants of complex traits and particularly epistatic interactions.

A complex genetic network involving a broad-spectrum locus and strain-specific loci controls resistance to different pathotypes of Aphanomyces euteiches in Medicago truncatula

A higher understanding of genetic and genomic bases of partial resistance in plants and their diversity regarding pathogen variability is required for a more durable management of resistance genetic factors in sustainable cropping systems. In this study, we investigated the diversity of genetic factors involved in partial resistance to Aphanomyces euteiches, a very damaging pathogen on pea and alfalfa, in Medicago truncatula. A mapping population of 178 recombinant inbred lines, from the cross F83005.5 (susceptible) and DZA045.5 (resistant), was used to identify quantitative trait loci for resistance to four A. euteiches reference strains belonging to the four main pathotypes currently known on pea and alfalfa. A major broad-spectrum genomic region, previously named AER1, was localized to a reduced 440 kb interval on chromosome 3 and was involved in complete or partial resistance, depending on the A. euteiches strain. We also identified 21 additive and/or epistatic genomic regions specific to one or two strains, several of them being anchored to the M. truncatula physical map. These results show that, in M. truncatula, a complex network of genetic loci controls partial resistance to different pea and alfalfa pathotypes of A. euteiches, suggesting a diversity of molecular mechanisms underlying partial resistance.

Partial resistance of Medicago truncatula to Aphanomyces euteiches is associated with protection of the root stele and is controlled by a major QTL rich in proteasome-related genes

A pathosystem between Aphanomyces euteiches, the causal agent of pea root rot disease, and the model legume Medicago truncatula was developed to gain insights into mechanisms involved in resistance to this oomycete. The F83005.5 French accession and the A17-Jemalong reference line, susceptible and partially resistant, respectively, to A. euteiches, were selected for further cytological and genetic analyses. Microscopy analyses of thin root sections revealed that a major difference between the two inoculated lines occurred in the root stele, which remained pathogen free in A17. Striking features were observed in A17 roots only, including i) frequent pericycle cell divisions, ii) lignin deposition around the pericycle, and iii) accumulation of soluble phenolic compounds. Genetic analysis of resistance was performed on an F7 population of 139 recombinant inbred lines and identified a major quantitative trait locus (QTL) near the top of chromosome 3. A second study, with near-isogenic line responses to A. euteiches confirmed the role of this QTL in expression of resistance. Fine-mapping allowed the identification of a 135-kb sequenced genomic DNA region rich in proteasome-related genes. Most of these genes were shown to be induced only in inoculated A17. Novel mechanisms possibly involved in the observed partial resistance are proposed.

Translating Medicago truncatula genomics to crop legumes

Genomic resources developed in the model legume, Medicago truncatula, have the potential to accelerate practical advances in crop legumes. M. truncatula is closely related to many economically important legumes, frequently displaying genome-scale synteny. Translating genome data from M. truncatula should be highly effective in marker development, gene discovery, and positional cloning in crop legumes. The M. truncatula genome sequence also provides valuable insights about gene families of practical importance, especially those that are legume-specific. The M. truncatula genome sequence should also simplify the assembly of next-generation sequence data in closely related taxa, especially alfalfa. Genomic resources, such as whole-genome arrays, make it possible to pursue detailed questions about gene expression in both M. truncatula and related crop species, while tagged mutant populations simplify the process of determining gene function.

AER1, a major gene conferring resistance to Aphanomyces euteiches in Medicago truncatula

Aphanomyces euteiches is a major soilborne oomycete pathogen that infects various legume species, including pea and alfalfa. The model legume Medicago truncatula has recently emerged as a valuable genetic system for understanding the genetic basis of resistance to A. euteiches in leguminous crops. The objective of this study was to identify genetic determinants of resistance to a broad host-range pea-infecting strain of A. euteiches in M. truncatula. Two M. truncatula segregating populations of 178 F(5) recombinant inbred lines and 200 F(3) families from the cross F83005.5 (susceptible) x DZA045.5 (resistant) were screened for resistance to A. euteiches. Phenotypic distributions observed suggested a dominant monogenic control of resistance. A major locus associated with resistance to A. euteiches, namely AER1, was mapped by bulk segregant analysis to a terminal end of chromosome 3 in M. truncatula and explained 88% of the phenotypic variation. AER1 was identified in a resistance-gene-rich region, where resistance gene analogs and genes associated with disease resistance phenotypes have been identified. Discovery of AER1 opens up new prospects for improving resistance to A. euteiches in cultivated legumes using a comparative genomics approach.

Can differences of nitrogen nutrition level among Medicago truncatula genotypes be assessed non-destructively?: Probing with a recombinant inbred lines population

The international consensus on Medicago truncatula as a model system has lead to the development of powerful approaches for dissecting the genetic and molecular bases of legume nitrogen nutrition. However, such approaches now come up against a poor knowledge of the phenotypic traits that should be used for the large-scale screening of the genotypic variability associated with nitrogen nutrition. This issue was unravelled in a previous report, in which an ecophysiological approach allowed a better understanding of the relationships between plant nitrogen nutrition and plant growth traits, for the model symbiotic association between M. truncatula cv. Jemalong and Rhizobium meliloti strain 2011. From this analysis, phenotypic traits were identified as potentially relevant for the large-scale screening of the genotypic variability. Here, by the phenotyping of a recombinant inbred lines population, we show that the proposed methodology provides a valuable support for assisting the detection of genetic variants affected for nitrogen uptake. Especially, the relative expansion rate of plant leaf area is identified as a good proxy for ranking genotypes according to their ability to uptake nitrogen in given environmental conditions. As leaf area can be measured non-destructively, such finding should pave the way for a more efficient evaluation of the genotypic variability.