Identification and Characterization of Circular Intronic RNAs Derived from Insulin Gene.

Circular RNAs (circRNAs) are a large family of noncoding RNAs that have emerged as novel regulators of gene expression. However, little is known about the function of circRNAs in pancreatic ß-cells. Here, transcriptomic analysis of mice pancreatic islet RNA-sequencing data identified 77 differentially expressed circRNAs between mice fed with a normal diet and a high-fat diet. Surprisingly, multiple circRNAs were derived from the intron 2 of the preproinsulin 2 (Ins2) gene and are termed as circular intronic (ci)-Ins2. The expression of ci-Ins2 transcripts in mouse pancreatic islets, and ßTC6 cells were confirmed by reverse transcription PCR, DNA sequencing, and RNase R treatment experiments. The level of ci-Ins2 was altered in ßTC6 cells upon exposure to elevated levels of palmitate and glucose. Computational analysis predicted the interaction of several RNA-binding proteins with ci-Ins2 and their flanking region, suggesting their role in the ci-Ins2 function or biogenesis. Additionally, bioinformatics analysis predicted the association of several microRNAs with ci-Ins2. Gene ontology and pathway analysis of genes targeted by miRNAs associated with ci-Ins2 suggested the regulation of several key biological processes. Together, our findings indicate that differential expression of circRNAs, especially ci-Ins2 transcripts, may regulate ß-cell function and may play a critical role in the development of diabetes.

Modulation of specific cell cycle phases in human embryonic stem cells by lncRNA RNA decoys.

Recent studies have shown that long noncoding RNAs (lncRNAs) are crucial regulators of human embryonic stem cells (hESCs). However, modes of actions of lncRNAs in hESCs are not well illustrated. Here, we predicted a regulatory network in hESCs in which lncRNAs interact with TFs and thereby control the expressions of downstream targets of TFs. The predicted network is comprised of 2289 3-motif subgraphs which are characterized by 3 nodes: (i) a lncRNA which is predicted to interact with (ii) a TF and (iii) a gene which is a target of TF and coexpressing with lncRNA. We performed functional annotation of the network by identifying hub nodes followed by pathway enrichment study, which unveiled an active G1-S cell cycle phase transition-specific subnetwork that encompasses 2 lncRNAs, MALAT1 and DANCR. Our analysis revealed that MALAT1 and DANCR might be playing key roles in G1-S phase transition by acting as RNA decoy via interacting with crucial stemness maintaining TFs. We predicted that MALAT1 possibly compete with DNMT1 and CDCA7 genes to bind to E2F1 thereby interrupting repression of DNMT1 and activation of CDCA7 by E2F1 in hESCs, whereas DANCR possibly competes with IPO7 gene to bind to MYC thereby interrupting MYC-mediated activation of IPO7 in hESCs. Both of these are conjectured to contribute to rapid G1-S phase transition aiding in stemness maintenance of hESCs. This study presents a crucial TF target cross talks mediated by lncRNAs in hESCs regulating its properties which needs further investigation.

Deciphering synergistic regulatory networks of microRNAs in hESCs and fibroblasts.

MicroRNAs (miRNAs) play crucial roles in pluripotency and differentiation of human Embryonic Stem Cells (hESCs). However, synergism among multiple miRNAs and their regulatory effects on stem cells are largely unknown. We investigated the synergistic regulations among miRNAs, which are differentially expressed in hESCs and fibroblasts (Fibs) to gain deeper insights into the regulatory mechanisms of miRNA-miRNA synergism in the differentiation of hESCs into Fibs. In this study, we identified miRNA targets incorporating differential expression profiles of miRNAs and genes in hESCs-Fibs with miRNA targeting features followed by enrichment analysis. We then built the active miRNA-miRNA synergistic network (MMSN) integrating miRNA-miRNA synergistic pairs, which identified 16 miRNA-miRNA closed communities. Further, topology assessment, community overlapping and functional studies revealed hsa-miR-873-5p as an important pluripotent miRNA which might be hindering hESCs differentiation into Fibs by targeting regulators of epithelial to mesenchymal transition (EMT). On the other hand, synergism among hsa-miR-98-5p, hsa-let-7i-5p, hsa-let-7g-5p, hsa-miR-30a-3p and hsa-miR-29b-2-5p is predicted to be highly essential for promoting differentiation of hESCs in to Fibs. This study deepens our understanding of miRNA synergism and its regulatory impacts on properties of stem cells.

