Chapter 11: Genome-Wide Association Studies, William S. Bush, Jason H. Moore

Dartmouth Scholarship

Genome-wide association studies (GWAS) have evolved over the last ten years into a powerful tool for investigating the genetic architecture of human disease. In this work, we review the key concepts underlying GWAS, including the architecture of common diseases, the structure of common human genetic variation, technologies for capturing genetic information, study designs, and the statistical methods used for data analysis. We also look forward to the future beyond GWAS.

Phylogeny And Evolutionary Patterns In The Dwarf Crayfish Subfamily (Decapoda: Cambarellinae), C. Pedraza-Lara, I. Doadrio, J. Breinholt, Keith A. Crandall

Computational Biology Institute

The Dwarf crayfish or Cambarellinae, is a morphologically singular subfamily of decapod crustaceans that contains only one genus, Cambarellus. Its intriguing distribution, along the river basins of the Gulf Coast of United States (Gulf Group) and into Central México (Mexican Group), has until now lacked of satisfactory explanation. This study provides a comprehensive sampling of most of the extant species of Cambarellus and sheds light on its evolutionary history, systematics and biogeography. We tested the impact of Gulf Group versus Mexican Group geography on rates of cladogenesis using a maximum likelihood framework, testing different models of birth/extinction of lineages. We …

Cross-Ontology Multi-Level Association Rule Mining In The Gene Ontology., Prashanti Manda, Seval Ozkan, Hui Wang, Fiona M. Mccarthy, Susan M. Bridges

BCoE Publications

The Gene Ontology (GO) has become the internationally accepted standard for representing function, process, and location aspects of gene products. The wealth of GO annotation data provides a valuable source of implicit knowledge of relationships among these aspects. We describe a new method for association rule mining to discover implicit co-occurrence relationships across the GO sub-ontologies at multiple levels of abstraction. Prior work on association rule mining in the GO has concentrated on mining knowledge at a single level of abstraction and/or between terms from the same sub-ontology. We have developed a bottom-up generalization procedure called Cross-Ontology Data Mining-Level by …

Phylogenetics Links Monster Larva To Deep-Sea Shrimp., Heather D Bracken-Grissom, Darryl L Felder, Nicole L Vollmer, Joel W Martin, Keith A Crandall

Computational Biology Institute

Mid-water plankton collections commonly include bizarre and mysterious developmental stages that differ conspicuously from their adult counterparts in morphology and habitat. Unaware of the existence of planktonic larval stages, early zoologists often misidentified these unique morphologies as independent adult lineages. Many such mistakes have since been corrected by collecting larvae, raising them in the lab, and identifying the adult forms. However, challenges arise when the larva is remarkably rare in nature and relatively inaccessible due to its changing habitats over the course of ontogeny. The mid-water marine species Cerataspis monstrosa (Gray 1828) is an armored crustacean larva whose adult identity …

Physiologically-Based Pharmacokinetic Modeling Of Acetaminophen Metabolism And Toxicity, David M. Ng, Ali Navid

STAR Program Research Presentations

Acetaminophen is a common analgesic and antipyretic. Metabolism of acetaminophen and acetaminophen-induced liver necrosis are predicted using physiologically-based pharmacokinetic (PBPK) modeling. Pharmacokinetic means the model determines where the drug is distributed in the body over time. Physiologically-based means the anatomy and physiology of the human body is reflected in the structure and functioning of the model. Acetaminophen is usually safe and effective when taken as recommended, but consumption at higher levels may lead to liver damage. Additionally, other factors such as alcoholic liver disease, smoking, and malnutrition affect the maximum safe dose of acetaminophen.

Classification Of Genomic Sequences By Latent Semantic Analysis, Samuel F. Way

Department of Electrical and Computer Engineering: Dissertations, Theses, and Student Research

Evolutionary distance measures provide a means of identifying and organizing related organisms by comparing their genomic sequences. As such, techniques that quantify the level of similarity between DNA sequences are essential in our efforts to decipher the genetic code in which they are written.

Traditional methods for estimating the evolutionary distance separating two genomic sequences often require that the sequences first be aligned before they are compared. Unfortunately, this preliminary step imposes great computational burden, making this class of techniques impractical for applications involving a large number of sequences. Instead, we desire new methods for differentiating genomic sequences that eliminate …

Systems Biology Approaches To Probe Gene Regulation In Bacteria, Diogo F. Troggian Veiga

