Enhanced perfusion following exposure to radiotherapy: A theoretical investigation.

Tumour angiogenesis leads to the formation of blood vessels that are structurally and spatially heterogeneous. Poor blood perfusion, in conjunction with increased hypoxia and oxygen heterogeneity, impairs a tumour's response to radiotherapy. The optimal strategy for enhancing tumour perfusion remains unclear, preventing its regular deployment in combination therapies. In this work, we first identify vascular architectural features that correlate with enhanced perfusion following radiotherapy, using in vivo imaging data from vascular tumours. Then, we present a novel computational model to determine the relationship between these architectural features and blood perfusion in silico. If perfusion is defined to be the proportion of vessels that support blood flow, we find that vascular networks with small mean diameters and large numbers of angiogenic sprouts show the largest increases in perfusion post-irradiation for both biological and synthetic tumours. We also identify cases where perfusion increases due to the pruning of hypoperfused vessels, rather than blood being rerouted. These results indicate the importance of considering network composition when determining the optimal irradiation strategy. In the future, we aim to use our findings to identify tumours that are good candidates for perfusion enhancement and to improve the efficacy of combination therapies.

Using a probabilistic approach to derive a two-phase model of flow-induced cell migration.

Interstitial fluid flow is a feature of many solid tumors. In vitro experiments have shown that such fluid flow can direct tumor cell movement upstream or downstream depending on the balance between the competing mechanisms of tensotaxis (cell migration up stress gradients) and autologous chemotaxis (downstream cell movement in response to flow-induced gradients of self-secreted chemoattractants). In this work we develop a probabilistic-continuum, two-phase model for cell migration in response to interstitial flow. We use a kinetic description for the cell velocity probability density function, and model the flow-dependent mechanical and chemical stimuli as forcing terms that bias cell migration upstream and downstream. Using velocity-space averaging, we reformulate the model as a system of continuum equations for the spatiotemporal evolution of the cell volume fraction and flux in response to forcing terms that depend on the local direction and magnitude of the mechanochemical cues. We specialize our model to describe a one-dimensional cell layer subject to fluid flow. Using a combination of numerical simulations and asymptotic analysis, we delineate the parameter regime where transitions from downstream to upstream cell migration occur. As has been observed experimentally, the model predicts downstream-oriented chemotactic migration at low cell volume fractions, and upstream-oriented tensotactic migration at larger volume fractions. We show that the locus of the critical volume fraction, at which the system transitions from downstream to upstream migration, is dominated by the ratio of the rate of chemokine secretion and advection. Our model also predicts that, because the tensotactic stimulus depends strongly on the cell volume fraction, upstream, tensotaxis-dominated migration occurs only transiently when the cells are initially seeded, and transitions to downstream, chemotaxis-dominated migration occur at later times due to the dispersive effect of cell diffusion.

A mathematical framework for the emergence of winners and losers in cell competition.

Cell competition is a process in multicellular organisms where cells interact with their neighbours to determine a "winner" or "loser" status. The loser cells are eliminated through programmed cell death, leaving only the winner cells to populate the tissue. Cell competition is context-dependent; the same cell type can win or lose depending on the cell type it is competing against. Hence, winner/loser status is an emergent property. A key question in cell competition is: how do cells acquire their winner/loser status? In this paper, we propose a mathematical framework for studying the emergence of winner/loser status based on a set of quantitative criteria that distinguishes competitive from non-competitive outcomes. We apply this framework in a cell-based modelling context, to both highlight the crucial role of active cell death in cell competition and identify the factors that drive cell competition.

Tumour induced osteopenia due to phosphaturic mesenchymal sinonasal tumour presenting with delayed onset insufficiency fractures.

We present a case of a 48-year-old female who presented with epistaxis. Magnetic resonance imaging (MRI) revealed a mass within the left nasal cavity which was revealed to be a phosphaturic mesenchymal sinonasal tumour. The patient defaulted treatment at this stage and later re-presented with pelvic and groin pain for which plain radiographs and computed tomography (CT) scan demonstrated diffuse osteopenia and multiple pelvic fractures of varying ages. MRI of the pelvis and both thighs revealed abnormal marrow signal of the bones and confirmed the presence of pelvic fractures. Multiple pseudo-fractures were seen at both femurs and scapula. The radiological findings along with abnormal biochemical markers were attributed to the paraneoplastic entity of tumour induced osteomalacia, in the context of unresected phosphaturic mesenchymal tumour. The tumour was resected, and patient showed complete reversal of the associated biochemical abnormalities. This case exemplifies that with early identification and complete resection of the causative tumour, the prognosis is excellent.

