Search

In rice, drought stress at reproductive stage drastically reduces yield, which in turn hampers farmer's efforts towards crop production. The majority of the rice varieties have resistance genes against several abiotic and biotic stresses. Therefore, the traditional landraces were studied to identify QTLs/candidate genes associated with drought tolerance.A high-density SNP-based genetic map was constructed using a Genotyping-by-sequencing (GBS) approach. The recombinant inbred lines (RILs) derived from crossing 'Banglami × Ranjit' were used for QTL analysis. A total map length of 1306.424 cM was constructed, which had an average inter-marker distance of 0.281 cM. The phenotypic evaluation of FThese QTLs/genes may be deployed for marker-assisted pyramiding to improve drought tolerance in the existing rice varieties.

Deepwater is an abiotic stress that limits rice cultivation worldwide due to recurrent floods. The miRNAs and lncRNAs are two non-coding RNAs emerging as major regulators of gene expressions under different abiotic stresses. However, the regulation of these two non-coding RNAs under deepwater stress in rice is still unexplored. In this study, small RNA-seq and RNA-seq from internode and node tissues were analyzed to predict deepwater stress responsive miRNAs and lncRNAs, respectively. Additionally, a competitive endogenous RNA (ceRNA) study revealed about 69 and 25 lncRNAs acting as endogenous target mimics (eTM) with the internode and node miRNAs, respectively. In ceRNA analyses, some of the key miRNAs such as miR1850.1, miR1848, and IN-nov-miR145 were upregulated while miR159e was downregulated, and their respective eTM lncRNAs and targets were found to have opposite expressions. Moreover, we have transiently expressed one module (IN-nov-miR145-Cc-TCONS_00011544-Os11g36430.3) in tobacco leaves. The integrated analysis has identified differentially expressed (DE) miRNAs, lncRNAs and their target genes, and the complex regulatory network, which might lead to stem elongation under deepwater stress. In this novel attempt to identify and characterize miRNAs and lncRNAs under deepwater stress in rice, we have provided, probably for the first time, a reference platform to study the interactions of these two non-coding RNAs with respective target genes through transient expression analyses.

Centella asiatica (L.) is one of the most valuable medicinal plants since prehistoric times. The pharmaceutical importance of this herb is due to the accumulation of large quantities of pentacyclic triterpenoid saponins, collectively known as centelloids synthesized by the isoprenoid biosynthesis pathway. Biosynthesis of triterpenoid in the plants proceeds via either of the two pathways, viz. Mevalonate (MVA) pathway (in the cytosol) or 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway (in plastid). In Centella, the pathway leading to the accumulation of triterpenoid is still not known or elucidated. Thus, to know whether the MVA or MEP pathway or a cross-talk between the pathway leads to the biosynthesis of triterpenoid, silencing the key regulatory gene using RNAi tool, of each of the pathway and then analyze a metabolite is an efficient approach. The key regulatory enzyme of the MVA pathway i.e. 3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) has already been successfully silenced using RNAi tool [1]. In the present study, the 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) a key regulatory enzyme in MEP pathway is silenced. The RNAi-DXR construct in pHANNIBAL vector was cloned into a binary vector pART27 and subsequently transformed into Agrobacterium strain AGL1. The transient analysis of the RNAi-CaDXR using semi-quantitative RT-PCR confirmed the silencing of the endogenous DXR gene in Nicotiana and further confirmed in Centella asiatica. The present study is the first step aimed to delineate the MEP pathway using RNAi silencing approach to elucidate its role in the accumulation of triterpenoid in this important medicinal plant.

The 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR; EC1.1.1.267), an NADPH-dependent reductase, plays a pivotal role in the methylerythritol 4-phosphate pathway (MEP), in the conversion of 1-deoxy-d-xylulose-5-phosphate (DXP) into MEP. Photochemical profiles, as well as pharmaceutical activities of Centella asiatica (L.), one of the most valuable medicinal plants, divulge the presence of secondary metabolites called Centellosides. Despite well-studied pharmaceutical activities, not much is known about the genes responsible for the synthesis of these compounds. In the present study, the full-length DXR gene sequence (JQ965955) of Centella submitted in NCBI was characterized using various bioinformatics tools and tissue specific differential expression studies were also carried out. The full-length CDNA of CaDXR contains an open reading frame (ORF) of 1425 bp which encodes a peptide of 474 amino acids. The molecular weight of this protein was found to be 51.5 kDa with isoelectric point of 6.33. The protein contains three conserved domain, namely NADPH (GSTGSIGT and LAAGSNV), substrate binding (LPADSEHSAI and NKGLEVIEAHY) and Cys-Ser-(Ala/Met/Val/Thr) cleavage-site domains. Phylogenetic studies of CaDXR sequence show close homology with DXR sequence of Angelica sinensis and Daucus carota subsp sativus as they all belong to Apiaceae family. In silico analysis predicted that CaDXR protein contains 21 α-helix and 11 β-sheets and further DXR protein model was validated by Ramachandran plot analysis. The results of molecular dynamics (MD) simulations unveil dynamic stability of the proposed model and docking studies suggest that the NDP cofactor tightly binds in the active site of the protein with a strong network of hydrogen and hydrophobic interactions. The expression studies by semi-RT followed by qRT—PCR suggests that CaDXR is differentially expressed in different tissues (with maximal expression in the node and lowest in the roots). Thus, characterization and structure–function analysis of DXR gene in Centella facilitate us to understand not only the functions of DXR gene but also regulatory mechanisms involved in the MEP pathway in C. asiatica plant at the molecular level.

