Non-stationary statistical modeling of circadian oscillations for research in cancer chronotherapy

Lead Research Organisation: University of Warwick
Department Name: Statistics

Abstract

Oscillators are found on all levels in Nature. One of the most prominent examples is the circadian rhythm which displays an endogenous, entrainable oscillation of about 24 hours. Chronotherapy is a form of personalized therapy where treatments are administered according to a schedule that corresponds to a person's rhythms in order to maximize effectiveness and minimize side effects of the therapy. There is now large evidence both in experimental models and in patients that the Circadian Timing System (CTS) of patients is affecting treatment tolerability and efficacy of anticancer drugs while the progression of the disease and treatment with anticancer drugs have disruptive effect on the CTS. Our projects focuses on the development of statistical methods that address the modeling of such time varying effects. An important aim is to develop a real-time prediction system of a patient's CTS state in order to achieve an improved personalized chronotherapy.

Technical Summary

There is now large evidence both in experimental models and in patients that the Circadian Timing System (CTS) of the patients is affecting treatment tolerability and efficacy of anticancer drugs while, on the other hand, the development of the disease and the treatment may have a disruptive effect on the CTS. For example, Costa et al. (2013) found a significant shift in period length during anticancer treatment suggesting that the patient's daily rhythm is altered while receiving anticancer drugs. Hence a personalized chronotherapy needs to take into account the potential effect of the treatment as well as the disease on the CTS of the patient. The aim of this project is develop a statistical methodology that is able to detect and quantify (in real time) such temporal variation in the circadian oscillation and that can be used to study how these are associated with changes in covariates (tumour growth, treatment with drugs and timing of their administration, external factors such as light, temperature, food and other potential Zeitgebers). We address this through non-stationary statistical modeling with the aim of developing a system for online prediction of a patient's CTS. Using spectrum and state-space approaches we will develop non-stationary statistical methods using data from all available biomarkers to construct a robust online prediction of the current state of a patient's CTS that can be used to individualize the timing of chronotherapy. Secondly, we will enhance a mechanistic mathematical modeling approach to gene expression and provide a statistical algorithm of fitting such a model to circadian gene expression data with the aim of gaining a better understanding of how the underlying gene regulatory processes might be affected by the disease and other covariates.

Planned Impact

This programme will develop some methodologies that will make it possible to analyze the impact of changing conditions on the circadian system in living organism. The PI has been involved for over 10 years in interdisciplinary research in systems biology and has extensive experience with the analysis of circadian data and the modeling gene expression. CI 1 Prof. Francis Levi is a world expert on chronotherapy in cancer who will direct the scientific content of our project while CI 2 will bring essential domain expertise concerning mathematical modeling in systems biology, particularly circadian systems.

We envisage an impact of this research on

(i) Statisticians working to provide appropriate statistical methods for systems biology and medicine and the statistical community working in the area of statistical signal processing.

(ii) Theoretical and experimental biologists working in the area of chronobiology.
Methods for modeling and estimating the time varying nature of circadian or other oscillations as a result of environmental or experimental conditions, with
online updating of information bases supported by appropriate inferential methodologies will continue to be developed and implemented in a wide range of applications.

(iii) There is large evidence of the Circadian Timing System (CTS) modulation of treatment tolerability and efficacy in experimental models and in cancer patients. The chronotherapy concept has thus obvious implications both for daily patient care and outcomes, and for drug development.

(iv) Our methods in assessing a patient's CTS in real time may also benefit researchers and patients of
chronotherapy in other fields, examples of this are the treatment of asthma, hypertension, and multiple types of depression.

