Evolutionary Bioinformatics

Phylogenetic methods unearth evolutionary history when supported by three starting points of reason: (1) the axiom begs the existence of a "model" of evolutionary change, (2) the axiom defines the historical ground plan (phylogeny) in which biological entities (taxa) evolve, and (3) the axiom demands identification of biological attributes (characters) with historical information. Axiom consequences are interlinked, making the retrodiction enterprise an endeavor of reciprocal fulfillment. In particular, establishing direction of evolutionary change (character polarization) roots phylogenies and enables testing the existence of historical memory (homology). Unfortunately, rooting phylogenies, especially the "tree of life," generally follow narratives instead of integrating empirical and theoretical knowledge of retrodictive exploration. This stems mostly from a focus on molecular sequence analysis and uncertainties about rooting methods. Here, we review available ...

Molecular phylogenetics relies on accurate identification of orthologous sequences among the taxa of interest. Most orthology inference programs available for use in phylogenomics rely on small sets of pre-defined orthologs from model organisms or phenetic approaches such as all-versus-all sequence comparisons followed by Markov graph-based clustering. Such approaches have high sensitivity but may erroneously include paralogous sequences. We developed PhyloTreePruner, a software utility that uses a phylogenetic approach to refine orthology inferences made using phenetic methods. PhyloTreePruner checks single-gene trees for evidence of paralogy and generates a new alignment for each group containing only sequences inferred to be orthologs. Importantly, PhyloTreePruner takes into account support values on the tree and avoids unnecessarily deleting sequences in cases where a weakly supported tree topology incorrectly indicates paralogy. A test of PhyloTreePruner on a dataset generated ...

Foot-and-mouth disease (FMD), an endemic disease of cloven-hoofed animals, causes an annual economic loss of US$60-150 million in Bangladesh. There is no cross-protection among the foot-and-mouth disease virus (FMDV) serotypes and vaccination escape mutation may happen. Peptide vaccine is a safer alternative. The aim of this study is to predict and map the B and T cell epitopes of VP1 proteins of FMDV serotypes O and A that were circulating in Bangladesh from 2011 to 2013. Using evolutionary and computational approach (BCPred, BepiPred, DiscoTope, ElliPro, and ProPred-I, IEDB analysis for MHC-I prediction), a total of 11 B and T cell epitopes were predicted. Also, the three-dimensional (3D) structure of VP1 protein showed that the predicted five epitopes residing on N- and C-termini can be considered as good vaccine candidates, and epitopes on the G-H loop can serve as receptor recognition sites for vaccine design. The scores of predicted epitopes of one method were cross-checked wi...

Microbiomes are ubiquitous and are found in the ocean, the soil, and in/on other living organisms. Changes in the microbiome can impact the health of the environmental niche in which they reside. In order to learn more about these communities, different approaches based on data from multiple omics have been pursued. Metagenomics produces a taxonomical profile of the sample, metatranscriptomics helps us to obtain a functional profile, and metabolomics completes the picture by determining which byproducts are being released into the environment. Although each approach provides valuable information separately, we show that, when combined, they paint a more comprehensive picture. We conclude with a review of network-based approaches as applied to integrative studies, which we believe holds the key to in-depth understanding of microbiomes.

DNA sequences accumulating in the International Nucleotide Sequence Databases (INSD) form a rich source of information for taxonomic and ecological meta-analyses. However, these databases include many erroneous entries, and the data itself is poorly annotated with metadata, making it difficult to target and extract entries of interest with any degree of precision. Here we describe the web-based workbench PlutoF, which is designed to bridge the gap between the needs of contemporary research in biology and the existing ...

Our study searched all available sequences of from NCBI (National Center for Biotechnology Information) and tested for their species resolution capability in single as well as in combination forms. A total of 28 loci were applied for analyses in the study. From the nuclear genome, the highest resolution was of , followed by 4, and 51. These 4 loci were found to be even better than the popular region ITS for identification. Among the chloroplast regions, the intergenic spacer -L gave the highest species resolution (76.7%), followed by K, L, C2, and 1. The divergence of , 18S, 1, A, L, and 2 was very low and should not be used as identifying markers for . In addition, 2-locus combinations could improve significantly the resolving capability for the genus, in which 14/36 data sets could be resolved completely (100%) with interspecies relationships. The indel information was also effective supporting data for molecular discrimination of species.

Recent advances in virus detection strategies and deep sequencing technologies have enabled the identification of a multitude of new viruses that persistently infect mosquitoes but do not infect vertebrates. These are usually referred to as insect-specific viruses (ISVs). These novel viruses have generated considerable interest in their modes of transmission, persistence in mosquito populations, the mechanisms that restrict their host range to mosquitoes, and their interactions with pathogens transmissible by the same mosquito. In this article, we discuss studies in our laboratory and others that demonstrate that many ISVs are efficiently transmitted directly from the female mosquito to their progeny via infected eggs, and, moreover, that persistent infection of mosquito cell cultures or whole mosquitoes with ISVs can restrict subsequent infection, replication, and transmission of some mosquito-borne viral pathogens. This suggests that some ISVs may act as natural regulators of arbo...

