Molecular mapping of drought-responsive QTLs during the reproductive stage of rice using a GBS (genotyping-by-sequencing) based SNP linkage map.

BACKGROUND: 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. METHODS AND RESULTS: 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 F6 and F7 RILs were performed under drought stress and control conditions. A total of 42 QTLs were identified under drought stress and control conditions for yield component traits explaining 1.95-13.36% of the total phenotypic variance (PVE). Among these, 19 QTLs were identified under drought stress conditions, whereas 23 QTLs were located under control conditions. A total of 4 QTLs explained a PVE ≥ 10% which are considered as the major QTLs. Moreover, bioinformatics analysis revealed the presence of 6 candidate genes, which showed differential expression under drought and control conditions. CONCLUSION: These QTLs/genes may be deployed for marker-assisted pyramiding to improve drought tolerance in the existing rice varieties.

Identification and analysis of miRNAs-lncRNAs-mRNAs modules involved in stem-elongation of deepwater rice (Oryza sativa L.).

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.

In silico characterization and differential expression analysis of 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR) of Centella asiatica.

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.

Genome-wide association studies using 50 K rice genic SNP chip unveil genetic architecture for anaerobic germination of deep-water rice population of Assam, India.

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.

Revealing shared differential co-expression profiles in rice infected by virus from reoviridae and sequiviridae group.

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.

Identification and validation of plant miRNA from NGS data-an experimental approach.

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.

An Insight Into Structure, Function, and Expression Analysis of 3-Hydroxy-3-Methylglutaryl-CoA Reductase of Cymbopogon winterianus.

Citronella (Cymbopogon winterianus) is one of the richest sources of high-value isoprenoid aromatic compounds used as flavour, fragrance, and therapeutic elements. These isoprenoid compounds are synthesized by 2 independent pathways: mevalonate pathway and 2-C-methyl-d-erythritol-4-phosphate pathway. Evidence suggests that 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) is a rate-controlling enzyme for the synthesis of variety of isoprenoids. This study reports the isolation, characterization, and tissue-specific expression analysis of HMGR from citronella. The modelled HMGR is a class I type of HMGR enzyme with 3-domain architecture. The active site comprises a cofactor (nicotinamide adenine dinucleotide phosphate) and the substrate-binding motifs. The real-time and quantitative reverse transcription-polymerase chain reaction results revealed equal expression level in both leaf sheath and root tissue. The results from our study shall be a valuable resource for future molecular intervention to alter the metabolic flux towards improvement of key active ingredient in this important medicinal plant.

Bruchid egg induced transcript dynamics in developing seeds of black gram (Vigna mungo).

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.

Genome wide transcriptome profiling reveals differential gene expression in secondary metabolite pathway of Cymbopogon winterianus.

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.

Molecular cloning, characterization and expression analysis of 3-hydroxy-3-methylglutaryl coenzyme A reductase gene from Centella asiatica L.

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.

Novel insights into structure-function mechanism and tissue-specific expression profiling of full-length dxr gene from Cymbopogon winterianus.

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.

Molecular recognition of avirulence protein (avrxa5) by eukaryotic transcription factor xa5 of rice (Oryza sativa L.): insights from molecular dynamics simulations.

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.

E-microsatellite markers for Centella asiatica (Gotu Kola) genome: validation and cross-transferability in Apiaceae family for plant omics research and development.

Abstract Centella asiatica (Gotu Kola) is a plant that grows in tropical swampy regions of the world and has important medicinal and culinary use. It is often considered as part of Ayurvedic medicine, traditional African medicine, and traditional Chinese medicine. The unavailability of genomics resources is significantly impeding its genetic improvement. To date, no attempt has been made to develop Expressed Sequence Tags (ESTs) derived Simple Sequence Repeat (SSR) markers (eSSRs) from the Centella genome. Hence, the present study aimed to develop eSSRs and their further experimental validation and cross-transferability of these markers in different genera of the Apiaceae family to which Centella belongs. An in-house pipeline was developed for the entire analyses by combining bioinformatics tools and perl scripts. A total of 4443 C. asiatica EST sequences from dbEST were processed, which generated 2617 nonredundant high quality EST sequences consisting 441 contigs and 2176 singletons. Out of 1776.5 kb of examined sequences, 417 (15.9%) ESTs containing 686 SSRs were detected with a density of one SSR per 2.59 kb. The gene ontology study revealed 282 functional domains involved in various processes, components, and functions, out of which 64 ESTs were found to have both SSRs and functional domains. Out of 603 designed EST-SSR primers, 18 pairs of primers were selected for validation based on the optimum parameter value. Reproducible amplification was obtained for six primer pairs in C. asiatica that were further tested for cross-transferability in nine other important genera/species of the Apiaceae family. Cross-transferability of the EST-SSR markers among the species were examined and Centella javanica showed highest transferability (83.3%). The study revealed six highly polymorphic EST-SSR primers with an average PIC value of 0.95. In conclusion, these EST-SSR markers hold a big promise for the genomics analysis of Centella asiatica, to facilitate comparative map-based analyses across other related species within the Apiaceae family, and future marker-assisted breeding programs. To the best of our knowledge, this is the first report of development of EST-SSRs in Centella asiatica by in silico approaches, which offers a veritable potential in further use in plant omics research and development.

Structural comparison, substrate specificity, and inhibitor binding of AGPase small subunit from monocot and dicot: present insight and future potential.

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.

Rediscovering medicinal plants' potential with OMICS: microsatellite survey in expressed sequence tags of eleven traditional plants with potent antidiabetic properties.

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.

Structure-based computational study of two disease resistance gene homologues (Hm1 and Hm2) in maize (Zea mays L.) with implications in plant-pathogen interactions.

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.

Computational identification and characterization of conserved miRNAs and their target genes in garlic (Allium sativum L.) expressed sequence tags.

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.

Structural analysis and molecular dynamics simulations of novel δ-endotoxin Cry1Id from Bacillus thuringiensis to pave the way for development of novel fusion proteins against insect pests of crops.

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.

In silico identification and characterization of conserved miRNAs and their target genes in sweet potato (Ipomoea batatas L.) expressed sequence tags (ESTs).

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.

Insights into the structure-function relationship of disease resistance protein HCTR in maize (Zea mays L.): a computational structural biology approach.

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.