Veterans Health Administration response to 2021 recall of Philips Respironics devices: A case study.

This case study describes, for the time frame of June 2021 through August 2022, the U.S. Veterans Health Administration (VHA) organizational response to a manufacturer's recall of positive airway pressure devices used in the treatment of sleep disordered breathing. VHA estimated it could take over a year for Veterans to receive replacement devices. Veterans awaiting a replacement faced a dilemma. They could continue using the recalled devices and bear the product safety risks that led to the recall, or they could stop using them and bear the risks of untreated sleep disordered breathing. Using a program monitoring approach, we report on the processes VHA put in place to respond to the recall. Specifically, we report on the strategic, service, and operational plans associated with VHA's response to the recall for Veterans needing replacement devices. In program monitoring, the strategic plan reflects the internal process objectives for the program. The service plan articulates how the delivery of services will intersect the customer journey. The operational plan describes how the program's resources and actions must support the service delivery plan. VHA's strategic plan featured a clinician-led, as opposed to primarily legal or administrative response to the recall. The recall response team also engaged with VHA's medical ethics service to articulate an ethical framework guiding the allocation of replacement devices under conditions of scarcity. This framework proposed allocating scarce devices to Veterans according to their clinical need. The service plan invited Veterans to schedule visits with sleep providers who could assess their clinical need and counsel them accordingly. The operational plan distributed devices according to clinical need as they became available. Monitoring our program processes in real time helped VHA launch and adapt its response to a recall affecting more than 700,000 Veterans.

An uncommon case of traumatic pisiform dislocation with triquetral fracture.

The pisiform is a sesamoid bone that acts as one of the key medial stabilizers of the wrist. We present a case of a 35-year-old gentleman who presented with medial wrist pain following a fall while rollerblading. Radiographs and Magnetic resonance imaging (MRI) revealed a rare combination of an acute pisiform dislocation with associated triquetral fracture. Subsequently, he was successfully treated with excision of the pisiform. Pisiform dislocation is an uncommon injury and can easily be missed in an acute emergency presentation. Therefore, it is important to be aware of the characteristic imaging appearance to avoid a delay in diagnosis and treatment.

A DNA-structured mathematical model of cell-cycle progression in cyclic hypoxia.

New experimental data have shown how the periodic exposure of cells to low oxygen levels (i.e., cyclic hypoxia) impacts their progress through the cell-cycle. Cyclic hypoxia has been detected in tumours and linked to poor prognosis and treatment failure. While fluctuating oxygen environments can be reproduced in vitro, the range of oxygen cycles that can be tested is limited. By contrast, mathematical models can be used to predict the response to a wide range of cyclic dynamics. Accordingly, in this paper we develop a mechanistic model of the cell-cycle that can be combined with in vitro experiments to better understand the link between cyclic hypoxia and cell-cycle dysregulation. A distinguishing feature of our model is the inclusion of impaired DNA synthesis and cell-cycle arrest due to periodic exposure to severely low oxygen levels. Our model decomposes the cell population into five compartments and a time-dependent delay accounts for the variability in the duration of the S phase which increases in severe hypoxia due to reduced rates of DNA synthesis. We calibrate our model against experimental data and show that it recapitulates the observed cell-cycle dynamics. We use the calibrated model to investigate the response of cells to oxygen cycles not yet tested experimentally. When the re-oxygenation phase is sufficiently long, our model predicts that cyclic hypoxia simply slows cell proliferation since cells spend more time in the S phase. On the contrary, cycles with short periods of re-oxygenation are predicted to lead to inhibition of proliferation, with cells arresting from the cell-cycle in the G2 phase. While model predictions on short time scales (about a day) are fairly accurate (i.e, confidence intervals are small), the predictions become more uncertain over longer periods. Hence, we use our model to inform experimental design that can lead to improved model parameter estimates and validate model predictions.