Clostridium haemolyticum is the causal agent of bacillary hemoglobinuria in cattle, goat, sheep, and ruminants. In this study, we report the first recorded human-infecting C. haemolyticum strain collected from an 18-year-old woman diagnosed with acute lymphoblastic leukemia. After failure of traditional techniques, only next-generation sequencing (NGS) technology in combination with bioinformatics, phylogenetic, and pathogenomics analyses revealed that our King Faisal Specialist Hospital and Research Center (KFSHRC) bacterial isolate belongs to C. haemolyticum species. KFSHRC isolate is composed of 1 chromosome and 4 plasmids. The total genome size is estimated to be 2.7 Mbp with a low GC content of 28.02%. Comparative pathogenomics analysis showed that C. haemolyticum KFSHRC isolate is a potential virulent pathogenic bacterium as it possesses the virulence factors necessary to establish an infection, acquire essential nutrients, resist antimicrobial agents, and tolerate hostile con...

Human respiratory syncytial virus (RSV) is a major viral pathogen associated with acute lower respiratory tract infections (ALRTIs) among hospitalized children. In this study, the genetic diversity of the RSV strains was investigated among nasopharyngeal aspirates (NPA) taken from children less than 5 years of age hospitalized with ALRTIs in Hospital Serdang, Malaysia. A total of 165 NPA samples were tested for the presence of RSV and other respiratory viruses from June until December 2009. RSV was found positive in 83 (50%) of the samples using reverse transcription polymerase chain reaction (RT-PCR). Further classification of 67 RSV strains showed that subgroups A and B comprised 11/67 (16.4%) and 56/67 (83.6%) of the strains, respectively. The second hypervariable region at the carboxyl-terminal of the G gene was amplified and sequenced in order to do phylogenetic study. The phylogenetic relationships of the samples were determined separately for subgroups A and B using neighbor ...

Temperature is an abiotic factor that affects various biological and physiological processes in fish. Temperature stress is known to increase the production of reactive oxygen species (ROS) that subsequently cause oxidative stress. Fish is known to evolve a system of antioxidant enzymes to reduce ROS toxicology. Glutathione peroxidase (GPx) family consists of key enzymes that protect fish from oxidative stress. In this study, full-length GPx1 cDNA (GenBank accession no. KY984468) of Tor tambroides was cloned and characterized by rapid amplification of cDNA ends (RACE). The 899-base-pair (bp) GPx1 cDNA includes a 576-bp open reading frame encoding for 191 amino acids, plus 28 bp of 5′-untranslated region (UTR) and 295 bp of 3′-UTR. Homology analysis revealed that GPx1 of T tambroides (Tor-GPx1) shared high similarity with GPx1 sequences of other fish species. The phylogenetic construction based on the amino acid sequence showed that Tor-GPx1 formed a clade with GPx1 sequences of vari...

A comprehensive phylogeny of species, i.e., a tree of life, has potential uses in a variety of contexts, including research, education, and public policy. Yet, accessing the tree of life typically requires special knowledge, complex software, or long periods of training. The Phylotastic project aims make it as easy to get a phylogeny of species as it is to get driving directions from mapping software. In prior work, we presented a design for an open system to validate and manage taxon names, find phylogeny resources, extract subtrees matching a user’s taxon list, scale trees to time, and integrate related resources such as species images. Here, we report the implementation of a set of tools that together represent a robust, accessible system for on-the-fly delivery of phylogenetic knowledge. This set of tools includes a web portal to execute several customizable workflows to obtain species phylogenies (scaled by geologic time and decorated with thumbnail images); more than 30 underl...

Glucose transporters (GLUT) are twelve-transmembrane spanning proteins that contain two pores capable of transporting glucose and dehydroascorbate in and out of cells. The mechanism by which transport is effected is unknown. An evolutionarily-based hypothesis for the mechanism of glucose transport is presented here based on reports that insulin has multiple binding sites for glucose. It is proposed that insulin-like peptides were incorporated as modular elements into transmembrane proteins during evolution, resulting in glucose transporting capacity. Homology searching reveals that all GLUT contain multiple copies of insulin-like regions. These regions map onto a model of GLUT in positions that define the glucose transport cores. This observation provides a mechanism for glucose transport involving the diffusion of glucose from one insulin-like glucose-binding region to another. It also suggests a mechanism by which glucose disregulation may occur in both type 1 and type 2 diabetes:...

