Multivariate point process models

Size: px
Start display at page:

Download "Multivariate point process models"

Transcription

1 Faculty of Science Multivariate point process models Niels Richard Hansen Department of Mathematical Sciences January 8, 200 Slide /20

2 Ideas and outline General aim: To build and implement a flexible (non-parametric), intensity based framework of multivariate point process models aka discretely marked point processes. Secondary aim: Development and implementation of model selection techniques LASSO, grouped LASSO, penalization based ideas with the multivariate point processes as a testbed. Joint with Patricia Bouret, Vincent Rivoirard. Main application: Organizational model of active transcription regulatory elements along genomes joint with Lisbeth Carstensen, Albin Sandelin, Ole Winther. Slide 2/20 Niels Richard Hansen Multivariate point process models January 8, 200

3 Genomic scales an e- se an c- s, ister re la- - n- in of ay ral s- ra- We focus on the distribution of point-like transcriptional regulatory elements at the meso-genomic scale. Figure from Zhang et al, Genome Res., 8-, 200 Slide /20 Niels Richard Hansen Multivariate point process models January 8, 200

4 Transcription regulator binding loci Attempt of a broad definition: Transcription regulators are proteins that modify, interact with or bind to the DNA, chromatin or other transcription regulators to either activate or repress the transcription of DNA. An active transcription regulatory loci is a loci on the genome where we observe the presence of a transcription regulator. Fact: Transcription regulators cluster in promoter regions and intergenic regions. Why? Is there a combined effect? Do they recruit each other...? With the different regulators as marks and with measurements of the active loci as points on the meso-genomic scale we use a multivariate point process model of the organization of active loci. Slide 4/20 Niels Richard Hansen Multivariate point process models January 8, 200

5 Obtaining and preprocessing data (ChIP-chip) TGGCA ACCGT GCAAC CGTTG cutoff actual binding sites Slide /20 Niels Richard Hansen Multivariate point process models January 8, 200

6 ENCODE affymetrix data sirt rara pu pol2 type brg cebpe ctcf type p00 h4kac4histone hk2mehistone hk2mehistone h4kac4histone p00 pol2 ctcf cebpe brg position pu rara sirt The 0 Affymetrix ChIP-chip measurements for the ENCODE pilot project segment ENm00. Histone modifications are not regarded as point-like, figure gives start positions. Slide 6/20 Niels Richard Hansen Multivariate point process models January 8, 200

7 Embryonic mouse stem cell data (ChIP-seq) Suz2 STAT type Sox2 Nanog Oct4 type Smad p00 Smad p00 Sox2 STAT Oct4 Suz2 Nanog.0e+0.0e+08.e+08 position Mouse, chromosome I: of active transcription factor binding loci measured by ChIP-seq for embryonic stem cells. 2 2 Chen et al, Cell, 06-, 2008 Slide /20 Niels Richard Hansen Multivariate point process models January 8, 200

8 The Hawkes process If t i <... < t ini denote the observed points for transcription regulator i we specify the intensity process for the occurrence of regulator j as λ j (t) = ϕ(β 0 + h i,j (t t ik )) i k:t ik <t for a fixed ϕ. We expand the h-function in a suitable basis; B i,j h i,j (s) = β i,j b B b(s) b= and the minus-log-likelihood in terms of β = (β 0, β ij b ) is directly expressible. Slide 8/20 Niels Richard Hansen Multivariate point process models January 8, 200

9 ppstat The current implementation in the R package ppstat offers Formula interface to model specification. Standard (e.g. spline) basis function expansions of linear filters. Inclusion of continuous time covariate effects and additive model specification. Currently only with ϕ = exp. Soon with a more flexible class of ϕ s including ϕ(x) = x +. Why is ϕ(x) = x + a problem? The minus-log-likelihood is a priory not twice differentiable, but it turns out to be so under certain regularity assumptions. Slide /20 Niels Richard Hansen Multivariate point process models January 8, 200

10 Some spline bases B(t) t B(t) t Slide 0/20 Niels Richard Hansen Multivariate point process models January 8, 200

11 university of copenhagen department of mathematical sciences Estimated multiplicative effects - ES cells E2f Zfx c.myc n.myc Nanog Sox2 Oct4 Suz2 STAT p00 Smad E2f Zfx c.myc n.myc Nanog Sox2 Oct4 Suz2 STAT p00 Smad position Slide /20 Niels Richard Hansen Multivariate point process models January 8,

12 Interaction terms The formula i k:t ik <t h i,j (t t ik ) that entered in the intensity is an additive model of linear, univariate filters. Question: Can we capture interaction effects? Suggestion: Introduce bivariate, linear filters h (i,r),j (t t ik, t t rl ) k:t ik <t l:t rl <t and expand them in a suitable bivariate basis. Slide 2/20 Niels Richard Hansen Multivariate point process models January 8, 200

13 Bivariate bases Many possible choices of 2d-bases: Thin plate splines h(s, t) = (s s 0 ) 2 + (t t 0 ) 2 log((s s 0 ) 2 + (t t 0 ) 2 ). Tensor products h(s, t) = h (s)h 2 (t). For illustration we consider an example with a spline based tensor product construction with 6 basis functions. Slide /20 Niels Richard Hansen Multivariate point process models January 8, 200

14 Bivariate spline tensor product bases value Sox Oct4 Slide 4/20 Niels Richard Hansen Multivariate point process models January 8, 200

15 Additive estimated Oct4-Sox2 effect on Nanog what joint oct4 sox2 Sox2 400 value Oct4 ĥ Oct4,Nanog (s) + ĥ Sox2,Nanog (t) Slide /20 Niels Richard Hansen Multivariate point process models January 8, 200

