Hildebrand Dijkstra

Quantitative and qualitative assessment of structural magnetic resonance imaging data in a two-center study

BackgroundMulti-center magnetic resonance imaging (MRI) studies present an opportunity to advance research by pooling data. However, brain measurements derived from MR-images are susceptible to differences in MR-sequence parameters. It is therefore necessary to determine whether there is an interaction between the sequence parameters and the effect of interest, and to minimise any such interaction by careful choice of acquisition parameters. As an exemplar of the issues involved in multi-center studies, we present data from a study in which we aimed to optimize a set of volumetric MRI-protocols to define a protocol giving data that are consistent and reproducible across two centers and over time.MethodsOptimization was achieved based on data quality and quantitative measures, in our case using FreeSurfer and Voxel Based Morphometry approaches. Our approach consisted of a series of five comparisons. Firstly, a single-center dataset was collected, using a range of candidate pulse-sequences and parameters chosen on the basis of previous literature. Based on initial results, a number of minor changes were implemented to optimize the pulse-sequences, and a second single-center dataset was collected. FreeSurfer data quality measures were compared between datasets in order to determine the best performing sequence(s), which were taken forward to the next stage of testing. We subsequently acquired short-term and long-term two-center reproducibility data, and quantitative measures were again assessed to determine the protocol with the highest reproducibility across centers. Effects of a scanner software and hardware upgrade on the reproducibility of the protocols at one of the centers were also evaluated.ResultsAssessing the quality measures from the first two datasets allowed us to define artefact-free protocols, all with high image quality as assessed by FreeSurfer. Comparing the quantitative test and retest measures, we found high within-center reproducibility for all protocols, but lower between-center reproducibility for some protocols than others. The upgrade showed no important effects.ConclusionsWe were able to determine (for the scanners used in this study) an optimised protocol, which gave the highest within- and between-center reproducibility of those assessed, and give details of this protocol here. More generally, we discuss some of the issues raised by multi-center studies and describe a methodical approach to take towards optimization and standardization, and recommend performing this kind of procedure to other investigators.

64 slice MDCT generally underestimates coronary calcium scores as compared to EBT: A phantom study

The objective of our study was the determination of the influence of the sequential and spiral acquisition modes on the concordance and deviation of the calcium score on 64-slice multi-detector computed tomography (MDCT) scanners in comparison to electron beam tomography (EBT) as the gold standard. Our methods and materials were an anthropomorphic cardio CT phantom with different calcium inserts scanned in sequential and spiral acquisition modes on three identical 64-slice MDCT scanners of manufacturer A and on three identical 64-slice MDCT scanners of manufacturer B and on an EBT system. Every scan was repeated 30 times with and 15 times without a small random variation in the phantom position for both sequential and spiral modes. Significant differences were observed between EBT and 64-slice MDCT data for all inserts, both acquisition modes, and both manufacturers of MDCT systems. High regression coefficients (0.90-0.98) were found between the EBT and 64-slice MDCT data for both scoring methods and both systems with high correlation coefficients (R2>0.94). System A showed more significant differences between spiral and sequential mode than system B. Almost no differences were observed in scanners of the same manufacturer for the Agatston score and no differences for the Volume score. The deviations of the Agatston and Volume scores showed regression dependencies approximately equal to the square root of the absolute score. The Agatston and Volume scores obtained with 64-slice MDCT imaging are highly correlated with EBT-obtained scores but are significantly underestimated (-10% to -2%) for both sequential and spiral acquisition modes. System B is more independent of acquisition mode to calcium score than system A. The Volume score shows no intramanufacturer dependency and its use is advocated versus the Agatston score. Using the same cut points for MDCT-based calcium scores as for EBT-based calcium scores can result in classifying individuals into a too low risk category. System information and scanprotocol is therefore needed for every calcium score procedure to ensure a correct clinical interpretation of the obtained calcium score results.

