A Linear Syst0ems Toolkit for Biology

Lead Research Organisation: University of Cambridge
Department Name: Plant Sciences

Abstract

Biology is complex; cells are made up of 1000s of proteins, a similar number of metabolites and tens of thousands of genes. A goal of biological research is to understand how this complexity brings about the functions of life. One way to achieve this goal is to understanding the connections between the 1000s of components that make up cells. Measuring the connections between all the components is challenging, particularly because cells are dynamical systems that are constantly changing. Accurate descriptions of the dynamical network interactions that take place in a cell are required to make the advances required for improved crops for food security and new medicines.
We have adapted a new tool set from Engineering to describe biological networks in a mathematical form. We make models of each of the connections which are used to predict how the system will change over time, which is very useful in discovering how cells respond to signals such as changes in temperature, hormones or drugs. Our new mathematical tool set allows researchers to identify and quantify the changes in a biological network, which can lead to the discovery of the gene(s) or pathways that are involved in responses to stresses or drugs and might underlie disease. Our new mathematical tool set will have wide utility in understanding a wide range of cellular systems, from the effects of drugs in humans to the response of a crop plant to environmental changes or attack by pests. Our development of a tool that measures how biological networks change is important for understanding biology, curing disease and improving crop plants to provide enhanced food security.
We propose to develop this so called Nu gap analysis as a practical tool for biologists. In our implementation, we identify and describe connections in biological systems using simple liner models. The Nu gap measures the difference between the mathematical descriptions of the connections obtained in different conditions, such as following a response to a drug, or an environmental stress.
To develop the Nu gap as a practical tool we will undertake a research programme that increases with complexity over time. This will permit rigorous testing, development and deployment of Nu gap analyses. First, we will perform theoretical analyses of the Nu gap on models derived from fabricated datasets designed specifically to assess the strengths and limitations of the Nu gap. This will inform as to where application of the toolset would be best, and conversely the situations where the Nu gap might be less informative. Having developed good theoretical understanding of the system, we will apply the Nu gap to real world data obtained by our laboratories. We will begin using data describing the circadian regulation of gene expression in the model plant Arabidopsis. A major goal will be to investigate the effect of a pharmacological and a genetic perturbation to the circadian system. Both profoundly affect the functioning of the circadian clock, but the mechanisms by which these affect the circadian clock is uncertain.
We will move from investigating the fundamental properties of the circadian clock in the model plant Arabidopsis to using linear modelling and Nu gap analyses to describe the circadian clock in a major crop, barley. The circadian clock regulates many important agronomic traits such as flowering time, seed set and cold tolerance. Our studies have the potential to inform breeders of useful gene targets. Recognising that biological systems are more than a series of interactions between genetic components we will extend our analysis to incorporate the physiology of the cell, such as changes in the concentration of calcium in the cytosol, which act as key regulators of signalling in stressful conditions.

Technical Summary

We have developed a tool set from Engineering to identify the causal connections in biological networks. We use linear time invariant (LTI) models to describe the dynamic relationships in biological systems based on analysis of time series datasets. The LTI models describe causal relationships in the network and have predictive power concerning the dynamical system. We have used LTI modelling successfully to describe the circadian clock of Arabidopsis.
We wish to build on this advance by developing a new approach based on the Nu gap metric, which identifies the causal changes in a network in response to stimulation or perturbation, (e.g. pharmacological agents or genetic mutation). The Nu gap identifies those connections that are altered by measuring the degree of change in the LTI models that describe those connections. Identification of those network connections that have changed in response to treatment allows follow up studies to be focused exclusively on the affected nodes that represent primary candidate gene targets.
Our preliminary studies demonstrate the utility of our approach. LTI modelling and Nu gap analyses of circadian transcriptomes of Arabidopsis has identified the target for the metabolite nicotinamide in the circadian clock and we have confirmed the mathematical prediction through experimentation.
We will develop LTI modelling with Nu gap analysis both theoretically and practically. We will apply Nu gap analysis to circadian datasets obtained both from model systems that are well understood and from crop plants in which circadian networks are not fully understood. We will use LTI modelling coupled with Nu gap analysis to investigate transcriptional responses to pharmacological and genetic. We will extend the utility of the LTI modelling with Nu gap analysis by incorporating non-transcriptional data. Theoretical developments will attempt to extend the power of Nu gap analysis to non-linear models.

