Chaste: an open source C++ library for computational physiology and biology

Chaste - Cancer, Heart And Soft Tissue Environment - is an open source C++ library for the computational simulation of mathematical models developed for physiology and biology. Code development has been driven by two initial applications: cardiac electrophysiology and cancer development. A large number of cardiac electrophysiology studies have been enabled and performed, including high performance computational investigations of defibrillation on realistic human cardiac geometries. New models for the initiation and growth of tumours have been developed. In particular, cell-based simulations have provided novel insight into the role of stem cells in the colorectal crypt. Chaste is constantly evolving and is now being applied to a far wider range of problems. The code provides modules for handling common scientific computing components, such as meshes and solvers for ordinary and partial differential equations (ODEs/PDEs). Re-use of these components avoids the need for researchers to "re-invent the wheel" with each new project, accelerating the rate of progress in new applications. Chaste is developed using industrially-derived techniques, in particular test-driven development, to ensure code quality, re-use and reliability. In this article we provide examples that illustrate the types of problems Chaste can be used to solve, which can be run on a desktop computer. We highlight some scientific studies that have used or are using Chaste, and the insights they have provided. The source code, both for specific releases and the development version, is available to download under an open source Berkeley Software Distribution (BSD) licence at http://www.cs.ox.ac.uk/chaste, together with details of a mailing list and links to documentation and tutorials

Thermal behaviour of Anopheles stephensi in response to infection with malaria and fungal entomopathogens

<p>Abstract</p> <p>Background</p> <p>Temperature is a critical determinant of the development of malaria parasites in mosquitoes, and hence the geographic distribution of malaria risk, but little is known about the thermal preferences of <it>Anopheles</it>. A number of other insects modify their thermal behaviour in response to infection. These alterations can be beneficial for the insect or for the infectious agent. Given current interest in developing fungal biopesticides for control of mosquitoes, <it>Anopheles stephensi </it>were examined to test whether mosquitoes showed thermally-mediated behaviour in response to infection with fungal entomopathogens and the rodent malaria, <it>Plasmodium yoelii</it>.</p> <p>Methods</p> <p>Over two experiments, groups of <it>An. stephensi </it>were infected with one of three entomopathogenic fungi, and/or <it>P. yoelii</it>. Infected and uninfected mosquitoes were released on to a thermal gradient (14 – 38°C) for "snapshot" assessments of thermal preference during the first five days post-infection. Mosquito survival was monitored for eight days and, where appropriate, oocyst prevalence and intensity was assessed.</p> <p>Results and conclusion</p> <p>Both infected and uninfected <it>An. stephensi </it>showed a non-random distribution on the gradient, indicating some capacity to behaviourally thermoregulate. However, chosen resting temperatures were not altered by any of the infections. There is thus no evidence that thermally-mediated behaviours play a role in determining malaria prevalence or that they will influence the performance of fungal biopesticides against adult <it>Anopheles</it>.</p

Projected WIMP sensitivity of the LUX-ZEPLIN dark matter experiment

LUX-ZEPLIN (LZ) is a next-generation dark matter direct detection experiment that will operate 4850 feet underground at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector with an active mass of 7 tonnes, LZ will search primarily for low-energy interactions with weakly interacting massive particles (WIMPs), which are hypothesized to make up the dark matter in our galactic halo. In this paper, the projected WIMP sensitivity of LZ is presented based on the latest background estimates and simulations of the detector. For a 1000 live day run using a 5.6-tonne fiducial mass, LZ is projected to exclude at 90% confidence level spin-independent WIMP-nucleon cross sections above 1.4×10-48 cm2 for a 40 GeV/c2 mass WIMP. Additionally, a 5σ discovery potential is projected, reaching cross sections below the exclusion limits of recent experiments. For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of 2.3×10-43 cm2 (7.1×10-42 cm2) for a 40 GeV/c2 mass WIMP is expected. With underground installation well underway, LZ is on track for commissioning at SURF in 2020

Simulations of events for the LUX-ZEPLIN (LZ) dark matter experiment

The LUX-ZEPLIN dark matter search aims to achieve a sensitivity to the WIMP-nucleon spin-independent cross-section down to (1–2)×10−12 pb at a WIMP mass of 40 GeV/c2. This paper describes the simulations framework that, along with radioactivity measurements, was used to support this projection, and also to provide mock data for validating reconstruction and analysis software. Of particular note are the event generators, which allow us to model the background radiation, and the detector response physics used in the production of raw signals, which can be converted into digitized waveforms similar to data from the operational detector. Inclusion of the detector response allows us to process simulated data using the same analysis routines as developed to process the experimental data

Projected sensitivity of the LUX-ZEPLIN experiment to the 0νββ decay of Xe 136

The LUX-ZEPLIN (LZ) experiment will enable a neutrinoless double β decay search in parallel to the main science goal of discovering dark matter particle interactions. We report the expected LZ sensitivity to Xe136 neutrinoless double β decay, taking advantage of the significant (&gt;600 kg) Xe136 mass contained within the active volume of LZ without isotopic enrichment. After 1000 live-days, the median exclusion sensitivity to the half-life of Xe136 is projected to be 1.06×1026 years (90% confidence level), similar to existing constraints. We also report the expected sensitivity of a possible subsequent dedicated exposure using 90% enrichment with Xe136 at 1.06×1027 years

