Genome-wide analysis of glutathione S-transferase gene family in chickpea suggests its role during seed development and abiotic stress.

Glutathione S-transferases (GSTs) are multifunctional proteins that help in oxidative stress metabolism and detoxification of xenobiotic compounds. Studies pertaining to GST gene family have been undertaken in various plant species, however no information is available with respect to GST genes in chickpea. In the current study, we identified a total of 51 GST encoding genes in chickpea (CaGST) genome. Phylogenetic analysis revealed that GST gene family can be divided into eleven distinct classes. Tau and phi were the major classes in chickpea and one third of the CaGST genes represented segmental duplication and purifying selection was common among these genes. Expression of many CaGST genes, in particular, members of tau class were found to be upregulated under abiotic stress conditions. In addition, CaGST genes displayed differential expression patterns across diverse organs/tissues, suggesting their roles in developmental processes. Many CaGST genes showed opposite expression pattern in small- and large-seeded chickpea cultivars during seed development. Higher expression of CaGST genes in small-seeded cultivar at maturation stages of seed development suggested their important role in seed development and seed size/weight determination in chickpea. Overall, these results provide a comprehensive information on GST gene family members in chickpea and is expected to provide a rational platform to explore versatile role of these genes in semi-arid legume crops.

Development of EST-based new SSR markers in seabuckthorn.

EST-based SSR markers were developed by screening a collection of 1584 clustered ESTs of seabuckthorn (Hippophae rhamnoides). PCR primers were designed for the amplification of 30 microsatellite loci. Two to five allelic bands were displayed by nine primer pairs in H. rhamnoides genotypes and by eleven primer pairs in H. salicifolia genotypes. None of the thirty primer pairs detected polymorphism in H. tibetana genotypes. Considering the high polymorphism detected in the tested genotypes and their direct origin from the genic regions, these EST-SSR markers hold immense promise in seabuckthorn genome analysis, molecular breeding and population genetics.

Updates on Genomic Resources in Chickpea for Crop Improvement.

In recent years, rapid advancement has been done in generation of genomic resources for the important legume crop chickpea. Here, we provide an update on important advancements made on availability of genomic resources for this crop. The availability of reference genome and transcriptome sequences, and resequencing of several accessions have enabled the discovery of gene space and molecular markers in chickpea. These resources have helped in elucidating evolutionary relationship and identification of quantitative trait loci for important agronomic traits. Gene expression in different tissues/organs during development and under abiotic/biotic stresses has been interrogated. In addition, single-base resolution DNA methylation patterns in different organs have been analyzed to understand gene regulation. Overall, we provide a consolidated overview of available genomic resources of chickpea that may help in fulfilling the promises for improvement of this important crop.

De novo transcriptome assembly and comprehensive expression profiling in Crocus sativus to gain insights into apocarotenoid biosynthesis.

Saffron (Crocus sativus L.) is commonly known as world's most expensive spice with rich source of apocarotenoids and possesses magnificent medicinal properties. To understand the molecular basis of apocarotenoid biosynthesis/accumulation, we performed transcriptome sequencing from five different tissues/organs of C. sativus using Illumina platform. After comprehensive optimization of de novo transcriptome assembly, a total of 105, 269 unique transcripts (average length of 1047 bp and N50 length of 1404 bp) were obtained from 206 million high-quality paired-end reads. Functional annotation led to the identification of many genes involved in various biological processes and molecular functions. In total, 54% of C. sativus transcripts could be functionally annotated using public databases. Transcriptome analysis of C. sativus revealed the presence of 16721 SSRs and 3819 transcription factor encoding transcripts. Differential expression analysis revealed preferential/specific expression of many transcripts involved in apocarotenoid biosynthesis in stigma. We have revealed the differential expression of transcripts encoding for transcription factors (MYB, MYB related, WRKY, C2C2-YABBY and bHLH) involved in secondary metabolism. Overall, these results will pave the way for understanding the molecular basis of apocarotenoid biosynthesis and other aspects of stigma development in C. sativus.

Genome-wide analysis of homeobox gene family in legumes: identification, gene duplication and expression profiling.

Homeobox genes encode transcription factors that are known to play a major role in different aspects of plant growth and development. In the present study, we identified homeobox genes belonging to 14 different classes in five legume species, including chickpea, soybean, Medicago, Lotus and pigeonpea. The characteristic differences within homeodomain sequences among various classes of homeobox gene family were quite evident. Genome-wide expression analysis using publicly available datasets (RNA-seq and microarray) indicated that homeobox genes are differentially expressed in various tissues/developmental stages and under stress conditions in different legumes. We validated the differential expression of selected chickpea homeobox genes via quantitative reverse transcription polymerase chain reaction. Genome duplication analysis in soybean indicated that segmental duplication has significantly contributed in the expansion of homeobox gene family. The Ka/Ks ratio of duplicated homeobox genes in soybean showed that several members of this family have undergone purifying selection. Moreover, expression profiling indicated that duplicated genes might have been retained due to sub-functionalization. The genome-wide identification and comprehensive gene expression profiling of homeobox gene family members in legumes will provide opportunities for functional analysis to unravel their exact role in plant growth and development.

Transcriptome analysis of Catharanthus roseus for gene discovery and expression profiling.