16 Estimated Oct4-Sox2 interaction effect on Nanog what joint oct4 sox2 Sox2 400 value Oct4 ĥ Oct4,Nanog (s) + ĥ Sox2,Nanog (t) + ĥ (Sox2,Oct4),Nanog (s, t) Slide 6/20 Niels Richard Hansen Multivariate point process models January 8, 200

17 Difference what joint oct4 sox2 Sox2 400 value Oct4 ĥ (Sox2,Oct4),Nanog (s, t) Slide /20 Niels Richard Hansen Multivariate point process models January 8, 200

18 What do we model? The point process model is a model of the organizational part of evolution. It is a model of which loci are active for a given cell line/organism/tissue under given experimental conditions. The driving stochastic mechanism is evolution we do not account for measurement noise or biological variation in an explicit way. Data are integrated over cells we do not observe single cell transcription factor bindings. Slide 8/20 Niels Richard Hansen Multivariate point process models January 8, 200

19 Which conclusions do the model facilitate? Working assumption: Organizational association Functional association It is a systemic analysis that facilitate phrasing and answering question of a systemic nature. Example: For two tissues find systemic differences in the organization of active regulatory loci. Slide /20 Niels Richard Hansen Multivariate point process models January 8, 200

20 Thanks to Lisbeth Carstensen, Albin Sandelin and Ole Winther and to all of you for attending this meeting. Slide 20/20 Niels Richard Hansen Multivariate point process models January 8, 200

Complete all warm up questions Focus on operon functioning we will be creating operon models on Monday

Complete all warm up questions Focus on operon functioning we will be creating operon models on Monday Complete all warm up questions Focus on operon functioning we will be creating operon models on Monday 1. What is the Central Dogma? 2. How does prokaryotic DNA compare to eukaryotic DNA? 3. How is DNA

More information

Chapter 15 Active Reading Guide Regulation of Gene Expression

Chapter 15 Active Reading Guide Regulation of Gene Expression Name: AP Biology Mr. Croft Chapter 15 Active Reading Guide Regulation of Gene Expression The overview for Chapter 15 introduces the idea that while all cells of an organism have all genes in the genome,

More information

REVIEW SESSION. Wednesday, September 15 5:30 PM SHANTZ 242 E

REVIEW SESSION. Wednesday, September 15 5:30 PM SHANTZ 242 E REVIEW SESSION Wednesday, September 15 5:30 PM SHANTZ 242 E Gene Regulation Gene Regulation Gene expression can be turned on, turned off, turned up or turned down! For example, as test time approaches,

More information

Geert Geeven. April 14, 2010

Geert Geeven. April 14, 2010 iction of Gene Regulatory Interactions NDNS+ Workshop April 14, 2010 Today s talk - Outline Outline Biological Background Construction of Predictors The main aim of my project is to better understand the

More information

Biology. Biology. Slide 1 of 26. End Show. Copyright Pearson Prentice Hall

Biology. Biology. Slide 1 of 26. End Show. Copyright Pearson Prentice Hall Biology Biology 1 of 26 Fruit fly chromosome 12-5 Gene Regulation Mouse chromosomes Fruit fly embryo Mouse embryo Adult fruit fly Adult mouse 2 of 26 Gene Regulation: An Example Gene Regulation: An Example

More information

Measuring TF-DNA interactions

Measuring TF-DNA interactions Measuring TF-DNA interactions How is Biological Complexity Achieved? Mediated by Transcription Factors (TFs) 2 Regulation of Gene Expression by Transcription Factors TF trans-acting factors TF TF TF TF

More information

12-5 Gene Regulation

12-5 Gene Regulation 12-5 Gene Regulation Fruit fly chromosome 12-5 Gene Regulation Mouse chromosomes Fruit fly embryo Mouse embryo Adult fruit fly Adult mouse 1 of 26 12-5 Gene Regulation Gene Regulation: An Example Gene

More information

Control of Gene Expression

Control of Gene Expression Control of Gene Expression Mechanisms of Gene Control Gene Control in Eukaryotes Master Genes Gene Control In Prokaryotes Epigenetics Gene Expression The overall process by which information flows from

More information

Alignment. Peak Detection

Alignment. Peak Detection ChIP seq ChIP Seq Hongkai Ji et al. Nature Biotechnology 26: 1293-1300. 2008 ChIP Seq Analysis Alignment Peak Detection Annotation Visualization Sequence Analysis Motif Analysis Alignment ELAND Bowtie

More information

FYTN05/TEK267 Chemical Forces and Self Assembly

FYTN05/TEK267 Chemical Forces and Self Assembly FYTN05/TEK267 Chemical Forces and Self Assembly FYTN05/TEK267 Chemical Forces and Self Assembly 1 Michaelis-Menten (10.3,10.4) M-M formalism can be used in many contexts, e.g. gene regulation, protein

More information

The Research Plan. Functional Genomics Research Stream. Transcription Factors. Tuning In Is A Good Idea

The Research Plan. Functional Genomics Research Stream. Transcription Factors. Tuning In Is A Good Idea Functional Genomics Research Stream The Research Plan Tuning In Is A Good Idea Research Meeting: March 23, 2010 The Road to Publication Transcription Factors Protein that binds specific DNA sequences controlling

More information

Lecture 18 June 2 nd, Gene Expression Regulation Mutations

