Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter September 21, 2016

A simulation framework for correlated count data of features subsets in high-throughput sequencing or proteomics experiments

  • Jochen Kruppa , Frank Kramer , Tim Beißbarth and Klaus Jung EMAIL logo


As part of the data processing of high-throughput-sequencing experiments count data are produced representing the amount of reads that map to specific genomic regions. Count data also arise in mass spectrometric experiments for the detection of protein-protein interactions. For evaluating new computational methods for the analysis of sequencing count data or spectral count data from proteomics experiments artificial count data is thus required. Although, some methods for the generation of artificial sequencing count data have been proposed, all of them simulate single sequencing runs, omitting thus the correlation structure between the individual genomic features, or they are limited to specific structures. We propose to draw correlated data from the multivariate normal distribution and round these continuous data in order to obtain discrete counts. In our approach, the required distribution parameters can either be constructed in different ways or estimated from real count data. Because rounding affects the correlation structure we evaluate the use of shrinkage estimators that have already been used in the context of artificial expression data from DNA microarrays. Our approach turned out to be useful for the simulation of counts for defined subsets of features such as individual pathways or GO categories.


The work was partially supported by the German Federal Ministry of Education and Research (BMBF) through the projects ELSA Genoperspektiv (grant number 01GP1402) and e:Med MMML-Demonstrators (grant number 031A428B). The authors would also like to thank the reviewers for further references.


Adler, A. S., M. L. McCleland, S. Yee, M. Yaylaoglu, S. Hussain, E. Cosino, E. Quinones, Z. Modrusan, S. Seshagiri, E. Torres, V. S. Chopra, B. Haley, Z. Zhang, E. M. Blackwood, M. Singh, M. Junttila, J-P. Stephan, J. Liu, G. Pau, E. R. Fearon, Z. Jiang and R. Firestein (2014): “An integrative analysis of colon cancer identifies an essential function for PRPF6 in tumor growth,” Genes. Dev., 28, 1068–1084.10.1101/gad.237206.113Search in Google Scholar PubMed PubMed Central

Allen, G. I. and Z. Liu (2012): “A log-linear graphical model for inferring genetic networks from high-throughput sequencing data,” IEEE Int. Conf. Bioinf. Biomed., 41–46. doi: 10.1109/BIBM.2012.6392619.Search in Google Scholar

Allen, G. I. and Z. Liu (2013): “A local poisson graphical model for inferring genetic networks from next generation sequencing data,” IEEE Trans. Nanobiosci., 12, 1–10.Search in Google Scholar

Anders, S. and W. Huber (2010): “Differential expression analysis for sequence count data,” Genome Biol., 11, R106.10.1186/gb-2010-11-10-r106Search in Google Scholar PubMed PubMed Central

Anders, S., P. T. Pyl and W. Huber (2015): “HTSeq – A Python framework to work with high-throughput sequencing data,” Bioinformatics, 31, 166–169.10.1093/bioinformatics/btu638Search in Google Scholar PubMed PubMed Central

Böhning, D., E. Dietz and P. Schlattmann (1999): “The zero-inflated poisson model and the decayed, missing and filled teeth index in dental epidemiology,” J. Royal. Stat. Soc., Series A, 162, 195–209.10.1111/1467-985X.00130Search in Google Scholar

Canale, A. and D. B. Dunson (2012): “Nonparametric Bayes modelling of count processes,” Biometrika, 100, 801–816.10.1093/biomet/ast037Search in Google Scholar

Choi, H., D. Fermin and A. I. Nesvizhskii (2008): “Significance analysis of spectral count data in label-free shotgun proteomics,” Mol. Cell. Proteomics, 7, 2373–2385.10.1074/mcp.M800203-MCP200Search in Google Scholar PubMed PubMed Central

Demir, E., M. P. Cary, S. Paley, K. Fukuda, C. Lemer, I. Vastrik, G. Wu, P. D’Eustachio, C. Schaefer, J. Luciano, F. Schacherer, I. Martinez-Flores, Z. Hu, V. Jimenez-Jacinto, G. Joshi-Tope, K. Kandasamy, A. C. Lopez-Fuentes, H. Mi, E. Pichler, I. Rodchenkov, A. Splendiani, S. Tkachev, J. Zucker, G. Gopinath, H. Rajasimha, R. Ramakrishnan, I. Shah, M. Syed, N. Anwar, O. Babur, M. Blinov, E. Brauner, D. Corwin, S. Donaldson, F. Gibbons, R. Goldberg, P. Hornbeck, A. Luna, P. Murray-Rust, E. Neumann, O. Ruebenacker, M. Samwald, M. van Iersel, S. Wimalaratne, K. Allen, B. Braun, M. Whirl-Carrillo, K. H. Cheung, K. Dahlquist, A. Finney, M. Gillespie, E. Glass, L. Gong, R. Haw, M. Honig, O. Hubaut, D. Kane, S. Krupa, M. Kutmon, J. Leonard, D. Marks, D. Merberg, V. Petri, A. Pico, D. Ravenscroft, L. Ren, N. Shah, M. Sunshine, R. Tang, R. Whaley, S. Letovksy, K. H. Buetow, A. Rzhetsky, V. Schachter, B. S. Sobral, U. Dogrusoz, S. McWeeney, M. Aladjem, E. Birney, J. Collado-Vides, S. Goto, M. Hucka, N. Le Novère, N. Maltsev, A. Pandey, P. Thomas, E. Wingender, P. D. Karp, C. Sander and G. D. Bader (2010): “The BioPAX community standard for pathway data sharing,” Nat. Biotechnol., 28, 935–942.10.1038/nbt.1666Search in Google Scholar PubMed PubMed Central

