Genetics and Genomics Commons^™

Computational Network Analysis Of The Anatomical And Genetic Organizations In The Mouse Brain, Shuiwang Ji

Computer Science Faculty Publications

Motivation: The mammalian central nervous system (CNS) generates high-level behavior and cognitive functions. Elucidating the anatomical and genetic organizations in the CNS is a key step toward understanding the functional brain circuitry. The CNS contains an enormous number of cell types, each with unique gene expression patterns. Therefore, it is of central importance to capture the spatial expression patterns in the brain. Currently, genome-wide atlas of spatial expression patterns in the mouse brain has been made available, and the data are in the form of aligned 3D data arrays. The sheer volume and complexity of these data pose significant challenges …

Minimum Description Length Measures Of Evidence For Enrichment, Zhenyu Yang, David R. Bickel

COBRA Preprint Series

In order to functionally interpret differentially expressed genes or other discovered features, researchers seek to detect enrichment in the form of overrepresentation of discovered features associated with a biological process. Most enrichment methods treat the p-value as the measure of evidence using a statistical test such as the binomial test, Fisher's exact test or the hypergeometric test. However, the p-value is not interpretable as a measure of evidence apart from adjustments in light of the sample size. As a measure of evidence supporting one hypothesis over the other, the Bayes factor (BF) overcomes this drawback of the p-value but lacks …

Appearance Based Stage Recognition Of Drosophila Embryos, Gopi Chand Nutakki

Masters Theses & Specialist Projects

Stages in Drosophila development denote the time after fertilization at which certain specific events occur in the developmental cycle. Stage information of a host embryo, as well as spatial information of a gene expression region is indispensable input for the discovery of the pattern of gene-gene interaction. Manual labeling of stages is becoming a bottleneck under the circumstance of high throughput embryo images. Automatic recognition based on the appearances of embryos is becoming a more desirable scheme. This problem, however, is very challenging due to severe variations of illumination and gene expressions. In this research thesis, we propose an appearance …

A Bayesian Shared Component Model For Genetic Association Studies, Juan J. Abellan, Carlos Abellan, Juan R. Gonzalez

COBRA Preprint Series

We present a novel approach to address genome association studies between single nucleotide polymorphisms (SNPs) and disease. We propose a Bayesian shared component model to tease out the genotype information that is common to cases and controls from the one that is specific to cases only. This allows to detect the SNPs that show the strongest association with the disease. The model can be applied to case-control studies with more than one disease. In fact, we illustrate the use of this model with a dataset of 23,418 SNPs from a case-control study by The Welcome Trust Case Control Consortium (2007) …

Minimum Description Length And Empirical Bayes Methods Of Identifying Snps Associated With Disease, Ye Yang, David R. Bickel

COBRA Preprint Series

The goal of determining which of hundreds of thousands of SNPs are associated with disease poses one of the most challenging multiple testing problems. Using the empirical Bayes approach, the local false discovery rate (LFDR) estimated using popular semiparametric models has enjoyed success in simultaneous inference. However, the estimated LFDR can be biased because the semiparametric approach tends to overestimate the proportion of the non-associated single nucleotide polymorphisms (SNPs). One of the negative consequences is that, like conventional p-values, such LFDR estimates cannot quantify the amount of information in the data that favors the null hypothesis of no disease-association.

We …

Computational Biology, Harvey Greenberg, Allen Holder

Mathematical Sciences Technical Reports (MSTR)

Computational biology is an interdisciplinary field that applies the techniques of computer science, applied mathematics, and statistics to address biological questions. OR is also interdisciplinary and applies the same mathematical and computational sciences, but to decision-making problems. Both focus on developing mathematical models and designing algorithms to solve them. Models in computational biology vary in their biological domain and can range from the interactions of genes and proteins to the relationships among organisms and species.

G-Lattices For An Unrooted Perfect Phylogeny, Monica Grigg

Mathematical Sciences Technical Reports (MSTR)

We look at the Pure Parsimony problem and the Perfect Phylogeny Haplotyping problem. From the Pure Parsimony problem we consider structures of genotypes called g-lattices. These structures either provide solutions or give bounds to the pure parsimony problem. In particular, we investigate which of these structures supports an unrooted perfect phylogeny, a condition that adds biological interpretation. By understanding which g-lattices support an unrooted perfect phylogeny, we connect two of the standard biological inference rules used to recreate how genetic diversity propagates across generations.