The medicinal plant, Catharanthus roseus, accumulates wide range of terpenoid indole alkaloids, which are well documented therapeutic agents. In this study, deep transcriptome sequencing of C. roseus was carried out to identify the pathways and enzymes (genes) involved in biosynthesis of these compounds. About 343 million reads were generated from different tissues (leaf, flower and root) of C. roseus using Illumina platform. Optimization of de novo assembly involving a two-step process resulted in a total of 59,220 unique transcripts with an average length of 1284 bp. Comprehensive functional annotation and gene ontology (GO) analysis revealed the representation of many genes involved in different biological processes and molecular functions. In total, 65% of C. roseus transcripts showed homology with sequences available in various public repositories, while remaining 35% unigenes may be considered as C. roseus specific. In silico analysis revealed presence of 11,620 genic simple sequence repeats (excluding mono-nucleotide repeats) and 1820 transcription factor encoding genes in C. roseus transcriptome. Expression analysis showed roots and leaves to be actively participating in bisindole alkaloid production with clear indication that enzymes involved in pathway of vindoline and vinblastine biosynthesis are restricted to aerial tissues. Such large-scale transcriptome study provides a rich source for understanding plant-specialized metabolism, and is expected to promote research towards production of plant-derived pharmaceuticals.

Optimization of de novo short read assembly of seabuckthorn (Hippophae rhamnoides L.) transcriptome.

Seabuckthorn (Hippophaerhamnoides L.) is known for its medicinal, nutritional and environmental importance since ancient times. However, very limited efforts have been made to characterize the genome and transcriptome of this wonder plant. Here, we report the use of next generation massive parallel sequencing technology (Illumina platform) and de novo assembly to gain a comprehensive view of the seabuckthorn transcriptome. We assembled 86,253,874 high quality short reads using six assembly tools. At our hand, assembly of non-redundant short reads following a two-step procedure was found to be the best considering various assembly quality parameters. Initially, ABySS tool was used following an additive k-mer approach. The assembled transcripts were subsequently subjected to TGICL suite. Finally, de novo short read assembly yielded 88,297 transcripts (> 100 bp), representing about 53 Mb of seabuckthorn transcriptome. The average length of transcripts was 610 bp, N50 length 1198 BP and 91% of the short reads uniquely mapped back to seabuckthorn transcriptome. A total of 41,340 (46.8%) transcripts showed significant similarity with sequences present in nr protein databases of NCBI (E-value < 1E-06). We also screened the assembled transcripts for the presence of transcription factors and simple sequence repeats. Our strategy involving the use of short read assembler (ABySS) followed by TGICL will be useful for the researchers working with a non-model organism's transcriptome in terms of saving time and reducing complexity in data management. The seabuckthorn transcriptome data generated here provide a valuable resource for gene discovery and development of functional molecular markers.

Expressed sequence tag based identification and expression analysis of some cold inducible elements in seabuckthorn (Hippophae rhamnoides L.).

A cDNA library was constructed from the mature leaves of seabuckthorn (Hippophae rhamnoides). Expressed Sequence Tags (ESTs) were generated by single pass sequencing of 4500 cDNA clones. We submitted 3412 ESTs to dbEST of NCBI. Clustering of these ESTs yielded 1665 unigenes comprising of 345 contigs and 1320 singletons. Out of 1665 unigenes, 1278 unigenes were annotated by similarity search while the remaining 387 unannotated unigenes were considered as organism specific. Gene Ontology (GO) analysis of the unigene dataset showed 691 unigenes related to biological processes, 727 to molecular functions and 588 to cellular component category. On the basis of similarity search and GO annotation, 43 unigenes were found responsive to biotic and abiotic stresses. To validate this observation, 13 genes that are known to be associated with cold stress tolerance from previous studies in Arabidopsis and 3 novel transcripts were examined by Real time RT-PCR to understand the change in expression pattern under cold/freeze stress. In silico study of occurrence of microsatellites in these ESTs revealed the presence of 62 Simple Sequence Repeats (SSRs), some of which are being explored to assess genetic diversity among seabuckthorn collections. This is the first report of generation of transcriptome data providing information about genes involved in managing plant abiotic stress in seabuckthorn, a plant known for its enormous medicinal and ecological value.

Isolation of good quality RNA from a medicinal plant seabuckthorn, rich in secondary metabolites.

Medicinal plants are being widely investigated owing to their ability to produce molecules of therapeutic significance. Isolation of good quality RNA is a tedious but primary step towards undertaking molecular biology experiments. However, medicinal plants are rich in secondary metabolites and not amenable to standard RNA isolation protocols involving Guanidine isothiocyanate (GITC). So an RNA isolation protocol from difficult samples (richer in secondary metabolites) is of highest desiderata. Here we propose a new protocol suitable for isolating RNA from plant tissues rich in secondary metabolites. To standard CTAB (Cetyl Trimethyl Ammonium Bromide) buffer, addition of 2% PVPP (polyvinyl polypyrrolidone) and 350 mM beta-mercaptoethanol was found useful. Use of glacial acetic acid (1M) along with ethanol for precipitation after phenolization and chloroform extraction enhanced the RNA yield. This is the first report of using glacial acetic acid in a CTAB based protocol for the precipitation of RNA. This protocol has been validated in medicinal plant Hippophae rhamnoides vern. seabuckthorn, where standard RNA isolation methods involving GITC and TRIZol extraction buffers failed. The RNA isolated by this method was of good quality as gauged by spectrophotometric readings and denaturing agarose gel electrophoresis. To the best of our knowledge, this RNA isolation protocol has never been published before. The RNA thus obtained could be suitably used for the downstream molecular procedures like Reverse Transcription Polymerase Chain Reaction (RT-PCR), Real Time-PCR, cDNA library construction, etc.