HGT-Gen: a tool for generating a phylogenetic tree with horizontal gene transfer

Horizontal gene transfer (HGT) is a common event in prokaryotic evolution. Therefore, it is very important to consider HGT in the study of molecular evolution of prokaryotes. This is true also for conducting computer simulations of their molecular phylogeny because HGT is known to be a serious disturbing factor for estimating their correct phylogeny. To the best of our knowledge, no existing computer program has generated a phylogenetic tree with HGT from an original phylogenetic tree. We developed a program called HGT-Gen that generates a phylogenetic tree with HGT on the basis of an original phylogenetic tree of a protein or gene. HGT-Gen converts an operational taxonomic unit or a clade from one place to another in a given phylogenetic tree. We have also devised an algorithm to compute the average length between any pair of branches in the tree. It defines and computes the relative evolutionary time to normalize evolutionary time for each lineage. The algorithm can generate an HGT between a pair of donor and acceptor lineages at the same evolutionary time. HGT-Gen is used with a sequence-generating program to evaluate the influence of HGT on the molecular phylogeny of prokaryotes in a computer simulation study.

Availability

The database is available for free at http://www.grl.shizuoka.ac.jp/˜thoriike/HGT-Gen.html     

GHOSTX: An Improved Sequence Homology Search Algorithm Using a Query Suffix Array and a Database Suffix Array

DNA sequences are translated into protein coding sequences and then further assigned to protein families in metagenomic analyses, because of the need for sensitivity. However, huge amounts of sequence data create the problem that even general homology search analyses using BLASTX become difficult in terms of computational cost. We designed a new homology search algorithm that finds seed sequences based on the suffix arrays of a query and a database, and have implemented it as GHOSTX. GHOSTX achieved approximately 131–165 times acceleration over a BLASTX search at similar levels of sensitivity. GHOSTX is distributed under the BSD 2-clause license and is available for download at http://www.bi.cs.titech.ac.jp/ghostx/. Currently, sequencing technology continues to improve, and sequencers are increasingly producing larger and larger quantities of data. This explosion of sequence data makes computational analysis with contemporary tools more difficult. We offer this tool as a potential solution to this problem.     

Mammalian Mitochondrial ncRNA Database

Mammalian Mitochondrial ncRNA is a web-based database, which provides specific information on non-coding RNA in mammals. This database includes easy searching, comparing with BLAST and retrieving information on predicted structure and its function about mammalian ncRNAs.

Availability

The database is available for free at http://www.iitm.ac.in/bioinfo/mmndb/     

PRD: A protein�CRNA interaction database

Although protein–RNA interactions (PRIs) are involved in various important cellular processes, compiled data on PRIs are still limited. This contrasts with protein–protein interactions, which have been intensively recorded in public databases and subjected to network level analysis. Here, we introduce PRD, an online database of PRIs, dispersed across several sources, including scientific literature. Currently, over 10,000 interactions have been stored in PRD using PSI-MI 2.5, which is a standard model for describing detailed molecular interactions, with an emphasis on gene level data. Users can browse all recorded interactions and execute flexible keyword searches against the database via a web interface. Our database is not only a reference of PRIs, but will also be a valuable resource for studying characteristics of PRI networks.

Availability

PRD can be freely accessed at http://pri.hgc.jp/     

A Novel Cholinergic Action of Alcohol and the Development of Tolerance to That Effect in Caenorhabditis elegans

