Tracer Kinetic Modeling Research Articles

Comparison of tracer kinetic models in differentiating malignant from normal prostate tissue using dynamic contrast-enhanced MRI

PurposeThe aim of this study was to evaluate the diagnostic value of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) derived kinetic parameters with high spatiotemporal resolution in discriminating malignant from normal prostate tissues.MethodsFifty patients with suspicious of malignant diseases in prostate were included in this study. Regions of interest (ROI) were manually delineated by experienced radiologists. Voxel-wise kinetic parameters were produced with the following tracer kinetic models (TKMs): Tofts model, extended Tofts model (ETM), Brix’s conventional two-compartment model (Brix), adiabatic tissue homogeneity model (ATH), and distributed parameter model (DP). The initial area under the signal-time curve (IAUC) with an uptake integral approach was also included. Mann–Whitney U test and receiver operating characteristic (ROC) curves were used to evaluate the capability of distinguishing tumor lesions from normal tissues. A p-value of 0.05 or less is considered statistically significant. ROI based parameters correlation analysis between DP and ETM were performed.Results624 lesions and 269 normal tissue ROIs were obtained. Thirty parameters were derived from the six kinetic models. Except for PS from Brix, statistically significant differences between lesions and normal tissues (P<0.05) were observed in other parameters.Ve from DP, ATH and Brix and PS from ATH have AUC values less than 0.6 in the ROC analysis. MTT, Vp and PS from DP, Ktrans from ETM and Tofts, E and PS from ATH, IAUC parameters and F from Brix have AUC values larger than 0.8. Ve and Vp from DP and ETM are correlated (r> 0.65). The correlation coefficient between Ktrans from ETM and PS from DP is 0.751.ConclusionMTT, Vp and PS from DP, Ktrans from ETM and Tofts, E and PS from ATH, F from Brix and IAUC parameters can be used to differentiate malignant lesions from normal tissues in the prostate.

Retinal Vascular Permeability in Diabetic Subjects without Retinopathy Compared with Mild Diabetic Retinopathy and Healthy Controls

ObjectivesTo investigate retinal vascular permeability mapping as a potential biomarker for diabetic retinopathy in subjects with diabetes with no signs of retinopathy and with mild non-proliferative retinopathy. DesignThis is a case-control study. SubjectsParticipants included seven healthy controls, twenty-two subjects with diabetes mellitus and no clinical signs of retinopathy (DMnoDR), and seven subjects with mild non-proliferative diabetic retinopathy (NPDR). MethodsAll participants underwent routine retinal fluorescein videoangiography (FVA). Each FVA dataset was analyzed with the dynamic tracer kinetic model (DTKM) method to estimate five parameters: extraction fraction, blood flow, arrival time, transit time, and rate constant defined via adiabatic solution. The DTKM method was based on indicator dilution theory, including sequential use of two prominent kinetic models: the plug flow model and the adiabatic approximation to the tissue homogeneity model. Main Outcome MeasuresExtraction fraction (E), i.e., the fluorescein dye leakage measured during one pass through surrounding retinal tissue, is extracted via DTKM method and directly relates to retinal vascular permeability. Thus, E represents the pre-clinical biomarker, retinal vascular permeability. ResultsThe three diagnostic groups were found to have significantly different permeability (p = 0.003). Despite having no clinical signs of retinopathy, the mean rank of average vascular extraction fractions (E) was significantly higher in DMnoDR subjects compared to healthy controls (p = 0.04), as was the mean rank of E for mild NPDR subjects (p = 0.002). The average E for mild NPDR, DMnoDR, and control subjects was 0.1 ± 0.04, 0.07 ± 0.04, and 0.04 ± 0.01, respectively. ConclusionsThe vascular permeability extracted from FVA datasets using the DTKM method is a promising biomarker for detecting preclinical retinal pathology in patients with diabetes. Longitudinal studies are ongoing to explore the ability of this biomarker to distinguish those subjects with diabetes who will progress to clinically apparent retinopathy from those who will not.

Cerebral perfusion metrics calculated directly from a hypoxia-induced step change in deoxyhemoglobin

