Non-invasive estimation of hepatic blood perfusion from H2 15O PET images using tissue-derived arterial and portal input functions

Quantitative Perfusion Imaging with Total-Body PET

Recently, PET systems with a long axial field of view have become the current state of the art. Total-body PET scanners enable unique possibilities for scientific research and clinical diagnostics, but this new technology also raises numerous challenges. A key advantage of total-body imaging is that having all the organs in the field of view allows studying biologic interaction of all organs simultaneously. One of the new, promising imaging techniques is total-body quantitative perfusion imaging. Currently, 15O-labeled water provides a feasible option for quantitation of tissue perfusion at the total-body level. This review summarizes the status of the methodology and the analysis and provides examples of preliminary findings on applications of quantitative parametric perfusion images for research and clinical work. We also describe the opportunities and challenges arising from moving from single-organ studies to modeling of a multisystem approach with total-body PET, and we discuss future directions for total-body imaging.

[15O]H2O PET: Potential or Essential for Molecular Imaging?

Imaging water pathways in the human body provides an excellent way of measuring accurately the blood flow directed to different organs. This makes it a powerful diagnostic tool for a wide range of diseases that are related to perfusion and oxygenation. Although water PET has a long history, its true potential has not made it into regular clinical practice. The article highlights the potential of water PET in molecular imaging and suggests its prospective role in becoming an essential tool for the 21st century precision medicine in different domains ranging from preclinical to clinical research and practice. The recent technical advances in high-sensitivity PET imaging can play a key accelerating role in empowering this technique, though there are still several challenges to overcome.

Hepatic Positron Emission Tomography: Applications in Metabolism, Haemodynamics and Cancer

Evaluating in vivo the metabolic rates of the human liver has been a challenge due to its unique perfusion system. Positron emission tomography (PET) represents the current gold standard for assessing non-invasively tissue metabolic rates in vivo. Here, we review the existing literature on the assessment of hepatic metabolism, haemodynamics and cancer with PET. The tracer mainly used in metabolic studies has been [¹⁸F]2-fluoro-2-deoxy-D-glucose (¹⁸F-FDG). Its application not only enables the evaluation of hepatic glucose uptake in a variety of metabolic conditions and interventions, but based on the kinetics of ¹⁸F-FDG, endogenous glucose production can also be assessed. 14(R,S)-[¹⁸F]fluoro-6-thia-Heptadecanoic acid (¹⁸F-FTHA), ¹¹C-Palmitate and ¹¹C-Acetate have also been applied for the assessment of hepatic fatty acid uptake rates (¹⁸F-FTHA and ¹¹C-Palmitate) and blood flow and oxidation (¹¹C-Acetate). Oxygen-15 labelled water (¹⁵O-H₂O) has been used for the quantification of hepatic perfusion. ¹⁸F-FDG is also the most common tracer used for hepatic cancer diagnostics, whereas ¹¹C-Acetate has also shown some promising applications in imaging liver malignancies. The modelling approaches used to analyse PET data and also the challenges in utilizing PET in the assessment of hepatic metabolism are presented.

A continuum thermomechanical model of in vivo electrosurgical heating of hydrated soft biological tissues

Radio-frequency (RF) heating of soft biological tissues during electrosurgical procedures is a fast process that involves phase change through evaporation and transport of intra- and extra-cellular water, and where variations in physical properties with temperature and water content play significant role. Accurately predicting and capturing these effects would improve the modeling of temperature change in the tissue allowing the development of improved instrument design and better understanding of tissue damage and necrosis. Previous models based on the Pennes' bioheat model neglect both evaporation and transport or consider evaporation through numerical correlations, however, do not account for changes in physical properties due to mass transport or phase change, nor capture the pressure increase due to evaporation within the tissue. While a porous media approach can capture the effects of evaporation, transport, pressure and changes in physical properties, the model assumes free diffusion of liquid and gas without a careful examination of assumptions on transport parameters in intact tissue resulting in significant under prediction of temperature. These different approaches have therefore been associated with errors in temperature prediction exceeding 20% when compared to experiments due to inaccuracies in capturing the effects of evaporation losses and transport. Here, we present a model of RF heating of hydrated soft tissue based on mixture theory where the multiphase nature of tissue is captured within a continuum thermomechanics framework, simultaneously considering the transport, deformation and phase change losses due to evaporation that occur during electrosurgical heating. The model predictions are validated against data obtained for in vivo ablation of porcine liver tissue at various power settings of the electrosurgical unit. The model is able to match the mean experimental temperature data with sharp gradients in the vicinity of the electrode during rapid low and high power ablation procedures with errors less than 7.9%. Additionally, the model is able to capture fast vaporization losses and the corresponding increase in pressure due to vapor buildup which have a significant effect on temperature prediction beyond 100 °C.

