Allelic diversity in the transcriptomes of contrasting rust-infected genotypes of Lathyrus sativus, a lasting resource for smart breeding

Functional and evolutionary study of MLO gene family in the regulation of Sclerotinia stem rot resistance in Brassica napus L

BACKGROUND

Oilseed rape (Brassica napus L.) is known as one of the most important oilseed crops cultivated around the world. However, its production continuously faces a huge challenge of Sclerotinia stem rot (SSR), a destructive disease caused by the fungus Sclerotinia sclerotiorum, resulting in huge yield loss annually. The SSR resistance in B. napus is quantitative and controlled by a set of minor genes. Identification of these genes and pyramiding them into a variety are a major strategy for SSR resistance breeding in B. napus.

RESULTS

Here, we performed a genome-wide association study (GWAS) using a natural population of B. napus consisting of 222 accessions to identify BnaA08g25340D (BnMLO2_2) as a candidate gene that regulates the SSR resistance. BnMLO2_2 was a member of seven homolog genes of Arabidopsis Mildew Locus O 2 (MLO2) and the significantly SNPs were mainly distributed in the promoter of BnMLO2_2, suggesting a role of BnMLO2_2 expression level in the regulation of SSR resistance. We expressed BnMLO2_2 in Arabidopsis and the transgenic plants displayed an enhanced SSR resistance. Transcriptome profiling of different tissues of B. napus revealed that BnMLO2_2 had the most expression level in leaf and silique tissues among all the 7 BnMLO2 members and also expressed higher in the SSR resistant accession than in the susceptible accession. In Arabidopsis, mlo2 plants displayed reduced resistance to SSR, whereas overexpression of MLO2 conferred plants an enhanced SSR resistance. Moreover, a higher expression level of MLO2 showed a stronger SSR resistance in the transgenic plants. The regulation of MLO2 in SSR resistance may be associated with the cell death. Collinearity and phylogenetic analysis revealed a large expansion of MLO family in Brassica crops.

CONCLUSION

Our study revealed an important role of BnMLO2 in the regulation of SSR resistance and provided a new gene candidate for future improvement of SSR resistance in B. napus and also new insights into understanding of MLO family evolution in Brassica crops.

Genome-wide characterization and expression analysis of the MLO gene family sheds light on powdery mildew resistance in Lagenaria siceraria

MLO (mildew locus O) genes play a vital role in plant disease defense system, especially powdery mildew (PM). Lagenaria siceraria is a distinct Cucurbitaceae crop, and PM is one of the most serious diseases threatening crop production and quality. Although MLOs have been exploited in many Cucurbitaceae species, genome-wide mining of MLO gene family in bottle gourd has not been surveyed yet. Here we identified 16 MLO genes in our recently assembled L. siceraria genome. A total of 343 unique MLO protein sequences from 20 species were characterized and compared to deduce a generally high level of purifying selection and the occurrence of regions related to candidate susceptibility factors in the evolutional divergence. LsMLOs were clustered in six clades containing seven conserved transmembrane domains and 10 clade-specific motifs along with deletion and variation. Three genes (LsMLO3, LsMLO6, and LsMLO13) in clade V showed high sequence identity with orthologues involved in PM susceptibility. The expression pattern of LsMLOs was tissue-specific but not cultivar-specific. Furthermore, it was indicated by qRT-PCR and RNA-seq that LsMLO3 and LsMLO13 were highly upregulated in response to PM stress. Subsequent sequence analysis revealed the structural deletion of LsMLO13 and a single nonsynonymous substitution of LsMLO3 in the PM-resistant genotype. Taken all together, it is speculated that LsMLO13 is likely a major PM susceptibility factor. The results of this study provide new insights into MLO family genes in bottle gourd and find a potential candidate S gene for PM tolerance breeding.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Association Mapping of Lathyrus sativus Disease Response to Uromyces pisi Reveals Novel Loci Underlying Partial Resistance

