Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

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

BMC Plant Biology volume 14, Article number: 376 (2014) Cite this article

Abstract

Background

Grass pea (Lathyrus sativus L.) is a valuable resource for potentially durable partial resistance to rust. To gain insight into the resistance mechanism and identify potential resistance genes, we generated the first comprehensive transcriptome assemblies from control and Uromyces pisi inoculated leafs of a susceptible and a partially rust-resistant grass pea genotype by RNA-seq.

Results

134,914 contigs, shared by both libraries, were used to analyse their differential expression in response to rust infection. Functional annotation grouped 60.4% of the contigs present in plant databases (37.8% of total) to 33 main functional categories, being “protein”, “RNA”, “signalling”, “transport” and “stress” the most represented. Transcription profiles revealed considerable differences in regulation of major phytohormone signalling pathways: whereas Salicylic and Abscisic Acid pathways were up-regulated in the resistant genotype, Jasmonate and Ethylene pathways were down-regulated in the susceptible one. As potential Resistance-genes we identified a mildew resistance locus O (MLO)-like gene, and MLO-related transcripts. Also, several pathogenesis-related genes were up-regulated in the resistant and exclusively down regulated in the susceptible genotype. Pathogen effectors identified in both inoculated libraries, as e.g. the rust Rtp1 transcript, may be responsible for the down-regulation of defence-related transcripts. The two genotypes contained 4,892 polymorphic contigs with SNPs unevenly distributed between different functional categories. Protein degradation (29.7%) and signalling receptor kinases (8.2%) were the most diverged, illustrating evolutionary adaptation of grass pea to the host/pathogens arms race.

Conclusions

The vast array of novel, resistance-related genomic information we present here provides a highly valuable resource for future smart breeding approaches in this hitherto under-researched, valuable legume crop.

Background

Rusts are among the most important diseases of legumes [1] and grass pea (Lathyrus sativus L.) is not an exception [2]-[4]. Rusts are caused by biotrophic fungi that keep infected host cells alive for their development. They form elaborate intracellular accommodation structures called haustoria, which maintain an intimate contact between fungal and plant cells over a prolonged period of time [5].

Rust in Lathyrus spp. is caused by Uromyces pisi (Pers.) Wint and U. viciae-fabae (Pers.) J. Schröt [6],[7], but and in addition to Lathyrus, U. pisi infects a broad range of other legumes too [7],[8]. Plants have developed multifaceted defence responses, many of which are induced only upon pathogen attack. These responses may include induction of pathogenesis related (PR) genes, the production of secondary metabolites (as e.g. phytoalexins), as well as the reinforcement of cell walls [9]. Associated with these responses may be the production of reactive oxygen species (ROS) and the induction of localized cell death (the hypersensitive response, HR) [10]. The induction of this basal plant defence machinery occurs upon the recognition of conserved molecules which are present in a variety of microbial species, but absent in the host. These pathogen associated molecular patterns (PAMPs) are molecular components highly conserved within a class of microbes, where they have essential functions for their fitness or survival [11]. These include, for example, fungal chitin, β-glucan and ergosterol. The specific virulence factors of the pathogen, known as fungal effectors, are recognized by corresponding resistance (R) genes of the host plant. Both rust-causing pathogens of Lathyrus are able to efficiently overcome R-gene based resistance [12]. To date, most fungal effectors identified are lineage-specific small secreted proteins (SSP) of unknown function [13],[14]. The U. viciae-fabae rust transferred protein 1 (Rtp1) was the first fungal effector visualized in the host cytoplasm and nucleus after in planta secretion by the rust fungus [15]. Rtp1 belongs to a family of cysteine protease inhibitors that are conserved in the rust species (order Pucciniales formerly known as Uredinales) [16].

Gene-for-gene resistance is associated with the activation of, for instance, the salicylic acid (SA)-dependent signalling pathway, leading to expression of defence-related genes like PR1, the production of ROS and finally to programmed cell death [17],[18]. Other phytohormones involved in plant/pathogen interaction are ethylene (ET) and jasmonates (JA). Plant defence responses appear specifically adapted to the attacking pathogen, with SA-dependent defences acting mainly against biotrophs, and JA- and ET-dependent responses acting mainly against necrotrophs [5],[19],[20].

Grass pea is a diploid species (2n = 14) with a genome size of approx. 8.2 Gbp [21]. Although grass pea is primarily self-pollinated, a 2 to 36% outcrossing rate was reported, depending on location and genotype [22]-[24]. Outcrossing is mainly driven by pollinators, and therefore can be minimized when grown in isolation [22]. There is a great potential for the expansion of grass pea in dry areas and zones that are becoming more drought-prone as a result of climate change [25]. Partial resistance to U. pisi has been reported in grass pea as a clear example of prehaustorial resistance, with no associated necrosis. This resistance is due to restriction of haustoria formation accompanied by frequent early abortion of the colonies, reduction in the number of haustoria per colony and decreased intercellular growth of infecting hyphae [26]. Though prehaustorial resistance is typical for non-hosts, it has also been implicated in host partial resistance [27],[28] and is common in resistance of major cool season grain legumes against rusts [1],[29]. Additionally, resistant Lathyrus genotypes may serve as a source of new and useful genetic traits in the breeding of related major legume crops such as peas, lentils and vetches. Cross-incompatibility has been reported between pea and L. sativus, but successful fusion of Pisum sativum and L. sativus protoplasts [30] creates new possibilities for gene transfer between these species. However, the slow progress in understanding the genetic control of important traits, such as disease resistance, in Lathyrus species hampered the development of modern cultivars or the introgression of their interesting traits into related species.

In economically important warm season legumes such as common bean and soybean, complete monogenically controlled resistances to rusts and associated rust resistance genes have been described together with closely linked markers for use in marker assisted backcrossing [29],[31]-[34]. By contrast, most rust resistances described so far in cool season food legumes are incomplete in nature and the genetic basis of resistance is largely unknown. Although QTL mapping studies confirmed the polygenic control of resistance as e.g. in pea [35], faba bean [36] and chickpea [37], no markers suitable for marker assisted selection (MAS) are available yet.

Genomic resources for grass pea are still scarce (e.g. in April 2014 the NCBI database contained only 178 EST sequences from L. sativus [38]), and the two linkage maps existing for grass pea do not contain sufficiently informative markers to bridge between them [4].

The advent of next-generation sequencing (NGS) technologies was an important breakthrough enabling the sensitive and quantitative high-throughput transcriptome analysis referred to as RNA-seq [39],[40]. RNA-seq discriminated between microbial and host transcriptomes, during plant-microbe interactions, using original or phylogenetically related genomes as a reference for transcript annotation [41]-[44]. RNA-seq gene expression patterns provided also information on complex regulatory networks and on variations in expressed genes, such as SNPs and SSRs, in an increasing number of non-model plants [45] and thus may be well suited to overcome the bottleneck of lacking genomic resources in Lathyrus.

Here we employed RNA-seq to study the response of L. sativus to U. pisi infection. We used MapMan and metabolic pathway analyses to interpret the results and assessed allelic diversity in transcripts as a source for genic markers for future (comparative) mapping studies. In addition, the expression of a set of selected genes was measured by qRT-PCR to validate the RNA-seq results.

To our knowledge, this is the first study on the global expression profiling of genes in grass pea/pathogen interaction using NGS. Our results will assist the elucidation of pathways and genes associated with resistance to rust in grass pea and related species. This approach may represent one of the initial steps towards the development of effective strategies for resistance breeding against such a quickly evolving pathogen.

Results

Contigs from the RNA-seq transcriptomes of resistant and susceptible L. sativusgenotypes

The RNA-seq libraries from control and inoculated leafs from the resistant genotype BGE015746 were united prior to assembly to generate a comprehensive data set enabling the generation of contigs of maximum length. They included 46,994,629 reads which were assembled into 105,288 contigs, ranging in size from 150 to 13,929 bp, with a mean contig length of 544 bp. The respective united library from the susceptible genotype BGE024709 comprised 72,566,465 reads which assembled in 119,870 contigs, with a size range of 150 to 15,658 bp and a mean contig length of 524 bp.

A reference assembly using both genotypes and treatments assembled in 134,914 contigs, ranging in size from 150 to 13,916 bp, with a mean contig length of 501 bp. The mapping and quantification of both genotypes’ libraries to the reference assembly allowed the analysis of their differential expression in response to U. pisi infection. 9,501 contigs were unique to the resistant and 15,645 contigs were unique to the susceptible genotype.

Redundancy of the reference assembly was checked using the clustering algorithm UCLUST, identifying only 49 (0.036%) transcripts with identity higher than 95%.

This Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GBSS00000000. The version described in this paper is the first version, GBSS01000000.

RNA-seq validation by quantitative RT-PCR assay

