title: Precise periodic components estimation for chronobiological signals through Bayesian Inference with sparsity enforcing prior.
abstract: The toxicity and efficacy of more than 30 anticancer agents present very high variations, depending on the dosing time. Therefore, the biologists studying the circadian rhythm require a very precise method for estimating the periodic component (PC) vector of chronobiological signals. Moreover, in recent developments, not only the dominant period or the PC vector present a crucial interest but also their stability or variability. In cancer treatment experiments, the recorded signals corresponding to different phases of treatment are short, from 7 days for the synchronization segment to 2 or 3 days for the after-treatment segment. When studying the stability of the dominant period, we have to consider very short length signals relative to the prior knowledge of the dominant period, placed in the circadian domain. The classical approaches, based on Fourier transform (FT) methods are inefficient (i.e., lack of precision) considering the particularities of the data (i.e., the short length). Another particularity of the signals considered in such experiments is the level of noise: such signals are very noisy and establishing the periodic components that are associated with the biological phenomena and distinguishing them from the ones associated with the noise are difficult tasks. In this paper, we propose a new method for the estimation of the PC vector of biomedical signals, using the biological prior informations and considering a model that accounts for the noise. The experiments developed in cancer treatment context are recording signals expressing a limited number of periods. This is a prior information that can be translated as the sparsity of the PC vector. The proposed method considers the PC vector estimation as an Inverse Problem (IP) using the general Bayesian inference in order to infer the unknown of our model, i.e. the PC vector but also the hyperparameters (i.e the variances). The sparsity prior information is modeled using a sparsity enforcing prior law. In this paper, we propose a Student'stdistribution, viewed as the marginal distribution of a bivariate normal-inverse gamma distribution. We build a general infinite Gaussian scale mixture (IGSM) hierarchical model where we assign prior distributions also for the hyperparameters. The expression of the joint posterior law of the unknown PC vector and hyperparameters is obtained via Bayes rule, and then, the unknowns are estimated via joint maximum a posteriori (JMAP) or posterior mean (PM). For the PM estimator, the expression of the posterior distribution is approximated by a separable one, via variational Bayesian approximation (VBA), using the Kullback-Leibler (KL) divergence. For the PM estimation, two possibilities are considered: an approximation with a partially separable distribution and an approximation with a fully separable one. Both resulting algorithms corresponding to the PM estimation and the one corresponding to the JMAP estimation are iterative algorithms. The algorithms are presented in detail and are compared with the ones corresponding to the Gaussian model. We examine the convergency of the algorithms and give simulation results to compare their performances. Finally, we show simulation results on synthetic and real data in cancer treatment applications. The real data considered in this paper examines the rest-activity patterns of KI/KI Per2::luc mouse, aged 10 weeks, singly housed in RealTime Biolumicorder (RT-BIO).
authors.affiliation: Laboratoire des signaux et systèmes (L2S), UMR 8506 CNRS-CentraleSupélec-Univ. Paris-Sud, CentraleSupélec, Plateau de Moulon, Gif-sur-Yvette, 91192 France ; Rythmes Biologiques et Cancers (RBC), UMR 776 INSERM-Univ. Paris-Sud, Campus CNRS, Villejuif, 94801 France.
authors.forename: Mircea
authors.initials: M
authors.lastname: Dumitru
authors.affiliation: Laboratoire des signaux et systèmes (L2S), UMR 8506 CNRS-CentraleSupélec-Univ. Paris-Sud, CentraleSupélec, Plateau de Moulon, Gif-sur-Yvette, 91192 France.
authors.forename: Ali
authors.initials: A
authors.lastname: Mohammad-Djafari
authors.affiliation: Laboratoire des signaux et systèmes (L2S), UMR 8506 CNRS-CentraleSupélec-Univ. Paris-Sud, CentraleSupélec, Plateau de Moulon, Gif-sur-Yvette, 91192 France ; Department of Molecular Biotechnology and Health Sciences, University of Turin, Turin, 10126 Italy.
authors.forename: Simona Baghai
authors.initials: SB
authors.lastname: Sain
keywords: Bayesian parameter estimation; Cancer treatment; Chronobiology; Circadian rhythm; Infinite Gaussian scale mixture (IGSM); Inverse problem; Joint maximum a posteriori (JMAP); Kullback-Leibler (KL) divergence; Normal-inverse gamma; Periodic component (PC) vector estimation; Posterior mean (PM); Sparsity enforcing; Variational Bayesian approximation (VBA)
references.citation: Chronobiol Int. 2009 Feb;26(2):258-81
references.pmid: 19212840
references.citation: BMC Syst Biol. 2011 Jun 24;5:97
references.pmid: 21702910
references.citation: Int J Cancer. 2012 Dec 1;131(11):2684-92
references.pmid: 22488038
references.citation: J Natl Cancer Inst. 2000 Jun 21;92(12):994-1000
references.pmid: 10861311
references.citation: Biostatistics. 2011 Oct;12(4):682-94
references.pmid: 21551122
references.citation: Bioinformatics. 2013 Dec 15;29(24):3174-80
references.pmid: 24058056
references.citation: Cancer. 2003 Jan 1;97(1):155-69
references.pmid: 12491517
references.citation: Annu Rev Pharmacol Toxicol. 2010;50:377-421
references.pmid: 20055686
references.citation: Clin Cancer Res. 2000 Aug;6(8):3038-45
references.pmid: 10955782
references.citation: Nat Rev Cancer. 2003 May;3(5):350-61
references.pmid: 12724733
references.citation: Chronobiol Int. 2011 Aug;28(7):586-600
references.pmid: 21859417
references.citation: PLoS One. 2011;6(6):e20393
references.pmid: 21674030
references.citation: Proc Natl Acad Sci U S A. 2009 Mar 31;106(13):5394-9
references.pmid: 19282474
references.citation: Biol Rev Camb Philos Soc. 2004 Aug;79(3):533-56
references.pmid: 15366762
references.citation: Annu Rev Pharmacol Toxicol. 2007;47:593-628
references.pmid: 17209800
references.citation: BMC Syst Biol. 2012 Jun 15;6:69
references.pmid: 22703599
references.citation: Annu Rev Physiol. 2010;72:517-49
references.pmid: 20148687
references.citation: Annu Rev Neurosci. 2012;35:445-62
references.pmid: 22483041
references.citation: Nat Rev Neurosci. 2003 Aug;4(8):649-61
references.pmid: 12894240
references.citation: Stat Appl Genet Mol Biol. 2014 Jun;13(3):227-73
references.pmid: 24864301
references.citation: Cancer Res. 2009 Jun 1;69(11):4700-7
references.pmid: 19470769
references.citation: Genes Dev. 2013 Jul 1;27(13):1526-36
references.pmid: 23824542
references.citation: Lancet Oncol. 2001 May;2(5):307-15
references.pmid: 11905786
references.citation: Cold Spring Harb Symp Quant Biol. 2011;76:39-47
references.pmid: 22179985
references.citation: Handb Exp Pharmacol. 2013;(217):261-88
references.pmid: 23604483
references.citation: Cancer Res. 2013 Dec 15;73(24):7176-88
references.pmid: 24154875
references.citation: J Natl Cancer Inst. 2002 May 1;94(9):690-7
references.pmid: 11983758
references.citation: Biostatistics. 2013 Sep;14(4):792-806
references.pmid: 23743206
references.citation: J Natl Cancer Inst. 2005 Apr 6;97(7):507-17
references.pmid: 15812076