Dependence of cell-free-layer width on rheological parameters: Combining empirical data on flow separation at microvascular bifurcations with geometrical considerations.

When blood flows through vessel networks, red blood cells (RBCs) are typically concentrated close to the vessel center line, creating a lubrication layer near the vessel wall. As RBCs bind oxygen, the width of this cell-free layer (CFL) impacts not only the blood rheology inside the vasculature, but also oxygen delivery to the tissues they perfuse and, hence, their function. Existing attempts to relate the width of the CFL to the rheological properties of the blood and the geometrical properties of the vessel are based on an analysis of the forces acting on RBCs suspended in the blood. However, the complexity of interactions in the blood makes this a challenging task. Here, we propose an alternative, two-step approach to derive such a functional relationship. First, we extend widely accepted empirical fits describing the minimum flow fraction needed for RBCs to enter a daughter vessel downstream of a microvascular bifurcation so that it depends not only on the diameter and discharge haematocrit of the parent vessel, but also on its average shear rate. Second, we propose a simple geometrical model for the minimum flow fraction based on the cross-sectional blood flow profile within the parent vessel upstream of the bifurcation-considering uniform, parabolic, and blunt velocity profiles-and derive the leading-order approximation to this model for small CFL widths. By equating the functional relationships obtained using these two approaches, we derive expressions relating the CFL width to the vessel diameter, discharge haematocrit, and mean shear rate. The resulting expressions are in good agreement with available in vivo data and represent a promising basis for future research.

Abnormal morphology biases hematocrit distribution in tumor vasculature and contributes to heterogeneity in tissue oxygenation.

Oxygen heterogeneity in solid tumors is recognized as a limiting factor for therapeutic efficacy. This heterogeneity arises from the abnormal vascular structure of the tumor, but the precise mechanisms linking abnormal structure and compromised oxygen transport are only partially understood. In this paper, we investigate the role that red blood cell (RBC) transport plays in establishing oxygen heterogeneity in tumor tissue. We focus on heterogeneity driven by network effects, which are challenging to observe experimentally due to the reduced fields of view typically considered. Motivated by our findings of abnormal vascular patterns linked to deviations from current RBC transport theory, we calculated average vessel lengths [Formula: see text] and diameters [Formula: see text] from tumor allografts of three cancer cell lines and observed a substantial reduction in the ratio [Formula: see text] compared to physiological conditions. Mathematical modeling reveals that small values of the ratio λ (i.e., [Formula: see text]) can bias hematocrit distribution in tumor vascular networks and drive heterogeneous oxygenation of tumor tissue. Finally, we show an increase in the value of λ in tumor vascular networks following treatment with the antiangiogenic cancer agent DC101. Based on our findings, we propose λ as an effective way of monitoring the efficacy of antiangiogenic agents and as a proxy measure of perfusion and oxygenation in tumor tissue undergoing antiangiogenic treatment.

Case Report of a Mediastinal Vascular Malformation Mimicking Esophageal Varices on Endoscopy, with Emphasis on Radiological Findings.

BACKGROUND Mediastinal vascular malformations are rare, and most patients are asymptomatic or present with unrelated symptoms. Imaging can be challenging to interpret, but plays an important role in diagnosis and prognostication. CASE REPORT We present the case of a 48-year-old man with history of intravenous drug abuse and incompletely treated pulmonary tuberculosis. A computed tomography (CT) scan done for respiratory symptoms showed an extensive soft-tissue mass in the mediastinum and upper abdomen, initially thought to represent tuberculous adenitis with possible esophageal involvement, which appeared variceal in nature on endoscopy. Further investigation with open mediastinal biopsy and magnetic resonance imaging (MRI) eventually led to the diagnosis of a low-flow venous mediastinal vascular malformation. The patient responded well to conservative management, with the malformation remaining stable on follow-up CT up to a decade later. CONCLUSIONS Radiologists should be aware of the rare but important differential diagnosis of a vascular malformation, particularly when an extensive infiltrative calcified mediastinal soft-tissue mass is encountered. Multi-modality imaging, particularly MRI, which can demonstrate typical features, is crucial for diagnosis and prognostication, thereby avoiding unnecessary invasive procedures and treatment.