Post-genomic era of biology is faced with a major challenge in deciphering the gene function out of enormous amount of data generated by NGS technology, in addition to large number of EST sequences, available in the public domain. Among many approaches, one important approach is to knock out the gene and analysis of the visible effect of loss of gene function. RNA interference (RNAi), as a reverse genetic approach currently in use for studies of gene function, holds a great promise in this context. The efficacy of RNAi as tool for functional genomics study has already been successfully demonstrated in Caenorhabditis elegans. With the availability of comprehensive resources of genomic sequence data and knowledge of the biological mechanism of RNAi, use of RNAi in functional genomics is quickly gaining space and popularity. Results from the study of transgene-induced RNAi suggests many variables that should be considered while designing experiment to decipher gene function.

North Eastern part of India such as Assam is inundated by flood every year where the farmers are forced to grow the traditional tall deep-water rice. Genetic improvement of this type of rice is slow because of insufficient knowledge about their genetic architecture and population structure. In the present investigation, the genetic diversity architecture of 94 deep-water rice genotypes of Assam and association mapping strategy was, for the first time, applied to determine the significant SNPs and genes for deep-water rice. These genotypes are known for their unique elongation ability under deep-water condition. The anaerobic germination (AG) related trait-associated genes identified here can provide affluent resources for rice breeding especially in flood-prone areas. We investigated the genome-wide association studies (GWAS) using 50 K rice genic SNP chip across 94 deep-water rice genotypes collected from different flood-prone districts/villages of Assam. Population structure and diversity analysis revealed that these genotypes were stratified into four sub-populations. Using GWAS approach, 20 significant genes were identified and found to be associated with AG-related traits. Of them, two most relevant genes (OsXDH1and SSXT) have been identified which explain phenotypic variability (R2 > 20%) in the population. These genes were located in Chr 3 (LOC_Os03g31550) which encodes for enzyme xanthine dehydrogenase 1(OsXDH1) and in Chr 12 (LOC_Os12g31350) which encodes for SSXT family protein. Both of these genes were found to be associated with anaerobic response index (increase in the coleoptile length under water in anaerobic condition with respect to control), respectively. Interestingly, OsXDH1is involved in purine catabolism pathway and acts as a scavenger of reactive oxygen species in plants, whereas SSXT is GRF1-interacting factor 3. These two candidate genes associated with AG of deep-water rice have been found to be reported for the first time. Thus, this study provides a greater resource for breeders not only for improvement of deep-water rice, but also for AG tolerant variety useful for direct-seeded rice in flood-affected areas.

Being the most important mode of propagation of plants, disease-free seed is a basic requirement of progressive agriculture. However, many plant pathogens such as bacteria, fungi, virus, virus-like agents, nematodes, etc. are transmitted through seeds to the next generation. Seed transmission enables the earliest possible infection, which is a determining factor of disease severity. Plant viruses are the most important plant pathogens causing up to 96% yield losses. Although, seed transmissibility is presently known to occur for around 18% of plant viruses only, it is estimated that one-third of the plant viruses will eventually prove to be seed transmitted in at least one host. Nevertheless, the most likely mechanism of plant viruses’ persistence between crop seasons is through the seeds, and even a seed transmission rate as low as 0.001 has the potential to initiate an epidemic. Therefore, study of the plant-pathogen interactions in seed-transmitted viruses is of utmost importance to help formulate preventive measures against such pathogens. Although the specific molecular mechanism of transmission of plant viruses through seeds is not completely understood till date, broadly, the mode of infection of the embryos (the future seeds) is known to be either direct (via mother plant) or indirect (via infected pollen); these mechanisms, however, may not be mutually exclusive. The modes of viral movement from infected maternal and paternal tissues to the embryo, the genetics of host-virus interaction, etc. are needed to be worked out for designing successful management strategies against seed-transmitted viruses.

miRNAs are class of endogenously initiated noncoding RNAs, which are most critical in gene expression and regulation at posttranscriptional level. They do so either by cleavage of the target mRNA or by translational repression. miRNAs are being given enough attention in recent years because of its role in myriad developmental processes including tumorogenesis and host-pathogen interaction. Advent of Next Generation Sequencing (NGS) technology and computational approach made it possible to pinpoint the precise role of miRNA and their target. Identification of miRNAs and their target has several approaches depending on efficiency, cost and time. The present review summarizes the developments in the field of plant miRNA w.r.t. to experimental approaches that are being followed to identify and validate the miRNAs and their targets.