Uromyces pisi ([Pers.] D.C.) Wint. is an important foliar biotrophic pathogen infecting grass pea (Lathyrus sativus L.), compromising their yield stability. To date, few efforts have been made to assess the natural variation in grass pea resistance and to identify the resistance loci operating against this pathogen, limiting its efficient breeding exploitation. To overcome this knowledge gap, the genetic architecture of grass pea resistance to U. pisi was investigated using a worldwide collection of 182 accessions through a genome-wide association approach. The response of the grass pea collection to rust infection under controlled conditions and at the seedling stage did not reveal any hypersensitive response but a continuous variation for disease severity, with the identification of promising sources of partial resistance. A panel of 5,651 high-quality single-nucleotide polymorphism (SNP) markers previously generated was used to test for SNP-trait associations, based on a mixed linear model accounting for population structure. We detected seven SNP markers significantly associated with U. pisi disease severity, suggesting that partial resistance is oligogenic. Six of the associated SNP markers were located in chromosomes 4 and 6, while the remaining SNP markers had no known chromosomal position. Through comparative mapping with the pea reference genome, a total of 19 candidate genes were proposed, encoding for leucine-rich repeat, NB-ARC domain, and TGA transcription factor family, among others. Results presented in this study provided information on the availability of partial resistance in grass pea germplasm and advanced our understanding of the molecular mechanisms of quantitative resistance to rust in grass pea. Moreover, the detected associated SNP markers constitute promising genomic targets for the development of molecular tools to assist disease resistance precision breeding.

Contrasting β-ODAP content correlates with stress gene expression in Lathyrus cultivars

Lathyrus sativus, commonly known as grass pea, is a nutrient-rich pulse crop with remarkable climate-resilient attributes. However, wide use of this nutritious crop is not adopted owing to the presence of a non-protein amino acid β-N-oxalyl-l-α,β-diaminopropionic acid (β-ODAP), which is neurotoxic if consumed in large quantities. We conducted a de novo transcriptomic profiling of two ODAP contrasting cultivars, Pusa-24 and its somaclonal variant Ratan, to understand the genetic changes leading to and associated with β-ODAP levels. Differential gene expression analysis showed that a variety of genes are downregulated in low β-ODAP cultivar Ratan and include genes involved in biotic/abiotic stress tolerance, redox metabolism, hormonal metabolism, and sucrose, and starch metabolism. Several genes related to chromatin remodeling are differentially expressed in cultivar Ratan. β-ODAP biosynthetic genes in these cultivars showed differential upregulation upon stress. ODAP content of these cultivars varied differentially upon stress and development. Physiological experiments indicate reduced relative water content and perturbed abscisic acid levels in the low ODAP cultivar. Altogether, our results suggest that the low ODAP cultivar may have a reduced stress tolerance. The dataset provides insight into the biological role of ODAP and will be helpful for hypothesis-driven experiments to understand ODAP biosynthesis and regulation.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

Grass pea natural variation reveals oligogenic resistance to Fusarium oxysporum f. sp. pisi

Grass pea (Lathyrus sativus L.) is an annual legume species, phylogenetically close to pea (Pisum sativum L.), that may be infected by Fusarium oxysporum f. sp. pisi (Fop), the causal agent of fusarium wilt in peas with vast worldwide yield losses. A range of responses varying from high resistance to susceptibility to this pathogen has been reported in grass pea germplasm. Nevertheless, the genetic basis of that diversity of responses is still unknown, hampering its breeding exploitation. To identify genomic regions controlling grass pea resistance to fusarium wilt, a genome-wide association study approach was applied on a grass pea worldwide collection of accessions inoculated with Fop race 2. Disease responses were scored in this collection that was also subjected to high-throughput based single nucleotide polymorphisms (SNP) screening through genotyping-by-sequencing. A total of 5,651 high-quality SNPs were considered for association mapping analysis, performed using mixed linear models accounting for population structure. Because of the absence of a fully assembled grass pea reference genome, SNP markers' genomic positions were retrieved from the pea's reference genome v1a. In total, 17 genomic regions were associated with three fusarium wilt response traits in grass pea, anticipating an oligogenic control. Seven of these regions were located on pea chromosomes 1, 6, and 7. The candidate genes underlying these regions were putatively involved in secondary and amino acid metabolism, RNA (regulation of transcription), transport, and development. This study revealed important fusarium wilt resistance favorable grass pea SNP alleles, allowing the development of molecular tools for precision disease resistance breeding.

