Genetic dissection of domestication traits in interspecific chickpea populations

Chickpea (Cicer arietinum) is a pulse crop that provides an integral source of nutrition for human consumption. The close wild relatives Cicer reticulatum and Cicer echinospermum harbor untapped genetic diversity that can be exploited by chickpea breeders to improve domestic varieties. Knowledge of genomic loci that control important chickpea domestication traits will expedite the development of improved chickpea varieties derived from interspecific crosses. Therefore, we set out to identify genomic loci underlying key chickpea domestication traits by both association and quantitative trait locus (QTL) mapping using interspecific F2 populations. Diverse phenotypes were recorded for various agronomic traits. A total of 11 high-confidence markers were detected on chromosomes 1, 3, and 7 by both association and QTL mapping; these were associated with growth habit, flowering time, and seed traits. Furthermore, we identified candidate genes linked to these markers, which advanced our understanding of the genetic basis of domestication traits and validated known genes such as the FLOWERING LOCUS gene cluster that regulates flowering time. Collectively, this study has elucidated the genetic basis of chickpea domestication traits, which can facilitate the development of superior chickpea varieties.

Genome-wide identification of Sclerotinia sclerotiorum small RNAs and their endogenous targets

BACKGROUND

Several phytopathogens produce small non-coding RNAs of approximately 18-30 nucleotides (nt) which post-transcriptionally regulate gene expression. Commonly called small RNAs (sRNAs), these small molecules were also reported to be present in the necrotrophic pathogen Sclerotinia sclerotiorum. S. sclerotiorum causes diseases in more than 400 plant species, including the important oilseed crop Brassica napus. sRNAs can further be classified as microRNAs (miRNAs) and short interfering RNAs (siRNAs). Certain miRNAs can activate loci that produce further sRNAs; these secondary sRNA-producing loci are called 'phased siRNA' (PHAS) loci and have only been described in plants. To date, very few studies have characterized sRNAs and their endogenous targets in S. sclerotiorum.

RESULTS

We used Illumina sequencing to characterize sRNAs from fungal mycelial mats of S. sclerotiorum spread over B. napus leaves. In total, eight sRNA libraries were prepared from in vitro, 12 h post-inoculation (HPI), and 24 HPI mycelial mat samples. Cluster analysis identified 354 abundant sRNA clusters with reads of more than 100 Reads Per Million (RPM). Differential expression analysis revealed upregulation of 34 and 57 loci at 12 and 24 HPI, respectively, in comparison to in vitro samples. Among these, 25 loci were commonly upregulated. Altogether, 343 endogenous targets were identified from the major RNAs of 25 loci. Almost 88% of these targets were annotated as repeat element genes, while the remaining targets were non-repeat element genes. Fungal degradome reads confirmed cleavage of two transposable elements by one upregulated sRNA. Altogether, 24 milRNA loci were predicted with both mature and milRNA* (star) sequences; these are both criteria associated previously with experimentally verified miRNAs. Degradome sequencing data confirmed the cleavage of 14 targets. These targets were related to repeat element genes, phosphate acetyltransferases, RNA-binding factor, and exchange factor. A PHAS gene prediction tool identified 26 possible phased interfering loci with 147 phasiRNAs from the S. sclerotiorum genome, suggesting this pathogen might produce sRNAs that function similarly to miRNAs in higher eukaryotes.

CONCLUSIONS

Our results provide new insights into sRNA populations and add a new resource for the study of sRNAs in S. sclerotiorum.

Virulence Profiles and Genome-Wide Association Study for Ascochyta lentis Isolates Collected from Australian Lentil-Growing Regions

Ascochyta lentis, the causal organism of Ascochyta blight (AB) of lentil (Lens culinaris), has been shown to produce an avirulence effector protein that mediates AB resistance in certain lentil cultivars. The two known forms of the effector protein were identified from a biparental mapping population between isolates that have reciprocal virulence on 'PBA Hurricane XT' and 'Nipper'. The effector AlAvr1-1 was described for the PBA Hurricane XT-avirulent isolate P94-24 and AlAvr1-2 characterized in the PBA Hurricane XT-virulent isolate AlKewell. Here, we performed a genome-wide association study to identify other loci associated with AB for a differential set of lentil cultivars from a diverse panel of isolates collected in the Australian lentil-growing regions from 2013 to 2020. The chromosome 3 AlAvr1 locus was strongly associated with the PBA Hurricane XT, 'Indianhead', and Nipper disease responses, but one other genomic region on chromosome 11 was also associated with the Nipper disease trait. Our results corroborate earlier work that identified the AlAvr1 locus for field-collected isolates that span the period before release and after widespread adoption of PBA Hurricane XT. A multiplex PCR assay was developed to differentiate the genes AlAvr1-1 and AlAvr1-2 to predict PBA Hurricane XT avirulence and pathotype designation in the diversity panel. Increasing numbers of the PBA Hurricane XT-virulent pathotype 2 isolates across that time indicate strong selection for isolates with the AlAvr1-2 allele. Furthermore, one other region of the A. lentis genome may contribute to the pathogen-host interaction for lentil AB.

The broad host range pathogen Sclerotinia sclerotiorum produces multiple effector proteins that induce host cell death intracellularly

Sclerotinia sclerotiorum is a broad host range necrotrophic fungal pathogen, which causes disease on many economically important crop species. S. sclerotiorum has been shown to secrete small effector proteins to kill host cells and acquire nutrients. We set out to discover novel necrosis-inducing effectors and characterize their activity using transient expression in Nicotiana benthamiana leaves. Five intracellular necrosis-inducing effectors were identified with differing host subcellular localization patterns, which were named intracellular necrosis-inducing effector 1-5 (SsINE1-5). We show for the first time a broad host range pathogen effector, SsINE1, that uses an RxLR-like motif to enter host cells. Furthermore, we provide preliminary evidence that SsINE5 induces necrosis via an NLR protein. All five of the identified effectors are highly conserved in globally sourced S. sclerotiorum isolates. Taken together, these results advance our understanding of the virulence mechanisms employed by S. sclerotiorum and reveal potential avenues for enhancing genetic resistance to this damaging fungal pathogen.

Field Pea (Pisum sativum) Germplasm Screening for Seedling Ascochyta Blight Resistance and Genome-Wide Association Studies Reveal Loci Associated with Resistance to Peyronellaea pinodes and Ascochyta koolunga