A model for quantitative correction of coronary calcium scores on multidetector, dual source, and electron beam computed tomography for influences of linear motion, calcification density, and temporal resolution: a cardiac phantom study.

PURPOSE The objective of this study is to quantify the influence of linear motion, calcification density, and temporal resolution on coronary calcium determination using multidetector computed tomography (MDCT), dual source CT (DSCT), and electron beam tomography (EBT) and to find a quantitative method which corrects for the influences of these parameters using a linear moving cardiac phantom. METHODS On a robotic arm with artificial arteries with four calcifications of increasing density, a linear movement was applied between 0 and 120 mm/s (step of 10 mm/s). The phantom was scanned five times on 64-slice MDCT, DSCT, and EBT using a standard acquisition protocol. The average Agatston, volume, and mass scores were determined for each velocity, calcification, and scanner. Susceptibility to motion was quantified using a cardiac motion susceptibility (CMS) index. Resemblance to EBT and physical volume and mass was quantified using aΔ index. RESULTS Increasing motion artifacts were observed at increasing velocities on all scanners, with increasing severity from EBT to DSCT to 64-slice MDCT. The calcium score showed a linear dependency on motion from which a correction factor could be derived. This correction factor showed a linear dependency on the mean calcification density with a good fit for all three scoring methods and all three scanners (0.73≤R2≤0.95). The slope and offset of this correction factor showed a linear dependency on temporal resolution with a good fit for all three scoring methods and all three scanners (0.83≤R2≤0.98). CMS was minimal for EBT and increasing values were observed for DSCT and highest values for 64-slice MDCT. CMS was minimal for mass score and increasing values were observed for volume score and highest values for Agatston score. For all densities and scoring methods DSCT showed on average the closest resemblance to EBT calcium scores. When using the correction factor, CMS index decreased on average by 15% and Δ index decreased by 35%. CONCLUSIONS Calcium scores determined on DSCT and 64-slice MDCT are highly susceptible to motion as compared to EBT. The mass score is less susceptible to motion compared to volume and Agatston score. Calcium scores determined on DSCT bear a closer resemblance to EBT obtained calcium scores than 64-slice MDCT. In addition, the calcium score is highly dependent on the average density of individual calcifications and the dependency of the calcium score on motion showed a linear behavior on calcification density. From these relations, a quantitative method could be derived which corrects the measured calcium score for the influence of linear motion, mean calcification density, and temporal resolution.

Semi-automated quantitative intravoxel incoherent motion analysis and its implementation in breast diffusion-weighted imaging

To optimize and validate intravoxel incoherent motion (IVIM) modeled diffusion‐weighted imaging (DWI) compared with the apparent diffusion coefficient (ADC) for semi‐automated analysis of breast lesions using a multi‐reader setup.

Quantitative DWI implemented after DCE-MRI yields increased specificity for BI-RADS 3 and 4 breast lesions.

To assess if specificity can be increased when semiautomated breast lesion analysis of quantitative diffusion‐weighted imaging (DWI) is implemented after dynamic contrast‐enhanced (DCE‐) magnetic resonance imaging (MRI) in the workup of BI‐RADS 3 and 4 breast lesions larger than 1 cm.

Threshold adjusted calcium scoring using CT is less susceptible to cardiac motion and more accurate.