Chaste: Cancer, Heart and Soft Tissue Environment.

Chaste (Cancer, Heart And Soft Tissue Environment) is an open source simulation package for the numerical solution of mathematical models arising in physiology and biology. To date, Chaste development has been driven primarily by applications that include continuum modelling of cardiac electrophysiology ('Cardiac Chaste'), discrete cell-based modelling of soft tissues ('Cell-based Chaste'), and modelling of ventilation in lungs ('Lung Chaste'). Cardiac Chaste addresses the need for a high-performance, generic, and verified simulation framework for cardiac electrophysiology that is freely available to the scientific community. Cardiac chaste provides a software package capable of realistic heart simulations that is efficient, rigorously tested, and runs on HPC platforms. Cell-based Chaste addresses the need for efficient and verified implementations of cell-based modelling frameworks, providing a set of extensible tools for simulating biological tissues. Computational modelling, along with live imaging techniques, plays an important role in understanding the processes of tissue growth and repair. A wide range of cell-based modelling frameworks have been developed that have each been successfully applied in a range of biological applications. Cell-based Chaste includes implementations of the cellular automaton model, the cellular Potts model, cell-centre models with cell representations as overlapping spheres or Voronoi tessellations, and the vertex model. Lung Chaste addresses the need for a novel, generic and efficient lung modelling software package that is both tested and verified. It aims to couple biophysically-detailed models of airway mechanics with organ-scale ventilation models in a package that is freely available to the scientific community. Chaste is designed to be modular and extensible, providing libraries for common scientific computing infrastructure such as linear algebra operations, finite element meshes, and ordinary and partial differential equation solvers. This infrastructure is used by libraries for specific applications, such as continuum mechanics, cardiac models, and cell-based models. The software engineering techniques used to develop Chaste are intended to ensure code quality, re-usability and reliability. Primary applications of the software include cardiac and respiratory physiology, cancer and developmental biology.

Cooperation, competition and antibiotic resistance in bacterial colonies.

Bacteria commonly live in dense and genetically diverse communities associated with surfaces. In these communities, competition for resources and space is intense, and yet we understand little of how this affects the spread of antibiotic-resistant strains. Here, we study interactions between antibiotic-resistant and susceptible strains using in vitro competition experiments in the opportunistic pathogen Pseudomonas aeruginosa and in silico simulations. Selection for intracellular resistance to streptomycin is very strong in colonies, such that resistance is favoured at very low antibiotic doses. In contrast, selection for extracellular resistance to carbenicillin is weak in colonies, and high doses of antibiotic are required to select for resistance. Manipulating the density and spatial structure of colonies reveals that this difference is partly explained by the fact that the local degradation of carbenicillin by ß-lactamase-secreting cells protects neighbouring sensitive cells from carbenicillin. In addition, we discover a second unexpected effect: the inducible elongation of cells in response to carbenicillin allows sensitive cells to better compete for the rapidly growing colony edge. These combined effects mean that antibiotic treatment can select against antibiotic-resistant strains, raising the possibility of treatment regimes that suppress sensitive strains while limiting the rise of antibiotic resistance. We argue that the detailed study of bacterial interactions will be fundamental to understanding and overcoming antibiotic resistance.

Reporting of CT cervical spine after office hours by radiology trainees-analysis of discrepancy rates and RADPEER scores.

