Baseline PET/CT imaging parameters for prediction of treatment outcome in Hodgkin and diffuse large B cell lymphoma: a systematic review.

PURPOSE: To systematically review the literature evaluating clinical utility of imaging metrics derived from baseline fluorine-18 fluorodeoxyglucose positron emission tomography/computed tomography (PET/CT) for prediction of progression-free (PFS) and overall survival (OS) in patients with classical Hodgkin lymphoma (HL) and diffuse large B cell lymphoma (DLBCL). METHODS: A search of MEDLINE/PubMed, Web of Science, Cochrane, Scopus and clinicaltrials.gov databases was undertaken for articles evaluating PET/CT imaging metrics as outcome predictors in HL and DLBCL. PRISMA guidelines were followed. Risk of bias was assessed using the Quality in Prognosis Studies (QUIPS) tool. RESULTS: Forty-one articles were included (31 DLBCL, 10 HL). Significant predictive ability was reported in 5/20 DLBCL studies assessing SUVmax (PFS: HR 0.13-7.35, OS: HR 0.83-11.23), 17/19 assessing metabolic tumour volume (MTV) (PFS: HR 2.09-11.20, OS: HR 2.40-10.32) and 10/13 assessing total lesion glycolysis (TLG) (PFS: HR 1.078-11.21, OS: HR 2.40-4.82). Significant predictive ability was reported in 1/4 HL studies assessing SUVmax (HR not reported), 6/8 assessing MTV (PFS: HR 1.2-10.71, OS: HR 1.00-13.20) and 2/3 assessing TLG (HR not reported). There are 7/41 studies assessing the use of radiomics (4 DLBCL, 2 HL); 5/41 studies had internal validation and 2/41 included external validation. All studies had overall moderate or high risk of bias. CONCLUSION: Most studies are retrospective, underpowered, heterogenous in their methodology and lack external validation of described models. Further work in protocol harmonisation, automated segmentation techniques and optimum performance cut-off is required to develop robust methodologies amenable for clinical utility.

Tissue microarray (TMA) use in post mortem neuropathology.

BACKGROUND: Tissue microarrays (TMAs), where each block (and thus section) contains multiple tissue cores from multiple blocks potentially allow more efficient use of tissue, reagents and time in neuropathology. NEW METHOD: The relationship between data from TMA cores and whole sections was investigated using 'virtual' TMA cores. This involved quantitative assessments of microglial pathology in white matter lesions and motor neuron disease, alongside qualitative TDP-43 inclusion status in motor neuron disease cases. Following this, a protocol was developed for TMA construction. RESULTS: For microglial pathology we found good concordance between virtual cores and whole sections for volume density using one 1.75 mm core (equivalent to a 2 mm core after accounting for peripheral tissue loss). More sophisticated microglial cell size and measures required two cores. Qualitative results of pTDP-43 pathology showed use of one 1.75 mm core gave a 100 % sensitivity and specificity within grey matter, and 88.3 % sensitivity and 100 % specificity within white matter. A method of producing the TMAs was suitable for immunohistochemistry both manually and by autostainer, with the minimal core loss from the microscope slide. COMPARISON WITH EXISTING METHODS: TMAs have been used infrequently in post mortem neuropathology research. However, we believe TMAs give comparable tissue assessment results and can be constructed, sectioned and stained with relative ease. CONCLUSIONS: We found TMAs could be used to assess both quantitative (microglial pathology) and qualitative pathology (TDP-43 proteinopathy) with greatly reduced quantities of tissue, time and reagents. These could be used for further work to improve data acquisition efficiency.

Highly integrated workflows for exploring cardiovascular conditions: Exemplars of precision medicine in Alzheimer's disease and aortic dissection.