Centella asiatica is an important medicinal plant and synthesizes the triterpenoid saponins via the isoprenoid pathway. The metabolic engineering approach is a promising tool to improve the key component which is limiting in many of the medicinal plants but, prior knowledge of the pathway leading to accumulation of the secondary metabolite is prerequisite before any molecular intervention. The present study aimed to knock-down 3-hydroxy-3-methylglutaryl-CoA reductases (HMGR), a key regulatory enzyme in the mevalonic (MVA) pathway using RNAi approach. Long term goal of this project is to delineate the biochemical pathway responsible for synthesis of key active ingredient in Centella. 3-Hydroxy-3-methylglutaryl-CoA reductases (HMGR) play an important role in catalyzing the first committed step of isoprenoid biosynthesis in the mevalonate (MVA) pathway in plants. The RNAi-HMGR construct was designed using the pHANNIBAL vector and subsequently cloned into a binary vector pART27. The transient analysis of the RNAi-CaHMGR using semi quantitative RT-PCR confirmed the silencing of the endogenous HMGR gene in Centella asiatica. Thus, the designing of the RNAi-HMGR construct and transient expression analysis of the construct gene is the first step to elucidate triterpenoid biosynthesis pathway that shall pave the way to intervene following metabolic engineering approach in future.

The 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR; EC1.1.1.267), an NADPH-dependent reductase, plays a pivotal role in the methylerythritol 4-phosphate pathway (MEP), in the conversion of 1-deoxy-d-xylulose-5-phosphate (DXP) into MEP. Photochemical profiles, as well as pharmaceutical activities of Centella asiatica (L.), one of the most valuable medicinal plants, divulge the presence of secondary metabolites called Centellosides. Despite well-studied pharmaceutical activities, not much is known about the genes responsible for the synthesis of these compounds. In the present study, the full-length DXR gene sequence (JQ965955) of Centella submitted in NCBI was characterized using various bioinformatics tools and tissue specific differential expression studies were also carried out. The full-length CDNA of CaDXR contains an open reading frame (ORF) of 1425 bp which encodes a peptide of 474 amino acids. The molecular weight of this protein was found to be 51.5 kDa with isoelectric point of 6.33. The protein contains three conserved domain, namely NADPH (GSTGSIGT and LAAGSNV), substrate binding (LPADSEHSAI and NKGLEVIEAHY) and Cys-Ser-(Ala/Met/Val/Thr) cleavage-site domains. Phylogenetic studies of CaDXR sequence show close homology with DXR sequence of Angelica sinensis and Daucus carota subsp sativus as they all belong to Apiaceae family. In silico analysis predicted that CaDXR protein contains 21 α-helix and 11 β-sheets and further DXR protein model was validated by Ramachandran plot analysis. The results of molecular dynamics (MD) simulations unveil dynamic stability of the proposed model and docking studies suggest that the NDP cofactor tightly binds in the active site of the protein with a strong network of hydrogen and hydrophobic interactions. The expression studies by semi-RT followed by qRT—PCR suggests that CaDXR is differentially expressed in different tissues (with maximal expression in the node and lowest in the roots). Thus, characterization and structure–function analysis of DXR gene in Centella facilitate us to understand not only the functions of DXR gene but also regulatory mechanisms involved in the MEP pathway in C. asiatica plant at the molecular level. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-021-02723-w.

Rice production is adversely constrained by drought especially if encountered during panicle initiation and flowering stage. Several traditional rice cultivars/landraces have inherent functioning mechanism of adaptation in diverse environmental adversities. The genetic and molecular characterization of these traditional tolerant germplasm at the omics level is critical to understand the schemas and molecular networks regulating stress tolerance and finally to identify candidate genes involved in stress tolerance. The current study attempts to unravel miRNA in response to drought stress in ARC-10372, a drought tolerant important rice genotype of Assam Rice Collection. Small RNA sequencing was performed that resulted a total of 45,742,638 good quality reads in both drought stressed and control samples. Upon annotation, 178 and 131 known miRNA were found to be expressed in control and drought samples, respectively. Out of these, 109 known miRNA belonging to 58 families were found to be differentially expressed. Nine (9) earlier reported drought responsive miRNAs were found to be expressed significantly in ARC-10372 and their expression pattern is similar to the earlier studies and, 10 miRNAs (miR531, miR827, miR8175 miR977, miR6300, miR1861, miR440, miR9773, miR3982 and miR1876) which were not reported earlier as drought responsive were found to be regulated significantly under drought stress in this local landrace. In silico identification revealed that target of most of these miRNAs are transcription factors such as NAC, MYB, MYC,ARF thereby further confirming their role in regulation during stress adaptation. Four new stress responsive miRNA targeted stress related transcript. The expression pattern of miRNA in sRNA sequencing data were validated by qPCR and most of them exhibited concordance with NGS data. In light of the reported studies, the present study for the first time gives a comparative aspect of drought regulated miRNome of leaf at panicle initiation stage in drought tolerant rice cultivar ARC-10372.

Differential co-expression is a cutting-edge approach to analyze gene expression data and identify both shared and divergent expression patterns. The availability of high-throughput gene expression datasets and efficient computational approaches have unfolded the opportunity to a systems level understanding of functional genomics of different stresses with respect to plants. We performed the meta-analysis of the available microarray data for reoviridae and sequiviridae infection in rice with the aim to identify the shared gene co-expression profile. The microarray data were downloaded from ArrayExpress and analyzed through a modified Weighted Gene Co-expression Network Analysis (WGCNA) protocol. WGCNA clustered the genes based on the expression intensities across the samples followed by identification of modules, eigengenes, principal components, topology overlap, module membership and module preservation. The module preservation analysis identified 4 modules; salmon (638 genes), midnightblue (584 genes), lightcyan (686 genes) and red (562 genes), which are highly preserved in both the cases. The networks in case of reoviridae infection showed neatly packed clusters whereas, in sequiviridae, the clusters were loosely connected which is due to the differences in the correlation values. We also identified 83 common transcription factors targeting the hub genes from all the identified modules. This study provides a coherent view of the comparative aspect of the expression of common genes involved in different virus infections which may aid in the identification of novel targets and development of new intervention strategy against the virus.