The Prp43 DExD/H-box protein is required for progression of the biochemically distinct pre-messenger RNA and ribosomal RNA (rRNA) maturation pathways. In Saccharomyces cerevisiae, the Spp382/Ntr1, Sqs1/Pfa1, and Pxr1/Gno1 proteins are implicated as cofactors necessary for Prp43 helicase activation during spliceosome dissociation (Spp382) and rRNA processing (Sqs1 and Pxr1). While otherwise dissimilar in primary sequence, these Prp43-binding proteins each contain a short glycine-rich G-patch motif required for function and thought to act in protein or nucleic acid recognition. Here yeast two-hybrid, domain-swap, and site-directed mutagenesis approaches are used to investigate G-patch domain activity and portability. Our results reveal that the Spp382, Sqs1, and Pxr1 G-patches differ in Prp43 two-hybrid response and in the ability to reconstitute the Spp382 and Pxr1 RNA processing factors. G-patch protein reconstitution did not correlate with the apparent strength of the Prp43 two-hybrid response, suggesting that this domain has function beyond that of a Prp43 tether. Indeed, while critical for Pxr1 activity, the Pxr1 G-patch appears to contribute little to the yeast two-hybrid interaction. Conversely, deletion of the primary Prp43 binding site within Pxr1 (amino acids 102–149) does not impede rRNA processing but affects small nucleolar RNA (snoRNA) biogenesis, resulting in the accumulation of slightly extended forms of select snoRNAs, a phenotype unexpectedly shared by the prp43 loss-of-function mutant. These and related observations reveal differences in how the Spp382, Sqs1, and Pxr1 proteins interact with Prp43 and provide evidence linking G-patch identity with pathway-specific DExD/H-box helicase activity.     

Larval host records of butterflies in Japan