Dissertations and Theses (Open Access)

Mechanisms that allow pathogens to colonize the host are not the product of isolated genes, but instead emerge from the concerted operation of regulatory networks. Therefore, identifying components and the systemic behavior of networks is necessary to a better understanding of gene regulation and pathogenesis. To this end, I have developed systems biology approaches to study transcriptional and post-transcriptional gene regulation in bacteria, with an emphasis in the human pathogen Mycobacterium tuberculosis (Mtb).

First, I developed a network response method to identify parts of the Mtb global transcriptional regulatory network utilized by the pathogen to counteract phagosomal stresses …

Gene Ontology Analysis Of Pairwise Genetic Associations In Two Genome-Wide Studies Of Sporadic Als, Nora Chung Kim, Peter C. Andrews, Folkert W. Asselbergs, H Robert Frost, Scott M. Williams, Brent T. Harris, Cynthia Read, Kathleen D. Askland, Jason H. Moore

Dartmouth Scholarship

It is increasingly clear that common human diseases have a complex genetic architecture characterized by both additive and nonadditive genetic effects. The goal of the present study was to determine whether patterns of both additive and nonadditive genetic associations aggregate in specific functional groups as defined by the Gene Ontology (GO).

Dramatic Shifts In Benthic Microbial Eukaryote Communities Following The Deepwater Horizon Oil Spill, Holly M. Bik, Kenneth M. Halanych, Jyotsna Sharma, W. Kelley Thomas

Hubbard Center for Genome Studies (HCGS)

Benthic habitats harbour a significant (yet unexplored) diversity of microscopic eukaryote taxa, including metazoan phyla, protists, algae and fungi. These groups are thought to underpin ecosystem functioning across diverse marine environments. Coastal marine habitats in the Gulf of Mexico experienced visible, heavy impacts following the Deepwater Horizon oil spill in 2010, yet our scant knowledge of prior eukaryotic biodiversity has precluded a thorough assessment of this disturbance. Using a marker gene and morphological approach, we present an intensive evaluation of microbial eukaryote communities prior to and following oiling around heavily impacted shorelines. Our results show significant changes in community structure, …

An Integrated Bioinformatics And Computational Biology Approach Identifies New Bh3-Only Protein Candidates, Robert G. Hawley, Yuzhong Chen, Irene Riz, Chen Zeng

Anatomy and Regenerative Biology Faculty Publications

FoxD4L1 is a forkhead transcription factor that expands the neural ectoderm by down-regulating genes that promote the onset of neural differentiation and up-regulating genes that maintain proliferative neural precursors in an immature state. We previously demonstrated that binding of Grg4 to an Eh-1 motif enhances the ability of FoxD4L1 to down-regulate target neural genes but does not account for all of its repressive activity. Herein we analyzed the protein sequence for additional interaction motifs and secondary structure. Eight conserved motifs were identified in the C-terminal region of fish and frog proteins. Extending the analysis to mammals identified a high scoring …

Sparse Integrative Clustering Of Multiple Omics Data Sets, Ronglai Shen, Sijian Wang, Qianxing Mo

Memorial Sloan-Kettering Cancer Center, Dept. of Epidemiology & Biostatistics Working Paper Series

High resolution microarrays and second-generation sequencing platforms are powerful tools to investigate genome-wide alterations in DNA copy number, methylation, and gene expression associated with a disease. An integrated genomic profiling approach measuring multiple omics data types simultaneously in the same set of biological samples would render an integrated data resolution that would not be available with any single data type. In a previous publication (Shen et al., 2009), we proposed a latent variable regression with a lasso constraint (Tibshirani, 1996) for joint modeling of multiple omics data types to identify common latent variables that can be used to cluster patient …

Structural Characterization Of Human Uch37, E. S. Burgie, C. Bingman, Ameet Soni, G. N. Phillips Jr.

Computer Science Faculty Works

Uch37 is a de-ubiquitylating enzyme that is functionally linked with the 26S proteasome via Rpn13, and is essential for metazoan development. Here, we report the X-ray crystal structure of full-length human Uch37 at 2.95 Å resolution. Uch37's catalytic domain is similar to those of all UCH enzymes characterized to date. The C-terminal extension is elongated, predominantly helical and contains coiled coil interactions. Additionally, we provide an initial characterization of Uch37's oligomeric state and identify a systematic error in previous analyses of Uch37 activity. Taken together, these data provide a strong foundation for further analysis of Uch37's several functions.