Liver blood dynamics after bariatric surgery: the effects of mixed-meal test and incretin infusions

AIMS/HYPOTHESIS

The mechanisms for improved glycemic control after bariatric surgery in subjects with type 2 diabetes (T2D) are not fully known. We hypothesized that dynamic hepatic blood responses to a mixed-meal are changed after bariatric surgery in parallel with an improvement in glucose tolerance.

METHODS

A total of ten morbidly obese subjects with T2D were recruited to receive a mixed-meal and a glucose-dependent insulinotropic polypeptide (GIP) infusion before and early after (within a median of less than three months) bariatric surgery, and hepatic blood flow and volume (HBV) were measured repeatedly with combined positron emission tomography/MRI. Ten lean non-diabetic individuals served as controls.

RESULTS

Bariatric surgery leads to a significant decrease in weight, accompanied with an improved β-cell function and glucagon-like peptide 1 (GLP-1) secretion, and a reduction in liver volume. Blood flow in portal vein (PV) was increased by 1.65-fold (P = 0.026) in response to a mixed-meal in subjects after surgery, while HBV decreased in all groups (P < 0.001). When the effect of GIP infusion was tested separately, no change in hepatic arterial and PV flow was observed, but HBV decreased as seen during the mixed-meal test.

CONCLUSIONS/INTERPRETATION

Early after bariatric surgery, PV flow response to a mixed-meal is augmented, improving digestion and nutrient absorption. GIP influences the post-prandial reduction in HBV thereby diverting blood to the extrahepatic sites.

Reconstruction of input functions from a dynamic PET image with sequential administration of 15O2 and [Formula: see text] for noninvasive and ultra-rapid measurement of CBF, OEF, and CMRO2

CBF, OEF, and CMRO2 images can be quantitatively assessed using PET. Their image calculation requires arterial input functions, which require invasive procedure. The aim of the present study was to develop a non-invasive approach with image-derived input functions (IDIFs) using an image from an ultra-rapid O2 and C15O2 protocol. Our technique consists of using a formula to express the input using tissue curve with rate constants. For multiple tissue curves, the rate constants were estimated so as to minimize the differences of the inputs using the multiple tissue curves. The estimated rates were used to express the inputs and the mean of the estimated inputs was used as an IDIF. The method was tested in human subjects ( n = 24). The estimated IDIFs were well-reproduced against the measured ones. The difference in the calculated CBF, OEF, and CMRO2 values by the two methods was small (<10%) against the invasive method, and the values showed tight correlations ( r = 0.97). The simulation showed errors associated with the assumed parameters were less than ∼10%. Our results demonstrate that IDIFs can be reconstructed from tissue curves, suggesting the possibility of using a non-invasive technique to assess CBF, OEF, and CMRO2.

Digital liver biopsy: Bio-imaging of fatty liver for translational and clinical research

The rapidly growing field of functional, molecular and structural bio-imaging is providing an extraordinary new opportunity to overcome the limits of invasive liver biopsy and introduce a "digital biopsy" for in vivo study of liver pathophysiology. To foster the application of bio-imaging in clinical and translational research, there is a need to standardize the methods of both acquisition and the storage of the bio-images of the liver. It can be hoped that the combination of digital, liquid and histologic liver biopsies will provide an innovative synergistic tri-dimensional approach to identifying new aetiologies, diagnostic and prognostic biomarkers and therapeutic targets for the optimization of personalized therapy of liver diseases and liver cancer. A group of experts of different disciplines (Special Interest Group for Personalized Hepatology of the Italian Association for the Study of the Liver, Institute for Biostructures and Bio-imaging of the National Research Council and Bio-banking and Biomolecular Resources Research Infrastructure) discussed criteria, methods and guidelines for facilitating the requisite application of data collection. This manuscript provides a multi-Author review of the issue with special focus on fatty liver.

