A knowledge-driven protocol for prediction of proteins of interest with an emphasis on biosynthetic pathways

This protocol describes a stepwise process to identify proteins of interest from a query proteome derived from NGS data. We implemented this protocol on Moringa oleifera transcriptome to identify proteins involved in secondary metabolite and vitamin biosynthesis and ion transport. This knowledge-driven protocol identifies proteins using an integrated approach involving sensitive sequence search and evolutionary relationships. We make use of functionally important residues (FIR) specific for the query protein family identified through its homologous sequences and literature. We screen protein hits based on the clustering with true homologues through phylogenetic tree reconstruction complemented with the FIR mapping. The protocol was validated for the protein hits through qRT-PCR and transcriptome quantification. Our protocol demonstrated a higher specificity as compared to other methods, particularly in distinguishing cross-family hits. This protocol was effective in transcriptome data analysis of M. oleifera as described in Pasha et al.•

Knowledge-driven protocol to identify secondary metabolite synthesizing protein in a highly specific manner.

•

Use of functionally important residues for screening of true hits.

•

Beneficial for metabolite pathway reconstruction in any (species, metagenomics) NGS data.

Genome sequencing of herb Tulsi (Ocimum tenuiflorum) unravels key genes behind its strong medicinal properties

BACKGROUND

Krishna Tulsi, a member of Lamiaceae family, is a herb well known for its spiritual, religious and medicinal importance in India. The common name of this plant is 'Tulsi' (or 'Tulasi' or 'Thulasi') and is considered sacred by Hindus. We present the draft genome of Ocimum tenuiflurum L (subtype Krishna Tulsi) in this report. The paired-end and mate-pair sequence libraries were generated for the whole genome sequenced with the Illumina Hiseq 1000, resulting in an assembled genome of 374 Mb, with a genome coverage of 61 % (612 Mb estimated genome size). We have also studied transcriptomes (RNA-Seq) of two subtypes of O. tenuiflorum, Krishna and Rama Tulsi and report the relative expression of genes in both the varieties.

RESULTS

The pathways leading to the production of medicinally-important specialized metabolites have been studied in detail, in relation to similar pathways in Arabidopsis thaliana and other plants. Expression levels of anthocyanin biosynthesis-related genes in leaf samples of Krishna Tulsi were observed to be relatively high, explaining the purple colouration of Krishna Tulsi leaves. The expression of six important genes identified from genome data were validated by performing q-RT-PCR in different tissues of five different species, which shows the high extent of urosolic acid-producing genes in young leaves of the Rama subtype. In addition, the presence of eugenol and ursolic acid, implied as potential drugs in the cure of many diseases including cancer was confirmed using mass spectrometry.

CONCLUSIONS

The availability of the whole genome of O.tenuiflorum and our sequence analysis suggests that small amino acid changes at the functional sites of genes involved in metabolite synthesis pathways confer special medicinal properties to this herb.

Comparative analyses of stress-responsive genes in Arabidopsis thaliana: insight from genomic data mining, functional enrichment, pathway analysis and phenomics

Biotic and abiotic stresses adversely affect agriculture by reducing crop growth and productivity worldwide. To investigate the abiotic stress-responsive genes in Arabidopsis thaliana, we compiled a dataset of stress signals and differentially upregulated genes (>= 2.5 fold change) from Stress-responsive transcription Factors DataBase (STIFDB) with additional set of stress signals and genes curated from PubMed and Gene Expression Omnibus. A dataset of 3091 genes differentially upregulated due to 14 different stress signals (abscisic acid, aluminum, cold, cold-drought-salt, dehydration, drought, heat, iron, light, NaCl, osmotic stress, oxidative stress, UV-B and wounding) were curated and used for the analysis. Details about stress-responsive enriched genes and their association with stress signals can be obtained from STIFDB2 database . The gene-stress-signal data were analyzed using an enrichment-based meta-analysis framework consisting of two different ontologies (Gene Ontology and Plant Ontology), biological pathway and functional domain annotations. We found several shared and distinct biological processes, cellular components and molecular functions associated with stress-responsive genes. Pathway analysis revealed that stress-responsive genes perturbed the pathways under the "Metabolic pathways" category. We also found several shared and stress-signal specific protein domains, suggesting functional mechanisms regulating stress-response. Phenomic characteristics of abiotic stress-responsive genes were ascertained for several stresses and found to be shared by multiple stresses in both anatomy and temporal categories of Plant Ontology. We found several constitutive stress-responsive genes that are differentially upregulated due to perturbation of different stress signals, for example a gene (AT1G68440) involved in phenylpropanoid metabolism and polyamine catabolism as responsive to seven different stress signals. We also performed structure-function prediction of five genes associated responsive to multiple abiotic stress signals. We envisage that results from our analysis that provide insight into functional repertoire, metabolic pathways and phenomic characteristics common and specifically associated with stress signals would help to understand abiotic stress regulome in Arabidopsis thaliana and may also help to develop an improved plant variety using molecular breeding and genetic engineering techniques that are rapidly stress-responsive and tolerant.

STIFDB2: an updated version of plant stress-responsive transcription factor database with additional stress signals, stress-responsive transcription factor binding sites and stress-responsive genes in Arabidopsis and rice

Understanding the principles of abiotic and biotic stress responses, tolerance and adaptation remains important in plant physiology research to develop better varieties of crop plants. Better understanding of plant stress response mechanisms and application of knowledge derived from integrated experimental and bioinformatics approaches are gaining importance. Earlier, we showed that compiling a database of stress-responsive transcription factors and their corresponding target binding sites in the form of Hidden Markov models at promoter, untranslated and upstream regions of stress-up-regulated genes from expression analysis can help in elucidating various aspects of the stress response in Arabidopsis. In addition to the extensive content in the first version, STIFDB2 is now updated with 15 stress signals, 31 transcription factors and 5,984 stress-responsive genes from three species (Arabidopsis thaliana, Oryza sativa subsp. japonica and Oryza sativa subsp. indica). We have employed an integrated biocuration and genomic data mining approach to characterize the data set of transcription factors and consensus binding sites from literature mining and stress-responsive genes from the Gene Expression Omnibus. STIFDB2 currently has 38,798 associations of stress signals, stress-responsive genes and transcription factor binding sites predicted using the Stress-responsive Transcription Factor (STIF) algorithm, along with various functional annotation data. As a unique plant stress regulatory genomics data platform, STIFDB2 can be utilized for targeted as well as high-throughput experimental and computational studies to unravel principles of the stress regulome in dicots and gramineae. STIFDB2 is available from the URL: http://caps.ncbs.res.in/stifdb2.