Ascochyta blight is a damaging disease that affects the stems, leaves, and pods of field pea (Pisum sativum) and impacts yield and grain quality. In Australia, field pea Ascochyta blight is primarily caused by the necrotrophic fungal species Peyronellaea pinodes and Ascochyta koolunga. In this study, we screened 1,276 Pisum spp. germplasm accessions in seedling disease assays with a mix of three isolates of P. pinodes and 641 accessions with three mixed isolates of A. koolunga (513 accessions were screened with both species). A selection of three P. sativum accessions with low disease scores for either pathogen, or in some cases both, were crossed with Australian field pea varieties PBA Gunyah and PBA Oura, and recombinant inbred line populations were made. Populations at the F3:4 and F4:5 generation were phenotyped for their disease response to P. pinodes and A. koolunga, and genotypes were determined using the diversity arrays technology genotyping method. Marker-trait associations were identified using a genome-wide association study approach. Trait-associated loci were mapped to the published P. sativum genome assembly, and candidate resistance gene analogues were identified in the corresponding genomic regions. One locus on chromosome 2 (LG1) was associated with resistance to P. pinodes, and the 8 Mb genomic region contains 156 genes, two of which are serine/threonine protein kinases, putatively contributing to the resistance trait. A second locus on chromosome 5 (LG3) was associated with resistance to A. koolunga, and the 35 Mb region contains 488 genes, of which five are potential candidate resistance genes, including protein kinases, a mitogen-activated protein kinase, and an ethylene-responsive protein kinase homolog.

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.

Identification of Sclerotinia stem rot resistance quantitative trait loci in a chickpea (Cicer arietinum) recombinant inbred line population

Sclerotinia stem rot (SSR), caused by Sclerotinia sclerotiorum , is one of the most economically devastating diseases in chickpea (Cicer arietinum L.). No complete resistance is available in chickpea to this disease, and the inheritance of partial resistance is not understood. Two hundred F7 recombinant inbred lines (RILs) derived from a cross between a partially resistant variety PBA HatTrick, and a highly susceptible variety Kyabra were characterised for their responses to SSR inoculation. Quantitative trait locus (QTL) analysis was conducted for the area under the disease progress curve (AUDPC) after RIL infection with S. sclerotiorum . Four QTLs on chromosomes, Ca4 (qSSR4-1, qSSR4-2), Ca6 (qSSR6-1) and Ca7 (qSSR7-1), individually accounted for between 4.2 and 15.8% of the total estimated phenotypic variation for the response to SSR inoculation. Candidate genes located in these QTL regions are predicted to be involved in a wide range of processes, including phenylpropanoid biosynthesis, plant-pathogen interaction, and plant hormone signal transduction. This is the first study investigating the inheritance of resistance to S. sclerotiorum in chickpea. Markers associated with the identified QTLs could be employed for marker-assisted selection in chickpea breeding.

A Trimethylguanosine Synthase1-like (TGS1) homologue is implicated in vernalisation and flowering time control

KEY MESSAGE

A plant-specific Trimethylguanosine Synthase1-like homologue was identified as a candidate gene for the efl mutation in narrow-leafed lupin, which alters phenology by reducing vernalisation requirement. The vernalisation pathway is a key component of flowering time control in plants from temperate regions but is not well understood in the legume family. Here we examined vernalisation control in the temperate grain legume species, narrow-leafed lupin (Lupinus angustifolius L.), and discovered a candidate gene for an ethylene imine mutation (efl). The efl mutation changes phenology from late to mid-season flowering and additionally causes transformation from obligate to facultative vernalisation requirement. The efl locus was mapped to pseudochromosome NLL-10 in a recombinant inbred line (RIL) mapping population developed by accelerated single seed descent. Candidate genes were identified in the reference genome, and a diverse panel of narrow-leafed lupins was screened to validate mutations specific to accessions with efl. A non-synonymous SNP mutation within an S-adenosyl-L-methionine-dependent methyltransferase protein domain of a Trimethylguanosine Synthase1-like (TGS1) orthologue was identified as the candidate mutation giving rise to efl. This mutation caused substitution of an amino acid within an established motif at a position that is otherwise highly conserved in several plant families and was perfectly correlated with the efl phenotype in F2 and F6 genetic population and a panel of diverse accessions, including the original efl mutant. Expression of the TGS1 homologue did not differ between wild-type and efl genotypes, supporting altered functional activity of the gene product. This is the first time a TGS1 orthologue has been associated with vernalisation response and flowering time control in any plant species.

Identification of Sources of Sclerotinia sclerotiorum Resistance in a Collection of Wild Cicer Germplasm

Sclerotinia sclerotiorum is an important fungal pathogen of chickpea (Cicer arietinum L.), and it can cause yield losses up to 100%. The wild progenitors are much more diverse than domesticated chickpea, and this study describes how this relates to S. sclerotiorum resistance. Initially, the pathogenicity of nine Australian S. sclerotiorum isolates was examined on three Cicer lines to develop a robust phenotyping assay, and significant differences in isolate aggressiveness were identified with six isolates being classed as highly aggressive and three as moderately aggressive. We identified two S. sclerotiorum isolates, CU8.20 and CU10.12, to be highly aggressive and moderately aggressive, respectively. A subsequent phenotyping assay was conducted using the two isolates to evaluate 86 wild Cicer accessions (Cicer reticulatum and Cicer echinospermum) and two C. arietinum varieties for resistance to S. sclerotiorum. A subset of 12 genotypes was further evaluated, and subsequently, two wild Cicer accessions with consistently high levels of resistance to S. sclerotiorum were examined using the initially characterized nine isolates. Wild Cicer accessions Karab_084 and Deste_063 demonstrated consistent partial resistance to S. sclerotiorum. There were significant differences in responses to S. sclerotiorum across wild Cicer collection sites. The Cermik, Karabahce, and Destek sites' responses to the aggressive isolate CU8.20 ranged from resistant to susceptible, highlighting an interaction between isolate genotype and chickpea collection site for sclerotinia stem rot resistance. This is the first evidence of partial stem resistance identified in wild Cicer germplasm, which can be adopted in chickpea breeding programs to enhance S. sclerotiorum resistance in future chickpea varieties.

fIdentification of B. napus small RNAs responsive to infection by a necrotrophic pathogen

BACKGROUND

Small RNAs are short non-coding RNAs that are key gene regulators controlling various biological processes in eukaryotes. Plants may regulate discrete sets of sRNAs in response to pathogen attack. Sclerotinia sclerotiorum is an economically important pathogen affecting hundreds of plant species, including the economically important oilseed B. napus. However, there are limited studies on how regulation of sRNAs occurs in the S. sclerotiorum and B. napus pathosystem.

RESULTS