Probabilistic Ensembles For Improved Inference In Protein-Structure Determination, Ameet Soni, J. Shavlik

Computer Science Faculty Works

Protein X-ray crystallography — the most popular method for determining protein structures — remains a laborious process requiring a great deal of manual crystallographer effort to interpret low-quality protein images. Automating this process is critical in creating a high-throughput protein-structure determination pipeline. Previously, our group developed ACMI, a probabilistic framework for producing protein-structure models from electron-density maps produced via X-ray crystallography. ACMI uses a Markov Random Field to model the three-dimensional (3D) location of each non-hydrogen atom in a protein. Calculating the best structure in this model is intractable, so ACMI uses approximate inference methods to estimate the optimal structure. …

Desert Hedgehog Is A Mammal-Specific Gene Expressed During Testicular And Ovarian Development In A Marsupial, William A. O'Hara

Master's Theses

Desert hedgehog (DHH) belongs to the hedgehog gene family that act as secreted intercellular signal transducers. DHH is an essential morphogen for normal testicular development and function in both mice and humans but is not present in the avian lineage. Like other hedgehog proteins, DHH signals through the patched (PTCH) receptors 1 and 2. Here we examine the expression and protein distribution of DHH, PTCH1 and PTCH2 in the developing testes of a marsupial mammal (the tammar wallaby) to determine whether DHH signalling is a conserved factor in gonadal development in all therian mammals.

Phylogenetic Search Through Partial Tree Mixing., Kenneth Sundberg, Mark Clement, Quinn Snell, Dan Ventura, Michael Whiting, Keith Crandall

Computational Biology Institute

BACKGROUND: Recent advances in sequencing technology have created large data sets upon which phylogenetic inference can be performed. Current research is limited by the prohibitive time necessary to perform tree search on a reasonable number of individuals. This research develops new phylogenetic algorithms that can operate on tens of thousands of species in a reasonable amount of time through several innovative search techniques.

RESULTS: When compared to popular phylogenetic search algorithms, better trees are found much more quickly for large data sets. These algorithms are incorporated in the PSODA application available at http://dna.cs.byu.edu/psoda

CONCLUSIONS: The use of Partial Tree Mixing …

Testing Phylogenetic Hypotheses Of The Subgenera Of The Freshwater Crayfish Genus Cambarus (Decapoda: Cambaridae)., Jesse W Breinholt, Megan L Porter, Keith A Crandall

Computational Biology Institute

BACKGROUND: The genus Cambarus is one of three most species rich crayfish genera in the Northern Hemisphere. The genus has its center of diversity in the Southern Appalachians of the United States and has been divided into 12 subgenera. Using Cambarus we test the correspondence of subgeneric designations based on morphology used in traditional crayfish taxonomy to the underlying evolutionary history for these crayfish. We further test for significant correlation and explanatory power of geographic distance, taxonomic model, and a habitat model to estimated phylogenetic distance with multiple variable regression.

METHODOLOGY/PRINCIPAL FINDINGS: We use three mitochondrial and one nuclear gene …

A Genomic Island In Salmonella Enterica Ssp. Salamae Provides New Insights On The Genealogy Of The Locus Of Enterocyte Effacement., P Scott Chandry, Simon Gladman, Sean C Moore, Torsten Seemann, Keith A Crandall, Narelle Fegan

Computational Biology Institute

The genomic island encoding the locus of enterocyte effacement (LEE) is an important virulence factor of the human pathogenic Escherichia coli. LEE typically encodes a type III secretion system (T3SS) and secreted effectors capable of forming attaching and effacing lesions. Although prominent in the pathogenic E. coli such as serotype O157:H7, LEE has also been detected in Citrobacter rodentium, E. albertii, and although not confirmed, it is likely to also be in Shigella boydii. Previous phylogenetic analysis of LEE indicated the genomic island was evolving through stepwise acquisition of various components. This study describes a new LEE region from two …

Mtbindingsim: Simulate Protein Binding To Microtubules, Julia T. Philip, Charles H. Pence, Holly V. Goodson

Faculty Publications

Summary: Many protein–protein interactions are more complex than can be accounted for by 1:1 binding models. However, biochemists have few tools available to help them recognize and predict the behaviors of these more complicated systems, making it difficult to design experiments that distinguish between possible binding models. MTBindingSim provides researchers with an environment in which they can rapidly compare different models of binding for a given scenario. It is written specifically with microtubule polymers in mind, but many of its models apply equally well to any polymer or any protein–protein interaction. MTBindingSim can thus both help in training intuition about …

