TinoTranscriptDB: A Database of Transcripts and Microsatellite Markers of Tinospora cordifolia, an Important Medicinal Plant.

Tinospora cordifolia, commonly known as "Giloe" in India, is a shrub belonging to the family Menispermaceae. It is an important medicinal plant known for its antipyretic, anti-inflammatory, antispasmodic, and antidiabetic properties and is used in the treatment of jaundice, gout, and rheumatism. Despite its economic importance, the limited information related to its genomic resources prohibits its judicious exploitation through molecular breeding or biotechnological approaches. In this study, we generated a meta-transcriptome assembly of 43,090 non-redundant transcripts by merging the RNASeq data obtained from Roche 454 GS-FLX, and Illumina platforms, and report the first transcriptome-based database for simple sequence repeats and transcription factors ("TinoTranscriptDB" (Tinospora cordifolia Transcriptome Database)). We annotated 26,716 (62%) of the total transcripts successfully from National Center for Biotechnology Information non-redundant protein (NCBI-NR), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-Prot, and Pfam databases. This database contains information of 2620 perfect simple sequence repeats (P-SSRs) with a relative abundance of 340.12 (loci/Mb), and relative density of 6309.29 (bp/Mb). Excluding mono-nucleotides, the most abundant SSR motifs were tri-nucleotides (54.31%), followed by di-nucleotides (37.51%), tetra-nucleotides (4.54%), penta-nucleotides (3.16%) and hexa-nucleotides (0.45%). Additionally, we also identified 4,311 transcription factors (TFs) and categorized them into 55 sub-families. This database is expected to fill the gap in genomic resource availability in T. cordifolia and thus accelerate molecular breeding and related functional and other applied studies aimed towards genetic improvements of T. cordifolia and related species.

Refinement of Draft Genome Assemblies of Pigeonpea (Cajanus cajan).

Genome assembly of short reads from large plant genomes remains a challenge in computational biology despite major developments in next generation sequencing. Of late several draft assemblies have been reported in sequenced plant genomes. The reported draft genome assemblies of Cajanus cajan have different levels of genome completeness, a large number of repeats, gaps, and segmental duplications. Draft assemblies with portions of genome missing are shorter than the referenced original genome. These assemblies come with low map accuracy affecting further functional annotation and the prediction of gene components as desired by crop researchers. Genome coverage, i.e., the number of sequenced raw reads mapped onto a certain location of the genome is an important quality indicator of completeness and assembly quality in draft assemblies. The present work aimed to improve the coverage in reported de novo sequenced draft genomes (GCA_000340665.1 and GCA_000230855.2) of pigeonpea, a legume widely cultivated in India. The two recently sequenced assemblies, A1 and A2 comprised 72% and 75% of the estimated coverage of the genome, respectively. We employed an assembly reconciliation approach to compare the draft assemblies and merge them, filling the gaps by employing an algorithm size sorting mate-pair library to generate a high quality and near complete assembly with enhanced contiguity. The majority of gaps present within scaffolds were filled with right-sized mate-pair reads. The improved assembly reduced the number of gaps than those reported in draft assemblies resulting in an improved genome coverage of 82.4%. Map accuracy of the improved assembly was evaluated using various quality metrics and for the presence of specific trait-related functional genes. Employed pair-end and mate-pair local libraries helped us to reduce gaps, repeats, and other sequence errors resulting in lengthier scaffolds compared to the two draft assemblies. We reported the prediction of putative host resistance genes against Fusarium wilt disease by their performance and evaluated them both in wet laboratory and field phenotypic conditions.

Comparative genomic analysis of monosporidial and monoteliosporic cultures for unraveling the complexity of molecular pathogenesis of Tilletia indica pathogen of wheat.

Tilletia indica (Ti) - a quarantined fungal pathogen of wheat and its pathogenesis is chiefly governed by pathogen effectors secreted inside the host plant. The de novo genome sequencing of several field isolates and stages available could be used for understanding the molecular pathogenesis. The presence of gaps and low coverage of assembled genomes poses a problem in accurate functional annotation of such functions. In the present study attempts were made to improve the Ti draft genome through reconciliation of globally available datasets of three highly virulent monoteliospore cultures of Ti field isolates. It has sequence depth of 107x and N50 scaffold size of 80,772 (more than 26 times as large as achieved in the draft assembly) with highest sequence contiguity, more accurate and nearly complete. Functional annotation revealed that Ti genome contains 9209 genes evolved with many expanded gene families and arranged mostly in a cluster. About 79% of Ti genes were orthologous to other basidiomycetes fungi, Around 7.93% proteins were having secretary signals and 6.66% were identified as highly virulent pathogenicity genes. Using improved Ti genome as a reference, the genomic variation was assessed with respect to repeats, SNPs/InDel, gene families and correct set of virulence associated genes during its life cycle. The comparative intra-species, inter-stage and inter-species genomic variation will have broader implications to understand the gene regulatory networks involved in growth, mating and virulence behaviour of Tilletia f. spp. and also for better appreciation of fungal biology and disease management.