Lecture 18 June 2 nd, Gene Expression Regulation Mutations Lecture 18 June 2 nd, 2016 Gene Expression Regulation Mutations From Gene to Protein Central Dogma Replication DNA RNA PROTEIN Transcription Translation RNA Viruses: genome is RNA Reverse Transcriptase

More information

Kernel Density Estimation

Kernel Density Estimation Kernel Density Estimation If Y {1,..., K} and g k denotes the density for the conditional distribution of X given Y = k the Bayes classifier is f (x) = argmax π k g k (x) k If ĝ k for k = 1,..., K are

More information

Comparing transcription factor regulatory networks of human cell types. The Protein Network Workshop June 8 12, 2015

Comparing transcription factor regulatory networks of human cell types. The Protein Network Workshop June 8 12, 2015 Comparing transcription factor regulatory networks of human cell types The Protein Network Workshop June 8 12, 2015 KWOK-PUI CHOI Dept of Statistics & Applied Probability, Dept of Mathematics, NUS OUTLINE

More information

Genome 541! Unit 4, lecture 2! Transcription factor binding using functional genomics

Genome 541! Unit 4, lecture 2! Transcription factor binding using functional genomics Genome 541 Unit 4, lecture 2 Transcription factor binding using functional genomics Slides vs chalk talk: I m not sure why you chose a chalk talk over ppt. I prefer the latter no issues with readability

More information

3.B.1 Gene Regulation. Gene regulation results in differential gene expression, leading to cell specialization.

3.B.1 Gene Regulation. Gene regulation results in differential gene expression, leading to cell specialization. 3.B.1 Gene Regulation Gene regulation results in differential gene expression, leading to cell specialization. We will focus on gene regulation in prokaryotes first. Gene regulation accounts for some of

More information

Computational Biology: Basics & Interesting Problems

Computational Biology: Basics & Interesting Problems Computational Biology: Basics & Interesting Problems Summary Sources of information Biological concepts: structure & terminology Sequencing Gene finding Protein structure prediction Sources of information

More information

Genome 541 Gene regulation and epigenomics Lecture 2 Transcription factor binding using functional genomics

Genome 541 Gene regulation and epigenomics Lecture 2 Transcription factor binding using functional genomics Genome 541 Gene regulation and epigenomics Lecture 2 Transcription factor binding using functional genomics I believe it is helpful to number your slides for easy reference. It's been a while since I took

More information

Genome 541! Unit 4, lecture 3! Genomics assays

Genome 541! Unit 4, lecture 3! Genomics assays Genome 541! Unit 4, lecture 3! Genomics assays Much easier to follow with slides. Good pace.! Having the slides was really helpful clearer to read and easier to follow the trajectory of the lecture.!!

More information

Computational Genomics. Systems biology. Putting it together: Data integration using graphical models

Computational Genomics. Systems biology. Putting it together: Data integration using graphical models 02-710 Computational Genomics Systems biology Putting it together: Data integration using graphical models High throughput data So far in this class we discussed several different types of high throughput

More information

System Biology - Deterministic & Stochastic Dynamical Systems

System Biology - Deterministic & Stochastic Dynamical Systems System Biology - Deterministic & Stochastic Dynamical Systems System Biology - Deterministic & Stochastic Dynamical Systems 1 The Cell System Biology - Deterministic & Stochastic Dynamical Systems 2 The

More information

Welcome to Class 21!

Welcome to Class 21! Welcome to Class 21! Introductory Biochemistry! Lecture 21: Outline and Objectives l Regulation of Gene Expression in Prokaryotes! l transcriptional regulation! l principles! l lac operon! l trp attenuation!

More information

Modelling gene expression dynamics with Gaussian processes

Modelling gene expression dynamics with Gaussian processes Modelling gene expression dynamics with Gaussian processes Regulatory Genomics and Epigenomics March th 6 Magnus Rattray Faculty of Life Sciences University of Manchester Talk Outline Introduction to Gaussian

More information

A segmentation-clustering problem for the analysis of array CGH data

A segmentation-clustering problem for the analysis of array CGH data A segmentation-clustering problem for the analysis of array CGH data F. Picard, S. Robin, E. Lebarbier, J-J. Daudin UMR INA P-G / ENGREF / INRA MIA 518 APPLIED STOCHASTIC MODELS AND DATA ANALYSIS Brest

More information

Newly made RNA is called primary transcript and is modified in three ways before leaving the nucleus:

Newly made RNA is called primary transcript and is modified in three ways before leaving the nucleus: m Eukaryotic mrna processing Newly made RNA is called primary transcript and is modified in three ways before leaving the nucleus: Cap structure a modified guanine base is added to the 5 end. Poly-A tail

More information

Learning Sequence Motif Models Using Expectation Maximization (EM) and Gibbs Sampling

Learning Sequence Motif Models Using Expectation Maximization (EM) and Gibbs Sampling Learning Sequence Motif Models Using Expectation Maximization (EM) and Gibbs Sampling BMI/CS 776 www.biostat.wisc.edu/bmi776/ Spring 009 Mark Craven craven@biostat.wisc.edu Sequence Motifs what is a sequence

More information

INTEGRATING EPIGENETIC PRIORS FOR IMPROVING COMPUTATIONAL IDENTIFICATION OF TRANSCRIPTION FACTOR BINDING SITES AFFAN SHOUKAT

INTEGRATING EPIGENETIC PRIORS FOR IMPROVING COMPUTATIONAL IDENTIFICATION OF TRANSCRIPTION FACTOR BINDING SITES AFFAN SHOUKAT INTEGRATING EPIGENETIC PRIORS FOR IMPROVING COMPUTATIONAL IDENTIFICATION OF TRANSCRIPTION FACTOR BINDING SITES AFFAN SHOUKAT A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILMENT

More information

Network Biology-part II

Network Biology-part II Network Biology-part II Jun Zhu, Ph. D. Professor of Genomics and Genetic Sciences Icahn Institute of Genomics and Multi-scale Biology The Tisch Cancer Institute Icahn Medical School at Mount Sinai New

More information

Deciphering regulatory networks by promoter sequence analysis

Deciphering regulatory networks by promoter sequence analysis Bioinformatics Workshop 2009 Interpreting Gene Lists from -omics Studies Deciphering regulatory networks by promoter sequence analysis Elodie Portales-Casamar University of British Columbia www.cisreg.ca

More information

Variable Selection in Structured High-dimensional Covariate Spaces

Variable Selection in Structured High-dimensional Covariate Spaces Variable Selection in Structured High-dimensional Covariate Spaces Fan Li 1 Nancy Zhang 2 1 Department of Health Care Policy Harvard University 2 Department of Statistics Stanford University May 14 2007

More information

Functional Genomics Research Stream. Lecture: May 5, 2009 The Road to Publication

Functional Genomics Research Stream. Lecture: May 5, 2009 The Road to Publication Functional Genomics Research Stream Lecture: May 5, 2009 The Road to Publication Final Lab Issues Lab Hours this Week: Wednesday 10 AM - 6 PM Thursday 10 AM - 6 PM One Hour Required (sign in, sign out)

More information

ChIP-seq. Expression Networks. Genomics: Lecture #14. Institut für Medizinische Genetik und Humangenetik Charité Universitätsmedizin Berlin.

ChIP-seq. Expression Networks. Genomics: Lecture #14. Institut für Medizinische Genetik und Humangenetik Charité Universitätsmedizin Berlin. Expression Networks Institut für Medizinische Genetik und Humangenetik Charité Universitätsmedizin Berlin Genomics: Lecture #14 Gene Expression Networks The ultimate goal of experiments is to measure genome

More information

Regulation of Gene Expression

Regulation of Gene Expression Chapter 18 Regulation of Gene Expression PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley with contributions from

More information

Latent Variable models for GWAs

Latent Variable models for GWAs Latent Variable models for GWAs Oliver Stegle Machine Learning and Computational Biology Research Group Max-Planck-Institutes Tübingen, Germany September 2011 O. Stegle Latent variable models for GWAs

More information

Introduction to Bioinformatics

Introduction to Bioinformatics CSCI8980: Applied Machine Learning in Computational Biology Introduction to Bioinformatics Rui Kuang Department of Computer Science and Engineering University of Minnesota kuang@cs.umn.edu History of Bioinformatics

More information

Regulation of gene Expression in Prokaryotes & Eukaryotes

Regulation of gene Expression in Prokaryotes & Eukaryotes Regulation of gene Expression in Prokaryotes & Eukaryotes 1 The trp Operon Contains 5 genes coding for proteins (enzymes) required for the synthesis of the amino acid tryptophan. Also contains a promoter

More information

Topic 4: Equilibrium binding and chemical kinetics

Topic 4: Equilibrium binding and chemical kinetics Topic 4: Equilibrium binding and chemical kinetics Outline: Applications, applications, applications use Boltzmann to look at receptor-ligand binding use Boltzmann to look at PolII-DNA binding and gene

More information

Lecture 8 Learning Sequence Motif Models Using Expectation Maximization (EM) Colin Dewey February 14, 2008

Lecture 8 Learning Sequence Motif Models Using Expectation Maximization (EM) Colin Dewey February 14, 2008 Lecture 8 Learning Sequence Motif Models Using Expectation Maximization (EM) Colin Dewey February 14, 2008 1 Sequence Motifs what is a sequence motif? a sequence pattern of biological significance typically

More information

Regularization of the Inverse Laplace Transform with Applications in Nuclear Magnetic Resonance Relaxometry Candidacy Exam

Regularization of the Inverse Laplace Transform with Applications in Nuclear Magnetic Resonance Relaxometry Candidacy Exam Regularization of the Inverse Laplace Transform with Applications in Nuclear Magnetic Resonance Relaxometry Candidacy Exam Applied Mathematics, Applied Statistics, & Scientific Computation University of

More information

Genome 541 Introduction to Computational Molecular Biology. Max Libbrecht

Genome 541 Introduction to Computational Molecular Biology. Max Libbrecht Genome 541 Introduction to Computational Molecular Biology Max Libbrecht Genome 541 units Max Libbrecht: Gene regulation and epigenomics Postdoc, Bill Noble s lab Yi Yin: Bayesian statistics Postdoc, Jay

More information

When is MLE appropriate

When is MLE appropriate When is MLE appropriate As a rule of thumb the following to assumptions need to be fulfilled to make MLE the appropriate method for estimation: The model is adequate. That is, we trust that one of the

More information

Bioinformatics 2 - Lecture 4

Bioinformatics 2 - Lecture 4 Bioinformatics 2 - Lecture 4 Guido Sanguinetti School of Informatics University of Edinburgh February 14, 2011 Sequences Many data types are ordered, i.e. you can naturally say what is before and what

More information

Sig2GRN: A Software Tool Linking Signaling Pathway with Gene Regulatory Network for Dynamic Simulation

Sig2GRN: A Software Tool Linking Signaling Pathway with Gene Regulatory Network for Dynamic Simulation Sig2GRN: A Software Tool Linking Signaling Pathway with Gene Regulatory Network for Dynamic Simulation Authors: Fan Zhang, Runsheng Liu and Jie Zheng Presented by: Fan Wu School of Computer Science and