For precision medicine to be implemented through the lens of in silico technology, it is imperative that biophysical research workflows offer insight into treatments that are specific to a particular illness and to a particular subject. The boundaries of precision medicine can be extended using multiscale, biophysics-centred workflows that consider the fundamental underpinnings of the constituents of cells and tissues and their dynamic environments. Utilising numerical techniques that can capture the broad spectrum of biological flows within complex, deformable and permeable organs and tissues is of paramount importance when considering the core prerequisites of any state-of-the-art precision medicine pipeline. In this work, a succinct breakdown of two precision medicine pipelines developed within two Virtual Physiological Human (VPH) projects are given. The first workflow is targeted on the trajectory of Alzheimer's Disease, and caters for novel hypothesis testing through a multicompartmental poroelastic model which is integrated with a high throughput imaging workflow and subject-specific blood flow variability model. The second workflow gives rise to the patient specific exploration of Aortic Dissections via a multi-scale and compliant model, harnessing imaging, computational fluid-dynamics (CFD) and dynamic boundary conditions. Results relating to the first workflow include some core outputs of the multiporoelastic modelling framework, and the representation of peri-arterial swelling and peri-venous drainage solution fields. The latter solution fields were statistically analysed for a cohort of thirty-five subjects (stratified with respect to disease status, gender and activity level). The second workflow allowed for a better understanding of complex aortic dissection cases utilising both a rigid-wall model informed by minimal and clinically common datasets as well as a moving-wall model informed by rich datasets.

Wall shear stress at the initiation site of cerebral aneurysms.

Hemodynamics are believed to play an important role in the initiation of cerebral aneurysms. In particular, studies have focused on wall shear stress (WSS), which is a key regulator of vascular biology and pathology. In line with the observation that aneurysms predominantly occur at regions of high WSS, such as bifurcation apices or outer walls of vascular bends, correlations have been found between the aneurysm initiation site and high WSS. The aim of our study was to analyze the WSS field at an aneurysm initiation site that was neither a bifurcation apex nor the outer wall of a vascular bend. Ten cases with aneurysms on the A1 segment of the anterior cerebral artery were analyzed and compared with ten controls. Aneurysms were virtually removed from the vascular models of the cases to mimic the pre-aneurysm geometry. Computational fluid dynamics (CFD) simulations were created to assess the magnitude, gradient, multidirectionality, and pulsatility of the WSS. To aid the inter-subject comparison of hemodynamic variables, we mapped the branch surfaces onto a two-dimensional parametric space. This approach made it possible to view the whole branch at once for qualitative evaluation. It also allowed us to empirically define a patch for quantitative analysis, which was consistent among subjects and encapsulated the aneurysm initiation sites in our dataset. To test the sensitivity of our results, CFD simulations were repeated with a second independent observer virtually removing the aneurysms and with a 20 % higher flow rate at the inlet. We found that branches harboring aneurysms were characterized by high WSS and high WSS gradients. Among all assessed variables, the aneurysm initiation site most consistently coincided with peaks of temporal variation in the WSS magnitude.

Improved hybrid/GPU algorithm for solving cardiac electrophysiology problems on Purkinje networks.

Cardiac Purkinje fibers provide an important pathway to the coordinated contraction of the heart. We present a numerical algorithm for the solution of electrophysiology problems across the Purkinje network that is efficient enough to be used in in silico studies on realistic Purkinje networks with physiologically detailed models of ion exchange at the cell membrane. The algorithm is on the basis of operator splitting and is provided with 3 different implementations: pure CPU, hybrid CPU/GPU, and pure GPU. Compared to our previous work, we modify the explicit gap junction term at network bifurcations to improve its mathematical consistency. Due to this improved consistency of the model, we are able to perform an empirical convergence study against analytical solutions. The study verified that all 3 implementations produce equivalent convergence rates, and shows that the algorithm produces equivalent result across different hardware platforms. Finally, we compare the efficiency of all 3 implementations on Purkinje networks of increasing spatial resolution using membrane models of increasing complexity. Both hybrid and pure GPU implementations outperform the pure CPU implementation, but their relative performance difference depends on the size of the Purkinje network and the complexity of the membrane model used.

Integration of Multi-Plane Tissue Doppler and B-Mode Echocardiographic Images for Left Ventricular Motion Estimation.

