Advancing the STMS genomic resources for defining new locations on the intraspecific genetic linkage map of chickpea (Cicer arietinum L.)

<p>Abstract</p> <p>Background</p> <p>Chickpea (<it>Cicer arietinum </it>L.) is an economically important cool season grain legume crop that is valued for its nutritive seeds having high protein content. However, several biotic and abiotic stresses and the low genetic variability in the chickpea genome have continuously hindered the chickpea molecular breeding programs. STMS (Sequence Tagged Microsatellite Sites) markers which are preferred for the construction of saturated linkage maps in several crop species, have also emerged as the most efficient and reliable source for detecting allelic diversity in chickpea. However, the number of STMS markers reported in chickpea is still limited and moreover exhibit low rates of both inter and intraspecific polymorphism, thereby limiting the positions of the SSR markers especially on the intraspecific linkage maps of chickpea. Hence, this study was undertaken with the aim of developing additional STMS markers and utilizing them for advancing the genetic linkage map of chickpea which would have applications in QTL identification, MAS and for <it>de novo </it>assembly of high throughput whole genome sequence data.</p> <p>Results</p> <p>A microsatellite enriched library of chickpea (enriched for <b>(</b>GT/CA)_nand (GA/CT)_nrepeats) was constructed from which 387 putative microsatellite containing clones were identified. From these, 254 STMS primers were designed of which 181 were developed as functional markers. An intraspecific mapping population of chickpea, [ICCV-2 (single podded) × JG-62 (double podded)] and comprising of 126 RILs, was genotyped for mapping. Of the 522 chickpea STMS markers (including the double-podding trait, screened for parental polymorphism, 226 (43.3%) were polymorphic in the parents and were used to genotype the RILs. At a LOD score of 3.5, eight linkage groups defining the position of 138 markers were obtained that spanned 630.9 cM with an average marker density of 4.57 cM. Further, based on the common loci present between the current map and the previously published chickpea intraspecific map, integration of maps was performed which revealed improvement of marker density and saturation of the region in the vicinity of <it>sfl </it>(double-podding) gene thereby bringing about an advancement of the current map.</p> <p>Conclusion</p> <p>An arsenal of 181 new chickpea STMS markers was reported. The developed intraspecific linkage map defined map positions of 138 markers which included 101 new locations.Map integration with a previously published map was carried out which revealed an advanced map with improved density. This study is a major contribution towards providing advanced genomic resources which will facilitate chickpea geneticists and molecular breeders in developing superior genotypes with improved traits.</p

Transcriptome Analysis in Chickpea (Cicer arietinum L.): Applications in Study of Gene Expression, Non-Coding RNA Prediction, and Molecular Marker Development

Extensive analyses of transcriptome have been carried out in chickpea, which is the third most important legume valued as a source of dietary protein and micronutrients. Over the last two decades, several laboratories have used a wide range of techniques encompassing expressed sequence tag (EST) analysis, serial analysis of gene expression (SAGE), microarray and next-generation sequencing (NGS) technologies for analysing the chickpea transcriptomes. However, chickpea transcriptome analysis witnessed significant progress with the advent of the NGS platforms. Gene expression analyses using NGS platforms were carried out in the vegetative and reproductive tissues such as shoot, root, mature leaf, flower bud, young pod, seed and nodule by various groups which resulted in identification of several tissue-specific transcripts. Some laboratories have utilized transcriptomics to explore the response of chickpea to abiotic and biotic stresses such as drought, salinity, heat, cold, Fusarium oxysporum and Ascochyta rabiei differentially expressed genes and also established crosstalk between biotic and abiotic stress responses. Transcriptome analysis has been utilized extensively to identify non-coding RNAs such as miRNAs and long intergenic non-coding (LINC) RNAs. Transcriptome analysis has facilitated the development of molecular markers such as simple sequence repeats (SSRs), single-nucleotide polymorphisms (SNPs) and potential intron polymorphisms (PIPs) that are being used to expedite the chickpea breeding programmes. The available chickpea transcriptomes will continue to serve as the foundation for devising strategies for chickpea improvement

High density linkage mapping of genomic and transcriptomic SNPs for synteny analysis and anchoring the genome sequence of chickpea

This study presents genome-wide discovery of SNPs through next generation sequencing of the genome of Cicer reticulatum. Mapping of the C. reticulatum sequenced reads onto the draft genome assembly of C. arietinum (desi chickpea) resulted in identification of 842,104 genomic SNPs which were utilized along with an additional 36,446 genic SNPs identified from transcriptome sequences of the aforementioned varieties. Two new chickpea Oligo Pool All (OPAs) each having 3,072 SNPs were designed and utilized for SNP genotyping of 129 Recombinant Inbred Lines (RILs). Using Illumina GoldenGate Technology genotyping data of 5,041 SNPs were generated and combined with the 1,673 marker data from previously published studies, to generate a high resolution linkage map. The map comprised of 6698 markers distributed on eight linkage groups spanning 1083.93 cM with an average inter-marker distance of 0.16 cM. Utility of the present map was demonstrated for improving the anchoring of the earlier reported draft genome sequence of desi chickpea by ~30% and that of kabuli chickpea by 18%. The genetic map reported in this study represents the most dense linkage map of chickpea , with the potential to facilitate efficient anchoring of the draft genome sequences of desi as well as kabuli chickpea varieties

An advanced draft genome assembly of a desi type chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is an important pulse legume crop. We previously reported a draft genome assembly of the desi chickpea cultivar ICC 4958. Here we report an advanced version of the ICC 4958 genome assembly (version 2.0) generated using additional sequence data and an improved genetic map. This resulted in 2.7-fold increase in the length of the pseudomolecules and substantial reduction of sequence gaps. The genome assembly covered more than 94% of the estimated gene space and predicted the presence of 30,257 protein-coding genes including 2230 and 133 genes encoding potential transcription factors (TF) and resistance gene homologs, respectively. Gene expression analysis identified several TF and chickpea-specific genes with tissue-specific expression and displayed functional diversification of the paralogous genes. Pairwise comparison of pseudomolecules in the desi (ICC 4958) and the earlier reported kabuli (CDC Frontier) chickpea assemblies showed an extensive local collinearity with incongruity in the placement of large sequence blocks along the linkage groups, apparently due to use of different genetic maps. Single nucleotide polymorphism (SNP)-based mining of intra-specific polymorphism identified more than four thousand SNPs differentiating a desi group and a kabuli group of chickpea genotypes

Development and use of genic molecular markers (GMMs) for construction of a transcript map of chickpea (Cicer arietinum L.)

A transcript map has been constructed by the development and integration of genic molecular markers (GMMs) including single nucleotide polymorphism (SNP), genic microsatellite or simple sequence repeat (SSR) and intron spanning region (ISR)-based markers, on an inter-specific mapping population of chickpea, the third food legume crop of the world and the first food legume crop of India. For SNP discovery through allele re-sequencing, primer pairs were designed for 688 genes/expressed sequence tags (ESTs) of chickpea and 657 genes/ESTs of closely related species of chickpea. High-quality sequence data obtained for 220 candidate genic regions on 2–20 genotypes representing 9 Cicer species provided 1,893 SNPs with an average frequency of 1/35.83 bp and 0.34 PIC (polymorphism information content) value. On an average 2.9 haplotypes were present in 220 candidate genic regions with an average haplotype diversity of 0.6326. SNP2CAPS analysis of 220 sequence alignments, as mentioned above, provided a total of 192 CAPS candidates. Experimental analysis of these 192 CAPS candidates together with 87 CAPS candidates identified earlier through in silico mining of ESTs provided scorable amplification in 173 (62.01%) cases of which predicted assays were validated in 143 (82.66%) cases (CGMM). Alignments of chickpea unigenes with Medicago truncatula genome were used to develop 121 intron spanning region (CISR) markers of which 87 yielded scorable products. In addition, optimization of 77 EST-derived SSR (ICCeM) markers provided 51 scorable markers. Screening of easily assayable 281 markers including 143 CGMMs, 87 CISRs and 51 ICCeMs on 5 parental genotypes of three mapping populations identified 104 polymorphic markers including 90 markers on the inter-specific mapping population. Sixty-two of these GMMs together with 218 earlier published markers (including 64 GMM loci) and 20 other unpublished markers could be integrated into this genetic map. A genetic map developed here, therefore, has a total of 300 loci including 126 GMM loci and spans 766.56 cM, with an average inter-marker distance of 2.55 cM. In summary, this is the first report on the development of large-scale genic markers including development of easily assayable markers and a transcript map of chickpea. These resources should be useful not only for genome analysis and genetics and breeding applications of chickpea, but also for comparative legume genomics

Integrated physical, genetic and genome map of chickpea (Cicer arietinum L.)

Physical map of chickpea was developed for the reference chickpea genotype (ICC 4958) using bacterial artificial chromosome (BAC) libraries targeting 71,094 clones (~12× coverage). High information content fingerprinting (HICF) of these clones gave high-quality fingerprinting data for 67,483 clones, and 1,174 contigs comprising 46,112 clones and 3,256 singletons were defined. In brief, 574 Mb genome size was assembled in 1,174 contigs with an average of 0.49 Mb per contig and 3,256 singletons represent 407 Mb genome. The physical map was linked with two genetic maps with the help of 245 BAC-end sequence (BES)-derived simple sequence repeat (SSR) markers. This allowed locating some of the BACs in the vicinity of some important quantitative trait loci (QTLs) for drought tolerance and reistance to Fusarium wilt and Ascochyta blight. In addition, fingerprinted contig (FPC) assembly was also integrated with the draft genome sequence of chickpea. As a result, ~965 BACs including 163 minimum tilling path (MTP) clones could be mapped on eight pseudo-molecules of chickpea forming 491 hypothetical contigs representing 54,013,992 bp (~54 Mb) of the draft genome. Comprehensive analysis of markers in abiotic and biotic stress tolerance QTL regions led to identification of 654, 306 and 23 genes in drought tolerance “QTL-hotspot” region, Ascochyta blight resistance QTL region and Fusarium wilt resistance QTL region, respectively. Integrated physical, genetic and genome map should provide a foundation for cloning and isolation of QTLs/genes for molecular dissection of traits as well as markers for molecular breeding for chickpea improvement

Dissecting the Root Nodule Transcriptome of Chickpea (Cicer arietinum L.).

A hallmark trait of chickpea (Cicer arietinum L.), like other legumes, is the capability to convert atmospheric nitrogen (N2) into ammonia (NH3) in symbiotic association with Mesorhizobium ciceri. However, the complexity of molecular networks associated with the dynamics of nodule development in chickpea need to be analyzed in depth. Hence, in order to gain insights into the chickpea nodule development, the transcriptomes of nodules at early, middle and late stages of development were sequenced using the Roche 454 platform. This generated 490.84 Mb sequence data comprising 1,360,251 reads which were assembled into 83,405 unigenes. Transcripts were annotated using Gene Ontology (GO), Cluster of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways analysis. Differential expression analysis revealed that a total of 3760 transcripts were differentially expressed in at least one of three stages, whereas 935, 117 and 2707 transcripts were found to be differentially expressed in the early, middle and late stages of nodule development respectively. MapMan analysis revealed enrichment of metabolic pathways such as transport, protein synthesis, signaling and carbohydrate metabolism during root nodulation. Transcription factors were predicted and analyzed for their differential expression during nodule development. Putative nodule specific transcripts were identified and enriched for GO categories using BiNGO which revealed many categories to be enriched during nodule development, including transcription regulators and transporters. Further, the assembled transcriptome was also used to mine for genic SSR markers. In conclusion, this study will help in enriching the transcriptomic resources implicated in understanding of root nodulation events in chickpea

Isolation of <i>Catharanthus roseus</i> (L.) G. Don Nuclei and Measurement of Rate of <i>Tryptophan decarboxylase</i> Gene Transcription Using Nuclear Run-On Transcription Assay

<div><p>Background</p><p>An accurate assessment of transcription ‘rate’ is often desired to describe the promoter activity. In plants, isolation of transcriptionally active nuclei and their subsequent use in nuclear run-on assays has been challenging and therefore limit an accurate measurement of gene transcription ‘rate’. <i>Catharanthus roseus</i> has emerged as a model medicinal plant as it exhibits an unsurpassed spectrum of chemodiversity, producing over 130 alkaloids through the terpenoid indole alkaloid (TIA) pathway and therefore serves as a ‘molecular hub’ to understand gene expression profiles.</p><p>Results</p><p>The protocols presented here streamline, adapt and optimize the existing methods of nuclear run-on assay for use in <i>C</i>. <i>roseus</i>. Here, we fully describe all the steps to isolate transcriptionally active nuclei from <i>C</i>. <i>roseus</i> leaves and utilize them to perform nuclear run-on transcription assay. Nuclei isolated by this method transcribed at a level consistent with their response to external stimuli, as transcription rate of TDC gene was found to be higher in response to external stimuli i.e. when seedlings were subjected to UV-B light or to methyl jasmonate (MeJA). However, the relative transcript abundance measured parallel through qRT-PCR was found to be inconsistent with the synthesis rate indicating that some post transcriptional events might have a role in transcript stability in response to stimuli.</p><p>Conclusions</p><p>Our study provides an optimized, efficient and inexpensive method of isolation of intact nuclei and nuclear ‘run-on’ transcription assay to carry out <i>in-situ</i> measurement of gene transcription rate in <i>Catharanthus roseus</i>. This would be valuable in investigating the transcriptional and post transcriptional response of other TIA pathway genes in <i>C</i>. <i>roseus</i>. Isolated nuclei may also provide a resource that could be used for performing the chip assay as well as serve as the source of nuclear proteins for <i>in-vitro</i> EMSA studies. Moreover, nascent nuclear run-on transcript could be further subjected to RNA-Seq for global nuclear run-on assay (GNRO-Seq) for genome wide <i>in-situ</i> measurement of transcription rate of plant genes.</p></div