A Perturbation Method For Inference On Regularized Regression Estimates, Jessica Minnier, Lu Tian, Tianxi Cai

Harvard University Biostatistics Working Paper Series

No abstract provided.

The Strength Of Statistical Evidence For Composite Hypotheses: Inference To The Best Explanation, David R. Bickel

COBRA Preprint Series

A general function to quantify the weight of evidence in a sample of data for one hypothesis over another is derived from the law of likelihood and from a statistical formalization of inference to the best explanation. For a fixed parameter of interest, the resulting weight of evidence that favors one composite hypothesis over another is the likelihood ratio using the parameter value consistent with each hypothesis that maximizes the likelihood function over the parameter of interest. Since the weight of evidence is generally only known up to a nuisance parameter, it is approximated by replacing the likelihood function with …

Powerful Snp Set Analysis For Case-Control Genome Wide Association Studies, Michael C. Wu, Peter Kraft, Michael P. Epstein, Deanne M. Taylor, Stephen J. Chanock, David J. Hunter, Xihong Lin

Harvard University Biostatistics Working Paper Series

No abstract provided.

An Entrepreneurial Approach To Librarianship, Flora Shrode, Jennifer Duncan, Wendy Holliday

Library Faculty & Staff Publications

Librarians from Utah State University explain recent efforts to encourage subject librarians to take a more holistic view of their roles. We are shifting from a traditional emphasis primarily on collection development and refocusing on natural connections between collections, instruction, liaison, and reference service. The poster provides background about Utah State University’s situation and explains our approach to analyzing local needs and culture to inform development of a new organizational structure. We describe our vision of subject librarianship, the process by which we assessed librarians’ ideas and goals for performing as subject librarians, and the actions we are taking to …

Quasispecies-Like Behavior Observed In Catalytic Rna Populations Evolving In A Test Tube, Carolina Diaz Arenas, Niles Lehman

Chemistry Faculty Publications and Presentations

Background: During the RNA World, molecular populations were probably very small and highly susceptible to the force of strong random drift. In conjunction with Muller's Ratchet, this would have imposed difficulties for the preservation of the genetic information and the survival of the populations. Mechanisms that allowed these nascent populations to overcome this problem must have been advantageous.

Results: Using continuous in vitro evolution experimentation with an increased mutation rate imposed by MnCl2, it was found that clonal 100-molecule populations of ribozymes clearly exhibit certain characteristics of a quasispecies. This is the first time this has been seen …

Sex, Mixability, And Modularity, Adi Livnat, Christos Papadimitriou, Nicholas Pippenger, Marcus W. Feldman

All HMC Faculty Publications and Research

The assumption that different genetic elements can make separate contributions to the same quantitative trait was originally made in order to reconcile biometry and Mendelism and ever since has been used in population genetics, specifically for the trait of fitness. Here we show that sex is responsible for the existence of separate genetic effects on fitness and, more generally, for the existence of a hierarchy of genetic evolutionary modules. Using the tools developed in the process, we also demonstrate that in terms of their fitness effects, separation and fusion of genes are associated with the increase and decrease of the …

Ancient Dna Analyses Exclude Humans As The Driving Force Behind Late Pleistocene Musk Ox (Ovibos Moschatus) Population Dynamics, Paula F. Campos, Eske Willerslev, Andrei Sher, Ludovic Orlando, Erik Axelsson, Alexei Tikhonov, Kim Aaris-Sorensen, Alex D. Greenwood, Ralf-Dietrich Kahlke, Pavel Kosintsev

Biological Sciences Faculty Publications

The causes of the late Pleistocene megafaunal extinctions are poorly understood. Different lines of evidence point to climate change, the arrival of humans, or a combination of these events as the trigger. Although many species went extinct, others, such as caribou and bison, survived to the present. The musk ox has an intermediate story: relatively abundant during the Pleistocene, it is now restricted to Greenland and the Arctic Archipelago. In this study, we use ancient DNA sequences, temporally unbiased summary statistics, and Bayesian analytical techniques to infer musk ox population dynamics throughout the late Pleistocene and Holocene. Our results reveal …

