Model of Host-Pathogen Interaction Dynamics Links In Vivo Optical Imaging and Immune Responses.

Tracking disease progression in vivo is essential for the development of treatments against bacterial infection. Optical imaging has become a central tool for in vivo tracking of bacterial population development and therapeutic response. For a precise understanding of in vivo imaging results in terms of disease mechanisms derived from detailed postmortem observations, however, a link between the two is needed. Here, we develop a model that provides that link for the investigation of Citrobacter rodentium infection, a mouse model for enteropathogenic Escherichia coli (EPEC). We connect in vivo disease progression of C57BL/6 mice infected with bioluminescent bacteria, imaged using optical tomography and X-ray computed tomography, to postmortem measurements of colonic immune cell infiltration. We use the model to explore changes to both the host immune response and the bacteria and to evaluate the response to antibiotic treatment. The developed model serves as a novel tool for the identification and development of new therapeutic interventions.

MEANS: python package for Moment Expansion Approximation, iNference and Simulation.

MOTIVATION: Many biochemical systems require stochastic descriptions. Unfortunately these can only be solved for the simplest cases and their direct simulation can become prohibitively expensive, precluding thorough analysis. As an alternative, moment closure approximation methods generate equations for the time-evolution of the system's moments and apply a closure ansatz to obtain a closed set of differential equations; that can become the basis for the deterministic analysis of the moments of the outputs of stochastic systems. RESULTS: We present a free, user-friendly tool implementing an efficient moment expansion approximation with parametric closures that integrates well with the IPython interactive environment. Our package enables the analysis of complex stochastic systems without any constraints on the number of species and moments studied and the type of rate laws in the system. In addition to the approximation method our package provides numerous tools to help non-expert users in stochastic analysis. AVAILABILITY AND IMPLEMENTATION: https://github.com/theosysbio/means CONTACTS: m.stumpf@imperial.ac.uk or e.lakatos13@imperial.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Early recognition of lung cancer by integrin targeted imaging in K-ras mouse model.

Non-small cell lung cancer is characterized by slow progression and high heterogeneity of tumors. Integrins play an important role in lung cancer development and metastasis and were suggested as a tumor marker; however their role in anticancer therapy remains controversial. In this work, we demonstrate the potential of integrin-targeted imaging to recognize early lesions in transgenic mouse model of lung cancer based on spontaneous introduction of mutated human gene bearing K-ras mutation. We conducted ex vivo and fluorescence molecular tomography-X-ray computed tomography (FMT-XCT) in vivo imaging and analysis for specific targeting of early lung lesions and tumors in rodent preclinical model for lung cancer. The lesions and tumors were characterized by histology, immunofluorescence and immunohistochemistry using a panel of cancer markers. Ex vivo, the integrin-targeted fluorescent signal significantly differed between wild type lung tissue and K-ras pulmonary lesions (PL) at all ages studied. The panel of immunofluorescence experiments demonstrated that PL, which only partially show cancer cell features were detected by αvß3-integrin targeted imaging. Human patient material analysis confirmed the specificity of target localization in different lung cancer types. Most importantly, small tumors in the lungs of 4-week-old animals could be noninvasively detected in vivo on the fluorescence channel of FMT-XCT. Our findings demonstrated αvß3-integrin targeted fluorescent imaging to specifically detect premalignant pleural lesions in K-ras mice. Integrin targeted imaging may find application areas in preclinical research and clinical practice, such as early lung cancer diagnostics, intraoperative assistance or therapy monitoring.

FMT-XCT: in vivo animal studies with hybrid fluorescence molecular tomography-X-ray computed tomography.

The development of hybrid optical tomography methods to improve imaging performance has been suggested over a decade ago and has been experimentally demonstrated in animals and humans. Here we examined in vivo performance of a camera-based hybrid fluorescence molecular tomography (FMT) system for 360° imaging combined with X-ray computed tomography (XCT). Offering an accurately co-registered, information-rich hybrid data set, FMT-XCT has new imaging possibilities compared to stand-alone FMT and XCT. We applied FMT-XCT to a subcutaneous 4T1 tumor mouse model, an Aga2 osteogenesis imperfecta model and a Kras lung cancer mouse model, using XCT information during FMT inversion. We validated in vivo imaging results against post-mortem planar fluorescence images of cryoslices and histology data. Besides offering concurrent anatomical and functional information, FMT-XCT resulted in the most accurate FMT performance to date. These findings indicate that addition of FMT optics into the XCT gantry may be a potent upgrade for small-animal XCT systems.