microRNAs (miRNAs) are a class of endogenous small RNAs (sRNAs) that play pivotal roles in plant development, abiotic stress response, and pathogen response. miRNAs have been extensively studied in plants, but rarely in Nicotiana benthamiana, despite its wide use in plant virology studies, particularly for studying N protein-tobacco mosaic virus (TMV) interactions. We report an efficient method using high-throughput sequencing and bioinformatics to identify genome-wide miRNAs in N. benthamiana. A total of 30 conserved miRNA families and 113 novel miRNAs belonging to 93 families were identified. Some miRNAs were clustered on chromosomes, and some were embedded in host gene introns. The predicted miRNA targets were involved in diverse biological processes, such as metabolism, signaling, and responses to stimuli. miRNA expression profiling revealed that most of them were differentially expressed during N-mediated immunity to TMV. This study provides a framework for further analysis of ...

In December, 2006, a group of 26 software developers from some of the most widely used life science programming toolkits and phylogenetic software projects converged on Durham, North Carolina, for a Phyloinformatics Hackathon, an intense five-day collaborative software coding event sponsored by the National Evolutionary Synthesis Center (NESCent). The goal was to help researchers to integrate multiple phylogenetic software tools into automated workflows. Participants addressed deficiencies in interoperability between programs by implementing “glue code” and improving support for phylogenetic data exchange standards (particularly NEXUS) across the toolkits. The work was guided by use-cases compiled in advance by both developers and users, and the code was documented as it was developed. The resulting software is freely available for both users and developers through incorporation into the distributions of several widely-used open-source toolkits. We explain the motivation for the hac...

Advances in high-throughput sequencing have generated a vast amount of transcriptomic data that are being increasingly used in phylogenetic reconstruction. However, processing the vast datasets for a huge number of genes and even identifying optimal analytical methodology are challenging. Through de novo sequenced and retrieved data from public databases, we identified 221 orthologous protein-coding genes to reconstruct the phylogeny of Ericales, an order characterized by rapid ancient radiation. Seven species representing different families in Ericales were used as in-groups. Both concatenation and coalescence methods yielded the same well-supported topology as previous studies, with only two nodes conflicting with previously reported relationships. The results revealed that a partitioning strategy could improve the traditional concatenation methodology. Rapidly evolving genes negatively affected the concatenation analysis, while slowly evolving genes slightly affected the coalesce...

The polyphasic approach is the most progressive system that has been suggested for distinguishing and phylogenetically classifying Cyanoprokaryota (Cyanobacteria/Cyanophyta). Several oscillatorialean genera ( Lyngbya, Phormidium, Plectonema, and Leptolyngbya) have problematic phylogenetic position and taxonomic state because of their heterogeneity and polyphyletic nature. To accurately resolve the phylogenetic relationship of some filamentous species ( Nodosilinea bijugata, Phormidium molle, Phormidium papyraceum), we have performed phylogenetic analyses based on 16S rRNA gene and the phycocyanin Operon (PC-IGS) by using maximum-likelihood (ML) tree inference methods. These analyses were combined with morphological re-evaluation. Our phylogenetic analyses support the taxonomic separation of genus Nodosilinea from the polyphyletic genus Leptolyngbya. Investigated Nodosilinea strains always formed a coherent genetic cluster supported with a high bootstrap value. The molecular phylogen...

The recent availability of the complete genome sequences of a large number of model organisms, together with the immense amount of data being produced by the new high-throughput technologies, means that we can now begin comparative analyses to understand the mechanisms involved in the evolution of the genome and their consequences in the study of biological systems. Phylogenetic approaches provide a unique conceptual framework for performing comparative analyses of all this data, for propagating information between different systems and for predicting or inferring new knowledge. As a result, phylogeny-based inference systems are now playing an increasingly important role in most areas of high throughput genomics, including studies of promoters (phylogenetic footprinting), interactomes (based on the presence and degree of conservation of interacting proteins), and in comparisons of transcriptomes or proteomes (phylogenetic proximity and co-regulation/co-expression). Here we review th...

Scientific workflows can be understood as arrangements of managed activities executed by different processing entities. It is a regular Bioinformatics approach applying workflows to solve problems in Molecular Biology, notably those related to sequence analyses. Due to the nature of the raw data and the in silico environment of Molecular Biology experiments, apart from the research subject, 2 practical and closely related problems have been studied: reproducibility and computational environment. When aiming to enhance the reproducibility of Bioinformatics experiments, various aspects should be considered. The reproducibility requirements comprise the data provenance, which enables the acquisition of knowledge about the trajectory of data over a defined workflow, the settings of the programs, and the entire computational environment. Cloud computing is a booming alternative that can provide this computational environment, hiding technical details, and delivering a more affordable, accessible, and configurable on-demand environment for researchers. Considering this specific scenario, we proposed a solution to improve the reproducibility of Bioinformatics workflows in a cloud computing environment using both Infrastructure as a Service (IaaS) and Not only SQL (NoSQL) database systems. To meet the goal, we have built 3 typical Bioinformatics workflows and ran them on 1 private and 2 public clouds, using different types of NoSQL database systems to persist the provenance data according to the Provenance Data Model (PROV-DM). We present here the results and a guide for the deployment of a cloud environment for Bioinformatics exploring the characteristics of various NoSQL database systems to persist provenance data.