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

Generation and analysis of expressed sequence tags (ESTs) for marker development in yam (Dioscorea alata L.)

BMC Genomics volume 12, Article number: 100 (2011) Cite this article

Abstract

Background

Anthracnose (Colletotrichum gloeosporioides) is a major limiting factor in the production of yam (Dioscorea spp.) worldwide. Availability of high quality sequence information is necessary for designing molecular markers associated with resistance. However, very limited sequence information pertaining to yam is available at public genome databases. Therefore, this collaborative project was developed for genetic improvement and germplasm characterization of yams using molecular markers. The current investigation is focused on studying gene expression, by large scale generation of ESTs, from one susceptible (TDa 95-0310) and two resistant yam genotypes (TDa 87-01091, TDa 95-0328) challenged with the fungus. Total RNA was isolated from young leaves of resistant and susceptible genotypes and cDNA libraries were sequenced using Roche 454 technology.

Results

A total of 44,757 EST sequences were generated from the cDNA libraries of the resistant and susceptible genotypes. Greater than 56% of ESTs were annotated using MapMan Mercator tool and Blast2GO search tools. Gene annotations were used to characterize the transcriptome in yam and also perform a differential gene expression analysis between the resistant and susceptible EST datasets. Mining for SSRs in the ESTs revealed 1702 unique sequences containing SSRs and 1705 SSR markers were designed using those sequences.

Conclusion

We have developed a comprehensive annotated transcriptome data set in yam to enrich the EST information in public databases. cDNA libraries were constructed from anthracnose fungus challenged leaf tissues for transcriptome characterization, and differential gene expression analysis. Thus, it helped in identifying unique transcripts in each library for disease resistance. These EST resources provide the basis for future microarray development, marker validation, genetic linkage mapping and QTL analysis in Dioscorea species.

Background

Yams (Dioscorea spp.) are the primary agricultural commodities and major staple crop in Africa. Yam tubers are nutritionally rich and a major source of dietary fiber, carbohydrates, vitamin C and essential minerals. Worldwide, yam consumption is 18 million tons http://www.IITA.org. In 2007, yam production was 52 million tons worldwide, of which Africa produced 96%, and Nigeria is the major producer (71%) with more than 37 million tons [1]. The consumer demand for yam is very high in sub-Saharan region of Africa, but the yam production is declining in this region due to factors including anthracnose disease caused by a fungus, Colletotrichum gloeosporioides[2], pests, and decline in soil fertility [3]. Yams are polyploid crop species and are propagated vegetatively from tubers (whole or setts).

The water yam (D. alata) is the most widely cultivated species and is highly susceptible to anthracnose disease [4, 5]. The genetic improvement of yam at IITA and CTCRI (India) concentrated on the development of disease resistant and high yielding varieties. Through classical breeding, it would be very difficult to develop a resistant cultivar due to constraints such as the long growth cycle (8-10 months), dioecious and poor flowering nature, polyploidy, vegetative propagation and heterozygous genetic background [6]. A large collection of germplasm with huge genetic variability is available, and it would be profitable to use candidate gene approach and trap the important trait information for disease and pest resistance.

Colletotrichum is a large genus of ascomycete fungi, containing many species which cause anthracnose or blight on a wide range of important crops and ornamental plants [7]. Two genes, clk1 which encodes a serine/threonine kinase in C. lindemuthianum[8], and cap20 which encodes a wall glycoprotein of C. gloeosporioides appressoria[9], have been shown to have a role in pathogenicity and virulence, respectively. Necrotrophy is clearly linked to the increased expression of plant cell wall degrading enzymes such as endo-polygalacturonases (endo- PG) and pectin lyases. Previous studies [10] showed that endo-PG, two forms of pectin lyase, a- and b-galactopyranosidase, a-arabinofuranosidase, and a protease are secreted into culture medium containing polypectate of bean cell walls. Pectin lyase activity was first observed 4 days after inoculation of beans with C. lindemuthianum, rising to maximum activity at 7 days, after which activity declined. Thus, the expression of pectin lyase activity correlates well with the onset of necrotrophy and the subsequent development of lesions [11]. Therefore the yam leaf material for current investigation was collected on 3rd and 7th day after inoculation with C. gloeosporioides fungus.

Initial genetic inheritance studies showed that resistance to yam anthracnose in water yam is dominant and quantitatively inherited [12]. A single major dominant locus controlling resistance in the breeding line TDa 95-0328 was tentatively designated Dcg-1, until allelism was investigated [6]. Lower number of molecular markers in yam limits the genetic mapping efforts. Previously developed SSRs from Dioscorea tokoro[13], a wild diploid yam species of East Asia and Japan were found not to be useful in cultivated yams D. rotundata and D. alata[4]. Most of the current molecular markers for the yam genome by Mignouna et al. are based on AFLP and RAPD [14, 15]. However, for effective gene discovery and marker-assisted breeding, it is important to develop more user-friendly, efficient, transportable and co-dominant markers such as simple sequence repeat (SSR) markers and single nulceotide polymorphism markers (SNPs). EST analysis is the most efficient and effective approach for the identification of candidate genes and also assist in new molecular markers such as EST-SSRs and SNPs [16].