The types of impact that we expect from this project are
- Impact on Knowledge, in particular on techniques (through beneficiaries in groups [i] and [ii]).
- Impact on Society, in particular impact on health and quality of life (through beneficiaries in group [iii]
and [iv] above).
- Impact on People, in particular on skills (through training the RA involved and facilitating the involvement of postdoctoral researchers in the workshops.

Publications

10 25 50
 
Description Time Series module for Doctoral Training Centre
Geographic Reach Local/Municipal/Regional 
Policy Influence Type Influenced training of practitioners or researchers
Impact improvement of statistical skills oriented towards analysing patient time series data obtained from sensors
 
Description IDEAs (Identification of DEterminants of Altered circadian rhythm) study.
Amount £400,000 (GBP)
Organisation Philips Research 
Sector Private
Country United States
Start 04/2017 
End 03/2019
 
Title R toolbox on time changing spectral analysis 
Description Statistical toolbox written in R to analyse and plot time changing spectral properties (for instance change in periodicities) of circadian data. 
Type Of Material Improvements to research infrastructure 
Provided To Others? No  
Impact Only just started to use it within the group. Further development of the toolbox might be necessary before providing it to other groups. 
 
Title R-codes to obtain results for homogeneous and harmonic HMMs 
Description Wearable computing devices allow collection of densely sampled real-time information on movement enabling researchers and medical experts to obtain objective and non-obtrusive records of actual activity of a subject in the real world over many days. Our interest is motivated by the use of activity data for evaluating and monitoring the circadian rhythmicity of subjects for research in chronobiology and chronotherapeutic healthcare. In order to translate the information from such high-volume data arising we propose the use of a hidden Markov modelling (HMM) approach which (i) naturally captures the notable square wave form observed in activity data along with heterogeneous ultradian variances over the circadian cycle of human activity, (ii) thresholds activity into different states in a probabilistic way while respecting time dependence and (iii) gives rise to circadian rhythm parameter estimates, based on probabilities of transitions between rest and activity, that are interpretable and of interest to circadian research. R-codes are available to obtain results for homogeneous and harmonic HMMs. 
Type Of Material Improvements to research infrastructure 
Year Produced 2018 
Provided To Others? Yes  
Impact The paper is only just published a month ago. However, the methodology/R-code is used (and further adapted) to analyse new data sets arising in the new collaborating project CIRCADIEM 
URL https://github.com/huang1010/hmms
 
Title Generated novel datasets by recruiting people (students, friends, family, neighbours) wearing 
Description Generated 19 new datasets by recruiting people (members of this team including PI and PDRA, students, friends, family, neighbours) to wear actigraph sensor located on chest as well as around wrist for 7 days. 
Type Of Material Database/Collection of data 
Provided To Others? No  
Impact Data sets will be used to develop statistical methodology 
 
Title Time varying Spectral analysis method 
Description Time varying spectral analysis method coded in R that detects changes in period, phase and amplitudes over time. Can be applied to periodic time series data (for actigraphy, gene expression, physiological data, etc to detect strength of circadian Rhythm) 
Type Of Material Data analysis technique 
Year Produced 2016 
Provided To Others? No  
Impact The code has been passed to company Altran for industrial exploitation 
 
Title Use of Hidden Markov Models for Actigraphy data 
Description Investigated the application of Hidden Markov models for modelling actigraphy data. Developed R-codes in addition to using pre-existing codes to fit these to actigraphy data. Developed new visualisation methods to analyse the data. 
Type Of Material Computer model/algorithm 
Provided To Others? No  
Impact The results are currently written up as a research paper 
 
Description Joined Chronotherapy group at Warwick 
Organisation University of Warwick
Department Warwick Medical School
Country United Kingdom 
Sector Academic/University 
PI Contribution PI Barbel Finkenstadt and PDRA Qi Huang joined Chronotherapy group at Warwick Medical School. There are regular meetings every 2 weeks where we have opportunity to regularly present our work work to chronotherapy researchers at the Medical School.
Collaborator Contribution We gain scientific insight from the group that is useful for our modelling approaches.
Impact Collaboration is multidisciplinary: Medicine, Biological Sciences, Mathematics, Statistics
Start Year 2015
 
Description Project CIRCADIEM 
Organisation National Institute of Health and Medical Research (INSERM)
Country France 
Sector Academic/University 
PI Contribution Statistical Methodology developed in project are used by the project CIRCADIEM which investigates altimetric data on shiftworking nurses at Paul Brousse Hospital, Villejuif, France
Collaborator Contribution New data sets on shiftworking individuals (about 160 nurses) will be available to us. These will help in further developing our methods and learning parameters to investigate circadian disruption for shift workers versus 'normal' day working individuals.
Impact Data sets are available for further research and the project will be listed as collaborators in planned grant application. As of October 2018 a PhD student in Statistics, Warwick, is working on the analysis of this data in co-supervision between Finkenstadt and Brettschneider
Start Year 2017
 
Description Inference for circadian pacemaking 
Form Of Engagement Activity Participation in an open day or visit at my research institution
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Other audiences
Results and Impact Offer-holder visit day, future students and their parents. Talk to make students enthusiastic about statistics and its interesting applications across disciplines
Year(s) Of Engagement Activity 2019
 
Description Invited Talk at Workshop "The Interplay of Stochastic and Deterministic Dynamics in Networks" (MBI Ohio) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Other audiences
Results and Impact Academic conference attracting some famous and well known Mathematicians in Mathematical Biology
Year(s) Of Engagement Activity 2016
URL https://mbi.osu.edu/event/?id=897
 
Description OxWaSP workshop: Talk on "Large Data Sets and Complex Models: A view from Systems Biology" 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach Local
Primary Audience Postgraduate students
Results and Impact Talk at Oxford/Warwick Doctoral Training centre Workshop. The Talk reached about 40 graduate students (some of them on search for PhD topic and the hope was to attract future PhD student) and about 10-15 other academics from Warwick and Oxford.
Year(s) Of Engagement Activity 2015
 
Description Symposium Medicine in the 4rth dimension 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Study participants or study members
Results and Impact The 2nd Warwick Quantitative Biomedicine Programme Symposium brought together leading scientists and clinicians to discuss systems biology approaches to medically relevant dynamic processes.
Year(s) Of Engagement Activity 2016
URL http://www2.warwick.ac.uk/fac/med/qbp/symposia/symposium2/
 
Description Talk at UK clock club summer 2016 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Postgraduate students
Results and Impact 20 min Talk by PDRA on Non-Stationary Spectral Analysis for Circadian Data at UK Clock club meeting on 4/7/2017 at Warwick Medical School
The talk was on this occasion also attended by international audience from US as it was in connection with a symposium held at WArwick.
Year(s) Of Engagement Activity 2016
URL http://www2.warwick.ac.uk/fac/med/news/events/ukclockclub/