More information

Inferring Transcriptional Regulatory Networks from High-throughput Data

Inferring Transcriptional Regulatory Networks from High-throughput Data Inferring Transcriptional Regulatory Networks from High-throughput Data Lectures 9 Oct 26, 2011 CSE 527 Computational Biology, Fall 2011 Instructor: Su-In Lee TA: Christopher Miles Monday & Wednesday 12:00-1:20

More information

scrna-seq Differential expression analysis methods Olga Dethlefsen NBIS, National Bioinformatics Infrastructure Sweden October 2017

scrna-seq Differential expression analysis methods Olga Dethlefsen NBIS, National Bioinformatics Infrastructure Sweden October 2017 scrna-seq Differential expression analysis methods Olga Dethlefsen NBIS, National Bioinformatics Infrastructure Sweden October 2017 Olga (NBIS) scrna-seq de October 2017 1 / 34 Outline Introduction: what

More information

From transcriptome to epigenome, and back? Mattia Pelizzola Center for Genomic Science Istituto Italiano di Tecnologia (IIT), Milano

From transcriptome to epigenome, and back? Mattia Pelizzola Center for Genomic Science Istituto Italiano di Tecnologia (IIT), Milano From transcriptome to epigenome, and back? Mattia Pelizzola Center for Genomic Science Istituto Italiano di Tecnologia (IIT), Milano m6a epitranscriptome m6a dynamics synthesis processing degradation m5c

More information

Mixture models for analysing transcriptome and ChIP-chip data

Mixture models for analysing transcriptome and ChIP-chip data Mixture models for analysing transcriptome and ChIP-chip data Marie-Laure Martin-Magniette French National Institute for agricultural research (INRA) Unit of Applied Mathematics and Informatics at AgroParisTech,

More information

ChIP seq peak calling. Statistical integration between ChIP seq and RNA seq

ChIP seq peak calling. Statistical integration between ChIP seq and RNA seq Institute for Computational Biomedicine ChIP seq peak calling Statistical integration between ChIP seq and RNA seq Olivier Elemento, PhD ChIP-seq to map where transcription factors bind DNA Transcription

More information

RNA Synthesis and Processing

RNA Synthesis and Processing RNA Synthesis and Processing Introduction Regulation of gene expression allows cells to adapt to environmental changes and is responsible for the distinct activities of the differentiated cell types that

More information

Eukaryotic Gene Expression

Eukaryotic Gene Expression Eukaryotic Gene Expression Lectures 22-23 Several Features Distinguish Eukaryotic Processes From Mechanisms in Bacteria 123 Eukaryotic Gene Expression Several Features Distinguish Eukaryotic Processes

More information

GAUSSIAN PROCESS REGRESSION

GAUSSIAN PROCESS REGRESSION GAUSSIAN PROCESS REGRESSION CSE 515T Spring 2015 1. BACKGROUND The kernel trick again... The Kernel Trick Consider again the linear regression model: y(x) = φ(x) w + ε, with prior p(w) = N (w; 0, Σ). The

More information

Expression QTLs and Mapping of Complex Trait Loci. Paul Schliekelman Statistics Department University of Georgia

Expression QTLs and Mapping of Complex Trait Loci. Paul Schliekelman Statistics Department University of Georgia Expression QTLs and Mapping of Complex Trait Loci Paul Schliekelman Statistics Department University of Georgia Definitions: Genes, Loci and Alleles A gene codes for a protein. Proteins due everything.

More information

Networks & pathways. Hedi Peterson MTAT Bioinformatics

Networks & pathways. Hedi Peterson MTAT Bioinformatics Networks & pathways Hedi Peterson (peterson@quretec.com) MTAT.03.239 Bioinformatics 03.11.2010 Networks are graphs Nodes Edges Edges Directed, undirected, weighted Nodes Genes Proteins Metabolites Enzymes

More information

L3.1: Circuits: Introduction to Transcription Networks. Cellular Design Principles Prof. Jenna Rickus

L3.1: Circuits: Introduction to Transcription Networks. Cellular Design Principles Prof. Jenna Rickus L3.1: Circuits: Introduction to Transcription Networks Cellular Design Principles Prof. Jenna Rickus In this lecture Cognitive problem of the Cell Introduce transcription networks Key processing network

More information

Introduction. Gene expression is the combined process of :

Introduction. Gene expression is the combined process of : 1 To know and explain: Regulation of Bacterial Gene Expression Constitutive ( house keeping) vs. Controllable genes OPERON structure and its role in gene regulation Regulation of Eukaryotic Gene Expression

More information

Domains of a- and b-globin Genes in the Context of the Structural Functional Organization of the Eukaryotic Genome

Domains of a- and b-globin Genes in the Context of the Structural Functional Organization of the Eukaryotic Genome ISSN 0006-2979, Biochemistry (Moscow), 2012, Vol. 77, No. 13, pp. 1409-1423. Pleiades Publishing, Ltd., 2012. Original Russian Text S. V. Razin, S. V. Ulianov, E. S. Ioudinkova, E. S. Gushchanskaya, A.

More information

Biostatistics-Lecture 16 Model Selection. Ruibin Xi Peking University School of Mathematical Sciences

Biostatistics-Lecture 16 Model Selection. Ruibin Xi Peking University School of Mathematical Sciences Biostatistics-Lecture 16 Model Selection Ruibin Xi Peking University School of Mathematical Sciences Motivating example1 Interested in factors related to the life expectancy (50 US states,1969-71 ) Per

More information