We identified different classes of sRNAs from B. napus using high throughput sequencing of replicated mock and infected samples at 24 h post-inoculation (HPI). Overall, 3999 sRNA loci were highly expressed, of which 730 were significantly upregulated during infection. These 730 up-regulated sRNAs targeted 64 genes, including disease resistance proteins and transcriptional regulators. A total of 73 conserved miRNA families were identified in our dataset. Degradome sequencing identified 2124 cleaved mRNA products from these miRNAs from combined mock and infected samples. Among these, 50 genes were specific to infection. Altogether, 20 conserved miRNAs were differentially expressed and 8 transcripts were cleaved by the differentially expressed miRNAs miR159, miR5139, and miR390, suggesting they may have a role in the S. sclerotiorum response. A miR1885-triggered disease resistance gene-derived secondary sRNA locus was also identified and verified with degradome sequencing. We also found further evidence for silencing of a plant immunity related ethylene response factor gene by a novel sRNA using 5'-RACE and RT-qPCR.

CONCLUSIONS

The findings in this study expand the framework for understanding the molecular mechanisms of the S. sclerotiorum and B. napus pathosystem at the sRNA level.

Modeling first order additive × additive epistasis improves accuracy of genomic prediction for sclerotinia stem rot resistance in canola

The fungus Sclerotinia sclerotiorum infects hundreds of plant species including many crops. Resistance to this pathogen in canola (Brassica napus L. subsp. napus) is controlled by numerous quantitative trait loci (QTL). For such polygenic traits, genomic prediction may be useful for breeding as it can capture many QTL at once while also considering nonadditive genetic effects. Here, we test application of common regression models to genomic prediction of S. sclerotiorum resistance in canola in a diverse panel of 218 plants genotyped at 24,634 loci. Disease resistance was scored by infection with an aggressive isolate and monitoring over 3 wk. We found that including first-order additive × additive epistasis in linear mixed models (LMMs) improved accuracy of breeding value estimation between 3 and 40%, depending on method of assessment, and correlation between phenotypes and predicted total genetic values by 14%. Bayesian models performed similarly to or worse than genomic relationship matrix-based models for estimating breeding values or overall phenotypes from genetic values. Bayesian ridge regression, which is most similar to the genomic relationship matrix-based approach in the amount of shrinkage it applies to marker effects, was the most accurate of this family of models. This confirms several studies indicating the highly polygenic nature of sclerotinia stem rot resistance. Overall, our results highlight the use of simple epistasis terms for prediction of breeding values and total genetic values for a complex disease resistance phenotype in canola.

Analysis of differentially expressed Sclerotinia sclerotiorum genes during the interaction with moderately resistant and highly susceptible chickpea lines

BACKGROUND

Sclerotinia sclerotiorum, the cause of Sclerotinia stem rot (SSR), is a host generalist necrotrophic fungus that can cause major yield losses in chickpea (Cicer arietinum) production. This study used RNA sequencing to conduct a time course transcriptional analysis of S. sclerotiorum gene expression during chickpea infection. It explores pathogenicity and developmental factors employed by S. sclerotiorum during interaction with chickpea.

RESULTS

During infection of moderately resistant (PBA HatTrick) and highly susceptible chickpea (Kyabra) lines, 9491 and 10,487 S. sclerotiorum genes, respectively, were significantly differentially expressed relative to in vitro. Analysis of the upregulated genes revealed enrichment of Gene Ontology biological processes, such as oxidation-reduction process, metabolic process, carbohydrate metabolic process, response to stimulus, and signal transduction. Several gene functional categories were upregulated in planta, including carbohydrate-active enzymes, secondary metabolite biosynthesis clusters, transcription factors and candidate secreted effectors. Differences in expression of four S. sclerotiorum genes on varieties with different levels of susceptibility were also observed.

CONCLUSION

These findings provide a framework for a better understanding of S. sclerotiorum interactions with hosts of varying susceptibility levels. Here, we report for the first time on the S. sclerotiorum transcriptome during chickpea infection, which could be important for further studies on this pathogen's molecular biology.

Identification of Novel Sources of Resistance to Ascochyta Blight in a Collection of Wild Cicer Accessions

Chickpea production is constrained worldwide by the necrotrophic fungal pathogen Ascochyta rabiei, the causal agent of Ascochyta blight (AB). To reduce the impact of this disease, novel sources of resistance are required in chickpea cultivars. Here, we screened a new collection of wild Cicer accessions for AB resistance and identified accessions resistant to multiple, highly pathogenic isolates. In addition to this, analyses demonstrated that some collection sites of C. echinospermum harbor predominantly resistant accessions, knowledge that can inform future collection missions. Furthermore, a genome-wide association study identified regions of the C. reticulatum genome associated with AB resistance and investigation of these regions identified candidate resistance genes. Taken together, these results can be utilized to enhance the resistance of chickpea cultivars to this globally yield-limiting disease.

A functional genomics approach to dissect spotted alfalfa aphid resistance in Medicago truncatula

Aphids are virus-spreading insect pests affecting crops worldwide and their fast population build-up and insecticide resistance make them problematic to control. Here, we aim to understand the molecular basis of spotted alfalfa aphid (SAA) or Therioaphis trifolii f. maculata resistance in Medicago truncatula, a model organism for legume species. We compared susceptible and resistant near isogenic Medicago lines upon SAA feeding via transcriptome sequencing. Expression of genes involved in defense and stress responses, protein kinase activity and DNA binding were enriched in the resistant line. Potentially underlying some of these changes in gene expression was the finding that members of the MYB, NAC, AP2 domain and ERF transcription factor gene families were differentially expressed in the resistant versus susceptible lines. A TILLING population created in the resistant cultivar was screened using exome capture sequencing and served as a reverse genetics tool to functionally characterise genes involved in the aphid resistance response. This screening revealed three transcription factors (a NAC, AP2 domain and ERF) as important regulators in the defence response, as a premature stop-codon in the resistant background led to a delay in aphid mortality and enhanced plant susceptibility. This combined functional genomics approach will facilitate the future development of pest resistant crops by uncovering candidate target genes that can convey enhanced aphid resistance.

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.

An RNAi supplemented diet as a reverse genetics tool to control bluegreen aphid, a major pest of legumes