3-Hydroxy-3-methylglutaryl-CoA reductases (HMGR) plays an important role in catalyzing the first committed step of isoprenoid biosynthesis in the mevelonic (MVA) pathway (catalyzes the conversion of HMG-CoA to MVA) in plants. The present manuscript reports the full length cDNA cloning of HMGR (CaHMGR, GenBank accession number: KJ939450.2) and its characterization from Centella asiatica. Sequence analysis indicated that the cDNA was of 1965 bp, which had an open reading frame of 1617 bp and encoded a protein containing 539 amino-acids with a mol wt of 57.9 kDa. A BLASTp search against non-redundant (nr) protein sequence showed that C. asiatica HMGR (CaHMGR) has 65-81% identity with HMGRs from different plant species and multi-alignment comparison analysis showed the presence of two motif each corresponding to HMG-CoA-binding and NADP(H)-binding. The Conserved Domain Database analysis predicted that CaHMGR belongs to Class I hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase. Three-dimensional modeling confirmed the novelty of CaHMGR with a spatial structure similar to Homo sapiens (PDB id: 1IDQ8_A). Tissue Expression analysis indicates that CaHMGR is ubiquitous albeit differentially expressed among different tissues analysed, Strong expression was recorded in the nodes and leaves and low in the roots. The present investigation confirmed that nodes are vital to terpenoid synthesis in C. asiatica. Thus, the cloning of full length CDS, characterization and structure-function analysis of HMGR gene in Centella facilitate to understand the HMGR's functions and regulatory mechanisms involved in mevalonate pathway in C. asiatica at genetic level.

Highlights • We sequenced the complete coding sequence of citronella dxr gene. • Analysis of CwDXR revealed the plastidial nature of the enzyme. • Comparative modelling of CwDXR was performed. • Molecular dynamics simulation and docking explained the mode of cofactor recognition. • Tissue-specific expression of CwDXR was performed using RT-PCR and qRT-PCR., The 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR; EC1.1.1.267), an NADPH-dependent reductase, plays a pivotal role in the methylerythritol 4-phosphate pathway (MEP), in the conversion of 1-deoxy-d-xylulose-5-phosphate (DXP) into MEP. The sheath and leaf of citronella (Cymbopogon winterianus) accumulates large amount of terpenes and sesquiterpenes with proven medicinal value and economic uses. Thus, sequencing of full length dxr gene and its characterization seems to be a valuable resource in metabolic engineering to alter the flux of isoprenoid active ingredients in plants. In this study, full length DXR from citronella was characterized through in silico and tissue-specific expression studies to explain its structure–function mechanism, mode of cofactor recognition and differential expression. The modelled DXR has a three-domain architecture and its active site comprised of a cofactor (NADPH) binding pocket and the substrate-binding pocket. Molecular dynamics simulation studies indicated that DXR model retained most of its secondary structure during 10 ns simulation in aqueous solution. The modelled DXR superimposes well with its closest structural homolog but subtle variations in the charge distribution over the cofactor recognition site were noticed. Molecular docking study revealed critical residues aiding tight anchoring NADPH within the active pocket of DXR. Tissue-specific differential expression analysis using semi-quantitative RT-PCR and qRT-PCR in various tissues of citronella plant revealed distinct differential expression of DXR. To our knowledge, this is the first ever report on DXR from the important medicinal plant citronella and further characterization of this gene will open up better avenues for metabolic engineering of secondary metabolite pathway genes from medicinal plants in the near future.

Black gram (Vigna mungo) seeds are a rich source of digestible proteins, however, during storage these seeds are severely damaged by bruchids (Callosobruchus spp.), reducing seed quality and yield losses. Most of the cultivated genotypes of black gram are susceptible to bruchids, however, few tolerant genotypes have also been identified but the mechanism of tolerance is poorly understood. We employed Suppression Subtractive Hybridization (SSH) to identify specifically, but rarely expressed bruchid egg induced genes in black gram. In this study, Suppression Subtractive Hybridization (SSH) library was constructed to study the genes involved in defense response in black gram against bruchid infestation. An EST library of 277 clones was obtained for further analyses. Based on CAP3 assembly, 134 unigenes were computationally annotated using Blast2GOPRO software. In all, 20 defense related genes were subject to quantitative PCR analysis (qPCR) out of which 12 genes showed up-regulation in developing seeds of the pods oviposited by bruchids. Few major defense genes like defensin, pathogenesis related protein (PR), lipoxygenase (LOX) showed high expression levels in the oviposited population when compared with the non-oviposited plants. This is the first report on defense related gene transcript dynamics during the bruchid-black gram interaction using SSH library. This library would be useful to clone defense related gene(s) such as defensin as represented in our library for crop improvement.

