title: Learning "graph-mer" motifs that predict gene expression trajectories in development.
abstract: A key problem in understanding transcriptional regulatory networks is deciphering what cis regulatory logic is encoded in gene promoter sequences and how this sequence information maps to expression. A typical computational approach to this problem involves clustering genes by their expression profiles and then searching for overrepresented motifs in the promoter sequences of genes in a cluster. However, genes with similar expression profiles may be controlled by distinct regulatory programs. Moreover, if many gene expression profiles in a data set are highly correlated, as in the case of whole organism developmental time series, it may be difficult to resolve fine-grained clusters in the first place. We present a predictive framework for modeling the natural flow of information, from promoter sequence to expression, to learn cis regulatory motifs and characterize gene expression patterns in developmental time courses. We introduce a cluster-free algorithm based on a graph-regularized version of partial least squares (PLS) regression to learn sequence patterns--represented by graphs of k-mers, or "graph-mers"--that predict gene expression trajectories. Applying the approach to wildtype germline development in Caenorhabditis elegans, we found that the first and second latent PLS factors mapped to expression profiles for oocyte and sperm genes, respectively. We extracted both known and novel motifs from the graph-mers associated to these germline-specific patterns, including novel CG-rich motifs specific to oocyte genes. We found evidence supporting the functional relevance of these putative regulatory elements through analysis of positional bias, motif conservation and in situ gene expression. This study demonstrates that our regression model can learn biologically meaningful latent structure and identify potentially functional motifs from subtle developmental time course expression data.
authors.affiliation: Department of Physics, Columbia University, New York, New York, United States of America.
authors.forename: Xuejing
authors.initials: X
authors.lastname: Li
authors.forename: Casandra
authors.initials: C
authors.lastname: Panea
authors.forename: Chris H
authors.initials: CH
authors.lastname: Wiggins
authors.forename: Valerie
authors.initials: V
authors.lastname: Reinke
authors.forename: Christina
authors.initials: C
authors.lastname: Leslie
mesh_terms: D000465:Algorithms; D000818:Animals; D017173:Caenorhabditis elegans / Q000235:genetics; D020869:Gene Expression Profiling / Q000379:methods*; D016018:Least-Squares Analysis; D008297:Male; D005810:Multigene Family; D015999:Multivariate Analysis; D009865:Oocytes; D025341:Principal Component Analysis; D011401:Promoter Regions, Genetic; D012044:Regression Analysis; D015203:Reproducibility of Results; D013094:Spermatozoa
references.citation: Cell. 2004 Apr 16;117(2):185-98
references.pmid: 15084257
references.citation: Mol Cells. 1999 Oct 31;9(5):535-41
references.pmid: 10597043
references.citation: Brief Bioinform. 2007 Jan;8(1):32-44
references.pmid: 16772269
references.citation: PLoS Comput Biol. 2009 Jan;5(1):e1000269
references.pmid: 19180174
references.citation: Pac Symp Biocomput. 2003;:29-40
references.pmid: 12603015
references.citation: Proc Int Conf Intell Syst Mol Biol. 1994;2:28-36
references.pmid: 7584402
references.citation: Science. 2000 Oct 27;290(5492):809-12
references.pmid: 11052945
references.citation: BMC Bioinformatics. 2003 Jan 13;4:2
references.pmid: 12525261
references.citation: Genome Biol. 2006;7(5):R36
references.pmid: 16686963
references.citation: Nature. 2008 Jan 31;451(7178):535-40
references.pmid: 18172436
references.citation: Nature. 2002 Dec 5;420(6915):520-62
references.pmid: 12466850
references.citation: Bioinformatics. 2007 Jun 15;23(12):1486-94
references.pmid: 17463025
references.citation: Development. 2004 Jan;131(2):311-23
references.pmid: 14668411
references.citation: Nat Genet. 2003 Jun;34(2):166-76
references.pmid: 12740579
references.citation: Mol Endocrinol. 1999 May;13(5):774-86
references.pmid: 10319327
references.citation: Mol Syst Biol. 2007;3:74
references.pmid: 17224918
references.citation: Nat Genet. 2001 Feb;27(2):167-71
references.pmid: 11175784
references.citation: Nat Biotechnol. 2005 Jan;23(1):137-44
references.pmid: 15637633
references.citation: Genome Res. 2003 Nov;13(11):2498-504
references.pmid: 14597658
references.citation: BMC Bioinformatics. 2007 Feb 01;8:35
references.pmid: 17270037
references.citation: Nucleic Acids Res. 2006;34(20):5730-9
references.pmid: 17041233
references.citation: PLoS Comput Biol. 2008 Nov;4(11):e1000224
references.pmid: 19008939
references.citation: Pac Symp Biocomput. 2000;:455-66
references.pmid: 10902193
references.citation: Bioinformatics. 2007 Feb 15;23(4):493-5
references.pmid: 17138590