Aphids are important agricultural pests causing major yield losses worldwide. Since aphids can rapidly develop resistance to chemical insecticides there is an urgent need to find alternative aphid pest management strategies. Despite the economic importance of bluegreen aphid (Acyrthosiphon kondoi), very few genetic resources are available to expand our current understanding and help find viable control solutions. An artificial diet is a desirable non-invasive tool to enable the functional characterisation of genes in bluegreen aphid and discover candidate target genes for future use in RNA interference (RNAi) mediated crop protection against aphids. To date no artificial diet has been developed for bluegreen aphid, so we set out to develop a suitable diet by testing and optimising existing diets. Here, we describe an artificial diet for rearing bluegreen aphid and also provide a proof of concept for the supplementation of the diet with RNAi molecules targeting the salivary gland transcript C002 and gap gene hunchback, resulting in bluegreen aphid mortality which has not yet been documented in this species. Managing this pest, for example via RNAi delivery through artificial feeding will be a major improvement to test bluegreen aphid candidate target genes for future pest control and gain significant insights into bluegreen aphid gene function.

A detailed in silico analysis of secondary metabolite biosynthesis clusters in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum

BACKGROUND

The broad host range pathogen Sclerotinia sclerotiorum infects over 400 plant species and causes substantial yield losses in crops worldwide. Secondary metabolites are known to play important roles in the virulence of plant pathogens, but little is known about the secondary metabolite repertoire of S. sclerotiorum. In this study, we predicted secondary metabolite biosynthetic gene clusters in the genome of S. sclerotiorum and analysed their expression during infection of Brassica napus using an existing transcriptome data set. We also investigated their sequence diversity among a panel of 25 previously published S. sclerotiorum isolate genomes.

RESULTS

We identified 80 putative secondary metabolite clusters. Over half of the clusters contained at least three transcriptionally coregulated genes. Comparative genomics revealed clusters homologous to clusters in the closely related plant pathogen Botrytis cinerea for production of carotenoids, hydroxamate siderophores, DHN melanin and botcinic acid. We also identified putative phytotoxin clusters that can potentially produce the polyketide sclerin and an epipolythiodioxopiperazine. Secondary metabolite clusters were enriched in subtelomeric genomic regions, and those containing paralogues showed a particularly strong association with repeats. The positional bias we identified was borne out by intraspecific comparisons that revealed putative secondary metabolite genes suffered more presence / absence polymorphisms and exhibited a significantly higher sequence diversity than other genes.

CONCLUSIONS

These data suggest that S. sclerotiorum produces numerous secondary metabolites during plant infection and that their gene clusters undergo enhanced rates of mutation, duplication and recombination in subtelomeric regions. The microevolutionary regimes leading to S. sclerotiorum secondary metabolite diversity have yet to be elucidated. Several potential phytotoxins documented in this study provide the basis for future functional analyses.

Additive and epistatic interactions between AKR and AIN loci conferring bluegreen aphid resistance and hypersensitivity in Medicago truncatula

Aphids, including the bluegreen aphid (BGA; Acyrthosiphon kondoi), are important pests in agriculture. Two BGA resistance genes have been identified in the model legume Medicago truncatula, namely AKR (Acyrthosiphon kondoi resistance) and AIN (Acyrthosiphon induced necrosis). In this study, progeny derived from a cross between a resistant accession named Jester and a highly susceptible accession named A20 were used to study the interaction between the AKR and AIN loci with respect to BGA performance and plant response to BGA infestation. These studies demonstrated that AKR and AIN have additive effects on the BGA resistance phenotype. However, AKR exerts dominant suppression epistasis on AIN-controlled macroscopic necrotic lesions. Nevertheless, both AKR and AIN condition production of H2O2 at the BGA feeding site. Electrical penetration graph analysis demonstrated that AKR prevents phloem sap ingestion, irrespective of the presence of AIN. Similarly, the jasmonic acid defense signaling pathway is recruited by AKR, irrespective of AIN. This research identifies an enhancement of aphid resistance through gene stacking, and insights into the interaction of distinct resistance genes against insect pests.

The role of jasmonate signalling in quinolizidine alkaloid biosynthesis, wounding and aphid predation response in narrow-leafed lupin

Quinolizidine alkaloids (QAs) are toxic secondary metabolites produced in lupin species that protect the plant against insects. They form in vegetative tissues and accumulate to a different extent in the grains: high levels in 'bitter' narrow-leafed lupin (NLL) and low levels in 'sweet' NLL. Grain QA levels vary considerably, and sometimes exceed the industry limit for food and feed purposes. We hypothesised that jasmonates regulate QA biosynthesis in response to environmental stresses such as wounding and aphid predation, which may explain non-genetic variability in grain QA levels. Methyl jasmonate (MeJA)-inducible genes were identified and verified in NLL. Exogenous MeJA application-induced expression of QA biosynthetic genes and QA levels for bitter, but not sweet NLL. Although MeJA-inducible genes responded to wounding, the expression of QA biosynthetic genes was not induced for bitter and sweet NLL. We assessed the effect of aphid predation on QA production for two cultivars - one moderately resistant and one susceptible to aphid predation. Although MeJA-inducible genes responded to aphid predation, no change in QA levels was found for either cultivar. These findings offer insights into the regulation of QA biosynthesis in bitter and sweet NLL and concludes that aphids are not a concern for increasing grain QAs in NLL cultivars.

A whole genome scan of SNP data suggests a lack of abundant hard selective sweeps in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum

The pathogenic fungus Sclerotinia sclerotiorum infects over 600 species of plant. It is present in numerous environments throughout the world and causes significant damage to many agricultural crops. Fragmentation and lack of gene flow between populations may lead to population sub-structure. Within discrete recombining populations, positive selection may lead to a ‘selective sweep’. This is characterised by an increase in frequency of a favourable allele leading to reduction in genotypic diversity in a localised genomic region due to the phenomenon of genetic hitchhiking. We aimed to assess whether isolates of S. sclerotiorum from around the world formed genotypic clusters associated with geographical origin and to determine whether signatures of population-specific positive selection could be detected. To do this, we sequenced the genomes of 25 isolates of S. sclerotiorum collected from four different continents–Australia, Africa (north and south), Europe and North America (Canada and the northen United States) and conducted SNP based analyses of population structure and selective sweeps. Among the 25 isolates, there was evidence for two major population clusters. One of these consisted of 11 isolates from Canada, the USA and France (population 1), and the other consisted of nine isolates from Australia and one from Morocco (population 2). The rest of the isolates were genotypic outliers. We found that there was evidence of outcrossing in these two populations based on linkage disequilibrium decay. However, only a single candidate selective sweep was observed, and it was present in population 2. This sweep was close to a Major Facilitator Superfamily transporter gene, and we speculate that this gene may have a role in nutrient uptake from the host. The low abundance of selective sweeps in the S. sclerotiorum genome contrasts the numerous examples in the genomes of other fungal pathogens. This may be a result of its slow rate of evolution and low effective recombination rate due to self-fertilisation and vegetative reproduction.

