Relating multi-sequence longitudinal intensity profiles and clinical covariates in new multiple sclerosis lesions

Structural magnetic resonance imaging (MRI) can be used to detect lesions in the brains of multiple sclerosis (MS) patients. The formation of these lesions is a complex process involving inflammation, tissue damage, and tissue repair, all of which are visible on MRI. Here we characterize the lesion formation process on longitudinal, multi-sequence structural MRI from 34 MS patients and relate the longitudinal changes we observe within lesions to therapeutic interventions. In this article, we first outline a pipeline to extract voxel level, multi-sequence longitudinal profiles from four MRI sequences within lesion tissue. We then propose two models to relate clinical covariates to the longitudinal profiles. The first model is a principal component analysis (PCA) regression model, which collapses the information from all four profiles into a scalar value. We find that the score on the first PC identifies areas of slow, long-term intensity changes within the lesion at a voxel level, as validated by two experienced clinicians, a neuroradiologist and a neurologist. On a quality scale of 1 to 4 (4 being the highest) the neuroradiologist gave the score on the first PC a median rating of 4 (95% CI: [4,4]), and the neurologist gave it a median rating of 3 (95% CI: [3,3]). In the PCA regression model, we find that treatment with disease modifying therapies (p-value < 0.01), steroids (p-value < 0.01), and being closer to the boundary of abnormal signal intensity (p-value < 0.01) are associated with a return of a voxel to intensity values closer to that of normal-appearing tissue. The second model is a function-on-scalar regression, which allows for assessment of the individual time points at which the covariates are associated with the profiles. In the function-on-scalar regression both age and distance to the boundary were found to have a statistically significant association with the profiles

Predicting the effects of climate change on water yield and forest production in the northeastern United States

Rapid and simultaneous changes in temperature, precipitation and the atmospheric concentration of CO2 are predicted to occur over the next century. Simple, well-validated models of ecosystem function are required to predict the effects of these changes. This paper describes an improved version of a forest carbon and water balance model (PnET-II) and the application of the model to predict stand- and regional-level effects of changes in temperature, precipitation and atmospheric CO2 concentration. PnET-II is a simple, generalized, monthly time-step model of water and carbon balances (gross and net) driven by nitrogen availability as expressed through foliar N concentration. Improvements from the original model include a complete carbon balance and improvements in the prediction of canopy phenology, as well as in the computation of canopy structure and photosynthesis. The model was parameterized and run for 4 forest/site combinations and validated against available data for water yield, gross and net carbon exchange and biomass production. The validation exercise suggests that the determination of actual water availability to stands and the occurrence or non-occurrence of soil-based water stress are critical to accurate modeling of forest net primary production (NPP) and net ecosystem production (NEP). The model was then run for the entire NewEngland/New York (USA) region using a 1 km resolution geographic information system. Predicted long-term NEP ranged from -85 to +275 g C m-2 yr-1 for the 4 forest/site combinations, and from -150 to 350 g C m-2 yr-1 for the region, with a regional average of 76 g C m-2 yr-1. A combination of increased temperature (+6*C), decreased precipitation (-15%) and increased water use efficiency (2x, due to doubling of CO2) resulted generally in increases in NPP and decreases in water yield over the region

Intraoperative Management of Robotic-Assisted Versus Open Radical Prostatectomy

Robotic-assisted laparoscopic radical prostatectomy was found to be a shorter procedure characterized by minimal blood loss, reduced fluid requirements, and shorter hospital stay compared with traditional open procedures

Cell and molecular transitions during efficient dedifferentiation

Dedifferentiation is a critical response to tissue damage, yet is not well understood, even at a basic phenomenological level. Developing Dictyostelium cells undergo highly efficient dedifferentiation, completed by most cells within 24 hr. We use this rapid response to investigate the control features of dedifferentiation, combining single cell imaging with high temporal resolution transcriptomics. Gene expression during dedifferentiation was predominantly a simple reversal of developmental changes, with expression changes not following this pattern primarily associated with ribosome biogenesis. Mutation of genes induced early in dedifferentiation did not strongly perturb the reversal of development. This apparent robustness may arise from adaptability of cells: the relative temporal ordering of cell and molecular events was not absolute, suggesting cell programmes reach the same end using different mechanisms. In addition, although cells start from different fates, they rapidly converged on a single expression trajectory. These regulatory features may contribute to dedifferentiation responses during regeneration

Phase protection of Fano-Feshbach resonances

Decay of bound states due to coupling with free particle states is a general phenomenon occurring at energy scales from MeV in nuclear physics to peV in ultracold atomic gases. Such a coupling gives rise to Fano-Feshbach resonances (FFR) that have become key to understanding and controlling interactions—in ultracold atomic gases, but also between quasiparticles, such as microcavity polaritons. Their energy positions were shown to follow quantum chaotic statistics. In contrast, their lifetimes have so far escaped a similarly comprehensive understanding. Here, we show that bound states, despite being resonantly coupled to a scattering state, become protected from decay whenever the relative phase is a multiple of π. We observe this phenomenon by measuring lifetimes spanning four orders of magnitude for FFR of spin–orbit excited molecular ions with merged beam and electrostatic trap experiments. Our results provide a blueprint for identifying naturally long-lived states in a decaying quantum system

Narrative-based computational modelling of the Gp130/JAK/STAT signalling pathway.

BACKGROUND: Appropriately formulated quantitative computational models can support researchers in understanding the dynamic behaviour of biological pathways and support hypothesis formulation and selection by "in silico" experimentation. An obstacle to widespread adoption of this approach is the requirement to formulate a biological pathway as machine executable computer code. We have recently proposed a novel, biologically intuitive, narrative-style modelling language for biologists to formulate the pathway which is then automatically translated into an executable format and is, thus, usable for analysis via existing simulation techniques. RESULTS: Here we use a high-level narrative language in designing a computational model of the gp130/JAK/STAT signalling pathway and show that the model reproduces the dynamic behaviour of the pathway derived by biological observation. We then "experiment" on the model by simulation and sensitivity analysis to define those parameters which dominate the dynamic behaviour of the pathway. The model predicts that nuclear compartmentalisation and phosphorylation status of STAT are key determinants of the pathway and that alternative mechanisms of signal attenuation exert their influence on different timescales. CONCLUSION: The described narrative model of the gp130/JAK/STAT pathway represents an interesting case study showing how, by using this approach, researchers can model biological systems without explicitly dealing with formal notations and mathematical expressions (typically used for biochemical modelling), nevertheless being able to obtain simulation and analysis results. We present the model and the sensitivity analysis results we have obtained, that allow us to identify the parameters which are most sensitive to perturbations. The results, which are shown to be in agreement with existing mathematical models of the gp130/JAK/STAT pathway, serve us as a form of validation of the model and of the approach itself