Resting cerebral perfusion metrics can be calculated from the MRI ΔR2* signal during the first passage of an intravascular bolus of a Gadolinium-based contrast agent (GBCA), or more recently, a transient hypoxia-induced change in the concentration of deoxyhemoglobin ([dOHb]). Conventional analysis follows a proxy process that includes deconvolution of an arterial input function (AIF) in a tracer kinetic model. We hypothesized that the step reduction in magnetic susceptibility accompanying a step decrease in [dOHb] that occurs when a single breath of oxygen terminates a brief episode of lung hypoxia permits direct calculation of relative perfusion metrics. The time course of the ΔR2* signal response enables both the discrimination of blood arrival times and the time course of voxel filling. We calculated the perfusion metrics implied by this step signal change in seven healthy volunteers and compared them to those from conventional analyses of GBCA and dOHb using their AIF and indicator dilution theory. Voxel-wise maps of relative cerebral blood flow and relative cerebral blood volume had a high spatial and magnitude congruence for all three analyses (r > 0.9) and were similar in appearance to published maps. The mean (SD) transit times (s) in grey and white matter respectively for the step response (7.4 (1.1), 8.05 (1.71)) were greater than those for GBCA (2.6 (0.45), 3.54 (0.83)) attributable to the nature of their respective calculation models. In conclusion we believe these calculations of perfusion metrics derived directly from ΔR2* have superior merit to calculations via AIF by virtue of being calculated from a direct signal rather than through a proxy model which encompasses errors inherent in designating an AIF and performing deconvolution calculations.

Multi-Tracer Studies of Brain Oxygen and Glucose Metabolism Using a Time-of-Flight Positron Emission Tomography - Computed Tomography Scanner.

The authors have developed a paradigm using positron emission tomography (PET) with multiple radiopharmaceutical tracers that combines measurements of cerebral metabolic rate of glucose (CMRGlc), cerebral metabolic rate of oxygen (CMRO2), cerebral blood flow (CBF), and cerebral blood volume (CBV), culminating in estimates of brain aerobic glycolysis (AG). These in vivo estimates of oxidative and non-oxidative glucose metabolism are pertinent to the study of the human brain in health and disease. The latest positron emission tomography-computed tomography (PET-CT) scanners provide time-of-flight (TOF) imaging and critical improvements in spatial resolution and reduction of artifacts. This has led to significantly improved imaging with lower radiotracer doses. Optimized methods for the latest PET-CT scanners involve administering a sequence of inhaled 15O-labeled carbon monoxide (CO) and oxygen (O2), intravenous 15O-labeled water (H2O), and 18F-deoxyglucose (FDG)-all within 2-h or 3-h scan sessions that yield high-resolution, quantitative measurements of CMRGlc, CMRO2, CBF, CBV, and AG. This methods paper describes practical aspects of scanning designed for quantifying brain metabolism with tracer kinetic models and arterial blood samples and provides examples of imaging measurements of human brain metabolism.

DIFFUSION MODEL-BASED POSTERIOR DISTRIBUTION PREDICTION FOR KINETIC PARAMETER ESTIMATION IN DYNAMIC PET.

Positron Emission Tomography (PET) is a valuable imaging method for studying molecular-level processes in the body, such as hyperphosphorylated tau (p-tau) protein aggregates, a hallmark of several neurodegenerative diseases including Alzheimer's disease. P-tau density and cerebral perfusion can be quantified from PET data using tracer kinetic modeling techniques. However, noise in PET images leads to uncertainty in the estimated kinetic parameters. This can be quantified in a Bayesian framework by the posterior distribution of kinetic parameters given PET measurements. Markov Chain Monte Carlo (MCMC) techniques can be employed to estimate the posterior distribution, although with significant computational needs. In this paper, we propose to leverage deep learning inference efficiency to infer the posterior distribution. A novel approach using denoising diffusion probabilistic model (DDPM) is introduced. The performance of the proposed method was evaluated on a [18F]MK6240 study and compared to an MCMC method. Our approach offered significant reduction in computation time (over 30 times faster than MCMC) and consistently predicted accurate (< 0.8 % mean error) and precise (< 5.77 % standard deviation error) posterior distributions.