Fischer, M., S. Zilkenat, R. G. Gerlach, S. Wagner and B. Y. Renard (2014): “Pre- and post-processing workflow for affinity purification mass spectrometry data,” J. Proteom. Res., 13, 2239–2249.10.1021/pr401249bSearch in Google Scholar PubMed

Frazee, A. C., G. Pertea, A. E. Jaffe, B. Langmead, S. L. Salzberg and J. T. Leek (2014): “Flexible isoform-level differential expression analysis with Ballgown,” bioRxiv reprint, doi: 10.1101/003665.Search in Google Scholar

Fröhlich, H., Ö. Sahin, D. Arlt, C. Bender and T. Beissbarth (2009): “Deterministic Effects Propagation Networks for reconstructing protein signaling networks from multiple interventions,” BMC Bioinform., 10, 322.10.1186/1471-2105-10-322Search in Google Scholar PubMed PubMed Central

Galati, J. C., K. A. Seaton, K. J. Lee, J. A. Simpson and J. B. Carlin (2014): “Rounding non-binary categorical variables following multivariate normal imputation: evaluation of simple methods and implications for practice,” J. Stat. Comput. Simul., 84, 798–811.10.1080/00949655.2012.727815Search in Google Scholar

Goeman, J. J., S.A. van de Geer, F. de Kort and H. C. van Houwelingen (2004): “A global test for groups of genes: testing association with a clinical outcome,” Bioinformatics, 20, 93–99.10.1093/bioinformatics/btg382Search in Google Scholar PubMed

Griebel, T., B. Zacher, P. Ribeca, E. Raineri, V. Lacroix, R. Guigó and M. Sammeth (2012): “Modelling and simulating generic RNA-Seq experiments with the flux simulator,” Nucleic Acids Res., 40, 10073–10083.10.1093/nar/gks666Search in Google Scholar PubMed PubMed Central

Higham, N. (2002): “Computing the nearest correlation matrix – a problem from finance,” IMA J. Numer. Anal., 22, 329–343.10.1093/imanum/22.3.329Search in Google Scholar

Horton, N. J., S. R. Lipsitz and M. Parzen (2003): “A potential for bias when rounding in multiple imputation,” Am. Stat., 57, 229–232.10.1198/0003130032314Search in Google Scholar

Jung, K., H. Dihazi, A. Bibi, G. H. Dihazi and T. Beissbarth (2014): “Adaption of the global test idea to proteomics data with missing values,” Bioinformatics, 30, 1424–1430.10.1093/bioinformatics/btu062Search in Google Scholar PubMed

Karlis, D. and L. Meligkotsidou (2005): “Multivariate Poisson regression with covariance structure,” Stat. Comput., 15, 255–265.10.1007/s11222-005-4069-4Search in Google Scholar

Kirk, P. D. W. and M. P. H. Stumpf (2009): “Gaussian process regression bootstrapping: exploring the effects of uncertainty in time course data,” Bioinformatics, 25, 1300–1306.10.1093/bioinformatics/btp139Search in Google Scholar PubMed PubMed Central

Kramer, F. (2014): “Integration of Pathway Data as Prior Knowledge into Methods for Network Reconstruction,” Dissertation, Georg-August-Universit at Göttingen.Search in Google Scholar

Kramer, F., M. Bayerlová, F. Klemm, A. Bleckmann and T. Beissbarth (2013): “rBiopaxParser – an R package to parse, modify and visualize BioPAX data,” Bioinformatics, 29, 520–522.10.1093/bioinformatics/bts710Search in Google Scholar PubMed

Kramer, F., M. Bayerlová and T. Beißbarth (2014): “R-based software for the integration of pathway data into bioinformatic algorithms,” Biology, 3, 85–100.10.3390/biology3010085Search in Google Scholar

Ledoit, O. and M. Wolf (2003): “Improved estimation of the covariance matrix of stock returns with an application to portfolio selection,” J. Empir. Financ., 10, 603–621.10.1016/S0927-5398(03)00007-0Search in Google Scholar