The smaller heparin fragments bind non-specifically through the IAPP sequence: An in silico study.

The diminishing ß-cell mass of pancreas in type II diabetes mellitus (TIIDM) is intricately linked with high fibrillation propensity of islet amyloid polypeptide (IAPP, aka amylin). IAPP is one of the most amyloidogenic peptide secreted by pancreatic ß-cells. In the autopsy of TIIDM patients, IAPP rich amyloid plaques are found containing different components of extracellular matrix (ECM), including heparin. For a positively charged IAPP, interaction with heparin which has accessible high density negatively charged functional groups is anticipated to moderate the fibrillation kinetics. Hence, the heparin has shown to affect the amyloidogenicity and cytotoxicity of IAPP depending upon its polymer length; short polymer inhibited the amyloidogenicity and longer fragment enhanced the propensity. Here using docking and molecular dynamic (MD) simulations studies, the work investigates key interactions between IAPP and different heparin fragments, those are likely involved in moderating IAPP fibrillation kinetics in presence of different length heparin fragments. The findings indicate that the heparin fragments of longer length, >dp7, predominantly interact with IAPP N- and C-termini, resulting in a stable complex with solvent accessible self-recognition element (SRE) of IAPP sequence. However, shorter fragment non-specifically binds through the IAPP sequence, including N-terminus residues and SRE sequences.

Functional efficiency of MerA protein among diverse mercury resistant bacteria for efficient use in bioremediation of inorganic mercury.

MerA protein of mer operon in mercury resistant bacteria influences transformation of Hg2+ to Hg0. Both in-silico and in-vivo studies have been carried out and MerA sequences, conserved motifs for mercury binding and NADPH (GCVPSK and LSCCA) varied widely in both Gram-positive and Gram-negative bacteria. As MerA-NADPH-FAD complex plays an important role in mercury volatilization, molecular interaction studies between MerA, NADPH, FAD and Hg2+ was carried out to study the efficiency of transformation of Hg2+ to Hg0 in mercury resistant bacteria. After the prediction of suitable models and molecular interaction analysis, the potential energies in the selected bacteria were as follows: Bacillus thuringiensis (NADPH: -5.15 kcal/mol and FAD: -9.63 kcal/mol), Pseudomonas aeruginosa (NADPH: -3.8 kcal/mol and FAD: -8.56 kcal/mol), Exiguobacterium sp. (NADPH: -3.37 kcal/mol and FAD: -8.42 kcal/mol), Vibrio sp. (NADPH: -3.3 kcal/mol and FAD: -6.7 kcal/mol) and Escherichia coli (NADPH: -3.28 kcal/mol and FAD: -5.69 kcal/mol). Additionally, the binding scores between MerA and Hg2+ followed the similar trend and found higher in B. thuringiensis (3.79) followed by P. aeruginosa (3.57), Exiguobacterium sp. (2.37), Vibrio sp. (1.47) and E. coli (1.07). ANOVA (2-way) result showed the significant (P < 0.05) variation among the energy values obtained after interaction studies. In-vivo analysis of expression of merA gene and Hg2+ removal efficiency also followed the same pattern with a highly significant correlation (P < 0.001) between the binding energy, binding score and expression pattern of merA gene as well as Hg2+ volatilization. Thus, the mercury removal efficiency of bacteria is genera specific which is correlated with the binding efficiency between MerA-NADPH complex and Hg2+ in mer operon mediated mercury resistant bacteria.

An integrative approach predicted co-expression sub-networks regulating properties of stem cells and differentiation.

