Direct estimation of kinetic parametric images for dynamic PET.

Dynamic positron emission tomography (PET) can monitor spatiotemporal distribution of radiotracer in vivo. The spatiotemporal information can be used to estimate parametric images of radiotracer kinetics that are of physiological and biochemical interests. Direct estimation of parametric images from raw projection data allows accurate noise modeling and has been shown to offer better image quality than conventional indirect methods, which reconstruct a sequence of PET images first and then perform tracer kinetic modeling pixel-by-pixel. Direct reconstruction of parametric images has gained increasing interests with the advances in computing hardware. Many direct reconstruction algorithms have been developed for different kinetic models. In this paper we review the recent progress in the development of direct reconstruction algorithms for parametric image estimation. Algorithms for linear and nonlinear kinetic models are described and their properties are discussed

Inverse Problems in data-driven multi-scale Systems Medicine: application to cancer physiology

Systems Medicine is an interdisciplinary framework involving reciprocal feedback between clinical investigation and mathematical modeling/analysis. Its aim is to improve the understanding of complex diseases by integrating knowledge and data across multiple levels of biological organization. This Thesis focuses on three inverse problems, arising from three kinds of data and related to cancer physiology, at different scales: tissues, cells, molecules. The general assumption of this piece of research is that cancer is associated toa path ological glucose consumption and, in fact, its functional behavior can be assessed by nuclear medicine experiments using [18F]-fluorodeoxyglucose (FDG) as a radioactive tracer mimicking the glucose properties. At tissue-scale, this Thesis considers the Positron Emission Tomography (PET) imaging technique, and deals with two distinct issues within compartmental analysis. First, this Thesis presents a compartmental approach, referred to as reference tissue model, for the estimation of FDG kinetics inside cancer tissues when the arterial blood input of the system is unknown. Then, this Thesis proposes an efficient and reliable method for recovering the compartmental kinetic parameters for each PET image pixel in the context of parametric imaging, exploiting information on the tissue physiology. Standard models in compartmental analysis assume that phosphorylation and dephosphorylation of FDG occur in the same intracellular cytosolic volume. Advances in cell biochemistry have shown that the appropriate location of dephosphorylation is the endoplasmic reticulum (ER). Therefore, at cell-scale, this Thesis formalizes a biochemically-driven compartmental model accounting for the specific role played by the ER, and applies it to the analysis of in vitro experiments on FDG uptake by cancer cell cultures obtained with a LigandTracer (LT) device. Finally, at molecule-scale, this Thesis provides a preliminary mathematical investigation of a chemical reaction network (CRN), represented by a huge Molecular Interaction Map (MIM), describing the biochemical interactions occurring between signaling proteins in specific pathways within a cancer cell. The main issue addressed in this case is the network parameterization problem, i.e. how to determine the reaction rate coefficients from protein concentration data

Investigation of the Effects of Image Signal-to-Noise Ratio on TSPO PET Quantification of Neuroinflammation

Neuroinflammation may be imaged using positron emission tomography (PET) and the tracer [11C]-PK11195. Accurate and precise quantification of 18 kilodalton Translocator Protein (TSPO) binding parameters in the brain has proven difficult with this tracer, due to an unfavourable combination of low target concentration in tissue, low brain uptake of the tracer and relatively high non-specific binding, all of which leads to higher levels of relative image noise. To address these limitations, research into new radioligands for the TSPO, with higher brain uptake and lower non-specific binding relative to [11C]-PK11195, is being conducted world-wide. However, factors other than radioligand properties are known to influence signal-to-noise ratio in quantitative PET studies, including the scanner sensitivity, image reconstruction algorithms and data analysis methodology. The aim of this thesis was to investigate and validate computational tools for predicting image noise in dynamic TSPO PET studies, and to employ those tools to investigate the factors that affect image SNR and reliability of TSPO quantification in the human brain. The feasibility of performing multiple (n≥40) independent Monte Carlo simulations for each dynamic [11C]-PK11195 frame- with realistic modelling of the radioactivity source, attenuation and PET tomograph geometries- was investigated. A Beowulf-type high performance computer cluster, constructed from commodity components, was found to be well suited to this task. Timing tests on a single desktop computer system indicated that a computer cluster capable of simulating an hour-long dynamic [11C]-PK11195 PET scan, with 40 independent repeats, and with a total simulation time of less than 6 weeks, could be constructed for less than 10,000 Australian dollars. A computer cluster containing 44 computing cores was therefore assembled, and a peak simulation rate of 2.84x105 photon pairs per second was achieved using the GEANT4 Application for Tomographic Emission (GATE) Monte Carlo simulation software. A simulated PET tomograph was developed in GATE that closely modelled the performance characteristics of several real-world clinical PET systems in terms of spatial resolution, sensitivity, scatter fraction and counting rate performance. The simulated PET system was validated using adaptations of the National Electrical Manufacturers Association (NEMA) quality assurance procedures within GATE. Image noise in dynamic TSPO PET scans was estimated by performing n=40 independent Monte Carlo simulations of an hour-long [11C]-PK11195 scan, and of an hour- long dynamic scan for a hypothetical TSPO ligand with double the brain activity concentration of [11C]-PK11195. From these data an analytical noise model was developed that allowed image noise to be predicted for any combination of brain tissue activity concentration and scan duration. The noise model was validated for the purpose of determining the precision of kinetic parameter estimates for TSPO PET. An investigation was made into the effects of activity concentration in tissue, radionuclide half-life, injected dose and compartmental model complexity on the reproducibility of kinetic parameters. Injecting 555 MBq of carbon-11 labelled TSPO tracer produced similar binding parameter precision to 185 MBq of fluorine-18, and a moderate (20%) reduction in precision was observed for the reduced carbon-11 dose of 370 MBq. Results indicated that a factor of 2 increase in frame count level (relative to [11C]-PK11195, and due for example to higher ligand uptake, injected dose or absolute scanner sensitivity) is required to obtain reliable binding parameter estimates for small regions of interest when fitting a two-tissue compartment, four-parameter compartmental model. However, compartmental model complexity had a similarly large effect, with the reduction of model complexity from the two-tissue compartment, four-parameter to a one-tissue compartment, two-parameter model producing a 78% reduction in coefficient of variation of the binding parameter estimates at each tissue activity level and region size studied. In summary, this thesis describes the development and validation of Monte Carlo methods for estimating image noise in dynamic TSPO PET scans, and analytical methods for predicting relative image noise for a wide range of tissue activity concentration and acquisition durations. The findings of this research suggest that a broader consideration of the kinetic properties of novel TSPO radioligands, with a view to selection of ligands that are potentially amenable to analysis with a simple one-tissue compartment model, is at least as important as efforts directed towards reducing image noise, such as higher brain uptake, in the search for the next generation of TSPO PET tracers