Applying dynamic contrast-enhanced MRI tracer kinetic models to differentiate benign and malignant soft tissue tumors

BackgroundTo explore the potential of different quantitative dynamic contrast-enhanced (qDCE)-MRI tracer kinetic (TK) models and qDCE parameters in discriminating benign from malignant soft tissue tumors (STTs).MethodsThis research included 92 patients (41females, 51 males; age range 16–86 years, mean age 51.24 years) with STTs. The qDCE parameters (Ktrans, Kep, Ve, Vp, F, PS, MTT and E) for regions of interest of STTs were estimated by using the following TK models: Tofts (TOFTS), Extended Tofts (EXTOFTS), adiabatic tissue homogeneity (ATH), conventional compartmental (CC), and distributed parameter (DP). We established a comprehensive model combining the morphologic features, time-signal intensity curve shape, and optimal qDCE parameters. The capacities to identify benign and malignant STTs was evaluated using the area under the curve (AUC), degree of accuracy, and the analysis of the decision curve.ResultsTOFTS-Ktrans, EXTOFTS-Ktrans, EXTOFTS-Vp, CC-Vp and DP-Vp demonstrated good diagnostic performance among the qDCE parameters. Compared with the other TK models, the DP model has a higher AUC and a greater level of accuracy. The comprehensive model (AUC, 0.936, 0.884–0.988) demonstrated superiority in discriminating benign and malignant STTs, outperforming the qDCE models (AUC, 0.899–0.915) and the traditional imaging model (AUC, 0.802, 0.712–0.891) alone.ConclusionsVarious TK models successfully distinguish benign from malignant STTs. The comprehensive model is a noninvasive approach incorporating morphological imaging aspects and qDCE parameters, and shows significant potential for further development.

Optimization-derived blood input function using a kernel method and its evaluation with total-body PET for brain parametric imaging

Dynamic PET allows quantification of physiological parameters through tracer kinetic modeling. For dynamic imaging of brain or head and neck cancer on conventional PET scanners with a short axial field of view, the image-derived input function (ID-IF) from intracranial blood vessels such as the carotid artery (CA) suffers from severe partial volume effects. Alternatively, optimization-derived input function (OD-IF) by the simultaneous estimation (SIME) method does not rely on an ID-IF but derives the input function directly from the data. However, the optimization problem is often highly ill-posed. We proposed a new method that combines the ideas of OD-IF and ID-IF together through a kernel framework. While evaluation of such a method is challenging in human subjects, we used the uEXPLORER total-body PET system that covers major blood pools to provide a reference for validation. MethodsThe conventional SIME approach estimates an input function using a joint estimation together with kinetic parameters by fitting time activity curves from multiple regions of interests (ROIs). The input function is commonly parameterized with a highly nonlinear model which is difficult to estimate. The proposed kernel SIME method exploits the CA ID-IF as a priori information via a kernel representation to stabilize the SIME approach. The unknown parameters are linear and thus easier to estimate. The proposed method was evaluated using 18F-fluorodeoxyglucose studies with both computer simulations and 20 human-subject scans acquired on the uEXPLORER scanner. The effect of the number of ROIs on kernel SIME was also explored. ResultsThe estimated OD-IF by kernel SIME showed a good match with the reference input function and provided more accurate estimation of kinetic parameters for both simulation and human-subject data. The kernel SIME led to the highest correlation coefficient (R = 0.97) and the lowest mean absolute error (MAE = 10.5 %) compared to using the CA ID-IF (R = 0.86, MAE = 108.2 %) and conventional SIME (R = 0.57, MAE = 78.7 %) in the human-subject evaluation. Adding more ROIs improved the overall performance of the kernel SIME method. ConclusionThe proposed kernel SIME method shows promise to provide an accurate estimation of the blood input function and kinetic parameters for brain PET parametric imaging.

