A big data pipeline: Identifying dynamic gene regulatory networks from time-course Gene Expression Omnibus data with applications to influenza infection

Abstract

A biological host response to an external stimulus or intervention such as a disease or infection is a dynamic process, which is regulated by an intricate network of many genes and their products. Understanding the dynamics of this gene regulatory network allows us to infer the mechanisms involved in a host response to an external stimulus, and hence aids the discovery of biomarkers of phenotype and biological function. In this article, we propose a modeling/analysis pipeline for dynamic gene expression data, called Pipeline4DGEData, which consists of a series of statistical modeling techniques to construct dynamic gene regulatory networks from the large volumes of high-dimensional time-course gene expression data that are freely available in the Gene Expression Omnibus repository. This pipeline has a consistent and scalable structure that allows it to simultaneously analyze a large number of time-course gene expression data sets, and then integrate the results across different studies. We apply the proposed pipeline to influenza infection data from nine studies and demonstrate that interesting biological findings can be discovered with its implementation.

Keywords

Time-course data differential equations gene regulatory network

Get full access to this article

View all access options for this article.

References

Edgar

Domrachev

Lash

. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 2002; 30: 207–210.

Barrett

Wilhite

Ledoux

et al.

NCBI GEO: archive for functional genomics data sets – update. Nucleic Acids Res 2013; 41: D991–D995.

NCBI. GEO Summary 2016, www.ncbi.nlm.nih.gov/geo/summary/ (accessed 26 November 2017).

De Jong

. Modeling and simulation of genetic regulatory systems: a literature review. J Comput Biol 2002; 9: 67–103.

Hecker

Lambeck

Toepfer

et al.

Gene regulatory network inference: data integration in dynamic models – a review. Biosystems 2009; 96: 86–103.

Hood

Friend

. Predictive, personalized, preventive, participatory (P4) cancer medicine. Nature Rev Clin Oncol 2011; 8: 184–187.

Khoury

Gwinn

Glasgow

et al.

A population approach to precision medicine. Am J Prevent Med 2012; 42: 639–645.

Steuer

Kurths

Daub

et al.

The mutual information: detecting and evaluating dependencies between variables. Bioinformatics 2002; 18(suppl 2): S231–S240.

Stuart

Segal

Koller

et al.

A gene-coexpression network for global discovery of conserved genetic modules. Science 2003; 302: 249–255.

10.

Thomas

. Boolean formalization of genetic control circuits. J Theor Biol 1973; 42: 563–585.

11.

Shmulevich

Dougherty

Kim

et al.

Probabilistic Boolean networks: a rule-based uncertainty model for gene regulatory networks. Bioinformatics 2002; 18: 261–274.

12.

Friedman

Linial

Nachman

et al.

Using Bayesian networks to analyze expression data. J Comput Biol 2000; 7: 601–620.

13.

Kim

Imoto

Miyano

. Inferring gene networks from time series microarray data using dynamic Bayesian networks. Brief Bioinform 2003; 4: 228–235.

14.

Perrin

Ralaivola

Mazurie

et al.

Gene networks inference using dynamic Bayesian networks. Bioinformatics 2003; 19(suppl 2): ii138–ii148.

15.

Holter

Maritan

Cieplak

et al.

Dynamic modeling of gene expression data. Proc Natl Acad Sci 2001; 98: 1693–1698.

16.

Sakamoto E and Iba H. Inferring a system of differential equations for a gene regulatory network by using genetic programming. In: Proceedings of the 2001 congress on evolutionary computation (vol. 1). Hong Kong, China: IEEE, 2001, pp.720–726.

17.

Cargill

Altshuler

Ireland

et al.

Characterization of single-nucleotide polymorphisms in coding regions of human genes. Nat Genet 1999; 22: 231–238.

18.

Affymetririx. Statistical algorithms description document. Affymetrix white papers, 2002.

19.

Irizarry

Hobbs

Collin

et al.

Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003; 4: 249–264.

20.

Irizarry

Gentleman

et al.

A model-based background adjustment for oligonucleotide expression arrays. J Am Stat Assoc 2004; 99: 909–917.

21.

Millenaar

Okyere

May

et al.

How to decide? Different methods of calculating gene expression from short oligonucleotide array data will give different results. BMC Bioinform 2006; 7: 137–137.

22.

Lim

Wang

Lefebvre

et al.

Comparative analysis of microarray normalization procedures: effects on reverse engineering gene networks. Bioinformatics 2007; 23: i282–i288.

23.

Green

Silverman

. Nonparametric regression and generalized linear models: a roughness penalty approach, Florida: CRC Press, 1993.

24.

Silverman

Ramsay

. Functional data analysis, New York: Springer, 2005.

25.

Golub

Heath

Wahba

. Generalized cross-validation as a method for choosing a good ridge parameter. Technometrics 1979; 21: 215–223.

