Improving In Vivo High-Resolution CT Imaging of the Tumour Vasculature in Xenograft Mouse Models through Reduction of Motion and Bone-Streak Artefacts

Performance evaluation of image co-registration methods in photoacoustic mesoscopy of the vasculature

Objective:The formation of functional vasculature in solid tumours enables delivery of oxygen and nutrients, and is vital for effective treatment with chemotherapeutic agents. Longitudinal characterisation of vascular networks can be enabled using mesoscopic photoacoustic imaging, but requires accurate image co-registration to precisely assess local changes across disease development or in response to therapy. Co-registration in photoacoustic imaging is challenging due to the complex nature of the generated signal, including the sparsity of data, artefacts related to the illumination/detection geometry, scan-to-scan technical variability, and biological variability, such as transient changes in perfusion. To better inform the choice of co-registration algorithms, we compared five open-source methods, in physiological and pathological tissues, with the aim of aligning evolving vascular networks in tumours imaged over growth at different time-points.Approach:Co-registration techniques were applied to 3D vascular images acquired with photoacoustic mesoscopy from murine ears and breast cancer patient-derived xenografts, at a fixed time-point and longitudinally. Images were pre-processed and segmented using an unsupervised generative adversarial network. To compare co-registration quality in different settings, pairs of fixed and moving intensity images and/or segmentations were fed into five methods split into the following categories: affine intensity-based using 1)mutual information (MI) or 2)normalised cross-correlation (NCC) as optimisation metrics, affine shape-based using 3)NCC applied to distance-transformed segmentations or 4)iterative closest point algorithm, and deformable weakly supervised deep learning-based using 5)LocalNet co-registration. Percent-changes in Dice coefficients, surface distances, MI, structural similarity index measure and target registration errors were evaluated.Main results:Co-registration using MI or NCC provided similar alignment performance, better than shape-based methods. LocalNet provided accurate co-registration of substructures by optimising subfield deformation throughout the volumes, outperforming other methods, especially in the longitudinal breast cancer xenograft dataset by minimising target registration errors.Significance:We showed the feasibility of co-registering repeatedly or longitudinally imaged vascular networks in photoacoustic mesoscopy, taking a step towards longitudinal quantitative characterisation of these complex structures. These tools open new outlooks for monitoring tumour angiogenesis at the meso-scale and for quantifying treatment-induced co-localised alterations in the vasculature.

Functional biomarkers derived from computed tomography and magnetic resonance imaging differentiate PDAC subgroups and reveal gemcitabine-induced hypo-vascularization

PURPOSE

Pancreatic ductal adenocarcinoma (PDAC) is a molecularly heterogeneous tumor entity with no clinically established imaging biomarkers. We hypothesize that tumor morphology and physiology, including vascularity and perfusion, show variations that can be detected by differences in contrast agent (CA) accumulation measured non-invasively. This work seeks to establish imaging biomarkers for tumor stratification and therapy response monitoring in PDAC, based on this hypothesis.

METHODS AND MATERIALS

Regional CA accumulation in PDAC was correlated with tumor vascularization, stroma content, and tumor cellularity in murine and human subjects. Changes in CA distribution in response to gemcitabine (GEM) were monitored longitudinally with computed tomography (CT) Hounsfield Units ratio (HUr) of tumor to the aorta or with magnetic resonance imaging (MRI) ΔR1 area under the curve at 60 s tumor-to-muscle ratio (AUC60r). Tissue analyses were performed on co-registered samples, including endothelial cell proliferation and cisplatin tissue deposition as a surrogate of chemotherapy delivery.

RESULTS

Tumor cell poor, stroma-rich regions exhibited high CA accumulation both in human (meanHUr 0.64 vs. 0.34, p < 0.001) and mouse PDAC (meanAUC60r 2.0 vs. 1.1, p < 0.001). Compared to the baseline, in vivo CA accumulation decreased specifically in response to GEM treatment in a subset of human (HUr -18%) and mouse (AUC60r -36%) tumors. Ex vivo analyses of mPDAC showed reduced cisplatin delivery (GEM: 0.92 ± 0.5 mg/g, vs. vehicle: 3.1 ± 1.5 mg/g, p = 0.004) and diminished endothelial cell proliferation (GEM: 22.3% vs. vehicle: 30.9%, p = 0.002) upon GEM administration.

CONCLUSION