Pulmonary CT perfusion robustly measures cardiac output in the context of multilevel pulmonary occlusion: a porcine study

BackgroundTo validate pulmonary computed tomography (CT) perfusion in a porcine model by invasive monitoring of cardiac output (CO) using thermodilution method.MethodsAnimals were studied at a single center, using a Swan-Ganz catheter for invasive CO monitoring as a reference. Fifteen pigs were included. Contrast-enhanced CT perfusion of the descending aorta and right and left pulmonary artery was performed. For variation purposes, a balloon catheter was inserted to block the contralateral pulmonary vascular bed; additionally, two increased CO settings were created by intravenous administration of catecholamines. Finally, stepwise capillary occlusion was performed by intrapulmonary arterial injection of 75-μm microspheres in four stages. A semiautomatic selection of AFs and a recirculation-aware tracer-kinetics model to extract the first-pass of AFs, estimating blood flow with the Stewart-Hamilton method, was implemented. Linear mixed models (LMM) were developed to calibrate blood flow calculations accounting with individual- and cohort-level effects.ResultsNine of 15 pigs had complete datasets. Strong correlations were observed between calibrated pulmonary (0.73, 95% confidence interval [CI] 0.6–0.82) and aortic blood flow measurements (0.82, 95% CI, 0.73–0.88) and the reference as well as agreements (± 2.24 L/min and ± 1.86 L/min, respectively) comparable to the state of the art, on a relatively wide range of right ventricle-CO measurements.ConclusionsCT perfusion validly measures CO using LMMs at both individual and cohort levels, as demonstrated by referencing the invasive CO.Relevance statementPossible clinical applications of CT perfusion for measuring CO could be in acute pulmonary thromboembolism or to assess right ventricular function to show impairment or mismatch to the left ventricle.Key points• CT perfusion measures flow in vessels.• CT perfusion measures cumulative cardiac output in the aorta and pulmonary vessels.• CT perfusion validly measures CO using LMMs at both individual and cohort levels, as demonstrated by using the invasive CO as a reference standard.Graphical

Unpaired deep learning for pharmacokinetic parameter estimation from dynamic contrast-enhanced MRI without AIF measurements

DCE-MRI provides information about vascular permeability and tissue perfusion through the acquisition of pharmacokinetic parameters. However, traditional methods for estimating these pharmacokinetic parameters involve fitting tracer kinetic models, which often suffer from computational complexity and low accuracy due to noisy arterial input function (AIF) measurements. Although some deep learning approaches have been proposed to tackle these challenges, most existing methods rely on supervised learning that requires paired input DCE-MRI and labeled pharmacokinetic parameter maps. This dependency on labeled data introduces significant time and resource constraints and potential noise in the labels, making supervised learning methods often impractical. To address these limitations, we present a novel unpaired deep learning method for estimating pharmacokinetic parameters and the AIF using a physics-driven CycleGAN approach. Our proposed CycleGAN framework is designed based on the underlying physics model, resulting in a simpler architecture with a single generator and discriminator pair. Crucially, our experimental results indicate that our method does not necessitate separate AIF measurements and produces more reliable pharmacokinetic parameters than other techniques.

PET-determined myocardial perfusion and flow in coronary artery disease characterization

Positron emission tomography (PET) myocardial perfusion imaging in conjunction with tracer-kinetic modeling enables the concurrent assessment of myocardial perfusion and regional myocardial blood flow (MBF) of the left ventricle in absolute terms in milliliters per gram per minute (mL/g/min). The non-invasive quantification of MBF during pharmacologically induced hyperemia, at rest, and corresponding myocardial flow reserve (MFR) opens a new avenue for the identification and characterization of classical or endogen type of coronary microvascular dysfunction (CMD) as functional substrate for microvascular angina in patients with non-obstructive coronary artery disease (CAD) and/or no CAD at all. Further, PET-MBF quantification expands the scope of conventional myocardial perfusion imaging from the identification of advanced, and flow-limiting, epicardial CAD to early stages of atherosclerosis and/or CMD. Adding MBF assessment to myocardial perfusion may also reliably unravel diffuse ischemia owing to significant left main stenosis and/or multivessel CAD, commonly confirmed by peak stress transient ischemic cavity dilation of the left ventricle during maximal vasomotor stress compared to rest on gated PET images. Owing to high spatial and contrast resolution in conjunction with photon-attenuation free myocardial perfusion PET images, PET is preferentially used for CAD detection in advanced obesity and women with pronounced breast habitus. With increasing clinical use of cardiac PET perfusion and MBF assessment, individualized, and image-guided cardiovascular treatment decisions in CAD patients is likely to ensue, while its translation into improved cardiovascular outcome remains to be investigated.