Leisch, F., A. Weingessel and K. Hornik (1998): “On the generation of correlated artificial binary data.” Working Papers SFB ‘Adaptive Information Systems and Modelling in Economics and Management Science’, 13. SFB Adaptive Information Systems and Modelling in Economics and Management Science, WU Vienna University of Economics and Business, Vienna.Search in Google Scholar

Li, B. and C. Dewey (2011): “RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome,” BMC Bioinform., 12, 323.10.1201/b16589-5Search in Google Scholar

Li, C.-S., J.-C. Lu, J. Park, K. Kim, P. A. Brinkley and J. P. Peterson (1999): “Multivariate zero-inflated poisson models and their applications,” Technometrics, 41, 29–38.10.1080/00401706.1999.10485593Search in Google Scholar

Liao, Y., G. K. Smyth and W. Shi (2014): “FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features,” Bioinformatics, 30, 923–930.10.1093/bioinformatics/btt656Search in Google Scholar PubMed

Liu, Z., F. Sun, J. Braun, D. P. B. McGovern and S. Piantadosi (2015): “Multilevel regularized regression for simultaneous taxa selection and network construction with metagenomic count data,” Bioinformatics, 31, 1067–1074.10.1093/bioinformatics/btu778Search in Google Scholar PubMed PubMed Central

Mansmann, U. and R. Meister (2006): “Testing differential gene expression in functional groups,” Methods Inf. Med., 44, 449–453.10.1055/s-0038-1633992Search in Google Scholar

Opgen-Rhein, R. and K. Strimmer (2007): “Accurate ranking of differentially expressed genes by a distribution-free shrinkage approach,” Statist. Appl. Genet. Mol. Biol., 6, 9.10.2202/1544-6115.1252Search in Google Scholar PubMed

R Core Team (2013): R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL in Google Scholar

Robinson, M. D., D. J. McCarthy and G. K. Smyth (2010): “edgeR: a Bioconductor package for differential expression analysis of digital gene expression data,” Bioinformatics, 26, 139–140.10.1093/bioinformatics/btp616Search in Google Scholar PubMed PubMed Central

Rustici G., N. Kolesnikov, M. Brandizi, T. Burdett, M. Dylag, I. Emam, A. Farne, E. Hastings, J. Ison, M. Keays, N. Kurbatova, J. Malone, R. Mani, A. Mupo, R. Pedro Pereira, E. Pilicheva, J. Rung, A. Sharma, Y. A. Tang, T. Ternent, A. Tikhonov, D. Welter, E. Williams, A. Brazma, H. Parkinson and U. Sarkans (2013): “ArrayExpress update – trends in database growth and links to data analysis tools,” Nucleic Acids. Res., 31, D987–D990.10.1093/nar/gks1174Search in Google Scholar PubMed PubMed Central

Schaefer, C. F., K. Anthony, S. Krupa, J. Buchoff, M. Day, T. Hannay and K. H. Buetow (2009): “PID: the pathway interaction database,” Nucleic Acids. Res., 37, D674–D679.10.1093/nar/gkn653Search in Google Scholar PubMed PubMed Central

Schäfer, J. and K. Strimmer (2005): “A shrinkage approach to large-scale covariance estimation and implications for functional genomics,” Statist. Appl. Genet. Mol. Biol., 4, 32.10.2202/1544-6115.1175Search in Google Scholar PubMed

Shi, P. and E. A. Valdez (2014): “Multivariate negative binomial models for insurance claim counts,” Insur. Math. Econ., 55, 18–29.10.1016/j.insmatheco.2013.11.011Search in Google Scholar

Shin, K. and R. Pasupathy (2007): “A method for fast generation of bivariate Poisson random vectors,” Proc 2007 Winter Simulation Conf, 472–479.10.1109/WSC.2007.4419637Search in Google Scholar

Yahav, I. and G. Shmueli (2012): “On generating multivariate Poisson data in management science applications,” Appl. Stoch. Model. Bus., 28, 91–102.10.1002/asmb.901Search in Google Scholar

Zhang, L. and B. K. Mallick (2013): “Inferring gene networks from discrete expression data,” Biostatistics, 14, 708–722.10.1093/biostatistics/kxt021Search in Google Scholar PubMed

Zhao, T. and H. Liu (2012): “The huge Package for High-Dimensional Undirected Graph Estimation in R,” J. Mach. Learn. Res., 13, 1059–1062.Search in Google Scholar

Zhou, H., J. Jin, Z. Haojun, Y. Bo, M. Wozniak and W. Limsoon (2012): “IntPath – an integrated pathway gene relationship database for model organisms and important pathogens,” BMC Syst. Biol., 6:Suppl 2, S2.10.1186/1752-0509-6-S2-S2Search in Google Scholar PubMed PubMed Central

Published Online: 2016-9-21
Published in Print: 2016-10-1

©2016 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 1.4.2023 from
Scroll to top button