Current Perspectives on Reducing the β-ODAP Content and Improving Potential Agronomic Traits in Grass Pea (Lathyrus sativus L.)

Grass pea is well-established as one of the most resilient and versatile crops that can thrive under extreme climatic circumstances such as cold, heat, drought, salt-affected soils, submergence, and excessive rainfall along with resistance to several diseases and pests. However, despite the awareness of its virtues, its cultivation globally has decreased recently owing to the presence of a neurotoxin, β-N-oxalyl-L-α, β-diaminopropionic acid (β-ODAP), in the seedlings and seeds of this legume, which has been reported to cause neurolathyrism, a non-reversible neurological disorder in humans and animals. Significant repositories of Lathyrus germplasm are available across countries that have provided access to a wide range of agro-morphological traits as well as the low β ODAP content. Efforts have been made worldwide to use these germplasms for the genetic enhancement of grass pea to make this food safe for human consumption. Efforts on molecular breeding of this crop are also lagging. However, during the last decade, the research scenario has changed with some efforts being made toward improving this climate resilient pulse in terms of genomic resources. Molecular markers have also been used to evaluate the interspecific diversity as well as the phylogenetic relationship among the species and mapping studies. Intron-targeted amplified polymorphic, genomic simple sequence repeat, resistance genes analogs, and disease resistance markers developed for other legume species have been successfully cross-amplified in grass pea. Transcriptomic studies have recently been undertaken on grass pea by deploying several second-generation sequencing techniques. In addition, a few studies have attempted to unveil the genes and the underlying mechanism conferring biotic and abiotic stress or regulating the pathway of β-ODAP in grass pea. Proteomics has accelerated the identification studies on differential proteomes in response to salinity and low-temperature stress conditions for unveiling the common signaling pathways involved in mitigating these abiotic stresses and in discovering differentially regulated proteins. In grass pea, a metabolomics approach has been used to identify the metabolic processes associated with β-ODAP synthesis. Genome sequencing of grass pea is under way which is expected to be vital for whole-genome re-sequencing and gene annotation toward the identification of genes with novel functions. Recently, a draft genome sequence of grass pea was developed, and some efforts are underway to re-sequence a diverse panel of grass pea comprising 384 germplasm lines. Owing to the scantiness of a successful transformation protocol, research on the application of modern approaches of genome editing like the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) or CRISPR-associated protein 9 (CRISPR/Cas9) system for the engineering of signaling pathways or regulatory mechanisms seeks immediate attention to reduce the β-ODAP content in seeds and to improve the potential agronomic traits in grass pea.

The MLO1 powdery mildew susceptibility gene in Lathyrus species: The power of high-density linkage maps in comparative mapping and synteny analysis

Powdery mildews are major diseases for a range of crops. The loss of function of specific Mildew Locus O (MLO) genes has long been associated with pre-haustorial plant resistance to powdery mildew and has proven to be durable in several species. Erysiphe pisi is the major causal agent of powdery mildew in pea (Pisum sativum L.) and in the closely related Lathyrus sativus L. and Lathyrus cicera L. PsMLO1 has been extensively studied in pea. However, no MLO gene family members have been isolated and characterized in Lathyrus species so far. In this study, MLO1 genes were isolated and characterized in L. sativus and L. cicera genotypes with varied levels of partial resistance against powdery mildew. Phylogenetic analyses confirmed that Lathyrus MLO1 belongs to Clade V, like all dicot MLO proteins associated with powdery mildew susceptibility. A L. sativus recombinant inbred line population (RIL) was genotyped by sequencing to develop a high-density L. sativus genetic linkage map. DNA sequence polymorphisms between the analyzed genotypes allowed the location of MLO1 in the newly developed L. sativus RIL genetic linkage map. Subsequent comparative mapping between L. sativus and L. cicera genetic maps and P. sativum, Lens culinaris Medik., and Medicago truncatula Gaertn. reference genomes revealed important aspects of the conservation of the MLO1 locus position and of the overall chromosomal rearrangements occurring during legume evolution, with relevance to legume disease resistance breeding programs.