Total-Body PET/CT Applications in Cardiovascular Diseases: A Perspective Document of the SNMMI Cardiovascular Council.

Digital PET/CT systems with a long axial field of view have become available and are emerging as the current state of the art. These new camera systems provide wider anatomic coverage, leading to major increases in system sensitivity. Preliminary results have demonstrated improvements in image quality and quantification, as well as substantial advantages in tracer kinetic modeling from dynamic imaging. These systems also potentially allow for low-dose examinations and major reductions in acquisition time. Thereby, they hold great promise to improve PET-based interrogation of cardiac physiology and biology. Additionally, the whole-body coverage enables simultaneous assessment of multiple organs and the large vascular structures of the body, opening new opportunities for imaging systemic mechanisms, disorders, or treatments and their interactions with the cardiovascular system as a whole. The aim of this perspective document is to debate the potential applications, challenges, opportunities, and remaining challenges of applying PET/CT with a long axial field of view to the field of cardiovascular disease.

Motion-compensated image reconstruction for improved kidney function assessment using dynamic contrast-enhanced MRI.

Accurately measuring renal function is crucial for pediatric patients with kidney conditions. Traditional methods have limitations, but dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides a safe and efficient approach for detailed anatomical evaluation and renal function assessment. However, motion artifacts during DCE-MRI can degrade image quality and introduce misalignments, leading to unreliable results. This study introduces a motion-compensated reconstruction technique for DCE-MRI data acquired using golden-angle radial sampling. Our proposed method achieves three key objectives: (1) identifying and removing corrupted data (outliers) using a Gaussian process model fitting with a -space center navigator, (2) efficiently clustering the data into motion phases and performing interphase registration, and (3) utilizing a novel formulation of motion-compensated radial reconstruction. We applied the proposed motion correction (MoCo) method to DCE-MRI data affected by varying degrees of motion, including both respiratory and bulk motion. We compared the outcomes with those obtained from the conventional radial reconstruction. Our evaluation encompassed assessing the quality of images, concentration curves, and tracer kinetic model fitting, and estimating renal function. The proposed MoCo reconstruction improved the temporal signal-to-noise ratio for all subjects, with a 21.8% increase on average, while total variation values of the aorta, right, and left kidney concentration were improved for each subject, with 32.5%, 41.3%, and 42.9% increases on average, respectively. Furthermore, evaluation of tracer kinetic model fitting indicated that the median standard deviation of the estimated filtration rate ( ), mean normalized root-mean-squared error (nRMSE), and chi-square goodness-of-fit of tracer kinetic model fit were decreased from 0.10 to 0.04, 0.27 to 0.24, and, 0.43 to 0.27, respectively. The proposed MoCo technique enabled more reliable renal function assessment and improved image quality for detailed anatomical evaluation in the case of bulk and respiratory motion during the acquisition of DCE-MRI.

Contrast-free dynamic susceptibility contrast using sinusoidal and bolus oxygenation challenges.