26.

Yao

Müller

Wang

. Functional data analysis for sparse longitudinal data. J Am Stat Assoc 2005; 100: 577–590.

27.

. More powerful significant testing for time course gene expression data using functional principal component analysis approaches. BMC Bioinform 2013; 14: 6–6.

28.

Carey

Gan

et al.

Correlation-based iterative clustering methods for time course data: the identification of temporal gene response modules for influenza infection in humans. Infect Dis Model. 2016; 1: 28–39, http://dxdoiorg/101016/jidm201607001.

29.

Huang

Sherman

Lempicki

. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 2009; 4: 44–57.

30.

Huang

Sherman

Lempicki

. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 2009; 37: 1–13.

31.

Liang

et al.

High-dimensional ODEs coupled with mixed-effects modeling techniques for dynamic gene regulatory network identification. J Am Stat Assoc 2011; 106: 1242–1258.

32.

Wu S, Liu ZP, Qiu X, et al. High-dimensional ordinary differential equation models for reconstructing genome-wide dynamic regulatory networks. In: Hu M, Liu Y and Lin J (eds) Topics in applied statistics: 2012 Symposium of the International Chinese Statistical Association, Springer Proceedings in Mathematics & Statistics. New York: Springer, 2013, pp.173–190.

33.

Hemker

. Numerical methods for differential equations in system simulation and in parameter estimation. In: Hemker HC and Hess B (eds) Anal Simul Biochem Syst. North Holland, Amsterdam: Elsevier, 1972; 25: 59–80.

34.

Bard Y. Nonlinear Parameter Estimation, Academic Press, New York, 1974.

35.

Voit

Almeida

. Decoupling dynamical systems for pathway identification from metabolic profiles. Bioinformatics 2004; 20: 1670–1681.

36.

Liang

. Parameter estimation for differential equation models using a framework of measurement error in regression models. J Am Stat Assoc 2008; 103: 1570–1583.

37.

Tibshirani

. Regression shrinkage and selection via the lasso. J R Stat Soc B 1996; 58: 267–288.

38.

Zou

. The adaptive lasso and its oracle properties. J Am Stat Assoc 2006; 101: 1418–1429.

39.

Fan

. Variable selection via nonconcave penalized likelihood and its oracle properties. J Am Stat Assoc 2001; 96: 1348–1360.

40.

Bates DM and Watts DG. Nonlinear regression: iterative estimation and linear approximations. In: Bates, Douglas M and Watts, Donald G (eds) Nonlinear regression analysis and its applications. John Wiley & Sons, Inc., 2008, pp. 32–66.

41.

Huber

Carey

Long

et al.

Graphs in molecular biology. BMC Bioinform 2007; 8: 1–1.

42.

Lee

Hsu

Sajdak

et al.

Coexpression analysis of human genes across many microarray data sets. Genome Res 2004; 14: 1085–1094.

43.

Konganti

Wang

Yang

et al.

SBEToolbox: a Matlab toolbox for biological network analysis. Evol Bioinform 2013; 9: 355–362.

44.

Nieminen

. On the centrality in a graph. Scand J Psychol 1974; 15: 332–336.

45.

Jeong

Mason

Barabási

et al.

Lethality and centrality in protein networks. Nature 2001; 411: 41–42.

46.

Koschützki

Schreiber

. Centrality analysis methods for biological networks and their application to gene regulatory networks. Gene Regul Syst Biol 2008; 2: 193–201.

47.

Bavelas

. Communication patterns in task-oriented groups. J Acoust Soci Am 1950; 22: 725–730.

48.

Hage

Harary

. Eccentricity and centrality in networks. Social Netw 1995; 17: 57–63.

49.

Barrenas

Chavali

Holme

et al.

Network properties of complex human disease genes identified through genome-wide association studies. PloS One 2009; 4: e8090–e8090.

50.

Hwang W, Cho Yr, Zhang A, et al. Bridging centrality: identifying bridging nodes in scale-free networks. In: Proceedings of the 12th ACM SIGKDD international conference on knowledge discovery and data mining, 2008, Las Vegas, Nevada, USA: ACM, pp.336–344.

51.

Hallinan JS. Cluster analysis of the p53 genetic regulatory network: topology and biology. In: Proceedings of the 2004 IEEE symposium on computational intelligence in bioinformatics and computational biology. CIBCB’04. La Jolla, CA, USA: IEEE, 2004, pp.1–8.

52.

Palumbo

Colosimo

Giuliani

et al.

Functional essentiality from topology features in metabolic networks: a case study in yeast. FEBS Lett 2005; 579: 4642–4646.

53.

Watts

Strogatz

. Collective dynamics of small-world networks. Nature 1998; 393: 440–442.

54.

Wasserman S and Faust K. Social network analysis: methods and applications (vol. 8). Cambridge, United Kingdom: Cambridge University Press, 1994.

