On the choice and number of microarrays for transcriptional regulatory network inference.
BACKGROUND
Transcriptional regulatory network inference (TRNI) from large compendia of DNA microarrays has become a fundamental approach for discovering transcription factor (TF)-gene interactions at the genome-wide level. In correlation-based TRNI, network edges can in principle be evaluated using standard statistical tests. However, while such tests nominally assume independent microarray experiments, we expect dependency between the experiments in microarray compendia, due to both project-specific factors (e.g., microarray preparation, environmental effects) in the multi-project compendium setting and effective dependency induced by gene-gene correlations. Herein, we characterize the nature of dependency in an Escherichia coli microarray compendium and explore its consequences on the problem of determining which and how many arrays to use in correlation-based TRNI.


RESULTS
We present evidence of substantial effective dependency among microarrays in this compendium, and characterize that dependency with respect to experimental condition factors. We then introduce a measure neff of the effective number of experiments in a compendium, and find that corresponding to the dependency observed in this particular compendium there is a huge reduction in effective sample size i.e., neff = 14.7 versus n = 376. Furthermore, we found that the neff of select subsets of experiments actually exceeded neff of the full compendium, suggesting that the adage 'less is more' applies here. Consistent with this latter result, we observed improved performance in TRNI using subsets of the data compared to results using the full compendium. We identified experimental condition factors that trend with changes in TRNI performance and neff , including growth phase and media type. Finally, using the set of known E. coli genetic regulatory interactions from RegulonDB, we demonstrated that false discovery rates (FDR) derived from neff -adjusted p-values were well-matched to FDR based on the RegulonDB truth set.


CONCLUSIONS
These results support utilization of neff as a potent descriptor of microarray compendia. In addition, they highlight a straightforward correlation-based method for TRNI with demonstrated meaningful statistical testing for significant edges, readily applicable to compendia from any species, even when a truth set is not available. This work facilitates a more refined approach to construction and utilization of mRNA expression compendia in TRNI.
Cosgrove
Elissa J
EJ
Department of Mathematics and Statistics, Boston University, Boston, MA, USA.
Gardner
Timothy S
TS
Kolaczyk
Eric D
ED
D004926:Escherichia coli / Q000378:metabolism; D005786:Gene Expression Regulation; D053263:Gene Regulatory Networks / Q000235:genetics*; D020411:Oligonucleotide Array Sequence Analysis; D012333:RNA, Messenger / Q000378:metabolism; D050436:Regulatory Elements, Transcriptional; D018401:Sample Size; D014157:Transcription Factors / Q000235:genetics / Q000378:metabolism
D012333:RNA, Messenger; D014157:Transcription Factors
Pac Symp Biocomput. 2000;:418-29
10902190
Nucleic Acids Res. 2008 Jan;36(Database issue):D866-70
17932051
J Biol Chem. 2002 Oct 4;277(40):37001-8
12121991
Bioinformatics. 2005 Mar;21(6):754-64
15479708
BMC Syst Biol. 2009 Apr 14;3:39
19366454
Ann Appl Stat. 2009 Jan 1;3(3):922-942
20563291
Bioinformatics. 2004 Jun 12;20(9):1436-46
14962916
Nucleic Acids Res. 2008 Jan;36(Database issue):D120-4
18158297
Genome Biol. 2006;7(4):R32
16626502
PLoS Biol. 2007 Jan;5(1):e8
17214507
Nature. 2007 Oct 18;449(7164):804-10
17943116
Nucleic Acids Res. 2009 Jan;37(Database issue):D885-90
18940857
Stat Appl Genet Mol Biol. 2005;4:Article34
16646853
BMC Bioinformatics. 2008 Jul 09;9:303
18613966
Proc AMIA Symp. 1999;:711-5
10566452
Bioinformatics. 2004 Dec 12;20(18):3565-74
15284096
Nucleic Acids Res. 2003 Feb 15;31(4):e15
12582260