PURPOSE: To review the errors made by radiology trainees in the reporting of cervical spine CTs (CCT) and to compare the discrepancy rates between the stages of training. METHODS: All CCTs reported by trainees after office hours between January 2015 and December 2015 were retrospectively reviewed by a team of five musculoskeletal consultants with experience ranging between 7 and 15 years. Discrepancies between the provisional report by the trainee and the findings by the musculoskeletal consultants were graded according to the RADPEER scoring system. Sensitivity and specificity of the trainees were assessed. RESULTS: Of 254 CCT provisional reports, there were 12 (4.7%) discrepancies, of which 5 (2.0%) discrepancies were likely to be clinically significant. We found a clinically significant difference between the stage of training of the trainee and RADPEER score (P = 0.023). The sensitivity and specificity of the senior radiology trainees were 97.0 and 98.1%, respectively, and that of the junior radiology trainees were 80 and 98.0% respectively (P = 0.039). Conditions misinterpreted as fractures include degenerative changes (n = 2) and nutrient vessel (n = 1). Other missed abnormalities include ossification of the posterior longitudinal ligament (n = 1), fracture of the foramen transversarium (n = 2), vertebral body fractures (n = 2), articular facet fractures (n = 2), and transverse process fractures (n = 2). CONCLUSION: Cervical spine CTs performed after office hours can be safely interpreted by senior radiology trainees to a reasonable degree, although a targeted intervention to improve diagnostic performance of junior radiology trainees may be of clinical benefit.

The importance of geometry in the corneal micropocket angiogenesis assay.

The corneal micropocket angiogenesis assay is an experimental protocol for studying vessel network formation, or neovascularization, in vivo. The assay is attractive due to the ease with which the developing vessel network can be observed in the same animal over time. Measurements from the assay have been used in combination with mathematical modeling to gain insights into the mechanisms of angiogenesis. While previous modeling studies have adopted planar domains to represent the assay, the hemispherical shape of the cornea and asymmetric positioning of the angiogenic source can be seen to affect vascular patterning in experimental images. As such, we aim to better understand: i) how the geometry of the assay influences vessel network formation and ii) how to relate observations from planar domains to those in the hemispherical cornea. To do so, we develop a three-dimensional, off-lattice mathematical model of neovascularization in the cornea, using a spatially resolved representation of the assay for the first time. Relative to the detailed model, we predict that the adoption of planar geometries has a noticeable impact on vascular patterning, leading to increased vessel 'merging', or anastomosis, in particular when circular geometries are adopted. Significant differences in the dynamics of diffusible aniogenesis simulators are also predicted between different domains. In terms of comparing predictions across domains, the 'distance of the vascular front to the limbus' metric is found to have low sensitivity to domain choice, while metrics such as densities of tip cells and vessels and 'vascularized fraction' are sensitive to domain choice. Given the widespread adoption and attractive simplicity of planar tissue domains, both in silico and in vitro, the differences identified in the present study should prove useful in relating the results of previous and future theoretical studies of neovascularization to in vivo observations in the cornea.

Microvessel Chaste: An Open Library for Spatial Modeling of Vascularized Tissues.

Spatial models of vascularized tissues are widely used in computational physiology. We introduce a software library for composing multiscale, multiphysics models for applications including tumor growth, angiogenesis, osteogenesis, coronary perfusion, and oxygen delivery. Composition of such models is time consuming, with many researchers writing custom software. Recent advances in imaging have produced detailed three-dimensional (3D) datasets of vascularized tissues at the scale of individual cells. To fully exploit such data there is an increasing need for software that allows user-friendly composition of efficient, 3D models of vascularized tissues, and comparison of predictions with in vivo or in vitro experiments and alternative computational formulations. Microvessel Chaste can be used to build simulations of vessel growth and adaptation in response to mechanical and chemical stimuli; intra- and extravascular transport of nutrients, growth factors and drugs; and cell proliferation in complex 3D geometries. In addition, it can be used to develop custom software for integrating modeling with experimental data processing workflows, facilitated by a comprehensive Python interface to solvers implemented in C++. This article links to two reproducible example problems, showing how the library can be used to build simulations of tumor growth and angiogenesis with realistic vessel networks.

Bayesian Calibration, Validation and Uncertainty Quantification for Predictive Modelling of Tumour Growth: A Tutorial.

In this work, we present a pedagogical tumour growth example, in which we apply calibration and validation techniques to an uncertain, Gompertzian model of tumour spheroid growth. The key contribution of this article is the discussion and application of these methods (that are not commonly employed in the field of cancer modelling) in the context of a simple model, whose deterministic analogue is widely known within the community. In the course of the example, we calibrate the model against experimental data that are subject to measurement errors, and then validate the resulting uncertain model predictions. We then analyse the sensitivity of the model predictions to the underlying measurement model. Finally, we propose an elementary learning approach for tuning a threshold parameter in the validation procedure in order to maximize predictive accuracy of our validated model.