Although modern ultrasound acquisition systems allow recording of 3D echocardiographic images, tracking anatomical structures from them is still challenging. In addition, since these images are typically created from information obtained across several cardiac cycles, it is not yet possible to acquire high-quality 3D images from patients presenting varying heart rhythms. In this paper, we propose a method to estimate the motion field from multi-plane echocardiographic images of the left ventricle, which are acquired simultaneously during a single cardiac cycle. The method integrates tri-plane B-mode and tissue Doppler images acquired at different rotation angles around the long axis of the left ventricle. It uses a diffeomorphic continuous spatio-temporal transformation model with a spherical data representation for a better interpolation in the circumferential direction. This framework allows exploiting the spatial relation among the acquired planes. In addition, higher temporal resolution of the transformation in the beam direction is achieved by uncoupling the estimation of the different components of the velocity field. The method was validated using a realistic synthetic dataset including healthy and ischemic cases, obtaining errors of 0.14 ± 0.09 mm for displacements, 0.96 ± 1.03% for longitudinal strain and 3.94 ± 4.38% for radial strain estimation. In addition, the method was also demonstrated on a healthy volunteer and two patients with ischemia.

Accuracy and reproducibility of patient-specific hemodynamic models of stented intracranial aneurysms: report on the Virtual Intracranial Stenting Challenge 2011.