Deoxygenation-based dynamic susceptibility contrast (dDSC) MRI uses respiratory challenges as a source of endogenous contrast as an alternative to gadolinium injection. These gas challenges induce T2*-weighted MRI signal losses, after which tracer kinetics modeling was applied to calculate cerebral perfusion. This work compares three gas challenges, desaturation (transient hypoxia), resaturation (transient normoxia), and SineO2 (sinusoidal modulation of end-tidal oxygen pressures) in a cohort of 10 healthy volunteers (age 37 ± 11 years; 60% female). Perfusion estimates consisted of cerebral blood flow (CBF), cerebral blood volume (CBV), and mean transit time (MTT). Calculations were computed using a traditional tracer kinetics model in the time domain for desaturation and resaturation and in the frequency domain for SineO2. High correlations and limits of agreement were observed among the three deoxygenation-based paradigms for CBV, although MTT and CBF estimates varied with the hypoxic stimulus. Cross-modality correlation with gadolinium DSC was lower, particularly for MTT, but on a par with agreement between the other perfusion references. Overall, this work demonstrated the feasibility and reliability of oxygen respiratory challenges to measure brain perfusion. Additional work is needed to assess the utility of dDSC in the diagnostic evaluation of various pathologies such as ischemic strokes, brain tumors, and neurodegenerative diseases.

Comparison of tracer kinetic models for 68Ga-PSMA-11 PET in intermediate-risk primary prostate cancer patients

Background68Ga-PSMA-11 positron emission tomography enables the detection of primary, recurrent, and metastatic prostate cancer. Regional radiopharmaceutical uptake is generally evaluated in static images and quantified as standard uptake values (SUVs) for clinical decision-making. However, analysis of dynamic images characterizing both tracer uptake and pharmacokinetics may offer added insights into the underlying tissue pathophysiology. This study was undertaken to evaluate the suitability of various kinetic models for 68Ga-PSMA-11 PET analysis. Twenty-three lesions in 18 patients were included in a retrospective kinetic evaluation of 55-min dynamic 68Ga-PSMA-11 pre-prostatectomy PET scans from patients with biopsy-demonstrated intermediate- to high-risk prostate cancer. Three kinetic models—a reversible one-tissue compartment model, an irreversible two-tissue compartment model, and a reversible two-tissue compartment model, were evaluated for their goodness of fit to lesion and normal reference prostate time-activity curves. Kinetic parameters obtained through graphical analysis and tracer kinetic modeling techniques were compared for reference prostate tissue and lesion regions of interest.ResultsSupported by goodness of fit and information loss criteria, the irreversible two-tissue compartment model optimally fit the time-activity curves. Lesions exhibited significant differences in kinetic rate constants (K1, k2, k3, Ki) and semiquantitative measures (SUV and %ID/kg) when compared with reference prostatic tissue. The two-tissue irreversible tracer kinetic model was consistently appropriate across prostatic zones.ConclusionsAn irreversible tracer kinetic model is appropriate for dynamic analysis of 68Ga-PSMA-11 PET images. Kinetic parameters estimated by Patlak graphical analysis or full compartmental analysis can distinguish tumor from normal prostate tissue.

Review of tracer kinetic models in evaluation of gliomas using dynamic contrast-enhanced imaging.

Glioma is the most common type of primary malignant tumor of the central nervous system (CNS), and is characterized by high malignancy, high recurrence rate and poor survival. Conventional imaging techniques only provide information regarding the anatomical location, morphological characteristics, and enhancement patterns. In contrast, advanced imaging techniques such as dynamic contrast-enhanced (DCE) MRI or DCE CT can reflect tissue microcirculation, including tumor vascular hyperplasia and vessel permeability. Although several studies have used DCE imaging to evaluate gliomas, the results of data analysis using conventional tracer kinetic models (TKMs) such as Tofts or extended-Tofts model (ETM) have been ambiguous. More advanced models such as Brix's conventional two-compartment model (Brix), tissue homogeneity model (TH) and distributed parameter (DP) model have been developed, but their application in clinical trials has been limited. This review attempts to appraise issues on glioma studies using conventional TKMs, such as Tofts or ETM model, highlight advancement of DCE imaging techniques and provides insights on the clinical value of glioma management using more advanced TKMs.

Assessing Perfusion in Steno-Occlusive Cerebrovascular Disease Using Transient Hypoxia-Induced Deoxyhemoglobin as a Dynamic Susceptibility Contrast Agent.