Regularized estimation and model selection in compartment models

Dynamic imaging series acquired in medical and biological research are often analyzed with the help of compartment models. Compartment models provide a parametric, nonlinear function of interpretable, kinetic parameters describing how some concentration of interest evolves over time. Aiming to estimate the kinetic parameters, this leads to a nonlinear regression problem. In many applications, the number of compartments needed in the model is not known from biological considerations but should be inferred from the data along with the kinetic parameters. As data from medical and biological experiments are often available in the form of images, the spatial data structure of the images has to be taken into account. This thesis addresses the problem of parameter estimation and model selection in compartment models. Besides a penalized maximum likelihood based approach, several Bayesian approaches-including a hierarchical model with Gaussian Markov random field priors and a model state approach with flexible model dimension-are proposed and evaluated to accomplish this task. Existing methods are extended for parameter estimation and model selection in more complex compartment models. However, in nonlinear regression and, in particular, for more complex compartment models, redundancy issues may arise. This thesis analyzes difficulties arising due to redundancy issues and proposes several approaches to alleviate those redundancy issues by regularizing the parameter space. The potential of the proposed estimation and model selection approaches is evaluated in simulation studies as well as for two in vivo imaging applications: a dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) study on breast cancer and a study on the binding behavior of molecules in living cell nuclei observed in a fluorescence recovery after photobleaching (FRAP) experiment

Measuring Inorganic Carbon Fluxes from Carbonate Mineral Weathering from Large River Basins: The Ohio River Basin

Rising atmospheric CO2 concentrations have motivated efforts to better quantify reservoirs and fluxes of Earth’s carbon. Of these fluxes from the atmosphere, one that has received relatively little attention is the atmospheric carbon sink associated with carbonate mineral dissolution. Osterhoudt (2014) and Salley (2016) explored new normalization techniques to improve and standardize a process for measuring this flux over large river basins. The present research extends this work to the 490,600 km2 Ohio River drainage basin and 11 subbasins. The study estimated the DIC flux leaving these basins between October 1, 2013, and September 30, 2014, based on secondary hydrogeochemical, geologic, and climatic data. The total annual DIC flux for the Ohio River basin was estimated to be 7.54 x 1012 g carbon (C). The time-volume normalized value of DIC flux for the Ohio basin was 3.36 x 108 g C/km3 day, where the km3 refers to the amount of water available during the year. This was within 71.4% agreement with the Barren River data (Salley, 2016) and within 63.9% agreement with the Green River data (Osterhoudt, 2014). In general, normalized DIC flux values of sub-basins containing at least modest amounts (more than 8%) of exposed carbonates (Tennessee, Cumberland, Green, Kentucky, Licking, Monongahela, and Allegheny) were in strong agreement with the normalized DIC flux of the Ohio River basin, whereas inclusion of basins with little or no near surface carbonates (Wabash, Great Miami, Scioto and Kanawha) yielded poor agreement. Regression analysis yielded strong agreement between DIC flux and the normalization parameters for the carbonate-bearing sub-basins (R2 = 0.97, p