An Intelligent Data-Centric Approach Toward Identification Of Conserved Motifs In Protein Sequences, Kathryn Dempsey Cooper, Benjamin Currall, Richard Hallworth, Hesham Ali

Interdisciplinary Informatics Faculty Proceedings & Presentations

The continued integration of the computational and biological sciences has revolutionized genomic and proteomic studies. However, efficient collaboration between these fields requires the creation of shared standards. A common problem arises when biological input does not properly fit the expectations of the algorithm, which can result in misinterpretation of the output. This potential confounding of input/output is a drawback especially when regarding motif finding software. Here we propose a method for improving output by selecting input based upon evolutionary distance, domain architecture, and known function. This method improved detection of both known and unknown motifs in two separate case studies. …

Genbank, Dennis A. Benson, Ilene Karasch-Mizrachi, David J. Lipman, James Ostell, Eric W. Sayers

Harold W. Manter Laboratory of Parasitology: Library Materials

GenBank(R) is a comprehensive database that contains publicly available nucleotide sequences for more than 380,000 organisms named at the genus level or lower, obtained primarily through submissions from individual laboratories and batch submissions from large-scale sequencing projects, including whole genome shotgun (WGS) and environmental sampling projects. Most submissions are made using the web-based BankIt or standalone Sequin programs, and accession numbers are assigned by GenBank staff upon receipt. Daily data exchange with the European Nucleotide Archive (ENA) and the DNA Data Bank of Japan (DDBJ) ensures worldwide coverage. GenBank is accessible through the NCBI Entrez retrieval system that integrates data …

A Decomposition Of The Pure Parsimony Problem, Allen Holder, Thomas M. Langley

Mathematical Sciences Technical Reports (MSTR)

We partially order a collection of genotypes so that we can represent the problem of inferring the least number of haplotypes in terms of substructures we call g-lattices. This representation allows us to prove that if the genotypes partition into chains with certain structure, then the NP-Hard problem can be solved efficiently. Even without the specified structure, the decomposition shows how to separate the underlying integer programming model into smaller models.

Fisheries Management Paper No.245 - Aquaculture Of Coral, Live Rocks And Associated Products. Aquaculture Policy, Department Of Fisheries, Western Australia

Fisheries Management Papers

Aquarium enthusiasts use corals, live rock and associated products (e.g. live sand) for display and to maintain healthy water quality conditions in their marine aquariums. There is a large and growing market for these products. While the wild harvest fisheries for aquarium products in Western Australia are sustainably managed, output is capped, so increased product is needed to meet the market demand. With a trend towards the capping of wild harvest and additional pressures, such as the introduction of Marine Parks and other marine protected areas, the aquaculture of these products is set for industry development and expansion.

The large …

Failure To Replicate A Genetic Association May Provide Important Clues About Genetic Architecture, Casey S. Greene, Nadia M. Penrod, Scott M. Williams, Jason H. Moore

Dartmouth Scholarship

Replication has become the gold standard for assessing statistical results from genome-wide association studies. Unfortunately this replication requirement may cause real genetic effects to be missed. A real result can fail to replicate for numerous reasons including inadequate sample size or variability in phenotype definitions across independent samples. In genome-wide association studies the allele frequencies of polymorphisms may differ due to sampling error or population differences. We hypothesize that some statistically significant independent genetic effects may fail to replicate in an independent dataset when allele frequencies differ and the functional polymorphism interacts with one or more other functional polymorphisms. To …

Resampling-Based Multiple Hypothesis Testing With Applications To Genomics: New Developments In The R/Bioconductor Package Multtest, Houston N. Gilbert, Katherine S. Pollard, Mark J. Van Der Laan, Sandrine Dudoit

U.C. Berkeley Division of Biostatistics Working Paper Series

The multtest package is a standard Bioconductor package containing a suite of functions useful for executing, summarizing, and displaying the results from a wide variety of multiple testing procedures (MTPs). In addition to many popular MTPs, the central methodological focus of the multtest package is the implementation of powerful joint multiple testing procedures. Joint MTPs are able to account for the dependencies between test statistics by effectively making use of (estimates of) the test statistics joint null distribution. To this end, two additional bootstrap-based estimates of the test statistics joint null distribution have been developed for use in the …