Identification and profiling of narrow-leafed lupin (Lupinus angustifolius) microRNAs during seed development

BACKGROUND

Whilst information regarding small RNAs within agricultural crops is increasing, the miRNA composition of the nutritionally valuable pulse narrow-leafed lupin (Lupinus angustifolius) remains unknown.

RESULTS

By conducting a genome- and transcriptome-wide survey we identified 7 Dicer-like and 16 Argonaute narrow-leafed lupin genes, which were highly homologous to their legume counterparts. We identified 43 conserved miRNAs belonging to 16 families, and 13 novel narrow-leafed lupin-specific miRNAs using high-throughput sequencing of small RNAs from foliar and root and five seed development stages. We observed up-regulation of members of the miRNA families miR167, miR399, miR156, miR319 and miR164 in narrow-leafed lupin seeds, and confirmed expression of miR156, miR166, miR164, miR1507 and miR396 using quantitative RT-PCR during five narrow-leafed lupin seed development stages. We identified potential targets for the conserved and novel miRNAs and were able to validate targets of miR399 and miR159 using 5' RLM-RACE. The conserved miRNAs are predicted to predominately target transcription factors and 93% of the conserved miRNAs originate from intergenic regions. In contrast, only 43% of the novel miRNAs originate from intergenic regions and their predicted targets were more functionally diverse.

CONCLUSION

This study provides important insights into the miRNA gene regulatory networks during narrow-leafed lupin seed development.

INDEL variation in the regulatory region of the major flowering time gene LanFTc1 is associated with vernalization response and flowering time in narrow-leafed lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) cultivation was transformed by 2 dominant vernalization-insensitive, early flowering time loci known as Ku and Julius (Jul), which allowed expansion into shorter season environments. However, reliance on these loci has limited genetic and phenotypic diversity for environmental adaptation in cultivated lupin. We recently predicted that a 1,423-bp deletion in the cis-regulatory region of LanFTc1, a FLOWERING LOCUS T (FT) homologue, derepressed expression of LanFTc1 and was the underlying cause of the Ku phenotype. Here, we surveyed diverse germplasm for LanFTc1 cis-regulatory variation and identified 2 further deletions of 1,208 and 5,162 bp in the 5' regulatory region, which overlap the 1,423-bp deletion. Additionally, we confirmed that no other polymorphisms were perfectly associated with vernalization responsiveness. Phenotyping and gene expression analyses revealed that Jul accessions possessed the 5,162-bp deletion and that the Jul and Ku deletions were equally capable of removing vernalization requirement and up-regulating gene expression. The 1,208-bp deletion was associated with intermediate phenology, vernalization responsiveness, and gene expression and therefore may be useful for expanding agronomic adaptation of lupin. This insertion/deletion series may also help resolve how the vernalization response is mediated at the molecular level in legumes.

Characterization of the genetic factors affecting quinolizidine alkaloid biosynthesis and its response to abiotic stress in narrow-leafed lupin (Lupinus angustifolius L.)

Quinolizidine alkaloids (QAs) are toxic secondary metabolites that complicate the end use of narrow-leafed lupin (NLL; Lupinus angustifolius L.) grain, as levels sometimes exceed the industry limit for its use as a food and feed source. The genotypic and environmental influences on QA production in NLL are poorly understood. Here, the expression of QA biosynthetic genes was analysed in vegetative and reproductive tissues of bitter (high QA) and sweet (low QA) accessions. It was demonstrated that sweet accessions are characterized by lower QA biosynthetic gene expression exclusively in leaf and stem tissues than bitter NLL, consistent with the hypothesis that QAs are predominantly produced in aerial tissues and transported to seeds, rather than synthesized within the seed itself. This analysis informed our identification of additional candidate genes involved in QA biosynthesis. Drought and temperature stress are two major abiotic stresses that often occur during NLL pod set. Hence, we assessed the effect of drought, increased temperature, and their combination, on QA production in three sweet NLL cultivars. A cultivar-specific response to drought and temperature in grain QA levels was observed, including the identification of a cultivar where alkaloid levels did not change with these stress treatments.

Exploring the genetic and adaptive diversity of a pan-Mediterranean crop wild relative: narrow-leafed lupin

KEY MESSAGE

This first pan-Mediterranean analysis of genetic diversity in wild narrow-leafed lupin revealed strong East-West genetic differentiation of populations, an historic eastward migration, and signatures of genetic adaptation to climatic variables. Most grain crops suffer from a narrow genetic base, which limits their potential for adapting to new challenges such as increased stresses associated with climate change. Plant breeders are returning to the wild ancestors of crops and their close relatives to broaden the genetic base of their crops. Understanding the genetic adaptation of these wild relatives will help plant breeders most effectively use available wild diversity. Here, we took narrow-leafed lupin (Lupinus angustifolius L.) as a model to understand adaptation in a wild crop ancestor. A set of 142 wild accessions of narrow-leafed lupin from across the Mediterranean basin were subjected to genotyping-by-sequencing using Diversity Arrays Technology. Phylogenetic, linkage disequilibrium and demographic analyses were employed to explore the history of narrow-leafed lupin within the Mediterranean region. We found strong genetic differentiation between accessions from the western and eastern Mediterranean, evidence of an historic West to East migration, and that eastern Mediterranean narrow-leafed lupin experienced a severe and recent genetic bottleneck. We showed that these two populations differ for flowering time as a result of local adaptation, with the West flowering late while the East flowers early. A genome-wide association study identified single nucleotide polymorphism markers associated with climatic adaptation. Resolving the origin of wild narrow-leafed lupin and how its migration has induced adaptation to specific regions of the Mediterranean serves as a useful resource not only for developing narrow-leafed lupin cultivars with greater resilience to a changing climate, but also as a model which can be applied to other legumes.

A comprehensive draft genome sequence for lupin (Lupinus angustifolius), an emerging health food: insights into plant-microbe interactions and legume evolution

Lupins are important grain legume crops that form a critical part of sustainable farming systems, reducing fertilizer use and providing disease breaks. It has a basal phylogenetic position relative to other crop and model legumes and a high speciation rate. Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is gaining popularity as a health food, which is high in protein and dietary fibre but low in starch and gluten-free. We report the draft genome assembly (609 Mb) of NLL cultivar Tanjil, which has captured >98% of the gene content, sequences of additional lines and a dense genetic map. Lupins are unique among legumes and differ from most other land plants in that they do not form mycorrhizal associations. Remarkably, we find that NLL has lost all mycorrhiza-specific genes, but has retained genes commonly required for mycorrhization and nodulation. In addition, the genome also provided candidate genes for key disease resistance and domestication traits. We also find evidence of a whole-genome triplication at around 25 million years ago in the genistoid lineage leading to Lupinus. Our results will support detailed studies of legume evolution and accelerate lupin breeding programmes.