Increased Liver Fatty Acid Uptake Is Partly Reversed and Liver Fat Content Normalized After Bariatric Surgery

OBJECTIVE

Changes in liver fatty acid metabolism are important in understanding the mechanisms of diabetes remission and metabolic changes after bariatric surgery.

RESEARCH DESIGN AND METHODS

Liver fatty acid uptake (LFU), blood flow, and fat content (LFC) were measured in 25 obese subjects before bariatric surgery and 6 months after using positron emission tomography/computed tomography and MRS; 14 lean individuals served as the control subjects.

RESULTS

The increased LFU in obese subjects was associated with body adiposity. LFU was reduced postoperatively but was still high compared with the control subjects. LFC was normalized. Liver blood flow (per unit volume) was higher in obese subjects than in the control subjects at baseline and was further increased postoperatively; however, the total organ blood flow was unchanged as the liver volume decreased.

CONCLUSIONS

The findings suggest that in a postoperative state, intrahepatic fatty acids are not stored in the liver but are used for oxidation to provide energy. Changes in perfusion may contribute to improved liver metabolism postoperatively.

Application of adaptive kinetic modelling for bias propagation reduction in direct 4D image reconstruction

Parametric imaging in thoracic and abdominal PET can provide additional parameters more relevant to the pathophysiology of the system under study. However, dynamic data in the body are noisy due to the limiting counting statistics leading to suboptimal kinetic parameter estimates. Direct 4D image reconstruction algorithms can potentially improve kinetic parameter precision and accuracy in dynamic PET body imaging. However, construction of a common kinetic model is not always feasible and in contrast to post-reconstruction kinetic analysis, errors in poorly modelled regions may spatially propagate to regions which are well modelled. To reduce error propagation from erroneous model fits, we implement and evaluate a new approach to direct parameter estimation by incorporating a recently proposed kinetic modelling strategy within a direct 4D image reconstruction framework. The algorithm uses a secondary more general model to allow a less constrained model fit in regions where the kinetic model does not accurately describe the underlying kinetics. A portion of the residuals then is adaptively included back into the image whilst preserving the primary model characteristics in other well modelled regions using a penalty term that trades off the models. Using fully 4D simulations based on dynamic [(15)O]H2O datasets, we demonstrate reduction in propagation-related bias for all kinetic parameters. Under noisy conditions, reductions in bias due to propagation are obtained at the cost of increased noise, which in turn results in increased bias and variance of the kinetic parameters. This trade-off reflects the challenge of separating the residuals arising from poor kinetic modelling fits from the residuals arising purely from noise. Nonetheless, the overall root mean square error is reduced in most regions and parameters. Using the adaptive 4D image reconstruction improved model fits can be obtained in poorly modelled regions, leading to reduced errors potentially propagating to regions of interest which the primary biologic model accurately describes. The proposed methodology, however, depends on the secondary model and choosing an optimal model on the residual space is critical in improving model fits.

Pancreatic metabolism, blood flow, and β-cell function in obese humans

CONTEXT

Glucolipotoxicity is believed to induce pancreatic β-cell dysfunction in obesity. Previously, it has not been possible to study pancreatic metabolism and blood flow in humans.

OBJECTIVE

The objective of the study was to investigate whether pancreatic metabolism and blood flow are altered in obesity using positron emission tomography (PET). In the preclinical part, the method was validated in animals.

DESIGN

This was a cross-sectional study.

SETTING

The study was conducted in a clinical research center.

PARTICIPANTS

Human studies consisted of 52 morbidly obese and 25 healthy age-matched control subjects. Validation experiments were done with rodents and pigs.

INTERVENTIONS