Multivariate moment closure techniques for stochastic kinetic models.

Stochastic effects dominate many chemical and biochemical processes. Their analysis, however, can be computationally prohibitively expensive and a range of approximation schemes have been proposed to lighten the computational burden. These, notably the increasingly popular linear noise approximation and the more general moment expansion methods, perform well for many dynamical regimes, especially linear systems. At higher levels of nonlinearity, it comes to an interplay between the nonlinearities and the stochastic dynamics, which is much harder to capture correctly by such approximations to the true stochastic processes. Moment-closure approaches promise to address this problem by capturing higher-order terms of the temporally evolving probability distribution. Here, we develop a set of multivariate moment-closures that allows us to describe the stochastic dynamics of nonlinear systems. Multivariate closure captures the way that correlations between different molecular species, induced by the reaction dynamics, interact with stochastic effects. We use multivariate Gaussian, gamma, and lognormal closure and illustrate their use in the context of two models that have proved challenging to the previous attempts at approximating stochastic dynamics: oscillations in p53 and Hes1. In addition, we consider a larger system, Erk-mediated mitogen-activated protein kinases signalling, where conventional stochastic simulation approaches incur unacceptably high computational costs.

Ex-vivo assessment and non-invasive in vivo imaging of internal hemorrhages in Aga2/+ mutant mice.

Mutations in type I collagen genes (COL1A1/2) typically lead to Osteogenesis imperfecta, the most common heritable cause of skeletal fractures and bone deformation in humans. Heterozygous Col1a1(Aga2/+), animals with a dominant mutation in the terminal C-propeptide domain of type I collagen develop typical skeletal hallmarks and internal hemorrhages starting from 6 day after birth. The disease progression for Aga2/+ mice, however, is not uniform differing between severe phenotype lethal at the 6-11th day of life, and moderate-to-severe one with survival to adulthood. Herein we investigated whether a new modality that combines X-ray computer tomography with fluorescence tomography in one hybrid system can be employed to study internal bleedings in relation to bone fractures and obtain insights into disease progression. The disease phenotype was characterized on Aga2/+ vs. wild type mice between 6 and 9 days postnatal. Anatomical and functional findings obtained in-vivo were contrasted to the ex-vivo appearance of the same tissues under cryo-slicing.

A general moment expansion method for stochastic kinetic models.

Moment approximation methods are gaining increasing attention for their use in the approximation of the stochastic kinetics of chemical reaction systems. In this paper we derive a general moment expansion method for any type of propensities and which allows expansion up to any number of moments. For some chemical reaction systems, more than two moments are necessary to describe the dynamic properties of the system, which the linear noise approximation is unable to provide. Moreover, also for systems for which the mean does not have a strong dependence on higher order moments, moment approximation methods give information about higher order moments of the underlying probability distribution. We demonstrate the method using a dimerisation reaction, Michaelis-Menten kinetics and a model of an oscillating p53 system. We show that for the dimerisation reaction and Michaelis-Menten enzyme kinetics system higher order moments have limited influence on the estimation of the mean, while for the p53 system, the solution for the mean can require several moments to converge to the average obtained from many stochastic simulations. We also find that agreement between lower order moments does not guarantee that higher moments will agree. Compared to stochastic simulations, our approach is numerically highly efficient at capturing the behaviour of stochastic systems in terms of the average and higher moments, and we provide expressions for the computational cost for different system sizes and orders of approximation. We show how the moment expansion method can be employed to efficiently quantify parameter sensitivity. Finally we investigate the effects of using too few moments on parameter estimation, and provide guidance on how to estimate if the distribution can be accurately approximated using only a few moments.

Preconditioning of the fluorescence diffuse optical tomography sensing matrix based on compressive sensing.