Resting brain tissue perfusion in cerebral steno-occlusive vascular disease can be assessed by MR imaging using gadolinium-based susceptibility contrast agents. Recently, transient hypoxia-induced deoxyhemoglobin has been investigated as a noninvasive MR imaging contrast agent. Here we present a comparison of resting perfusion metrics using transient hypoxia-induced deoxyhemoglobin and gadolinium-based contrast agents in patients with known cerebrovascular steno-occlusive disease. Twelve patients with steno-occlusive disease underwent DSC MR imaging using a standard bolus of gadolinium-based contrast agent compared with transient hypoxia-induced deoxyhemoglobin generated in the lungs using an automated gas blender. A conventional multi-slice 2D gradient echo sequence was used to acquire the perfusion data and analyzed using a standard tracer kinetic model. MTT, relative CBF, and relative CBV maps were generated and compared between contrast agents. The spatial distributions of the perfusion metrics generated with both contrast agents were consistent. Perfusion metrics in GM and WM were not statistically different except for WM MTT. Cerebral perfusion metrics generated with noninvasive transient hypoxia-induced changes in deoxyhemoglobin are very similar to those generated using a gadolinium-based contrast agent in patients with cerebrovascular steno-occlusive disease.

Exploring the potential of Physics-Informed Neural Networks to extract vascularization data from DCE-MRI in the presence of diffusion

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is widely used to assess tissue vascularization, particularly in oncological applications. However, the most widely used pharmacokinetic (PK) models do not account for contrast agent (CA) diffusion between neighboring voxels, which can limit the accuracy of the results, especially in cases of heterogeneous tumors. To address this issue, previous works have proposed algorithms that incorporate diffusion phenomena into the formulation. However, these algorithms often face convergence problems due to the ill-posed nature of the problem.In this work, we present a new approach to fitting DCE-MRI data that incorporates CA diffusion by using Physics-Informed Neural Networks (PINNs). PINNs can be trained to fit measured data obtained from DCE-MRI while ensuring the mass conservation equation from the PK model. We compare the performance of PINNs to previous algorithms on different 1D cases inspired by previous works from literature. Results show that PINNs retrieve vascularization parameters more accurately from diffusion-corrected tracer-kinetic models. Furthermore, we demonstrate the robustness of PINNs compared to other traditional algorithms when faced with noisy or incomplete data. Overall, our results suggest that PINNs can be a valuable tool for improving the accuracy of DCE-MRI data analysis, particularly in cases where CA diffusion plays a significant role.

Quantitative functional imaging with CT perfusion: technical considerations, kinetic modeling, and applications

CT perfusion (CTP)-derived quantitative maps of hemodynamic parameters have found important clinical applications in stroke, cancer, and cardiovascular disease. Blood flow, blood volume, transit time, and other perfusion parameters are sensitive markers of pathophysiology with impaired perfusion. This review summarizes the basic principles of CTP including image acquisition, tracer kinetic modeling, deconvolution algorithms, and diagnostic interpretation. The focus is on practical and theoretical considerations for accurate quantitative parametric imaging. Recommended CTP scan parameters to maintain CT number accuracy and optimize radiation dose versus image noise are first reviewed. Tracer kinetic models, which describe how injected contrast material is distributed between blood and the tissue microenvironment by perfusion and bidirectional passive exchange, are then derived. Deconvolution algorithms to solve for hemodynamic parameters of kinetic models are discussed and their quantitative accuracy benchmarked. The applications and diagnostic interpretation of CTP in stroke, cancer, and cardiovascular disease are summarized. Finally, we conclude with a discussion of future directions for CTP research, including radiation dose reduction, new opportunities with novel CT hardware, and emerging diagnostic applications.

Efficient estimation of pharmacokinetic parameters from breast dynamic contrast-enhanced MRI based on a convolutional neural network for predicting molecular subtypes