PET and magnetic resonance imaging studies using a glucose analog ([(18)F]fluoro-2-deoxy-d-glucose), a palmitate analog [14(R,S)-[(18)F]fluoro-6-thia-heptadecanoic acid], and radiowater ([(15)O]H2O) were performed. In animals, a comparison between ex vivo and in vivo data was performed.

MAIN OUTCOME MEASURES

Pancreatic glucose/fatty acid (FA) uptake, fat accumulation, and blood flow parameters of β-cell function were measured.

RESULTS

PET proved to be a feasible method to measure pancreatic metabolism. Compared with healthy participants, obese participants had elevated pancreatic FA uptake (P < .0001), more fat accumulation (P = .0001), lowered glucose uptake both during fasting and euglycemic hyperinsulinemia, and blunted blood flow (P < .01) in the pancreas. Blood flow, FA uptake, and fat accumulation were negatively associated with multiple markers of β-cell function.

CONCLUSIONS

Obesity leads to changes in pancreatic energy metabolism with a substrate shift from glucose to FAs. In morbidly obese humans, impaired pancreatic blood flow may contribute to β-cell dysfunction and in the pathogenesis of type 2 diabetes.

Rapid quantitative CBF and CMRO(2) measurements from a single PET scan with sequential administration of dual (15)O-labeled tracers

Positron emission tomography (PET) with (15)O tracers provides essential information in patients with cerebral vascular disorders, such as cerebral blood flow (CBF), oxygen extraction fraction (OEF), and metabolic rate of oxygen (CMRO(2)). However, most of techniques require an additional C(15)O scan for compensating cerebral blood volume (CBV). We aimed to establish a technique to calculate all functional images only from a single dynamic PET scan, without losing accuracy or statistical certainties. The technique was an extension of previous dual-tracer autoradiography (DARG) approach, but based on the basis function method (DBFM), thus estimating all functional parametric images from a single session of dynamic scan acquired during the sequential administration of H(2)(15)O and (15)O(2). Validity was tested on six monkeys by comparing global OEF by PET with those by arteriovenous blood sampling, and tested feasibility on young healthy subjects. The mean DBFM-derived global OEF was 0.57±0.06 in monkeys, in an agreement with that by the arteriovenous method (0.54±0.06). Image quality was similar and no significant differences were seen from DARG; 3.57%±6.44% and 3.84%±3.42% for CBF, and -2.79%±11.2% and -6.68%±10.5% for CMRO(2). A simulation study demonstrated similar error propagation between DBFM and DARG. The DBFM method enables accurate assessment of CBF and CMRO(2) without additional CBV scan within significantly shortened examination period, in clinical settings.

Bringing physiology into PET of the liver

Several physiologic features make interpretation of PET studies of liver physiology an exciting challenge. As with other organs, hepatic tracer kinetics using PET is quantified by dynamic recording of the liver after the administration of a radioactive tracer, with measurements of time-activity curves in the blood supply. However, the liver receives blood from both the portal vein and the hepatic artery, with the peak of the portal vein time-activity curve being delayed and dispersed compared with that of the hepatic artery. The use of a flow-weighted dual-input time-activity curve is of importance for the estimation of hepatic blood perfusion through initial dynamic PET recording. The portal vein is inaccessible in humans, and methods of estimating the dual-input time-activity curve without portal vein measurements are being developed. Such methods are used to estimate regional hepatic blood perfusion, for example, by means of the initial part of a dynamic (18)F-FDG PET/CT recording. Later, steady-state hepatic metabolism can be assessed using only the arterial input, provided that neither the tracer nor its metabolites are irreversibly trapped in the prehepatic splanchnic area within the acquisition period. This is used in studies of regulation of hepatic metabolism of, for example, (18)F-FDG and (11)C-palmitate.

Hepatic blood perfusion measured by 3-minute dynamic 18F-FDG PET in pigs

