Learning "graph-mer" motifs that predict gene expression trajectories in development.
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.
Li
Xuejing
X
Department of Physics, Columbia University, New York, New York, United States of America.
Panea
Casandra
C
Wiggins
Chris H
CH
Reinke
Valerie
V
Leslie
Christina
C
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
Cell. 2004 Apr 16;117(2):185-98
15084257
Mol Cells. 1999 Oct 31;9(5):535-41
10597043
Brief Bioinform. 2007 Jan;8(1):32-44
16772269
PLoS Comput Biol. 2009 Jan;5(1):e1000269
19180174
Pac Symp Biocomput. 2003;:29-40
12603015
Proc Int Conf Intell Syst Mol Biol. 1994;2:28-36
7584402
Science. 2000 Oct 27;290(5492):809-12
11052945
BMC Bioinformatics. 2003 Jan 13;4:2
12525261
Genome Biol. 2006;7(5):R36
16686963
Nature. 2008 Jan 31;451(7178):535-40
18172436
Nature. 2002 Dec 5;420(6915):520-62
12466850
Bioinformatics. 2007 Jun 15;23(12):1486-94
17463025
Development. 2004 Jan;131(2):311-23
14668411
Nat Genet. 2003 Jun;34(2):166-76
12740579
Mol Endocrinol. 1999 May;13(5):774-86
10319327
Mol Syst Biol. 2007;3:74
17224918
Nat Genet. 2001 Feb;27(2):167-71
11175784
Nat Biotechnol. 2005 Jan;23(1):137-44
15637633
Genome Res. 2003 Nov;13(11):2498-504
14597658
BMC Bioinformatics. 2007 Feb 01;8:35
17270037
Nucleic Acids Res. 2006;34(20):5730-9
17041233
PLoS Comput Biol. 2008 Nov;4(11):e1000224
19008939
Pac Symp Biocomput. 2000;:455-66
10902193
Bioinformatics. 2007 Feb 15;23(4):493-5
17138590