As of 10th October 2010, there were only 31 EST sequences stored at GenBank (National Center for Biotechnology and Information, http://www.ncbi.nlm.nih.gov/) for the genus Dioscorea. To understand the transcriptome of yam and for new marker discovery to assist in yam crop improvement programs, generation of ESTs would be extremely desirable. With the objective of characterizing transcriptome in yam, initial efforts were focused on successful isolation of mRNA, cDNA library construction and sequencing using Sanger's method [17]. Two hundred EST sequences were generated using this approach, but reported little functional significance during BLAST search analysis. To enrich and characterize expressed gene sequence information especially during the anthracnose infection for candidate gene identification and also to utilize the sequence information for new marker discovery, a collaborative effort to generate large number of ESTs was initiated between VSU, USA and IITA, Africa. For this purpose, two anthracnose resistant germplasm lines 87-01091, 95-00328 and a susceptible line 95-0310 were selected for cDNA library construction from leaf tissues inoculated with the C. gloeosporioides fungus. The Roche 454 pyrosequencing technology was used to generate EST sequence information (Agencourt Biosciences, MA).

Results and Discussion

Generation of ESTs

Assembly process generated 15,196 ESTs in TDa 95-0328; 15,984 ESTs in TDa 95-0310, and 13,577 ESTs in TDa 87-01091 with average sequence lengths of 426, 411 and 524 bases, respectively (Table 1). All the EST sequences generated were submitted to GenBank at NCBI and were assigned Genbank Accession numbers HO809681-HO825421 with dbEST id from 71421255 to 71436995 for library TDa 95-0310; HO825422-HO840419 with dbEST ids from 71436996 to 71451993 for library TDa 95-0328; HO850622 to HO864016 with dbEST ids from 71462196 to 71475590 for library TDa 87-01091. These sequences contributed to >99% of ESTs currently available in yam at Genbank and would serve as the major resource for the yam research community.

Table 1 List of Yam Genotypes used, average EST length and total number of ESTs generated in the present study.
Full size table

Transcriptome Analysis

EST sequences obtained from each individual library were first analyzed using Blast2GO http://blast2go.org[18]. Blast2GO provides the gene annotation information for the ESTs and also helps in comparison of two EST datasets to find expression differences. BLAST search analysis across all three libraries revealed highest homology with the sequences from Oryza sativa, Vitis vinifera, Populus trichocarpa and other crop species. Mapping and annotation steps in Blast2GO program resulted in providing GO IDs (Gene Ontology identifiers) to individual ESTs in each dataset of the three libraries (Additional file 1). A total of 9,133 (57.86%) of 15,787 ESTs were annotated with 4,074 non-redundant (nr) GO IDs in the susceptible TDa 95-0310 germplasm line, while in resistant TDa 95-0328 line, 8,985 (59.75%) of 15,038 ESTs were annotated with 3,839 nr GO IDs and in other resistant line TDa 87-01091 library, 4,577 (34.08%) of 13,429 ESTs were annotated with 2,465 nr GO IDs. Approximately 50% of ESTs could not be annotated and remained unknown using Blas2GO. To complement and comprehensively annotate the Yam transcriptome, all the EST sequences were submitted to Mercator tool of MapMan database. Mercator tool generated functional predictions by searching 6 reference databases (3 BLAST based and 2 based on reverse position-specific BLAST and InterProScan) and assigned MapMan Bin Ids (additional file 2). Still 44% of ESTs lacked significant homology with genes in other crops in each of the three datasets and these may include novel metabolic genes in yam. A comparative analysis was done using the Bin ID information from each library to identify unique ESTs in each library and also the common genes shared with other libraries. The distribution of ESTs with known function is presented as a Venn diagram in Figure 1. A large subset of 6295 ESTs representing 660 unique Bin IDs were shared across three libraries and these may account for housekeeping genes involved in general cellular metabolism. The two resistant lines TDa 95-0328 and TDa 87-01091 had 115 and 180 unique ESTs, respectively, which may account for the tolerance against C. gloeosporioides fungus.

Figure 1
figure 1

Venn Diagram describing distribution of ESTs across three libraries, TDa 95-0310, TDa 95-0328, and TDa 87-01091.

Full size image

In TDa 95-0310 susceptible genotype, there were 71 unique sequences which predominantly matched genes for protein synthesis and nucleotide synthesis (Figure 2). However, there were several sequences annotated as protein kinases and were known to have a definite role in disease signaling [19]. Another interesting feature is the presence of genes for synthesis of cell wall precursors and secondary metabolites which might be highly induced due to severe biotic stress. There were also genes for recycling NAD to maintain glycolysis and substrate level phosphorylation in the absence of oxygen in damaged leaf of susceptible genotype [20]. Removal of the toxic compounds like proline is very essential for plant survival under stress [21]. There were unique genes for proline degradation in susceptible genotype suggesting the presence of a possible mechanism of detoxification during pathogenesis.

Figure 2
figure 2

Distribution of unique ESTs into MapMan functional categories in the three yam genotypes.

Full size image

Common sequences between TDa 95-0310 and TDa 95-0328, were annotated as 122 genes representing cell wall precursors, co-factor and vitamin metabolism, tetrapyrrole synthesis, C1-metabolism, lipid metabolism, amino acid metabolism, secondary metabolism, hormone metabolism, major and minor CHO metabolism, nucleotide metabolism, genes for regulation of transcription, amino acid activation, protein synthesis, protein degradation, energy metabolism, and genes for sugar and light signaling. On the other hand, the common genes between TDa 95-0328 and TDa 87-01091 were cell wall proteins, sulfate assimilation related genes, phosphatidylserine decarboxylase, hydroxymethylpyrimidine kinase, L-Galactono-1,4-lactone dehydrogenase, genes for dicer degradation, GeBP like transcription factor, and myo inositol oxygenase. Similarly, TDa 87-01091 and TDa 95-0310 shared sequences related to dihydroneopterin aldolase, Glycine cleavage H protein, phosphomannose isomerase, C-lectin, alcohol dehydrogenase, and malate dehydrogenase genes. Although these genes have specific role in metabolism, regulation and signaling, the significance of their sharing between yam genotypes cannot be speculated and is beyond the scope of this paper.

Functional Analysis

Using the GO ID information, the ESTs were classified into subsets based on molecular function and is represented as a pie diagram for each library (Figure 3). A major subset of ESTs (> 60%) across all libraries was linked to binding and catalytic activity, where the remaining groups involved transporter, signal transduction, transcription factors, secondary metabolism and antioxidant activities. Based on biological process classification, a major proportion of the genes were involved in oxidation-reduction activity across all libraries and these may involve genes in energy metabolisms and other housekeeping activities. ESTs were also grouped under other biological processes such as photosynthesis, membrane transport, transcription, translation, protein folding, and others. However, a significant proportion of the expressed genes (~5-10%) were biotic stress and other stress responsive genes (Figure 3).

Figure 3
figure 3

GO distribution in TDa 95-0310, TDa 95-0328 and TDa 87-01091 based on molecular function.

Full size image

Differential gene expression analysis

From the EST datasets, it is interesting to identify genes involved in defense mechanisms that were differentially expressed in susceptible and resistant yam genotypes. Using the Bin Ids of ESTs in each library, the expression of various ESTs were compared from one library to the other. These results were summarized in Figure 2 and additional file 3. The resistant genotypes had unique ESTs related to carbohydrate metabolism, cell wall biogenesis, lipid and amino acid metabolism, secondary and hormone metabolism, transcription factors, protein synthesis, and signaling proteins (Figure 2). Even though, the EST datasets is not comprehensive of the entire yam transcriptome, this analysis would shed some light on the gene expression differences during anthracnose infection. Further investigation of these differentially expressed genes using microarray platforms or RNA-seq or real-time quantitative PCR experiments, would assist in elucidating the underlying molecular mechanism for anthracnose tolerance in yam.

Informative defense related genes during anthracnose infection

The differential gene expression analysis provided the gene information that might be potentially over-expressed or under-expressed compared to other libraries. These results were summated in additional file 3.

Multiple pathogenesis-related genes and host defense related genes including signaling genes were identified from the ESTs in the resistant lines and they include nicotiana lesion inducing proteins, erwinia induced proteins, cysteine proteases, cysteine protease inhibitors, peroxidases, extensin, wall associated kinases, brassinosteroid receptor kinases, ethylene responsive genes, phosphatidylinositol 4-kinase, leucine rich repeat (LRR) genes, serine-threonine protein kinases and others signaling proteins. Signal transduction in plants, leading to the expression of defense genes, is initiated by leucine rich repeat (LRR)-type membrane receptors containing intrinsic kinase activity as in the case of FLS2 [22].

There were unique genes involved in isoprenoid and flavonoids biosynthesis and fatty acid synthesis and elongase activity in TDa 87-01091 compared to TDa 95-0310 and TDa 95-0328. There were hits specifically showing homologies to peroxisome biogenesis protein pex1 (GO: 0009851, GO: 0017111) and perxisome assembly factor (GO: 0016558) in library 87 and 328. This confirms previous reports of plant peroxisomes role in the biosynthesis of the signaling molecules like jasmonic acid, β-oxidation of indole butyric acid (IBA), and sulphur and polyamine metabolism. Moreover, evidence is emerging from recent studies that peroxisomes have important functions in specific defense mechanisms, conferring resistance against pathogen attack [23].

There were homologies to 12-oxo-phytodienoic acid reductases (OPRs), enzymes of the octadecanoid pathway, which convert linolenic acid to a phytohormone jasmonic acid. Fifty nine significant matches to auxin responsive proteins (GO: 0006417, 9725, 9734, 46983), auxin signaling f-box 3 (GO: 0002237) and auxin response transcription factors (GO: 0045449) were observed in TDa 87-01091 and TDa 95-0328, confirming their role in host defense against pathogens.

Actually Plant peroxidases (POXs) transduce the extracellular signals into the redox signals that eventually stimulate the intracellular Ca2+ signaling required for induction of defense responses [24]. There were significant hits to different plant peroxidases in resistant libraries confirming their requirement for defense against Colletotrichum fungus.

The hits to pectate lyase were few, one each in TDa 95-0328 (GO: 30570) and TDa 87-01091 (GO: 0016829) compared to four in TDa 95/0310. The contribution of a single Pectate lyase gene (pel) to the pathogenic abilities of C. magna confirmed that Pectate lyase is a pathogenicity factor required for the penetration and colonization of Colletotrichum species [25]. Cell-wall-degrading enzymes (CWDEs), like pectolytic enzymes, are considered to play a role in the pathogenesis of bacteria and fungi on their hosts [26]. The peptidoglycan-binding LysM domain-containing protein and beta-1,4-glucanases may help in fungal cell wall degradation [27] as there were significant hits observed in the resistant genotypes. The disruption of a single enzyme may be complemented by the activity of other CWDEs [28]. There was also homology to chitin binding, chitinase activity and chitin catabolic/metabolic process and were abundant in 87-01091, a genotype resistant to C. gloeosporioides strain SGG.

The activation of the defense response, after pathogen inoculation, is the most important factor for the success of host resistance [29]. The EST data from our study also suggest strength of the expression of ubiquitous early acting genes for establishment of the resistance response of the resistant genotype 87-01091.

Marker discovery

The EST sequence information generated in this study also serves as a major resource to generate new molecular markers. Simple sequence repeat motifs detected in EST sequences are usually termed as EST-SSRs. These are very informative markers as they exist within the gene sequence. A total of 44,254 ESTs from three genotypes, with an average length of 500 bp were used to evaluate for the presence of SSR motifs. There are a total of 1702 EST sequences containing SSRs identified from three libraries. In this work, SSRs were considered for primer design that fitted the following criteria: a minimum pattern length of 12 bp, excluding polyA and polyT repeat, at least 7 repeat units in case of di-nucleotide and at least 5 repeat units for tri-, tetra-, penta- and hexa-nucleotide SSRs. There are about a total of 3459 good quality repeats containing di (289), tri (1271), tetra (66), penta (19), hexa (53) and hepta (3) repeats which were used for designing primers (Table 1). The default settings for Primer3 program: optimum temperature of 63oC and an optimum primer size of 24 bases were selected for designing the primers. A total of 1705 primer pairs were designed from 1702 EST sequences. A subset were tested for screening genomic DNA of 3 genotypes using simple 1.2% agarose gel electrophoresis for polymorphism detection and found polymorphic. These SSR primers are summarized in additional file 4.

Identification of homologous SNPs is very challenging in a polyploid species such as yam. HaploSNPer generates putative SNP haplotype for each contig using CAP3 generated .ace file. Therefore, we can verify whether a SNP is heterozygous within each genotype or homozygous at specific loci. If they are heterozygous within each genotype, basically they may be copies from homeologous genomes or paralogous copies. Even though multiple sequence clusters were identified with single nucleotide variations, most of them were polymorphic within a germplasm line at many loci. Usually they will be heterozygous always and exhibit no polymorphism between the genotypes. Selection of such SNPs usually results in false positive assays. However, we were able to detect 104 candidate SNPs between libraries 310 and 328 that are homologous within each genotype. These SNPs can be used for further genotyping and generating genetic maps in yam.

Conclusion

Despite the economic importance of yam, very little genomic or transcriptome sequence data is publicly available. This is the first large scale EST generation attempt made in yam with the objective of providing a comprehensive annotated transcriptome dataset that will be publicly available, and also used as a resource for novel gene and new marker discovery for crop improvement in yam. We have successfully sequenced transcriptomes of two resistant and one susceptible lines of yam under the conditions of anthracnose infection. Approximately 56% of total ESTs were annotated and analyzed for functional characterization and differential expression of genes for tolerance to anthracnose disease. We also used this dataset as a resource to design new SSR and SNP markers. They serve as a useful tool in identifying genetic variation in the current cultivars, wild relatives and also assist in generating genetic linkage maps. The SSR markers generated are currently being evaluated at IITA for use in their yam improvement program. The identified SNP markers will be validated to identify QTL regions through the ongoing research on yam at IITA, USDA-ARS and VSU. The markers generated could be used on F1 population generated at IITA to develop first generation EST-SSR based QTL map for yam.

Methods

The D. alata genotypes were selected based on their scoring (0-6 scale) to leaf damage by [30] anthracnose for resistance (0-17.5% mean leaf area damage). The fully expanded, pathogen challenged, young leaves from yam genotypes with differential resistance/susceptibility to FGS (fast-growing salmon) and SGG (slow-growing grey) strains of anthracnose (Table 1) were harvested at IITA. The pathogen inoculations were done at IITA following the standardized inoculation procedures [31] similar to the method described by Green et al. [32]. A spore suspension of C. gleosporoides in sterile distilled water was prepared from 7-10 days old single-spore cultures (based on virulence on D. alata) and standardised to 106 spores per ml with tween 80 (Merck)(1.2% v/v) added as wetting agent. Young leaves of both the genotypes were inoculated on both sides and the green leaf tissues next to the infestation was collected on 4th day and 7th day after infestation.

These leaf samples were freeze dried and supplied to VSU for genomic studies in 2007. Extraction of high quality RNA is very challenging from yam tissues due to high polyphenolic content and other residues. Total RNA was isolated from these leaf samples using the optimized RNA isolation protocol [17] and was quantified using Spectrophotometer (Bio-Rad Laboratoires, CA, USA). The quality was further confirmed by running a 1.2% Formaldehyde Agarose Gel Electrophoresis. The cDNA library was constructed using Clontech's Creator SMART cDNA library construction kit with the pDNR-lib vector (Clonetech Laboratories, CA, USA). The total RNA was reverse transcribed to cDNA using Powerscript Reverse Transcriptase, using kit primers SMART IV Oligonucleotide and CDS III/3' PCR primer. The cDNA was PCR-amplified using the Advantage 2 PCR kit, using the SMART 5' PCR III primer and CDS III/3' PCR primer, using between 18 and 26 cycles according to the manufacturer's instructions. The libraries from all the genotypes were arrayed in 96 well plates and stored at -80°C.

To assess the quality of the library, electro competent cells, DH10B E. coli bacterial cells (T1 phage resistant) were transformed with the vector containing cDNA inserts (Gene pulser, Bio-Rad Laboratoires, CA, USA) grown in liquid suspension for two hours and then plated onto LB agar plates with 30 ug/ml chloramphenicol. The independent colonies were collected separately and arrayed into 96 well plates for Sanger sequencing. Preliminary sequencing efforts revealed very little information about disease resistance due to less number of hits to functional genes; we focused our efforts to next generation sequencing methods. The cDNA library was constructed for two resistant and one susceptible line. The cDNA quality was tested by random cloning of the sequences using creator SMART cDNA library construction kit (Clonetech Laboratories, CA) and individual clones were sequenced using Sanger sequencing method. Agencourt Bioscience Corporation (MA, USA) was identified as a preferred partner to sequence the cDNA libraries from pathogen challenged leaf tissues using Roche 454 pyrosequencing technologies. High quality fragment cDNA library was constructed and 454 sequence data was generated by Agencourt Biosciences (MA, USA). Raw sequence data generated from three cDNA libraries using 454 GS FLX (Titanium chemistry) instrument was assembled into contigs to generate EST dataset for each library (Agencourt Biosciences, MA). Sequencing of these libraries from the challenged leaves of resistant and susceptible genotypes generated a total of 195,937 (TDa 87-01091), 217,868(TDa 95-00328) and 399,987(TDa 95-0310) cDNA raw reads, respectively.

Sequence assembly, analysis and annotation

  1. a.

    Preliminary sequencing data analysis: The Sanger sequences were analyzed using NCBI BLAST tool. The vector sequences from each sequence were trimmed manually, Sequences less than 150 base pairs were removed and then ESTs were identified based on NCBI "blastn" similarities, Blast hits considered significant cover at least 150 bases with 80-96% identity. All ESTs sharing significant similarities were clustered together. ESTs with no significant similarities to any other ESTs were given their own ID number and referred as singletons.

  2. b.

    454 sequence data analysis: The sequencing agency provided high quality raw reads as well as aligned EST sequences for three libraries. These sequences were analyzed using Blast2GO Program. BLAST hits with an e-value of 10-5 or less, which corresponds approximately to a 60-bp contiguous perfect match in the data set, were considered to be successful hits against the transcriptome. The aligned EST Sequences with BlastX hits were mapped and annotated according to gene ontology terms (GO) using the program Blast2Go [33]. The distribution of genes in each ontology categories was examined and the percentages of unique sequences in each of the assigned GO term namely, biological process, molecular function, and cellular component, were computed and presented. Using the available annotated EST information in each library, a Fisher exact test was performed at p-value 0.05 in Blast2GO program to compare the expression level of various ESTs from one library to the other [18].

  3. c.

    Mapman analysis: For data analysis by MapMan software version 3.5.0 BETA [34] http://gabi.rzpd.de/projects/MapMan/, the EST sequences from yam libraries were uploaded into Mercator tools and the mapping file was generated. Functional predictions and Bin Ids were generated by searching a variety of reference databases (currently 6 are available: 3 BLAST-based, 2 reverse position-specific BLAST based and InterProScan) and subsequently evaluating and compiling the search results for each input gene to propose a functional Bin based on the manually curated binning of the reference database entries. List of ESTs represented in the mapping file and their corresponding Bin identifiers are listed in additional file 2.

  4. d.

    BLAST search against Colletotrichum genes: The total yam EST sequences were blasted against 962 nucleotide sequences available at NCBI for C. gloeosporioides (as of October 17, 2010). There were not any significant hits from the yam ESTs (e-value <10-15) except a few that belong to 18 s ribosomal unit of Colletotrichum and the rest being the calmodulin, actin and EF1-alpha genes. The alignment length of these hits is also less than 50 bp.

  5. e.

    BLAST search against Plant Vs fungal gene data base: We find very few significant similarities with the fungal genes 1,20,8, respectively in TDa 95-0310, TDa 95-0328 and TDa 87-01091 with their alignment length ranging from 637-917 at eValue = 0 and BIT score more than 100. There were 12,296 genes unique to yam not showing hits to either plant or fungal genes while 31,561 ESTs showing similarities to plant genes. These 29 ESTS were eliminated for final EST analysis using Mapman.

Bioinformatics mining of microsatellites

The total ESTs more than 150 base pairs were searched for microsatellites using Schroeder's SSR finder software [35]. In this work, SSRs were considered for primer design that fitted the following criteria: a minimum pattern length of 12 bp, excluding polyA and polyT repeat, at least 7 repeat units in case of di-nucleotide and at least 5 repeat units for tri-, tetra-, penta- and hexa-nucleotide SSRs. The default settings for Primer3 input were optimum temperature of 63oC and an optimum primer size of 24 bases.

Identification of SNPs

EST sequences were pooled based on the GO ID information captured during annotation process using a relational database. These homo or orthologous sequences were submitted to the web-based HaploSNPar program to generate haplotype clusters and identify candidate SNPs with default parameters [35].

Abbreviations

RNA:

Ribonucleic Acid

DNA:

Deoxy Ribonucleic Acid

cDNA:

Complementary DNA

BLAST:

Basic Local Alignment Search Tool

bp:

base pairs

EST:

expressed sequence tag

GO:

gene ontology

SSR:

Simple Sequence Repeat

SNP:

Single Nucleotide Polymorphism

References

  1. FAO 2008. accessed on September 3, 2010, [http://www.fao.org]

  2. Abang MM, Winter S, Green KR, Hoffmann P, Mignouna HD, Wolf GA: Molecular identification of Colletotrichum gloeosporioides causing anthracnose of yam in Nigeria. Plant Pathol. 2002, 51: 63-71. 10.1046/j.0032-0862.2001.00655.x.

    Article  Google Scholar 

  3. Brunt AA, Jackson GVH, Frison EA: FAO/IBPGR Technical Guidelines for the Safe Movement of Yam Germplasm. Food and Agriculture Organization of the United Nations. 1989, Rome/International Board for Plant Genetic Resources, Rome

    Google Scholar 

  4. McDonald FD, Alleyne AT, Ogarro LW, Delauney AJ: Yam anthracnose in the English-speaking islands of the Eastern Caribbean-successes and research advances in disease management. Trop Agric. 1998, 75: 53-57.

    Google Scholar 

  5. Winch JE, Newhook FJ, Jackson GVH, Cole JS: Studies of Colletotrichum gloeosporioides disease on yam, Dioscorea alata, in Solomon Islands. Plant Pathol. 1984, 33: 467-477. 10.1111/j.1365-3059.1984.tb02870.x.

    Article  Google Scholar 

  6. Mignouna HD, Abang MM, Asiedu R: Harnessing modern biotechnology for tropical tuber crop improvement: Yam (Dioscorea spp.) molecular breeding. African Journal of Biotechnology. 2003, 2 (12): 478-485.

    Article  CAS  Google Scholar 

  7. Bailey JA, Jeger MJ: Colletotrichum: Biology, Pathology and Control. 1992, CAB International, Wallingford

    Google Scholar 

  8. Dufresne M, Bailey JA, Dron M, Langin T: Clk1, a serine/threonine protein kinase-encoding gene, is involved in pathogenicity of C. lindemuthianum on common bean. Mol Plant Microbe Interact. 1998, 11: 99-108. 10.1094/MPMI.1998.11.2.99.

    Article  CAS  Google Scholar 

  9. Hwang C-H, Kolattukudy PE: Isolation and characterisation of genes expressed uniquely during appressorium formation by Colletotrichum gloeosporioides conidia induced by the host surface wax. Mol Gen Genet. 1995, 247: 282-294. 10.1007/BF00293196.

    Article  CAS  Google Scholar 

  10. Wijesundera RLC, Bailey JA, Byrde RJW, Fielding AH: Cell wall degrading enzymes of Colletotrichum lindemuthianum: their role in the development of bean anthracnose. Physiological and Molecular Plant Pathology. 1989, 34: 403-413. 10.1016/0885-5765(89)90067-2.

    Article  CAS  Google Scholar 

  11. Sarah EP, Bleddyn HH, Richard JO, Jonathan RG: Colletotrichum: A Model Genus for Studies on Pathology and Fungal-Plant Interactions. Fungal genetics and Biology. 1999, 27 (2-3): 186-196. 10.1006/fgbi.1999.1143.

    Article  Google Scholar 

  12. Abang MM, Green KR, Asiedu R: Inheritance of resistance in water yam (Dioscorea alata) to anthracnose (Colletotrichum gloesoporioides). Theoret Appl Genet. 2001, 103: 52-55. 10.1007/s001220000531.

    Article  Google Scholar 

  13. Terauchi R, Konuma A: Microsatellite polymorphism in Dioscorea tokoro , a wild yam species. Genome. 1994, 37: 794-801.

    Article  CAS  Google Scholar 

  14. Mignouna HD, Mank RA, Ellis THN, an den Bosch N, Asiedu R, Abang MM, Peleman J: A genetic linkage map of water yam (Discorea alata L.) based on AFLP markers and QTL analysis for anthracose resistance. Theoret Appl Genet. 2002, 105: 726-735. 10.1007/s00122-002-0912-6.

    Article  CAS  Google Scholar 

  15. Abang MM, Onasanya A, Asiedu R: Identification and application of RAPD markers for anthracnose resistance in water yam (Dioscorea alata). Ann Appl Biol. 2002, 141: 61-66. 10.1111/j.1744-7348.2002.tb00195.x.

    Article  Google Scholar 

  16. Morgante M, Olivieri AM: PCR-amplified microsatellites as markers in plant genetics. Plant J. 1993, 3 (1): 175-182. 10.1111/j.1365-313X.1993.tb00020.x.

    Article  CAS  Google Scholar 

  17. Satya SN, Mohamed IA, Asiedu RS, Mignouna RH: A Standardized RNA Isolation Protocol for Yam (Dioscorea alata L) cDNA Library Construction. Virginia Journal of sciences. 2009, 60 (4): 171-180.

    Google Scholar 

  18. Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36 (10): 3420-3435. 10.1093/nar/gkn176.

    Article  CAS  Google Scholar 

  19. Lee SW, Han SW, Bartley LE, Ronald PC: From the Academy: Colloquium review. Unique characteristics Xanthomonas oryzae pv. oryzae AvrXa21 and implications for plant innate immunity. Proceeding of the National Academy of Sciences USA. 2006, 103: 18395-400. 10.1073/pnas.0605508103.

    Article  CAS  Google Scholar 

  20. Kottapalli KR, Rakwal R, Satoh K, Shibato J, Kottapalli P, Iwahashi H, Kikuchi S: Transcriptional profiling of indica rice cultivar IET8585 (Ajaya) infected with bacterial leaf blight pathogen Xanthomonas oryzae pv oryzae. Plant Physiol Biochem. 2007, 45: 834-850. 10.1016/j.plaphy.2007.07.013.

    Article  CAS  Google Scholar 

  21. Larkindale J, Vierling E: Core genome responses involved in acclimation to high temperature. Plant Physiol. 2008, 146: 748-761. 10.1104/pp.107.112060.

    Article  CAS  Google Scholar 

  22. Hirt Heribert: A new Blue print for plant pathogen resistance. News and Views. 2010, [http://www.heribert-hirt.info/pdf/prr80_BlueprintPlantPathogenResistance.pdf]

    Google Scholar 

  23. Kaur N, Reumann S, Hu J: Peroxisome biogenesis and function. The Arabidopsis book. The American Society of Plant Biologists. 2009, 1-41.

    Google Scholar 

  24. Kawano T: Roles of the reactive oxygen species-generating peroxidase reactions in plant defense and growth induction. Plant cell reports. 2003, 21 (9): 829-837.

    Article  CAS  Google Scholar 

  25. Nir Y, Stanley F, Amos D, Noel TK, Dov P: Expression of Pectate Lyase from Colletotrichum gloeosporioides in C. magna Promotes Pathogenicity. The American Pathological Society. 2000, 13 (8): 887-891.

    Google Scholar 

  26. Annis LS, Goodwin HP: Recent advances in the molecular genetics of plant cell wall-degrading enzymes produced by plant pathogenic fungi. Eur J Plant Pathol. 1997, 103: 1-14. 10.1023/A:1008656013255.

    Article  CAS  Google Scholar 

  27. Hamer JE, Holden DW: Linking approaches in the study of fungal pathogenesis: a commentary. Fungal Genet Biol. 1997, 21: 11-16. 10.1006/fgbi.1997.0964.

    Article  Google Scholar 

  28. Hedrick SA, Bell JN, T Lamb CJ: Chitinase cDNA cloning and mRNA induction by fungal elicitor, wounding and infection. Plant Physiol. 1988, 86: 182-186. 10.1104/pp.86.1.182.

    Article  CAS  Google Scholar 

  29. Abang MM: Genetic diversity of Colletotrichum gloeosporioides Penz. causing anthracnose disease of yam (Dioscorea spp.) in Nigeria. Bibliotheca Mycologica. 2003, 197: J. Cramer in der Gebr. Borntraeger Science Publishers Berlin, Stuttgart (ISBN 3-443-59099-3 paperback). 139 pp

    Google Scholar 

  30. Mignouna HD, Abang MM, Green KR, Asiedu R: Inheritance of resistance in water yam (Dioscorea alata) to anthracnose (Colletotrichum gloeosporioides). Theor Appl Genet. 2001, 103: 52-55. 10.1007/s001220000531.

    Article  Google Scholar 

  31. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.

    Article  CAS  Google Scholar 

  32. Green KR, Abang MM, Iloba C: A rapid bioassay for screening yam germplasm for response to anthracnose. Tropical Sci. 2000, 40: 132-8.

    Google Scholar 

  33. Thimm O, Blaesing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller 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-39. 10.1111/j.1365-313X.2004.02016.x.

    Article  CAS  Google Scholar 

  34. Schroeder Steven: University of Missouri-Columbia, Maize Mapping Project 2003-01-01., [http://www.maizemap.org/bioinformatics/SSRFINDER/README.ssrfinder], SchroederSG@missouri.edu

  35. Tang1 Jifeng, Vosman Ben, Voorrips Roeland, van der Linden C Gerard, Leunissen Jack: QualitySNP: a pipeline for detecting single nucleotide polymorphisms and insertions/deletions in EST data from diploid and polyploid species. BMC Bioinformatics. 2006, 7: 438-10.1186/1471-2105-7-438.

    Article  Google Scholar 

Download references

Acknowledgements

This project was supported by a USAID linkage grant through IITA awarded to Mohamed Ali, Mignouna HD Robert and Asiedu Robert. Authors like to acknowledge Dr. Siva Kumpatla and Dr. Teklu Andebrahan for helpful hints during execution of the experiment, Pradeep K Sripathi for continued support and Linda Ballard for BLAST search to fungal gene database. Authors acknowledge the support from MapMan team for annotation and generating the mapping file.

Author information

Authors and Affiliations

  1. USDA-ARS, Stoneville, MS, USA

    Satya S Narina & Brian E Scheffler

  2. Dow AgroSciences, Indianapolis, IN, USA

    Ramesh Buyyarapu

  3. Department of Plant and Soil Science, Texas Tech University, Lubbock, Texas, 79409, USA

    Kameswara Rao Kottapalli

  4. International Institute for Tropical Agriculture (IITA), Oyo Road, PMB 5320, Ibadan, Nigeria

    Alieu M Sartie & Asiedu Robert

  5. USDA-NIFA, Washington, D.C., USA

    Mohamed I Ali

  6. African Agricultural Technology Foundation, Nairobi, 00100, Kenya

    Mignouna JD Hodeba

  7. Virginia State University, Petersburg, VA23806, USA

    Brian L Sayre

Authors
  1. Satya S Narina
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ramesh Buyyarapu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Kameswara Rao Kottapalli
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Alieu M Sartie
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Mohamed I Ali
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Asiedu Robert
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Mignouna JD Hodeba
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Brian L Sayre
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Brian E Scheffler
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Satya S Narina or Brian E Scheffler.

Additional information

Authors' contributions

SN conducted the major part of this study, analyzed the data and prepared the manuscript. RB participated in data analysis, manuscript preparation and revisions. KRK analyzed the data using MapMan software, manuscript preparation and reviewed the manuscript. MA supervised the study and reviewed the manuscript. AR and AS supplied the challenged leaf material and participated in manuscript revisions. MH, BLS and BS reviewed the manuscript. All the authors read and approved the final manuscript.

Electronic supplementary material

12864_2010_9900_MOESM1_ESM.XLS

Additional file 1:Gene Ontology (GO IDs) identifiers for three libraries. Mapping and annotation steps in Blast2GO program resulted in providing GO IDs to individual ESTs in each dataset of the three libraries, TDa 95-0310, TDa 95-0328 and TDa 87-01091. (XLS 7 MB)

12864_2010_9900_MOESM2_ESM.XLS

Additional file 2:Yam sequences annotated using Mercator tool of MapMan database. Yam ESTs for three libraries with MapMan annotations were presented in sheet 1, Mapping file generated from Mercator tool with Bin Codes for each yam EST matching to specific genes on sheet 2 and unique Bin IDs identified in each library related to anthracnose were listed on sheet 3 of the file. (XLS 13 MB)

12864_2010_9900_MOESM3_ESM.XLS

Additional file 3:Differential gene expression among the three libraries. The summary of differential gene expression analysis providing the gene information that might be potentially over-expressed or under-expressed compared to other libraries. The sheet 1 of the file has expression differences between TDa 87-01091 compared with TDa 95-0310 and the sheet 2 has expression differences between TDa 95/0328 compared with TDa 95/0310 while sheet 3 has expression differences between TDa 87-01091 compared with TDa 95-0328. (XLS 306 KB)

12864_2010_9900_MOESM4_ESM.XLS

Additional file 4:Simple Sequence Repeat (SSR) primers generated from the project. The SSR primers generated from library TDa 95-310, TDa 95-0328 and TDa 87-01091 were respectively presented on sheet 1, 2 and 3 of the file. (XLS 992 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Narina, S.S., Buyyarapu, R., Kottapalli, K.R. et al. Generation and analysis of expressed sequence tags (ESTs) for marker development in yam (Dioscorea alata L.). BMC Genomics 12, 100 (2011). https://doi.org/10.1186/1471-2164-12-100

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-12-100

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us