Validation studies are prerequisites for computational fluid dynamics (CFD) simulations to be accepted as part of clinical decision-making. This paper reports on the 2011 edition of the Virtual Intracranial Stenting Challenge. The challenge aimed to assess the reproducibility with which research groups can simulate the velocity field in an intracranial aneurysm, both untreated and treated with five different configurations of high-porosity stents. Particle imaging velocimetry (PIV) measurements were obtained to validate the untreated velocity field. Six participants, totaling three CFD solvers, were provided with surface meshes of the vascular geometry and the deployed stent geometries, and flow rate boundary conditions for all inlets and outlets. As output, they were invited to submit an abstract to the 8th International Interdisciplinary Cerebrovascular Symposium 2011 (ICS'11), outlining their methods and giving their interpretation of the performance of each stent configuration. After the challenge, all CFD solutions were collected and analyzed. To quantitatively analyze the data, we calculated the root-mean-square error (RMSE) over uniformly distributed nodes on a plane slicing the main flow jet along its axis and normalized it with the maximum velocity on the slice of the untreated case (NRMSE). Good agreement was found between CFD and PIV with a NRMSE of 7.28%. Excellent agreement was found between CFD solutions, both untreated and treated. The maximum difference between any two groups (along a line perpendicular to the main flow jet) was 4.0 mm/s, i.e. 4.1% of the maximum velocity of the untreated case, and the average NRMSE was 0.47% (range 0.28-1.03%). In conclusion, given geometry and flow rates, research groups can accurately simulate the velocity field inside an intracranial aneurysm-as assessed by comparison with in vitro measurements-and find excellent agreement on the hemodynamic effect of different stent configurations.

Leptin may play a role in bone microstructural alterations in obese children.

CONTEXT: Bone mass is low and fracture risk is higher in obese children. Hormonal changes in relation to skeletal microstructure and biomechanics have not been studied in obese children. OBJECTIVE: The objective of the study was to ascertain the relationships of obesity-related changes in hormones with skeletal microstructure and biomechanics. DESIGN: High resolution peripheral quantitative computed tomography (HR-pQCT) was used to compare three-dimensional cortical and trabecular microstructure and biomechanics at load-bearing and nonload bearing sites in obese and lean children. The relationship between leptin, adiponectin, testosterone, estrogen, osteocalcin and sclerostin and skeletal microstructure was also determined. SETTING: The study was conducted at a tertiary pediatric endocrine unit in the United Kingdom. PARTICIPANTS: Obese and lean children were matched by gender and pubertal stage. RESULTS: Radial cortical porosity (mean difference -0.01 [95% CI: -0.02, -0.004], P = .003) and cortical pore diameter (mean difference -0.005 mm [95% CI: -0.009, -0.001], P = .011) were lower in obese children. Tibial trabecular thickness was lower (mean difference -0.009 mm [95% CI: -0.014, -0.004], P = .003), and trabecular number was higher (mean difference 0.23 mm(-1) [95% CI: 0.08, 0.38], P = .004) in obese children. At the radius, fat mass percentage negatively correlated with cortical porosity (r = -0.57, P < .001) and pore diameter (r = -0.38, P = .02) and negatively correlated with trabecular thickness (r = -0.62, P < .001) and trabecular von Mises stress (r = -0.39, P = .019) at the tibia. No difference was observed in the other biomechanical parameters of the radius and tibia. Leptin was higher in obese children (805.3 ± 440.6 pg/ml vs 98.1 ± 75.4 pg/ml, P < .001) and was inversely related to radial cortical porosity (r = 0.60, 95% CI: [-0.80, -0.30], P < .001), radial cortical pore diameter (r = 0.51, 95% CI [-0.75, -0.16], P = .002), tibial trabecular thickness (r = 0.55, 95% CI: [-0.78, -0.21], P = .001) and tibial trabecular von Mises stress (r = -0.39, 95% CI: -0.65, 0.04, P = .02). CONCLUSION: Childhood obesity alters radial and tibial microstructure. Leptin may direct these changes. Despite this, the biomechanical properties of the radius and tibia do not adapt sufficiently in obese children to withstand the increased loading potential from a fall. This may explain the higher incidence of fracture in obese children.

Improved myocardial motion estimation combining tissue Doppler and B-mode echocardiographic images.

We propose a technique for myocardial motion estimation based on image registration using both B-mode echocardiographic images and tissue Doppler sequences acquired interleaved. The velocity field is modeled continuously using B-splines and the spatiotemporal transform is constrained to be diffeomorphic. Images before scan conversion are used to improve the accuracy of the estimation. The similarity measure includes a model of the speckle pattern distribution of B-mode images. It also penalizes the disagreement between tissue Doppler velocities and the estimated velocity field. Registration accuracy is evaluated and compared to other alternatives using a realistic synthetic dataset, obtaining mean displacement errors of about 1 mm. Finally, the method is demonstrated on data acquired from six volunteers, both at rest and during exercise. Robustness is tested against low image quality and fast heart rates during exercise. Results show that our method provides a robust motion estimate in these situations.

Approximating hemodynamics of cerebral aneurysms with steady flow simulations.

Computational fluid dynamics (CFD) simulations can be employed to gain a better understanding of hemodynamics in cerebral aneurysms and improve diagnosis and treatment. However, introduction of CFD techniques into clinical practice would require faster simulation times. The aim of this study was to evaluate the use of computationally inexpensive steady flow simulations to approximate the aneurysm's wall shear stress (WSS) field. Two experiments were conducted. Experiment 1 compared for two cases the time-averaged (TA), peak systole (PS) and end diastole (ED) WSS field between steady and pulsatile flow simulations. The flow rate waveform imposed at the inlet was varied to account for variations in heart rate, pulsatility index, and TA flow rate. Consistently across all flow rate waveforms, steady flow simulations accurately approximated the TA, but not the PS and ED, WSS field. Following up on experiment 1, experiment 2 tested the result for the TA WSS field in a larger population of 20 cases covering a wide range of aneurysm volumes and shapes. Steady flow simulations approximated the space-averaged WSS with a mean error of 4.3%. WSS fields were locally compared by calculating the absolute error per node of the surface mesh. The coefficient of variation of the root-mean-square error over these nodes was on average 7.1%. In conclusion, steady flow simulations can accurately approximate the TA WSS field of an aneurysm. The fast computation time of 6 min per simulation (on 64 processors) could help facilitate the introduction of CFD into clinical practice.

Benchmarking framework for myocardial tracking and deformation algorithms: an open access database.

In this paper we present a benchmarking framework for the validation of cardiac motion analysis algorithms. The reported methods are the response to an open challenge that was issued to the medical imaging community through a MICCAI workshop. The database included magnetic resonance (MR) and 3D ultrasound (3DUS) datasets from a dynamic phantom and 15 healthy volunteers. Participants processed 3D tagged MR datasets (3DTAG), cine steady state free precession MR datasets (SSFP) and 3DUS datasets, amounting to 1158 image volumes. Ground-truth for motion tracking was based on 12 landmarks (4 walls at 3 ventricular levels). They were manually tracked by two observers in the 3DTAG data over the whole cardiac cycle, using an in-house application with 4D visualization capabilities. The median of the inter-observer variability was computed for the phantom dataset (0.77 mm) and for the volunteer datasets (0.84 mm). The ground-truth was registered to 3DUS coordinates using a point based similarity transform. Four institutions responded to the challenge by providing motion estimates for the data: Fraunhofer MEVIS (MEVIS), Bremen, Germany; Imperial College London - University College London (IUCL), UK; Universitat Pompeu Fabra (UPF), Barcelona, Spain; Inria-Asclepios project (INRIA), France. Details on the implementation and evaluation of the four methodologies are presented in this manuscript. The manually tracked landmarks were used to evaluate tracking accuracy of all methodologies. For 3DTAG, median values were computed over all time frames for the phantom dataset (MEVIS=1.20mm, IUCL=0.73 mm, UPF=1.10mm, INRIA=1.09 mm) and for the volunteer datasets (MEVIS=1.33 mm, IUCL=1.52 mm, UPF=1.09 mm, INRIA=1.32 mm). For 3DUS, median values were computed at end diastole and end systole for the phantom dataset (MEVIS=4.40 mm, UPF=3.48 mm, INRIA=4.78 mm) and for the volunteer datasets (MEVIS=3.51 mm, UPF=3.71 mm, INRIA=4.07 mm). For SSFP, median values were computed at end diastole and end systole for the phantom dataset(UPF=6.18 mm, INRIA=3.93 mm) and for the volunteer datasets (UPF=3.09 mm, INRIA=4.78 mm). Finally, strain curves were generated and qualitatively compared. Good agreement was found between the different modalities and methodologies, except for radial strain that showed a high variability in cases of lower image quality.

Personalization of a cardiac electromechanical model using reduced order unscented Kalman filtering from regional volumes.

Patient-specific cardiac modeling can help in understanding pathophysiology and therapy planning. However it requires to combine functional and anatomical data in order to build accurate models and to personalize the model geometry, kinematics, electrophysiology and mechanics. Personalizing the electromechanical coupling from medical images is a challenging task. We use the Bestel-Clément-Sorine (BCS) electromechanical model of the heart, which provides reasonable accuracy with a reasonable number of parameters (14 for each ventricle) compared to the available clinical data at the organ level. We propose a personalization strategy from cine MRI data in two steps. We first estimate global parameters with an automatic calibration algorithm based on the Unscented Transform which allows to initialize the parameters while matching the volume and pressure curves. In a second step we locally personalize the contractilities of all AHA (American Heart Association) zones of the left ventricle using the reduced order unscented Kalman filtering on Regional Volumes. This personalization strategy was validated synthetically and tested successfully on eight healthy and three pathological cases.

Understanding the mechanisms amenable to CRT response: from pre-operative multimodal image data to patient-specific computational models.

This manuscript describes our recent developments towards better understanding of the mechanisms amenable to cardiac resynchronization therapy response. We report the results from a full multimodal dataset corresponding to eight patients from the euHeart project. The datasets include echocardiography, MRI and electrophysiological studies. We investigate two aspects. The first one focuses on pre-operative multimodal image data. From 2D echocardiography and 3D tagged MRI images, we compute atlas based dyssynchrony indices. We complement these indices with presence and extent of scar tissue and correlate them with CRT response. The second one focuses on computational models. We use pre-operative imaging to generate a patient-specific computational model. We show results of a fully automatic personalized electromechanical simulation. By case-per-case discussion of the results, we highlight the potential and key issues of this multimodal pipeline for the understanding of the mechanisms of CRT response and a better patient selection.

Generating anatomical models of the heart and the aorta from medical images for personalized physiological simulations.

The anatomy and motion of the heart and the aorta are essential for patient-specific simulations of cardiac electrophysiology, wall mechanics and hemodynamics. Within the European integrated project euHeart, algorithms have been developed that allow to efficiently generate patient-specific anatomical models from medical images from multiple imaging modalities. These models, for instance, account for myocardial deformation, cardiac wall motion, and patient-specific tissue information like myocardial scar location. Furthermore, integration of algorithms for anatomy extraction and physiological simulations has been brought forward. Physiological simulations are linked closer to anatomical models by encoding tissue properties, like the muscle fibers, into segmentation meshes. Biophysical constraints are also utilized in combination with image analysis to assess tissue properties. Both examples show directions of how physiological simulations could provide new challenges and stimuli for image analysis research in the future.

RADStation3G: a platform for cardiovascular image analysis integrating PACS, 3D+t visualization and grid computing.

RADStation3G is a software platform for cardiovascular image analysis and surgery planning. It provides image visualization and management in 2D, 3D and 3D+t; data storage (images or operational results) in a PACS (using DICOM); and exploitation of patients' data such as images and pathologies. Further, it provides support for computationally expensive processes with grid technology. In this article we first introduce the platform and present a comparison with existing systems, according to the platform's modules (for cardiology, angiology, PACS archived enriched searching and grid computing), and then RADStation3G is described in detail.

Intra-aneurysmal pressure and flow changes induced by flow diverters: relation to aneurysm size and shape.

BACKGROUND AND PURPOSE: Effects of blood flow modification by flow diverters are observed to lead often to aneurysm thrombosis and reverse remodeling. For this process, to further understand the potential roles of intra-aneurysmal blood pressure changes and aneurysm morphologies, 23 patients were studied by numeric simulation. MATERIALS AND METHODS: 3D imaging of aneurysms of different sizes and shapes, all located at the supraclinoid segment of the ICA (n=23), was prepared for CFD simulations. Hemodynamic variables were calculated for conditions before and after virtual FD implantation, reconstituting a vessel wall scaffold across the aneurysm neck. WSS, velocity, residence time, turnover time, and intra-aneurysmal pressure were assessed statistically. RESULTS: After placement of FDs, significant reductions inside the aneurysm were observed for most hemodynamic variables (P<.01) except mean intra-aneurysmal pressures. For minimum/maximum intra-aneurysmal pressure values, small but significant changes were found; however, they were considered too small to be of relevance. CONCLUSIONS: Calculations in 23 cases did not reveal significant intra-aneurysmal mean or peak pressure changes, indicating a minor role of pressure changes in the rare event of secondary ruptures after FD use. Other hemodynamic variables (WSS and velocity) exhibited more significant changes, indicating their role in intra-aneurysmal thrombus formation. Size-dependent, significantly higher reduction in WSS (P=.069) and velocity (P=.013) was observed in small aneurysms compared with larger ones. When it came to shape, there were significantly higher reductions in WSS (P=.055) and velocity (P=.065) and a significantly higher increase in turnover time in fusiform aneurysms compared with saccular aneurysms.

Cardiac motion estimation by joint alignment of tagged MRI sequences.

Image registration has been proposed as an automatic method for recovering cardiac displacement fields from tagged Magnetic Resonance Imaging (tMRI) sequences. Initially performed as a set of pairwise registrations, these techniques have evolved to the use of 3D+t deformation models, requiring metrics of joint image alignment (JA). However, only linear combinations of cost functions defined with respect to the first frame have been used. In this paper, we have applied k-Nearest Neighbors Graphs (kNNG) estimators of the α-entropy (H(α)) to measure the joint similarity between frames, and to combine the information provided by different cardiac views in an unified metric. Experiments performed on six subjects showed a significantly higher accuracy (p<0.05) with respect to a standard pairwise alignment (PA) approach in terms of mean positional error and variance with respect to manually placed landmarks. The developed method was used to study strains in patients with myocardial infarction, showing a consistency between strain, infarction location, and coronary occlusion. This paper also presents an interesting clinical application of graph-based metric estimators, showing their value for solving practical problems found in medical imaging.

Deployment of self-expandable stents in aneurysmatic cerebral vessels: comparison of different computational approaches for interventional planning.

In the last few years, there has been a growing focus on faster computational methods to support clinicians in planning stenting procedures. This study investigates the possibility of introducing computational approximations in modelling stent deployment in aneurysmatic cerebral vessels to achieve simulations compatible with the constraints of real clinical workflows. The release of a self-expandable stent in a simplified aneurysmatic vessel was modelled in four different initial positions. Six progressively simplified modelling approaches (based on Finite Element method and Fast Virtual Stenting--FVS) have been used. Comparing accuracy of the results, the final configuration of the stent is more affected by neglecting mechanical properties of materials (FVS) than by adopting 1D instead of 3D stent models. Nevertheless, the differences showed are acceptable compared to those achieved by considering different stent initial positions. Regarding computational costs, simulations involving 1D stent features are the only ones feasible in clinical context.