Herbal medicines and traditionally used medicinal plants present an untapped potential for novel molecular target discovery using systems science and OMICS biotechnology driven strategies. Since up to 40% of the world's poor people have no access to government health services, traditional and folk medicines are often the only therapeutics available to them. In this vein, North East (NE) India is recognized for its rich bioresources. As part of the Indo-Burma hotspot, it is regarded as an epicenter of biodiversity for several plants having myriad traditional uses, including medicinal use. However, the improvement of these valuable bioresources through molecular breeding strategies, for example, using genic microsatellites or Simple Sequence Repeats (SSRs) or Expressed Sequence Tags (ESTs)-derived SSRs has not been fully utilized in large scale to date. In this study, we identified a total of 47,700 microsatellites from 109,609 ESTs of 11 medicinal plants (pineapple, papaya, noyontara, bitter orange, bermuda brass, ratalu, barbados nut, mango, mulberry, lotus, and guduchi) having proven antidiabetic properties. A total of 58,159 primer pairs were designed for the non-redundant 8060 SSR-positive ESTs and putative functions were assigned to 4483 unique contigs. Among the identified microsatellites, excluding mononucleotide repeats, di-/trinucleotides are predominant, among which repeat motifs of AG/CT and AAG/CTT were most abundant. Similarity search of SSR containing ESTs and antidiabetic gene sequences revealed 11 microsatellites linked to antidiabetic genes in five plants. GO term enrichment analysis revealed a total of 80 enriched GO terms widely distributed in 53 biological processes, 17 molecular functions, and 10 cellular components associated with the 11 markers. The present study therefore provides concrete insights into the frequency and distribution of SSRs in important medicinal resources. The microsatellite markers reported here markedly add to the genetic stock for cross transferability in these plants and the literature on biomarkers and novel drug discovery for common chronic diseases such as diabetes.

The theoretical three-dimensional structure of a novel δ-endotoxin Cry1Id (81 kDa) belonging to Cry1I class, toxic to many of the lepidopteran pests has been investigated through comparative modeling. Molecular dynamics (MD) simulations was carried out to characterize its structural and dynamical features at 10 ns in explicit solvent using the GROMACS version 4.5.4. Finally the simulated model was validated by the SAVES, WHAT IF, MetaMQAP, ProQ, ModFOLD and MolProbity servers. Despite low sequence identity with its structural homologs, Cry1Id not only resembles the previously reported Cry structures but also shares the common five conserved blocks of amino acid residues. Although the domain II of Cry1Id superpose well with its closest structural homolog Cry8Ea1, variation of amino acids and length in the apical loop2 of domain II was observed. In this work, we have hypothesized that the variations in apical loop2 might be the sole factor for providing variable surface accessibility to Cry1Id protein that could be important in receptor recognition. MD simulation showed the proposed endotoxin retains its stable conformation in aqueous solution. The result from this study is expected to aid in the development hybrid Cry proteins and new potent fusion proteins with novel specificities against different insect pests for improved pest management of crop plants.

The disease resistance gene Hm1 of maize encodes a NADPH-dependent reductase enzyme, HC-toxin reductase (HCTR) that detoxifies the HC toxin secreted by the race specific fungus Cochliobolus carbonum race 1. HCTR enzyme shares 29.6% sequence identity with dihydroflavonol reductase (DFR) of grape, a key enzyme involved in flavonoid biosynthesis. Here we report the comparative modelling, molecular dynamics simulation and docking studies to explain the structure–function relationship and the mode of cofactor (NADPH) binding in HCTR enzyme at the molecular level. The nucleotide binding domain of modelled HCTR adopts a classic Rossmann fold and possesses a consensus glycine rich GxGxxG motif. Molecular simulation studies suggested that HCTR model retained stability throughout the simulation in aqueous solution. HCTR model showed considerable structural identities with the cofactor binding site of DFR, but significant difference in the catalytic site might be the reason of functional divergence between these families of proteins. Similarly electrostatic surface potential analysis of both HCTR and DFR revealed profound variations in the charge distribution over the substrate binding site, which can be correlated with the sequence variability and may suggest distinct substrate-binding patterns and differences in the catalytic mechanism. Docking results indicated Phe19, Gly21, Arg40, Thr90, Gly208, Arg218, Glu221 and Thr222 are important residues for cofactor (NADPH) binding through strong hydrogen bonding and electrostatic interactions. Alanine scanning and analysis of docking energies of mutant proteins suggested that Phe19, and Arg40 are two critical residues for the cofactor binding. The result from the present study is expected to pave the way for exploration of similar genes in other economically important crop varieties.

Actin-1 is a ubiquitous protein belonging to the reproductive class of Actin family in Arabidopsis thaliana. Thisprotein is involved in the formation of filaments that are major components of the cytoskeleton. Despite the im-portance of this protein, very little information is availa ble regarding its structure and function in plants. In thisstudy, analysis of the protein sequence was done an d comparative model of Actin-1 was constructed (UNIPROTID: P0CJ46) from Arabidopsis thaliana using the crystal structure of Dictyostelium discoideum actin (PDB ID:1NLV-A) as template employing Modelle r version 9.9. The stable structure was generated by 5 nanosecond mole-cular dynamics simulation steps usin g GROMOS43A1 96 force field that characterized its structural and dynamicfeature. The biochemical function of the final simulated structure was also investigated using PROFUNC. Themolecular simulation study suggested that the modeled Acti n-1 protein retain its stable conformation in aqueoussolution. The predicted binding sites in the modeled Actin- 1 protein are very informative for further protein-ligandinteraction study.

