Computational Biology Commons^™

Testing The Independence Hypothesis Of Accepted Mutations For Pairs Of Adjacent Amino Acids In Protein Sequences, Jyotsna Ramanan, Peter Revesz

School of Computing: Faculty Publications

Evolutionary studies usually assume that the genetic mutations are independent of each other. However, that does not imply that the observed mutations are independent of each other because it is possible that when a nucleotide is mutated, then it may be biologically beneficial if an adjacent nucleotide mutates too. With a number of decoded genes currently available in various genome libraries and online databases, it is now possible to have a large-scale computer-based study to test whether the independence assumption holds for pairs of adjacent amino acids. Hence the independence question also arises for pairs of adjacent amino acids within …

Isolation And Comparative Genomic Analysis Of Final Third Of Satis Genome, Kelly Hartigan, Nicole Curnutt, Matthew Mcdermut

Undergraduate Research Symposium Posters

A highly novel Streptomyces phage, Satis, was isolated from a direct environmental sample collected from outside Danforth House on the Washington University campus. Satis infects bacterial species Streptomyces lividans producing pinpoint, cloudy plaques less than 1mm in diameter. Electron microscope data shows rare atypical physical features. Rather than the common octahedral capsid shape, Satis has a prolate head with visible cross-linked hexagonal protein structure and average measurements of 285 nm by 47 nm with a long, flexible tail measuring 268 nm. Upon sequencing, it was found that Satis contains the longest phage genome discovered to date through the SEA-PHAGE program …

Statistical Methods For Two Problems In Cancer Research: Analysis Of Rna-Seq Data From Archival Samples And Characterization Of Onset Of Multiple Primary Cancers, Jialu Li

Dissertations and Theses (Open Access)

My dissertation is focused on quantitative methodology development and application for two important topics in translational and clinical cancer research.

The first topic was motivated by the challenge of applying transcriptome sequencing (RNA-seq) to formalin-fixation and paraffin-embedding (FFPE) tumor samples for reliable diagnostic development. We designed a biospecimen study to directly compare gene expression results from different protocols to prepare libraries for RNA-seq from human breast cancer tissues, with randomization to fresh-frozen (FF) or FFPE conditions. To comprehensively evaluate the FFPE RNA-seq data quality for expression profiling, we developed multiple computational methods for assessment, such as the uniformity and continuity …

Phylogenetic Studies Of The Madagascan Freshwater Crabs (Potamoidea, Potamonautidae, Deckeniinae), Rainee Stevens

All NMU Master's Theses

Relationships within the Madagascan freshwater crab fauna were examined based on unstudied specimens from all parts of Madagascar. This study allowed the examination of the validity of existing genera and species, and identified potential new taxa. The new specimens provided phylogenetic and evolutionary data, as well as new insights into the distribution patterns of the Malagasy freshwater crab fauna. In addition, the large number of new localities from previously unsurveyed areas of Madagascar allowed the construction of updated distribution maps.

The present analysis included 62 unidentified specimens plus 13 identified species of Malagasy freshwater crab taxa that had already been …

Discovery And Validation Of Information Theory-Based Transcription Factor And Cofactor Binding Site Motifs., Ruipeng Lu, Eliseos J Mucaki, Peter K Rogan

Biochemistry Publications

Data from ChIP-seq experiments can derive the genome-wide binding specificities of transcription factors (TFs) and other regulatory proteins. We analyzed 765 ENCODE ChIP-seq peak datasets of 207 human TFs with a novel motif discovery pipeline based on recursive, thresholded entropy minimization. This approach, while obviating the need to compensate for skewed nucleotide composition, distinguishes true binding motifs from noise, quantifies the strengths of individual binding sites based on computed affinity and detects adjacent cofactor binding sites that coordinate with the targets of primary, immunoprecipitated TFs. We obtained contiguous and bipartite information theory-based position weight matrices (iPWMs) for 93 sequence-specific TFs, …

Transcriptome Approach For Identifying Potential Biomarkers For Endocrine Disruption Due To Crude Oil Exposure Using Killifish (Fundulus Heteroclitus)., Frank J. Zadlock Iv

Seton Hall University Dissertations and Theses (ETDs)