Quinolizidine Alkaloid Biosynthesis in Lupins and Prospects for Grain Quality Improvement

Quinolizidine alkaloids (QAs) are toxic secondary metabolites found within the genus Lupinus, some species of which are commercially important grain legume crops including Lupinus angustifolius (narrow-leafed lupin, NLL), L. luteus (yellow lupin), L. albus (white lupin), and L. mutabilis (pearl lupin), with NLL grain being the most largely produced of the four species in Australia and worldwide. While QAs offer the plants protection against insect pests, the accumulation of QAs in lupin grain complicates its use for food purposes as QA levels must remain below the industry threshold (0.02%), which is often exceeded. It is not well understood what factors cause grain QA levels to exceed this threshold. Much of the early work on QA biosynthesis began in the 1970-1980s, with many QA chemical structures well-characterized and lupin cell cultures and enzyme assays employed to identify some biosynthetic enzymes and pathway intermediates. More recently, two genes associated with these enzymes have been characterized, however, the QA biosynthetic pathway remains only partially elucidated. Here, we review the research accomplished thus far concerning QAs in lupin and consider some possibilities for further elucidation and manipulation of the QA pathway in lupin crops, drawing on examples from model alkaloid species. One breeding strategy for lupin is to produce plants with high QAs in vegetative tissues while low in the grain in order to confer insect resistance to plants while keeping grain QA levels within industry regulations. With the knowledge achieved on alkaloid biosynthesis in other plant species in recent years, and the recent development of genomic and transcriptomic resources for NLL, there is considerable scope to facilitate advances in our knowledge of QAs, leading to the production of improved lupin crops.

Narrow-Leafed Lupin (Lupinus angustifolius) β1- and β6-Conglutin Proteins Exhibit Antifungal Activity, Protecting Plants against Necrotrophic Pathogen Induced Damage from Sclerotinia sclerotiorum and Phytophthora nicotianae

Vicilins (7S globulins) are seed storage proteins and constitute the main protein family in legume seeds, particularly in narrow-leafed lupin (Lupinus angustifolius L.; NLL), where seven vicilin genes, called β1- to β7-conglutin have been identified. Vicilins are involved in germination processes supplying amino acids for seedling growth and plant development, as well as in some cases roles in plant defense and protection against pathogens. The roles of NLL β-conglutins in plant defense are unknown. Here the potential role of five NLL β-conglutin family members in protection against necrotrophic fungal pathogens was investigated and it was demonstrated that recombinant purified 6xHis-tagged β1- and β6-conglutin proteins exhibited the strongest in vitro growth inhibitory activity against a range of necrotrophic fungal pathogens compared to β2, β3, and β4 conglutins. To examine activity in vivo, two representative necrotrophic pathogens, the fungus Sclerotinia sclerotiorum and oomycete Phytophthora nicotianae were used. Transient expression of β1- and β6-conglutin proteins in Nicotiana benthamiana leaves demonstrated in vivo growth suppression of both of these pathogens, resulting in low percentages of hyphal growth and elongation in comparison to control treated leaves. Cellular studies using β1- and β6-GFP fusion proteins showed these conglutins localized to the cell surface including plasmodesmata. Analysis of cellular death following S. sclerotiorum or P. nicotianae revealed both β1- and β6-conglutins suppressed pathogen induced cell death in planta and prevented pathogen induced suppression of the plant oxidative burst as determined by protein oxidation in infected compared to mock-inoculated leaves.

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.

Comparative genomics and prediction of conditionally dispensable sequences in legume-infecting Fusarium oxysporum formae speciales facilitates identification of candidate effectors

BACKGROUND

Soil-borne fungi of the Fusarium oxysporum species complex cause devastating wilt disease on many crops including legumes that supply human dietary protein needs across many parts of the globe. We present and compare draft genome assemblies for three legume-infecting formae speciales (ff. spp.): F. oxysporum f. sp. ciceris (Foc-38-1) and f. sp. pisi (Fop-37622), significant pathogens of chickpea and pea respectively, the world's second and third most important grain legumes, and lastly f. sp. medicaginis (Fom-5190a) for which we developed a model legume pathosystem utilising Medicago truncatula.

RESULTS

Focusing on the identification of pathogenicity gene content, we leveraged the reference genomes of Fusarium pathogens F. oxysporum f. sp. lycopersici (tomato-infecting) and F. solani (pea-infecting) and their well-characterised core and dispensable chromosomes to predict genomic organisation in the newly sequenced legume-infecting isolates. Dispensable chromosomes are not essential for growth and in Fusarium species are known to be enriched in host-specificity and pathogenicity-associated genes. Comparative genomics of the publicly available Fusarium species revealed differential patterns of sequence conservation across F. oxysporum formae speciales, with legume-pathogenic formae speciales not exhibiting greater sequence conservation between them relative to non-legume-infecting formae speciales, possibly indicating the lack of a common ancestral source for legume pathogenicity. Combining predicted dispensable gene content with in planta expression in the model legume-infecting isolate, we identified small conserved regions and candidate effectors, four of which shared greatest similarity to proteins from another legume-infecting ff. spp.

CONCLUSIONS

We demonstrate that distinction of core and potential dispensable genomic regions of novel F. oxysporum genomes is an effective tool to facilitate effector discovery and the identification of gene content possibly linked to host specificity. While the legume-infecting isolates didn't share large genomic regions of pathogenicity-related content, smaller regions and candidate effector proteins were highly conserved, suggesting that they may play specific roles in inducing disease on legume hosts.

Characterization and mapping of LanrBo: a locus conferring anthracnose resistance in narrow-leafed lupin (Lupinus angustifolius L.)

A novel and highly effective source of anthracnose resistance in narrow-leafed lupin was identified. Resistance was shown to be governed by a single dominant locus. Molecular markers have been developed, which can be used for selecting resistant genotypes in lupin breeding. A screening for anthracnose resistance of a set of plant genetic resources of narrow-leafed lupin (Lupinus angustifolius L.) identified the breeding line Bo7212 as being highly resistant to anthracnose (Colletotrichum lupini). Segregation analysis indicated that the resistance of Bo7212 is inherited by a single dominant locus. The corresponding resistance gene was given the designation LanrBo. Previously published molecular anchor markers allowed us to locate LanrBo on linkage group NLL-11 of narrow-leafed lupin. Using information from RNAseq data obtained with inoculated resistant vs. susceptible lupin entries as well as EST-sequence information from the model genome Lotus japonicus, additional SNP and EST markers linked to LanrBo were derived. A bracket of two LanrBo-flanking markers allows for precise marker-assisted selection of the novel resistance gene in narrow-leafed lupin breeding programs.

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.