In PDAC, CA accumulation, which is related to tumor vascularization and perfusion, inversely correlates with tumor cellularity. The standard of care GEM treatment results in decreased CA accumulation, which impedes drug delivery. Further investigation is warranted into potentially detrimental effects of GEM in combinatorial therapy regimens.

Vessel network extraction and analysis of mouse pulmonary vasculature via X-ray micro-computed tomographic imaging

In this work, non-invasive high-spatial resolution three-dimensional (3D) X-ray micro-computed tomography (μCT) of healthy mouse lung vasculature is performed. Methodologies are presented for filtering, segmenting, and skeletonizing the collected 3D images. Novel methods for the removal of spurious branch artefacts from the skeletonized 3D image are introduced, and these novel methods involve a combination of distance transform gradients, diameter-length ratios, and the fast marching method (FMM). These new techniques of spurious branch removal result in the consistent removal of spurious branches without compromising the connectivity of the pulmonary circuit. Analysis of the filtered, skeletonized, and segmented 3D images is performed using a newly developed Vessel Network Extraction algorithm to fully characterize the morphology of the mouse pulmonary circuit. The removal of spurious branches from the skeletonized image results in an accurate representation of the pulmonary circuit with significantly less variability in vessel diameter and vessel length in each generation. The branching morphology of a full pulmonary circuit is characterized by the mean diameter per generation and number of vessels per generation. The methods presented in this paper lead to a significant improvement in the characterization of 3D vasculature imaging, allow for automatic separation of arteries and veins, and for the characterization of generations containing capillaries and intrapulmonary arteriovenous anastomoses (IPAVA).

Functional Parameters Derived from Magnetic Resonance Imaging Reflect Vascular Morphology in Preclinical Tumors and in Human Liver Metastases

Purpose: Tumor vessels influence the growth and response of tumors to therapy. Imaging vascular changes in vivo using dynamic contrast-enhanced MRI (DCE-MRI) has shown potential to guide clinical decision making for treatment. However, quantitative MR imaging biomarkers of vascular function have not been widely adopted, partly because their relationship to structural changes in vessels remains unclear. We aimed to elucidate the relationships between vessel function and morphology in vivo Experimental Design: Untreated preclinical tumors with different levels of vascularization were imaged sequentially using DCE-MRI and CT. Relationships between functional parameters from MR (iAUC, K trans, and BATfrac) and structural parameters from CT (vessel volume, radius, and tortuosity) were assessed using linear models. Tumors treated with anti-VEGFR2 antibody were then imaged to determine whether antiangiogenic therapy altered these relationships. Finally, functional-structural relationships were measured in 10 patients with liver metastases from colorectal cancer.Results: Functional parameters iAUC and K trans primarily reflected vessel volume in untreated preclinical tumors. The relationships varied spatially and with tumor vascularity, and were altered by antiangiogenic treatment. In human liver metastases, all three structural parameters were linearly correlated with iAUC and K trans For iAUC, structural parameters also modified each other's effect.Conclusions: Our findings suggest that MR imaging biomarkers of vascular function are linked to structural changes in tumor vessels and that antiangiogenic therapy can affect this link. Our work also demonstrates the feasibility of three-dimensional functional-structural validation of MR biomarkers in vivo to improve their biological interpretation and clinical utility. Clin Cancer Res; 24(19); 4694-704. ©2018 AACR.

Estimating oxygen distribution from vasculature in three-dimensional tumour tissue

Regions of tissue which are well oxygenated respond better to radiotherapy than hypoxic regions by up to a factor of three. If these volumes could be accurately estimated, then it might be possible to selectively boost dose to radio-resistant regions, a concept known as dose-painting. While imaging modalities such as 18F-fluoromisonidazole positron emission tomography (PET) allow identification of hypoxic regions, they are intrinsically limited by the physics of such systems to the millimetre domain, whereas tumour oxygenation is known to vary over a micrometre scale. Mathematical modelling of microscopic tumour oxygen distribution therefore has the potential to complement and enhance macroscopic information derived from PET. In this work, we develop a general method of estimating oxygen distribution in three dimensions from a source vessel map. The method is applied analytically to line sources and quasi-linear idealized line source maps, and also applied to full three-dimensional vessel distributions through a kernel method and compared with oxygen distribution in tumour sections. The model outlined is flexible and stable, and can readily be applied to estimating likely microscopic oxygen distribution from any source geometry. We also investigate the problem of reconstructing three-dimensional oxygen maps from histological and confocal two-dimensional sections, concluding that two-dimensional histological sections are generally inadequate representations of the three-dimensional oxygen distribution.