The purpose of this paper is to investigate calcium scoring on computed tomography (CT) using an adjusted threshold depending on the maximum Hounsfield value within the calcification (HU(peak)). The volume of 19 calcifications was retrospectively determined on 64-slice multidetector CT and dual source CT (DSCT) at different thresholds and the threshold associated with the physical volume was determined. In addition, approximately 10 000 computer simulations were done simulating the same process for calcifications with mixed density. Using these data a relation between the HU(peak) and the threshold could be established. Hereafter, this relation was assessed by scanning six calcifications in a phantom at 40-110 beats per minute using DSCT. The influence of motion was determined and the measured calcium scores were compared to the physical volumes and mass. A positive linear correlation was found between the scoring threshold and the HU(peak) of the calcifications both for the phantom measurements as for the computer simulations. Using this relation the individual threshold for each calcification could be calculated. Calcium scores of the moving calcifications determined with an adjusted threshold were approximately 30% less susceptible to cardiac motion compared to standard calcium scoring. Furthermore, these scores approximated the physical volume and mass at least 10% better than the standard calcium scores. The threshold in calcium scoring should be adjusted for each individual calcification based on the HU(peak) of the calcification. Calcium scoring using an adjusted threshold is less susceptible to cardiac motion and more accurate compared to the physical values.

Diminished liver microperfusion in Fontan patients : A biexponential DWI study

It has been demonstrated that hepatic apparent diffusion coefficients (ADC) are decreasing in patients with a Fontan circulation. It remains however unclear whether this is a true decrease of molecular diffusion, or rather reflects decreased microperfusion due to decreased portal blood flow. The purpose of this study was therefore to differentiate diffusion and microperfusion using intravoxel incoherent motion (IVIM) modeled diffusion-weighted imaging (DWI) for different liver segments in patients with a Fontan circulation, compare to a control group, and relate with liver function, chronic hepatic congestion and hepatic disease. For that purpose, livers of 59 consecutively included patients with Fontan circulation (29 men; mean-age, 19.1 years) were examined (Oct 2012─Dec 2013) with 1.5T MRI and DWI (b = 0,50,100,250,500,750,1500,1750 s/mm2). IVIM (Dslow, Dfast, ffast) and ADC were calculated for eight liver segments, compared to a control group (19 volunteers; 10 men; mean-age, 32.9 years), and correlated to follow-up duration, clinical variables, and laboratory measurements associated with liver function. The results demonstrated that microperfusion was reduced (p<0.001) in Fontan livers compared to controls with ─38.1% for Dfast and ─32.6% for ffast. Molecular diffusion (Dslow) was similar between patients and controls, while ADC was significantly lower (─14.3%) in patients (p<0.001). ADC decreased significantly with follow-up duration after Fontan operation (r = ─0.657). Dslow showed significant inverse correlations (r = ─0.591) with follow-up duration whereas Dfast and ffast did not. From these results it was concluded that the decreasing ADC values in Fontan livers compared with controls reflect decreases in hepatic microperfusion rather than any change in molecular diffusion. However, with the time elapsed since the Fontan operation molecular diffusion and ADC decreased while microperfusion remained stable. This indicates that after Fontan operation initial blood flow effects on the liver are followed by intracellular changes preceding the formation of fibrosis and cirrhosis.

The Fontan circulation and the liver: A magnetic resonance diffusion-weighted imaging study.

BACKGROUND Patients with a Fontan circulation tend to develop liver fibrosis, liver cirrhosis and even hepatocellular carcinoma. The aim of this study is to use the magnetic resonance technique diffusing-weighted imaging (DWI) for detecting liver fibrosis/cirrhosis in Fontan patients and to establish whether DWI results are associated with functional aspects of the Fontan circulation. METHODS In a cross-sectional study, 59 Fontan patients were evaluated by liver DWI. The association between apparent diffusion coefficients (ADC) and patient characteristics, laboratory measurements and functional aspects of the Fontan circulation (NYHA class, maximum oxygen uptake during exercise and cardiac index) was assessed. RESULTS Liver ADC values were low (0.82×10(-3)±0.11×10(-3) mm2/s) compared with literature values for healthy volunteers and correlated negatively with calculated liver fibrosis/cirrhosis scores (Fib-4 score, p=0.019; AST/ALT ratio, p=0.009) and gamma-glutamyl transferase (p=0.001). Furthermore, ADC values correlated negatively with follow-up duration (p<0.001) and positively with cardiac index (p=0.019). No correlation between ADC values and exercise tests was found. In multivariable analysis, the ADC values were independently correlated with follow-up duration after Fontan completion. CONCLUSIONS The results of the current study suggest that progressive liver damage due to chronic congestion and potential hypoperfusion is reflected in the liver ADC values in Fontan patients. This study highlights that liver damage in the context of the Fontan circulation might be far more common than previously thought, and that the implementation of liver assessment in the routine follow-up of Fontan patients is recommendable.