Tea is the most popular non-alcoholic and healthy beverage across the world. The understanding of the genetic organization and molecular biology of tea plant, which is very poorly understood at present, is required for quantum increase in productivity and efficient use of germplasm for either cultivation or breeding program. Single-pass sequencing of randomly selected cDNA clones is the most widely accepted technique for gene identification and cloning. In the present study, a good quality cDNA library was constructed and preliminary analysis of ESTs was carried out. The titers of unamplified and amplified libraries were 1.4 × 106 pfu/ml and 5.27 × 108 pfu/ml respectively. A total of 210 cDNA clones from the constructed cDNA library were sequenced and analyzed. A total of 84 high quality Expressed Sequence Tags (ESTs) were generated, among which 71 ESTs had significant homology with sequences in NCBI non-redundant protein database by BLAST X analysis. About 80% ESTs had poly (A) tail at 3′ end indicating that the cDNAs were full length. The database-matched ESTs were classified into putative cellular roles, viz. energy-related category (corresponding to 20% of total BLAST X matched ESTs), Transcription (14.2%), protein synthesis (14.2%) cell growth and division (8.6%), cell structure (5.7%), signal transduction (5.7%), transporters (2.9%), disease and defenses (2.9%), secondary metabolism (2.9%) and gene regulation (2.9%). This study provides an overview of the mRNA expression profile and first hand information of gene sequence expressed in tender leaves and apical buds of tea plant.

Superoxide dismutases (SODs), members of the metalloenzymes family are most effective intracellular enzymatic antioxidant in aerobic organisms. These enzymes provide the first line of defense in plants against the toxic effects of elevated levels of reactive oxygen species (ROS) generated during various environmental stresses. The availability of high-throughput computational tools has provided better opportunities to characterize the protein features and determine their function. In the present study an attempt was made to gain an insight into the structure and evolution of subunits of SODs (Cu-Zn, Mn and Fe SODs) of rice. The 3-Dimensional structures of SODs were modeled based on available X-ray crystal structures and further validated. The primary sequence, secondary and tertiary structure analysis revealed Mn and Fe SOD to be structurally homologous while Cu-Zn SOD is un-related to either of them. Comparative structural study also revealed former two were dominated by α-helices followed by β-strands in contrast; Cu-Zn SOD dominated by β-strands. Molecular phylogeny indicated a common evolutionary origin of Mn and Fe SOD while Cu-Zn SOD may have evolved separately.

Simple Sequence Repeats (SSRs) developed from Expressed Sequence Tags (ESTs), known as EST-SSRs are most widely used and potentially valuable source of gene based markers for their high levels of crosstaxon portability, rapid and less expensive development. The EST sequence information in the publicly available databases is increasing in a faster rate. The emerging computational approach provides a better alternative process of development of SSR markers from the ESTs than the conventional methods. In the present study, 12,851 EST sequences of Camellia sinensis, downloaded from National Center for Biotechnology Information (NCBI) were mined for the development of Microsatellites. 6148 (4779 singletons and 1369 contigs) non redundant EST sequences were found after preprocessing and assembly of these sequences using various computational tools. Out of total 3822.68 kb sequence examined, 1636 (26.61%) EST sequences containing 2371 SSRs were detected with a density of 1 SSR/1.61 kb leading to development of 245 primer pairs. These mined EST-SSR markers will help further in the study of variability, mapping, evolutionary relationship in Camellia sinensis. In addition, these developed SSRs can also be applied for various studies across species.

Advances in transcriptome sequencing provide fast, cost-effective and reliable approach to generate large expression datasets especially suitable for non-model species to identify putative genes, key pathway and regulatory mechanism. Citronella (Cymbopogon winterianus) is an aromatic medicinal grass used for anti-tumoral, antibacterial, anti-fungal, antiviral, detoxifying and natural insect repellent properties. Despite of having number of utilities, the genes involved in terpenes biosynthetic pathway is not yet clearly elucidated. The present study is a pioneering attempt to generate an exhaustive molecular information of secondary metabolite pathway and to increase genomic resources in Citronella. Using high-throughput RNA-Seq technology, root and leaf transcriptome was analysed at an unprecedented depth (11.7 Gb). Targeted searches identified majority of the genes associated with metabolic pathway and other natural product pathway viz. antibiotics synthesis along with many novel genes. Terpenoid biosynthesis genes comparative expression results were validated for 15 unigenes by RT-PCR and qRT-PCR. Thus the coverage of these transcriptome is comprehensive enough to discover all known genes of major metabolic pathways. This transcriptome dataset can serve as important public information for gene expression, genomics and function genomics studies in Citronella and shall act as a benchmark for future improvement of the crop.