The killifish, Fundulus heteroclitus, is a common fish model in aquatic toxicology. However, little is known in this organism about how endocrine disrupting compounds, (EDCs) including polycyclic aromatic hydrocarbons (PAHs), impact the reproductive status from a molecular standpoint. The objective of this project is to apply high throughput sequencing to F. heteroclitus testes to examine the molecular mechanisms that impair reproduction when exposed to crude oil. First, a crude oil exposure experiment was conducted. Following exposure, semi-quantitative PCR was performed to detect changes in gene expression of the following genes known to respond to EDCs: gonadal aromatase (CYP19a), vitellogenin …

P08. Unravelling Organelle Genome Evolution Architecture Using Rna-Sequencing Data, Matheus Sanita Lima, David Roy Smith

Western Research Forum

Background: Mitochondria genomes vary from 11 Mb to 6 kb, while plastids can vary from 1 Mb to 30 kb. Non-coding DNA accounts for most of this size variation, but the mechanistic and evolutionary reasons for that are still unknown. Next generation sequencing has generated unprecedented amounts of genomic and transcriptomic data that can be used for organelle genome evolution studies. However, most of these data is used only for the study of cell nucleus. Therefore, I decided to use these untapped data source to investigate the transcription of organelle genomes in plastid-bearing protists.

Methods: I mapped the transcriptomes over …

An Out-Of-Core Gpu Based Dimensionality Reduction Algorithm For Big Mass Spectrometry Data And Its Application In Bottom-Up Proteomics, Muaaz Awan, Fahad Saeed

Parallel Computing and Data Science Lab Technical Reports

Modern high resolution Mass Spectrometry instruments can generate millions of spectra in a single systems biology experiment. Each spectrum consists of thousands of peaks but only a small number of peaks actively contribute to deduction of peptides. Therefore, pre-processing of MS data to detect noisy and non-useful peaks are an active area of research. Most of the sequential noise reducing algorithms are impractical to use as a pre-processing step due to high time-complexity. In this paper, we present a GPU based dimensionality-reduction algorithm, called G-MSR, for MS2 spectra. Our proposed algorithm uses novel data structures which optimize the memory and …

Characterization Of A Large Vertebrate Genome And Homomorphic Sex Chromosomes In The Axolotl, Ambystoma Mexicanum, Melissa Keinath

Theses and Dissertations--Biology