Planned Impact

WHO WILL BENEFIT?
(1) Academic scientists interested in circadian rhythms, crop biology, stress physiology of plants and systems biologists in all organisms
(2) Industrial scientists interested in generating crop varieties with enhanced stress tolerance and those using systems tools for gene and drug discovery.
(3) Research staff.
(4) The general public.

HOW WILL THEY BENEFIT?
(1) We will ensure wide dissemination and use of the Nu gap in systems biology. (a) We will organise a one day training course in Nu gap techniques at Cambridge for doctoral and pre-doctoral scientists. We have requested funds to support travel and subsistence for up to 20 participants. We will advertise the Nu gap training course through our contact networks, the Cambridge Networks Network, UKPSF, GARNet and Engineering and Systems Biology message boards. Materials associated with the training course will be made available via our websites. (b) We will present our research at the following conferences; IEE Conference on Decision and Control, the International Conference on Systems Biology and other appropriate meetings. These will also be used to advertise the Nu gap training course. (c) We will ensure maximum impact by publishing our research in a timely manner. The applicants have a track record of publishing in high impact journals and widespread dissemination.

(2) Industrial scientists will benefit because we will develop a new STEM tool that allows identification of the causal changes in biological systems that occur in response to stimulation. We envisage that this will have utility in both the agricultural and pharmaceutical industries. In agriculture the Nu gap might be implemented to identify candidate genes for breeding programmes. In the pharmaceutical industry Nu gap could be exceptionally powerful in identifying drug targets and may offer considerable advantage compared to correlative tools currently in use. We will use our current relationships with industrial partners at Bayer Crop Science and Microsoft Research to attract industrially-based scientists to the Nu gap training course.

(3) The PDRA will gain considerable benefit from being employed on the project. This will include training in circadian systems and linear modelling tools. The training in the specialist control theory approaches will place the PDRA in a good position for a further career in academia or the pharmaceutical, agricultural, financial or engineering industries. PDRAs from the Webb laboratory have had excellent career advancement. All BBSRC-funded PDRAs in the Webb laboratory have obtained publications in Science or Nature and eight former members of the Webb laboratory have obtained Faculty positions. The PDRA will gain considerable experience on helping develop and deliver the Nu gap training course and associated material.

(4) The general public will benefit from outreach activities at the Department of Plant Sciences, Cambridge. During Science Week numerous interactive and more formal displays on aspects of plant biology and research are presented and 7,000 visit the Plant Sciences displays on 'Science Saturday' which will include dissemination of findings from this project. We take every opportunity to publicise our findings, Dr Webb has appeared on Radio interviews (e.g. BBC Farming Today) and his recent findings have been summarised in media outlets as diverse as the Financial Times and Comedy Central's Colbert Report. It is hoped in the long term that the public will benefit from food security generated from the novel agricultural products that arise from our findings. Whilst recognising that in any field of study the translation rate from laboratory finding to industrial product is always low, we make every effort with our industrial partners (Bayer Cropscience) to translate our findings for public benefit. We are currently registering IP on one of our discoveries
 
Description We have been developed a new tools that allows analysing data sets from biological systems to identify how biological systems respond to biological events. We performed our tool development using data from the analysis of circadian rhythms in crops and the model species Arabidopsis. We have developed an approach using linear models and metric that measures the difference in models called the nu gap. We have applied this to data from Arabidopsis treated with a chemical and also different varieties of barley to find the changes in the circadian clock that result in different behaviours. We have found that the two circadian clock genes PRR7 and PRR9 regulate circadian period. We also found that calcium ions regulate circadian period. We have also used our tools to generate the first mathematical model of the barley circadian clock. This was made possible by our new mathematical tools. In further investigations of the role of the calcium ion in regulating circadian period we found that calcium is sensed by a Calmodulin like 24, a protein of unknown function, and this appears to regulate the evening circadian clock protein TOC1. Lastly we have extended our analysis of the circadian regulation of calcium signalling, to show that the circadian clock regulates calcium signals also in the chloroplast..