ADP-glucose pyrophosphorylase (AGPase) is the first rate limiting enzyme of starch biosynthesis pathway and has been exploited as the target for greater starch yield in several plants. The structure-function analysis and substrate binding specificity of AGPase have provided enormous potential for understanding the role of specific amino acid or motifs responsible for allosteric regulation and catalytic mechanisms, which facilitate the engineering of AGPases. We report the three-dimensional structure, substrate, and inhibitor binding specificity of AGPase small subunit from different monocot and dicot crop plants. Both monocot and dicot subunits were found to exploit similar interactions with the substrate and inhibitor molecule as in the case of their closest homologue potato tuber AGPase small subunit. Comparative sequence and structural analysis followed by molecular docking and electrostatic surface potential analysis reveal that rearrangements of secondary structure elements, substrate, and inhibitor binding residues are strongly conserved and follow common folding pattern and orientation within monocot and dicot displaying a similar mode of allosteric regulation and catalytic mechanism. The results from this study along with site-directed mutagenesis complemented by molecular dynamics simulation will shed more light on increasing the starch content of crop plants to ensure the food security worldwide.

The avirulence gene avrxa5 of bacterial blight pathogen Xanthomonas oryzae pv. oryzae (Xoo) recognized by the resistant rice lines having corresponding resistance (xa5) gene in a gene–for–gene manner. We used a combinatorial approach involving protein–protein docking, molecular dynamics (MD) simulations and binding free energy calculations to gain novel insights into the gene-for-gene mechanism that governs the direct interaction of R-Avr protein. From the best three binding poses predicted by molecular docking, MD simulations were performed to explore the dynamic binding mechanism of xa5 and avrxa5. Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) techniques were employed to calculate the binding free energy and to uncover the thriving force behind the molecular recognition of avrxa5 by eukaryotic transcription factor xa5. Binding free energy analysis revealed van der Waals term as the most constructive component that favors the xa5 and avrxa5 interaction. In addition, hydrogen bonds (H-bonds) and essential electrostatic interactions analysis highlighted amino acid residues Lys54/Asp870, Lys56/Ala868, Lys56/Ala866, Lys56/Glu871, Ile59/His862, Gly61/Phe858, His62/Arg841, His62/Leu856, Ser101/Ala872 and Ser105/Asp870 plays pivotal role for the energetically stability of the R-Avr complex. Insights gained from the present study are expected to unveil the molecular mechanisms that define the transcriptional activator mediated transcriptome modification in host plants.

The NADPH-dependent HC-toxin reductases (HCTR1 and 2) encoded by enzymatic class of disease resistance homologous genes (Hm1 and Hm2) protect maize by detoxifying a cyclic tetrapeptide, HC-toxin, secreted by the fungus Cochliobolus carbonum race 1(CCR1). Unlike the other classes' resistance (R) genes, HCTR-mediated disease resistance is an inimitable mechanism where the avirulence (Avr) component from CCR1 is not involved in toxin degradation. In this study, we attempted to decipher cofactor (NADPH) recognition and mode of HC-toxin binding to HCTRs through molecular docking, molecular dynamics (MD) simulations and binding free energy calculation methods. The rationality and the stability of docked complexes were validated by 30-ns MD simulation. The binding free energy decomposition of enzyme-cofactor complex was calculated to find the driving force behind cofactor recognition. The overall binding free energies of HCTR1-NADPH and HCTR2-NADPH were found to be −616.989 and −16.9749 kJ mol−1 respectively. The binding free energy decomposition revealed that the binding of NADPH to the HCTR1 is mainly governed by van der Waals and nonpolar interactions, whereas electrostatic terms play dominant role in stabilizing the binding mode between HCTR2 and NADPH. Further, docking analysis of HC-toxin with HCTR-NADPH complexes showed a distinct mode of binding and the complexes were stabilized by a strong network of hydrogen bond and hydrophobic interactions. This study is the first in silico attempt to unravel the biophysical and biochemical basis of cofactor recognition in enzymatic class of R genes in cereal crop maize.

The endogenous small non-coding micro RNAs (miRNAs), which are typically ~21–24 nt nucleotides, play a crucial role in regulating the intrinsic normal growth of cells and development of the plants as well as in maintaining the integrity of genomes. These small non-coding RNAs function as the universal specificity factors in post-transcriptional gene silencing. Discovering miRNAs, identifying their targets, and further inferring miRNA functions is a routine process to understand normal biological processes of miRNAs and their roles in the development of plants. Comparative genomics based approach using expressed sequence tags (EST) and genome survey sequences (GSS) offer a cost-effective platform for identification and characterization of miRNAs and their target genes in plants. Despite the fact that sweet potato (Ipomoea batatas L.) is an important staple food source for poor small farmers throughout the world, the role of miRNA in various developmental processes remains largely unknown. In this paper, we report the computational identification of miRNAs and their target genes in sweet potato from their ESTs. Using comparative genomics-based approach, 8 potential miRNA candidates belonging to miR168, miR2911, and miR156 families were identified from 23 406 ESTs in sweet potato. A total of 42 target genes were predicted and their probable functions were illustrated. Most of the newly identified miRNAs target transcription factors as well as genes involved in plant growth and development, signal transduction, metabolism, defense, and stress response. The identification of miRNAs and their targets is expected to accelerate the pace of miRNA discovery, leading to an improved understanding of the role of miRNA in development and physiology of sweet potato, as well as stress response.