Prokaryotic Regulation

Prokaryotic Regulation Prokaryotic Regulation Control of transcription initiation can be: Positive control increases transcription when activators bind DNA Negative control reduces transcription when repressors bind to DNA regulatory

More information

Whole-genome analysis of GCN4 binding in S.cerevisiae

Whole-genome analysis of GCN4 binding in S.cerevisiae Whole-genome analysis of GCN4 binding in S.cerevisiae Lillian Dai Alex Mallet Gcn4/DNA diagram (CREB symmetric site and AP-1 asymmetric site: Song Tan, 1999) removed for copyright reasons. What is GCN4?

More information

Co-ordination occurs in multiple layers Intracellular regulation: self-regulation Intercellular regulation: coordinated cell signalling e.g.

Co-ordination occurs in multiple layers Intracellular regulation: self-regulation Intercellular regulation: coordinated cell signalling e.g. Gene Expression- Overview Differentiating cells Achieved through changes in gene expression All cells contain the same whole genome A typical differentiated cell only expresses ~50% of its total gene Overview

More information

Chapter 18 Lecture. Concepts of Genetics. Tenth Edition. Developmental Genetics

Chapter 18 Lecture. Concepts of Genetics. Tenth Edition. Developmental Genetics Chapter 18 Lecture Concepts of Genetics Tenth Edition Developmental Genetics Chapter Contents 18.1 Differentiated States Develop from Coordinated Programs of Gene Expression 18.2 Evolutionary Conservation

More information

MODEL-BASED APPROACHES FOR THE DETECTION OF BIOLOGICALLY ACTIVE GENOMIC REGIONS FROM NEXT GENERATION SEQUENCING DATA. Naim Rashid

MODEL-BASED APPROACHES FOR THE DETECTION OF BIOLOGICALLY ACTIVE GENOMIC REGIONS FROM NEXT GENERATION SEQUENCING DATA. Naim Rashid MODEL-BASED APPROACHES FOR THE DETECTION OF BIOLOGICALLY ACTIVE GENOMIC REGIONS FROM NEXT GENERATION SEQUENCING DATA Naim Rashid A dissertation submitted to the faculty of the University of North Carolina

More information

Genetic Networks. Korbinian Strimmer. Seminar: Statistical Analysis of RNA-Seq Data 19 June IMISE, Universität Leipzig

Genetic Networks. Korbinian Strimmer. Seminar: Statistical Analysis of RNA-Seq Data 19 June IMISE, Universität Leipzig Genetic Networks Korbinian Strimmer IMISE, Universität Leipzig Seminar: Statistical Analysis of RNA-Seq Data 19 June 2012 Korbinian Strimmer, RNA-Seq Networks, 19/6/2012 1 Paper G. I. Allen and Z. Liu.

More information

GENE REGULATION AND PROBLEMS OF DEVELOPMENT