Objective. Tracer kinetic models allow for estimating pharmacokinetic (PK) parameters, which are related to pathological characteristics, from breast dynamic contrast-enhanced magnetic resonance imaging. However, existing tracer kinetic models subject to inaccuracy are time-consuming for PK parameters estimation. This study aimed to accurately and efficiently estimate PK parameters for predicting molecular subtypes based on convolutional neural network (CNN). Approach. A CNN integrating global and local features (GL-CNN) was trained using synthetic data where known PK parameters map was used as the ground truth, and subsequently used to directly estimate PK parameters (volume transfer constant K trans and flux rate constant K ep) map. The accuracy assessed by the peak signal-to-noise ratio (PSNR), structural similarity (SSIM), and concordance correlation coefficient (CCC) was compared between the GL-CNN and Tofts-based PK parameters in synthetic data. Radiomic features were calculated from the PK parameters map in 208 breast tumors. A random forest classifier was constructed to predict molecular subtypes using a discovery cohort (n = 144). The diagnostic performance evaluated on a validation cohort (n = 64) using the area under the receiver operating characteristic curve (AUC) was compared between the GL-CNN and Tofts-based PK parameters. Main results. The average PSNR (48.8884), SSIM (0.9995), and CCC (0.9995) between the GL-CNN-based K trans map and ground truth were significantly higher than those between the Tofts-based K trans map and ground truth. The GL-CNN-based K trans obtained significantly better diagnostic performance (AUCs = 0.7658 and 0.8528) than the Tofts-based K trans for luminal B and HER2 tumors. The GL-CNN method accelerated the computation by speed approximately 79 times compared to the Tofts method for the whole breast of all patients. Significance. Our results indicate that the GL-CNN method can be used to accurately and efficiently estimate PK parameters for predicting molecular subtypes.

Total-Body Multiparametric PET Quantification of 18F-FDG Delivery and Metabolism in the Study of Coronavirus Disease 2019 Recovery.

Conventional whole-body static 18F-FDG PET imaging provides a semiquantitative evaluation of overall glucose metabolism without insight into the specific transport and metabolic steps. Here we demonstrate the ability of total-body multiparametric 18F-FDG PET to quantitatively evaluate glucose metabolism using macroparametric quantification and assess specific glucose delivery and phosphorylation processes using microparametric quantification for studying recovery from coronavirus disease 2019 (COVID-19). Methods: The study included 13 healthy subjects and 12 recovering COVID-19 subjects within 8 wk of confirmed diagnosis. Each subject had a 1-h dynamic 18F-FDG scan on the uEXPLORER total-body PET/CT system. Semiquantitative SUV and the SUV ratio relative to blood (SUVR) were calculated for different organs to measure glucose utilization. Tracer kinetic modeling was performed to quantify the microparametric blood-to-tissue 18F-FDG delivery rate [Formula: see text] and the phosphorylation rate k 3, as well as the macroparametric 18F-FDG net influx rate ([Formula: see text]). Statistical tests were performed to examine differences between healthy subjects and recovering COVID-19 subjects. The effect of COVID-19 vaccination was also investigated. Results: We detected no significant difference in lung SUV but significantly higher lung SUVR and [Formula: see text] in COVID-19 recovery, indicating improved sensitivity of kinetic quantification for detecting the difference in glucose metabolism. A significant difference was also observed in the lungs with the phosphorylation rate k 3 but not with [Formula: see text], which suggests that glucose phosphorylation, rather than glucose delivery, drives the observed difference of glucose metabolism. Meanwhile, there was no or little difference in bone marrow 18F-FDG metabolism measured with SUV, SUVR, and [Formula: see text] but a significantly higher bone marrow [Formula: see text] in the COVID-19 group, suggesting a difference in glucose delivery. Vaccinated COVID-19 subjects had a lower lung [Formula: see text] and a higher spleen [Formula: see text] than unvaccinated COVID-19 subjects. Conclusion: Higher lung glucose metabolism and bone marrow glucose delivery were observed with total-body multiparametric 18F-FDG PET in recovering COVID-19 subjects than in healthy subjects, implying continued inflammation during recovery. Vaccination demonstrated potential protection effects. Total-body multiparametric PET of 18F-FDG can provide a more sensitive tool and more insights than conventional whole-body static 18F-FDG imaging to evaluate metabolic changes in systemic diseases such as COVID-19.