In Silico Prediction Of Protein-Protein Interactions Between Cronobacter Sakazakii And Human Host, Sean Held

Theses

The genus Cronobacter accommodates 16 biogroups of the emerging opportunistic pathogen known formerly as Enterobacter sakazakii. Cronobacter are occasional contaminants of milk powder and, consequently, powdered infant formula (PIT). Since it is a pathogen that is transmitted in infant formula it represents a significant health risk to neonates. These bacteria are opportunistic pathogens (can only cause disease when the host's resistance is low) and are linked with life-threatening infections in neonates. Some of the clinical symptoms of Cronohacter infection include necrotizing enterocolitis, bacteremia, and meningitis, with case fatality rates of 50-80% being reported (Healy et al. 2010). Infant formula manufacturers …

Dynamic Herbal Monographs For A Digital World, Niamh O'Brien

Theses

Post analysis of a worldwide survey of Medical Herbalists, 93% of respondents were in favour of an online system which could update monographs dynamically. 63% of respondents suggested that some current monographs are out of date and lack certain practicalities in areas such as : Interactions, Dosage and Safety. Research into gaining optimal responses from surveys led to a 78% response rate. Survey analysis resulted in a modem up-to-date monograph template being created and each of the aforementioned information systems tested against same. Testing involved the generation of XML, HTML, PHP and OWL languages for encoding documents to allow for …

On The Global Stability Of A Generalized Cholera Epidemiological Model, Yuanji Cheng, Jin Wang, Xiuxiang Yang

Mathematics & Statistics Faculty Publications

In this paper, we conduct a careful global stability analysis for a generalized cholera epidemiological model originally proposed in [J. Wang and S. Liao, A generalized cholera model and epidemic/endemic analysis, J. Biol. Dyn. 6 (2012), pp. 568-589]. Cholera is a water-and food-borne infectious disease whose dynamics are complicated by the multiple interactions between the human host, the pathogen, and the environment. Using the geometric approach, we rigorously prove the endemic global stability for the cholera model in three-dimensional (when the pathogen component is a scalar) and four-dimensional (when the pathogen component is a vector) systems. This work unifies the …

Efficient And Robust Rna-Seq Process For Cultured Bacteria And Complex Community Transcriptomes, Georgia Giannoukos, Dawn M. Ciulla, Katherine Huang, Brian J. Haas, Jacques Izard, Joshua Z. Levin, Jonathan Livny, Ashlee M. Earl, Dirk Gevers, Doyle V. Ward, Chad Nusbaum, Bruce W. Birren, Andreas Gnirke

Department of Food Science and Technology: Faculty Publications

We have developed a process for transcriptome analysis of bacterial communities that accommodates both intact and fragmented starting RNA and combines efficient rRNA removal with strand-specific RNA-seq. We applied this approach to an RNA mixture derived from three diverse cultured bacterial species and to RNA isolated from clinical stool samples. The resulting expression profiles were highly reproducible, enriched up to 40-fold for non-rRNA transcripts, and correlated well with profiles representing undepleted total RNA.

Modeling Protein Expression And Protein Signaling Pathways, Donatello Telesca, Peter Muller, Steven Kornblau, Marc Suchard, Yuan Ji

COBRA Preprint Series

High-throughput functional proteomic technologies provide a way to quantify the expression of proteins of interest. Statistical inference centers on identifying the activation state of proteins and their patterns of molecular interaction formalized as dependence structure. Inference on dependence structure is particularly important when proteins are selected because they are part of a common molecular pathway. In that case inference on dependence structure reveals properties of the underlying pathway. We propose a probability model that represents molecular interactions at the level of hidden binary latent variables that can be interpreted as indicators for active versus inactive states of the proteins. The …

A Study Of Correlations Between The Definition And Application Of The Gene Ontology, Yuji Mo

Department of Computer Electronics and Engineering: Dissertations, Theses, and Student Research

When using the Gene Ontology (GO), nucleotide and amino acid sequences are annotated by terms in a structured and controlled vocabulary organized into relational graphs. The usage of the vocabulary (GO terms) in the annotation of these sequences may diverge from the relations defined in the ontology. We measure the consistency of the use of GO terms by comparing GO's defined structure to the terms' application. To do this, we first use synthetic data with different characteristics to understand how these characteristics influence the correlation values determined by various similarity measures. Using these results as a baseline, we found that …

Planning Combinatorial Disulfide Cross-Links For Protein Fold Determination, Fei Xiong, Alan M Friedman, Chris Bailey-Kellogg

Dartmouth Scholarship

Fold recognition techniques take advantage of the limited number of overall structural organizations, and have become increasingly effective at identifying the fold of a given target sequence. However, in the absence of sufficient sequence identity, it remains difficult for fold recognition methods to always select the correct model. While a native-like model is often among a pool of highly ranked models, it is not necessarily the highest-ranked one, and the model rankings depend sensitively on the scoring function used. Structure elucidation methods can then be employed to decide among the models based on relatively rapid biochemical/biophysical experiments.

Additive Functions In Boolean Models Of Gene Regulatory Network Modules, Christian Darabos, Ferdinando Ferdinando Di Cunto, Marco Tomassini, Jason H. Moore

Dartmouth Scholarship

Gene-on-gene regulations are key components of every living organism. Dynamical abstract models of genetic regulatory networks help explain the genome’s evolvability and robustness. These properties can be attributed to the structural topology of the graph formed by genes, as vertices, and regulatory interactions, as edges. Moreover, the actual gene interaction of each gene is believed to play a key role in the stability of the structure. With advances in biology, some effort was deployed to develop update functions in Boolean models that include recent knowledge. We combine real-life gene interaction networks with novel update functions in a Boolean model. We …

Transcriptomic Characterization Of A Synergistic Genetic Interaction During Carpel Margin Meristem Development In Arabidopsis Thaliana, April N. Wynn, Elizabeth E. Rueschhoff, Robert G. Franks

Biological Sciences Research

In flowering plants the gynoecium is the female reproductive structure. In Arabidopsis thalianaovules initiate within the developing gynoecium from meristematic tissue located along the margins of the floral carpels. When fertilized the ovules will develop into seeds. SEUSS (SEU) and AINTEGUMENTA (ANT) encode transcriptional regulators that are critical for the proper formation of ovules from the carpel margin meristem (CMM). The synergistic loss of ovule initiation observed in the seu ant double mutant suggests that SEU and ANT share overlapping functions during CMM development. However the molecular mechanism underlying this synergistic interaction is unknown. Using …

Gc-Content Normalization For Rna-Seq Data, Davide Risso, Katja Schwartz, Gavin Sherlock, Sandrine Dudoit

U.C. Berkeley Division of Biostatistics Working Paper Series

Background: Transcriptome sequencing (RNA-Seq) has become the assay of choice for high-throughput studies of gene expression. However, as is the case with microarrays, major technology-related artifacts and biases affect the resulting expression measures. Normalization is therefore essential to ensure accurate inference of expression levels and subsequent analyses thereof.

Results: We focus on biases related to GC-content and demonstrate the existence of strong sample-specific GC-content effects on RNA-Seq read counts, which can substantially bias differential expression analysis. We propose three simple within-lane gene-level GC-content normalization approaches and assess their performance on two different RNA-Seq datasets, involving different species and experimental designs. …

Magnetosome Genes In The Gammaproteobacterium Strain Bw-2, Lucero Rivera, Denis Trubitsyn, Dennis A. Bazylinski

Undergraduate Research Opportunities Program (UROP)

Magnetotactic bacteria (MTB) biomineralize intracellular nanometer-sized, magnetic crystals surrounded by a lipid bilayer membrane known as magnetosomes. These crystals, which consist of magnetite (Fe3O4) or greigite (Fe3S4), causes the cell to align along the geomagnetic field lines as they swim, a phenomenon known as magnetotaxis. Strain BW-2 is a magnetite-producing magnetotactic bacterium isolated from Badwater Basin, Death Valley National Park (California) and is one of only two species of MTB that are known to phylogenetically belong to the Gammaproteobacteria class of the Proteobacteria phylum. The biomineralization of magnetite in magnetotactic bacteria is mediated by a series of genes that include …

Multiple Testing Of Local Maxima For Detection Of Peaks In Chip-Seq Data, Armin Schwartzman, Andrew Jaffe, Yulia Gavrilov, Clifford A. Meyer

Harvard University Biostatistics Working Paper Series

No abstract provided.