The Arabidopsis KH-Domain RNA-Binding Protein ESR1 Functions in Components of Jasmonate Signalling, Unlinking Growth Restraint and Resistance to Stress

Glutathione S-transferases (GSTs) play important roles in the protection of cells against toxins and oxidative damage where one Arabidopsis member, GSTF8, has become a commonly used marker gene for early stress and defense responses. A GSTF8 promoter fragment fused to the luciferase reporter gene was used in a forward genetic screen for Arabidopsis mutants with up-regulated GSTF8 promoter activity. This identified the esr1-1 (enhanced stress response 1) mutant which also conferred increased resistance to the fungal pathogen Fusarium oxysporum. Through positional cloning, the ESR1 gene was found to encode a KH-domain containing RNA-binding protein (At5g53060). Whole transcriptome sequencing of esr1-1 identified altered expression of genes involved in responses to biotic and abiotic stimuli, hormone signaling pathways and developmental processes. In particular was an overall significant enrichment for jasmonic acid (JA) mediated processes in the esr1-1 down-regulated dataset. A subset of these genes were tested for MeJA inducibility and we found the expression of some but not all were reduced in esr1-1. The esr1-1 mutant was not impaired in other aspects of JA-signalling such as JA- sensitivity or development, suggesting ESR1 functions in specific components of the JA-signaling pathway. Examination of salicylic acid (SA) regulated marker genes in esr1-1 showed no increase in basal or SA induced expression suggesting repression of JA-regulated genes is not due to antagonistic SA-JA crosstalk. These results define new roles for KH-domain containing proteins with ESR1 unlinking JA-mediated growth and defense responses.

Analysis of conglutin seed storage proteins across lupin species using transcriptomic, protein and comparative genomic approaches

BACKGROUND

The major proteins in lupin seeds are conglutins that have primary roles in supplying carbon, sulphur and nitrogen and energy for the germinating seedling. They fall into four families; α, β, γ and δ. Interest in these conglutins is growing as family members have been shown to have beneficial nutritional and pharmaceutical properties.

RESULTS

An in-depth transcriptome and draft genome from the narrow-leafed lupin (NLL; Lupinus angustifolius) variety, Tanjil, were examined and 16 conglutin genes were identified. Using RNAseq data sets, the structure and expression of these 16 conglutin genes were analysed across eight lupin varieties from five lupin species. Phylogenic analysis suggest that the α and γ conglutins diverged prior to lupin speciation while β and δ members diverged both prior and after speciation. A comparison of the expression of the 16 conglutin genes was performed, and in general the conglutin genes showed similar levels of RNA expression among varieties within species, but quite distinct expression patterns between lupin species. Antibodies were generated against the specific conglutin families and immunoblot analyses were used to compare the levels of conglutin proteins in various tissues and during different stages of seed development in NLL, Tanjil, confirming the expression in the seed. This analysis showed that the conglutins were expressed highly at the mature seed stage, in all lupin species, and a range of polypeptide sizes were observed for each conglutin family.

CONCLUSIONS

This study has provided substantial information on the complexity of the four conglutin families in a range of lupin species in terms of their gene structure, phylogenetic relationships as well as their relative RNA and protein abundance during seed development. The results demonstrate that the majority of the heterogeneity of conglutin polypeptides is likely to arise from post-translational modification from a limited number of precursor polypeptides rather than a large number of different genes. Overall, the results demonstrate a high degree of plasticity for conglutin expression during seed development in different lupin species.

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.

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.

Characterization and genetic dissection of resistance to spotted alfalfa aphid (Therioaphis trifolii) in Medicago truncatula

Aphids cause significant yield losses in agricultural crops worldwide. Medicago truncatula, a model legume, cultivated pasture species in Australia and close relative of alfalfa (Medicago sativa), was used to study the defence response against Therioaphis trifolii f. maculate [spotted alfalfa aphid (SAA)]. Aphid performance and plant damage were compared among three accessions. A20 is highly susceptible, A17 has moderate resistance, and Jester is strongly resistant. Subsequent analyses using A17 and A20, reciprocal F1s and an A17×A20 recombinant inbred line (RIL) population revealed that this moderate resistance is phloem mediated and involves antibiosis and tolerance but not antixenosis. Electrical penetration graph analysis also identified a novel waveform termed extended potential drop, which occurred following SAA infestation of M. truncatula. Genetic dissection using the RIL population revealed three quantitative trait loci on chromosomes 3, 6, and 7 involved in distinct modes of aphid defence including antibiosis and tolerance. An antibiosis locus resides on linkage group 3 (LG3) and is derived from A17, whereas a plant tolerance and antibiosis locus resides on LG6 and is derived from A20, which exhibits strong temporary tolerance. The loci identified reside in regions harbouring classical resistance genes, and introgression of these loci in current medic cultivars may help provide durable resistance to SAA, while elucidation of their molecular mechanisms may provide valuable insight into other aphid-plant interactions.

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.

Development of genomic resources for the narrow-leafed lupin (Lupinus angustifolius): construction of a bacterial artificial chromosome (BAC) library and BAC-end sequencing

Background

Lupinus angustifolius L, also known as narrow-leafed lupin (NLL), is becoming an important grain legume crop that is valuable for sustainable farming and is becoming recognised as a potential human health food. Recent interest is being directed at NLL to improve grain production, disease and pest management and health benefits of the grain. However, studies have been hindered by a lack of extensive genomic resources for the species.

Results

A NLL BAC library was constructed consisting of 111,360 clones with an average insert size of 99.7 Kbp from cv Tanjil. The library has approximately 12 × genome coverage. Both ends of 9600 randomly selected BAC clones were sequenced to generate 13985 BAC end-sequences (BESs), covering approximately 1% of the NLL genome. These BESs permitted a preliminary characterisation of the NLL genome such as organisation and composition, with the BESs having approximately 39% G:C content, 16.6% repetitive DNA and 5.4% putative gene-encoding regions. From the BESs 9966 simple sequence repeat (SSR) motifs were identified and some of these are shown to be potential markers.

Conclusions