GENE REGULATION AND PROBLEMS OF DEVELOPMENT GENE REGULATION AND PROBLEMS OF DEVELOPMENT By Surinder Kaur DIET Ropar Surinder_1998@ yahoo.in Mob No 9988530775 GENE REGULATION Gene is a segment of DNA that codes for a unit of function (polypeptide,

More information

CSEP 590A Summer Lecture 4 MLE, EM, RE, Expression

CSEP 590A Summer Lecture 4 MLE, EM, RE, Expression CSEP 590A Summer 2006 Lecture 4 MLE, EM, RE, Expression 1 FYI, re HW #2: Hemoglobin History Alberts et al., 3rd ed.,pg389 2 Tonight MLE: Maximum Likelihood Estimators EM: the Expectation Maximization Algorithm

More information

CSEP 590A Summer Tonight MLE. FYI, re HW #2: Hemoglobin History. Lecture 4 MLE, EM, RE, Expression. Maximum Likelihood Estimators

CSEP 590A Summer Tonight MLE. FYI, re HW #2: Hemoglobin History. Lecture 4 MLE, EM, RE, Expression. Maximum Likelihood Estimators CSEP 59A Summer 26 Lecture 4 MLE, EM, RE, Expression FYI, re HW #2: Hemoglobin History 1 Alberts et al., 3rd ed.,pg389 2 Tonight MLE: Maximum Likelihood Estimators EM: the Expectation Maximization Algorithm

More information

Chapter 20. Initiation of transcription. Eukaryotic transcription initiation

Chapter 20. Initiation of transcription. Eukaryotic transcription initiation Chapter 20. Initiation of transcription Eukaryotic transcription initiation 2003. 5.22 Prokaryotic vs eukaryotic Bacteria = one RNA polymerase Eukaryotes have three RNA polymerases (I, II, and III) in

More information

10-810: Advanced Algorithms and Models for Computational Biology. microrna and Whole Genome Comparison

10-810: Advanced Algorithms and Models for Computational Biology. microrna and Whole Genome Comparison 10-810: Advanced Algorithms and Models for Computational Biology microrna and Whole Genome Comparison Central Dogma: 90s Transcription factors DNA transcription mrna translation Proteins Central Dogma:

More information

Honors Biology Reading Guide Chapter 11

Honors Biology Reading Guide Chapter 11 Honors Biology Reading Guide Chapter 11 v Promoter a specific nucleotide sequence in DNA located near the start of a gene that is the binding site for RNA polymerase and the place where transcription begins

More information

Notation. Pattern Recognition II. Michal Haindl. Outline - PR Basic Concepts. Pattern Recognition Notions

Notation. Pattern Recognition II. Michal Haindl. Outline - PR Basic Concepts. Pattern Recognition Notions Notation S pattern space X feature vector X = [x 1,...,x l ] l = dim{x} number of features X feature space K number of classes ω i class indicator Ω = {ω 1,...,ω K } g(x) discriminant function H decision

More information

Carri-Lyn Mead Thursday, January 13, 2005 Terry Fox Laboratory, Dr. Dixie Mager

Carri-Lyn Mead Thursday, January 13, 2005 Terry Fox Laboratory, Dr. Dixie Mager Investigating Trends in Transposable Element Insertion within Regulatory Regions Carri-Lyn Mead cmead@bcgsc.ca Thursday, January 13, 2005 Terry Fox Laboratory, Dr. Dixie Mager Outline Transposable Element

More information

Multiple Choice Review- Eukaryotic Gene Expression

Multiple Choice Review- Eukaryotic Gene Expression Multiple Choice Review- Eukaryotic Gene Expression 1. Which of the following is the Central Dogma of cell biology? a. DNA Nucleic Acid Protein Amino Acid b. Prokaryote Bacteria - Eukaryote c. Atom Molecule

More information

Slide 1 / 7. Free Response

Slide 1 / 7. Free Response Slide 1 / 7 Free Response Slide 2 / 7 Slide 3 / 7 1 The above diagrams illustrate the experiments carried out by Griffith and Hershey and Chaserespectively. Describe the hypothesis or conclusion that each

More information

Bi 8 Lecture 11. Quantitative aspects of transcription factor binding and gene regulatory circuit design. Ellen Rothenberg 9 February 2016

Bi 8 Lecture 11. Quantitative aspects of transcription factor binding and gene regulatory circuit design. Ellen Rothenberg 9 February 2016 Bi 8 Lecture 11 Quantitative aspects of transcription factor binding and gene regulatory circuit design Ellen Rothenberg 9 February 2016 Major take-home messages from λ phage system that apply to many

More information

CS-E5880 Modeling biological networks Gene regulatory networks

CS-E5880 Modeling biological networks Gene regulatory networks CS-E5880 Modeling biological networks Gene regulatory networks Jukka Intosalmi (based on slides by Harri Lähdesmäki) Department of Computer Science Aalto University January 12, 2018 Outline Modeling gene

More information

An Introduction to the spls Package, Version 1.0

An Introduction to the spls Package, Version 1.0 An Introduction to the spls Package, Version 1.0 Dongjun Chung 1, Hyonho Chun 1 and Sündüz Keleş 1,2 1 Department of Statistics, University of Wisconsin Madison, WI 53706. 2 Department of Biostatistics

More information

Big Idea 3: Living systems store, retrieve, transmit and respond to information essential to life processes. Tuesday, December 27, 16

Big Idea 3: Living systems store, retrieve, transmit and respond to information essential to life processes. Tuesday, December 27, 16 Big Idea 3: Living systems store, retrieve, transmit and respond to information essential to life processes. Enduring understanding 3.B: Expression of genetic information involves cellular and molecular

More information

Developmental Biology Lecture Outlines

Developmental Biology Lecture Outlines Developmental Biology Lecture Outlines Lecture 01: Introduction Course content Developmental Biology Obsolete hypotheses Current theory Lecture 02: Gametogenesis Spermatozoa Spermatozoon function Spermatozoon

More information

Topic 4 - #14 The Lactose Operon

Topic 4 - #14 The Lactose Operon Topic 4 - #14 The Lactose Operon The Lactose Operon The lactose operon is an operon which is responsible for the transport and metabolism of the sugar lactose in E. coli. - Lactose is one of many organic

More information

Statistics for high-dimensional data: Group Lasso and additive models

Statistics for high-dimensional data: Group Lasso and additive models Statistics for high-dimensional data: Group Lasso and additive models Peter Bühlmann and Sara van de Geer Seminar für Statistik, ETH Zürich May 2012 The Group Lasso (Yuan & Lin, 2006) high-dimensional

More information

Regularization in Cox Frailty Models

Regularization in Cox Frailty Models Regularization in Cox Frailty Models Andreas Groll 1, Trevor Hastie 2, Gerhard Tutz 3 1 Ludwig-Maximilians-Universität Munich, Department of Mathematics, Theresienstraße 39, 80333 Munich, Germany 2 University

More information

Ultra High Dimensional Variable Selection with Endogenous Variables

Ultra High Dimensional Variable Selection with Endogenous Variables 1 / 39 Ultra High Dimensional Variable Selection with Endogenous Variables Yuan Liao Princeton University Joint work with Jianqing Fan Job Market Talk January, 2012 2 / 39 Outline 1 Examples of Ultra High

More information

Second-Order Analysis of Spatial Point Processes

Second-Order Analysis of Spatial Point Processes Title Second-Order Analysis of Spatial Point Process Tonglin Zhang Outline Outline Spatial Point Processes Intensity Functions Mean and Variance Pair Correlation Functions Stationarity K-functions Some

More information

Cell Division: Mitosis and Meiosis

Cell Division: Mitosis and Meiosis Cell Division: Mitosis and Meiosis Omnis cellula e cellula Cell Div. In Eukaryotic cells Chromosome structure Biparental inheritance and diploidy Mitosis: Karyokinesis and cytokinesis 5 stages: PPMAT Regulating

More information

BIS &003 Answers to Assigned Problems May 23, Week /18.6 How would you distinguish between an enhancer and a promoter?

BIS &003 Answers to Assigned Problems May 23, Week /18.6 How would you distinguish between an enhancer and a promoter? Week 9 Study Questions from the textbook: 6 th Edition: Chapter 19-19.6, 19.7, 19.15, 19.17 OR 7 th Edition: Chapter 18-18.6 18.7, 18.15, 18.17 19.6/18.6 How would you distinguish between an enhancer and

More information

Linear regression. Example 3.6.1: Relation between abrasion loss and hardness of rubber tires.

Linear regression. Example 3.6.1: Relation between abrasion loss and hardness of rubber tires. Linear regression Example 3.6.1: Relation between abrasion loss and hardness of rubber tires. X i the abrasion loss 30 observations, i = 1,..., 30. Sample space R 30. y i the hardness, i = 1,..., 30. Model:

More information

16 CONTROL OF GENE EXPRESSION

16 CONTROL OF GENE EXPRESSION 16 CONTROL OF GENE EXPRESSION Chapter Outline 16.1 REGULATION OF GENE EXPRESSION IN PROKARYOTES The operon is the unit of transcription in prokaryotes The lac operon for lactose metabolism is transcribed

More information

Regulation of Gene Expression

Regulation of Gene Expression Chapter 18 Regulation of Gene Expression Edited by Shawn Lester PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley

More information

A Bayesian Perspective on Residential Demand Response Using Smart Meter Data

A Bayesian Perspective on Residential Demand Response Using Smart Meter Data A Bayesian Perspective on Residential Demand Response Using Smart Meter Data Datong-Paul Zhou, Maximilian Balandat, and Claire Tomlin University of California, Berkeley [datong.zhou, balandat, tomlin]@eecs.berkeley.edu

More information

Gene regulation: From biophysics to evolutionary genetics

Gene regulation: From biophysics to evolutionary genetics Gene regulation: From biophysics to evolutionary genetics Michael Lässig Institute for Theoretical Physics University of Cologne Thanks Ville Mustonen Johannes Berg Stana Willmann Curt Callan (Princeton)

More information

Missing Data Issues in the Studies of Neurodegenerative Disorders: the Methodology

Missing Data Issues in the Studies of Neurodegenerative Disorders: the Methodology Missing Data Issues in the Studies of Neurodegenerative Disorders: the Methodology Sheng Luo, PhD Associate Professor Department of Biostatistics & Bioinformatics Duke University Medical Center sheng.luo@duke.edu

More information

Introduction to Gaussian Processes

Introduction to Gaussian Processes Introduction to Gaussian Processes Neil D. Lawrence GPSS 10th June 2013 Book Rasmussen and Williams (2006) Outline The Gaussian Density Covariance from Basis Functions Basis Function Representations Constructing

More information

Machine learning - HT Maximum Likelihood

Machine learning - HT Maximum Likelihood Machine learning - HT 2016 3. Maximum Likelihood Varun Kanade University of Oxford January 27, 2016 Outline Probabilistic Framework Formulate linear regression in the language of probability Introduce

More information

Unlocking RNA-seq tools for zero inflation and single cell applications using observation weights

Unlocking RNA-seq tools for zero inflation and single cell applications using observation weights Unlocking RNA-seq tools for zero inflation and single cell applications using observation weights Koen Van den Berge, Ghent University Statistical Genomics, 2018-2019 1 The team Koen Van den Berge Fanny

More information

Models of transcriptional regulation

Models of transcriptional regulation Models of transcriptional regulation We have already discussed four simple mechanisms of transcriptional regulation, nuclear exclusion nuclear concentration modification of bound activator redirection

More information

Lecture 2: Read about the yeast MAT locus in Molecular Biology of the Gene. Watson et al. Chapter 10. Plus section on yeast as a model system Read

Lecture 2: Read about the yeast MAT locus in Molecular Biology of the Gene. Watson et al. Chapter 10. Plus section on yeast as a model system Read Lecture 2: Read about the yeast MAT locus in Molecular Biology of the Gene. Watson et al. Chapter 10. Plus section on yeast as a model system Read chapter 22 and chapter 10 [section on MATing type gene

More information

The Eukaryotic Genome and Its Expression. The Eukaryotic Genome and Its Expression. A. The Eukaryotic Genome. Lecture Series 11

The Eukaryotic Genome and Its Expression. The Eukaryotic Genome and Its Expression. A. The Eukaryotic Genome. Lecture Series 11 The Eukaryotic Genome and Its Expression Lecture Series 11 The Eukaryotic Genome and Its Expression A. The Eukaryotic Genome B. Repetitive Sequences (rem: teleomeres) C. The Structures of Protein-Coding

More information

Topographic Mapping and Dimensionality Reduction of Binary Tensor Data of Arbitrary Rank

Topographic Mapping and Dimensionality Reduction of Binary Tensor Data of Arbitrary Rank Topographic Mapping and Dimensionality Reduction of Binary Tensor Data of Arbitrary Rank Peter Tiňo, Jakub Mažgút, Hong Yan, Mikael Bodén Topographic Mapping and Dimensionality Reduction of Binary Tensor

More information

Warm-Up. Explain how a secondary messenger is activated, and how this affects gene expression. (LO 3.22)

Warm-Up. Explain how a secondary messenger is activated, and how this affects gene expression. (LO 3.22) Warm-Up Explain how a secondary messenger is activated, and how this affects gene expression. (LO 3.22) Yesterday s Picture The first cell on Earth (approx. 3.5 billion years ago) was simple and prokaryotic,

More information

Bayes methods for categorical data. April 25, 2017

Bayes methods for categorical data. April 25, 2017 Bayes methods for categorical data April 25, 2017 Motivation for joint probability models Increasing interest in high-dimensional data in broad applications Focus may be on prediction, variable selection,

More information