Joint Multiple Testing Procedures For Graphical Model Selection With Applications To Biological Networks, Houston N. Gilbert, Mark J. Van Der Laan, Sandrine Dudoit

U.C. Berkeley Division of Biostatistics Working Paper Series

Gaussian graphical models have become popular tools for identifying relationships between genes when analyzing microarray expression data. In the classical undirected Gaussian graphical model setting, conditional independence relationships can be inferred from partial correlations obtained from the concentration matrix (= inverse covariance matrix) when the sample size n exceeds the number of parameters p which need to estimated. In situations where n < p, another approach to graphical model estimation may rely on calculating unconditional (zero-order) and first-order partial correlations. In these settings, the goal is to identify a lower-order conditional independence graph, sometimes referred to as a ‘0-1 graphs’. For either choice of graph, model selection may involve a multiple testing problem, in which edges in a graph are drawn only after rejecting hypotheses involving (saturated or lower-order) partial correlation parameters. Most multiple testing procedures applied in previously proposed graphical model selection algorithms rely on standard, marginal testing methods which do not take into account the joint distribution of the test statistics derived from (partial) correlations. We propose and implement a multiple testing framework useful when testing for edge inclusion during graphical model selection. Two features of our methodology include (i) a computationally efficient and asymptotically valid test statistics joint null distribution derived from influence curves for correlation-based parameters, and (ii) the application of empirical Bayes joint multiple testing procedures which can effectively control a variety of popular Type I error rates by incorpo- rating joint null distributions such as those described here (Dudoit and van der Laan, 2008). Using a dataset from Arabidopsis thaliana, we observe that the use of more sophisticated, modular approaches to multiple testing allows one to identify greater numbers of edges when approximating an undirected graphical model using a 0-1 graph. Our framework may also be extended to edge testing algorithms for other types of graphical models (e.g., for classical undirected, bidirected, and directed acyclic graphs).

One Rna Plays Three Roles To Provide Catalytic Activity To A Group I Intron Lacking An Endogenous Internal Guide Sequence, Nilesh Vaidya, Niles Lehman

Chemistry Faculty Publications and Presentations

Catalytic RNA molecules possess simultaneously a genotype and a phenotype. However, a single RNA genotype has the potential to adopt two or perhaps more distinct phenotypes as a result of differential folding and/or catalytic activity. Such multifunctionality would be particularly significant if the phenotypes were functionally inter-related in a common biochemical pathway. Here, this phenomenon is demonstrated by the ability of the Azoarcus group I ribozyme to function when its canonical internal guide sequence (GUG) has been removed from the 5? end of the molecule, and added back exogenously in trans. The presence of GUG triplets in noncovalent fragments of …

Sparse Linear Discriminant Analysis For Simultaneous Testing For The Significance Of A Gene Set/Pathway And Gene Selection, Michael C. Wu, Lingson Zhang, Zhaoxi Wang, David C. Christiani, Xihong Lin

Harvard University Biostatistics Working Paper Series

No abstract provided.

High Resolution Single Molecule Optical Localization Of Multiple Fluorophores On Dna Origami Constructs Fluorophores On Dna Origami Constructs, Anuradha Rajulapati

Theses, Dissertations and Capstones

An ongoing challenge in the development of nanoelectronics and nanophotonics is the nondestructive, high-resolution localization in space of single molecules and multi-molecular assemblies. The apparent barrier to the use of optical microscopy at the sub-100 nm scale is the well known Abbe Limit, the diffraction limit to resolution. This laboratory has adapted a technique called Single-molecule high resolution imaging with photobleaching (SHRIMP). We have developed methods for utilizing SHRIMP for the determination of the separation of two fluorophores in single DNA origami constructs.

DNA Origami is extremely useful because it can address nanocomponents down to 2nm separation. It should be …

Understanding The Physical Properties That Control Protein Crystallization By Analysis Of Largescale Experimental Data, W. Nicholson Price Ii, Yang Chen, Samuel K. Handelman, Helen Neely, Philip Manor, Richard Karlin, Rajesh Nair, Jinfeng Liu, Michael Baran, John Everett, Saichiu N. Tong, Farhad Forouhar, Swarup S. Swaminathan, Thomas Acton, Rong Xiao, Joseph R. Luft, Angela Lauricella, George T. Detitta, Burkhard Rost, Gaetano T. Montelione, John T. Hunt