Lethal and Pre-Lethal Effects of a Fungal Biopesticide Contribute to Substantial and Rapid Control of Malaria Vectors

Rapidly emerging insecticide resistance is creating an urgent need for new active ingredients to control the adult mosquitoes that vector malaria. Biopesticides based on the spores of entomopathogenic fungi have shown considerable promise by causing very substantial mortality within 7–14 days of exposure. This mortality will generate excellent malaria control if there is a high likelihood that mosquitoes contact fungi early in their adult lives. However, where contact rates are lower, as might result from poor pesticide coverage, some mosquitoes will contact fungi one or more feeding cycles after they acquire malaria, and so risk transmitting malaria before the fungus kills them. Critics have argued that ‘slow acting’ fungal biopesticides are, therefore, incapable of delivering malaria control in real-world contexts. Here, utilizing standard WHO laboratory protocols, we demonstrate effective action of a biopesticide much faster than previously reported. Specifically, we show that transient exposure to clay tiles sprayed with a candidate biopesticide comprising spores of a natural isolate of Beauveria bassiana, could reduce malaria transmission potential to zero within a feeding cycle. The effect resulted from a combination of high mortality and rapid fungal-induced reduction in feeding and flight capacity. Additionally, multiple insecticide-resistant lines from three key African malaria vector species were completely susceptible to fungus. Thus, fungal biopesticides can block transmission on a par with chemical insecticides, and can achieve this where chemical insecticides have little impact. These results support broadening the current vector control paradigm beyond fast-acting chemical toxins

Measurement of the gamma ray background in the Davis cavern at the Sanford Underground Research Facility

Deep underground environments are ideal for low background searches due to the attenuation of cosmic rays by passage through the earth. However, they are affected by backgrounds from γ-rays emitted by 40K and the 238U and 232Th decay chains in the surrounding rock. The LUX-ZEPLIN (LZ) experiment will search for dark matter particle interactions with a liquid xenon TPC located within the Davis campus at the Sanford Underground Research Facility, Lead, South Dakota, at the 4850-foot level. In order to characterise the cavern background, in-situ γ-ray measurements were taken with a sodium iodide detector in various locations and with lead shielding. The integral count rates (0–3300 keV) varied from 596 Hz to 1355 Hz for unshielded measurements, corresponding to a total flux from the cavern walls of 1.9 ± 0.4 γ cm−2s−1. The resulting activity in the walls of the cavern can be characterised as 220 ± 60 Bq/kg of 40K, 29 ± 15 Bq/kg of 238U, and 13 ± 3 Bq/kg of 232Th

Brain Cortical Mapping by Simultaneous Recording of Functional Near Infrared Spectroscopy and Electroencephalograms from the Whole Brain During Right Median Nerve Stimulation

To investigate relationships between hemodynamic responses and neural activities in the somatosensory cortices, hemodynamic responses by near infrared spectroscopy (NIRS) and electroencephalograms (EEGs) were recorded simultaneously while subjects received electrical stimulation in the right median nerve. The statistical significance of the hemodynamic responses was evaluated by a general linear model (GLM) with the boxcar design matrix convoluted with Gaussian function. The resulting NIRS and EEGs data were stereotaxically superimposed on the reconstructed brain of each subject. The NIRS data indicated that changes in oxy-hemoglobin concentration increased at the contralateral primary somatosensory (SI) area; responses then spread to the more posterior and ipsilateral somatosensory areas. The EEG data indicated that positive somatosensory evoked potentials peaking at 22 ms latency (P22) were recorded from the contralateral SI area. Comparison of these two sets of data indicated that the distance between the dipoles of P22 and NIRS channels with maximum hemodynamic responses was less than 10 mm, and that the two topographical maps of hemodynamic responses and current source density of P22 were significantly correlated. Furthermore, when onset of the boxcar function was delayed 5–15 s (onset delay), hemodynamic responses in the bilateral parietal association cortices posterior to the SI were more strongly correlated to electrical stimulation. This suggests that GLM analysis with onset delay could reveal the temporal ordering of neural activation in the hierarchical somatosensory pathway, consistent with the neurophysiological data. The present results suggest that simultaneous NIRS and EEG recording is useful for correlating hemodynamic responses to neural activity

Simulations of events for the LUX-ZEPLIN (LZ) dark matter experiment

The LUX-ZEPLIN dark matter search aims to achieve a sensitivity to the WIMP-nucleon spin-independent cross-section down to (1–2)x10-12 pb at a WIMP mass of 40  GeV/c2. This paper describes the simulations framework that, along with radioactivity measurements, was used to support this projection, and also to provide mock data for validating reconstruction and analysis software. Of particular note are the event generators, which allow us to model the background radiation, and the detector response physics used in the production of raw signals, which can be converted into digitized waveforms similar to data from the operational detector. Inclusion of the detector response allows us to process simulated data using the same analysis routines as developed to process the experimental data