Ribosomal DNA localization on Lathyrus species chromosomes by FISH

BACKGROUND

Fluorescence In Situ Hybridization (FISH) played an essential role to locate the ribosomal RNA genes on the chromosomes that offered a new tool to study the chromosome structure and evolution in plant. The 45S and 5S rRNA genes are independent and localized at one or more loci per the chromosome complement, their positions along chromosomes offer useful markers for chromosome discriminations. In the current study FISH has been performed to locate 45S and 5S rRNA genes on the chromosomes of nine Lathyrus species belong to five different sections, all have chromosome number 2n=14, Lathyrus gorgoni Parl, Lathyrus hirsutus L., Lathyrus amphicarpos L., Lathyrus odoratus L., Lathyrus sphaericus Retz, Lathyrus incospicuus L, Lathyrus paranensis Burkart, Lathyrus nissolia L., and Lathyrus articulates L.

RESULTS

The revealed loci of 45S and 5S rDNA by FISH on metaphase chromosomes of the examined species were as follow: all of the studied species have one 45S rDNA locus and one 5S rDNA locus except L. odoratus L., L. amphicarpos L. and L. sphaericus Retz L. have two loci of 5S rDNA. Three out of the nine examined species have the loci of 45S and 5S rRNA genes on the opposite arms of the same chromosome (L. nissolia L., L. amphicarpos L., and L. incospicuus L.), while L. hirsutus L. has both loci on the same chromosome arm. The other five species showed the loci of the two types of rDNA on different chromosomes.

CONCLUSION

The detected 5S and 45S rDNA loci in Lathyrus could be used as chromosomal markers to discriminate the chromosome pairs of the examined species. FISH could discriminate only one chromosome pair out of the seven pairs in three species, in L. hirsutus L., L. nissolia L. and L. incospicuus L., and two chromosome pairs in five species, in L. paranensis Burkart, L. odoratus L., L. amphicarpos L., L. gorgoni Parl. and L. articulatus L., while it could discriminate three chromosome pairs in L. sphaericus Retz. these results could contribute into the physical genome mapping of Lathyrus species and the evolution of rDNA patterns by FISH in the coming studies in future.

Grass pea (Lathyrus sativus L.): orphan crop, nutraceutical or just plain food?

Although grass pea is an environmentally successful robust legume with major traits of interest for food and nutrition security, the genetic potential of this orphan crop has long been neglected. Grass pea (Lathyrus sativus L.) is a Neolithic plant that has survived millennia of cultivation and has spread over three continents. It is a robust legume crop that is considered one of the most resilient to climate changes and to be survival food during drought-triggered famines. The hardy penetrating root system allows the cultivation of grass pea in various soil types, including marginal ones. As an efficient nitrogen fixer, it meets its own nitrogen requirements and positively benefits subsequent crops. However, already in ancient India and Greece, overconsumption of the seeds and a crippling neurological disorder, later coined neurolathyrism, had been linked. Overemphasis of their suspected toxic properties has led to disregard the plant's exceptionally positive agronomic properties and dietary advantages. In normal socio-economic and environmental situations, in which grass pea is part of a balanced diet, neurolathyrism is virtually non-existent. The etiology of neurolathyrism has been oversimplified and the deficiency in methionine in the diet has been overlooked. In view of the global climate change, this very adaptable and nutritious orphan crop deserves more attention. Grass pea can become a wonder crop if the double stigma on its reputation as a toxic plant and as food of the poor can be disregarded. Additionally, recent research has exposed the potential of grass pea as a health-promoting nutraceutical. Development of varieties with an improved balance in essential amino acids and diet may be relevant to enhance the nutritional value without jeopardizing the multiple stress tolerance of this promising crop.

