- Research article
- Open Access
The first set of EST resource for gene discovery and marker development in pigeonpea (Cajanus cajanL.)
© Raju et al; licensee BioMed Central Ltd. 2010
- Received: 21 August 2009
- Accepted: 11 March 2010
- Published: 11 March 2010
Pigeonpea (Cajanus cajan (L.) Millsp) is one of the major grain legume crops of the tropics and subtropics, but biotic stresses [Fusarium wilt (FW), sterility mosaic disease (SMD), etc.] are serious challenges for sustainable crop production. Modern genomic tools such as molecular markers and candidate genes associated with resistance to these stresses offer the possibility of facilitating pigeonpea breeding for improving biotic stress resistance. Availability of limited genomic resources, however, is a serious bottleneck to undertake molecular breeding in pigeonpea to develop superior genotypes with enhanced resistance to above mentioned biotic stresses. With an objective of enhancing genomic resources in pigeonpea, this study reports generation and analysis of comprehensive resource of FW- and SMD- responsive expressed sequence tags (ESTs).
A total of 16 cDNA libraries were constructed from four pigeonpea genotypes that are resistant and susceptible to FW ('ICPL 20102' and 'ICP 2376') and SMD ('ICP 7035' and 'TTB 7') and a total of 9,888 (9,468 high quality) ESTs were generated and deposited in dbEST of GenBank under accession numbers GR463974 to GR473857 and GR958228 to GR958231. Clustering and assembly analyses of these ESTs resulted into 4,557 unique sequences (unigenes) including 697 contigs and 3,860 singletons. BLASTN analysis of 4,557 unigenes showed a significant identity with ESTs of different legumes (23.2-60.3%), rice (28.3%), Arabidopsis (33.7%) and poplar (35.4%). As expected, pigeonpea ESTs are more closely related to soybean (60.3%) and cowpea ESTs (43.6%) than other plant ESTs. Similarly, BLASTX similarity results showed that only 1,603 (35.1%) out of 4,557 total unigenes correspond to known proteins in the UniProt database (≤ 1E-08). Functional categorization of the annotated unigenes sequences showed that 153 (3.3%) genes were assigned to cellular component category, 132 (2.8%) to biological process, and 132 (2.8%) in molecular function. Further, 19 genes were identified differentially expressed between FW- responsive genotypes and 20 between SMD- responsive genotypes. Generated ESTs were compiled together with 908 ESTs available in public domain, at the time of analysis, and a set of 5,085 unigenes were defined that were used for identification of molecular markers in pigeonpea. For instance, 3,583 simple sequence repeat (SSR) motifs were identified in 1,365 unigenes and 383 primer pairs were designed. Assessment of a set of 84 primer pairs on 40 elite pigeonpea lines showed polymorphism with 15 (28.8%) markers with an average of four alleles per marker and an average polymorphic information content (PIC) value of 0.40. Similarly, in silico mining of 133 contigs with ≥ 5 sequences detected 102 single nucleotide polymorphisms (SNPs) in 37 contigs. As an example, a set of 10 contigs were used for confirming in silico predicted SNPs in a set of four genotypes using wet lab experiments. Occurrence of SNPs were confirmed for all the 6 contigs for which scorable and sequenceable amplicons were generated. PCR amplicons were not obtained in case of 4 contigs. Recognition sites for restriction enzymes were identified for 102 SNPs in 37 contigs that indicates possibility of assaying SNPs in 37 genes using cleaved amplified polymorphic sequences (CAPS) assay.
The pigeonpea EST dataset generated here provides a transcriptomic resource for gene discovery and development of functional markers associated with biotic stress resistance. Sequence analyses of this dataset have showed conservation of a considerable number of pigeonpea transcripts across legume and model plant species analysed as well as some putative pigeonpea specific genes. Validation of identified biotic stress responsive genes should provide candidate genes for allele mining as well as candidate markers for molecular breeding.
- Simple Sequence Repeat Marker
- Fusarium Wilt
- Polymorphic Information Content
- Cleave Amplify Polymorphic Sequence
- Cleave Amplify Polymorphic Sequence Marker
Pigeonpea (Cajanus cajan (L.) Millsp) is one of the major grain legume crops of the tropical and subtropical regions of the world . It is the only cultivated food crop of the Cajaninae sub-tribe and has a diploid genome with 11 pairs of chromosomes (2n = 2× = 22) and a genome size estimated to be 858 Mbp . The genus Cajanus comprises 32 species most of which are found in India, Australia and one is native to West Africa. Pigeonpea is a major food legume crop in South Asia and East Africa with India as the largest producer (3.5 Mha) followed by Myanmar (0.54 Mha) and Kenya (0.20 Mha) . It plays an important role in food security, balanced diet and alleviation of poverty because of its diverse usages as a food; fodder and fuel wood . Several abiotic (e.g. drought, salinity and water-logging) and biotic (e.g. diseases like Fusarium wilt, sterility mosaic and pod borer insects) stresses, are serious challenges for sustainable pigeonpea production to meet the demands of the resource poor people of several African and Asian countries.
Fusarium wilt (FW) caused by Fusarium udum is an important biotic constraint in pigeonpea production in the Indian subcontinent, which results in 16-47% crop losses . The fungus enters the host vascular system at the root tips through wounds or invasion made by nematodes, leading to progressive chlorosis of leaves, branches, wilting and collapse of the root system . In India alone, the loss due to this disease is estimated to be US $71 million and the percentage of disease incidence varies from 5.3 to 22.6% .
Sterility mosaic disease (SMD) caused by pigeonpea sterility mosaic virus (PPSMV) is one of the wide-spread diseases of pigeonpea, which is transmitted by an eriophyid mite (Aceria cajani Channabasavanna). The disease is characterized by the symptoms like bushy and pale green appearance of plants followed by reduction in size, increase in number of secondary and mosaic mottling of leaves and finally partial or complete cessation of reproductive structures. Some parts of the plant may show disease symptoms and other parts may remain unaffected .
Due to the above mentioned factors combined with limited water resources to the fields in the semi-arid tropic regions, where the crop is grown, the productivity has remained stagnant at around 0.7 t/ha during the past two decades . With the advent of genomic tools such as molecular markers, genetic maps, etc., conventional plant breeding has been facilitated greatly and improved genotypes/varieties with enhanced resistance/tolerance to biotic/abiotic stresses have been developed in several crop species [9, 10]. In case of pigeonpea, however, a very limited number of genomic tools are available so far [11, 12]. For instance, 156 microsatellite or simple sequence repeat (SSR) markers [13–16], 908 expressed sequence tags (ESTs), at the time of undertaking the study, were available in pigeonpea. For enhancing the genomic resources in pigeonpea, transcriptome sequencing to generate ESTs should be a fast approach. ESTs, which are generated by large-scale single pass sequencing of randomly picked cDNA clones, have been cost - effective and valuable resource for efficient and rapid identification of novel genes and development of molecular markers . Further, ESTs have been employed in bioinformatic analyses to identify the genes that are differentially expressed in various tissues, cell types, or developmental stages of the same or different genotypes [18, 19].
In view of above facts, this study was undertaken to obtain a comprehensive resource of FW- and SMD-responsive ESTs in pigeonpea with the following objectives: (i) generation of FW- and SMD- responsive ESTs, (ii) functional annotation of assembled unigenes, (iii) in silico identification of putative FW- and SMD- responsive genes, and (iv) development of novel SSR and SNP markers in pigeonpea.
Generation of FW- and SMD- responsive ESTs
Pigeonpea EST assembly
Comparison of pigeonpea unigenes with other plant EST databases
BLASTN analyses of pigeonpea unigenes against legume and model plant ESTs
High quality ESTs generated
Pigeonpea (Cajanus cajan) (908)
Chickpea (Cicer arietinum) (7,097)
Soybean (Glycine max) (880,561)
Cowpea (Vigna unguiculata) (183,757)
Medicago (Medicago truncatula) (249,625)
Lotus (Lotus japonicus) (183,153)
Common bean (Phaseolus vulgaris) (83,448)
Significant similarity with ESTs of at least
one legume species
Significant similarity across legume ESTs
No similarity with legume species
Model plant ESTs
Arabidopsis (Arabidopsis thaliana) (1,527,298)
Rice (Oryza sativa) (1,240,613)
Poplar (Poplus alba) (418,223)
Significant similarity with ESTs of at least one
Model plant species
Significant similarity across ESTs of all model plant
Significant similarity with ESTs of at least one
plant species analyzed
Significant similarity across ESTs of all plant
No similarity with ESTs of any plant species
Functional categorization of pigeonpea unigenes
In silicoexpression analysis
In Cluster I, 11.3% (8) of total genes were grouped and further sub divided into two groups with each sharing 2.8% (2) and 8.5% (6) genes, respectively. Similarly, Cluster II and Cluster III accounted for 4.2% (3) and 15.5% (11) genes and the largest Cluster IV, included 69.0% (49) of total genes with three sub groups IVa, IVb and IVc each sharing 14.0% (10), 10% (7) and 45% (32) of genes, respectively. Cluster analysis also showed high level expression of genes related to chloroplast/photosystem related proteins (22.5%), developmental proteins (19.7%), cellular proteins (15.4%), metabolic proteins (14.0%), defence/stimulus responsive proteins (4.3%), protein specific binding proteins (2.8%) and few uncharacterized proteins (19.8%).
EST based markers can assay the functional genetic variation compared to other class of genetic markers and hence were targeted for marker development . The unigene set based on generated ESTs in this study as well as the ones available in public domain was used for development of simple sequence repeats (SSR) and single nucleotide polymorphism (SNP) markers.
Identification and development of genic microsatellite markers
Features of SSRs identified in ESTs
SSR database mining
Total number of sequences examined
Total length of examined sequences (bp)
Number of ESTs containing SSRs
Number of identified SSRs
Number of sequences containing more than 1 SSR
Number of SSRs present in compound formation
Frequency of SSR
Distribution of SSRs
In general, mononucleotide SSRs are not included for primer designing and synthesis. However, as only a very limited number of SSR markers are currently available for pigeonpea in public domain and in a separate study some mononucleotide SSRs were found polymorphic , primer pairs were designed for 383 SSRs including mononucleotide SSRs. A total of 94 primer pairs were considered for validation after excluding the primers for monomeric SSR motifs and compound SSRs with mononucleotide repeats. However based on repeat number criteria, such as 5 minimum for di-, tri-, tetra-, penta-nucleotides, primer pairs were synthesized for 84 SSRs. The details of newly developed pigeonpea EST-SSR primers along with corresponding SSR motif, primer sequence, annealing temperature and product size are provided in Additional file 12.
Characteristics of pigeonpea EST-SSR markers
Product size (bp)
No. of alleles
SNP discovery and identification of CAPS markers
Summary of SNPs and CAPS markers identified from pigeonpea ESTs
Total number of contigs examined (UG-IV)
Number of contigs containing ≥ 5 ESTs
Number of contigs containing SNPs
Total length of 37 contigs (bp)
Total number of identified SNPs in 37 contigs
Average SNP frequency
Total number of contigs containing CAPS convertible SNPs
In order to perform cost-effective and robust genotyping assay for the detected 102 SNPs in 37 contigs, efforts were made to identify the restriction enzymes that can be used to assay SNPs via cleaved amplified polymorphic sequence (CAPS) assay. Results indicated that SNPs present in 37 contigs can be evaluated by using CAPS assay (Table 4).
Plants are known to have developed integrated defence mechanisms against fungal and viral infections by altering spatial and temporal transcriptional changes. The EST approach was successfully utilized in identification of disease-responsive genes from various tissues and growth stages in chickpea , Lathyrus , soybean , rice  and ginseng . Many earlier studies have shown that resistant genotypes have efficient mechanisms for stress perception and enhanced expression of defence-responsive genes, which maintain cellular survival and recovery . Hence, the present study was undertaken to identify catalog of defence related genes in response to FW and SMD infection in pigeonpea by generating ESTs from different stress challenged tissues at various time intervals.
Generation of cDNA libraries and unigene assemblies
Roots provide a structural and physiological support for plant interactions with the soil environment by conducting transport of water, ions and nutrients. Plants are encountered with many biotic stress factors which includes bacterial, fungal and viral infection. Roots and leaves are the primary sites of infection by these organisms. Therefore, a total of 16 cDNA libraries were generated at different time intervals to specifically target the roots infected with Fusarium udum and leaves infected with SMD. In total 5,680 high quality ESTs were generated from FW- and similarly 3,788 high quality ESTs from SMD- challenged genotypes. Earlier, at the time of analysis in November 2008, the public domain consisted of only 908 ESTs for pigeonpea. Thus the present study contributes approximately 10-fold increase in the pigeonpea EST resource and an addition of 4,557 pigeonpea unigenes (UG-III).
Functional annotation of pigeonpea unigenes
Homology searches (BLASTN and BLASTX) against other plant ESTs and functional characterization was done for all the defined unigene datasets (UG-I, UG-II, UG-III and UG-IV). Of the 5,085 unigenes (UG-IV) assembled from all the pigeonpea ESTs, 3,280 (64.5%) had significant identity with ESTs of at least one plant species analyzed, 299 (5.8%) unigenes showed significant identity with ESTs of all analyzed plant species in the study, while 41 (0.8%) were found to be novel to pigeonpea. A high significant identity was observed with soybean (56.3%), and the least percentage of similarity was observed with chickpea (22.7%) (Table 1). A similar BLASTN results were observed for the remaining three unigenes sets (UG-I, UG-II and UG-III) against the ESTs of plant species surveyed. Comparative analysis of newly defined UG-III dataset (4,557) with 908 public domain pigeonpea ESTs showed that only 508 (11.1%) shared identity and indicated that our EST sequencing study identified 4,049 (88.9%) new set of pigeonpea unigenes. Relatively, very low similarity of 36.5% with Lotus and 42.3% with Medicago was observed compared to soybean and cowpea than other legume species. These observations are in accordance with phylogenetic relationships of legumes .
The pigeonpea ESTs showed higher similarity to legume ESTs databases (22.7-56.3%) of the legume species than monocot species (27.3-33.4%). Comparative analysis of pigeonpea ESTs with monocot species like rice (27.3%) showed that the percentage of significance is much lower compared to any other legume species, inspite of larger EST repository. This is clearly attributed to phylogenetic divergence between dicots and monocots in course of evolution. These comparisons also indicate that several unigenes that were absent in analysed non-legumes but present in all legume species may be specifically confined to legumes.
BLASTX analyses indicated that those ESTs without significant identity to any other protein sequences in the existing database may be novel and involved in plant defence responses. Hence, this novel EST collection represented a significant addition to the existing pigeonpea EST resources and provides valuable information for further predictions/validation of gene functions in pigeonpea.
A comprehensive comparison of functionally categorized unigenes of all the four unigenes data sets (UG-I, UG-II, UG-III and UG-IV) showed a similar distribution. A large number of unigenes were involved in cell part, organelle, binding, organelle part, metabolic and cellular process among the significantly annotated ones. These observations are consistent with the earlier reported functional categorization studies in rice , soybean , barley  and tall fescue . However, the sequences encoding activities related to categories such as biological regulation and response to stimulus are 28 and 20 incase of FW-responsive ESTs compared to 0 and 2 in case of SMD-responsive ESTs. This was possibly due to the fact that the ESTs generated from FW- challenged root libraries were most abundantly involved in stimulus to pathogenesis and ESTs derived from SMD stress are chloroplast binding proteins. Earlier studies such as Lee and colleagues , Ablett and colleagues , also reported that photosynthesis-related proteins were the most prevalent from aerial parts of the plant, which would help to make energy related activities such as cell division, growth, elongation and development. Similarly in this study, photosynthesis related genes were identified in larger proportion (30%) in SMD-responsive cDNA libraries derived from leaf tissues.
In silicodifferential gene expression
The invasion of pathogen not only results in expression of novel genes/transcripts, but also in altering the abundances of different ESTs resulting in induction or repression. This was evident from differential expression of 19 genes between FW-responsive genotypes and 20 genes between SMD-responsive genotypes. It is however, important to mention that in silico method of gene expression is not the ideal method to identify the differentially expressed genes. Nevertheless, as large scale EST data were generated from FW- responsive and SMD- responsive genotypes, an effort like some earlier studies [18–21, 34] was made to identify some putative genes differentially expressed in FW- and SMD- resistant and sensitive genotypes. Validation of these candidate genes by Northern analysis or real-time quantitative PCR analysis is essential before these candidate genes are deployed in some other studies.
Significant number of unigene sequences related to proteins like kinases, phosphatases, peroxidases, ribonucleases, endochitinases, glucanases and hormones like Abscisic acid responsive (ABR) genes were identified to be differentially expressed and are known to play a vital role in defence mechanism. For example, the cell wall degrading enzymes like endochitinases (EC: 18.104.22.168) implicate a major defence mechanism against pathogen . Similarly, kinases play a major role in the plant's recognition to pathogen [38, 39]. For instance, chitinase protein (UniProt ID: P23472), a class of pathogenesis related (PR) proteins with bi-functional role in lysozyme/chitinase activity involved in random hydrolysation of N-aetyl-beta-D-glucosaminide-beta linkages in chitin and chitodextrins during systemic acquired resistance (SAR), was expressed at higher concentrations in FW-responsive resistant genotype ('ICPL 20102') compared to susceptible genotype ('ICP 2376'). The high expression levels of chitinase in resistant genotype indicate the effectiveness within a narrow range of pathogenesis [40, 41].
Similarly, the protein coding for ABA-responsive protein (ABR18) (UniProt ID: Q06930), which is involved in stimulus mechanism and cell localization etc. during plant development and one of the vital roles is in defence mechanism during biotic stress signaling. This gene was identified to be expressed relatively higher in SMD-resistant pigeonpea genotype 'ICP 7035' compared to the susceptible genotype 'TTB 7'. During pathogen infection ABA inhibits the transcription of a basic β-1, 3-glucanase (EC: 22.214.171.124) that can degrade the β-1, 3-glucan callose, forming a physical barrier to viral spread through plasmodesmata. This down regulation of β-1, 3-glucanase by ABA can be termed as a resistance factor in plant pathogen interactions . In our study, significant expression signals were observed in SMD resistant genotype 'ICP 7035' during viral infection. This positive correlation between the ABA levels and disease resistance was reported in plant species like common bean , rice  and tobacco . Different enzymes like methyltransferases (HMT3) (UniProt ID: Q8LAX0) and dehydrogenases (G3PC) (UniProt ID: P34921) are putatively involved in synthesis of lignin in cell walls. These enzymes also play a major role in defence against pathogen interaction [46, 47].
An in silico hierarchical clustering analysis of 71 differentially expressed and genes across 16 cDNA libraries using HCE V 2.0 was done to infer potential relation between the co-expressed genes. The profiles of some of the interesting gene families and genes that could play an important role in stress stimulus were explained.
In Cluster I, of the 8 contigs, 6 were identified to be highly expressed in FW- challenged libraries of susceptible genotype. The cluster includes genes encoding proteins involved in mitochondrial DNA (mtDNA) stability, Na+/H+ exchanger and a few uncharacterized proteins. The sub cluster Ia includes two genes, which are highly expressed in 'ICPL 20102' libraries (15 DAI and 25 DAI). The genes connected in sub cluster Ib are highly expressed in 'ICP 2376' libraries (6 DAI and 15 DAI). One of the putative proteins TAR1 (Transcript Antisense to Ribosomal RNA), a mitochondrial protein is known to be involved in regulation and respiratory metabolism. Over- expression of this protein suppresses the respiration-deficient petite phenotype of a point mutation in mitochondrial RNA polymerase that affects mitochondrial gene expression and mtDNA stability. This dysfunction of mitochondria might occur in response to biotic or abiotic stress . The over-expression of these genes was observed only in 15 DAI libraries of both the FW- responsive genotypes. And their immediate disappearance in the later stages of infection in resistant genotype libraries and continued expression in susceptible genotypes supports a hypothesis that continual expression of this protein may lead to mitochondrial dysfunction and subsequent cell degeneracy.
Another protein Na+/H+ exchanger 7 (UniProt ID: Q8BLV3) is an ubiquitous ion transporter that serves multiple cell physiological processes such as intracellular pH homeostasis and electro neutral exchange of protons for Na+ and K+ across endomembranes. Biochemical studies suggest that Na+/H+ exchangers in the plasma membrane of plant cells contribute to cellular sodium homeostasis during salt stress, although the above protein is expressed in high salt stressed plants, it may also be expressed during biotic stress. Its high expression in susceptible genotype 'ICPL 2376' at 6 and 15 DAI libraries shows the severity of stress during fungal pathogenesis. And in the remaining FW- and SMD- challenged libraries this shows normal expression.
Cluster II genes include hevamine (EC: 126.96.36.199) and Leucoanthocyanidin dioxygenase (EC: 188.8.131.52) specifically expressed in 'ICPL 20102' 30 DAI library. The important protein hevamine represents a new class of polysaccharide-hydrolyzing (βα)8 barrel enzyme belonging to families of plant chitinases and lysozymes, which are vital for plant defence against pathogenic bacteria and fungi. Recent results indicate that these enzymes may be involved not only in defence-related process or general stress response but also in growth and development processes . The high expression of these proteins in the late phase of Fusarium infection indicates their prolonged defensive role against fungal pathogenesis.
The genes connected in Cluster III are highly expressed in FW-challenged 'ICP 2376' - 25 and 30 DAI libraries. The genes related to endoprotease activity, beta-glucan synthesis, carboxypeptidases, alkaloid biosynthesis, secologanin biosynthesis, beta-galactosidase activity, microtubule-stabilizing activity, nucleic acid binding protein, ribonucleoprotein and a few uncharacterized proteins were constituted in this cluster. For instance, antifungal class proteins such as beta-glucanases (EC: 184.108.40.206) and Hansenula mrakii killer toxin-resistant proteins (UniProt ID: P41809) located in epidermal leaf cells are believed to be involved in cell differentiation and defence against fungal pathogens . Plants deficient in these enzymes generated by antisense transformation showed markedly reduced resistance to viral and fungal infection. Similarly, another class of proteins carboxypeptidases (EC: 3.4.16 - 3.4.18) with diverse functions ranging from catabolism to regulating biological processes, including function as a defence against pathogen attack .
The majority of genes (49) segregated in Cluster IV were highly expressed in SMD-responsive cDNA libraries derived from leaf tissues. In total Cluster IV genes showed high gene expression in SMD derived libraries. As expected, photosynthesis related transcripts were abundantly represented in sub clusters IVa, IVb and IVc, and these include putative transcripts like ribosomal proteins, mitochondrial proteins, chloroplast precursor proteins, photosystem I and II reaction centre proteins. The observed expression pattern of photosynthesis related proteins in this study is also consistent with experimental observations in barley .
Overall, the differentially expressed genes are involved in diverse pathways, displaying complex expression patterns. The different clusters based on monitoring of gene expression patterns propose that various pathways in response to biotic stress exist in pigeonpea and their interaction can lead to differential stress tolerance. The uncharacterized class of transcripts co-expressed could be repository of novel proteins and further characterization of these may reveal their significant role in plant stress responses .
Development of functional markers
One of our primary goals of our research programme is to develop molecular markers based on expressed sequences and screen them for polymorphism. During the last decade, microsatellites or SSRs have proven to be useful markers in plant genetic research and have been used for marker-assisted breeding purposes. The presence of SSRs in the coding region suggests their importance as functional or gene based markers [1, 11, 53]. Unfortunately, development of microsatellite markers is expensive, labor intensive, and time consuming if they are being developed from genomic libraries . The data mining of microsatellites markers from EST data can be a cost effective option. The cost of mining EST libraries is far lower than other traditional methods, and SSR development from ESTs has been successful in EST data mining [22, 23, 53–56]. SSR motifs with repeats more than eight for di-nucleotides, six for tri-nucleotides, and five for tetra-nucleotides were considered. Dimeric repeat motifs (40) were relatively abundant than trimeric repeats (33). In addition to this, tetra-, penta- and hexameric repeat motifs were considerably less represented. A total of 94 SSR markers have been synthesized and characterized for polymorphism survey. However, there are some distant contrasts in frequency and distribution of SSRs in ESTs and in genomic survey sequences (GSSs). In general di-nucleotide SSRs of all repeat lengths are more common in GSSs and tri-nucleotide SSRs are common in the ESTs [22, 23, 56, 57]. As against these reports, in our findings we observed that di-nucleotide repeats are more abundant than tri-nucleotide repeat motifs [58, 59]. However this observation is not unexpected as the frequency and distribution of SSR depends on several factors such as size of dataset, tools and criteria used for SSR discovery .
In this study, a total of 15 polymorphic EST-SSRs primer pairs were validated and used for diversity study on forty pigeonpea genotypes representing 32 cultivated (C. cajan) and 8 wild species (six C. scarabaeoides and two C. platycarpus). All markers detected at least one allele in all genotypes tested, suggesting transferability of all markers across the Cajanus genus. In addition to high transferability, EST-SSRs are good candidates for the development of conserved orthologous sequence (COS) markers for genetic analysis and breeding of different species . However, EST-SSRs were reported to be less polymorphic than genomic SSRs in crop plants due to greater DNA sequence conservation in transcribed regions [22, 60]. For instance, the 15 SSR loci provided only 60 alleles with an average of 4 alleles per loci and an average 0.43 PIC value. Similar kind of diversity features were observed in earlier SSR based diversity studies in pigeonpea [14, 15].
EST-SSR profiles obtained on 40 pigeonpea genotypes were used to compute pair-wise genetic distances among different genotypes to construct a dendrogram based on UPGMA clustering. The neighbor joining tree grouped 40 pigeonpea genotypes into three major clusters (Figure 10). The Cluster I comprising 32 genotypes (cultivated) is the largest cluster followed by Cluster III containing six wild genotypes (C. scarabaeoides) and Cluster II is the smallest cluster with two wild genotypes belonging to C. platycarpus species revealing clear segregation of the cultivated and the wild species. Less genetic variation was detected with in cultivated species, with only nine markers detecting polymorphism and a total of 35 alleles. The low genetic variability amongst cultivars when compared with the wild species genotypes suggests that natural and artificial selection has contributed to the selection of specific alleles and to changes of allelic frequencies at specific loci as reported by Odeny and colleagues . The distinctness of C. platycarpus with C. scarabaeoides accessions observed in this study correlate well with earlier studies . It is also important to note that 'ICPL 20097' and 'ICP 2376' genotypes were found closely related with high genetic similarity as both of these genotypes belong to the same geographic region. In conclusion, EST-SSR markers developed in this study complement the currently available or ongoing efforts on development of genomic SSRs that will be a valuable resource for linkage map development and marker assisted selection in pigeonpea .
SNPs and indels are an essentially inexhaustible resource of polymorphic markers for use in the high-resolution genetic map development of traits and for association studies. Although a variety of molecular markers are available SNPs are comparatively advantageous because of their abundance and amenability to high throughput approaches . In addition, SNPs also offer several advantages like high-throughput and cost-effective genotyping  and identification of functional/gene-based markers for complex trait through linkage map development or association genetics [9, 10, 64, 65]. Although SNP discovery was a cost effective task in past, advances in next generation sequencing technologies have made SNP discovery cheaper and faster . However in case for a given species, ESTs are available from more than one genotype, in silico mining of ESTs is still a very inexpensive and fast approach for SNP discovery [17, 64] and therefore we used this approach for mining SNPs in this study.
By using in silico mining approach in a total of 871 contigs coming from 10,376 ESTs (9,888 generated in this study and 908 available in public domain), a total of 102 potential SNPs were identified in 37 contigs that were consisted of ≥ 5 ESTs. Smaller contigs were not considered for SNP mining as these contigs are prone to errors due to lack of read depth as reported by Wang and colleagues . Sequence analysis of PCR products for a subset of 6 out of 10 contigs confirmed the occurrence of SNPs in all the cases. As PCR products could not be generated for remaining four contigs, the presence of SNPs could not be confirmed in those cases. Furthermore, as SNP genotyping is another important criteria in breeding programmes, identification of CAPS markers for 37 contigs will facilitate SNP genotyping even in low tech laboratories .
This study has contributed a new and significant set of 9,888 ESTs that together with 908 public domain ESTs provides a unigene set of 5,085 sequences for pigeonpea. Detailed analysis of these datasets have provided several important features of pigeonpea transcriptome such as conserved genes (across legumes and model plant species) as well as possible pigeonpea specific genes, assignment of pigeonpea genes to different GO categories, identification of differentially expressed genes in response to FW- and SMD- stresses, etc. In terms of applied aspect of developed resource in breeding, this study has demonstrated development and application of gene-based molecular markers i.e SSRs, SNPs and CAPS. In summary, it is anticipated that this study is a significant contribution to enhance genomic resources in a so called orphan legume crop that will eventually impact pigeonpea breeding [11, 12].
Four pigeonpea genotypes namely 'ICPL 20102' (FW- resistant), 'ICP 2376' (FW- susceptible), 'ICP 7035' (resistant to SMD) and 'TTB 7' (highly susceptible to SMD) were used for constructing the cDNA libraries and generating the ESTs. Seeds of two genotypes ('ICPL 20102', 'ICP 2376') were procured from Legume Pathology section at ICRISAT and for the remaining two genotypes ('ICP 7035' and 'TTB 7') were obtained from Dr. M Byregowda, University of Agricultural Sciences, Bangalore, India.
A total of 40 genotypes including 32 genotypes from cultivated species (C. cajan) and 8 genotypes from 2 wild species (C. platycarpus and C. scarabaeoides) were used for validation and diversity analysis with new set of EST-SSR markers. These genotypes were obtained from Pigeonpea Breeding (Dr. KB Saxena) and Genebank (Dr. HD Upadhyaya) and have been listed in Additional file 13.
Inoculation treatment for FW and SMD
Seeds of FW-tolerant ('ICPL 20102') and FW-susceptible ('ICP 2376') were germinated in 15-inch deep polythene covers filled with sterile soil and sand (1:1) in a glass house at 23 ± 3°C under 80% relative humidity. The root, being the primary target of the pathogen Fusarium udum and the possible site of the initial defence response, was selected as the tissue of study. Ten days old seedlings were uprooted from pots and the root system was thoroughly washed in running tap water and rinsed with distilled water. Seedlings of each genotype were inoculated by immersing the roots for 2 min in fungal inoculum (Fusarium udum culture). The spore suspension at 6 × 105 conidia/ml, was made by adding fungal spores from several culture plates (Fusarium was grown on potato dextrose media supplemented with 0.25 μg/ml tetracycline). Immediately following the inoculation, the seedlings were transplanted to sterilized sand and soil mixture (1:1) in pots and were transferred to glass house.
In order to capture the genes expressed in resistant and susceptible genotypes at different time periods after inoculation, six stages i.e. 6, 10, 15, 20, 25 and 30 days after inoculation (DAI) were selected arbitrarily to construct the cDNA libraries. From 6-15 DAI, chlorosis symptoms were observed on the leaves and aerial parts of plant material, indicated the severity of Fusarium wilt disease. Furthermore, from each stage of days after inoculation, shoot and root section cuttings were made and observed the fungal penetration into the vascular tissues. Based on microscopic observations at different stages, initial symptoms of fungal infection were noticed in root vascular tissue at 15 and 20 DAI stages. Beyond 20 DAI, though the fungus penetrates deeper into the vascular tissues of susceptible and resistant varieties to some extent, the susceptible variety shows complete Fusarium symptoms where as the resistant variety prevents most of the attacking fungus from reaching maturity and developing symptoms.
For SMD study, highly susceptible ('TTB 7') and resistant ('ICP 7035') pigeonpea genotypes that are parents of a mapping population segregating for resistance to SMD were chosen. Forty seeds from each accession were sown in plastic bags filled with sterilized soil and were maintained in a glass house under optimal physiological conditions as described above. Ten days after sowing, the aerial parts of the seedlings were stapled with mosaic virus infected leaves. The viral disease is caused by pigeonpea sterility mosaic virus (PPSMV) and transmitted by an eriophyid mite Aceria cajani Channabasavanna . The disease slowly spreads into the vascular tissues from the aerial parts through mite population which is characterized by a bushy and pale green appearance of plants. Based on the severity of disease symptoms, leaves with visible SMD lesions were harvested at 45 and 60 days after sowing (DAS) stages for construction of cDNA libraries.
cDNA library construction
Root and leaf tissue samples were collected from FW- and SMD- responsive genotypes at different time-points till the infection stage reached stagnant phase. RNA was isolated from the above two tissue samples according to the protocol described by Schmitt and colleagues . RNA quality was assessed using formamide gel electrophoresis and poly (A)+ RNA was isolated with poly (A) tract mRNA isolation system IV (Promega, Madison, WI, USA) as described by the manufacturers. Double-strand cDNA was constructed using Super SMART™ PCR cDNA Synthesis kit (Clontech®, Mountain View, CA, USA) as described in the manufacturer's instructions. The resulting cDNA was size fractioned on 1.2% agarose gel. cDNA fractions containing fragments greater than 500 bp were selected for library construction. Subsequently, the cDNA was ligated into pGEM® Easy vector (Promega®, Madison, WI, USA) and ligation was allowed to proceed overnight at 14°C. The resulting plasmids were electroporated using One Shot® Top 10 Electrocomp™ cells (Invitrogen, Carlsbad, CA, USA). The transformants were spread on LB Agar plates containing 100 mg/ml ampicillin for direct picking. Based on blue/white screening, recombinant clones were picked into Nunc-Immuno™ 96 MicroWell™ Plates (Nunc™, Roskilde, Denmark) containing LB broth with 100 μg/ml ampicillin and grown for overnight at 37°C on a rotary shaker at 220 rpm. Glycerol stocks in 96-well format were prepared by combining 38 μl of 60% (v/v) glycerol with 150 μl of culture and frozen at -80°C.
EST sequencing, editing and assembly
Clones were randomly selected and on an average of 500 clones were sequenced per library in case of FW-response study and 1000 clones per library in case of SMD-responsive study. The plasmid DNA from these clones (i.e. colonies) was extracted using a 96-well alkaline lysis method prior to sequencing . Plasmid DNA sequencing was performed by commercial DNA sequencing service provider (Macrogen Inc., Korea) using the standard M13 forward primer.
The FASTA files containing the raw sequences were edited by the software Sequencher™ 4.0 (Gene Codes Corporation, Ann Arbor. MI, USA) to remove the vector sequences. The vector screened sequences were subjected to EST trimmer , to trim poly-A ends and low quality sequences. High quality sequences of >100 bp were selected for further sequence analysis. ESTs were clustered and aligned into contigs and singletons using the CAP3 program .
In order to assess the number of unique and overlapping transcripts among the 16 libraries, four data sets were generated; those derived from libraries constructed from of FW-responsive genotypes (UG-I); those derived from libraries constructed from of SMD-responsive genotypes (UG-II); combined dataset of FW- and SMD- responsive ESTs (UG-III); and also from public domain sequences with total generated ESTs in this study (UG-IV). In addition to the above assembly of unigene sets, CAP3 analysis was also performed to libraries derived from FW- resistant genotype, FW- susceptible genotype, SMD- resistant genotype and from SMD- susceptible genotype individually.
Homology search and functional annotation
The unigene sequences were also characterized for nucleotide homology search against the EST datasets of selected legume species [pigeonpea (Cajanus cajan)-908, chickpea (Cicer arietinum)-7,097, soybean (Glycine max)-880,561, Medicago (Medicago truncatula)-249,625, common bean (Phaseolus vulgaris)-83,448, cowpea (Vigna unguiculata)-183,757 and Lotus (Lotus japonicus)-183,153] and selected model plant species [rice (Oryza sativa)-1,240,613, Arabidopsis (Arabidopsis thaliana)-1,527,298 and poplar (Populus alba)-418,223] available at National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov) using BLASTN algorithm . A match was considered significant at E-value ≤ 1E-05.
Each unigene dataset was subjected to BLASTX analysis against the non-redundant protein database of UniProt to deduce a putative function. Sequence similarity was considered as significant at E-value ≤ 1E-08. Each unigene was assigned a putative cellular function based on the significant database hit with lowest e-value. Subsequently, unigenes that showed a significant BLASTX hit were used for functional annotation based on Gene Ontology categories from UniProt database (UniProt-GO). This process allowed assignment of unigenes to the GO functional categories of biological process, cellular component and molecular function. Distribution of unigenes was further investigated in terms of their assignment to sub-categories of the main GO categories.
In silico expression and hierarchical clustering
In order to identify the differentially expressed genes in FW- and SMD- responsive genotypes, 389 contigs coming from FW-responsive genotypes and 328 contigs coming SMD-responsive genotypes were analyzed by using IDEG.6 web interface tool [73, 74]. The IDEG.6 web tool allows running six different statistical analyses for the detection of differentially expressed genes in multiple tag experiments. For pair-wise comparisons, the Audic and Claverie test, Fisher exact test and chi-square tests (Χ2) were used and in multiple comparisons R- statistics test, Greller and Tobin test and chi-square tests (Χ2) were used [73, 74].
Further, gene expression analysis was performed with hierarchical clustering expression (HCE) version 2.0 beta software  using transcript abundance data from UG-III set that includes 697 contigs derived from both the stress responsive libraries. As a pre-requisite for HCE analysis, all 697 contigs were subjected to R statistics (R>8) and only those contigs (71) were selected that have; (i) minimum 5 ESTs, and (ii) differential abundance of ESTs coming from different libraries. The matrix file developed based on the frequency of ESTs to each of 71 contigs was used as input file for above mentioned HCE tool.
Identification and development of SSR markers
A total of 5,085 unigenes (unigene set, UG-IV) developed based on 9,888 ESTs generated in this study and 908 public domain ESTs were searched with a Perl script program, MISA (MIcroSAtellite) [23, 76] for identification and localization of SSRs. The SSR motifs, with repeat units more than five times in di-, tri-, tetra-, penta- and hexa- nucleotides were considered as SSR search criteria in MISA script. The Primer3 programme  was used for designing the primer pairs for SSRs and custom synthesized by MWG (MWG-Biotech AG, Bangalore, India).
The primer pairs for SSRs were tested for their utility as potential genetic markers on 40 elite genotypes of pigeonpea (Additional file 13). PCR amplifications were performed in 5 μl reactions containing 5 ng of genomic DNA, 1× SE-Taq DNA polymerase buffer (including 1.5 mM MgCl2), 2 mM dNTPs, 10 pmol of each primer and 0.1 U Taq DNA polymerase (SibEnzyme, Novosibirsk, Russia) with the following touch down profile; 3 min at 95°C; 5 cycles of 20 sec at 94°C, 20 sec at 60°C minus 1°C/cycle, 30 sec at 72°C; 40 cycles of 20 sec at 94°C, 20 sec at 56°C, 30 sec at 72°C; and 20 min at 72°C for final extension. PCR products were separated on 6% non-denaturing polyacrylamide gels for 3 h at 600 V and visualized by silver staining. The polymorphism information content (PIC) of individual EST-SSR markers was calculated by using the standard formula . Only data from polymorphic SSR loci were used for diversity analysis. Genetic similarities between any two genotypes were estimated according to Nei and Li . All 40 genotypes were clustered with the Unweighted Pair Group Method using arithmetic average (UPGMA) in the SAHN procedure of the NTSYS-PC v2.10t .
SNP detection and their conversion into CAPS
All 871 contigs obtained from the collection of 5,085 unigenes (UG-IV) were searched for putative SNP/indels by using an integrated pipeline for large scale SNP discovery [80, 81]. The pipeline utilized the CAP3 output files as input to detect SNPs/indels based on the nucleotide redundancy in the multiple sequence alignments. The auto SNP pipeline generated text file includes contig ID, number of sequences in the contig ID, consensus length, number of SNPs, mutation type and SNP frequency. The threshold for identification of SNPs was based on the number of sequences (≥ 5) in each consensus sequence and two or more sequences from different genotype. In order to verify the SNPs at sequence level, the PCR amplicons of all four genotypes were sequenced using the corresponding forward and reverse primers for a set of 10 contigs (see Additional file 14). The amplicons were purified and further sequencing was done as described . The sequenced data along with the sequences of ESTs (that provided the SNPs initially) were aligned and analyzed using BioEdit programme http://www.mbio.ncsu.edu/BioEdit/bioedit.html.
For converting SNPs into cleaved amplified polymorphic sequence (CAPS) markers, SNPs present in 37 contigs were analyzed to identify the recognition site for any of commercially available 725 restriction enzymes  by using integrated SNP2CAPS pipeline .
Authors are thankful to Indo-USAgricultural Knowledge Initiative (Indo-USA KI) supported by Indian Council of Agricultural Research (ICAR), Government of India and SP2-Leader Discretionary Grant from Generation Challenge Program http://www.generationcp.org for the financial support to undertake this study. Thanks are also due to Department of Biotechnology (DBT), Government of India for sponsoring a Post-Doctoral Fellowship to NLR. Authors are thankful to Dr. K.B. Saxena and Dr. H.D. Upadhayaya of ICRISAT for providing the seeds/DNA of some genotypes used in this study. Thanks are also due to Mr. A. Bhanu Prakash and Ms. Spurthi Nayak for their help in data analysis and discussions. Authors are also thankful to three anonymous reviewers for their valuable suggestions on the first version of the MS that helped in improvement of the MS.
- Varshney RK, Hoisington DA, Upadhyaya HD, Gaur PM, Nigam SN,Saxena KB, Vadez V, Sethy NK, Bhatia S, Aruna R, Channabyre Gowda MV,Singh NK: Molecular genetics and breeding of grain legume crops forthe semi-arid tropics. Genomic-Assisted crop improvement SpringerNetherlands, The NetherlandsVarshney RK, Tuberosa R 2007, 207-241 full_text.View ArticleGoogle Scholar
- Greilhuber J, Obermayer R: Genome size variation in Cajanus cajan (Fabaceae): a reconsideration. Plant Syst Evol. 1998, 212: 135-141. 10.1007/BF00985225.View ArticleGoogle Scholar
- FAOSTAT 2006. [http://faostat.fao.org]
- Rao SC, Coleman SW, Mayeux HS: Forage production and nutritive value of selected pigeonpea ecotypes in the southern great plains. Crop Sci. 2002, 42: 1259-1263.View ArticleGoogle Scholar
- Prasad P, Eswara Reddy NP, Anandam RJ, Lakshmikantha Reddy G: Isozymes variability among Fusarium udum resistant cultivars of pigeonpea (Cajanus cajan (L.) (Millsp). Acta Physiol Plant. 2003, 25: 221-228. 10.1007/s11738-003-0001-8.View ArticleGoogle Scholar
- Butler EJ: The wilt disease of pigeonpea and pepper. Agricult J India. 1906, 1: 25-26.Google Scholar
- Kannaiyan J, Nene YL, Reddy MV, Ryan JG, Raju TN: Prevalence of pigeonpea disease and associated crop losses in Asia, Africa and America. Trop Pest Manage. 1984, 30: 62-71. 10.1080/09670878409370853.View ArticleGoogle Scholar
- Kumar PL, Jones AT, Reddy DVR: A novel mite transmitted virus with a divided RNA genome closely associated with Pigeonpea sterility mosaic disease. Phytopathol. 2003, 93: 81-91. 10.1094/PHYTO.2003.93.1.71.View ArticleGoogle Scholar
- Varshney RK, Graner A, Sorrells E: Genomics-assisted breeding for crop improvement. Trends Plant Sci. 2005, 10: 621-630. 10.1016/j.tplants.2005.10.004.PubMedView ArticleGoogle Scholar
- Varshney RK, Hoisington DA, Tyagi AK: Advances in cereal genomics and applications in crop breeding. Trends Biotechnol. 2006, 11: 490-499. 10.1016/j.tibtech.2006.08.006.View ArticleGoogle Scholar
- Varshney RK, Close TJ, Singh NK, Hoisington DA, Cook DR: Orphan legume crops enter the genomics era. Curr Opin Plant Biol. 2009, 12: 202-210. 10.1016/j.pbi.2008.12.004.PubMedView ArticleGoogle Scholar
- Varshney RK, Penmetsa RV, Dutta S, Kulwal PL, Saxena RK, Datta S, Sharma TR, Rosen B, Carrasquilla-Garcia N, Farmer AD, Dubey A, Saxena KB, Gao J, Fakrudin B, Singh MN, Singh BP, Wanjari KB, Yuan M, Srivastava RK, Kilian A, Upadhyaya HD, Mallikarjuna N, Town CD, Bruening GE, He G, May GD, McCombie R, Jackson SA, Singh NK, Cook DR: Pigeonpea genomics initiative (PGI): an international effort to improve crop productivity of pigeonpea (Cajanus cajanL.). Mol Breed. 2009.Google Scholar
- Burns MJ, Edwards KJ, Newbury HJ, Ford-Lloyd BV, Baggott CD: Development of simple sequence repeat (SSR) markers for the assessment of gene flow and genetic diversity in pigeonpea (Cajanus cajan). Mol Ecol Notes. 2001, 1: 283-285. 10.1046/j.1471-8278.2001.00109.x.View ArticleGoogle Scholar
- Odeny DA, Jayashree B, Ferguson M, Hoisington D, Crouch J, Gebhardt C: Development, characterization and utilization of microsatellite markers in pigeonpea. Plant Breed. 2007, 126: 130-136. 10.1111/j.1439-0523.2007.01324.x.View ArticleGoogle Scholar
- Saxena RK, Prathima C, Saxena K, Hoisington DA, Singh NK, Varshney RK: Novel SSR markers for polymorphism detection in pigeonpea (Cajanusspp.). Plant Breed. 2009.Google Scholar
- Odeny DA, Jayashree B, Gebhardt C, Crouch J: New microsatellite markers for pigeonpea (Cajanus cajan (L.) millsp.). BMC Res Notes. 2009, 2: 35-10.1186/1756-0500-2-35.PubMedPubMed CentralView ArticleGoogle Scholar
- Sreenivasulu N, Kavi Kishor PB, Varshney RK, Altschmied L: Mining functional information from cereal genomes -the utility of expressed sequence tags. Curr Sci. 2002, 83: 965-973.Google Scholar
- Ogihara Y, Mochida K, Nemoto Y, Murai K, Yamazaki Y, Shin IT, Kohara Y: Correlated clustering and virtual display of gene expression patterns in the wheat life cycle by large-scale statistical analyses of expressed sequence tags. Plant J. 2003, 33: 1001-1011. 10.1046/j.1365-313X.2003.01687.x.PubMedView ArticleGoogle Scholar
- Ronning CM, Stegalkina SS, Ascenzi RA, Bougri O, Hart AL, Utterbach TR, Vanaken SE, Riedmuller SB, White JA, Cho J: Comparative analyses of potato expressed sequence tag libraries. Plant Physiol. 2003, 131: 419-429. 10.1104/pp.013581.PubMedPubMed CentralView ArticleGoogle Scholar
- Ewing RM, Kahla AB, Poirot O, Lopez F, Audic S, Claverie JM: Large-scale statistical analyses of rice ESTs reveal correlated patterns of gene expression. Genome. 1999, 9: 950-959. 10.1101/gr.9.10.950.View ArticleGoogle Scholar
- Stekel DJ, Git Y, Falciani F: The comparison of gene expression from multiple cDNA libraries. Genome Res. 2000, 10: 2055-2061. 10.1101/gr.GR-1325RR.PubMedPubMed CentralView ArticleGoogle Scholar
- Varshney RK, Graner A, Sorrells E: Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 2005, 23: 48-55. 10.1016/j.tibtech.2004.11.005.PubMedView ArticleGoogle Scholar
- Thiel T, Michalek W, Varshney RK, Graner A: Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003, 106: 411-422.PubMedGoogle Scholar
- Temnykh S, DeClerck G, Lukashova A, Lipocic L, Cartinhour S, McCouch S: Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): frequency, length variation, transposon associations, and genetic marker potential. Genome Res. 2001, 11: 1441-1452. 10.1101/gr.184001.PubMedPubMed CentralView ArticleGoogle Scholar
- Kantety RV, La Rota M, Matthews DE, Sorrells ME: Data mining for simple sequence repeats in expressed sequence tags from barley, maize, rice, sorghum and wheat. Plant Mol Biol. 2002, 48: 501-510. 10.1023/A:1014875206165.PubMedView ArticleGoogle Scholar
- Coram TE, Pang ECK: Isolation and analysis of candidate ascochyta blight defence genes in chickpea. Part I. Generation and analysis of an expressed sequence tag (EST) library. Physiol Mol Plant Path. 2005, 66: 192-200. 10.1016/j.pmpp.2005.08.003.View ArticleGoogle Scholar
- Skiba B, Ford R, Pang CK: 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 Path. 2005, 66: 55-67. 10.1016/j.pmpp.2005.04.005.View ArticleGoogle Scholar
- Iqbal MJ, Yaegashi S, Ahsan R, Shopinski KL, Lightfoot DA: Root response to Fusarium solani f. sp. glycines: temporal accumulation of transcripts in partially resistant and susceptible soybean. Theor Appl Genet. 2005, 110: 1429-1438. 10.1007/s00122-005-1969-9.PubMedView ArticleGoogle Scholar
- Jantasuriyarat C, Gowda M, Haller K, Hatfield J, Lu G, Stahlberg E, Zhou B, Li H, Kim H, Yu Y, Dean RA, Wing RA, Soderlund C, Wang GL: Large-Scale identification of expressed sequence tags involved in rice and rice blast fungus interaction. Plant Physiol. 2005, 138: 105-115. 10.1104/pp.104.055624.PubMedPubMed CentralView ArticleGoogle Scholar
- Goswami R, Punja ZK: Molecular and biochemical characterization of defense responses in ginseng (Panax quinquefolius) roots challenged with Fusarium equiseti. Physiol Mol Plant Path. 2008, 72: 10-20. 10.1016/j.pmpp.2008.04.006.View ArticleGoogle Scholar
- Reddy PCO, Sairanganayakulu G, Thippeswamy M, Reddy PS, Reddy MK, Sudhakar C: Identification of stress-induced genes from the drought tolerant semi-arid legume crop horsegram (Macrotyloma uniflorum (Lam.) Verdc.) through analysis of subtracted expressed sequence tags. Plant Sci. 2008, 175: 372-384. 10.1016/j.plantsci.2008.05.012.View ArticleGoogle Scholar
- Wojciechowski MF, Sanderson MJ, Steele KP, Liston A: Molecular phylogeny of the "temperate herbaceous tribes" of papilionoid legumes: a supertree approach. Advances in Legume Systematics. Edited by: Herendeen P, Bruneau A. 2000, Royal Botanic Gardens: Kew, 9: 277-298.Google Scholar
- Alkharouf N, Khan R, Matthews B: Analysis of expressed sequence tags from roots of resistant soybean infected by the soybean cyst nematode. Genome. 2004, 47: 380-388.PubMedView ArticleGoogle Scholar
- Zhang H, Sreenivasulu N, Weschke W, Stein N, Rudd S, Radchuk V, Potokina E, Scholz U, Schweizer P, Zierold U, Langridge P, Varshney RK, Wobus U, Graner A: Large-scale analysis of the barley transcriptome based on expressed sequence tags. Plant J. 2004, 40: 276-290. 10.1111/j.1365-313X.2004.02209.x.PubMedView ArticleGoogle Scholar
- Mian MAR, Zhang Y, Wang Z, Zhang J, Cheng X, Chen L, Chekhovskiy K, Dai X, Mao C, Cheng F, Zhao X, He J, Scott AD, Town CD, May GD: Analysis of tall fescue ESTs representing different abiotic stresses, tissue types and developmental stages. BMC Plant Biol. 2008, 8: 27-10.1186/1471-2229-8-27.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee MC, Lee YJ, Lee MH, Nam HG, Cho TJ, Hahn TR, Cho MJ, Sohn U: Large-scale analysis of expressed genes from the leaf of oilseed rape (Brassica napus L.). Plant Cell Rep. 1998, 17: 930-936. 10.1007/s002990050512.View ArticleGoogle Scholar
- Ablett E, Seaton G, Scott K, Shelton D, Graham MW, Baver-stock P, Lee LS, Henry R: Analysis of grape ESTs: global gene expression patterns in leaf and berry. Plant Sci. 2000, 159: 87-95. 10.1016/S0168-9452(00)00335-6.PubMedView ArticleGoogle Scholar
- Yu I, Fengler KA, Clough SJ, Bent AF: Identification of Arabidopsis mutants exhibiting an altered hypersensitive response in gene-for-gene disease resistance. Microbe Interact. 2000, 13: 277-286. 10.1094/MPMI.2000.13.3.277.View ArticleGoogle Scholar
- Piffanelli P, Zhou R, Casais C, Orme J, Jarosch B, Schaffrath U, Collins NC, Panstruga R, Schulze-Lefeert P: The barley MLO modulator of defence and cell death is responsive to biotic and abiotic stress stimuli. Plant Physiol. 2002, 129: 1076-1085. 10.1104/pp.010954.PubMedPubMed CentralView ArticleGoogle Scholar
- de las Mercedes Dana M, Pintor-Toro JA, Cubero B: Transgenic tobacco plants over expressing chitinases of fungal origin show enhanced resistance to biotic and abiotic stress agents. Plant Physiol. 2006, 142: 722-730. 10.1104/pp.106.086140.PubMedView ArticleGoogle Scholar
- Ferreira RB, Monteiro S, Freitas R, Santos CN, Chen Z, Batista LM, Duarte J, Borges A, Teixeira AR: The role of plant defence proteins in fungal pathogenesis. Mol Plant Pathol. 2007, 8: 677-700. 10.1111/j.1364-3703.2007.00419.x.PubMedView ArticleGoogle Scholar
- Mauch-Mani B, Mauch F: The role of abscisic acid in plant-pathogen interactions. Curr Opin Biotech. 2005, 8: 409-414. 10.1016/j.pbi.2005.05.015.Google Scholar
- Mayek-Perez NO, Garcia-Espinosa R, Lopez-Castaneda C, Acosta-Gallegos JA, Simpson J: Water relations, histopathology and growth of common bean (Phaseolus vulgaris L.) during pathogenesis of macrophomina phaseolina under drought stress. Physiol Mol Plant Pathol. 2002, 60: 185-195. 10.1006/pmpp.2001.0388.View ArticleGoogle Scholar
- Koga H, Dohi K, Mori M: Abscisic acid and low temperatures suppress the whole plant-specific resistance reaction of rice plants to the infection with Magnaporthe grisea. Physiol Mol Plant Pathol. 2004, 65: 3-9. 10.1016/j.pmpp.2004.11.002.View ArticleGoogle Scholar
- Whenham RJ, Fraser RSS, Brown LP, Payne JA: Tobacco mosaic virus-induced increases in abscisic acid concentration in tobacco leaves: intracellular location in light and darkgreen areas, and relationship to symptom development. Planta. 1986, 168: 592-598. 10.1007/BF00392281.PubMedView ArticleGoogle Scholar
- Jennings DB, Ehrenshaft M, Pharr DM, Williamson JD: Roles for mannitol and mannitol dehydrogenase in active oxygen-mediated plant defense. Proc Natl Acad Sci USA. 1998, 95: 15129-15133. 10.1073/pnas.95.25.15129.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu CJ, Deavours BE, Richard SB, Ferrer JL, Blount JW, Huhman D, Dixon RA, Noel JP: Structural basis for dual functionality of isoflavonoid O-methyltransferases in the evolution of plant defense responses. Plant Cell. 2006, 18: 3656-3669. 10.1105/tpc.106.041376.PubMedPubMed CentralView ArticleGoogle Scholar
- Day DA, Millar AH, Whelan J: Plant mitochondria: from genome to function. Springer publishers, The Netherlands 2004.View ArticleGoogle Scholar
- Scheltinga ATV, Kalk KH, Beintema JJ, Dijkstra B: Crystal structures of Hevamine, a plant defence protein with chitinases and lysozymes activity, and its complex with an inhibitor. Structure. 1994, 2: 1181-1189. 10.1016/S0969-2126(94)00120-0.View ArticleGoogle Scholar
- Grenier J, Potvin C, Asselin A: Barley pathogenesis-related proteins with fungal cell wall lytic activity inhibit the growth of yeasts. Plant Physiol. 1993, 103: 1277-1283. 10.1104/pp.103.4.1277.PubMedPubMed CentralView ArticleGoogle Scholar
- De DN: Plant cell vacuoles.Collingwood, Australia: CSIRO Publishing 2000, 288.Google Scholar
- Mehta PA, Sivaprakash K, Parani M, Venkataraman G, Paida AK: Generation and analysis of expressed sequence tags from the salt-tolerant mangrove species Avicennia marina (Forsk) Vierh. Theor Appl Genet. 2005, 100: 416-424. 10.1007/s00122-004-1801-y.View ArticleGoogle Scholar
- Kota R, Varshney RK, Thiel T, Dehmer KJ, Graner A: Generation and comparison of EST-derived SSRs and SNPs in barley (Hordeum vulagre L.). Hereditas. 2001, 135: 145-151. 10.1111/j.1601-5223.2001.00145.x.PubMedView ArticleGoogle Scholar
- Gupta PK, Varshney RK: The development and use of microsatellite markers for genetic analysis and plant breeding with emphasis on bread wheat. Euphytica. 2000, 113: 163-185. 10.1023/A:1003910819967.View ArticleGoogle Scholar
- Shultz JL, Kazi S, Bashir R, Afzal JA, Lightfoot DA: The development of BAC-end sequence-based microsatellite markers and placement in the physical and genetic maps of soybean. Theor Appl Genet. 2007, 114: 1081-1090. 10.1007/s00122-007-0501-9.PubMedView ArticleGoogle Scholar
- Varshney RK, Thiel T, Stein N, Langridge P, Graner A: In silico analysis of frequency and distribution of microsatellites in ESTs of some cereal species. Cell Mol Biol Lett. 2002, 7: 537-546.PubMedGoogle Scholar
- Luo M, Dang P, Guo BZ, He G, Holbrook CC, Bausher MG, Lee RD: Generation of expressed sequence tags (ESTs) for gene discovery and marker development in cultivated peanut. Crop Sci. 2005, 45: 346-353.View ArticleGoogle Scholar
- Yu JK, Sun Q, Rota ML, Edwards H, Tefera H, Sorrells ME: Expressed sequence tag analysis in tef (Eragrostis tef (Zucc) Trotter). Genome. 2006, 49: 365-372. 10.1139/G05-118.PubMedView ArticleGoogle Scholar
- Quilang J, Wang S, Li P, Abernathy J, Peatman E, Wang Y, Wang L, Shi Y, Wallace R, Guo X, Liu Z: Generation and analysis of ESTs from the eastern oyster, Crassostrea virginica Gmelin and identification of microsatellite and SNP markers. BMC Genomics. 2007, 8: 157-10.1186/1471-2164-8-157.PubMedPubMed CentralView ArticleGoogle Scholar
- Scott KD, Eggler P, Seaton G, Rossetto M, Ablett EM, Lee LS, Henry RJ: Analysis of SSRs derived from grape ESTs. Theor Appl Genet. 2000, 100: 723-726. 10.1007/s001220051344.View ArticleGoogle Scholar
- Sivaramakrishnan S, Seetha K, Reddy LJ: Diversity in selected wild and cultivated species of pigeonpea using RFLP of mtDNA. Euphytica. 2002, 125: 21-28. 10.1023/A:1015759318497.View ArticleGoogle Scholar
- Kota R, Varshney RK, Prasad M, Zhang H, Stein N, Graner A: EST-derived single nucleotide polymorphism markers for assembling genetic and physical maps of the barley genome. Funct Integr Genomic. 2008, 8: 223-233. 10.1007/s10142-007-0060-9.View ArticleGoogle Scholar
- Varshney RK, Dubey A: Novel genomic tools and modern genetic and breeding approaches for crop improvement. J Plant Biochem Biotechnol. 2009, 18: 127-138.View ArticleGoogle Scholar
- Rafalski A: Applications of single nucleotide polymorphisms in crop genetics. Curr Opin Plant Biol. 2002, 5: 94-100. 10.1016/S1369-5266(02)00240-6.PubMedView ArticleGoogle Scholar
- Moreno-Vazquez S, Ochoa OE, Faber N, Chao S, Jacobs JM, Maison-Neuve B, Kesseli RV, Michelimore RW: SNP-based codominant markers for a recessive gene conferring resistance to corky root rot (Rhizomonas suberifaciens) in lettuce (Lactuca sativa). Genome. 2003, 46: 1059-1069. 10.1139/g03-073.PubMedView ArticleGoogle Scholar
- Varshney RK, Nayak SN, May GD, Jackson SA: Next-generation sequencing technologies and their implications for crop genetics and breeding. Trends Biotechnol. 2009, 27: 522-530. 10.1016/j.tibtech.2009.05.006.PubMedView ArticleGoogle Scholar
- Wang S, Sha Z, Sonstegard TS, Liu H, Xu P, Somridhivej B, Peatman E, Kucuktas H, Liu Z: Quality assessment parameters for EST-derived SNPs from catfish. BMC Genomics. 2008, 9: 450-10.1186/1471-2164-9-450.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmitt ME, Brown TA, Trumpower BL: A rapid and simple method for preparation of RNA from Saccharomyces cerevisiae. Nucl Acids Res. 1990, 18: 3091-3092. 10.1093/nar/18.10.3091.PubMedPubMed CentralView ArticleGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular cloning. Cold Spring Harbor Laboratory Press, Plainview, 1989, I-III.Google Scholar
- EST trimmer. [http://pgrc.ipk-gatersleben.de/misa/download/est_trimmer.pl]
- Huang X, Madan A: CAP3: a DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMedPubMed CentralView ArticleGoogle Scholar
- NCBI EST database. [http://www.ncbi.nlm.nih.gov/dbEST]
- IDEG.6 analysis tool. [http://telethon.bio.unipd.it/bioinfo/IDEG6_form]
- Romualdi C, Bortoluzzi S, D'alessi F, Danieli GA: IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics. 2003, 12: 159-162.PubMedView ArticleGoogle Scholar
- Hierarchical Clustering Explorer 2.0. [http://www.cs.umd.edu/hcil/hce/hce2.html]
- SSR identification tool. [http://pgrc.ipk-gatersleben.de/misa/]
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Humana Press, Totowa, NJ, 365-386. Source code, [http://fokker.wi.mit.edu/primer3/]Google Scholar
- Nei M, Li WH: Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc Natl Acad Sci USA. 1979, 76: 5269-5273. 10.1073/pnas.76.10.5269.PubMedPubMed CentralView ArticleGoogle Scholar
- Rohlf FJ: NTSYS-pc. Numerical taxonomy and multivariate analysis system, Version 2.10. Exeter Software, New York 2002.Google Scholar
- Jayashree B, Hanspal MS, Srinivasan R, Vigneshwaran R, Varshney RK, Sputhi N, Eshwar K, Ramesh N, Chandra S, Hoisngton DA: An integrated pipeline of open source software adapted for multi-CPU architectures: use in the large-scale identification of single nucleotide polymorphisms. Comp Funct Genomics. 2007, Article ID 35604-Google Scholar
- Jayashree B, Bhanuprakash A, Jami A, Reddy SP, Nayak S, Varshney RK: Perl module and PISE wrappers for the integrated analysis of sequence data and SNP features. BMC Res Notes. 2009, 2: 92-10.1186/1756-0500-2-92.PubMedPubMed CentralView ArticleGoogle Scholar
- Restriction enzyme data base. [http://rebase.neb.com/]
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.