Improved Draft Genome Sequence of a Monoteliosporic Culture of the Karnal Bunt (Tilletia indica) Pathogen of Wheat.

Karnal bunt of wheat is an internationally quarantined fungal pathogen disease caused by Tilletia indica and affects the international commercial seed trade of wheat. We announce here the first improved draft genome assembly of a monoteliosporic culture of the Tilletia indica fungus, consisting of 787 scaffolds with an approximate total genome size of 31.83 Mbp, which is more accurate and near to complete than the previous version.

RNA-seq of Rice Yellow Stem Borer Scirpophaga incertulas Reveals Molecular Insights During Four Larval Developmental Stages.

The yellow stem borer (YSB), Scirpophaga incertulas, is a prominent pest in rice cultivation causing serious yield losses. The larval stage is an important stage in YSB, responsible for maximum infestation. However, limited knowledge exists on the biology and mechanisms underlying the growth and differentiation of YSB. To understand and identify the genes involved in YSB development and infestation, so as to design pest control strategies, we performed de novo transcriptome analysis at the first, third, fifth, and seventh larval developmental stages employing Illumina Hi-seq. High-quality reads (HQR) of ∼229 Mb were assembled into 24,775 transcripts with an average size of 1485 bp. Genes associated with various metabolic processes, i.e., detoxification mechanism [CYP450, GSTs, and carboxylesterases (CarEs)], RNA interference (RNAi) machinery (Dcr-1, Dcr-2, Ago-1, Ago-2, Sid-1, Sid-2, Sid-3, and Sid-1-related gene), chemoreception (CSPs, GRs, OBPs, and ORs), and regulators [transcription factors (TFs) and hormones] were differentially regulated during the developmental stages. Identification of stage-specific transcripts made it possible to determine the essential processes of larval development. Comparative transcriptome analysis revealed that YSB has not evolved much with respect to the detoxification mechanism, but showed the presence of distinct RNAi machinery. The presence of strong specific visual recognition coupled with chemosensory mechanisms supports the monophagous nature of YSB. Designed expressed sequenced tags-simple-sequence repeats (EST-SSRs) will facilitate accurate estimation of the genetic diversity of YSB. This is the first report on characterization of the YSB transcriptome and the identification of genes involved in key processes, which will help researchers and industry to devise novel pest control strategies. This study also opens up a new avenue to develop next-generation resistant rice using RNAi or genome editing approaches.

Draft genome sequence of Karnal bunt pathogen (Tilletia indica) of wheat provides insights into the pathogenic mechanisms of quarantined fungus.

Karnal bunt disease in wheat is caused by hemibiotrophic fungus, Tilletia indica that has been placed as quarantine pest in more than 70 countries. Despite its economic importance, little knowledge about the molecular components of fungal pathogenesis is known. In this study, first time the genome sequence of T. indica has been deciphered for unraveling the effectors' functions of molecular pathogenesis of Karnal bunt disease. The T. indica genome was sequenced employing hybrid approach of PacBio Single Molecule Real Time (SMRT) and Illumina HiSEQ 2000 sequencing platforms. The genome was assembled into 10,957 contigs (N50 contig length 3 kb) with total size of 26.7 Mb and GC content of 53.99%. The number of predicted putative genes were 11,535, which were annotated with Gene Ontology databases. Functional annotation of Karnal bunt pathogen genome and classification of identified effectors into protein families revealed interesting functions related to pathogenesis. Search for effectors' genes using pathogen host interaction database identified 135 genes. The T. indica genome sequence and putative genes involved in molecular pathogenesis would further help in devising novel and effective disease management strategies including development of resistant wheat genotypes, novel biomarkers for pathogen detection and new targets for fungicide development.

De novo transcriptome sequencing facilitates genomic resource generation in Tinospora cordifolia.