Biotic and Abiotic Constraints in Mungbean Production-Progress in Genetic Improvement

Mungbean [Vigna radiata (L.) R. Wilczek var. radiata] is an important food and cash legume crop in Asia. Development of short duration varieties has paved the way for the expansion of mungbean into other regions such as Sub-Saharan Africa and South America. Mungbean productivity is constrained by biotic and abiotic factors. Bruchids, whitefly, thrips, stem fly, aphids, and pod borers are the major insect-pests. The major diseases of mungbean are yellow mosaic, anthracnose, powdery mildew, Cercospora leaf spot, halo blight, bacterial leaf spot, and tan spot. Key abiotic stresses affecting mungbean production are drought, waterlogging, salinity, and heat stress. Mungbean breeding has been critical in developing varieties with resistance to biotic and abiotic factors, but there are many constraints still to address that include the precise and accurate identification of resistance source(s) for some of the traits and the traits conferred by multi genes. Latest technologies in phenotyping, genomics, proteomics, and metabolomics could be of great help to understand insect/pathogen-plant, plant-environment interactions and the key components responsible for resistance to biotic and abiotic stresses. This review discusses current biotic and abiotic constraints in mungbean production and the challenges in genetic improvement.

First genetic linkage map of Lathyrus cicera based on RNA sequencing-derived markers: Key tool for genetic mapping of disease resistance

The Lathyrus cicera transcriptome was analysed in response to rust (Uromyces pisi) infection to develop novel molecular breeding tools with potential for genetic mapping of resistance in this robust orphan legume species. One RNA-seq library each was generated from control and rust-inoculated leaves from two L. cicera genotypes with contrasting quantitative resistance, de novo assembled into contigs and sequence polymorphisms were identified. In toto, 19,224 SNPs differentiate the susceptible from the partially resistant genotype's transcriptome. In addition, we developed and tested 341 expressed E-SSR markers from the contigs, of which 60.7% varied between the two L. cicera genotypes. A first L. cicera linkage map was created using part of the developed markers in a RIL population from the cross of the two genotypes. This map contains 307 markers, covered 724.2 cM and is organised in 7 major and 2 minor linkage groups, with an average mapping interval of 2.4 cM. The genic markers also enabled us to compare their position in L. cicera map with the physical position of the same markers mapped on Medicago truncatula genome, highlighting a high macrosyntenic conservation between both species. This study provides a large new set of genic polymorphic molecular markers with potential for mapping rust resistances. It represents the first step towards genomics-assisted precision breeding in L. cicera.

An RNA Sequencing Transcriptome Analysis of Grasspea (Lathyrus sativus L.) and Development of SSR and KASP Markers

Grasspea (Lathyrus sativus L., 2n = 14) has great agronomic potential because of its ability to survive under extreme conditions, such as drought and flood. However, this legume is less investigated because of its sparse genomic resources and very slow breeding process. In this study, 570 million quality-filtered and trimmed cDNA sequence reads with total length of over 82 billion bp were obtained using the Illumina NextSeq^TM 500 platform. Approximately two million contigs and 142,053 transcripts were assembled from our RNA-Seq data, which resulted in 27,431 unigenes with an average length of 1,250 bp and maximum length of 48,515 bp. The unigenes were of high-quality. For example, the stay-green (SGR) gene of grasspea was aligned with the SGR gene of pea with high similarity. Among these unigenes, 3,204 EST-SSR primers were designed, 284 of which were randomly chosen for validation. Of these validated unigenes, 87 (30.6%) EST-SSR primers produced polymorphic amplicons among 43 grasspea accessions selected from different geographical locations. Meanwhile, 146,406 SNPs were screened and 50 SNP loci were randomly chosen for the kompetitive allele-specific PCR (KASP) validation. Over 80% (42) SNP loci were successfully transformed to KASP markers. Comparison of the dendrograms according to the SSR and KASP markers showed that the different marker systems are partially consistent with the dendrogram constructed in our study.