The NLL BAC library and BAC-end sequences are powerful resources for genetic and genomic research on lupin. These resources will provide a robust platform for future high-resolution mapping, map-based cloning, comparative genomics and assembly of whole-genome sequencing data for the species.

Identification of potential early regulators of aphid resistance in Medicago truncatula via transcription factor expression profiling

*Resistance to aphids has been identified in a number of plant species, yet the molecular mechanisms underlying aphid resistance remain largely unknown. *Using high-throughput quantitative real-time PCR technology, the transcription profiles of 752 putative Medicago truncatula transcription factor genes were analysed in a pair of susceptible and resistant closely related lines of M. truncatula following 6 and 12 h of bluegreen aphid (Acyrthosiphon kondoi) infestation. *Eighty-two transcription factor genes belonging to 30 transcription factor families were responsive to bluegreen aphid infestation. More transcription factor genes were responsive in the resistant interaction than in the susceptible interaction; of the 36 genes that were induced at 6 and/or 12 h, 32 were induced only in the resistant interaction. Bluegreen aphid-induced expression of a subset of these genes was correlated with the presence of AKR, a single dominant gene conferring resistance to bluegreen aphids. Similar transcription factor expression patterns of this subset were associated with bluegreen aphid resistance in other M. truncatula genetic backgrounds, as well as with resistance to pea aphid (Acyrthosiphon pisum). *Our results suggest that these transcription factors are among the early aphid-responsive genes in resistant plants, and may play important roles in resistance to multiple aphid species.

Two independent resistance genes in the Medicago truncatula cultivar jester confer resistance to two different aphid species of the genus Acyrthosiphon

In recent years the biology of resistance to aphids and other sap-sucking insects has been studied in detail, whereas the genetic basis underlying this resistance is still poorly understood. Genetic resistance to Acyrthosiphon kondoi Shinji (bluegreen aphid; BGA) has been identified in Medicago truncatula Gaertn and backcrossed into susceptible cultivars. One of these M. truncatula cultivars named Jester also has good resistance to an Australian biotype of pea aphid (PA; A. pisum Harris). Until now it has been unclear whether resistance to each aphid species of the genus Acyrthosiphon is conferred by distinct genes or the same gene termed AKR for A. kondoi resistance. Infestation of the progenitors of the cultivar Jester with both aphid species revealed that resistance to BGA came from a different donor than resistance to PA, demonstrating that resistance to these aphid species is mediated by different resistance genes. However, an interaction between these genes for resistance to Acyrthosiphon species remains a possibility, given that PA resistance was not one of the parameters selected for in the creation of Jester. The identification of resistance to the model aphid, PA, and a closely related aphid BGA in the same genetic background of the model legume M. truncatula makes this system an attractive model for the study of plant-aphid interactions, as well as R gene specificity and evolution.

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.

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.

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.

Agroinfection-based high-throughput screening reveals specific recognition of INF elicitins in Solanum

SUMMARY We adapted and optimized the use of the Agrobacterium tumefaciens binary PVX expression system (PVX agroinfection) to screen Solanum plants for response to pathogen elicitors and applied the assay to identify a total of 11 clones of Solanum huancabambense and Solanum microdontum, out of 31 species tested, that respond to the elicitins INF1, INF2A and INF2B of Phytophthora infestans. Prior to this study, response to INF elicitins was only known in Nicotiana spp. within the Solanaceae. The identified S. huancabambense and S. microdontum clones also exhibited hypersensitivity-like cell death following infiltration with purified recombinant INF1, INF2A and INF2B, thereby validating the screening protocol. Comparison of INF elicitin activity revealed that Nicotiana plants responded to significantly lower concentrations than Solanum, suggesting variable levels of sensitivity to INF elicitins. We exploited natural variation in response to INF elicitins in the identified Solanum accessions to evaluate the relationship between INF recognition and late blight resistance. Interestingly, several INF-responsive Solanum plants were susceptible to P. infestans. Also, an S. microdontum xSolanum tuberosum (potato) population that segregates for INF response was generated but failed to identify a measurable contribution of INF response to resistance. These results suggest that in Solanum, INF elicitins are recognized as general elicitors and do not have a measurable contribution to disease resistance.

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

The world's oldest and largest Medicago truncatula collection is housed at the South Australian Research and Development Institute (SARDI). We used six simple sequence repeat (SSR) loci to analyse the genetic diversity and relationships between randomly selected individuals from 192 accessions in the core collection. M. truncatula is composed of three subspecies (ssp.): ssp. truncatula, ssp. longeaculeata, and ssp. tricycla. Analysis at the level of six SSR loci supports the concept of ssp. tricycla, all the samples of which showed unique alleles at two loci. Contingency Chi-squared tests were significant between ssp. tricycla and ssp. truncatula at four loci, suggesting a barrier to gene flow between these subspecies. In accessions defined as ssp. longeaculeata, no unique allelic distribution or diagnostic sizes were observed, suggesting this apparent ssp. is a morphological variant of ssp. truncatula. The data also suggest M. truncatula that exhibits unusually wide genotype dispersal throughout its native Mediterranean region, possibly due to animal and trade-related movements. Our results showed the collection to be highly diverse, exhibiting an average of 25 SSR alleles per locus, with over 90% of individuals showing discrete genotypes. The rich diversity of the SARDI collection provides an invaluable resource for studying natural allelic variation of M. truncatula. To efficiently exploit the variation in the SARDI collection, we have defined a subset of accessions (n = 61) that maximises the diversity.

A cDNA-AFLP based strategy to identify transcripts associated with avirulence in Phytophthora infestans

Expression profiling using cDNA-AFLP is commonly used to display the transcriptome of a specific tissue or developmental stage. Here, cDNA-AFLP was used to identify transcripts in a segregating F1 population of Phytophthora infestans, the oomycete pathogen that causes late blight. To find transcripts derived from putative avirulence (Avr) genes germinated cyst cDNA from F1 progeny with defined avirulence phenotypes was pooled and used in a bulked segregant analysis (BSA). Over 30,000 transcript derived fragments (TDFs) were screened resulting in 99 Avr-associated TDFs as well as TDFs with opposite pattern. With 142 TDF sequences homology searches and database mining was carried out. cDNA-AFLP analysis on individual F1 progeny revealed 100% co-segregation of four TDFs with particular AVR phenotypes and this was confirmed by RT-PCR. Two match the same P. infestans EST with unknown sequence and this is a likely candidate for Avr4. The other two are associated with the Avr3b-Avr10-Avr11 locus. This combined cDNA-AFLP/BSA strategy is an efficient approach to identify Avr-associated transcriptome markers that can complement positional cloning.