Tinospora cordifolia is known for its medicinal properties owing to the presence of useful constituents such as terpenes, glycosides, steroids, alkaloids, and flavonoids belonging to secondary metabolism origin. However, there is little information available pertaining to critical genomic elements (ESTs, molecular markers) necessary for judicious exploitation of its germplasm. We employed 454 GS-FLX pyrosequencing of entire transcripts and altogether ∼25 K assembled transcripts or Expressed sequence tags (ESTs) were identified. As the interest in T. cordifolia is primarily due to its secondary metabolite constituents, the ESTs pertaining to terpenoids biosynthetic pathway were identified in the present study. Additionally, several ESTs were assigned to different transcription factor families. To validate our transcripts dataset, the novel EST-SSR markers were generated to assess the genetic diversity among germplasm of T. cordifolia. These EST-SSR markers were found to be polymorphic and the dendrogram based on dice similarity index revealed three distinct clustering of accessions. The present study demonstrates effectiveness in using both NEWBLER and MIRA sequence read assembler software for enriching transcript-dataset and thus enables better exploitation of EST resources for mining candidate genes and designing molecular markers.

Expression of CspE by a psychrotrophic bacterium Enterobacter ludwigii PAS1, isolated from Indian Himalayan soil and in silico protein modelling, prediction of conserved residues and active sites.

Proteome analysis of Enterobacter ludwigii PAS1 provide a powerful set of tool to study the cold shock proteins along with that combination of bioinformatics is useful for interpretation of comparative results from many species. There is a considerable interest in the use of psychrotrophic bacteria for nitrogen fixation, especially at hilly regions, thus better understanding of cold adaptation mechanisms too. The psychrotrophic E. ludwigii PAS1 grown at 30 and 4 °C, isolated from Himalaya soil was undertaken for proteomic responses during optimal and cold shock conditions. Comparative proteomic analyses using two-dimensional gel electrophoresis (2-DE) and MALDI-TOF/TOF MS revealed the presence of Cold shock protein E (CspE). Three-dimensional structure of CspE of E. ludwigii PAS1 divulge the presence of five antiparallel ß-sheets forming a ß-barrel structure with surface exposed aromatic and basic residues that were responsible for nucleic acid binding and also reveals the presence of highly conserved nucleic acid-binding motifs RNP1 and RNP2 in Csp family.

Metabolomics in agriculture.

Metabolome refers to the complete set of metabolites synthesized through a series of multiple enzymatic steps from various biochemical pathways processing the information encrypted in the plant genome. Knowledge about synthesis and regulation of various plant metabolic substances has improved substantially with availability of Omics data originating from sequencing of plant genomes. Metabolic profiling of crops is increasingly becoming popular in assessing plant phenotypes and genetic diversity. Metabolic compositional changes vividly reflect the changes occurring during plant growth, development, and in response to stress. Hence, study of plant metabolic pathways, the interconnections between them in context of systems biology is increasingly becoming popular in identification of candidate genes. The present article reviews recent developments in analysis of plant metabolomics, available bioinformatics techniques and databases employed for comparative pathway analysis, metabolic QTLs, and their application in plants.

LOX genes in blast fungus (Magnaporthe grisea) resistance in rice.

Plant Lipoxygenases (LOX) are known to play major role in plant immunity by providing front-line defense against pathogen-induced injury. To verify this, we isolated a full-length OsLOX3 gene and also 12 OsLOX cDNA clones from Oryza sativa indica (cultivar Pusa Basmati 1). We have examined the role played by LOXs in plant development and during attack by blast pathogen Magnaporthe grisea. Gene expression, promoter region analysis, and biochemical and protein structure analysis of isolated OsLOX3 revealed significant homology with LOX super family. Protein sequence comparison of OsLOXs revealed high levels of homology when compared with japonica rice (up to100%) and Arabidopsis (up to 64%). Isolated LOX3 gene and 12 OsLOX cDNAs contained the catalytic LOX domains much required for oxygen binding and synthesis of oxylipins. Amino acid composition, protein secondary structure, and promoter region analysis (with abundance of motifs CGTCA and TGACG) support the role of OsLOX3 gene in providing resistance to diseases in rice plants. OsLOX3 gene expression analysis of root, shoot, flag leaf, and developing and mature seed revealed organ specific patterns during rice plant development and gave evidence to association between tissue location and physiological roles played by individual OsLOXs. Increased defense activity of oxylipins was observed as demonstrated by PCR amplification of OsLOX3 gene and upon inoculation with virulent strains of M. grisea and ectopic application of methyl jasmonate in the injured leaf tissue in adult rice plants.

Validation of arsenic resistance in Bacillus cereus strain AG27 by comparative protein modeling of arsC gene product.

