Comprehensive Transcriptome Assembly of Chickpea (Cicer arietinum L.) Using Sanger and Next Generation Sequencing Platforms: Development and Applications

• Published in: PloS one 2014-01, Vol.9 (1), p.e86039-e86039
• Main Authors: Kudapa, Himabindu, Azam, Sarwar, Sharpe, Andrew G, Taran, Bunyamin, Li, Rong, Deonovic, Benjamin, Cameron, Connor, Farmer, Andrew D, Cannon, Steven B, Varshney, Rajeev K, Lukens, Lewis
• Format: Article
• Language: English
• Subjects: Agriculture; Alfalfa; Assembly; Bioinformatics; Biology; Breeding; Cajanus cajan; chickpeas; Cicer - genetics; Cicer - metabolism; Cicer arietinum; Comparative analysis; Conservation; Contig Mapping; Crops; Datasets; Expressed sequence tags; Gene expression; Gene Expression Profiling; Gene mapping; Gene polymorphism; Genes; Genes, Plant; Genetic aspects; Genetic Linkage; Genetic Markers; Genetics; Genomes; Genomics; genotype; Genotypes; Glycine max; High-Throughput Nucleotide Sequencing; Introns; Legumes; linkage groups; Markers; Medicago; Medicago truncatula - genetics; Microsatellite Repeats; Molecular Sequence Annotation; Phaseolus vulgaris; Plant Proteins - genetics; Plant Proteins - metabolism; Polymorphism; Salinity; Sequence Analysis, DNA; Simple sequence repeats; Soybeans; Synteny; Transcription; Transcriptome
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0086039

A comprehensive transcriptome assembly of chickpea has been developed using 134.95 million Illumina single-end reads, 7.12 million single-end FLX/454 reads and 139,214 Sanger expressed sequence tags (ESTs) from >17 genotypes. This hybrid transcriptome assembly, referred to as C icer a rietinum T ranscriptome A ssembly version 2 (CaTA v2, available at http://data.comparative-legumes.org/transcriptomes/cicar/lista_cicar-201201), comprising 46,369 transcript assembly contigs (TACs) has an N50 length of 1,726 bp and a maximum contig size of 15,644 bp. Putative functions were determined for 32,869 (70.8%) of the TACs and gene ontology assignments were determined for 21,471 (46.3%). The new transcriptome assembly was compared with the previously available chickpea transcriptome assemblies as well as to the chickpea genome. Comparative analysis of CaTA v2 against transcriptomes of three legumes - Medicago , soybean and common bean, resulted in 27,771 TACs common to all three legumes indicating strong conservation of genes across legumes. CaTA v2 was also used for identification of simple sequence repeats (SSRs) and intron spanning regions (ISRs) for developing molecular markers. ISRs were identified by aligning TACs to the Medicago genome, and their putative mapping positions at chromosomal level were identified using transcript map of chickpea. Primer pairs were designed for 4,990 ISRs, each representing a single contig for which predicted positions are inferred and distributed across eight linkage groups. A subset of randomly selected ISRs representing all eight chickpea linkage groups were validated on five chickpea genotypes and showed 20% polymorphism with average polymorphic information content (PIC) of 0.27. In summary, the hybrid transcriptome assembly developed and novel markers identified can be used for a variety of applications such as gene discovery, marker-trait association, diversity analysis etc., to advance genetics research and breeding applications in chickpea and other related legumes.