Nucleic Acids Research Advance Access published online on April 11, 2007
Nucleic Acids Research, doi:10.1093/nar/gkm130
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Computational Biology |
Error-pooling-based statistical methods for identifying novel temporal replication profiles of human chromosomes observed by DNA tiling arrays
1Department of Statistics, Seoul National University, Korea, 2Department of Biochemistry and Molecular Genetics, University of Virginia School of Medicine, USA and 3Division of Biostatistics and Epidemiology, University of Virginia School of Medicine, USA
*To whom correspondence should be addressed. Tel: +1-434 982-1033; Fax: +1-434 924-8437; Email: jaeklee{at}virginia.edu
Received February 6, 2007. Revised March 20, 2007. Accepted February 16, 2007.
Statistical analysis on tiling array data is extremely challenging due to the astronomically large number of sequence probes, high noise levels of individual probes and limited number of replicates in these data. To overcome these difficulties, we first developed statistical error estimation and weighted ANOVA modeling approaches to high-density tiling array data, especially the former based on an advanced error-pooling method to accurately obtain heterogeneous technical error of small-sample tiling array data. Based on these approaches, we analyzed the high-density tiling array data of the temporal replication patterns during cell-cycle S phase of synchronized HeLa cells on human chromosomes 21 and 22. We found many novel temporal replication patterns, identifying about 26% of over 1 million tiling array sequence probes with significant differential replication during the four 2-h time periods of S phase. Among these differentially replicated probes, 126 941 sequence probes were matched to 417 known genes. The majority of these genes were found to be replicated within one or two consecutive time periods, while the others were replicated at two non-consecutive time periods. Also, coding regions found to be more differentially replicated in particular time periods than noncoding regions in the gene-poor chromosome 21 (25% differentially replicated among genic probes versus 18.6% among intergenic probes), while such a phenomenon was less prominent in gene-rich chromosome 22. A rigorous statistical testing for local proximity of differentially replicated genic and intergenic probes was performed to identify significant stretches of differentially replicated sequence regions. From this analysis, we found that adjacent genes were frequently replicated at different time periods, potentially implying the existence of quite dense replication origins. Evaluating the conditional probability significance of identified gene ontology terms on chromosomes 21 and 22, we detected some over-represented molecular functions and biological processes among these differentially replicated genes, such as the ones relevant to hydrolase, transferase and receptor-binding activities. Some of these results were confirmed showing >70% consistency with cDNA microarray data that were independently generated in parallel with the tiling arrays. Thus, our improved analysis approaches specifically designed for high-density tiling array data enabled us to reliably and sensitively identify many novel temporal replication patterns on human chromosomes.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  Y. Kim, S. Bekiranov, J. K. Lee, and T. Park Double error shrinkage method for identifying protein binding sites observed by tiling arrays with limited replication Bioinformatics, October 1, 2009; 25(19): 2486 - 2491. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C. Murie and R. Nadon A correction for estimating error when using the Local Pooled Error Statistical Test Bioinformatics, August 1, 2008; 24(15): 1735 - 1736. [Abstract] [Full Text] [PDF]
|
|