To validate the RNA-seq results, expression levels of a set of 9 selected genes were analysed by qRT-PCR. Genes were selected by their level of expression and transcript count, in order to represent a broad range of expression profiles. Further, the number of their transcripts differed between inoculated and control samples by log2 ratios ranging from −6.39 to 4.70 at q-values < 0.05. Their read count numbers were generally higher than 100, with exception of contig a45744;151, “mitochondrial chaperone BCS1”, with 2 counts in the resistant control and 36 counts in the resistant inoculated line, and contig a32859;123 “seed maturation protein”, with 3 counts in the susceptible inoculated line and 104 counts in the resistant inoculated line (Table 1). The best housekeeping genes for normalization suggested by the geNorm software were, for the resistant genotype samples, “β-tubulin” (a6507;507) and “photosystem I P700 apoprotein A2” (a160;902), and “O-methyltransferase” (a5102;390), for the susceptible genotype. A good correlation (R = 0.82 for the resistant and R = 0.80 for the susceptible genotypes) was observed between the log2 fold changes measured by RNA-seq and qRT-PCR (Figure 1).

Table 1 Log2 fold expression results for RNA-seq and qRT-PCR experiments
Full size table
Figure 1
figure 1

Correlation between RNA-seq and qRT-PCR. The relative expression levels obtained by RNA-seq using DEGseq and by qRT-PCR using the ΔΔCt method. Pearson’s correlation coefficient (R) between relative expression levels is shown above the trendline.

Full size image

Differential gene expression in resistant and susceptible L. sativusgenotypes during infection

Differentially expressed contigs were grouped by expression patterns based on up- or down-regulation (log2 ≥ 2 or log2 ≤ −2; respectively, q-value ≤ 0.05) after inoculation. Within each expression pattern group, comparisons were performed between genotypes. Expression patterns were grouped in eight response types, according to their up- or down-regulation, in susceptible and resistant genotypes, respectively. The number of contigs and description of each group is summarized in Table 2. Most representative groups are group F (contigs down regulated in both genotypes) and H (contigs down-regulated only in the resistant genotype) with 2,516 and 1,606 contigs respectively, followed by group A that includes 814 contigs up-regulated in both resistant and susceptible genotypes upon infection. A detailed list with all the identified contigs, their description and expression pattern groups can be found in Additional file 1.

Table 2 Classification of contigs according to their differential expression in the susceptible and resistant genotype upon infection with U. pisi
Full size table

As depicted in Figure 2, from the 134,914 contigs that could be identified and quantified, 68,889 were shared among all libraries. Of these, 974 contigs were up-regulated and 5,203 contigs down-regulated in the resistant genotype BGE015746 and 772 contigs up- and 4,617 down-regulated in the susceptible genotype BGE024709. Furthermore, from the 5,807 contigs only present in the resistant genotype’s libraries, 132 were up- and 485 down-regulated (inoculated vs. control). From the 7,938 contigs only found in the susceptible genotype’s libraries, 134 were up- and 689 down-regulated.

Figure 2
figure 2

Venn diagram of the number of unique and shared contigs between the two genotypes and its expression. In black boxes the number of up (log2 fold ≥ 2) and down (log2 fold ≤ −2) regulated contigs in the inoculated condition versus control. Resistant genotype: BGE015746, susceptible genotype: BGE024709.

Full size image

Annotation

From the 134,914 contigs detected in all libraries, 50,937 (37.75%) contigs could be matched via BLAST to entries in plant databases and 961 (0.71%) matched only to fungal databases. The latter contigs were present only in the inoculated libraries. Also, 4,558 contigs were absent in control samples and found exclusively in fungal databases, or with a higher bit-score in fungal databases than in plant databases and thus, most probably correspond to U. pisi sequences.

As indicated in Figure 3, BLAST produced hits mainly to other legume species with frequencies in the order Medicago truncatula (26,728; 19.81%), Glycine max (11,436; 8.48%), P. sativum (1,538; 1.14%) and Lotus japonicus (921; 0.68%). Vitis vinifera (2,409; 1.79%), Populus trichocarpa (656; 0.49%) and the model Arabidopsis thaliana (607; 0.45%) were the best matching non-legume species. BLAST hits from L. sativus comprised only 0.02% (33 contigs) of the total illustrating the scarcity of Lathyrus entries in the data bases.

Figure 3
figure 3

Number of contigs that could be BLASTed to different plant species.

Full size image