Rice (Oryza sativa L.), a model plant belonging to the family Poaceae, is a staple food for a majority of the people worldwide. Grown in the tropical and subtropical regions of the world, this important cereal crop is under constant and serious threat from both biotic and abiotic stresses. Among the biotic threats, Xanthomonas oryzae pv. oryzae, causing the damaging bacterial blight disease in rice, is a prominent pathogen. The xa5 gene in the host plant rice confers race-specific resistance to this pathogen. This recessive gene belongs to the Xa gene family of rice and encodes a gamma subunit of transcription factor IIA (TFIIAγ). In view of the importance of this gene in conferring resistance to the devastating disease, we reconstructed the phylogenetic relationship of this gene, developed a three-dimensional protein model, followed by long-term molecular dynamics simulation studies to gain a better understanding of the evolution, structure, and function of xa5. The modeled structure was found to fit well with the small subunit of TFIIA from human, suggesting that it may also act as a small subunit of TFIIA in rice. The model had a stable conformation in response to the atomic flexibility and interaction, when subjected to MD simulation at 20 nano second in aqueous solution. Further structural analysis of xa5 indicated that the protein retained its basic transcription factor function, suggesting that it might govern a novel pathway responsible for bacterial blight resistance. Future molecular docking studies of xa5 underway with its corresponding avirulence gene is expected to shed more direct light into plant-pathogen interactions at the molecular level and thus pave the way for richer agriproteomic insights.

The endogenous small non-coding functional microRNAs (miRNAs) are short in size, range from ~21 to 24 nucleotides in length, play a pivotal role in gene expression in plants and animals by silencing genes either by destructing or blocking of translation of homologous mRNA. Although various high-throughput, time consuming and expensive techniques like forward genetics and direct cloning are employed to detect miRNAs in plants but comparative genomics complemented with novel bioinformatic tools pave the way for efficient and cost-effective identification of miRNAs through homologous sequence search with previously known miRNAs. In this study, an attempt was made to identify and characterize conserved miRNAs in garlic expressed sequence tags (ESTs) through computational means. For identification of novel miRNAs in garlic, a total 3227 known mature miRNAs of plant kingdom Viridiplantae were searched for homology against 21,637 EST sequences resulting in identification of 6 potential miRNA candidates belonging to 6 different miRNA families. The psRNATarget server predicted 33 potential target genes and their probable functions for the six identified miRNA families in garlic. Most of the garlic miRNA target genes seem to encode transcription factors as well as genes involved in stress response, metabolism, plant growth and development. The results from the present study will shed more light on the understanding of molecular mechanisms of miRNA in garlic which may aid in the development of novel and precise techniques to understand some post-transcriptional gene silencing mechanism in response to stress tolerance.

Glycoside hydrolase family 19 chitinases (EC 3.2.1.14) widely distributed in plants, bacteria and viruses catalyse the hydrolysis of chitin and play a major role in plant defense mechanisms and development. Rice possesses several classes of chitinase, out of which a single structure of class I has been reported in PDB to date. In the present study an attempt was made to gain more insight into the structure, function and evolution of class I, II and IV chitinases of GH family 19 from rice. The three-dimensional structures of chitinases were modelled and validated based on available X-ray crystal structures. The structural study revealed that they are highly α-helical and bilobed in nature. These enzymes are single or multi domain and multi-functional in which chitin-binding domain (CBD) and catalytic domain (CatD) are present in class I and IV whereas class II lacks CBD. The CatD possesses a catalytic triad which is thought to be involved in catalytic process. Loop III, which is common in all three classes of chitinases, reflects that it may play a significant role in their function. Our study also confirms that the absence and presence of different loops in GH family 19 of rice may be responsible for various sized products. Molecular phylogeny revealed chitinases in monocotyledons and dicotyledons differed from each other forming two different clusters and may have evolved differentially. More structural study of this enzyme from different plants is required to enhance the knowledge of catalytic mechanism and substrate binding.

UNLABELLED With the advent of high-throughput sequencing technology, sequences from many genomes are being deposited to public databases at a brisk rate. Open access to large amount of expressed sequence tag (EST) data in the public databases has provided a powerful platform for simple sequence repeat (SSR) development in species where sequence information is not available. SSRs are markers of choice for their high reproducibility, abundant polymorphism and high inter-specific transferability. The mining of SSRs from ESTs requires different high-throughput computational tools that need to be executed individually which are computationally intensive and time consuming. To reduce the time lag and to streamline the cumbersome process of SSR mining from ESTs, we have developed a user-friendly, web-based EST-SSR pipeline "EST-SSR-MARKER PIPELINE (ESMP)". This pipeline integrates EST pre-processing, clustering, assembly and subsequently mining of SSRs from assembled EST sequences. The mining of SSRs from ESTs provides valuable information on the abundance of SSRs in ESTs and will facilitate the development of markers for genetic analysis and related applications such as marker-assisted breeding. AVAILABILITY The database is available for free at http://bioinfo.aau.ac.in/ESMP.

Microsatellites are the markers of choice due to their high abundance reproducibility, degree of polymorphism and co-dominant nature. These are mainly used for studying the genetic variability in different species and Marker assisted selection. Expressed Sequence Tags (ESTs) serve as the main resource for Simple Sequence Repeats (SSRs). The computational approach for detecting SSRs and developing SSR markers from EST-SSRs is preferred over the conventional methods as it reduces time and cost to a great extent. The available EST sequence databases, various web interfaces and standalone tools provide the platform for an easy analysis of the EST sequences leading to the development of potential EST-SSR Markers. This paper is an overview of in silico approach to develop SSR Markers from the EST sequence using some of the most efficient tools that are available freely for academic purpose.