Using Japanese literature, we created a consolidated list of host records of butterflies in Japan. The list used the host records described in eight major illustrated reference books, two checklists, and 14 other pieces of literature. The presence of larvae on plants, the observation of larvae eating plants or insects in the field were considered as host records. We collected all species recorded in Japan. Scientific, family, and Japanese names of butterflies were consolidated using the BINRAN database (http://binran.lepimages.jp/). Scientific and Japanese names of host plants were based on the YList database (http://ylist.info/). If scientific names of host plants were not found in YList, we used scientific names based on The Plant List (http://www.theplantlist.org/). Family names of host plants were based on the Catalogue of Life database (http://www.catalogueoflife.org/). Scientific, family, and Japanese names of host insects were based on the MOKUROKU database (http://konchudb.agr.agr.kyushu-u.ac.jp/mokuroku/) for Hymenoptera and the catalogue of the Paraneoptera of Japan published by the Entomological Society of Japan for Hemiptera. We also provided the references of each host record and the original names described in the referred literature. Two datasets, HostDB and ReferenceDB, were created to include 3600 records of butterfly larval hosts in Japan, along with scientific and Japanese names of each species and a literature list. These datasets will be useful for basic and applied biological studies of butterflies. Data files are stored in the Ecological Research Data Archives (http://db.cger.nies.go.jp/JaLTER/ER_DataPapers/) and available from http://hostbj.lepumus.net/. These datasets are published under the Creative Commons License Attribution-ShareAlike 4.0 (CC BY-SA, https://creativecommons.org/licenses/by-sa/4.0/).     

A grid layout algorithm for automatic drawing of biochemical networks

MOTIVATION: Visualization is indispensable in the research of complex biochemical networks. Available graph layout algorithms are not adequate for satisfactorily drawing such networks. New methods are required to visualize automatically the topological architectures and facilitate the understanding of the functions of the networks. RESULTS: We propose a novel layout algorithm to draw complex biochemical networks. A network is modeled as a system of interacting nodes on squared grids. A discrete cost function between each node pair is designed based on the topological relation and the geometric positions of the two nodes. The layouts are produced by minimizing the total cost. We design a fast algorithm to minimize the discrete cost function, by which candidate layouts can be produced efficiently. A simulated annealing procedure is used to choose better candidates. Our algorithm demonstrates its ability to exhibit cluster structures clearly in relatively compact layout areas without any prior knowledge. We developed Windows software to implement the algorithm for CADLIVE. AVAILABILITY: All materials can be freely downloaded from http://kurata21.bio.kyutech.ac.jp/grid/grid_layout.htm; http://www.cadlive.jp/ SUPPLEMENTARY INFORMATION: http://kurata21.bio.kyutech.ac.jp/grid/grid_layout.htm; http://www.cadlive.jp/     

bPE toolkit: toolkit for computational protein engineering

We present a computational toolkit consisting of five utility tools, for performing basic operations on a protein structure file in PDB format. The toolkit consists of five different programs which can be integrated as part of a pipeline for computational protein structure characterization or as a standalone analysis package. The programs include tools for chirality check for amino acids (ProChiral), contact map generation (CoMa), data redundancy (DaRe), hydrogen bond potential energy (HyPE) and electrostatic interaction energy (EsInE). All programs in the toolkit can be accessed and downloaded through the following link: http://www.iitg.ac.in/bpetoolkit/.     

CAB-Align: A Flexible Protein Structure Alignment Method Based on the Residue-Residue Contact Area

Proteins are flexible, and this flexibility has an essential functional role. Flexibility can be observed in loop regions, rearrangements between secondary structure elements, and conformational changes between entire domains. However, most protein structure alignment methods treat protein structures as rigid bodies. Thus, these methods fail to identify the equivalences of residue pairs in regions with flexibility. In this study, we considered that the evolutionary relationship between proteins corresponds directly to the residue–residue physical contacts rather than the three-dimensional (3D) coordinates of proteins. Thus, we developed a new protein structure alignment method, contact area-based alignment (CAB-align), which uses the residue–residue contact area to identify regions of similarity. The main purpose of CAB-align is to identify homologous relationships at the residue level between related protein structures. The CAB-align procedure comprises two main steps: First, a rigid-body alignment method based on local and global 3D structure superposition is employed to generate a sufficient number of initial alignments. Then, iterative dynamic programming is executed to find the optimal alignment. We evaluated the performance and advantages of CAB-align based on four main points: (1) agreement with the gold standard alignment, (2) alignment quality based on an evolutionary relationship without 3D coordinate superposition, (3) consistency of the multiple alignments, and (4) classification agreement with the gold standard classification. Comparisons of CAB-align with other state-of-the-art protein structure alignment methods (TM-align, FATCAT, and DaliLite) using our benchmark dataset showed that CAB-align performed robustly in obtaining high-quality alignments and generating consistent multiple alignments with high coverage and accuracy rates, and it performed extremely well when discriminating between homologous and nonhomologous pairs of proteins in both single and multi-domain comparisons. The CAB-align software is freely available to academic users as stand-alone software at http://www.pharm.kitasato-u.ac.jp/bmd/bmd/Publications.html.     

VMD DisRg: New User-Friendly Implement for calculation distance and radius of gyration in VMD program

Molecular dynamic simulation is a practical and powerful technique for analysis of protein structure. Several programs have been developed to facilitate the mentioned investigation, under them the visual molecular dynamic or VMD is the most frequently used programs. One of the beneficial properties of the VMD is its ability to be extendable by designing new plug-in. We introduce here a new facility of the VMD for distance analysis and radius of gyration of biopolymers such as protein and DNA.

Availability

The database is available for free at http://trc.ajums.ac.ir/HomePage.aspx/?TabID/=12618/&Site/=trc.ajums.ac/&Lang/=fa-IR     

Comparative Analysis of Pdf-Mediated Circadian Behaviors Between Drosophila melanogaster and D. virilis

PAR proteins (partitioning defective) are major regulators of cell polarity and asymmetric cell division. One of the par genes, par-1, encodes a Ser/Thr kinase that is conserved from yeast to mammals. In Caenorhabditis elegans, par-1 governs asymmetric cell division by ensuring the polar distribution of cell fate determinants. However the precise mechanisms by which PAR-1 regulates asymmetric cell division in C. elegans remain to be elucidated. We performed a genomewide RNAi screen and identified six genes that specifically suppress the embryonic lethal phenotype associated with mutations in par-1. One of these suppressors is mpk-1, the C. elegans homolog of the conserved mitogen activated protein (MAP) kinase ERK. Loss of function of mpk-1 restored embryonic viability, asynchronous cell divisions, the asymmetric distribution of cell fate specification markers, and the distribution of PAR-1 protein in par-1 mutant embryos, indicating that this genetic interaction is functionally relevant for embryonic development. Furthermore, disrupting the function of other components of the MAPK signaling pathway resulted in suppression of par-1 embryonic lethality. Our data therefore indicates that MAP kinase signaling antagonizes PAR-1 signaling during early C. elegans embryonic polarization.ASYMMETRIC cell division, a process in which a mother cell divides in two different daughter cells, is a fundamental mechanism to achieve cell diversity during development. We use the early embryo of Caenorhabditis elegans as a model system to study asymmetric cell division. The C. elegans one-cell embryo divides asymmetrically along its anteroposterior axis, generating two cells of different sizes and fates: the larger anterior daughter cell will generate somatic tissues while the smaller posterior daughter cell will generate the germline (Sulston et al. 1983).A group of proteins called PAR proteins (partitioning defective) is required for asymmetric cell division in C. elegans (Kemphues et al. 1988). Depletion of any of the seven par genes (par-1 to -6 and pkc-3) leads to defects in asymmetric cell division and embryonic lethality (Kemphues et al. 1988; Kirby et al. 1990; Tabuse et al. 1998; Hung and Kemphues 1999; Hao et al. 2006). PAR-3 and PAR-6 are conserved proteins that contain PDZ-domains and form a complex with PKC-3 (Etemad-Moghadam et al. 1995; Izumi et al. 1998; Tabuse et al. 1998; Hung and Kemphues 1999). This complex becomes restricted to the anterior cortex of the embryo in response to spatially defined actomyosin contractions occurring in the embryo upon fertilization (Goldstein and Hird 1996; Munro et al. 2004). The posterior cortex of the embryo that becomes devoid of the anterior PAR proteins is occupied by the RING protein PAR-2 and the Ser/Thr kinase PAR-1 (Guo and Kemphues 1995; Boyd et al. 1996; Cuenca et al. 2003). Once polarized, the anterior and posterior PAR proteins mutually exclude each other from their respective cortices (Etemad-Moghadam et al. 1995; Boyd et al. 1996; Cuenca et al. 2003; Hao et al. 2006). Loss of function of the gene par-1, as opposed to loss of most other par genes, results in embryos that exhibit only subtle effects on the polarized cortical domains occupied by the other PAR proteins (Cuenca et al. 2003). However defects in this gene are associated with a more symmetric division in size, an aberrant distribution of cell fate specification markers, altered cell fates of the daughter cells of the embryo, and ultimately embryonic lethality (Kemphues et al. 1988; Guo and Kemphues 1995).PAR-1 controls asymmetric cell division and cell fate specification by regulating the localization of the two highly similar CCCH-type zinc-finger proteins MEX-5 and MEX-6 (referred to as MEX-5/6). MEX-5 and MEX-6 are 70% identical in their amino acid sequence and fulfill partially redundant functions in the embryo (Schubert et al. 2000). In wild-type animals, endogenous MEX-5 and GFP fusions of MEX-6 localize primarily to the anterior of the embryo while both proteins are evenly distributed in par-1 mutant embryos (Schubert et al. 2000; Cuenca et al. 2003). This suggests that in wild-type animals, PAR-1 acts in part by restricting MEX-5 and MEX-6 to the anterior of the embryo. The precise mechanism of this regulation is not known, but an elegant study performed for MEX-5 indicates that differential protein mobility in the anterior and posterior cytoplasm of the one-cell embryo contributes to this asymmetry (Tenlen et al. 2008). While increased mobility in the posterior of the one-cell embryo correlates with a par-1- and par-4-dependent phosphorylation on MEX-5, the kinase directly phosphorylating MEX-5 remains to be identified (Tenlen et al. 2008).Some of the phenotypes associated with loss of par-1 function are dependent on the function of mex-5 and mex-6. First, loss of function of par-1 leads to a decreased stability and aberrant localization of the posterior cell fate specification marker PIE-1, a protein that is usually inherited by the posterior daughter cell in wild-type animals and ensures the correct specification of the germline (Mello et al. 1996; Seydoux et al. 1996). This decreased stability is dependent on mex-5/6 function as PIE-1 levels are restored, albeit with symmetrical distribution, in mex-6(RNAi); mex-5(RNAi); par-1(b274) embryos (Schubert et al. 2000; Cuenca et al. 2003; Derenzo et al. 2003). Second, embryos lacking par-1 function exhibit decreased amounts of P granules in the one-cell embryo, while these markers are present in mex-6(pk440); mex-5(zu199); par-1(RNAi) embryos of comparable age (Cheeks et al. 2004). Third, in par-1(RNAi) one-cell embryos the posterior cortical domain occupied by the polarity protein PAR-2 is extended anteriorly, when compared to wild-type embryos (Cuenca et al. 2003). This anterior extension is rescued in embryos deficient for both par-1 and mex-5/6 (Cuenca et al. 2003). Taken together, these results indicate that par-1 acts in the embryo—at least in part—by regulating the localization and/or activity of the proteins MEX-5 and MEX-6. However, it remains unclear whether other proteins can modulate PAR-1 function to affect MEX-5/6 activity.To gain insight into the mechanisms of par-1 function in the embryo, we sought to identify genes that act together with par-1 during embryonic development. We performed an RNAi-based screen for genetic interactors of the temperature-sensitive allele par-1(zu310), using the embryonic lethal phenotype of this mutant as a readout. This method has proven successful in previous screens to identify genes involved in early embryonic processes (Labbé et al. 2006; O''Rourke et al. 2007). We were able to identify six genes that, upon disruption of their function, suppress the embryonic lethal phenotype of par-1 mutant embryos. One of these genes is mpk-1, the C. elegans homolog of the highly conserved MAP kinase ERK. Closer analysis subsequently showed that reduction of function of mpk-1 not only increases viability of par-1 mutant embryos, but also reverts several polarity phenotypes associated with loss of function of par-1. Our data indicate that mpk-1 antagonizes par-1 activity to regulate polarization and asymmetric cell divisions in the early embryo.     

Construction of EST Database for Comparative Gene Studies of Acanthamoeba

The genus Acanthamoeba can cause severe infections such as granulomatous amebic encephalitis and amebic keratitis in humans. However, little genomic information of Acanthamoeba has been reported. Here, we constructed Acanthamoeba expressed sequence tags (EST) database (Acanthamoeba EST DB) derived from our 4 kinds of Acanthamoeba cDNA library. The Acanthamoeba EST DB contains 3,897 EST generated from amebae under various conditions of long term in vitro culture, mouse brain passage, or encystation, and downloaded data of Acanthamoeba from National Center for Biotechnology Information (NCBI) and Taxonomically Broad EST Database (TBestDB). The almost reported cDNA/genomic sequences of Acanthamoeba provide stand alone BLAST system with nucleotide (BLAST NT) and amino acid (BLAST AA) sequence database. In BLAST results, each gene links for the significant information including sequence data, gene orthology annotations, relevant references, and a BlastX result. This is the first attempt for construction of Acanthamoeba database with genes expressed in diverse conditions. These data were integrated into a database (http://www.amoeba.or.kr).     

PathVisio 3: An Extendable Pathway Analysis Toolbox

PathVisio is a commonly used pathway editor, visualization and analysis software. Biological pathways have been used by biologists for many years to describe the detailed steps in biological processes. Those powerful, visual representations help researchers to better understand, share and discuss knowledge. Since the first publication of PathVisio in 2008, the original paper was cited more than 170 times and PathVisio was used in many different biological studies. As an online editor PathVisio is also integrated in the community curated pathway database WikiPathways.Here we present the third version of PathVisio with the newest additions and improvements of the application. The core features of PathVisio are pathway drawing, advanced data visualization and pathway statistics. Additionally, PathVisio 3 introduces a new powerful extension systems that allows other developers to contribute additional functionality in form of plugins without changing the core application.PathVisio can be downloaded from http://www.pathvisio.org and in 2014 PathVisio 3 has been downloaded over 5,500 times. There are already more than 15 plugins available in the central plugin repository. PathVisio is a freely available, open-source tool published under the Apache 2.0 license (http://www.apache.org/licenses/LICENSE-2.0). It is implemented in Java and thus runs on all major operating systems. The code repository is available at http://svn.bigcat.unimaas.nl/pathvisio. The support mailing list for users is available on https://groups.google.com/forum/#!forum/wikipathways-discuss and for developers on https://groups.google.com/forum/#!forum/wikipathways-devel.

This is a PLOS Computational Biology software article.

     

PHProteomicDB： A Module for Two-dimensional Gel Electrophoresis Database Creation on Personal Web Sites

PHProteomicDB is a PHP-written module to help researchers in proteomics to share two-dimenslonal gel electrophoresis data using personal web sites. No technical or PHP knowledge is necessary except a few basics about web site management. PHProteomicDB has a user-friendly administration interface to enter and update data. It creates web pages on the fly displaying gel characteristics, gel pictures, and numbered gel spots with their related identifications pointing to their reference pages in protein databanks. The module is freely available at http：//www.huvec.com/index.php3？rub=Download.