-
The ACROBAT 2022 Challenge: Automatic Registration Of Breast Cancer Tissue
Authors:
Philippe Weitz,
Masi Valkonen,
Leslie Solorzano,
Circe Carr,
Kimmo Kartasalo,
Constance Boissin,
Sonja Koivukoski,
Aino Kuusela,
Dusan Rasic,
Yanbo Feng,
Sandra Sinius Pouplier,
Abhinav Sharma,
Kajsa Ledesma Eriksson,
Stephanie Robertson,
Christian Marzahl,
Chandler D. Gatenbee,
Alexander R. A. Anderson,
Marek Wodzinski,
Artur Jurgas,
Niccolò Marini,
Manfredo Atzori,
Henning Müller,
Daniel Budelmann,
Nick Weiss,
Stefan Heldmann
, et al. (16 additional authors not shown)
Abstract:
The alignment of tissue between histopathological whole-slide-images (WSI) is crucial for research and clinical applications. Advances in computing, deep learning, and availability of large WSI datasets have revolutionised WSI analysis. Therefore, the current state-of-the-art in WSI registration is unclear. To address this, we conducted the ACROBAT challenge, based on the largest WSI registration…
▽ More
The alignment of tissue between histopathological whole-slide-images (WSI) is crucial for research and clinical applications. Advances in computing, deep learning, and availability of large WSI datasets have revolutionised WSI analysis. Therefore, the current state-of-the-art in WSI registration is unclear. To address this, we conducted the ACROBAT challenge, based on the largest WSI registration dataset to date, including 4,212 WSIs from 1,152 breast cancer patients. The challenge objective was to align WSIs of tissue that was stained with routine diagnostic immunohistochemistry to its H&E-stained counterpart. We compare the performance of eight WSI registration algorithms, including an investigation of the impact of different WSI properties and clinical covariates. We find that conceptually distinct WSI registration methods can lead to highly accurate registration performances and identify covariates that impact performances across methods. These results establish the current state-of-the-art in WSI registration and guide researchers in selecting and developing methods.
△ Less
Submitted 29 May, 2023;
originally announced May 2023.
-
A survey of open questions in adaptive therapy: bridging mathematics and clinical translation
Authors:
Jeffrey West,
Fred Adler,
Jill Gallaher,
Maximilian Strobl,
Renee Brady-Nicholls,
Joel S. Brown,
Mark Robertson-Tessi,
Eunjung Kim,
Robert Noble,
Yannick Viossat,
David Basanta,
Alexander R. A. Anderson
Abstract:
Adaptive therapy is a dynamic cancer treatment protocol that updates (or "adapts") treatment decisions in anticipation of evolving tumor dynamics. This broad term encompasses many possible dynamic treatment protocols of patient-specific dose modulation or dose timing. Adaptive therapy maintains high levels of tumor burden to benefit from the competitive suppression of treatment-sensitive subpopula…
▽ More
Adaptive therapy is a dynamic cancer treatment protocol that updates (or "adapts") treatment decisions in anticipation of evolving tumor dynamics. This broad term encompasses many possible dynamic treatment protocols of patient-specific dose modulation or dose timing. Adaptive therapy maintains high levels of tumor burden to benefit from the competitive suppression of treatment-sensitive subpopulations on treatment-resistant subpopulations. This evolution-based approach to cancer treatment has been integrated into several ongoing or planned clinical trials, including treatment of metastatic castrate resistant prostate cancer, ovarian cancer, and BRAF-mutant melanoma. In the previous few decades, experimental and clinical investigation of adaptive therapy has progressed synergistically with mathematical and computational modeling. In this work, we discuss 11 open questions in cancer adaptive therapy mathematical modeling. The questions are split into three sections: 1) the necessary components of mathematical models of adaptive therapy 2) design and validation of dosing protocols, and 3) challenges and opportunities in clinical translation.
△ Less
Submitted 18 October, 2022;
originally announced October 2022.
-
Local contact inhibition leads to universal principles of cell population growth
Authors:
Gregory J. Kimmel,
Jeffrey West,
Mehdi Damaghi,
Alexander R. A. Anderson,
Philipp M. Altrock
Abstract:
Cancer cell population dynamics often exhibit remarkably replicable, universal laws despite their underlying heterogeneity. Mechanistic explanations of universal cell population growth remain partly unresolved to this day, whereby population feedback between the microscopic and mesoscopic configurations can lead to macroscopic growth laws. We here present a unification under density-dependent birt…
▽ More
Cancer cell population dynamics often exhibit remarkably replicable, universal laws despite their underlying heterogeneity. Mechanistic explanations of universal cell population growth remain partly unresolved to this day, whereby population feedback between the microscopic and mesoscopic configurations can lead to macroscopic growth laws. We here present a unification under density-dependent birth events via contact inhibition. We consider five classical tumor growth laws: exponential, generalized logistic, Gompertz, radial growth, and fractal growth, which can be seen as manifestations of a single microscopic model. Our theory is substantiated by agent based simulations and can explain growth curve differences in experimental data from in vitro cancer cell population growth. Thus, our framework offers a possible explanation for the large number of mean-field laws that can adequately capture seemingly unrelated cancer or microbial growth dynamics.
△ Less
Submitted 23 August, 2021;
originally announced August 2021.
-
A comparative study between discrete and continuum models for the evolution of competing phenotype-structured cell populations in dynamical environments
Authors:
Aleksandra Ardaševa,
Robert A. Gatenby,
Alexander R. A. Anderson,
Helen M. Byrne,
Philip K. Maini,
Tommaso Lorenzi
Abstract:
Deterministic continuum models formulated in terms of non-local partial differential equations for the evolutionary dynamics of populations structured by phenotypic traits have been used recently to address open questions concerning the adaptation of asexual species to periodically fluctuating environmental conditions. These deterministic continuum models are usually defined on the basis of popula…
▽ More
Deterministic continuum models formulated in terms of non-local partial differential equations for the evolutionary dynamics of populations structured by phenotypic traits have been used recently to address open questions concerning the adaptation of asexual species to periodically fluctuating environmental conditions. These deterministic continuum models are usually defined on the basis of population-scale phenomenological assumptions and cannot capture adaptive phenomena that are driven by stochastic variability in the evolutionary paths of single individuals. In this paper, we develop a stochastic individual-based model for the coevolution between two competing phenotype-structured cell populations that are exposed to time-varying nutrient levels and undergo spontaneous, heritable phenotypic variations with different probabilities. The evolution of every cell is described by a set of rules that result in a discrete-time branching random walk on the space of phenotypic states. We formally show that the deterministic continuum counterpart of this model comprises a system of non-local partial differential equations for the cell population density functions coupled with an ordinary differential equation for the nutrient concentration. We compare the individual-based model and its continuum analogue, focussing on scenarios whereby the predictions of the two models differ. Our results clarify the conditions under which significant differences between the two models can emerge due to stochastic effects associated with small population levels. These differences arise in the presence of low probabilities of phenotypic variation, and become more apparent when the two populations are characterised by less fit initial mean phenotypes and smaller initial levels of phenotypic heterogeneity.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
Evolutionary dynamics of competing phenotype-structured populations in periodically fluctuating environments
Authors:
Aleksandra Ardaševa,
Robert A. Gatenby,
Alexander R. A. Anderson,
Helen M. Byrne,
Philip K. Maini,
Tommaso Lorenzi
Abstract:
Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to s…
▽ More
Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to support a deeper understanding of the adaptive role of spontaneous phenotypic variations in fluctuating environments, we consider a system of non-local partial differential equations modelling the evolutionary dynamics of two competing phenotype-structured populations in the presence of periodically oscillating nutrient levels. The two populations undergo spontaneous phenotypic variations at different rates. The phenotypic state of each individual is represented by a continuous variable, and the phenotypic landscape of the populations evolves in time due to variations in the nutrient level. Exploiting the analytical tractability of our model, we study the long-time behaviour of the solutions to obtain a detailed mathematical depiction of evolutionary dynamics. The results suggest that when nutrient levels undergo small and slow oscillations, it is evolutionarily more convenient to rarely undergo spontaneous phenotypic variations. Conversely, under relatively large and fast periodic oscillations in the nutrient levels, which bring about alternating cycles of starvation and nutrient abundance, higher rates of spontaneous phenotypic variations confer a competitive advantage. We discuss the implications of our results in the context of cancer metabolism.
△ Less
Submitted 23 August, 2019; v1 submitted 28 May, 2019;
originally announced May 2019.
-
Cardiac re-entry dynamics & self-termination in DT-MRI based model of Human Foetal Heart
Authors:
Irina V. Biktasheva,
Richard A. Anderson,
Arun V. Holden,
Eleftheria Pervolaraki,
Fengcai Wen
Abstract:
The effect of heart geometry and anisotropy on cardiac re-entry dynamics and self-termination is studied here in anatomically realistic computer simulations of human foetal heart. 20 weeks of gestational age human foetal heart isotropic and anisotropic anatomy models from diffusion tensor MRI data sets are used in the computer simulations. The fibre orientation angles of the heart were obtained fr…
▽ More
The effect of heart geometry and anisotropy on cardiac re-entry dynamics and self-termination is studied here in anatomically realistic computer simulations of human foetal heart. 20 weeks of gestational age human foetal heart isotropic and anisotropic anatomy models from diffusion tensor MRI data sets are used in the computer simulations. The fibre orientation angles of the heart were obtained from the DT-MRI primary eigenvalues. In a spatially homogeneous electrophysiological mono domain model with the DT-MRI based heart geometries, we initiate simplified Fitz-Hugh-Nagumo kinetics cardiac re-entry at a prescribed location in a 2D slice, and in the full 3D anatomy model. In a slice of the heart, the MRI based fibre anisotropy changes the re-entry dynamics from pinned to anatomical re-entry. In the full 3D MRI based model, the foetal heart fibre anisotropy changes the re-entry dynamics from a persistent re-entry to the re-entry self-termination.
△ Less
Submitted 19 August, 2017; v1 submitted 13 April, 2017;
originally announced April 2017.
-
The Sensitivity of the Advanced LIGO Detectors at the Beginning of Gravitational Wave Astronomy
Authors:
D. V. Martynov,
E. D. Hall,
B. P. Abbott,
R. Abbott,
T. D. Abbott,
C. Adams,
R. X. Adhikari,
R. A. Anderson,
S. B. Anderson,
K. Arai,
M. A. Arain,
S. M. Aston,
L. Austin,
S. W. Ballmer,
M. Barbet,
D. Barker,
B. Barr,
L. Barsotti,
J. Bartlett,
M. A. Barton,
I. Bartos,
J. C. Batch,
A. S. Bell,
I. Belopolski,
J. Bergman
, et al. (239 additional authors not shown)
Abstract:
The Laser Interferometer Gravitational Wave Observatory (LIGO) consists of two widely separated 4 km laser interferometers designed to detect gravitational waves from distant astrophysical sources in the frequency range from 10 Hz to 10 kHz. The first observation run of the Advanced LIGO detectors started in September 2015 and ended in January 2016. A strain sensitivity of better than…
▽ More
The Laser Interferometer Gravitational Wave Observatory (LIGO) consists of two widely separated 4 km laser interferometers designed to detect gravitational waves from distant astrophysical sources in the frequency range from 10 Hz to 10 kHz. The first observation run of the Advanced LIGO detectors started in September 2015 and ended in January 2016. A strain sensitivity of better than $10^{-23}/\sqrt{\text{Hz}}$ was achieved around 100 Hz. Understanding both the fundamental and the technical noise sources was critical for increasing the observable volume in the universe. The average distance at which coalescing binary black hole systems with individual masses of 30 $M_\odot$ could be detected was 1.3 Gpc. Similarly, the range for binary neutron star inspirals was about 75 Mpc. With respect to the initial detectors, the observable volume of Universe increased respectively by a factor 69 and 43. These improvements allowed Advanced LIGO to detect the gravitational wave signal from the binary black hole coalescence, known as GW150914.
△ Less
Submitted 10 February, 2018; v1 submitted 1 April, 2016;
originally announced April 2016.
-
The impact of cellular characteristics on the evolution of shape homeostasis
Authors:
Philip Gerlee,
David Basanta,
Alexander R. A. Anderson
Abstract:
The importance of individual cells in a developing multicellular organism is well known but precisely how the individual cellular characteristics of those cells collectively drive the emergence of robust, homeostatic structures is less well understood. For example cell communication via a diffusible factor allows for information to travel across large distances within the population, and cell pola…
▽ More
The importance of individual cells in a developing multicellular organism is well known but precisely how the individual cellular characteristics of those cells collectively drive the emergence of robust, homeostatic structures is less well understood. For example cell communication via a diffusible factor allows for information to travel across large distances within the population, and cell polarisation makes it possible to form structures with a particular orientation, but how do these processes interact to produce a more robust and regulated structure? In this study we investigate the ability of cells with different cellular characteristics to grow and maintain homeostatic structures. We do this in the context of an individual-based model where cell behaviour is driven by an intra-cellular network that determines the cell phenotype. More precisely, we investigated evolution with 96 different permutations of our model, where cell motility, cell death, long-range growth factor (LGF), short-range growth factor (SGF) and cell polarisation were either present or absent. The results show that LGF has the largest positive impact on the fitness of the evolved solutions. SGF and polarisation also contribute, but all other capabilities essentially increase the search space, effectively making it more difficult to achieve a solution. By perturbing the evolved solutions, we found that they are highly robust to both mutations and wounding. In addition, we observed that by evolving solutions in more unstable environments they produce structures that were more robust and adaptive. In conclusion, our results suggest that robust collective behaviour is most likely to evolve when cells are endowed with long range communication, cell polarisation, and selection pressure from an unstable environment.
△ Less
Submitted 17 February, 2017; v1 submitted 8 December, 2015;
originally announced December 2015.
-
Characterization of the LIGO detectors during their sixth science run
Authors:
The LIGO Scientific Collaboration,
The Virgo Collaboration,
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amador. Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson,
W. G. Anderson
, et al. (846 additional authors not shown)
Abstract:
In 2009-2010, the Laser Interferometer Gravitational-wave Observa- tory (LIGO) operated together with international partners Virgo and GEO600 as a network to search for gravitational waves of astrophysical origin. The sensitiv- ity of these detectors was limited by a combination of noise sources inherent to the instrumental design and its environment, often localized in time or frequency, that cou…
▽ More
In 2009-2010, the Laser Interferometer Gravitational-wave Observa- tory (LIGO) operated together with international partners Virgo and GEO600 as a network to search for gravitational waves of astrophysical origin. The sensitiv- ity of these detectors was limited by a combination of noise sources inherent to the instrumental design and its environment, often localized in time or frequency, that couple into the gravitational-wave readout. Here we review the performance of the LIGO instruments during this epoch, the work done to characterize the de- tectors and their data, and the effect that transient and continuous noise artefacts have on the sensitivity of LIGO to a variety of astrophysical sources.
△ Less
Submitted 18 November, 2014; v1 submitted 28 October, 2014;
originally announced October 2014.
-
Searching for stochastic gravitational waves using data from the two co-located LIGO Hanford detectors
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amado. Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson
, et al. (852 additional authors not shown)
Abstract:
Searches for a stochastic gravitational-wave background (SGWB) using terrestrial detectors typically involve cross-correlating data from pairs of detectors. The sensitivity of such cross-correlation analyses depends, among other things, on the separation between the two detectors: the smaller the separation, the better the sensitivity. Hence, a co-located detector pair is more sensitive to a gravi…
▽ More
Searches for a stochastic gravitational-wave background (SGWB) using terrestrial detectors typically involve cross-correlating data from pairs of detectors. The sensitivity of such cross-correlation analyses depends, among other things, on the separation between the two detectors: the smaller the separation, the better the sensitivity. Hence, a co-located detector pair is more sensitive to a gravitational-wave background than a non-co-located detector pair. However, co-located detectors are also expected to suffer from correlated noise from instrumental and environmental effects that could contaminate the measurement of the background. Hence, methods to identify and mitigate the effects of correlated noise are necessary to achieve the potential increase in sensitivity of co-located detectors. Here we report on the first SGWB analysis using the two LIGO Hanford detectors and address the complications arising from correlated environmental noise. We apply correlated noise identification and mitigation techniques to data taken by the two LIGO Hanford detectors, H1 and H2, during LIGO's fifth science run. At low frequencies, 40 - 460 Hz, we are unable to sufficiently mitigate the correlated noise to a level where we may confidently measure or bound the stochastic gravitational-wave signal. However, at high frequencies, 460-1000 Hz, these techniques are sufficient to set a $95%$ confidence level (C.L.) upper limit on the gravitational-wave energy density of Ω(f)<7.7 x 10^{-4} (f/ 900 Hz)^3, which improves on the previous upper limit by a factor of $\sim 180$. In doing so, we demonstrate techniques that will be useful for future searches using advanced detectors, where correlated noise (e.g., from global magnetic fields) may affect even widely separated detectors.
△ Less
Submitted 2 December, 2014; v1 submitted 22 October, 2014;
originally announced October 2014.
-
Multimessenger Search for Sources of Gravitational Waves and High-Energy Neutrinos: Results for Initial LIGO-Virgo and IceCube
Authors:
The IceCube Collaboration,
The LIGO Scientific Collaboration,
The Virgo Collaboration,
M. G. Aartsen,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
D. Altmann,
T. Anderson,
C. Arguelles,
T. C. Arlen,
J. Auffenberg,
X. Bai,
S. W. Barwick,
V. Baum,
J. J. Beatty,
J. Becker Tjus,
K. -H. Becker,
S. BenZvi,
P. Berghaus,
D. Berley,
E. Bernardini,
A. Bernhard
, et al. (1166 additional authors not shown)
Abstract:
We report the results of a multimessenger search for coincident signals from the LIGO and Virgo gravitational-wave observatories and the partially completed IceCube high-energy neutrino detector, including periods of joint operation between 2007-2010. These include parts of the 2005-2007 run and the 2009-2010 run for LIGO-Virgo, and IceCube's observation periods with 22, 59 and 79 strings. We find…
▽ More
We report the results of a multimessenger search for coincident signals from the LIGO and Virgo gravitational-wave observatories and the partially completed IceCube high-energy neutrino detector, including periods of joint operation between 2007-2010. These include parts of the 2005-2007 run and the 2009-2010 run for LIGO-Virgo, and IceCube's observation periods with 22, 59 and 79 strings. We find no significant coincident events, and use the search results to derive upper limits on the rate of joint sources for a range of source emission parameters. For the optimistic assumption of gravitational-wave emission energy of $10^{-2}$\,M$_\odot$c$^2$ at $\sim 150$\,Hz with $\sim 60$\,ms duration, and high-energy neutrino emission of $10^{51}$\,erg comparable to the isotropic gamma-ray energy of gamma-ray bursts, we limit the source rate below $1.6 \times 10^{-2}$\,Mpc$^{-3}$yr$^{-1}$. We also examine how combining information from gravitational waves and neutrinos will aid discovery in the advanced gravitational-wave detector era.
△ Less
Submitted 9 October, 2014; v1 submitted 3 July, 2014;
originally announced July 2014.
-
Methods and results of a search for gravitational waves associated with gamma-ray bursts using the GEO600, LIGO, and Virgo detectors
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen,
R. A. Anderson,
S. B. Anderson
, et al. (868 additional authors not shown)
Abstract:
In this paper we report on a search for short-duration gravitational wave bursts in the frequency range 64 Hz-1792 Hz associated with gamma-ray bursts (GRBs), using data from GEO600 and one of the LIGO or Virgo detectors. We introduce the method of a linear search grid to analyse GRB events with large sky localisation uncertainties such as the localisations provided by the Fermi Gamma-ray Burst Mo…
▽ More
In this paper we report on a search for short-duration gravitational wave bursts in the frequency range 64 Hz-1792 Hz associated with gamma-ray bursts (GRBs), using data from GEO600 and one of the LIGO or Virgo detectors. We introduce the method of a linear search grid to analyse GRB events with large sky localisation uncertainties such as the localisations provided by the Fermi Gamma-ray Burst Monitor (GBM). Coherent searches for gravitational waves (GWs) can be computationally intensive when the GRB sky position is not well-localised, due to the corrections required for the difference in arrival time between detectors. Using a linear search grid we are able to reduce the computational cost of the analysis by a factor of O(10) for GBM events. Furthermore, we demonstrate that our analysis pipeline can improve upon the sky localisation of GRBs detected by the GBM, if a high-frequency GW signal is observed in coincidence. We use the linear search grid method in a search for GWs associated with 129 GRBs observed satellite-based gamma-ray experiments between 2006 and 2011. The GRBs in our sample had not been previously analysed for GW counterparts. A fraction of our GRB events are analysed using data from GEO600 while the detector was using squeezed-light states to improve its sensitivity; this is the first search for GWs using data from a squeezed-light interferometric observatory. We find no evidence for GW signals, either with any individual GRB in this sample or with the population as a whole. For each GRB we place lower bounds on the distance to the progenitor, assuming a fixed GW emission energy of $10^{-2} M_{\odot}c^{2}$, with a median exclusion distance of 0.8 Mpc for emission at 500 Hz and 0.3 Mpc at 1 kHz. The reduced computational cost associated with a linear search grid will enable rapid searches for GWs associated with Fermi GBM events in the Advanced detector era.
△ Less
Submitted 1 July, 2014; v1 submitted 5 May, 2014;
originally announced May 2014.
-
Bridging scales in cancer progression: Mapping genotype to phenotype using neural networks
Authors:
Philip Gerlee,
Eunjung Kim,
Alexander R. A. Anderson
Abstract:
In this review we summarize our recent efforts in trying to understand the role of heterogeneity in cancer progression by using neural networks to characterise different aspects of the mapping from a cancer cells genotype and environment to its phenotype. Our central premise is that cancer is an evolving system subject to mutation and selection, and the primary conduit for these processes to occur…
▽ More
In this review we summarize our recent efforts in trying to understand the role of heterogeneity in cancer progression by using neural networks to characterise different aspects of the mapping from a cancer cells genotype and environment to its phenotype. Our central premise is that cancer is an evolving system subject to mutation and selection, and the primary conduit for these processes to occur is the cancer cell whose behaviour is regulated on multiple biological scales. The selection pressure is mainly driven by the microenvironment that the tumour is growing in and this acts directly upon the cell phenotype. In turn, the phenotype is driven by the intracellular pathways that are regulated by the genotype. Integrating all of these processes is a massive undertaking and requires bridging many biological scales (i.e. genotype, pathway, phenotype and environment) that we will only scratch the surface of in this review. We will focus on models that use neural networks as a means of connecting these different biological scales, since they allow us to easily create heterogeneity for selection to act upon and importantly this heterogeneity can be implemented at different biological scales. More specifically, we consider three different neural networks that bridge different aspects of these scales and the dialogue with the micro-environment, (i) the impact of the micro-environment on evolutionary dynamics, (ii) the mapping from genotype to phenotype under drug-induced perturbations and (iii) pathway activity in both normal and cancer cells under different micro-environmental conditions.
△ Less
Submitted 28 April, 2014;
originally announced April 2014.
-
Search for gravitational waves associated with gamma-ray bursts detected by the InterPlanetary Network
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
F. Acernese,
K. Ackley,
C. Adams,
T. Adams,
P. Addesso,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
A. Alemic,
B. Allen,
A. Allocca,
D. Amariutei,
M. Andersen,
R. A. Anderson,
S. B. Anderson
, et al. (879 additional authors not shown)
Abstract:
We present the results of a search for gravitational waves associated with 223 gamma-ray bursts (GRBs) detected by the InterPlanetary Network (IPN) in 2005-2010 during LIGO's fifth and sixth science runs and Virgo's first, second and third science runs. The IPN satellites provide accurate times of the bursts and sky localizations that vary significantly from degree scale to hundreds of square degr…
▽ More
We present the results of a search for gravitational waves associated with 223 gamma-ray bursts (GRBs) detected by the InterPlanetary Network (IPN) in 2005-2010 during LIGO's fifth and sixth science runs and Virgo's first, second and third science runs. The IPN satellites provide accurate times of the bursts and sky localizations that vary significantly from degree scale to hundreds of square degrees. We search for both a well-modeled binary coalescence signal, the favored progenitor model for short GRBs, and for generic, unmodeled gravitational wave bursts. Both searches use the event time and sky localization to improve the gravitational-wave search sensitivity as compared to corresponding all-time, all-sky searches. We find no evidence of a gravitational-wave signal associated with any of the IPN GRBs in the sample, nor do we find evidence for a population of weak gravitational-wave signals associated with the GRBs. For all IPN-detected GRBs, for which a sufficient duration of quality gravitational-wave data is available, we place lower bounds on the distance to the source in accordance with an optimistic assumption of gravitational-wave emission energy of $10^{-2}M_{\odot}c^2$ at 150 Hz, and find a median of 13 Mpc. For the 27 short-hard GRBs we place 90% confidence exclusion distances to two source models: a binary neutron star coalescence, with a median distance of 12Mpc, or the coalescence of a neutron star and black hole, with a median distance of 22 Mpc. Finally, we combine this search with previously published results to provide a population statement for GRB searches in first-generation LIGO and Virgo gravitational-wave detectors, and a resulting examination of prospects for the advanced gravitational-wave detectors.
△ Less
Submitted 17 April, 2014; v1 submitted 26 March, 2014;
originally announced March 2014.
-
The evolution of carrying capacity in constrained and expanding tumour cell populations
Authors:
Philip Gerlee,
Alexander R. A. Anderson
Abstract:
Cancer cells are known to modify their micro-environment such that it can sustain a larger population, or, in ecological terms, they construct a niche which increases the carrying capacity of the population. It has however been argued that niche construction, which benefits all cells in the tumour, would be selected against since cheaters could reap the benefits without paying the cost. We have in…
▽ More
Cancer cells are known to modify their micro-environment such that it can sustain a larger population, or, in ecological terms, they construct a niche which increases the carrying capacity of the population. It has however been argued that niche construction, which benefits all cells in the tumour, would be selected against since cheaters could reap the benefits without paying the cost. We have investigated the impact of niche specificity on tumour evolution using an individual based model of breast tumour growth, in which the carrying capacity of each cell consists of two components: an intrinsic, subclone-specific part and a contribution from all neighbouring cells. Analysis of the model shows that the ability of a mutant to invade a resident population depends strongly on the specificity. When specificity is low selection is mostly on growth rate, while high specificity shifts selection towards increased carrying capacity. Further, we show that the long-term evolution of the system can be predicted using adaptive dynamics. By comparing the results from a spatially structured vs.\ well-mixed population we show that spatial structure restores selection for carrying capacity even at zero specificity, which a poses solution to the niche construction dilemma. Lastly, we show that an expanding population exhibits spatially variable selection pressure, where cells at the leading edge exhibit higher growth rate and lower carrying capacity than those at the centre of the tumour.
△ Less
Submitted 14 August, 2015; v1 submitted 4 February, 2014;
originally announced February 2014.
-
Application of a Hough search for continuous gravitational waves on data from the 5th LIGO science run
Authors:
The LIGO Scientific Collaboration,
The Virgo Collaboration,
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amador Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson,
W. G. Anderson
, et al. (855 additional authors not shown)
Abstract:
We report on an all-sky search for periodic gravitational waves in the frequency range $\mathrm{50-1000 Hz}$ with the first derivative of frequency in the range $-8.9 \times 10^{-10}$ Hz/s to zero in two years of data collected during LIGO's fifth science run. Our results employ a Hough transform technique, introducing a $χ^2$ test and analysis of coincidences between the signal levels in years 1…
▽ More
We report on an all-sky search for periodic gravitational waves in the frequency range $\mathrm{50-1000 Hz}$ with the first derivative of frequency in the range $-8.9 \times 10^{-10}$ Hz/s to zero in two years of data collected during LIGO's fifth science run. Our results employ a Hough transform technique, introducing a $χ^2$ test and analysis of coincidences between the signal levels in years 1 and 2 of observations that offers a significant improvement in the product of strain sensitivity with compute cycles per data sample compared to previously published searches. Since our search yields no surviving candidates, we present results taking the form of frequency dependent, 95$%$ confidence upper limits on the strain amplitude $h_0$. The most stringent upper limit from year 1 is $1.0\times 10^{-24}$ in the $\mathrm{158.00-158.25 Hz}$ band. In year 2, the most stringent upper limit is $\mathrm{8.9\times10^{-25}}$ in the $\mathrm{146.50-146.75 Hz}$ band. This improved detection pipeline, which is computationally efficient by at least two orders of magnitude better than our flagship Einstein$@$Home search, will be important for "quick-look" searches in the Advanced LIGO and Virgo detector era.
△ Less
Submitted 17 March, 2014; v1 submitted 11 November, 2013;
originally announced November 2013.
-
Constraints on cosmic strings from the LIGO-Virgo gravitational-wave detectors
Authors:
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amador Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson,
W. G. Anderson,
K. Arai,
M. C. Araya
, et al. (852 additional authors not shown)
Abstract:
Cosmic strings can give rise to a large variety of interesting astrophysical phenomena. Among them, powerful bursts of gravitational waves (GWs) produced by cusps are a promising observational signature. In this Letter we present a search for GWs from cosmic string cusps in data collected by the LIGO and Virgo gravitational wave detectors between 2005 and 2010, with over 625 days of live time. We…
▽ More
Cosmic strings can give rise to a large variety of interesting astrophysical phenomena. Among them, powerful bursts of gravitational waves (GWs) produced by cusps are a promising observational signature. In this Letter we present a search for GWs from cosmic string cusps in data collected by the LIGO and Virgo gravitational wave detectors between 2005 and 2010, with over 625 days of live time. We find no evidence of GW signals from cosmic strings. From this result, we derive new constraints on cosmic string parameters, which complement and improve existing limits from previous searches for a stochastic background of GWs from cosmic microwave background measurements and pulsar timing data. In particular, if the size of loops is given by the gravitational backreaction scale, we place upper limits on the string tension $Gμ$ below $10^{-8}$ in some regions of the cosmic string parameter space.
△ Less
Submitted 7 April, 2014; v1 submitted 9 October, 2013;
originally announced October 2013.
-
First Searches for Optical Counterparts to Gravitational-wave Candidate Events
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amador Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson,
W. G. Anderson
, et al. (883 additional authors not shown)
Abstract:
During the LIGO and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight su…
▽ More
During the LIGO and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight such GW candidates. We present the methods used to analyze the image data as well as the transient search results. No optical transient was identified with a convincing association with any of these candidates, and none of the GW triggers showed strong evidence for being astrophysical in nature. We compare the sensitivities of these observations to several model light curves from possible sources of interest, and discuss prospects for future joint GW-optical observations of this type.
△ Less
Submitted 21 October, 2013; v1 submitted 8 October, 2013;
originally announced October 2013.
-
A directed search for continuous Gravitational Waves from the Galactic Center
Authors:
The LIGO Scientific Collaboration,
The Virgo Collaboration,
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amador Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson,
W. G. Anderson
, et al. (850 additional authors not shown)
Abstract:
We present the results of a directed search for continuous gravitational waves from unknown, isolated neutron stars in the Galactic Center region, performed on two years of data from LIGO's fifth science run from two LIGO detectors. The search uses a semi-coherent approach, analyzing coherently 630 segments, each spanning 11.5 hours, and then incoherently combining the results of the single segmen…
▽ More
We present the results of a directed search for continuous gravitational waves from unknown, isolated neutron stars in the Galactic Center region, performed on two years of data from LIGO's fifth science run from two LIGO detectors. The search uses a semi-coherent approach, analyzing coherently 630 segments, each spanning 11.5 hours, and then incoherently combining the results of the single segments. It covers gravitational wave frequencies in a range from 78 to 496 Hz and a frequency-dependent range of first order spindown values down to -7.86 x 10^-8 Hz/s at the highest frequency. No gravitational waves were detected. We place 90% confidence upper limits on the gravitational wave amplitude of sources at the Galactic Center. Placing 90% confidence upper limits on the gravitational wave amplitude of sources at the Galactic Center, we reach ~3.35x10^-25 for frequencies near 150 Hz. These upper limits are the most constraining to date for a large-parameter-space search for continuous gravitational wave signals.
△ Less
Submitted 27 September, 2013; v1 submitted 24 September, 2013;
originally announced September 2013.
-
Search for long-lived gravitational-wave transients coincident with long gamma-ray bursts
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amador Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson,
W. G. Anderson
, et al. (854 additional authors not shown)
Abstract:
Long gamma-ray bursts (GRBs) have been linked to extreme core-collapse supernovae from massive stars. Gravitational waves (GW) offer a probe of the physics behind long GRBs. We investigate models of long-lived (~10-1000s) GW emission associated with the accretion disk of a collapsed star or with its protoneutron star remnant. Using data from LIGO's fifth science run, and GRB triggers from the swif…
▽ More
Long gamma-ray bursts (GRBs) have been linked to extreme core-collapse supernovae from massive stars. Gravitational waves (GW) offer a probe of the physics behind long GRBs. We investigate models of long-lived (~10-1000s) GW emission associated with the accretion disk of a collapsed star or with its protoneutron star remnant. Using data from LIGO's fifth science run, and GRB triggers from the swift experiment, we perform a search for unmodeled long-lived GW transients. Finding no evidence of GW emission, we place 90% confidence level upper limits on the GW fluence at Earth from long GRBs for three waveforms inspired by a model of GWs from accretion disk instabilities. These limits range from F<3.5 ergs cm^-2 to $F<1200 ergs cm^-2, depending on the GRB and on the model, allowing us to probe optimistic scenarios of GW production out to distances as far as ~33 Mpc. Advanced detectors are expected to achieve strain sensitivities 10x better than initial LIGO, potentially allowing us to probe the engines of the nearest long GRBs.
△ Less
Submitted 9 December, 2013; v1 submitted 24 September, 2013;
originally announced September 2013.
-
A filter-flow perspective of hematogenous metastasis offers a non-genetic paradigm for personalized cancer therapy
Authors:
Jacob G. Scott,
Alexander G. Fletcher,
Philip K. Maini,
Alexander R. A. Anderson,
Philip Gerlee
Abstract:
Research into mechanisms of hematogenous metastasis has largely become genetic in focus, attempting to understand the molecular basis of `seed-soil' relationships. Preceeding this biological mechanism is the physical process of dissemination of circulating tumour cells (CTCs). We utilize a `filter-flow' paradigm to show that assumptions about CTC dynamics strongly affect metastatic efficiency: wit…
▽ More
Research into mechanisms of hematogenous metastasis has largely become genetic in focus, attempting to understand the molecular basis of `seed-soil' relationships. Preceeding this biological mechanism is the physical process of dissemination of circulating tumour cells (CTCs). We utilize a `filter-flow' paradigm to show that assumptions about CTC dynamics strongly affect metastatic efficiency: without data on CTC dynamics, any attempt to predict metastatic spread in individual patients is impossible.
△ Less
Submitted 19 September, 2013;
originally announced September 2013.
-
Gravitational waves from known pulsars: results from the initial detector era
Authors:
J. Aasi,
J. Abadie,
B. P. Abbott,
R. Abbott,
T. Abbott,
M. R. Abernathy,
T. Accadia,
F. Acernese,
C. Adams,
T. Adams,
R. X. Adhikari,
C. Affeldt,
M. Agathos,
N. Aggarwal,
O. D. Aguiar,
P. Ajith,
B. Allen,
A. Allocca,
E. Amador Ceron,
D. Amariutei,
R. A. Anderson,
S. B. Anderson,
W. G. Anderson,
K. Arai,
M. C. Araya
, et al. (871 additional authors not shown)
Abstract:
We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produ…
▽ More
We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.
△ Less
Submitted 24 April, 2014; v1 submitted 16 September, 2013;
originally announced September 2013.
-
Exploiting ecological principles to better understand cancer progression and treatment
Authors:
David Basanta,
Alexander R. A. Anderson
Abstract:
A small but growing number of people are finding interesting parallels between ecosystems as studied by ecologists (think of a Savanna or the Amazon rain forest or a Coral reef) and tumours1-3. The idea of viewing cancer from an ecological perspective has many implications but fundamentally, it means that we should not see cancer just as a group of mutated cells. A more useful definition of cancer…
▽ More
A small but growing number of people are finding interesting parallels between ecosystems as studied by ecologists (think of a Savanna or the Amazon rain forest or a Coral reef) and tumours1-3. The idea of viewing cancer from an ecological perspective has many implications but fundamentally, it means that we should not see cancer just as a group of mutated cells. A more useful definition of cancer is to consider it a disruption in the complex balance of many interacting cellular and microenvironmental elements in a specific organ. This perspective means that organs undergoing carcinogenesis should be seen as sophisticated ecosystems in homeostasis that cancer cells can disrupt. It also makes cancer seem even more complex but may ultimately provides isights that make it more treatable. Here we discuss how ecological principles can be used to better understand cancer progression and treatment, using several mathematical and computational models to illustrate our argument.
△ Less
Submitted 9 May, 2013;
originally announced May 2013.
-
Evolution of intratumoral phenotypic heterogeneity: the role of trait inheritance
Authors:
Jill Gallaher,
Alexander R. A. Anderson
Abstract:
A tumor can be thought of as an ecosystem, which critically means that we cannot just consider it as a collection of mutated cells but more as a complex system of many interacting cellular and microenvironmental elements. At its simplest, a growing tumor with increased proliferation capacity must compete for space as a limited resource. Hypercellularity leads to a contact-inhibited core with a com…
▽ More
A tumor can be thought of as an ecosystem, which critically means that we cannot just consider it as a collection of mutated cells but more as a complex system of many interacting cellular and microenvironmental elements. At its simplest, a growing tumor with increased proliferation capacity must compete for space as a limited resource. Hypercellularity leads to a contact-inhibited core with a competitive proliferating rim. Evolution and selection occurs, and an individual cell's capacity to survive and propagate is determined by its combination of traits and interaction with the environment. With heterogeneity in phenotypes, the clone that will dominate is not always obvious as there are both local interactions and global pressures. Several combinations of phenotypes can coexist, changing the fitness of the whole.
To understand some aspects of heterogeneity in a growing tumor we build an off-lattice agent based model consisting of individual cells with assigned trait values for proliferation and migration rates. We represent heterogeneity in these traits with frequency distributions and combinations of traits with density maps. How the distributions change over time is dependent on how traits are passed on to progeny cells, which is our main inquiry. We bypass the translation of genetics to behavior by focussing on the functional end result of inheritance of the phenotype combined with the environmental influence of limited space.
△ Less
Submitted 2 May, 2013;
originally announced May 2013.
-
Ovarian volume throughout life: a validated normative model
Authors:
Thomas W. Kelsey,
Sarah K. Dodwell,
A. Graham Wilkinson,
Tine Greve,
Claus Y. Andersen,
Richard A. Anderson,
W. Hamish B. Wallace
Abstract:
The measurement of ovarian volume has been shown to be a useful indirect indicator of the ovarian reserve in women of reproductive age, in the diagnosis and management of a number of disorders of puberty and adult reproductive function, and is under investigation as a screening tool for ovarian cancer. To date there is no normative model of ovarian volume throughout life. By searching the publishe…
▽ More
The measurement of ovarian volume has been shown to be a useful indirect indicator of the ovarian reserve in women of reproductive age, in the diagnosis and management of a number of disorders of puberty and adult reproductive function, and is under investigation as a screening tool for ovarian cancer. To date there is no normative model of ovarian volume throughout life. By searching the published literature for ovarian volume in healthy females, and using our own data from multiple sources (combined n = 59,994) we have generated and robustly validated the first model of ovarian volume from conception to 82 years of age. This model shows that 69% of the variation in ovarian volume is due to age alone. We have shown that in the average case ovarian volume rises from 0.7 mL (95% CI 0.4 -- 1.1 mL) at 2 years of age to a peak of 7.7 mL (95% CI 6.5 -- 9.2 mL) at 20 years of age with a subsequent decline to about 2.8mL (95% CI 2.7 -- 2.9 mL) at the menopause and smaller volumes thereafter. Our model allows us to generate normal values and ranges for ovarian volume throughout life. This is the first validated normative model of ovarian volume from conception to old age; it will be of use in the diagnosis and management of a number of diverse gynaecological and reproductive conditions in females from birth to menopause and beyond.
△ Less
Submitted 23 April, 2013;
originally announced April 2013.
-
Senescent fibroblasts can drive melanoma initiation and progression
Authors:
Eunjung Kim,
Vito Rebecca,
Inna V. Fedorenko,
Jane L. Messina,
Rahel Mathew,
Silvya S. Maria-Engler,
David Basanta,
Keiran S. M. Smalley,
Alexander R. A. Anderson
Abstract:
Skin is one of the largest human organ systems whose primary purpose is the protection of deeper tissues. As such, the skin must maintain a homeostatic balance in the face of many microenvironmental and genetic perturbations. At its simplest, skin homeostasis is maintained by the balance between skin cell growth and death such that skin architecture is preserved. This study presents a hybrid multi…
▽ More
Skin is one of the largest human organ systems whose primary purpose is the protection of deeper tissues. As such, the skin must maintain a homeostatic balance in the face of many microenvironmental and genetic perturbations. At its simplest, skin homeostasis is maintained by the balance between skin cell growth and death such that skin architecture is preserved. This study presents a hybrid multiscale mathematical model of normal skin (vSkin). The model focuses on key cellular and microenvironmental variables that regulate homeostatic interactions among keratinocytes, melanocytes and fibroblasts, key components of the skin. The model recapitulates normal skin structure, and is robust enough to withstand physical as well as biochemical perturbations. Furthermore, the vSkin model revealed the important role of the skin microenvironment in melanoma initiation and progression. Our experiments showed that dermal fibroblasts, which are an important source of growth factors in the skin, adopt a phenotype that facilitates cancer cell growth and invasion when they become senescent. Based on these experimental results, we incorporated senescent fibroblasts into vSkin model and showed that senescent fibroblasts transform the skin microenvironment and enhance the growth and invasion of normal melanocytes as well as early stage melanoma cells. These predictions are consistent with our experimental results as well as clinical observations. Our co-culture experiments show that the senescent fibroblasts promote the growth and invasion of non-tumorigenic melanoma cells. We also observed increased proteolytic activity in stromal fields adjacent to melanoma lesions in human histology. Collectively, senescent fibroblasts create a pro-oncogenic environment that synergizes with mutations to drive melanoma initiation and progression and should therefore be considered as a potential future therapeutic target.
△ Less
Submitted 3 April, 2013;
originally announced April 2013.
-
A Markov chain model of evolution in asexually reproducing populations: insight and analytical tractability in the evolutionary process
Authors:
Daniel Nichol,
Peter Jeavons,
Robert Bonomo,
Philip K. Maini,
Jerome L. Paul,
Robert A. Gatenby,
Alexander R. A. Anderson,
Jacob G. Scott
Abstract:
The evolutionary process has been modelled in many ways using both stochastic and deterministic models. We develop an algebraic model of evolution in a population of asexually reproducing organisms in which we represent a stochastic walk in phenotype space, constrained to the edges of an underlying graph representing the genotype, with a time-homogeneous Markov Chain. We show its equivalence to a…
▽ More
The evolutionary process has been modelled in many ways using both stochastic and deterministic models. We develop an algebraic model of evolution in a population of asexually reproducing organisms in which we represent a stochastic walk in phenotype space, constrained to the edges of an underlying graph representing the genotype, with a time-homogeneous Markov Chain. We show its equivalence to a more standard, explicit stochastic model and show the algebraic model's superiority in computational efficiency. Because of this increase in efficiency, we offer the ability to simulate the evolution of much larger populations in more realistic genotype spaces. Further, we show how the algebraic properties of the Markov Chain model can give insight into the evolutionary process and allow for analysis using familiar linear algebraic methods.
△ Less
Submitted 17 January, 2013;
originally announced January 2013.
-
Intrinsic cell factors that influence tumourigenicity in cancer stem cells - towards hallmarks of cancer stem cells
Authors:
Jacob G. Scott,
Prakash Chinnaiyan,
Alexander R. A. Anderson,
Anita Hjelmeland,
David Basanta
Abstract:
Since the discovery of a cancer initiating side population in solid tumours, studies focussing on the role of so-called cancer stem cells in cancer initiation and progression have abounded. The biological interrogation of these cells has yielded volumes of information about their behaviour, but there has, as of yet, not been many actionable generalised theoretical conclusions. To address this poin…
▽ More
Since the discovery of a cancer initiating side population in solid tumours, studies focussing on the role of so-called cancer stem cells in cancer initiation and progression have abounded. The biological interrogation of these cells has yielded volumes of information about their behaviour, but there has, as of yet, not been many actionable generalised theoretical conclusions. To address this point, we have created a hybrid, discrete/continuous computational cellular automaton model of a generalised stem-cell driven tissue and explored the phenotypic traits inherent in the inciting cell and the resultant tissue growth. We identify the regions in phenotype parameter space where these initiating cells are able to cause a disruption in homeostasis, leading to tissue overgrowth and tumour formation. As our parameters and model are non-specific, they could apply to any tissue cancer stem-cell and do not assume specific genetic mutations. In this way, our model suggests that targeting these phenotypic traits could represent generalizable strategies across cancer types and represents a first attempt to identify the hallmarks of cancer stem cells.
△ Less
Submitted 20 August, 2013; v1 submitted 16 January, 2013;
originally announced January 2013.
-
A mathematical model of tumor self-seeding reveals secondary metastatic deposits as drivers of primary tumor growth
Authors:
Jacob G Scott,
David Basanta,
Alexander R. A. Anderson,
Philip Gerlee
Abstract:
Two models of circulating tumor cell (CTC) dynamics have been proposed to explain the phenomenon of tumor 'self-seeding', whereby CTCs repopulate the primary tumor and accelerate growth: Primary Seeding, where cells from a primary tumor shed into the vasculature and return back to the primary themselves; and Secondary Seeding, where cells from the primary first metastasize in a secondary tissue an…
▽ More
Two models of circulating tumor cell (CTC) dynamics have been proposed to explain the phenomenon of tumor 'self-seeding', whereby CTCs repopulate the primary tumor and accelerate growth: Primary Seeding, where cells from a primary tumor shed into the vasculature and return back to the primary themselves; and Secondary Seeding, where cells from the primary first metastasize in a secondary tissue and form microscopic secondary deposits, which then shed cells into the vasculature returning to the primary. These two models are difficult to distinguish experimentally, yet the differences between them is of great importance to both our understanding of the metastatic process and also for designing methods of intervention. Therefore we developed a mathematical model to test the relative likelihood of these two phenomena in the subset of tumours whose shed CTCs first encounter the lung capillary bed, and show that Secondary Seeding is several orders of magnitude more likely than Primary seeding. We suggest how this difference could affect tumour evolution, progression and therapy, and propose several possible methods of experimental validation.
△ Less
Submitted 25 February, 2013; v1 submitted 23 May, 2012;
originally announced May 2012.
-
A validated model of serum anti-Mullerian hormone from conception to menopause
Authors:
T W Kelsey,
P Wright,
S M Nelson,
R A Anderson,
W H B Wallace
Abstract:
Background: Anti-Mullerian hormone (AMH) is a product of growing ovarian follicles. The concentration of AMH in blood may also reflect the non-growing follicle (NGF) population, i.e. the ovarian reserve, and be of value in predicting reproductive lifespan. A full description of AMH production up to the menopause has not been reported. Methodology/Principal Findings: By searching the published lite…
▽ More
Background: Anti-Mullerian hormone (AMH) is a product of growing ovarian follicles. The concentration of AMH in blood may also reflect the non-growing follicle (NGF) population, i.e. the ovarian reserve, and be of value in predicting reproductive lifespan. A full description of AMH production up to the menopause has not been reported. Methodology/Principal Findings: By searching the published literature for AMH concentrations in healthy pre-menopausal females, and using our own data (combined n = 3,260) we have generated and robustly validated the first model of AMH concentration from conception to menopause. This model shows that 34% of the variation in AMH is due to age alone. We have shown that AMH peaks at age 24.5 years, followed by a decline to the menopause. We have also shown that there is a neonatal peak and a potential pre-pubertal peak. Our model allows us to generate normative data at all ages. Conclusions/Significance: These data highlight key inflection points in ovarian follicle dynamics. This first validated model of circulating AMH in healthy females describes a transition period in early adulthood, after which AMH reflects the progressive loss of the NGF pool. The existence of a neonatal increase in gonadal activity is confirmed for females. An improved understanding of the relationship between circulating AMH and age will lead more accurate assessment of ovarian reserve for the individual woman.
△ Less
Submitted 7 June, 2011;
originally announced June 2011.
-
A Spatial Model of Tumor-Host Interaction: Application of Chemotherapy
Authors:
Peter Hinow,
Philip Gerlee,
Lisa J. McCawley,
Vito Quaranta,
Madalina Ciobanu,
Shizhen Wang,
Jason M. Graham,
Bruce P. Ayati,
Jonathan Claridge,
Kristin R. Swanson,
Mary Loveless,
Alexander R. A. Anderson
Abstract:
In this paper we consider chemotherapy in a spatial model of tumor growth. The model, which is of reaction-diffusion type, takes into account the complex interactions between the tumor and surrounding stromal cells by including densities of endothelial cells and the extra-cellular matrix. When no treatment is applied the model reproduces the typical dynamics of early tumor growth. The initially…
▽ More
In this paper we consider chemotherapy in a spatial model of tumor growth. The model, which is of reaction-diffusion type, takes into account the complex interactions between the tumor and surrounding stromal cells by including densities of endothelial cells and the extra-cellular matrix. When no treatment is applied the model reproduces the typical dynamics of early tumor growth. The initially avascular tumor reaches a diffusion limited size of the order of millimeters and initiates angiogenesis through the release of vascular endothelial growth factor (VEGF) secreted by hypoxic cells in the core of the tumor. This stimulates endothelial cells to migrate towards the tumor and establishes a nutrient supply sufficient for sustained invasion. To this model we apply cytostatic treatment in the form of a VEGF-inhibitor, which reduces the proliferation and chemotaxis of endothelial cells. This treatment has the capability to reduce tumor mass, but more importantly, we were able to determine that inhibition of endothelial cell proliferation is the more important of the two cellular functions targeted by the drug. Further, we considered the application of a cytotoxic drug that targets proliferating tumor cells. The drug was treated as a diffusible substance entering the tissue from the blood vessels. Our results show that depending on the characteristics of the drug it can either reduce the tumor mass significantly or in fact accelerate the growth rate of the tumor. This result seems to be due to complicated interplay between the stromal and tumor cell types and highlights the importance of considering chemotherapy in a spatial context.
△ Less
Submitted 9 April, 2009; v1 submitted 6 October, 2008;
originally announced October 2008.
-
Stability Analysis of a Hybrid Cellular Automaton Model of Cell Colony Growth
Authors:
P. Gerlee,
A. R. A Anderson
Abstract:
Cell colonies of bacteria, tumour cells and fungi, under nutrient limited growth conditions, exhibit complex branched growth patterns. In order to investigate this phenomenon we present a simple hybrid cellular automaton model of cell colony growth. In the model the growth of the colony is limited by a nutrient that is consumed by the cells and which inhibits cell division if it falls below a ce…
▽ More
Cell colonies of bacteria, tumour cells and fungi, under nutrient limited growth conditions, exhibit complex branched growth patterns. In order to investigate this phenomenon we present a simple hybrid cellular automaton model of cell colony growth. In the model the growth of the colony is limited by a nutrient that is consumed by the cells and which inhibits cell division if it falls below a certain threshold. Using this model we have investigated how the nutrient consumption rate of the cells affects the growth dynamics of the colony. We found that for low consumption rates the colony takes on a Eden-like morphology, while for higher consumption rates the morphology of the colony is branched with a fractal geometry. These findings are in agreement with previous results, but the simplicity of the model presented here allows for a linear stability analysis of the system. By observing that the local growth of the colony is proportional to the flux of the nutrient we derive an approximate dispersion relation for the growth of the colony interface. This dispersion relation shows that the stability of the growth depends on how far the nutrient penetrates into the colony. For low nutrient consumption rates the penetration distance is large, which stabilises the growth, while for high consumption rates the penetration distance is small, which leads to unstable branched growth. When the penetration distance vanishes the dispersion relation is reduced to the one describing Laplacian growth without ultra-violet regularisation. The dispersion relation was verified by measuring how the average branch width depends on the consumption rate of the cells and shows good agreement between theory and simulations.
△ Less
Submitted 8 May, 2007;
originally announced May 2007.
-
Computational Methods and Results for Structured Multiscale Models of Tumor Invasion
Authors:
Bruce P. Ayati,
Glenn F. Webb,
Alexander R. A. Anderson
Abstract:
We present multiscale models of cancer tumor invasion with components at the molecular, cellular, and tissue levels. We provide biological justifications for the model components, present computational results from the model, and discuss the scientific-computing methodology used to solve the model equations. The models and methodology presented in this paper form the basis for developing and tre…
▽ More
We present multiscale models of cancer tumor invasion with components at the molecular, cellular, and tissue levels. We provide biological justifications for the model components, present computational results from the model, and discuss the scientific-computing methodology used to solve the model equations. The models and methodology presented in this paper form the basis for developing and treating increasingly complex, mechanistic models of tumor invasion that will be more predictive and less phenomenological. Because many of the features of the cancer models, such as taxis, aging and growth, are seen in other biological systems, the models and methods discussed here also provide a template for handling a broader range of biological problems.
△ Less
Submitted 14 April, 2005; v1 submitted 16 March, 2005;
originally announced March 2005.