Image reconstruction in fluorescence diffuse optical tomography (FDOT) is a highly ill-posed inverse problem due to a large number of unknowns and limited measurements. In FDOT, the fluorophore distribution is often sparse in the imaging domain, since most fluorophores are designed to accumulate in relatively small regions. Compressive sensing theory has shown that sparse signals can be recovered exactly from only a small number of measurements when the forward sensing matrix is sufficiently incoherent. In this Letter, we present a method of preconditioning the FDOT forward matrix to reduce its coherence. The reconstruction results using real data obtained from a phantom experiment show visual and quantitative improvements due to preconditioning in conjunction with convex relaxation and greedy-type sparse signal recovery algorithms.

Revisiting the normalized Born approximation: effects of scattering.

The normalized Born approximation has been suggested as a ratiometric method in fluorescence molecular tomography (FMT) applications, to account for heterogeneity variations. The method enabled practical inversions, as it offered fluorescence reconstruction accuracy over a wide range of absorption heterogeneity, while also accounting for unknown experimental factors, such as the various system gains and losses. Yet it was noted that scattering variations affect the robustness and accuracy. Herein we decompose the effects of absorption and scattering and capitalize on the recent development of hybrid FMT/x-ray computed tomography imaging methods to proposed amendments to the method, which improve the overall accuracy of the approach.

Performance evaluation of adaptive meshing algorithms for fluorescence diffuse optical tomography using experimental data.

Fluorescence diffuse optical tomography (FDOT) is a computationally demanding imaging problem. The discretizations of FDOT forward and inverse problems pose a trade-off between the accuracy and the computational efficiency of the image reconstruction. To address this trade-off, we analyzed the effect of discretization on the accuracy of FDOT imaging and proposed novel adaptive meshing algorithms for FDOT in a series of studies. In this Letter, we apply these new adaptive meshing algorithms to FDOT imaging using real data from a phantom experiment to demonstrate the practical advantages of our algorithms in FDOT image reconstruction.

Imaging performance of a hybrid x-ray computed tomography-fluorescence molecular tomography system using priors.

PURPOSE: The performance is studied of two newly introduced and previously suggested methods that incorporate priors into inversion schemes associated with data from a recently developed hybrid x-ray computed tomography and fluorescence molecular tomography system, the latter based on CCD camera photon detection. The unique data set studied attains accurately registered data of high spatially sampled photon fields propagating through tissue along 360 degrees projections. METHODS: Approaches that incorporate structural prior information were included in the inverse problem by adding a penalty term to the minimization function utilized for image reconstructions. Results were compared as to their performance with simulated and experimental data from a lung inflammation animal model and against the inversions achieved when not using priors. RESULTS: The importance of using priors over stand-alone inversions is also showcased with high spatial sampling simulated and experimental data. The approach of optimal performance in resolving fluorescent biodistribution in small animals is also discussed. CONCLUSIONS: Inclusion of prior information from x-ray CT data in the reconstruction of the fluorescence biodistribution leads to improved agreement between the reconstruction and validation images for both simulated and experimental data.

Patch-based anisotropic diffusion scheme for fluorescence diffuse optical tomography--part 2: image reconstruction.

Fluorescence diffuse optical tomography (fDOT) provides 3D images of fluorescence distributions in biological tissue, which represent molecular and cellular processes. The image reconstruction problem is highly ill-posed and requires regularisation techniques to stabilise and find meaningful solutions. Quadratic regularisation tends to either oversmooth or generate very noisy reconstructions, depending on the regularisation strength. Edge preserving methods, such as anisotropic diffusion regularisation (AD), can preserve important features in the fluorescence image and smooth out noise. However, AD has limited ability to distinguish an edge from noise. We propose a patch-based anisotropic diffusion regularisation (PAD), where regularisation strength is determined by a weighted average according to the similarity between patches around voxels within a search window, instead of a simple local neighbourhood strategy. However, this method has higher computational complexity and, hence, we wavelet compress the patches (PAD-WT) to speed it up, while simultaneously taking advantage of the denoising properties of wavelet thresholding. Furthermore, structural information can be incorporated into the image reconstruction with PAD-WT to improve image quality and resolution. In this case, the weights used to average voxels in the image are calculated using the structural image, instead of the fluorescence image. The regularisation strength depends on both structural and fluorescence images, which guarantees that the method can preserve fluorescence information even when it is not structurally visible in the anatomical images. In part 1, we tested the method using a denoising problem. Here, we use simulated and in vivo mouse fDOT data to assess the algorithm performance. Our results show that the proposed PAD-WT method provides high quality and noise free images, superior to those obtained using AD.