There is an unmet clinical need for an imaging method for quantification of hepatic blood perfusion. The purpose of the present study was to develop and validate a PET method using blood-to-cell clearance (K(1)) of (18)F-FDG, 3-O-(11)C-methylglucose ((11)C-MG), or 2-(18)F-fluoro-2-deoxy-D-galactose ((18)F-FDGal) as a measure of hepatic blood perfusion without the need for portal venous blood samples. We aimed to make the method as simple as possible with the prospect of future application to clinical studies. For this purpose, we examined the possibility of using a 3-min data acquisition and a model-derived dual input calculated from measurements of radioactivity concentrations in a peripheral artery.

METHODS

Pigs (40 kg) underwent dynamic PET of the liver with (18)F-FDG, (11)C-MG, or (18)F-FDGal with simultaneous measurements of time-activity curves in blood sampled from a femoral artery and the portal vein (PV); blood flow rates were measured in the hepatic artery (HA) and PV by transit-time flow meters. Two input functions were compared: A measured dual input and a model-derived dual input, the latter with the PV time-activity curve estimated from the measured arterial time-activity curve and a previously validated 1-parametric PV model. (K(1)) was estimated for each tracer by fitting compartmental models to the data, comparing 3-min and 60-min data acquisitions and the 2 dual-input time-activity curves.

RESULTS

Agreement between (K(1)) estimated using the measured and the model-derived dual input was good for all 3 tracers. For (18)F-FDG and (11)C-MG, (K(1)) (3-min data acquisition, model-derived dual input, and 1-tissue compartmental model) correlated to the measured blood perfusion (P = 0.01 and P = 0.07, respectively). For (18)F-FDGal, the correlation was not significant.

CONCLUSION

A simplified method for quantification of hepatic blood perfusion using 3-min dynamic (18)F-FDG PET or (11)C-MG PET with blood sampling from only a peripheral artery was developed. Parametric (K(1)) images were constructed and showed homogeneous blood perfusion in these normal livers.

Non-invasive estimation of hepatic glucose uptake from [18F]FDG PET images using tissue-derived input functions

PURPOSE

The liver is perfused through the portal vein and hepatic artery. Quantification of hepatic glucose uptake (HGU) using PET requires the use of an input function for both the hepatic artery and portal vein. The former can be generally obtained invasively, but blood withdrawal from the portal vein is not practical in humans. The aim of this study was to develop and validate a new technique to obtain quantitative HGU by estimating the input function from PET images.

METHODS

Normal pigs (n = 12) were studied with [18F]FDG PET, in which arterial and portal blood time-activity curves (TAC) were determined invasively to serve as reference measurements. The present technique consisted of two characteristics, i.e. using a model input function and simultaneously fitting multiple liver tissue TACs from images by minimizing the residual sum of square between the tissue TACs and fitted curves. The input function was obtained from the parameters determined from the fitting. The HGU values were computed by the estimated and measured input functions and compared between the methods.

RESULTS

The estimated input functions were well reproduced. The HGU values, ranging from 0.005 to 0.02 ml/min per ml, were not significantly different between the two methods (r = 0.95, p < 0.001). A Bland-Altman plot demonstrated a small overestimation by the image-derived method with a bias of 0.00052 ml/min per g for HGU.

CONCLUSION

The results presented demonstrate that the input function can be estimated directly from the PET image, supporting the fully non-invasive assessment of liver glucose metabolism in human studies.

Liver fat content in type 2 diabetes: relationship with hepatic perfusion and substrate metabolism

OBJECTIVE

Hepatic steatosis is common in type 2 diabetes. It is causally linked to the features of the metabolic syndrome, liver cirrhosis, and cardiovascular disease. Experimental data have indicated that increased liver fat may impair hepatic perfusion and metabolism. The aim of the current study was to assess hepatic parenchymal perfusion, together with glucose and fatty acid metabolism, in relation to hepatic triglyceride content.

RESEARCH DESIGN AND METHODS

Fifty-nine men with well controlled type 2 diabetes and 18 age-matched healthy normoglycemic men were studied using positron emission tomography to assess hepatic tissue perfusion, insulin-stimulated glucose, and fasting fatty acid metabolism, respectively, in relation to hepatic triglyceride content, quantified by proton magnetic resonance spectroscopy. Patients were divided into two groups with hepatic triglyceride content below (type 2 diabetes-low) or above (type 2 diabetes-high) the median of 8.6%.