The differentiation of human Embryonic Stem Cells (hESCs) is accompanied by the formation of different intermediary cells, gradually losing its stemness and acquiring differentiation. The precise mechanisms underlying hESCs integrity and its differentiation into fibroblast (Fib) are still elusive. Here, we aimed to assess important genes and co-expression sub-networks responsible for stemness, early differentiation of hESCs into embryoid bodies (EBs) and its lineage specification into Fibs. To achieve this, we compared transcriptional profiles of hESCs-EBs and EBs-Fibs and obtained differentially expressed genes (DEGs) exclusive to hESCs-EBs (early differentiation), EBs-Fibs (late differentiation) and common DEGs in hESCs-EBs and EBs-Fibs. Then, we performed gene set enrichment analysis (GSEA) followed by overrepresentation study and identified key genes for each gene category. The regulations of these genes were studied by integrating ChIP-Seq data of core transcription factors (TFs) and histone methylation marks in hESCs. Finally, we identified co-expression sub-networks from key genes of each gene category using k-clique sub-network extraction method. Our study predicted seven genes edicting core stemness properties forming a co-expression network. From the pathway analysis of sub-networks of hESCs-EBs, we hypothesize that FGF2 is contributing to pluripotent transcription network of hESCs in association with DNMT3B and JARID2 thereby facilitating cell proliferation. On the contrary, FGF2 is found to promote cell migration in Fibs along with DDR2, CAV1, DAB2, and PARVA. Moreover, our study identified three k-clique sub-networks regulating TGF-ß signaling pathway thereby promoting EBs to Fibs differentiation by: (i) modulating extracellular matrix involving ITGB1, TGFB1I1 and GBP1, (ii) regulating cell cycle remodeling involving CDKN1A, JUNB and DUSP1 and (iii) helping in epithelial to mesenchymal transition (EMT) involving THBS1, INHBA and LOX. This study put forward the unexplored genes and co-expression sub-networks regulating stemness and different stages of differentiation of hESCs which will undoubtedly add to the comprehensive understanding of hESCs biology.

Unraveling novel TF-miRNA regulatory crosstalk in metastasis of Soft Tissue Sarcoma.

Cancer metastasis is a disease of extreme clinical relevance, as it is responsible for more than 90% of cancer-associated mortality. The molecular mechanism and critical regulators involved in this complex multi-stage process of metastasis is poorly deciphered in soft tissue sarcomas (STS), a heterogeneous group of rare tumors with high metastatic potential. Therefore, we aimed at identifying miRNA and transcription factor (TF) regulatory networks and paths in STS metastasis. We integrated mRNA and miRNA expression profiles with curated regulations (TF→gene, TF→miRNA, miRNA→gene) from different databases and constructed a potentially active regulatory sub-network in STS metastasis. From functional and topological analysis, we found nine novel regulators of Notch signaling sub-network which are conjectured to play critical role in metastasis of STS. This illustrated that the sub-network is promising for identification of critical regulators. Further analysis deploying our developed tool 'RiNAcyc' and computing coverage ratio of known STS associated genes and miRNAs identified a 15 node active path. This potential path highlights the crucial role of BMP2, hsa-miR-24, AP2 and MYC as the up-stream regulators of the path and hsa-miR-215 and TYMS as potential indicator of chemotherapeutic benefit in STS metastasis.

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.

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.

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.

Molecular phylogeny, homology modeling, and molecular dynamics simulation of race-specific bacterial blight disease resistance protein (xa5) of rice: a comparative agriproteomics approach.

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.

A comparative proteomic approach to analyse structure, function and evolution of rice chitinases: a step towards increasing plant fungal resistance.

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.

Mining for SSRs and FDMs from expressed sequence tags of Camellia sinensis.

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.

An approach to delineate primers for a group of poorly conserved sequences incorporating the common motif region.

Glutathione synthetase (gshB) has previously been reported to confer tolerance to acidic soil condition in Rhizobium species. Cloning the gene coding for this enzyme necessitates the designing of proper primer sets which in turn depends on the identification of high quality sequence similarity in multiple global alignments. In this experiment, a group of homologous gene sequences related to gshB gene (accession no: gi-86355669:327589-328536) of Rhizobium etli CFN 42, were extracted from NCBI nucleotide sequence databases using BLASTN and were analyzed for designing degenerate primers. However, the T-coffee multiple global alignment results did not show any block of conserved region for the above sequence set to design the primers. Therefore, we attempted to identify the location of common motif region based on multiple local alignments employing the MEME algorithm supported with MAST and Primer3. The results revealed some common motif regions that enabled us to design the primer sets for related gshB gene sequences. The result will be validated in wet lab.

ESMP: A high-throughput computational pipeline for mining SSR markers from ESTs.

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.