Quantitative Myocardial Perfusion with Dynamic Contrast-Enhanced Imaging in MRI and CT: Theoretical Models and Current Implementation

Technological advances in magnetic resonance imaging (MRI) and computed tomography (CT), including higher spatial and temporal resolution, have made the prospect of performing absolute myocardial perfusion quantification possible, previously only achievable with positron emission tomography (PET). This could facilitate integration of myocardial perfusion biomarkers into the current workup for coronary artery disease (CAD), as MRI and CT systems are more widely available than PET scanners. Cardiac PET scanning remains expensive and is restricted by the requirement of a nearby cyclotron. Clinical evidence is needed to demonstrate that MRI and CT have similar accuracy for myocardial perfusion quantification as PET. However, lack of standardization of acquisition protocols and tracer kinetic model selection complicates comparison between different studies and modalities. The aim of this overview is to provide insight into the different tracer kinetic models for quantitative myocardial perfusion analysis and to address typical implementation issues in MRI and CT. We compare different models based on their theoretical derivations and present the respective consequences for MRI and CT acquisition parameters, highlighting the interplay between tracer kinetic modeling and acquisition settings.

Intermodel Agreement of Myocardial Blood Flow Estimation From Stress-Rest Myocardial Perfusion Magnetic Resonance Imaging in Patients With Coronary Artery Disease

ObjectivesThe aim of this study was to assess the intermodel agreement of different magnetic resonance myocardial perfusion models and evaluate their correspondence to stenosis diameter. Materials and MethodsIn total, 260 myocardial segments were analyzed from rest and adenosine stress first-pass myocardial perfusion magnetic resonance images (1.5 T, 0.050 ± 0.005 mmol/kg body weight gadolinium; 122 segments in rest, 138 in stress) in 10 patients with suspected or known coronary artery disease. Signal intensity curves were calculated per myocardial segment, of which the contours were traced with QMASS MR V.7.6 (Medis, Leiden, the Netherlands), and exported to Matlab. Myocardial blood flow quantification was performed with distributed parameter, extended Toft, Patlak, and Fermi parametric models (in-house programs; Matlab R2013a; Mathworks Inc, Natick, MA). Modeling was applied after the signal intensity curves were corrected for spatial magnetic field inhomogeneity and contrast saturation. Overall and grouped perfusion values based on presence of coronary stenosis (>50% diameter reduction) at coronary computed tomography angiography at second generation dual-source computed tomography were compared between the perfusion models. ResultsRest and stress myocardial perfusion estimates for all models were significantly related to each other (P < 0.001). The highest correlation coefficients were found between the extended Toft and Fermi models (R = 0.89−0.91) and low correlation coefficients between the distributed parameter and Patlak models (R = 0.66−0.68). The models resulted in significantly different perfusion estimates in stress (P = 0.03), but not in rest (P = 0.74). The differences in perfusion estimates in stress were caused by differences between the distributed parameter and Patlak models and between the Patlak and Fermi models (both P < 0.001). Significantly lower perfusion estimates were found for myocardial segments subtended by coronary arteries with versus without significant stenosis, but only for estimations produced by the extended Toft model (P = 0.04) and Fermi model (P = 0.01). There were no significant differences in rest perfusion values between models. ConclusionsQuantitative myocardial perfusion values in stress depend on the modeling method used to calculate the perfusion estimate. The difference in myocardial perfusion estimate with or without stenosis in the subtending coronary artery is most pronounced when the extended Toft or Fermi model is used.