RESULTS

Type 2 diabetes-high patients had the highest BMI and A1C and lowest whole-body insulin sensitivity (ANOVA, all P < 0.001). Compared with control subjects and type 2 diabetes-low patients, type 2 diabetes-high patients had the lowest hepatic parenchymal perfusion (P = 0.004) and insulin-stimulated hepatic glucose uptake (P = 0.013). The observed decrease in hepatic fatty acid influx rate constant, however, only reached borderline significance (P = 0.088). In type 2 diabetic patients, hepatic parenchymal perfusion (r = -0.360, P = 0.007) and hepatic fatty acid influx rate constant (r = -0.407, P = 0.007) correlated inversely with hepatic triglyceride content. In a pooled analysis, hepatic fat correlated with hepatic glucose uptake (r = -0.329, P = 0.004).

CONCLUSIONS

In conclusion, type 2 diabetic patients with increased hepatic triglyceride content showed decreased hepatic parenchymal perfusion and hepatic insulin mediated glucose uptake, suggesting a potential modulating effect of hepatic fat on hepatic physiology.

Assessment of tissue perfusion by contrast-enhanced ultrasound

Contrast-enhanced ultrasound (CEUS) with microbubble contrast agents is a new imaging technique for quantifying tissue perfusion. CEUS presents several advantages over other imaging techniques in assessing tissue perfusion, including the use of microbubbles as blood-pool agents, portability, availability and absence of exposure to radiation or nuclear tracers. Dedicated software packages are necessary to quantify the echo-signal intensity and allow the calculation of the degree of tissue contrast enhancement based on the accurate distinction between microbubble backscatter signals and native tissue background. The measurement of organ transit time after microbubble injection and the analysis of tissue reperfusion kinetics represent the two fundamental methods for the assessment of tissue perfusion by CEUS. Transit time measurement has been shown to be feasible and has started to become accepted as a clinical tool, especially in the liver. The loudness of audio signals from spectral Doppler analysis is used to generate time-intensity curves to follow the wash-in and wash-out of the microbubble bolus. Tissue perfusion may be quantified also by analysing the replenishment kinetics of the volume of microbubbles after their destruction in the imaged slice. This allows to obtain semiquantitative parameters related to local tissue perfusion, especially in the heart, brain, and kidneys.

Tracer input for kinetic modelling of liver physiology determined without sampling portal venous blood in pigs

Purpose

Quantification of hepatic tracer kinetics by PET requires measurement of tracer input from the hepatic artery (HA) and portal vein (PV). We wished to develop a method for estimating dual tracer input without the necessity to sample PV blood.

Methods

Pigs weighing 40 kg were given bolus doses of C¹⁵O (CO), 2-[¹⁸F]fluoro-2-deoxy-D-glucose (FDG), [¹¹C]-methylglucose (MG), 2-[¹⁸F]fluoro-2-deoxy-D-galactose (FDGal) or H₂¹⁵O (H₂O). Tracer concentration 3-min time courses were measured in the femoral artery and PV by blood sampling. Blood flow was measured in the HA and PV using flow-meters. A model for transfer of tracer through the splanchnic circulation was used to estimate values of a tracer-specific model parameter β. Tracer-specific mean values of β were used to estimate tracer concentration time courses in the PV from the measured arterial concentration. A model-derived dual-input was calculated using the mean HA flow fraction (0.25) and validated by comparison of the use of the measured dual-input and a kinetic model with a fixed ”true” K₁^true, i.e. clearance of tracer from blood to liver cells.

Results

The rank order of the means of β was CO < FDG ≈ MG < FDGal < H₂O, reflecting their different splanchnic mean transit times. Estimated K₁^est was not significantly different from “true” K₁^true.

Conclusion

The hepatic dual tracer input, which is of great importance for the assessment of processes such as transfer across the plasma-hepatocyte membrane or hepatic blood perfusion, can be well approximated in pigs without the necessity to sample PV blood and measure hepatic blood flow; only arterial blood sampling is needed.