Comparing individual-based approaches to modelling the self-organization of multicellular tissues.

The coordinated behaviour of populations of cells plays a central role in tissue growth and renewal. Cells react to their microenvironment by modulating processes such as movement, growth and proliferation, and signalling. Alongside experimental studies, computational models offer a useful means by which to investigate these processes. To this end a variety of cell-based modelling approaches have been developed, ranging from lattice-based cellular automata to lattice-free models that treat cells as point-like particles or extended shapes. However, it remains unclear how these approaches compare when applied to the same biological problem, and what differences in behaviour are due to different model assumptions and abstractions. Here, we exploit the availability of an implementation of five popular cell-based modelling approaches within a consistent computational framework, Chaste (http://www.cs.ox.ac.uk/chaste). This framework allows one to easily change constitutive assumptions within these models. In each case we provide full details of all technical aspects of our model implementations. We compare model implementations using four case studies, chosen to reflect the key cellular processes of proliferation, adhesion, and short- and long-range signalling. These case studies demonstrate the applicability of each model and provide a guide for model usage.

Cell morphology drives spatial patterning in microbial communities.

The clearest phenotypic characteristic of microbial cells is their shape, but we do not understand how cell shape affects the dense communities, known as biofilms, where many microbes live. Here, we use individual-based modeling to systematically vary cell shape and study its impact in simulated communities. We compete cells with different cell morphologies under a range of conditions and ask how shape affects the patterning and evolutionary fitness of cells within a community. Our models predict that cell shape will strongly influence the fate of a cell lineage: we describe a mechanism through which coccal (round) cells rise to the upper surface of a community, leading to a strong spatial structuring that can be critical for fitness. We test our predictions experimentally using strains of Escherichia coli that grow at a similar rate but differ in cell shape due to single amino acid changes in the actin homolog MreB. As predicted by our model, cell types strongly sort by shape, with round cells at the top of the colony and rod cells dominating the basal surface and edges. Our work suggests that cell morphology has a strong impact within microbial communities and may offer new ways to engineer the structure of synthetic communities.

Towards More Comprehensive Projections of Urban Heat-Related Mortality: Estimates for New York City under Multiple Population, Adaptation, and Climate Scenarios.

BACKGROUND: High temperatures have substantial impacts on mortality and, with growing concerns about climate change, numerous studies have developed projections of future heat-related deaths around the world. Projections of temperature-related mortality are often limited by insufficient information to formulate hypotheses about population sensitivity to high temperatures and future demographics. OBJECTIVES: The present study derived projections of temperature-related mortality in New York City by taking into account future patterns of adaptation or demographic change, both of which can have profound influences on future health burdens. METHODS: We adopted a novel approach to modeling heat adaptation by incorporating an analysis of the observed population response to heat in New York City over the course of eight decades. This approach projected heat-related mortality until the end of the 21st century based on observed trends in adaptation over a substantial portion of the 20th century. In addition, we incorporated a range of new scenarios for population change until the end of the 21st century. We then estimated future heat-related deaths in New York City by combining the changing temperature-mortality relationship and population scenarios with downscaled temperature projections from the 33 global climate models (GCMs) and two Representative Concentration Pathways (RCPs). RESULTS: The median number of projected annual heat-related deaths across the 33 GCMs varied greatly by RCP and adaptation and population change scenario, ranging from 167 to 3,331 in the 2080s compared with 638 heat-related deaths annually between 2000 and 2006. CONCLUSIONS: These findings provide a more complete picture of the range of potential future heat-related mortality risks across the 21st century in New York City, and they highlight the importance of both demographic change and adaptation responses in modifying future risks. Citation: Petkova EP, Vink JK, Horton RM, Gasparrini A, Bader DA, Francis JD, Kinney PL. 2017. Towards more comprehensive projections of urban heat-related mortality: estimates for New York City under multiple population, adaptation, and climate scenarios. Environ Health Perspect 125:47-55; http://dx.doi.org/10.1289/EHP166.