Changes in the structure, content and morphology of chromosomes accumulate over evolutionary time and contribute to cell, developmental and organismal biology. The axolotl (Ambystoma mexicanum) is an important model for studying these changes because: 1) it provides important phylogenetic perspective for reconstructing the evolution of vertebrate genomes and amphibian karyotypes, 2) its genome has evolved to a large size (~10X larger than human) but has maintained gene orders, and 3) it possesses potentially young sex chromosomes that have not undergone extensive differentiation in the structure that is typical of many other vertebrate sex chromosomes (e.g. mammalian XY chromosomes …

Identification Of Novel Sleep Related Genes From Large Scale Phenotyping Experiments In Mice, Shreyas Joshi

Theses and Dissertations--Biology

Humans spend a third of their lives sleeping but very little is known about the physiological and genetic mechanisms controlling sleep. Increased data from sleep phenotyping studies in mouse and other species, genetic crosses, and gene expression databases can all help improve our understanding of the process. Here, we present analysis of our own sleep data from the large-scale phenotyping program at The Jackson Laboratory (JAX), to identify the best gene candidates and phenotype predictors for influencing sleep traits.

The original knockout mouse project (KOMP) was a worldwide collaborative effort to produce embryonic stem (ES) cell lines with one of …

Genomic Perspectives On Amphibian Evolution Across Multiple Phylogenetic Scales, Paul Michael Hime

Theses and Dissertations--Biology

Genomes provide windows into the evolutionary histories of species. The recent accessibility of genome-scale data in non-model organisms and the proliferation of powerful statistical models are now providing unprecedented opportunities to uncover evolutionary relationships and to test hypotheses about the processes that generate and maintain biodiversity. This dissertation work reveals shallow-scale species boundaries and population genetic structure in two imperiled groups of salamanders and demonstrates that the number and information content of genomic regions used in species delimitation exert strong effects on the resulting inferences. Genome scans are employed to test hypotheses about the mechanisms of genetic sex determination in …

Red Snapper Distribution On Natural Habitats And Artificial Structures In The Northern Gulf Of Mexico, Mandy Karnauskas, John F. Walter Iii, Matthew D. Campbell, Adam G. Pollack, J. Marcus Drymon, Sean P. Powers

University Faculty and Staff Publications

In 2011, an intensive, multiple-gear, fishery-independent survey was carried out in the northern Gulf of Mexico (GOM) to collect comprehensive age and length information on Red Snapper Lutjanus campechanus. Based on this synoptic survey, we produced a spatial map of Red Snapper relative abundance that integrates both gear selectivity effects and ontogenetically varying habitat usage. Our methodology generated a spatial map of Red Snapper at a 10-km2 grid resolution that is consistent with existing knowledge of the species: Red Snapper occurred in relatively high abundances at depths of 50–90 m along the coasts of Texas and Louisiana and in smaller, …

K-Mer Analysis Pipeline For Classification Of Dna Sequences From Metagenomic Samples, Russell Kaehler

Graduate Student Theses, Dissertations, & Professional Papers

Biological sequence datasets are increasing at a prodigious rate. The volume of data in these datasets surpasses what is observed in many other fields of science. New developments wherein metagenomic DNA from complex bacterial communities is recovered and sequenced are producing a new kind of data known as metagenomic data, which is comprised of DNA fragments from many genomes. Developing a utility to analyze such metagenomic data and predict the sample class from which it originated has many possible implications for ecological and medical applications. Within this document is a description of a series of analytical techniques used to process …

Comparative Population Genomics And Speciation Of Snakes Across The North American Deserts, Edward A. Myers

Dissertations, Theses, and Capstone Projects

Understanding the process of speciation is of central interest to evolutionary biologists. Speciation can be studied using a phylogeographic approach, by identifying regions that promote lineage divergence, addressing whether speciation has occurred with gene flow, and when extended to multiple taxa, addressing if the same patterns of speciation are shared across codistributed groups with different ecologies. Here I examine the comparative phylogeographic histories and population genomics of thirteen snake taxa that are widely distributed and co-occur across the arid southwest of North America. I first quantify the degree to which these species groups have a shared history of population divergence …

Modeling And Analysis Of Germ Layer Formations Using Finite Dynamical Systems, Alexander Garza, Megan Eberle, Eric A. Eager

Spora: A Journal of Biomathematics

The development of an embryo from a fertilised egg to a multicellular organism proceeds through numerous steps, with the formation of the three germ layers (endoderm, mesoderm, ectoderm) being one of the first. In this paper we study the mesendoderm (the tissue that collectively gives rise to both mesoderm and endoderm) gene regulatory network for two species, \textit{Xenopus laevis} and the axolotl (\textit{Ambystoma mexicanum}) using Boolean networks. We find that previously-established bistability found in these networks can be reproduced using this Boolean framework, provided that some assumptions used in previously-published differential equations models are relaxed. We conclude by discussing our …

Data Development And Analysis Pathways For Marine Mammals And Turtles: Creating A User Interface, Sarina Fernandez, Warren Asfazadour, Eric Archer, Lisa Komoroske

STAR Program Research Presentations

A major obstacle in genetic research is developing streamlined methods for analyzing large amounts of data. The statistical computer programming language R provides users with the ability to develop packages containing specific functions in order to create more accessible data analysis pipelines. However, writing code in R can still be intimidating to those with little to no coding experience. Fortunately, the R package shiny provides a framework for developing web applications based on R functions. Using shiny, we developed a user-friendly web application containing functions of the R package strataG. The strataG package contains several functions for summarizing genetic data …

Development Of An In Silico Kir Genotyping Algorithm And Its Application To Population And Cancer Immunogenetic Analyses, Howard Rosoff

Dissertations and Theses (Open Access)

Gene content determination and variant calling in the complex KIR genomic region are useful for immune system function analysis, pathogenesis and disease risk factor elucidation, immunotherapy development, evolutionary investigations, and human migration modeling. Sequence-specific oligonucleotide and sequence-specific primer PCR methods are the de facto standards for KIR presence/absence identification, but the current platforms are unsuitable for SNP calling, impractical for KIR typing large cohorts of DNA samples, and inapplicable for typing repositories in which sequence data, but not cells or cell analytes, are available. Alternative typing methods, such as in silico sequence-based typing, can address the problems associated with amplicon-based …

Incremental Phylogenetics By Repeated Insertions: An Evolutionary Tree Algorithm, Peter Revesz, Zhiqiang Li

School of Computing: Faculty Publications

We introduce the idea of constructing hypothetical evolutionary trees using an incremental algorithm that inserts species one-by-one into the current evolutionary tree. The method of incremental phylogenetics by repeated insertions lead to an algorithm that can be used on DNA, RNA and amino acid sequences. According to experimental results on both synthetic and biological data, the new algorithm generates more accurate evolutionary trees than the UPGMA and the Neighbor Joining algorithms.

Ten Simple Rules For Taking Advantage Of Git And Github, Yasset Perez-Riverol, Laurent Gatto, Rui Wang, Timo Sachsenberg, Julian Uszkoreit, Felipe Da Veiga Leprevost, Christian Fufezan, Tobias Ternent, Stephen J. Eglen, Daniel S. Katz, Tom J. Pollard, Alexander Konovalov, Robert M. Flight, Kai Blin, Juan Antonio Vizcaíno

Molecular and Cellular Biochemistry Faculty Publications

No abstract provided.

Identification, Characterization, And Life Cycle Of Intein-Associated Homing Endonucleases, Joshua J. Skydel

Honors Scholar Theses

Inteins are molecular parasites that have been identified in unicellular organisms from the three domains of life. The intein self-excises following translation of the host gene, and therefore incurs a fitness cost for its carrier. The symbiotic state of the intein to its host is dependent on the presence or absence of a homing endonuclease domain, which facilitates horizontal transfer of the molecule. Identification of this domain provides information on the evolutionary history of the intein, as well as patterns of horizontal gene transfer in microbial communities. I have therefore developed Hidden Markov Models (HMMs) to identify homing endonuclease domains …

Using Hadoop To Identify False Positives In Bacterial Strain Typing From Dna Fingerprints, Colin C. Adams

Computer Science and Software Engineering

Pyroprinting is a novel technique used by the Department of Biological Sciences to obtain “fingerprints” from the DNA of E. coli isolates in order to categorize them into strains. To determine the number of false positives that occur in the pyroprinting process, isolates with the same pyroprints needed to be sequenced to see if their underlying alleles match. If they do match, this shows they are indeed the same strain and are a true positive. If the alleles don’t match, they are different strains and are a false positive. To do this 100 isolates with nucleotide identifiers were sequenced. Over …

Phylogenetic Analysis Of Human Cytomegalovirus Pus27 And Pus28: Ascertaining An Independent Or Linked Evolutionary History, Jessica A. Scarborough

Undergraduate Honors Theses

Human cytomegalovirus (HCMV) is a widespread pathogen that is particularly skilled at evading immune detection and defense mechanisms, largely due to extensive co-evolution with its host’s immune system. One aspect of this co-evolution involves the acquisition of four virally encoded GPCR chemokine receptor homologs, products of the US27, US28, UL33 and UL78 genes. G protein-coupled receptors (GPCR) are the largest family of cell surface proteins, found in organisms from yeast to humans. In this research, phylogenetic analysis was used to investigate the origins of the US27 and US28 genes, which are adjacent in the viral genome. The results indicate that …

Investigating Metastatic Lineage In Colorectal Cancer By Single Cell Dna Sequencing, Marco Leung

Dissertations and Theses (Open Access)

Metastasis is the primary cause of human cancer deaths. Patients with metastatic colorectal cancer (mCRC) show only an 11% 5-year survival rate, compared to those without local or distant metastases (92% 5-year survival rate). Understanding the CRC tumor evolution may provide valuable insights on how to improve treatment in patients with mCRC. However, the genomic basis of metastasis has been difficult to study, in part due to the extensive intratumor heterogeneity at both the primary and metastatic tumor sites, and the low frequency of subclones with metastatic potential. Previous studies have applied conventional bulk next-generation sequencing (NGS) methods, which have …

Accurate Mutation Annotation And Functional Prediction Enhance The Applicability Of -Omics Data In Precision Medicine, Tenghui Chen

Dissertations and Theses (Open Access)

Clinical sequencing has been recognized as an effective approach for enhancing the accuracy and efficiency of cancer patient management and therefore achieve the goals of personalized therapy. However, the accuracy of large scale sequencing data in clinics has been constrained by many different aspects, such as clinical detection, annotation and interpretation of the variants that are observed in clinical sequencing data. In my Ph.D thesis work, I mainly investigated how to comprehensively and efficiently apply high dimensional -omics data to enhance the capability of precision cancer medicine. Following this motivation, my dissertation has been focused on two important topics in …

Detecting Gene-Gene Interactions Using A Permutation-Based Random Forest Method, Jing Li, James D. Malley, Angeline S. Andrew, Margaret R. Karagas, Jason H. Moore

Dartmouth Scholarship

Identifying gene-gene interactions is essential to understand disease susceptibility and to detect genetic architectures underlying complex diseases. Here, we aimed at developing a permutation-based methodology relying on a machine learning method, random forest (RF), to detect gene-gene interactions. Our approach called permuted random forest (pRF) which identified the top interacting single nucleotide polymorphism (SNP) pairs by estimating how much the power of a random forest classification model is influenced by removing pairwise interactions.

Fastpop: A Rapid Principal Component Derived Method To Infer Intercontinental Ancestry Using Genetic Data, Yafang Li, Jinyoung Byun, Guoshuai Cai, Xiangjun Xiao, Younghun Han, Olivier Cornelis, James E. Dinulos, Joe Dennis, Douglas Easton, Ivan Gorlov, Michael F. Seldin, Christopher I. Amos

Dartmouth Scholarship

Identifying subpopulations within a study and inferring intercontinental ancestry of the samples are important steps in genome wide association studies. Two software packages are widely used in analysis of substructure: Structure and Eigenstrat. Structure assigns each individual to a population by using a Bayesian method with multiple tuning parameters. It requires considerable computational time when dealing with thousands of samples and lacks the ability to create scores that could be used as covariates. Eigenstrat uses a principal component analysis method to model all sources of sampling variation. However, it does not readily provide information directly relevant to ancestral origin; the …

Phagephisher: A Pipeline For The Discovery Of Covert Viral Sequences In Complex Genomic Datasets, Thomas Hatzopoulos, Siobhan C. Watkins, Catherine Putonti

Bioinformatics Faculty Publications

Obtaining meaningful viral information from large sequencing datasets presents unique challenges distinct from prokaryotic and eukaryotic sequencing efforts. The difficulties surrounding this issue can be ascribed in part to the genomic plasticity of viruses themselves as well as the scarcity of existing information in genomic databases. The open-source software PhagePhisher (http://www.putonti-lab.com/phagephisher) has been designed as a simple pipeline to extract relevant information from complex and mixed datasets, and will improve the examination of bacteriophages, viruses, and virally related sequences, in a range of environments. Key aspects of the software include speed and ease of use; PhagePhisher can be used with …

Hpcnmf: A High-Performance Toolbox For Non-Negative Matrix Factorization, Karthik Devarajan, Guoli Wang

COBRA Preprint Series

Non-negative matrix factorization (NMF) is a widely used machine learning algorithm for dimension reduction of large-scale data. It has found successful applications in a variety of fields such as computational biology, neuroscience, natural language processing, information retrieval, image processing and speech recognition. In bioinformatics, for example, it has been used to extract patterns and profiles from genomic and text-mining data as well as in protein sequence and structure analysis. While the scientific performance of NMF is very promising in dealing with high dimensional data sets and complex data structures, its computational cost is high and sometimes could be critical for …

Models For Hsv Shedding Must Account For Two Levels Of Overdispersion, Amalia Magaret

UW Biostatistics Working Paper Series

We have frequently implemented crossover studies to evaluate new therapeutic interventions for genital herpes simplex virus infection. The outcome measured to assess the efficacy of interventions on herpes disease severity is the viral shedding rate, defined as the frequency of detection of HSV on the genital skin and mucosa. We performed a simulation study to ascertain whether our standard model, which we have used previously, was appropriately considering all the necessary features of the shedding data to provide correct inference. We simulated shedding data under our standard, validated assumptions and assessed the ability of 5 different models to reproduce the …

Resolving Gnetum Evolutionary History, Angela Mcfadden

All Master's Theses

Gnetum are non-flowering seed plants of the tropics, indigenous to South America, Africa, and Asia. This group of about 40 species is fascinating to botanists because it shares distinctive morphological characteristics with flowering plants, such as broad leaves, woody stems, and flower-like strobili. There are still questions surrounding the relationships within the genus of Gnetum. With that in mind, I focused my work on generating phylogenetic hypotheses, using two molecular data sets: a concatenation of over 60 different chloroplast genes (66,815 base pairs), and the whole chloroplast genome (128,772 base pairs). This allowed me to compare the two phylogenies …