Phosphorelay of non-orthodox two component systems functions through a bi-molecular mechanism in vivo: the case of ArcB.

Two-component systems play a central part in bacterial signal transduction. Phosphorelay mechanisms have been linked to more robust and ultra-sensitive signalling dynamics. The molecular machinery that facilitates such a signalling is, however, only understood in outline. In particular the functional relevance of the dimerization of a non-orthodox or hybrid histidine kinase along which the phosphorelay takes place has been a subject of debate. We use a combination of molecular and genetic approaches, coupled to mathematical and statistical modelling, to demonstrate that the different possible intra- and inter-molecular mechanisms of phosphotransfer are formally non-identifiable in Escherichia coli expressing the ArcB non-orthodox histidine kinase used in anoxic redox control. In order to resolve this issue we further analyse the mathematical model in order to identify discriminatory experiments, which are then performed to address cis- and trans-phosphorelay mechanisms. The results suggest that exclusive cis- and trans-mechanisms will not be operating, instead the functional phosphorelay is likely to build around a sequence of allosteric interactions among the domain pairs in the histidine kinase. This is the first detailed mechanistic analysis of the molecular processes involved in non-orthodox two-component signalling and our results suggest strongly that dimerization facilitates more discriminatory proof-reading of external signals, via these allosteric reactions, prior to them being further processed.

Deep-tissue reporter-gene imaging with fluorescence and optoacoustic tomography: a performance overview.

PURPOSE: A primary enabling feature of near-infrared fluorescent proteins (FPs) and fluorescent probes is the ability to visualize deeper in tissues than in the visible. The purpose of this work is to find which is the optimal visualization method that can exploit the advantages of this novel class of FPs in full-scale pre-clinical molecular imaging studies. PROCEDURES: Nude mice were stereotactically implanted with near-infrared FP expressing glioma cells to from brain tumors. The feasibility and performance metrics of FPs were compared between planar epi-illumination and trans-illumination fluorescence imaging, as well as to hybrid Fluorescence Molecular Tomography (FMT) system combined with X-ray CT and Multispectral Optoacoustic (or Photoacoustic) Tomography (MSOT). RESULTS: It is shown that deep-seated glioma brain tumors are possible to visualize both with fluorescence and optoacoustic imaging. Fluorescence imaging is straightforward and has good sensitivity; however, it lacks resolution. FMT-XCT can provide an improved rough resolution of ∼1 mm in deep tissue, while MSOT achieves 0.1 mm resolution in deep tissue and has comparable sensitivity. CONCLUSIONS: We show imaging capacity that can shift the visualization paradigm in biological discovery. The results are relevant not only to reporter gene imaging, but stand as cross-platform comparison for all methods imaging near infrared fluorescent contrast agents.

Fluorescence background subtraction technique for hybrid fluorescence molecular tomography/x-ray computed tomography imaging of a mouse model of early stage lung cancer.

The ability to visualize early stage lung cancer is important in the study of biomarkers and targeting agents that could lead to earlier diagnosis. The recent development of hybrid free-space 360-deg fluorescence molecular tomography (FMT) and x-ray computed tomography (XCT) imaging yields a superior optical imaging modality for three-dimensional small animal fluorescence imaging over stand-alone optical systems. Imaging accuracy was improved by using XCT information in the fluorescence reconstruction method. Despite this progress, the detection sensitivity of targeted fluorescence agents remains limited by nonspecific background accumulation of the fluorochrome employed, which complicates early detection of murine cancers. Therefore we examine whether x-ray CT information and bulk fluorescence detection can be combined to increase detection sensitivity. Correspondingly, we research the performance of a data-driven fluorescence background estimator employed for subtraction of background fluorescence from acquisition data. Using mice containing known fluorochromes ex vivo, we demonstrate the reduction of background signals from reconstructed images and sensitivity improvements. Finally, by applying the method to in vivo data from K-ras transgenic mice developing lung cancer, we find small tumors at an early stage compared with reconstructions performed using raw data. We conclude with the benefits of employing fluorescence subtraction in hybrid FMT-XCT for early detection studies.