From the 4,558 contigs that were absent in control samples and found exclusively in fungal databases, or with a higher bit-score in fungal databases, 20 contigs from the 49 accessions described in UniProtKB/Swiss-Prot and UniProtKB/TrEMBL as U. viciae-fabae, were identified (see list in Additional file 2). None of these 20 contigs were significantly differentially expressed between the two inoculated genotypes. For example, among the eight contigs out of the 20 without a plant database hit, five were homologous to “invertase 1”, and the three others to “rust transferred protein – Rtp1”, “amino acid transporter” and “putative permease”. Six other contigs absent in control samples and found exclusively in fungal databases or with a higher bit-score in fungal databases, were homologous to housekeeping genes that can be found throughout different kingdoms (three “tubulin beta chain”, two “succinate dehydrogenase” and one “plasma membrane (H+) ATPase”.

Functional annotation of the contigs via Mercator and MapMan, depicted in Figure 4, grouped 60.4% of them into 33 main functional categories, of which the categories “protein” (11.0%), “RNA” (8.0%), “signalling” (6.7%), “transport” (5.4%) and “stress” (4.2%) were most crowded. A total of 39.4% could not be assigned to any functional category.

Figure 4
figure 4

Percentage of contigs assigned in each main functional category.

Full size image

Analysis of functional categories, within each expression pattern group, identified differences among the functions present within each group. Comparisons were also performed among the different expression profiles in each category (Figure 5). Transcripts included in the functional categories “stress” and “protein” were present at a higher percentage in up-regulated expression pattern groups (“stress” in A, B and C; “protein” in A, B and E), while the functional category “cell” was present at higher percentage in down-regulated expression pattern groups (F, G and H). The most prominent functional category in group C contigs up-regulated in both genotypes, with a higher expression in the resistant genotype was “cell wall”. “Lipid metabolism” and “DNA” were also over-represented. However, also the down-regulated groups F, G, and H contained a considerable number of contigs from the “cell wall” category. In group B, joining contigs up-regulated in both genotypes with a higher expression in the susceptible genotype, the categories “secondary metabolism” and “hormone metabolism” were over-represented. Interestingly, the functional category “signalling” was over-represented in contigs up- regulated only in the susceptible genotype, as in group D.

Figure 5
figure 5

Percentage of contigs assigned in each functional category for each expression pattern group. A - contigs up-regulated in both resistant and susceptible genotypes similarly; B - contigs up-regulated in both genotypes, with a higher expression in the susceptible genotype; C - contigs up-regulated in both genotypes, with a higher expression in the resistant genotype; D - contigs up-regulated only in the susceptible genotype; E - contigs up-regulated only in the resistant genotype; F - contigs down-regulated in both resistant and susceptible genotypes similarly; G - contigs down-regulated only in the susceptible genotype; H - contigs down-regulated only in the resistant genotype.

Full size image

Biotic stress related proteins

In order to restrict the number of analysed contigs to the ones probably more directly related to resistance, we focused mostly on contigs up-regulated at a higher ratio, or exclusively, in the resistant genotype (groups C and E), contigs exclusively down-regulated in the susceptible genotype (group G) and contigs exclusively down-regulated in the resistant genotype (group H).

From the subcategory “stress.biotic”, two contigs in group E corresponded to the well- studied mildew resistance locus O (MLO) gene which was first identified in barley, conferring resistance to powdery mildew [46]. Also, from a total of 25 “MLO-like” contigs, 12 were differentially expressed. One of these was down-regulated in the resistant genotype (group H). These might be related to MLO susceptibility genes, as reported by several previous studies [47]-[49]. Interestingly, in the susceptible genotype, one “PREDICTED: beta glucosidase 12-like”, identified by Mercator as “PENETRATION 2”, required for MLO-mediated resistance and belonging to the functional category “secondary metabolism”, was down-regulated (group G). Group G also contained one “acidic endochitinase” and two LRR proteins, one TIR-NBS-LRR and one containing LRR and NB-ARC domains. In group C, a pathogenesis related protein 1 (PR-1) contig was identified.

The subcategory “stress.abiotic” contained, i.a., genes involved in response to heat that also respond to biotic stresses. For example, in group C and E, we identified one “DNAJ heat shock protein” in both groups, three “heat shock protein 70 family” (group E) and one “18.1 kDa class 1 heat shock protein” (group E). Group G, however, contained one “DNAJ homolog subfamily B member” and one “double Clp-N motif-containing P-loop nucleoside triphosphate hydrolases superfamily protein”.

Several contigs related to secondary metabolism were exclusively up-regulated in the resistant genotype (group E). These comprised a “reticuline oxidase-like protein” involved in alkaloid biosynthesis, an “isoflavone 2’hydroxylase”, functioning in the isoflavonoid biosynthesis pathway, a “dihydroflavonol-4-reductase”, with roles in the flavonoid and brassinosteroid metabolic pathway and an “AMP-dependent CoA ligase”, acting in the JA and lignin biosynthesis pathways. In group G, 17 contigs were related to secondary metabolism including four involved in the flavonoid pathway, two in the isoprenoid/terpenoid pathway and one “WAX 2-like” involved in wax biosynthesis.

PTI (pathogen-associated molecular pattern triggered immunity) relies on an efficient signalling network in order to control the infection [50]. Receptor kinases are important for the plant’s pathogen recognition and their expression may be constitutively expressed or up-regulated in resistant genotypes or down-regulated in susceptible genotypes in response to effectors from the pathogen. Receptor kinases and kinases exclusively up-regulated in the resistant genotype and contained in group E may be part of such signalling cascades. These included one protein kinase with thaumatin (PR-5) domain, six “DUF 26”, one “CRINKLY4”, one “FERONIA receptor like kinase”, and also a MAP kinase “MAPKKK5” and a G-protein “zinc finger (Ran-binding) family protein”. In contrast, the down-regulation of such transcripts in the susceptible genotype (group G) may contribute to susceptibility. Here we identified three “DUF 26”, three LRR (“NIK1”, “RKF1” and “PXY”), two G-proteins (“guanine nucleotide-binding protein” and “dynamin-related protein 1E-like”), two MAP kinases (“PAS domain-containing protein tyrosine kinase family protein”) and three genes involved in calcium signalling (“calcium-transporting ATPase”, “calmodulin-binding heat-shock protein” and “calmodulin-domain protein kinase 9”). Interestingly, calmodulin also plays a role in the MLO response, where the lack of a calmodulin binding site decreases its defence response [51].

The “cell wall” category contained seven cellulose synthase contigs: one in group E “IRREGULAR XYLEM 3 (IRX3)” and four in group C (three “IRX1” and one “CESA1”). In group G, we identified two cellulose synthase “IRX14” and two “pectinesterase inhibitor” contigs.

From the genes normally associated with defence response, only one “endo-beta-1 3-glucanase” was identified in group C, while two others “endo-beta-1 3-glucanase” were detected in group G. Also in group G, we identified two “peroxidase” and two “glutathione S-transferase” genes.

SNPs in resistance pathways

In the 68,889 contigs present in both the susceptible and the resistant genotypes, we identified 2,634 contigs containing Single Nucleotide Polymorphisms (SNPs) discriminating between their respective alleles. The number of SNPs in functional (MapMan) categories varied considerably. The categories “RNA regulation of transcription” (9.5%) and “protein.degradation” (8.9%) contained by far the most SNPs, followed by the protein-related categories “protein.postranslational modification” (4.3%) and “protein.synthesis” (3.2%). Other categories including the most SNP-containing contigs were “signalling.receptor kinases” (2.5%), “protein.targeting” (2.4%) and the stress related categories, “stress.biotic” (1.8%) and “stress.abiotic” (1.6%) (Figure 6).

Figure 6
figure 6

Percentage of contigs containing SNPs between the resistant and susceptible genotypes in each Mercator mapping functional sub-category. FA: fatty acid; metabol: metabolism; misc.: miscellaneous; PS: photosynthesis; red.: redox.

Full size image

EST-SRR development

200 EST-SSR potential polymorphic markers between the two genotypes were designed. EST-SSRs were identified by the Phobos software [52], using as search parameters, perfect SSRs with a repeat unit lenght of two to six nucleotides. Polymorphisms between the resistant and susceptible genotypes were manually identified and flanked by primer pair using the Primer3 software [53]. To validate the EST-SSR sequences, 40 primer pairs were randomly selected for PCR amplification to confirm the presence of size polymorphism between the two accessions. PCR reactions were conducted twice in order to confirm the results. From the total 40 EST-SSR tested, 25 (62.5%) primer pairs successfully amplified polymorphic fragments between the two accessions. 6 (15.0%) primer pairs amplified monomorphic fragments and 5 (12.5%) produce a very complex pattern. The remaining 4 (10.0%) primer pairs were not able to produce any fragments.

Discussion

Lathyrus spp. is a potential source of resistance to several pathogens [4],[25] and especially L. sativus provides resistance to several fungal and bacterial diseases [26],[38],[54],[55]. However, the lack of genetic and/or genomic information was a barrier to further identify resistance-related genes and to use them in breeding.

In the present study we therefore attempted to improve this unfavourable situation by identifying ESTs and SNPs, potentially involved in resistance, that may be used in future smart breeding approaches. We describe for the first time a high-throughput transcriptome assembly of grass pea/pathogen interaction, using genotypes contrasting in response to rust infection, to unravel the involved partial resistance mechanism and associated resistant genes.

Our study has identified a large number of differentially expressed genes corresponding to biological categories that are thought to be most relevant in grass pea response to rust. A limitation of our study is the fact that only a single pooled sample was investigated for each genotype and condition. Although the biological variance could not be assessed in the bulked approach, the large number of individual samples in the pool is likely to level out many of possible outliers. Nevertheless, the validation of twelve genes by RT-qPCR, using three biological replicates, provided a good correlation with RNA-seq results.

Another motive that could also be influencing our results is that we used different cDNA synthesis primers, oligo(dT) for the RNA-seq and poly(A) for RT-qPCR, what might yield different quantities of poly-adenylated and non-adenylated transcripts.

Our study was severely hampered by the low number of annotated sequences, which is due to the lack of a reference genomic sequence for Lathyrus. Nevertheless, we could annotate between 34% and 46% of differentially expressed contigs to hits in plant databases, depending on the genotype and the infection status of the plants. We further developed new gene-based molecular tools as e.g. expressed sequence tags, gene-based simple sequence repeats (EST-SSR) and SNP-based markers. Moreover, we identified a number of U. pisi effectors in the infected tissues though the overall low number of observed fungal transcripts probably reflects the low quantity of fungal structures in early-infected leaves [56]. Thus, our present study will help to overcome the problems we encountered in previous work, where the transfer of molecular markers from close related species had a very low rate of success (18% for pea EST-SSRs and 6% for pea genomic SSRs, [57]) Therefore, the present RNA-seq libraries will boost the availability of specific EST-SSRs and SNP-based markers that will be equally important for future development of more effective grass pea resistance breeding approaches.

The high amplification rate of the developed EST-SSRs validates the quality of the RNA-seq data. The few primers that failed to produce amplification products or produced amplicons with an unexpected pattern may be caused by the location of the respective primers across splice regions or the presence of a large intron, since genomic regions are absent from cDNA. In addition also primers could be derived from chimeric cDNA clones [58].

Besides the novel markers, the deep insights into pathogenesis-related mechanisms provided by this study are of particular interest. The most interesting pathogenesis-related protein that we identified, the “MLO-like protein” is involved in signalling in response to biotic stress. MLO was described for the first time in barley, where it conferred partial resistance to powdery mildew by inducing the thickening of the cell wall at fungal penetration sites [46]. Two “MLO-like” contigs were up-regulated exclusively in the resistant genotype (group E), and perhaps related to this, also in group E and more strongly up-regulated in the resistant genotype (group C), we identified cellulose biosynthesis genes. The exclusively resistance-up-regulated group E contained one “IRREGULAR XYLEM 3” (IRX3) gene and three “IRX1”. Additionally, one “cellulose synthase 1” (CESA1) was stronger up-regulated in the resistant genotype than in the susceptible one (group C). Consistent with the assumed importance of MLO signalling for rust resistance, some genes important for MLO function as e.g. “calmodulin”, involved in calcium signalling as a prerequisite for MLO function [51], were down-regulated in the susceptible genotype (group G). Another gene involved in MLO resistance and down-regulated in the susceptible genotype, is the glycosyl hydrolase “PENETRATION 2” gene [59]. Since rust resistance in grass pea is of prehaustorial type we consider MLO as a candidate R-gene. In order to confirm this assumption, callose deposition, as a potentially durable resistance mechanism against rusts, should be further investigated in rust-resistant and susceptible grass pea genotypes. The down-regulation of several “MLO-like” contigs in response to infection in both, susceptible and resistant genotypes, does not necessarily contradict our assumption, since several MLO orthologs were demonstrated to function as susceptibility genes [47]-[49].

Plant responses to biotic stressors are, i.a., controlled by phytohormones as e.g. salicylic acid (SA), abscisic acid (ABA), jasmonates (JA) and ethylene (ET). Differences in expression of hormone-related genes of the susceptible and resistant genotype, in response to the pathogen, also occurred in our gene expression patterns. For example, plant resistance to biotrophic pathogens is mainly controlled by the SA pathway [18] and the importance of SA in the induction of systemic acquired resistance in legumes against rust fungi has been reported [60],[61]. In our study an inducer of the SA pathway, the “ethylene response factor 5” (ERF5) gene, which at the same time inhibits the JA and ET biosynthesis pathways [62], was exclusively up-regulated in the resistance genotype (group E), whereas two Apetala2/Ethylene Responsive Factor (AP2/ERF) transcription factor genes, important for the regulation of defence responses [63], were down- regulated in the susceptible genotype (group G).

ABA regulates defence responses through its effects on callose deposition and production of ROS intermediates [18], activating also stomata closure as a barrier against pathogen infection [64]. In the resistant genotype, the transcript for “9-cis-epoxycarotenoid dioxygenase 2”, a key regulator of ABA biosynthesis in response to drought [65], and involved in the crosstalk between ABA and SA signalling in plant-pathogen interactions [66], was up-regulated (group E), whereas several transcripts engaged in ABA, auxin and JA signalling, were down-regulated in the susceptible genotype (group G). This is consistent with a susceptible response to a biotroph attack [18].

The image emerging from the transcription profiles, of the resistant and susceptible genotype, further highlights that pathogenesis related (PR) proteins are key players in Lathyrus-rust interactions since several PR genes were mainly up-regulated in the resistant genotype, after inoculation. Among these were two chitinases (PR-3 and PR-9) involved in the degradation of the fungal cell wall [9] and a thaumatin (PR-5) gene, which causes an increase of the permeability of fungal membranes by pore-forming mechanisms [67]. In group E, we found a “pathogenesis related protein 1” (PR-1) and a “protein kinase-coding resistance protein”, a receptor kinase with a thaumatin domain (PR5K), presumably involved in thaumatin signaling and described previously as delaying infection [68]. Another important PR-gene, an “acidic endochitinase” (PR3), was down regulated exclusively in the susceptible genotype. Genes involved in secondary metabolism were also detected. Legumes utilize flavonoids, notably isoflavones and isoflavanones, for defence against pathogens and as signalling molecules, with a number of phenylpropanoids having antimicrobial activity and restricting pathogen growth and disease symptoms [69]. In group G, we identified a “reticuline oxidase-like protein”, up regulated in non- race-specific resistance to stripe rust in wheat [70], an “isoflavone 2′-hydroxylase” from the isoflavonoid pathway [71] and a “dihydroflavonol-4-reductase” catalysing the first enzymatic step in anthocyanin biosynthesis, in the flavonoid pathway [72].

Also, exclusively down regulated in the susceptible genotype, we found some genes important for defence response within the miscellaneous category, like “endo-beta-1 3-glucanase”, “glutathione S-transferase” (GST) and “peroxidase”. In plants, beta-glycosidases serve a number of diverse and important functions, including bioactivation of defence compounds, cell wall degradation in endosperm during seed germination, activation of phytohormones, and lignifications [73]. GSTs are detoxification-related proteins, protecting cells from secondary metabolites produced in response to pathogen attack, including peroxidases [74]. Finally, peroxidases function as resistance factors against invading fungi, inhibiting hyphal elongation, and when H2O2 is present, causing oxidative burst [75].

Effectors are expected to be excellent targets for the control of pathogens, but, unlike effectors from some other plant pathogens, relatively little is known about rust effectors [76]. In this study, several unigenes were identified in fungal databases. The most known rust effector identified was “rust transferred protein 1” (Rtp1). This effector aggregates into amyloid-like filaments in vitro [77]. Immunoelectron microscopy localized this effector to the extrahaustorial matrix protuberances extending into the host cytoplasm, although the exact role for this protein remains to be discovered [78]. Other Uromyces effectors identified in this study were “succinate dehydrogenase”, “invertase” and “permease”. From the total potential rust transcripts identified, a selection of effector proteins could be used as probes to identify the target host proteins as a first step in the development of effector-driven legume breeding, maximizing the durability of resistance against the quickly evolving rust pathogens [79].

From the DE contigs obtained in the present study, 2,634 presented SNPs between the resistant and the susceptible lines. The MapMan software aided in the functional categorisation of SNPs, revealing that the categories “RNA regulation of transcription” (9.5%) and “protein.degradation” (8.9%) contained by far the most SNPs. Within these categories, ubiquitins were most polymorphic (5.6%). Ubiquitins tag proteins for proteasome degradation and play a central role in signalling pathways [80]. Especially ubiquitin “E3 RING” and “SCF F-BOX” contigs contained a large number of SNPs (1.7% and 1.8% respectively). E3 RING and SCF F-BOX proteins are involved in several aspects of plant immunity ranging from pathogen recognition to both PTI to effector-triggered immunity (ETI). From the differentially expressed contigs identified as containing SNPs, we found one E3 RING, “PREDICTED: RING-H2 finger protein ATL2-like”, down regulated in the susceptible genotype. Four other functional categories, “RNA.regulation of transcription” (9.5%), “protein.postranslational modification” (4.3%), “protein.synthesis” (3.2%) and “signalling.receptor kinases” (2.5%), also contained significant numbers of polymorphisms. Especially the “signalling.receptor kinases” category may be of particular interest for further studies since receptor kinases recognize pathogen effectors and their rapid evolution, reflected by large numbers of polymorphisms, may represent plants adaptation to a rapidly changing spectrum of pathogens in the arms race between them and their hosts [81].

The large number of SNPs that we identified will be instrumental for the development of linkage and high-throughput association mapping approaches and for the expansion of our previous diversity studies in Lathyrus [57],[82].

Our results provide an overview of gene expression profiles of contrasting L. sativus genotypes inoculated with rust, offering a valuable set of sequence data for candidate rust resistant gene discovery.

Conclusions

Our transcriptome analysis provided comprehensive insight into the molecular mechanisms underlying prehaustorial rust resistance in L. sativus.

The differences in resistance between the two L. sativus genotypes investigated appear to be mainly due to the activation of the SA pathway and several pathogenesis related genes, including the ones regulated by MLO. The fastest-evolving pathways differentiating between the two genotypes are the general RNA’s regulation of transcription, followed by the Ubiquitin-26S proteasome system and having also as most mutated receptor-based signalling genes and biotic and abiotic stress related genes. The detected polymorphic SNPs will allow the development of new gene-based molecular tools. Altogether, 51 genes were identified as potential resistance genes, prioritizing them as specific targets for future functional studies on grass pea/rust interactions. Besides a plethora of pathogenesis-related host genes, 4,558 transcripts, including putative effectors, were also identified for the rust fungus U. pisi. As a consequence of the newly developed wider array of genetic and genomic resources, future work will focus on high throughput mapping of the genetic basis of disease resistance in L. sativus and eventual comparative mapping with other legume species, contributing all to an improved exploitation of this under used highly potential legume species.

Methods

Plant and fungal material, inoculation and RNA isolation

The two L. sativus genotypes, BGE015746 and BGE024709 analysed in the present work were kindly provided by the Plant Genetic Resources Centre (CRF-INIA), Madrid, Spain. Seeds were multiplied in insect proof cages in order to minimize outcrossing. Evaluation for their resistance against U. pisi demonstrated that BGE024709 is susceptible to rust, whereas BGE015746 displays partial resistance [26]. Upon infection, both genotypes present well-formed pustules, with no associated chlorosis or necrosis. They contrast, however, in disease severity (DS), i.e. the percentage of leaf area covered by the fungus. Whereas the partial resistant genotype has a DS = 9%, the susceptible one has a DS = 30%.

The U. pisi monosporic isolate UpCo-01 from the fungal collection of the Institute for Sustainable Agriculture-CSIC (Córdoba, Spain) was used for the experiment. Inoculum was multiplied on plants of the susceptible P. sativum cv. Messire before use.

Twenty-four plants per each genotype and treatment (inoculated/control) were used. Two-week-old L. sativus seedlings were inoculated by dusting all the plants at the same time with 2 mg of spores per plant, diluted in pure talk (1:10), with the help of a small manual dusting device in a complete random experiment. Inoculated and control plants were incubated for 24 h at 20°C, in complete darkness, and 100% relative humidity, then transferred to a growth chamber and kept at 20 ± 2°C under 14 h light (150 μmol m−2 s−1) and 10 h dark.

RNA was extracted from inoculated and non-inoculated fresh leaves collected 37 hours after inoculation. The material was immediately frozen in liquid nitrogen and stored at −80°C. RNA was isolated using the GeneJET Plant RNA Purification Mini Kit (Thermo Scientific, Vilnius, Lithuania), according to the manufacturer’s instructions. The extracted RNA was treated with Turbo DNase I (Ambion, Austin, TX, USA), and RNA quantification was carried out using the NanoDrop device (Thermo Scientific, Passau, Germany).

Sequencing and quantification

For each of the 4 combinations, genotype and treatment (BGE015746 control, BGE015746 inoculated, BGE024709 control and BGE024709 inoculated) total RNA from 24 plants was extracted and pooled in equal amounts for sequencing. Three RNA-seq libraries (one for each genotype and one reference assembly, including all genotypes and treatments) were generated by GenXPro GmbH, Germany, using a proprietary protocol. In short, mRNA was captured from 20 μg of total RNA using Oligo dT(25) beads (Dynabeads; life Technologies). The purified mRNA was randomly fragmented in a Zn2+ solution to obtain approximately 250 bp long RNA fragments. cDNA was synthesized by reverse transcription starting from 6(N) random hexamer oligonucleotides followed by second strand synthesis. Barcoded Y-adapters were ligated to the cDNA and the library was amplified with 10 cycles of PCR. The libraries were sequenced on an Illumina Hiseq2000 machine. After Illumina paired-end sequencing, raw sequence reads were passed through quality filtering, thereby also removing sequencing adapter primers and cDNA synthesis primers. All high-quality reads were assembled using the Trinity RNA-Seq de novo assembly (Version: trinityrnaseq_r2011-11-26). In order to minimize the redundancy, CAP3 software [83] was also used with overlap length cutoff of 30 bp and overlap percent identity cutoff of 75%. Redundancy was tested using the clustering algorithm UCLUST ([84], available at http://drive5.com/usearch/manual/uclust_algo.html). The resulting contigs were annotated via BLASTX to publically available databases (ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz, nr, plants only). To identify fungal transcripts, an additional BLASTX to public fungal databases (http://www.ebi.ac.uk/uniprot, UniProtKB/Swiss-Prot and UniProtKB/TrEMBL) was performed. The sequenced reads were mapped with novoalign software (V2.07.14; http://www.novocraft.com/) to the own assembled contigs. RPKM was calculated as the normalized transcript expression value [85]. Our obtained counts were subsequently passed through DEGSeq to calculate the differential gene expression (R package version 1.16.0) [86].

SNP detection

SNPs were discovered between the two L. sativus genotypes using Joint-SNV-Mix [87]. The mappings from the transcriptome analysis were also analysed by Joint-SNV-Mix and the output was furthermore processed by GenXPro’s in-house software to detect SNPs discriminating the bulks. SNP calling was performed taking in account only the inoculated samples. A minimum coverage of 15 reads in each genotype in the inoculated condition was needed to call a SNP. Polymorphic contigs and their respective SNPs are listed in Additional file 3.

EST-SSR development and genotyping

EST-SSRs were selected in silico by identifying the polymorphic SSRs between the two L. sativus genotypes. Identification of the SSRs was done using Phobos plug-in [52] for the Geneious software [88], using as search parameters, perfect SSRs with a repeat unit length of two to six nucleotides. Length polymorphisms were manually identified by aligning SSR-containing contigs of one genotype against the whole library of the other genotype. Primers were designed using Primer3 plug-in [53] for the Geneious software, using as parameters a melting temperature from 59 to 63°C, a GC content of 50 to 60% and a primer size raging from 18 to 24 nucleotides. The developed EST-SSR markers are listed in Additional file 4.

PCR reactions for the EST-SSRs genotyping were conducted using the M13 tail labelling strategy described by Schuelke [89] in a total volume of 10 μl containing 10 ng of template DNA, 0.04 μM of M13(−21) tagged forward primer, 0.16 μM of IRD700 or IRD 800 M13(−21) and 0.16 μM of reverse primer, 0.2 mM of each dNTP, 1.5 mM of MgCl2, and 0.2 unit of Taq DNA polymerase (Promega, Madison, USA). The amplification reaction consisted of a denaturing step of 5 min at 94°C, followed by 30 cycles of 30 s at 94°C, 45 s at 56°C, 45 s at 72°C, and 8 cycles of 30 s at 94°C, 45 s at 53°C, 45 s at 72°C. The reaction was terminated at 72°C for 10 min.

SSR fragments were resolved with 6.5% polyacrylamide gel using a LI-COR 4300 DNA Analyzer (Lincoln, NE, USA).

Quantitative RT-PCR assay

1 μg of total RNA from each of three randomly chosen plants per genotype, per treatment (inoculated/control), was reverse transcribed in duplicates, using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA) following manufacturer’s instructions. Three independent reverse-transcription reactions (RT) were performed for each cDNA sample in a total of nine samples per genotype, per treatment. For all genes studied, the product of each of these reactions was analysed in technical duplicates, in a total of six technical replicates per treatment. qRT-PCR reactions were performed with a iQ™5 Real-Time PCR Detection System (Bio-Rad, Munich, Germany). Primers were designed using the Primer3 software [53]. Primer sequences can be found in Additional file 5. For data analysis, the Genex software package (MultiD, Goteborg, Sweden), including the geNorm software [90] was used.

Contig annotation and data analysis

In order to classify the contigs into functional categories, the Mercator pipeline for automated sequence annotation ([91], available at http://mapman.gabipd.org/web/guest/app/mercator) was employed. The mapping file was created excluding contigs without BLAST hit in previous analyses and accessing the following, manually curated databases: Arabidopsis TAIR proteins (release 10), SwissProt/UniProt Plant Proteins (PPAP), TIGR5 rice proteins (ORYZA), Clusters of orthologous eukaryotic genes database (KOG), Conserved domain database (CDD) and InterPro scan (IPR). The Mercator mapping file was then employed for pathway analysis by the MapMan software ([92], available at http://mapman.gabipd.org/web/guest/mapman).

Differentially expressed contigs were identified by comparing their expression in leaves of the resistant genotype BGE015746, control vs. inoculated, and of the susceptible BGE024709, control vs. inoculated, using DEGseq [86]. In cases where a particular transcript reacted in the same way in both genotypes, the total transcript count before and after inoculation was compared, allowing the identification of basal genotypic differences between the two genotypes.

Availability of supporting data

The raw RNA-seq data supporting the result of this article is available in the Sequence Read Archive (SRA), with accession numbers SRS686331, SRS687370, SRS687371 and SRS687373. This Transcriptome Shotgun Assembly (TSA) project has been deposited at DDBJ/EMBL/GenBank under the accession GBSS00000000. The version described in this paper is the first version, GBSS01000000.

Additional files

References

  1. Sillero JC, Fondevilla S, Davidson J, Patto MCV, Warkentin TD, Thomas J, Rubiales D: Screening techniques and sources of resistance to rusts and mildews in grain legumes. Euphytica. 2006, 147 (1–2): 255-272. 10.1007/s10681-006-6544-1.

    Article  Google Scholar 

  2. Duke J: Handbook of Legumes of World Economic Importance. New York,USA: Plennun Press; 1981.

    Book  Google Scholar 

  3. Campbell CG: Promoting the conservation and use of underutilized andneglected crops. Grass pea. Lathyrus sativus L., vol. 18. Rome, Italy:International Plant Genetic Resources Institute; 1997.

    Google Scholar 

  4. Vaz Patto MC, Skiba B, Pang ECK, Ochatt SJ, Lambein F, Rubiales D: Lathyrus improvement for resistance against biotic and abiotic stresses: From classical breeding to marker assisted selection. Euphytica. 2006, 147 (1–2): 133-147. 10.1007/s10681-006-3607-2.

    Article  Google Scholar 

  5. O’Connell RJ, Panstruga R: Tête à tête inside a plant cell: establishing compatibility between plants and biotrophic fungi and oomycetes. New Phytol. 2006, 171 (4): 699-718. 10.1111/j.1469-8137.2006.01829.x.

    Article  PubMed  Google Scholar 

  6. Barilli E, Satovic Z, Sillero JC, Rubiales D, Torres AM: Phylogenetic analysis of Uromyces species infecting grain and forage legumes by sequence analysis of nuclear ribosomal internal transcribed spacer region. J Phytopathol. 2011, 159 (3): 137-145. 10.1111/j.1439-0434.2010.01736.x.

    Article  CAS  Google Scholar 

  7. Barilli E, Moral A, Sillero JC, Rubiales D: Clarification on rust species potentially infecting pea (Pisum sativum L.) crop and host range of Uromyces pisi (Pers.) Wint. Crop Prot. 2012, 37: 65-70. 10.1016/j.cropro.2012.01.019.

    Article  Google Scholar 

  8. Rubiales D, Sillero JC, Emeran AA: Response of vetches (Vicia spp.) to specialized forms of Uromyces vicia-fabae and to Uromyces pisi . Crop Prot. 2013, 46: 38-43. 10.1016/j.cropro.2012.12.011.

    Article  Google Scholar 

  9. Stintzi A, Heitz T, Prasad V, Wiedemann-Merdinoglu S, Kauffmann S, Geoffroy P, Legrand M, Fritig B: Plant ‘pathogenesis-related’ proteins and their role in defense against pathogens. Biochimie. 1993, 75 (8): 687-706. 10.1016/0300-9084(93)90100-7.

    Article  CAS  PubMed  Google Scholar 

  10. Zurbriggen MD, Carrillo N, Hajirezaei MR: ROS signaling in the hypersensitive response: when, where and what for?. Plant Signal Behav. 2010, 5 (4): 393-396. 10.4161/psb.5.4.10793.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Medzhitov R, Janeway CA: Innate immunity: the virtues of a nonclonal system of recognition. Cell. 1997, 91 (3): 295-298. 10.1016/S0092-8674(00)80412-2.

    Article  CAS  PubMed  Google Scholar 

  12. McDonald B, Linde C: The population genetics of plant pathogens and breeding strategies for durable resistance. Euphytica. 2002, 124 (2): 163-180. 10.1023/A:1015678432355.

    Article  CAS  Google Scholar 

  13. Stergiopoulos I, de Wit PJGM: Fungal effector proteins. Annu Rev Phytopathol. 2009, 47 (1): 233-263. 10.1146/annurev.phyto.112408.132637.

    Article  CAS  PubMed  Google Scholar 

  14. Schmidt SM, Panstruga R: Pathogenomics of fungal plant parasites: what have we learnt about pathogenesis?. Curr Opin Plant Biol. 2011, 14 (4): 392-399. 10.1016/j.pbi.2011.03.006.

    Article  CAS  PubMed  Google Scholar 

  15. Kemen E, Kemen AC, Rafiqi M, Hempel U, Mendgen K, Hahn M, Voegele RT: Identification of a protein from rust fungi transferred from haustoria into infected plant cells. Mol Plant-Microbe Interact. 2005, 18 (11): 1130-1139. 10.1094/MPMI-18-1130.

    Article  CAS  PubMed  Google Scholar 

  16. Pretsch K, Kemen A, Kemen E, Geiger M, Mendgen K, Voegele R: The rust transferred proteins-a new family of effector proteins exhibiting protease inhibitor function. Mol Plant Pathol. 2013, 14 (1): 96-107. 10.1111/j.1364-3703.2012.00832.x.

    Article  CAS  PubMed  Google Scholar 

  17. Pieterse CMJ, Van der Does D, Zamioudis C, Leon-Reyes A, Van Wees SCM: Hormonal modulation of plant immunity. Annu Rev Cell Dev Biol. 2012, 28: 489-521. 10.1146/annurev-cellbio-092910-154055.

    Article  CAS  PubMed  Google Scholar 

  18. Bari R, Jones JG: Role of plant hormones in plant defence responses. Plant Mol Biol. 2009, 69 (4): 473-488. 10.1007/s11103-008-9435-0.

    Article  CAS  PubMed  Google Scholar 

  19. Glazebrook J: Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu Rev Phytopathol. 2005, 43: 205-227. 10.1146/annurev.phyto.43.040204.135923.

    Article  CAS  PubMed  Google Scholar 

  20. Smith JL, De Moraes CM, Mescher MC: Jasmonate- and salicylate-mediated plant defense responses to insect herbivores, pathogens and parasitic plants. Pest Manag Sci. 2009, 65 (5): 497-503. 10.1002/ps.1714.

    Article  CAS  PubMed  Google Scholar 

  21. Bennett MD, Leitch IJ, et al: Angiosperm DNA C-values database (release 8.0, Dec. 2012). 2012 [http://www.kew.org/cvalues/]

    Google Scholar 

  22. Chowdhury MA, Slinkard AE: Natural outcrossing in grasspea. J Hered. 1997, 88 (2): 154-156. 10.1093/oxfordjournals.jhered.a023076.

    Article  Google Scholar 

  23. Rahman MM, Kumar J, Rahman MA, Ali Afzal M: Natural outcrossing in Lathyrus sativus L. Indian J Genet Plant Breed. 1995, 55 (2): 204-207.

    Google Scholar 

  24. Gutiérrez-Marcos JF, Vaquero F, de Miera LE S, Vences FJ: High genetic diversity in a world-wide collection of Lathyrus sativus L. revealed by isozymatic analysis. Plant Genet Res. 2006, 4 (03): 159-171. 10.1079/PGR2006115.

    Article  Google Scholar 

  25. Hillocks R, Maruthi M: Grass pea (Lathyrus sativus): Is there a case for further crop improvement?. Euphytica. 2012, 186 (3): 647-654. 10.1007/s10681-012-0702-4.

    Article  Google Scholar 

  26. Vaz Patto MC, Rubiales D: Identification and characterization of partial resistance to rust in a germplasm collection of Lathyrus sativus L. Plant Breed. 2009, 128 (5): 495-500. 10.1111/j.1439-0523.2008.01601.x.

    Article  Google Scholar 

  27. Rubiales D, Niks RE: Characterization of Lr34, a major gene conferring nonhypersensitive resistance to wheat leaf rust. Plant Dis. 1995, 79 (12): 1208-1212. 10.1094/PD-79-1208.

    Article  Google Scholar 

  28. Sillero JC, Rubiales D: Histological characterization of resistance to Uromyces viciae-fabae in faba bean. Phytopathology. 2002, 92 (3): 294-299. 10.1094/PHYTO.2002.92.3.294.

    Article  CAS  PubMed  Google Scholar 

  29. Rubiales D, Castillejo MA, Madrid E, Barilli E, Rispail N: Legume breeding for rust resistance: lessons to learn from the model Medicago truncatula . Euphytica. 2011, 180 (1): 89-98. 10.1007/s10681-011-0367-4.

    Article  Google Scholar 

  30. Durieu P, Ochatt SJ: Efficient intergeneric fusion of pea (Pisum sativum L.) and grass pea (Lathyrus sativus L.) protoplasts. J Exp Bot. 2000, 51 (348): 1237-1242. 10.1093/jexbot/51.348.1237.

    Article  CAS  PubMed  Google Scholar 

  31. Faleiro F, Ragagnin VA, Moreira MA, de Barros EG: Use of molecular markers to accelerate the breeding of common bean lines resistant to rust and anthracnose. Euphytica. 2004, 138 (3): 213-218. 10.1023/B:EUPH.0000047080.80405.72.

    Article  CAS  Google Scholar 

  32. Miklas P, Kelly J, Beebe S, Blair M: Common bean breeding for resistance against biotic and abiotic stresses: From classical to MAS breeding. Euphytica. 2006, 147 (1–2): 105-131. 10.1007/s10681-006-4600-5.

    Article  CAS  Google Scholar 

  33. Hyten DL, Hartman GL, Nelson RL, Frederick RD, Concibido VC, Narvel JM, Cregan PB: Map location of the Rpp1 locus that confers resistance to soybean rust in soybean. Crop Sci. 2007, 47 (2): 837-840. 10.2135/cropsci2006.07.0484.

    Article  CAS  Google Scholar 

  34. Garcia A, Calvo ES, de Souza Kiihl RA, Harada A, Hiromoto DM, Vieira LG: Molecular mapping of soybean rust (Phakopsora pachyrhizi) resistance genes: discovery of a novel locus and alleles. Theor Appl Genet. 2008, 117 (4): 545-553. 10.1007/s00122-008-0798-z.

    Article  CAS  PubMed  Google Scholar 

  35. Barilli E, Satovic Z, Rubiales D, Torres AM: Mapping of quantitative trait loci controlling partial resistance against rust incited by Uromyces pisi (Pers.) Wint. in a Pisum fulvum L. intraspecific cross. Euphytica. 2010, 175 (2): 151-159. 10.1007/s10681-010-0141-z.

    Article  CAS  Google Scholar 

  36. Rai R, Singh AK, Singh BD, Joshi AK, Chand R, Srivastava CP: Molecular mapping for resistance to pea rust caused by Uromyces fabae (Pers.) de-Bary. Theor Appl Genet. 2011, 123 (5): 803-813. 10.1007/s00122-011-1628-2.

    Article  PubMed  Google Scholar 

  37. Madrid E, Rubiales D, Moral A, Moreno MT, Millán T, Gil J, Rubio J: Mechanism and molecular markers associated with rust resistance in a chickpea interspecific cross (Cicer arietinum × Cicer reticulatum). Eur J Plant Pathol. 2008, 121 (1): 43-53. 10.1007/s10658-007-9240-7.

    Article  CAS  Google Scholar 

  38. Skiba B, Ford R, Pang ECK: Construction of a cDNA library of Lathyrus sativus inoculated with Mycosphaerella pinodes and the expression of potential defence-related expressed sequence tags (ESTs). Physiol Mol Plant Pathol. 2005, 66 (1–2): 55-67. 10.1016/j.pmpp.2005.04.005.

    Article  CAS  Google Scholar 

  39. Simon SA, Zhai J, Nandety RS, McCormick KP, Zeng J, Mejia D, Meyers BC: Short-read sequencing technologies for transcriptional analyses. Annu Rev Plant Biol. 2009, 60: 305-333. 10.1146/annurev.arplant.043008.092032.

    Article  CAS  PubMed  Google Scholar 

  40. Metzker ML: Sequencing technologies - the next generation. Nat Rev Genet. 2010, 11 (1): 31-46. 10.1038/nrg2626.

    Article  CAS  PubMed  Google Scholar 

  41. Kemen E, Gardiner A, Schultz-Larsen T, Kemen AC, Balmuth AL, Robert-Seilaniantz A, Bailey K, Holub E, Studholme DJ, Maclean D, Jones JD: Gene gain and loss during evolution of obligate parasitism in the white rust pathogen of Arabidopsis thaliana . PLoS Biol. 2011, 9 (7): e1001094-10.1371/journal.pbio.1001094.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Fernandez D, Tisserant E, Talhinhas P, Azinheira H, Vieira ANA, Petitot A-S, Loureiro A, Poulain J, Da Silva C, Silva MDOC, Duplessis S: 454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemileia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant–rust interaction. Mol Plant Pathol. 2012, 13 (1): 17-37. 10.1111/j.1364-3703.2011.00723.x.

    Article  CAS  PubMed  Google Scholar 

  43. Tisserant E, Kohler A, Dozolme-Seddas P, Balestrini R, Benabdellah K, Colard A, Croll D, Da Silva C, Gomez SK, Koul R, Ferrol N, Fiorilli V, Formey D, Franken P, Helber N, Hijri M, Lanfranco L, Lindquist E, Liu Y, Malbreil M, Morin E, Poulain J, Shapiro H, van Tuinen D, Waschke A, Azcón-Aguilar C, Bécard G, Bonfante P, Harrison MJ, Küster H, et al: The transcriptome of the arbuscular mycorrhizal fungus Glomus intraradices (DAOM 197198) reveals functional tradeoffs in an obligate symbiont. New Phytol. 2012, 193 (3): 755-769. 10.1111/j.1469-8137.2011.03948.x.

    Article  CAS  PubMed  Google Scholar 

  44. Westermann AJ, Gorski SA, Vogel J: Dual RNA-seq of pathogen and host. Nat Rev Microbiol. 2012, 10 (9): 618-630. 10.1038/nrmicro2852.

    Article  CAS  PubMed  Google Scholar 

  45. Strickler SR, Bombarely A, Mueller LA: Designing a transcriptome next-generation sequencing project for a nonmodel plant species. Am J Bot. 2012, 99 (2): 257-266. 10.3732/ajb.1100292.

    Article  CAS  PubMed  Google Scholar 

  46. Jørgensen IH: Discovery, characterization and exploitation of Mlo powdery mildew resistance in barley. Euphytica. 1992, 63 (1–2): 141-152. 10.1007/BF00023919.

    Article  Google Scholar 

  47. Kim DS, Hwang BK: The pepper MLO gene, CaMLO2, is involved in the susceptibility cell-death response and bacterial and oomycete proliferation. Plant J. 2012, 72 (5): 843-855. 10.1111/tpj.12003.

    Article  CAS  PubMed  Google Scholar 

  48. Zheng Z, Nonomura T, Appiano M, Pavan S, Matsuda Y, Toyoda H, Wolters A-MA, Visser RGF, Bai Y: Loss of function in Mlo orthologs reduces susceptibility of pepper and tomato to powdery mildew disease caused by Leveillula taurica . PLoS One. 2013, 8 (7): e70723-10.1371/journal.pone.0070723.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. McGrann GR, Stavrinides A, Russell J, Corbitt MM, Booth A, Chartrain L, Thomas WT, Brown JK: A trade off between mlo resistance to powdery mildew and increased susceptibility of barley to a newly important disease, Ramularia leaf spot. J Exp Bot. 2014, 65 (4): 1025-1037. 10.1093/jxb/ert452.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Nicaise V, Roux M, Zipfel C: Recent advances in PAMP-triggered immunity against bacteria: pattern recognition receptors watch over and raise the alarm. Plant Physiol. 2009, 150 (4): 1638-1647. 10.1104/pp.109.139709.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Kim MC, Panstruga R, Elliott C, Muller J, Devoto A, Yoon HW, Park HC, Cho MJ, Schulze-Lefert P: Calmodulin interacts with MLO protein to regulate defence against mildew in barley. Nature. 2002, 416 (6879): 447-451. 10.1038/416447a.

    Article  CAS  PubMed  Google Scholar 

  52. Mayer C, et al: Phobos 3.3.11. 2010 [http://www.rub.de/spezzoo/cm/cm_phobos.htm] ,

    Google Scholar 

  53. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG: Primer3-new capabilities and interfaces. Nucleic Acids Res. 2012, 40 (15): e115-10.1093/nar/gks596.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Vaz Patto MC, Fernández-Aparicio M, Moral A, Rubiales D: Characterization of resistance to powdery mildew (Erysiphe pisi) in a germplasm collection of Lathyrus sativus . Plant Breed. 2006, 125: 308-310. 10.1111/j.1439-0523.2006.01220.x.

    Article  Google Scholar 

  55. Martín-Sanz A, Palomo J, Pérez de la Vega M, Caminero C: Characterization of Pseudomonas syringae pv. syringae isolates associated with bacterial blight in Lathyrus spp. and sources of resistance. Eur J Plant Pathol. 2012, 134: 205-216. 10.1007/s10658-012-9980-x.

    Article  Google Scholar 

  56. Hacquard S, Veneault-Fourrey C, Delaruelle C, Frey P, Martin F, Duplessis S: Validation of Melampsora larici-populina reference genes for in planta RT-quantitative PCR expression profiling during time-course infection of poplar leaves. Physiol Mol Plant Pathol. 2011, 75 (3): 106-112. 10.1016/j.pmpp.2010.10.003.

    Article  CAS  Google Scholar 

  57. Almeida NF, Leitão ST, Caminero C, Torres AM, Rubiales D, Vaz Patto MC: Transferability of molecular markers from major legumes to Lathyrus spp. for their application in mapping and diversity studies. Mol Biol Rep. 2014, 41 (1): 269-283. 10.1007/s11033-013-2860-4.

    Article  CAS  PubMed  Google Scholar 

  58. Varshney RK, Graner A, Sorrells ME: Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 2005, 23 (1): 48-55. 10.1016/j.tibtech.2004.11.005.

    Article  CAS  PubMed  Google Scholar 

  59. Consonni C, Bednarek P, Humphry M, Francocci F, Ferrari S, Harzen A, Ver Loren van Themaat E, Panstruga R: Tryptophan-derived metabolites are required for antifungal defense in the Arabidopsis mlo2 mutant. Plant Physiol. 2010, 152 (3): 1544-1561. 10.1104/pp.109.147660.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Sillero JC, Rojas-Molina MM, Avila CM, Rubiales D: Induction of systemic acquired resistance against rust, ascochyta blight and broomrape in faba bean by exogenous application of salicylic acid and benzothiadiazole. Crop Prot. 2012, 34: 65-69. 10.1016/j.cropro.2011.12.001.

    Article  CAS  Google Scholar 

  61. Barilli E, Sillero JC, Rubiales D: Induction of systemic acquired resistance in pea against rust (Uromyces pisi) by exogenous application of biotic and abiotic inducers. J Phytopathol. 2010, 158 (1): 30-34. 10.1111/j.1439-0434.2009.01571.x.

    Article  CAS  Google Scholar 

  62. Son GH, Wan J, Kim HJ, Nguyen XC, Chung WS, Hong JC, Stacey G: Ethylene-Responsive Element-Binding Factor 5, ERF5, is involved in chitin-induced innate immunity response. Mol Plant-Microbe Interact. 2011, 25 (1): 48-60. 10.1094/MPMI-06-11-0165.

    Article  Google Scholar 

  63. Gutterson N, Reuber TL: Regulation of disease resistance pathways by AP2/ERF transcription factors. Curr Opin Plant Biol. 2004, 7 (4): 465-471. 10.1016/j.pbi.2004.04.007.

    Article  CAS  PubMed  Google Scholar 

  64. Melotto M, Underwood W, Koczan J, Nomura K, He SY: Plant stomata function in innate immunity against bacterial invasion. Cell. 2006, 126 (5): 969-980. 10.1016/j.cell.2006.06.054.

    Article  CAS  PubMed  Google Scholar 

  65. Qin X, Zeevaart JAD: The 9-cis-epoxycarotenoid cleavage reaction is the key regulatory step of abscisic acid biosynthesis in water-stressed bean. Proc Natl Acad Sci U S A. 1999, 96 (26): 15354-15361. 10.1073/pnas.96.26.15354.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Cao FY, Yoshioka K, Desveaux D: The roles of ABA in plant-pathogen interactions. J Plant Res. 2011, 124 (4): 489-499. 10.1007/s10265-011-0409-y.

    Article  CAS  PubMed  Google Scholar 

  67. Selitrennikoff CP: Antifungal proteins. Appl Environ Microbiol. 2001, 67 (7): 2883-2894. 10.1128/AEM.67.7.2883-2894.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  68. Guo Z, Bonos S, Meyer W, Day P, Belanger F: Transgenic creeping bentgrass with delayed dollar spot symptoms. Mol Breed. 2003, 11 (2): 95-101. 10.1023/A:1022458101221.

    Article  CAS  Google Scholar 

  69. Rojas-Molina MD, Rubiales D, Prats E, Sillero JC: Effects of phenylpropanoid and energetic metabolism inhibition on faba bean resistance mechanisms to rust. Phytopathology. 2007, 97 (1): 60-65. 10.1094/PHYTO-97-0060.

    Article  CAS  Google Scholar 

  70. Chen X, Coram T, Huang X, Wang M, Dolezal A: Understanding molecular mechanisms of durable and non-durable resistance to stripe rust in wheat using a transcriptomics approach. Curr Genomics. 2013, 14 (2): 111-126. 10.2174/1389202911314020004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  71. Liu CJ, Huhman D, Sumner LW, Dixon RA: Regiospecific hydroxylation of isoflavones by cytochrome p450 81E enzymes from Medicago truncatula . Plant J. 2003, 36 (4): 471-484. 10.1046/j.1365-313X.2003.01893.x.

    Article  CAS  PubMed  Google Scholar 

  72. Shimada S, Takahashi K, Sato Y, Sakuta M: Dihydroflavonol 4-reductase cDNA from non-anthocyanin-producing species in the Caryophyllales. Plant Cell Physiol. 2004, 45 (9): 1290-1298. 10.1093/pcp/pch156.

    Article  CAS  PubMed  Google Scholar 

  73. Morant AV, Bjarnholt N, Kragh ME, Kjærgaard CH, Jørgensen K, Paquette SM, Piotrowski M, Imberty A, Olsen CE, Møller BL, Bak S: The β-glucosidases responsible for bioactivation of hydroxynitrile glucosides in Lotus japonicus . Plant Physiol. 2008, 147 (3): 1072-1091. 10.1104/pp.107.109512.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Marrs KA: The functions and regulation of glutathione S-transferases in plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 127-158. 10.1146/annurev.arplant.47.1.127.

    Article  CAS  PubMed  Google Scholar 

  75. Ghosh M: Antifungal properties of haem peroxidase from Acorus calamus . Ann Bot. 2006, 98 (6): 1145-1153. 10.1093/aob/mcl205.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. Link TI, Lang P, Scheffler BE, Duke MV, Graham MA, Cooper B, Tucker ML, van de Mortel M, Voegele RT, Mendgen K, Baum TJ, Whitham SA: The haustorial transcriptomes of Uromyces appendiculatus and Phakopsora pachyrhizi and their candidate effector families. Mol Plant Pathol. 2014, 15 (4): 379-393. 10.1111/mpp.12099.

    Article  CAS  PubMed  Google Scholar 

  77. Kemen E, Kemen A, Ehlers A, Voegele R, Mendgen K: A novel structural effector from rust fungi is capable of fibril formation. Plant J. 2013, 75 (5): 767-780. 10.1111/tpj.12237.

    Article  CAS  PubMed  Google Scholar 

  78. Giraldo MC, Valent B: Filamentous plant pathogen effectors in action. Nat Rev Microbiol. 2013, 11 (11): 800-814. 10.1038/nrmicro3119.

    Article  CAS  PubMed  Google Scholar 

  79. Vleeshouwers VG, Raffaele S, Vossen JH, Champouret N, Oliva R, Segretin ME, Rietman H, Cano LM, Lokossou A, Kessel G, Pel MA, Kamoun S: Understanding and exploiting late blight resistance in the age of effectors. Annu Rev Phytopathol. 2011, 49: 507-531. 10.1146/annurev-phyto-072910-095326.

    Article  CAS  PubMed  Google Scholar 

  80. Marino D, Peeters N, Rivas S: Ubiquitination during plant immune signaling. Plant Physiol. 2012, 160 (1): 15-27. 10.1104/pp.112.199281.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  81. Karasov TL, Horton MW, Bergelson J: Genomic variability as a driver of plant–pathogen coevolution?. Curr Opin Plant Biol. 2014, 18: 24-30. 10.1016/j.pbi.2013.12.003.

    Article  PubMed Central  PubMed  Google Scholar 

  82. Vaz Patto MC, Rubiales D: Lathyrus diversity: available resources with relevance to crop improvement – L. sativus and L. cicera as case studies. Ann Bot. 2014, 113 (6): 895-908. 10.1093/aob/mcu024.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  83. Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  84. Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26 (19): 2460-2461. 10.1093/bioinformatics/btq461.

    Article  CAS  PubMed  Google Scholar 

  85. Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18 (9): 1509-1517. 10.1101/gr.079558.108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  86. Wang L, Feng Z, Wang X, Zhang X: DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 26 (1): 136-138. 10.1093/bioinformatics/btp612.

    Article  PubMed  Google Scholar 

  87. Roth A, Ding J, Morin R, Crisan A, Ha G, Giuliany R, Bashashati A, Hirst M, Turashvili G, Oloumi A, Marra MA, Aparicio S, Shah SP: JointSNVMix: a probabilistic model for accurate detection of somatic mutations in normal/tumour paired next-generation sequencing data. Bioinformatics. 2012, 28 (7): 907-913. 10.1093/bioinformatics/bts053.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  88. Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, Duran C, Field M, Heled J, Kearse M, Markowitz S, Moir R, Stones-Havas S, Sturrock S, Thierer T, Wilson A, et al: 2011. Geneious v5.4, Available from http://www.geneious.com/,

    Google Scholar 

  89. Schuelke M: An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000, 18 (2): 233-234. 10.1038/72708.

    Article  CAS  PubMed  Google Scholar 

  90. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): RESEARCH0034-10.1186/gb-2002-3-7-research0034.

    Article  PubMed Central  PubMed  Google Scholar 

  91. Lohse M, Nagel A, Herter T, May P, Schroda M, Zrenner R, Tohge T, Fernie AR, Stitt M, Usadel B: Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data. Plant Cell Environ. 2014, 37 (5): 1250-1258. 10.1111/pce.12231.

    Article  CAS  PubMed  Google Scholar 

  92. Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37 (6): 914-939. 10.1111/j.1365-313X.2004.02016.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

The authors thank the CRF-INIA, Madrid, Spain, for supplying the genotypes. This work was supported by the project LEGRESIST (GEN2006-27798-C6-5-/VEG) and by Fundação para a Ciência e a Tecnologia through Grants #PEst-OE/EQB/LA0004/2011 and #PTDC/AGR-GPL/103285/2008. NFA and MCVP were supported by Fundação para a Ciência e a Tecnologia (SFRH/BD/44357/2008 and Research Contracts by the Ciência 2008 program respectively).

Author information

Authors and Affiliations

  1. Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa, Av. da República, Oeiras, 2780-157, Portugal

    Nuno Felipe Almeida, Susana Trindade Leitão & Maria Carlota Vaz Patto

  2. Frankfurt am Main, GenXPro GmbH, D-60438, Germany

    Nicolas Krezdorn, Björn Rotter & Peter Winter

  3. Institute for Sustainable Agriculture, CSIC, Córdoba, E-14080, Spain

    Diego Rubiales

Authors
  1. Nuno Felipe Almeida
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Susana Trindade Leitão
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Nicolas Krezdorn
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Björn Rotter
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Peter Winter
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Diego Rubiales
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Maria Carlota Vaz Patto
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Nuno Felipe Almeida.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

NFA carried out the inoculations and sample processing for sequencing and RT-qPCR, designed the RT-qPCR primers, performed the data analysis, developed the EST-SSRs primers and drafted the manuscript. NFA and STL participated in the results generation on RT-qPCR and data analysis. STL performed the EST-SSR genotyping. NK and BR performed the RNA sequencing and bioinformatics data processing. PW and DR revised the manuscript critically. MCVP coordinated the study and participated to the drafting of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12870_2014_376_MOESM1_ESM.xlsx

Additional file 1: Lists of genes in expression pattern groups mapped to the reference assembly. RPKM: reads per kilobase per million. (XLSX 2 MB)

12870_2014_376_MOESM2_ESM.xlsx

Additional file 2: List of contigs, mapped to the reference assembly, present only in the inoculated condition with a higher bit-score in fungal databases than in plant databases. (XLSX 587 KB)

Additional file 3: List of detected polymorphic SNPs between grass pea accessions BGE015746 and BGE024709. (XLSX 5 MB)

12870_2014_376_MOESM4_ESM.xlsx

Additional file 4: List of detected polymorphic EST-SSRs between grass pea genotypes BGE015746 and BGE024709. (XLSX 26 KB)

Additional file 5: Contig information and primer sequences for qRT-PCR. (XLSX 10 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Almeida, N.F., Leitão, S.T., Krezdorn, N. et al. Allelic diversity in the transcriptomes of contrasting rust-infected genotypes of Lathyrus sativus, a lasting resource for smart breeding. BMC Plant Biol 14, 376 (2014). https://doi.org/10.1186/s12870-014-0376-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-014-0376-2

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us