We have further developed our mathematical tools, and we are applying those to circadian regulation of genes in Marchantia.
Exploitation Route We hope our new tool sets will be used by others to identify how drugs and genetic differences affect complex behaviours.
Sectors Agriculture, Food and Drink,Pharmaceuticals and Medical Biotechnology

 
Description We have produced a new tool (DYDE) for producing dynamic gene network models from time series transcriptional data. We have used this tool to make models of the plant circadian oscillator in Arabidopsis, and for the first time in Barley. This is the first mathematical model of a cereal circadian clock. Using DYDE we were able to describe the effects of chemical perturbations on the circadian oscillator. We also used this model to examine the regulation of the barley circadian transcriptome in day and night cycles, and describe for the first time the structure of the wheat circadian clock. The new model of the barley circadian oscillator and the DYDE tool allowed us to define a different structure of the circadian clock in cereals and the effect of genetic lesions on the circadian system. These data will provide an important tool for breeders to understand the impact of breeding on the circadian system
First Year Of Impact 2020
Sector Agriculture, Food and Drink
 
Title We developed a new method for measuring where in a biological network systems change occurs 
Description We have developed a method and reported code that allows time series data to be analysed to create a predicative network and allows for the way that network changes under stress to be determined. 
Type Of Material Technology assay or reagent 
Year Produced 2019 
Provided To Others? Yes  
Impact Using this tool we were able to identify how metabolites regulates circadian period. 
URL https://journals.plos.org/ploscompbiol/article?rev=2&id=10.1371/journal.pcbi.1006674
 
Title DyDE Dynamic determination of expression 
Description We have developed a tool to determine the structure of biological networks based on analysis and modelling of the dynamic control of gene expression. 
Type Of Material Computer model/algorithm 
Year Produced 2017 
Provided To Others? No  
Impact We have a paper in submission at the moment 
 
Title Models and Data for Mombaerts et al., "Dynamical Differential Expression (DyDE) Reveals the Period Control Mechanisms of the Arabidopsis Circadian Oscillator", PLOS Computational Biology, 2019. 
Description This dataset consists of experimental data that reveals the period control mechanisms of the Arabidopsis Circadian Oscillator (as described in the associated publication) and is provided here in an XLSX file. MATLAB code for performing the dynamic differential expression (DyDe) described in the publication are provided as .ZIP 
Type Of Material Database/Collection of data 
Year Produced 2018 
Provided To Others? Yes  
 
Description Collaboration with University of Luxembourg 
Organisation University of Luxembourg
Country Luxembourg 
Sector Academic/University 
PI Contribution We provide experimental data and insight in to the regulation of the circadian oscillator in plants and biological interpretation of the data
Collaborator Contribution Professor Jorge Goncalves and his team are providing experitise in linear mathematical modelling
Impact A research student in each of the Webb and Goncalves lab has received support and training from the PIs and PDRA.
Start Year 2017
 
Description Collaboration with the Max Planck Institute Cologne 
Organisation Max Planck Society
Department Max Planck Institute for Plant Breeding Research
Country Germany 
Sector Academic/University 
PI Contribution We are providing biological insight in to circadian systems and linear modelling toolkits
Collaborator Contribution The partners are providing data about the circadian regulation of barley genes
Impact Manuscript in preparation
Start Year 2015
 
Description Participation in Science Sunday 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach Local
Primary Audience Public/other audiences
Results and Impact The laboratory prepared an interactive display that described circadian rhythms for the Science Open Day. This usually attracts around 2000 people to the Plant Sciences display area.
Year(s) Of Engagement Activity 2016
 
Description Science Saturday 2016 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach Regional
Primary Audience Public/other audiences
Results and Impact Three PDRA supported by the BBSRC will develop a display for science Saturday concerning circadian biology
Year(s) Of Engagement Activity 2016