55.

Cai

Borenstein

Petrov

. Broker genes in human disease. Genome Biol Evol 2010; 2: 815–825.

56.

Huang

Chen

. Graph theory and stability analysis of protein complex interaction networks. IET Syst Biol 2016; 10: 64–75.

57.

Grewal

Singh

Chand

. Effect of aggregation operators on network-based disease gene prioritization: a case study on blood disorders. IEEE/ACM Trans Comput Biol Bioinform 2016; 99: 1–1.

58.

(WHO) WHO. Influenza (Seasonal) fact-sheet, www.who.int/mediacentre/factsheets/fs211/en/ (2016, accessed 26 November 2017).

59.

Woods

McClain

Chen

et al.

A host transcriptional signature for presymptomatic detection of infection in humans exposed to influenza H1N1 or H3N2. PLoS One 2013; 8: e52198–e52198.

60.

Djavani

Crasta

Zhang

et al.

Gene expression in primate liver during viral hemorrhagic fever. Virol J 2009; 6: 1–1.

61.

Shapira

Gat-Viks

Shum

et al.

A physical and regulatory map of host-influenza interactions reveals pathways in H1N1 infection. Cell 2009; 139: 1255–1267.

62.

Jin

Pan

et al.

Characterizing and controlling the inflammatory network during influenza A virus infection. Sci Rep 2014; 4: 3799–3799.

63.

Huang

Zaas

Rao

et al.

Temporal dynamics of host molecular responses differentiate symptomatic and asymptomatic influenza a infection. PLoS Genet 2011; 7: e1002234–e1002234.

64.

Thakar

Hartmann

Marjanovic

et al.

Comparative analysis of anti-viral transcriptomics reveals novel effects of influenza immune antagonism. BMC Immunol 2015; 16: 46–46.

65.

Niu

Chasman

Eisfeld

et al.

Multi-task consensus clustering of genome-wide transcriptomes from related biological conditions. Bioinformatics 2016; 32: 1509–1517.

66.

Mitchell

Eisfeld

Sims

et al.

A network integration approach to predict conserved regulators related to pathogenicity of influenza and SARS-CoV respiratory viruses. PloS One 2013; 8: e69374–e69374.

67.

Itkes AV. Oligoadenylate and cyclic AMP: interrelation and mutual regulation. In: Müller WEG. and Schröder HC (eds) Biological Response Modifiers – Interferons, Double-Stranded RNA and 2’,5’-Oligoadenylates. Progress in Molecular and Subcellular Biology, vol 14. Springer, Berlin, Heidelberg, 1994, pp. 209–221.

68.

Cao

Wei

Zhao

et al.

DnaJA1/Hsp40 is co-opted by influenza A virus to enhance its viral RNA polymerase activity. J Virol 2014; 88: 14078–14089.

69.

Bohn

BeMiller

. (1 3)-β-D-Glucans as biological response modifiers: a review of structure-functional activity relationships. Carbohydr Polym 1995; 28: 3–14.

70.

LoVerme

La Rana

Russo

et al.

The search for the palmitoylethanolamide receptor. Life Sci 2005; 77: 1685–1698.

71.

Yuan

. Drugs to cure avian influenza infection – multiple ways to prevent cell death. Cell Death Dis 2013; 4: e835–e835.

72.

Davidson

Crotta

McCabe

et al.

Pathogenic potential of interferon αβ in acute influenza infection. Nature Commun 2014; 5: 3864–3864.

73.

Dubois

Terrier

Rosa-Calatrava

. Influenza viruses and mRNA splicing: doing more with less. MBio 2014; 5: e00070–14.

74.

Dudek

Luig

Pauli

et al.

The clinically approved proteasome inhibitor PS-341 efficiently blocks influenza A virus and vesicular stomatitis virus propagation by establishing an antiviral state. J Virol 2010; 84: 9439–9451.

75.

Turpin

Luke

Jones

et al.

Influenza virus infection increases p53 activity: role of p53 in cell death and viral replication. J Virol 2005; 79: 8802–8811.

76.

Husain

Harrod

. Enhanced acetylation of alpha-tubulin in influenza A virus infected epithelial cells. FEBS Lett 2011; 585: 128–132.

77.

Han

Chen

et al.

Influenza virus A/Beijing/501/2009 (H1N1) NS1 interacts with β-tubulin and induces disruption of the microtubule network and apoptosis on A549 cells. PLoS One 2012; 7: e48340–e48340.

78.

Jackson

Dowling

Spiesman

et al.

Transmission of the common cold to volunteers under controlled conditions: I. The common cold as a clinical entity. AMA Arch Intern Med 1958; 101: 267–278.

79.

Hubert

Arabie

. Comparing partitions. J Classif 1985; 2: 193–218.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.55 MB