Laplace deconvolution on the basis of time domain data and its application to Dynamic Contrast Enhanced imaging

In the present paper we consider the problem of Laplace deconvolution with noisy discrete non-equally spaced observations on a finite time interval. We propose a new method for Laplace deconvolution which is based on expansions of the convolution kernel, the unknown function and the observed signal over Laguerre functions basis (which acts as a surrogate eigenfunction basis of the Laplace convolution operator) using regression setting. The expansion results in a small system of linear equations with the matrix of the system being triangular and Toeplitz. Due to this triangular structure, there is a common number $m$ of terms in the function expansions to control, which is realized via complexity penalty. The advantage of this methodology is that it leads to very fast computations, produces no boundary effects due to extension at zero and cut-off at $T$ and provides an estimator with the risk within a logarithmic factor of the oracle risk. We emphasize that, in the present paper, we consider the true observational model with possibly nonequispaced observations which are available on a finite interval of length $T$ which appears in many different contexts, and account for the bias associated with this model (which is not present when $T\rightarrow\infty$). The study is motivated by perfusion imaging using a short injection of contrast agent, a procedure which is applied for medical assessment of micro-circulation within tissues such as cancerous tumors. Presence of a tuning parameter $a$ allows to choose the most advantageous time units, so that both the kernel and the unknown right hand side of the equation are well represented for the deconvolution. The methodology is illustrated by an extensive simulation study and a real data example which confirms that the proposed technique is fast, efficient, accurate, usable from a practical point of view and very competitive.Comment: 36 pages, 9 figures. arXiv admin note: substantial text overlap with arXiv:1207.223

Astigmatism and Pseudoaccommodation in Pseudophakic Eyes

noAdvanced IOLs with circumferential zones of different power provide pseudoaccommodation. We investigated the potential for power variation with meridian, namely astigmatism, to provide pseudo-accommodation. With appropriate power and axis orientations, acceptable pseudo-accommodation can be achieved

Collagen microarchitecture mechanically controls myofibroblast differentiation.

Altered microarchitecture of collagen type I is a hallmark of wound healing and cancer that is commonly attributed to myofibroblasts. However, it remains unknown which effect collagen microarchitecture has on myofibroblast differentiation. Here, we combined experimental and computational approaches to investigate the hypothesis that the microarchitecture of fibrillar collagen networks mechanically regulates myofibroblast differentiation of adipose stromal cells (ASCs) independent of bulk stiffness. Collagen gels with controlled fiber thickness and pore size were microfabricated by adjusting the gelation temperature while keeping their concentration constant. Rheological characterization and simulation data indicated that networks with thicker fibers and larger pores exhibited increased strain-stiffening relative to networks with thinner fibers and smaller pores. Accordingly, ASCs cultured in scaffolds with thicker fibers were more contractile, expressed myofibroblast markers, and deposited more extended fibronectin fibers. Consistent with elevated myofibroblast differentiation, ASCs in scaffolds with thicker fibers exhibited a more proangiogenic phenotype that promoted endothelial sprouting in a contractility-dependent manner. Our findings suggest that changes of collagen microarchitecture regulate myofibroblast differentiation and fibrosis independent of collagen quantity and bulk stiffness by locally modulating cellular mechanosignaling. These findings have implications for regenerative medicine and anticancer treatments

A Semi-parametric Technique for the Quantitative Analysis of Dynamic Contrast-enhanced MR Images Based on Bayesian P-splines

Dynamic Contrast-enhanced Magnetic Resonance Imaging (DCE-MRI) is an important tool for detecting subtle kinetic changes in cancerous tissue. Quantitative analysis of DCE-MRI typically involves the convolution of an arterial input function (AIF) with a nonlinear pharmacokinetic model of the contrast agent concentration. Parameters of the kinetic model are biologically meaningful, but the optimization of the non-linear model has significant computational issues. In practice, convergence of the optimization algorithm is not guaranteed and the accuracy of the model fitting may be compromised. To overcome this problems, this paper proposes a semi-parametric penalized spline smoothing approach, with which the AIF is convolved with a set of B-splines to produce a design matrix using locally adaptive smoothing parameters based on Bayesian penalized spline models (P-splines). It has been shown that kinetic parameter estimation can be obtained from the resulting deconvolved response function, which also includes the onset of contrast enhancement. Detailed validation of the method, both with simulated and in vivo data, is provided

Data-driven modelling of biological multi-scale processes

Biological processes involve a variety of spatial and temporal scales. A holistic understanding of many biological processes therefore requires multi-scale models which capture the relevant properties on all these scales. In this manuscript we review mathematical modelling approaches used to describe the individual spatial scales and how they are integrated into holistic models. We discuss the relation between spatial and temporal scales and the implication of that on multi-scale modelling. Based upon this overview over state-of-the-art modelling approaches, we formulate key challenges in mathematical and computational modelling of biological multi-scale and multi-physics processes. In particular, we considered the availability of analysis tools for multi-scale models and model-based multi-scale data integration. We provide a compact review of methods for model-based data integration and model-based hypothesis testing. Furthermore, novel approaches and recent trends are discussed, including computation time reduction using reduced order and surrogate models, which contribute to the solution of inference problems. We conclude the manuscript by providing a few ideas for the development of tailored multi-scale inference methods.Comment: This manuscript will appear in the Journal of Coupled Systems and Multiscale Dynamics (American Scientific Publishers

Cylindrical illumination with angular coupling for whole-prostate photoacoustic tomography

Current diagnosis of prostate cancer relies on histological analysis of tissue samples acquired by biopsy, which could benefit from real-time identification of suspicious lesions. Photoacoustic tomography has the potential to provide real-time targets for prostate biopsy guidance with chemical selectivity, but light delivered from the rectal cavity has been unable to penetrate to the anterior prostate. To overcome this barrier, a urethral device with cylindrical illumination is developed for whole-prostate imaging, and its performance as a function of angular light coupling is evaluated with a prostate-mimicking phantom

Aerospace medicine and biology: A continuing bibliography with indexes, supplement 162, January 1977

This bibliography lists 189 reports, articles, and other documents introduced into the NASA scientific and technical information system in December 1976

Quantification of venous blood signal contribution to BOLD functional activation in the auditory cortex at 3 T

Most modern techniques for functional magnetic resonance imaging (fMRI) rely on blood-oxygen-level-dependent (BOLD) contrast as the basic principle for detecting neuronal activation. However, the measured BOLD effect depends on a transfer function related to neurophysiological changes accompanying electrical neural activation. The spatial accuracy and extension of the region of interest are determined by vascular effect, which introduces incertitude on real neuronal activation maps. Our efforts have been directed towards the development of a new methodology that is capable of combining morphological, vascular and functional information; obtaining new insight regarding foci of activation; and distinguishing the nature of activation on a pixel-by-pixel basis. Six healthy volunteers were studied in a parametric auditory functional experiment at 3 T; activation maps were overlaid on a high-resolution brain venography obtained through a novel technique. The BOLD signal intensities of vascular and nonvascular activated voxels were analyzed and compared: it was shown that nonvascular active voxels have lower values for signal peak (Pb10−7) and area (Pb10−8) with respect to vascular voxels. The analysis showed how venous blood influenced the measured BOLD signals, supplying a technique to filter possible venous artifacts that potentially can lead to misinterpretation of fMRI results. This methodology, although validated in the auditory cortex activation, maintains a general applicability to any cortical fMRI study, as the basic concepts on which it relies on are not limited to this cortical region. The results obtained in this study can represent the basis for new methodologies and tools that are capable of adding further characterization to the BOLD signal properties