The ars gene system provides arsenic resistance to a variety of microorganisms and can be chromosomal or plasmid-borne. The arsC gene, which codes for an arsenate reductase is essential for arsenate resistance and transforms arsenate into arsenite, which is extruded from the cell. Therefore, arsC gene from Bacillus cereus strain AG27 isolated from soil was amplified, cloned and sequenced. The strain exhibited a minimum inhibitory concentration of 40 and 35 mM to sodium arsenate and sodium arsenite, respectively. Homology of the sequence, when compared with available database using BLASTn search showed that 300 bp amplicons obtained possess partial arsC gene sequence which codes for arsenate reductase, an enzyme involved in the reduction of arsenate to arsenite which is then effluxed out of the cell, thereby indicating the presence of efflux mechanism of resistance in strain. The efflux mechanism was further confirmed by atomic absorption spectroscopy and scanning electron microscopy studies. Moreover, three dimensional structure of modeled arsC from Bacillus cereus strain shares significant structural similarity with arsenate reductase protein of B.subtilis, consisting of, highly similar overall fold with single α/ß domain containing a central four stranded, parallel, open-twisted ß-sheet flanked by α-helices on both sides. The structure harbors the arsenic binding motif AB loop or P-loop that is highly conserved in arsenate reductase family.

In silico analysis of sequential, structural and functional diversity of wheat cystatins and its implication in plant defense.

Phytocystatins constitute a multigene family that regulates the activity of endogenous and/or exogenous cysteine proteinases. Cereal crops like wheat are continuously threatened by a multitude of pathogens, therefore cystatins offer to play a pivotal role in deciding the plant response. In order to study the need of having diverse specificities and activities of various cystatins, we conducted comparative analysis of six wheat cystatins (WCs) with twelve rice, seven barley, one sorghum and ten corn cystatin sequences employing different bioinformatics tools. The obtained results identified highly conserved signature sequences in all the cystatins considered. Several other motifs were also identified, based on which the sequences could be categorized into groups in congruence with the phylogenetic clustering. Homology modeling of WCs revealed 3D structural topology so well shared by other cystatins. Protein-protein interaction of WCs with papain supported the notion that functional diversity is a con-sequence of existing differences in amino acid residues in highly conserved as well as relatively less conserved motifs. Thus there is a significant conservation at the sequential and structural levels; however, concomitant variations maintain the functional diversity in this protein family, which constantly modulates itself to reciprocate the diversity while counteracting the cysteine proteinases.

Basal expression studies of cystatins during specific growth stages of wheat spikes for defining their possible role in differential and stage dependent immunity against Karnal bunt (Tilletia indica).

Two genotypes showing differential immunity against Karnal bunt (Tilletia indica) were used to investigate the role of three members of cystatin gene family in growth stage dependent immunity in wheat (Triticum aestivum L.). Three members of cystatin gene family (WC1, WC2, and WC4) were cloned and sequenced. Analysis of sequenced data showed that there was 76-99% nucleotide and protein sequence identity between different genes of the wheat cystatin. In silico amino acid sequence analysis revealed the presence of a conserved signature pattern of residues and also the functional domains were presumed to be actively involved in imparting cysteine protease inhibition capability. The semi-quantitative and quantitative levels of these members were measured by means of RT-PCR, northern blotting, western blotting, and by ELISA techniques. The members of cystatin gene family were expressed in both resistant (HD 29) and susceptible genotypes (WH 542); however, the expression level was significantly (P < 0.001) higher in resistant compared to susceptible genotype at all the stages of wheat spikes. The patterns of expression of WC2, WC4 were similar except in the levels in S(1) and S(2) stages as it remained constant (P > 0.05) in contrary to WC1 family whose expression gradually increased from S(v) to S(2) stage. According to the intensity of the detected band in RT PCR, northern blot and western blot, WC1 family seems to be expressed more than the other gene families. The immunoassay results further showed that WC1 protein was abundantly expressed in resistant genotype and high expression was observed at the S2 stage as compared to susceptible genotype (P < 0.001) suggesting that low level of expression of WC1 in S2 stage is responsible for KB infection. The results of the present study clearly indicate the role of cystatin gene family in differential and stage dependent immunity against KB.

Exploration of Csp genes from temperate and glacier soils of the Indian Himalayas and in silico analysis of encoding proteins.

The metagenomic Csp library was constructed from the temperate and glacier soils of central Himalaya, India followed by polymerase chain reaction (PCR) amplification. The library was further screened for low-temperature adaptation, and the positive recombinants were sorted out by determining changes in the melting temperature (Tm). A homology search of cloned sequence showed their identity with the Csp genes of Pseudomonas fluorescens, Psychrobacter cryohalolentis K5, and Shewanella spp MR-4. Amino acid sequence analysis annotated the presence of conserved aromatic and basic amino acids as well as RNA binding motifs from the cold shock domain. Furthermore, a PROSITE scan showed a moderate identity of less than 60% with the known cold shock-inducible proteins (ribosomal proteins, rbfA, DEAD-box helicases), cold acclimation protein, and temperature-induced protein (SRP1/TIP1). This study highlighted the prevalence of Csp genes from cold Himalayan environments that can be explored for tailor-made crop constructions in future.