Genome-wide identification and comparison of legume MLO gene family

MLO proteins are highly conserved proteins with seven trans-membrane domains. Specific MLO genes have been linked to plant disease susceptibility. Others are involved in plant reproduction and in root thigmomorphogenesis. Functions of the remaining MLOs are still unknown. Here we performed a genome-wide survey of the MLO family in eight legume species from different clades of the Papillionoideae sub-family. A total of 118 MLO sequences were identified and characterized. Their deduced protein sequences shared the characteristics of MLO proteins. The total number of MLO genes per legume species varied from 13 to 20 depending on the species. Legume MLOs were evenly distributed over their genomes and tended to localize within syntenic blocks conserved across legume genomes. Phylogenetic analysis indicated that these sequences clustered in seven well-defined clades. Comparison of MLO protein sequences revealed 34 clade-specific motifs in the variable regions of the proteins. Comparative analyses of the MLO family between legume species also uncovered several evolutionary differences between the tropical legume species from the Phaseoloid clades and the other legume species. Altogether, this study provides interesting new features on the evolution of the MLO family. It also provides valuable clues to identify additional MLO genes from non-sequenced species.

Identification of Genes in a Partially Resistant Genotype of Avena sativa Expressed in Response to Puccinia coronata Infection

Cultivated oat (Avena sativa), an important crop in many countries, can suffer significant losses through infection by the fungus Puccinia coronata, the causal agent of crown rust disease. Understanding the molecular basis of existing partial resistance to this disease might provide targets of interest for crop improvement programs. A suppressive subtractive hybridization (SSH) library was constructed using cDNA from the partially resistant oat genotype MN841801-1 after inoculation with the pathogen. A total of 929 genes returned a BLASTx hit and were annotated under different GO terms, including 139 genes previously described as participants in mechanisms related to the defense response and signal transduction. Among these were genes involved in pathogen recognition, cell-wall modification, oxidative burst/ROS scavenging, and abscisic acid biosynthesis, as well genes related to inducible defense responses mediated by salicylic and jasmonic acid (although none of which had been previously reported involved in strong responses). These findings support the hypothesis that basal defense mechanisms are the main systems operating in oat partial resistance to P. coronata. When the expression profiles of 20 selected genes were examined at different times following inoculation with the pathogen, the partially resistant genotype was much quicker in mounting a response than a susceptible genotype. Additionally, a number of genes not previously described in oat transcriptomes were identified in this work, increasing our molecular knowledge of this crop.

Lathyrus sativus transcriptome resistance response to Ascochyta lathyri investigated by deepSuperSAGE analysis

Lathyrus sativus (grass pea) is a temperate grain legume crop with a great potential for expansion in dry areas or zones that are becoming more drought-prone. It is also recognized as a potential source of resistance to several important diseases in legumes, such as ascochyta blight. Nevertheless, the lack of detailed genomic and/or transcriptomic information hampers further exploitation of grass pea resistance-related genes in precision breeding. To elucidate the pathways differentially regulated during ascochyta-grass pea interaction and to identify resistance candidate genes, we compared the early response of the leaf gene expression profile of a resistant L. sativus genotype to Ascochyta lathyri infection with a non-inoculated control sample from the same genotype employing deepSuperSAGE. This analysis generated 14.387 UniTags of which 95.7% mapped to a reference grass pea/rust interaction transcriptome. From the total mapped UniTags, 738 were significantly differentially expressed between control and inoculated leaves. The results indicate that several gene classes acting in different phases of the plant/pathogen interaction are involved in the L. sativus response to A. lathyri infection. Most notably a clear up-regulation of defense-related genes involved in and/or regulated by the ethylene pathway was observed. There was also evidence of alterations in cell wall metabolism indicated by overexpression of cellulose synthase and lignin biosynthesis genes. This first genome-wide overview of the gene expression profile of the L. sativus response to ascochyta infection delivered a valuable set of candidate resistance genes for future use in precision breeding.