Law Faculty Scholarship

Crystallization is the most serious bottleneck in high-throughput protein-structure determination by diffraction methods. We have used data mining of the large-scale experimental results of the Northeast Structural Genomics Consortium and experimental folding studies to characterize the biophysical properties that control protein crystallization. This analysis leads to the conclusion that crystallization propensity depends primarily on the prevalence of well-ordered surface epitopes capable of mediating interprotein interactions and is not strongly influenced by overall thermodynamic stability. We identify specific sequence features that correlate with crystallization propensity and that can be used to estimate the crystallization probability of a given construct. Analyses of …

Multiple Imputation To Correct For Measurement Error In Admixture Estimates In Genetic Structured Association Testing, Miguel A. Padilla, Jamin Divers, Laura K. Vaughan, David B. Allison, Hemant K. Tiwari

Psychology Faculty Publications

Objectives: Structured association tests ( SAT), like any statistical model, assumes that all variables are measured without error. Measurement error can bias parameter estimates and confound residual variance in linear models. It has been shown that admixture estimates can be contaminated with measurement error causing SAT models to suffer from the same afflictions. Multiple imputation (MI) is presented as a viable tool for correcting measurement error problems in SAT linear models with emphasis on correcting measurement error contaminated admixture estimates. Methods: Several MI methods are presented and compared, via simulation, in terms of controlling Type I error rates for both …

The Strength Of Statistical Evidence For Composite Hypotheses With An Application To Multiple Comparisons, David R. Bickel

COBRA Preprint Series

The strength of the statistical evidence in a sample of data that favors one composite hypothesis over another may be quantified by the likelihood ratio using the parameter value consistent with each hypothesis that maximizes the likelihood function. Unlike the p-value and the Bayes factor, this measure of evidence is coherent in the sense that it cannot support a hypothesis over any hypothesis that it entails. Further, when comparing the hypothesis that the parameter lies outside a non-trivial interval to the hypotheses that it lies within the interval, the proposed measure of evidence almost always asymptotically favors the correct hypothesis …

Sheep Updates 2008 - Part 3, L. J. E. Karlsson, J. C. Greeff, L. Slocombe, K. Jones, N. Underwood, Fred Wilkinson, W. D. Hoffman, W. A. Mckiernan, V. H. Oddy, M. J. Mcphee, B. L. Mcintyre, P. F. Parnell, R. A. Clark, J. Timms, G. Griffith, C. Mulholland, P. Hyland, Danielle England, Fiona Jones, John Lucey, Martin Staines, Richard Morris, Megan Abrahams, Caroline Peek, Dennis Van Gool, Daniel Gardiner, Kari-Lee Falconer, Sandra Prosser, Mario D'Antuono, John Young, Andrew Thompson, Chris Oldham, Brown Besier, Angus Campbell, Ralph Behrendt

Sheep Updates

This session covers fiveteen papers from different authors:

CONTROLLING FLY STRIKE

1. Breeding for Blowfly Resistance - Indicatoe Traits, LJE Karlsson, JC Greeff, L Slocombe, Department of Agriculture & Food, Western Australia

2.A practical method to select for breech strike resistance in non-pedigreed Merino flocks, LJE Karlsson, JC Greeff, L Slocombe, K. Jones, N. Underwood, Department of Agriculture & Food, Western Australia

3. Twice a year shearing - no mulesing, Fred Wilkinson, Producer, Brookton WA

BEEF

4. Commercial testing of a new tool for prediction of fatness in beef cattle, WD Hoffman^A, WA McKiernan^A, VH Oddy …

Estimation And Testing For The Effect Of A Genetic Pathway On A Disease Outcome Using Logistic Kernel Machine Regression Via Logistic Mixed Models, Dawei Liu, Debashis Ghosh, Xihong Lin

Harvard University Biostatistics Working Paper Series

No abstract provided.

A Powerful And Flexible Multilocus Association Test For Quantitative Traits, Lydia Coulter Kwee, Dawei Liu, Xihong Lin, Debashis Ghosh, Michael P. Epstein

Harvard University Biostatistics Working Paper Series

No abstract provided.