Cardioprotective C-kit⁺ bone marrow cells attenuate apoptosis after acute myocardial infarction in mice - in-vivo assessment with fluorescence molecular imaging.

Cardiomyocyte loss via apoptosis plays a crucial role in ventricular remodeling following myocardial infarction (MI). Cell-based therapy approaches using bone marrow derived c-kit⁺ pluripotent cells may attenuate apoptosis following ischemic injury. We therefore thought to examine the early course of apoptosis following myocardial infarction - in-vivo - and non-invasively determine the effect of c-kit⁺ bone marrow cells on post-MI remodeling. We studied apoptosis in wild-type Kit(+/+) , c-kit mutant Kit(W)/Kit(W-v) and Kit(W)/Kit(W-v) mice after cell therapy with bone-marrow derived c-kit⁺ cells after ischemia-reperfusion injury. Mice were followed by hybrid Fluorescence Molecular Tomography/X-ray Computed Tomography (FMT-XCT) at 6h, 24h and 7 days after ischemia-reperfusion injury using an Annexin V-based fluorescent nanosensor targeting phosphatidylserine. Kit(W)/Kit(W-v) mice showed increased and prolonged apoptosis compared to control Kit(+/+) mice while c-kit cell therapy was able to attenuate the altered apoptosis rates. Increased apoptosis was accompanied by severe decline in heart function, determined by cardiac Magnetic Resonance Imaging, and cell therapy was able to rescue the animals from deleterious heart failure. Post-mortem cryoslicing and immunohistochemistry localized the fluorescence signal of the Annexin V sensor within the infarcted myocardium. Flow cytometry of digested infarct specimens identified apoptotic cardiomyocytes as the major source for the in-vivo Annexin V signal. In-vivo molecular imaging using hybrid FMT-XCT reveals increased cardiomyocyte apoptosis in Kit(W)/Kit(W-v) mice and shows that c-kit⁺ cardioprotective cells are able to attenuate post-MI apoptosis and rescue mice from progressive heart failure.

Improving limited-projection-angle fluorescence molecular tomography using a co-registered x-ray computed tomography scan.

We examine the improvement in imaging performance, such as axial resolution and signal localization, when employing limited-projection-angle fluorescence molecular tomography (FMT) together with x-ray computed tomography (XCT) measurements versus stand-alone FMT. For this purpose, we employed living mice, bearing a spontaneous lung tumor model, and imaged them with FMT and XCT under identical geometrical conditions using fluorescent probes for cancer targeting. The XCT data was employed, herein, as structural prior information to guide the FMT reconstruction. Gold standard images were provided by fluorescence images of mouse cryoslices, providing the ground truth in fluorescence bio-distribution. Upon comparison of FMT images versus images reconstructed using hybrid FMT and XCT data, we demonstrate marked improvements in image accuracy. This work relates to currently disseminated FMT systems, using limited projection scans, and can be employed to enhance their performance.

Hybrid system for simultaneous fluorescence and x-ray computed tomography.

A hybrid imaging system for simultaneous fluorescence tomography and X-ray computed tomography (XCT) of small animals has been developed and presented. The system capitalizes on the imaging power of a 360 ( degrees )-projection free-space fluorescence tomography system, implemented within a microcomputed tomography scanner. Image acquisition is based on techniques that automatically adjust a series of imaging parameters to offer a high dynamic range dataset. Image segmentation further allows the incorporation of structural priors in the optical reconstruction problem to improve the imaging performance. The functional system characteristics are showcased, and images from a brain imaging study are shown, which are reconstructed using XCT-derived priors into the optical forward problem.

Fast automatic segmentation of anatomical structures in x-ray computed tomography images to improve fluorescence molecular tomography reconstruction.

The recent development of hybrid imaging scanners that integrate fluorescence molecular tomography (FMT) and x-ray computed tomography (XCT) allows the utilization of x-ray information as image priors for improving optical tomography reconstruction. To fully capitalize on this capacity, we consider a framework for the automatic and fast detection of different anatomic structures in murine XCT images. To accurately differentiate between different structures such as bone, lung, and heart, a combination of image processing steps including thresholding, seed growing, and signal detection are found to offer optimal segmentation performance. The algorithm and its utilization in an inverse FMT scheme that uses priors is demonstrated on mouse images.