Quantitative Bioluminescence Tomography-guided System for Conformal Irradiation In Vivo

Although cone-beam CT (CBCT) has been used to guide irradiation for pre-clinical radiotherapy(RT) research, it is limited to localize soft tissue target especially in a low imaging contrast environment. Knowledge of target shape is a fundamental need for RT. Without such information to guide radiation, normal tissue can be irradiated unnecessarily, leading to experimental uncertainties. Recognition of this need led us to develop quantitative bioluminescence tomography (QBLT), which provides strong imaging contrast to localize optical targets. We demonstrated its capability of guiding conformal RT using an orthotopic bioluminescent glioblastoma (GBM) model. With multi-projection and multi-spectral bioluminescence imaging and a novel spectral derivative method, our QBLT system is able to reconstruct GBM with localization accuracy <1mm. An optimal threshold was determined to delineate QBLT reconstructed gross target volume (GTV_{QBLT}), which provides the best overlap between the GTV_{QBLT} and CBCT contrast labeled GBM (GTV), used as the ground truth for the GBM volume. To account for the uncertainty of QBLT in target localization and volume delineation, we also innovated a margin design; a 0.5mm margin was determined and added to GTV_{QBLT} to form a planning target volume (PTV_{QBLT}), which largely improved tumor coverage from 75% (0mm margin) to 98% and the corresponding variation (n=10) of the tumor coverage was significantly reduced. Moreover, with prescribed dose 5Gy covering 95% of PTV_{QBLT}, QBLT-guided 7-field conformal RT can irradiate 99.4 \pm 1.0% of GTV vs. 65.5 \pm 18.5% with conventional single field irradiation (n=10). Our QBLT-guided system provides a unique opportunity for researchers to guide irradiation for soft tissue targets and increase rigorous and reproducibility of scientific discovery

Bioluminescence tomography with Gaussian prior

Parameterizing the bioluminescent source globally in Gaussians provides several advantages over voxel representation in bioluminescence tomography. It is mathematically unique to recover Gaussians [Med. Phys. 31(8), 2289 (2004)] and practically sufficient to approximate various shapes by Gaussians in diffusive medium. The computational burden is significantly reduced since much fewer unknowns are required. Besides, there are physiological evidences that the source can be modeled by Gaussians. The simulations show that the proposed model and algorithm significantly improves accuracy and stability in the presence of Gaussian or non- Gaussian sources, noisy data or the optical background mismatch. It is also validated through in vivo experimental data

Cerenkov Luminescence Tomography for In Vivo Radiopharmaceutical Imaging

Cerenkov luminescence imaging (CLI) is a cost-effective molecular imaging tool for biomedical applications of radiotracers. The introduction of Cerenkov luminescence tomography (CLT) relative to planar CLI can be compared to the development of X-ray CT based on radiography. With CLT, quantitative and localized analysis of a radiopharmaceutical distribution becomes feasible. In this contribution, a feasibility study of in vivo radiopharmaceutical imaging in heterogeneous medium is presented. Coupled with a multimodal in vivo imaging system, this CLT reconstruction method allows precise anatomical registration of the positron probe in heterogeneous tissues and facilitates the more widespread application of radiotracers. Source distribution inside the small animal is obtained from CLT reconstruction. The experimental results demonstrated that CLT can be employed as an available in vivo tomographic imaging of charged particle emitters in a heterogeneous medium

A Born-type approximation method for bioluminescence tomography

In this paper, we present a Born-type approximation method for bioluminescence tomography (BLT), which is to reconstruct an internal bioluminescent source from the measured bioluminescent signal on the external surface of a small animal. Based on the diffusion approximation for the photon propagation in biological tissue, this BLT method utilizes the Green function to establish a linear relationship between the measured bioluminescent signal and the internal bioluminescent source distribution. The Green function can be modified to describe a heterogeneous medium with an arbitrary boundary using the Born approximation. The BLT reconstruction is formulated in a linear least-squares optimization framework with simple bounds constraint. The performance of this method is evaluated in numerical simulation and phantom experiments

Modeling pre-metastatic lymphvascular niche in the mouse ear sponge assay.

Lymphangiogenesis, the formation of new lymphatic vessels, occurs in primary tumors and in draining lymph nodes leading to pre-metastatic niche formation. Reliable in vivo models are becoming instrumental for investigating alterations occurring in lymph nodes before tumor cell arrival. In this study, we demonstrate that B16F10 melanoma cell encapsulation in a biomaterial, and implantation in the mouse ear, prevents their rapid lymphatic spread observed when cells are directly injected in the ear. Vascular remodeling in lymph nodes was detected two weeks after sponge implantation, while their colonization by tumor cells occurred two weeks later. In this model, a huge lymphangiogenic response was induced in primary tumors and in pre-metastatic and metastatic lymph nodes. In control lymph nodes, lymphatic vessels were confined to the cortex. In contrast, an enlargement and expansion of lymphatic vessels towards paracortical and medullar areas occurred in pre-metastatic lymph nodes. We designed an original computerized-assisted quantification method to examine the lymphatic vessel structure and the spatial distribution. This new reliable and accurate model is suitable for in vivo studies of lymphangiogenesis, holds promise for unraveling the mechanisms underlying lymphatic metastases and pre-metastatic niche formation in lymph nodes, and will provide new tools for drug testing

Luciferase expression allows bioluminescence imaging but imposes limitations on the orthotopic mouse (4T1) model of breast cancer

Funding Information: Experiments on the 4T1 and 4Tluc2D6 mouse models of breast cancer were supported by the Russian Scientific Foundation, grant 14-14-00882. MRI measurements were carried out on ClinScan 7T located at Center for Collective Usage (CKP) “Medical nanobiotechologies”, located in Russian National Research Medical University. Experiments on the optimization of protocols for DNA immunization were supported by the Russian Scientific Foundation grant 15-15-30039. Optimization of tumor challenge after DNA immunization was supported by the Russian Fund for Basic Research grant 17-04-00583. Participants were trained in the immunization and tumor challenge experiments in the frame of the European Union Twinning project VACTRAIN, grant agreement #692293, and Swedish Institute PI project 19806/2016. Maria Isaguliants and Stefan Petkov were supported by VACTRAIN, and Maria Isaguliants, also by BALTINFECT, grant agreement #316275. Athina Kilpeläinen was supported by the individual study grant of the Swedish Institute #19061/2014. Patrik Hort is gratefully acknowledged for the language editing. Natalia Belikova is gratefully acknowledged for help with the quantification of protein expression based on the exponential calibration curves. Publisher Copyright: © 2017 Nature Publishing Group. All rights reserved.Implantation of reporter-labeled tumor cells in an immunocompetent host involves a risk of their immune elimination. We have studied this effect in a mouse model of breast cancer after the orthotopic implantation of mammary gland adenocarcinoma 4T1 cells genetically labelled with luciferase (Luc). Mice were implanted with 4T1 cells and two derivative Luc-expressing clones 4T1luc2 and 4T1luc2D6 exhibiting equal in vitro growth rates. In vivo, the daughter 4T1luc2 clone exhibited nearly the same, and 4T1luc2D6, a lower growth rate than the parental cells. The metastatic potential of 4T1 variants was assessed by magnetic resonance, bioluminescent imaging, micro-computed tomography, and densitometry which detected 100-μm metastases in multiple organs and bones at the early stage of their development. After 3-4 weeks, 4T1 generated 11.4 ? 2.1, 4T1luc2D6, 4.5 ? 0.6; and 4T1luc2, 〈1 metastases per mouse, locations restricted to lungs and regional lymph nodes. Mice bearing Luc-expressing tumors developed IFN-? Response to the dominant CTL epitope of Luc. Induced by intradermal DNA-immunization, such response protected mice from the establishment of 4T1luc2-tumors. Our data show that natural or induced cellular response against the reporter restricts growth and metastatic activity of the reporter-labelled tumor cells. Such cells represent a powerful instrument for improving immunization technique for cancer vaccine applications.publishersversionPeer reviewe

Early diagnosis of bladder cancer by photoacoustic imaging of tumor-targeted gold nanorods

Detection and removal of bladder cancer lesions at an early stage is crucial for preventing tumor relapse and progression. This study aimed to develop a new technological platform for the visualization of small and flat urothelial lesions of high-grade bladder carcinoma in situ (CIS). We found that the integrin alpha 581, overexpressed in bladder cancer cell lines, murine orthotopic bladder cancer and human bladder CIS, can be exploited as a receptor for targeted delivery of GNRs functionalized with the cyclic CphgisoDGRG peptide (Iso4). The GNRs@Chit-Iso4 was stable in urine and selectively recognized alpha 581 positive neoplastic urothelium, while low frequency ultrasound-assisted shaking of intravesically instilled GNRs@Chit-Iso4 allowed the distribution of nanoparticles across the entire volume of the bladder. Photoacoustic imaging of GNRs@Chit-Iso4 bound to tumor cells allowed for the detection of neoplastic lesions smaller than 0.5 mm that were undetectable by ultrasound imaging and bioluminescence

Improvements in the robustness and accuracy of bioluminescence tomographic reconstructions of distributed sources within small animals

High quality three-dimensional bioluminescence tomographic (BLT) images, if available, would constitute a major advance and provide much more useful information than the two-dimensional bioluminescence images that are frequently used today. To-date, high quality BLT images have not been available, largely because of the poor quality of the data being input into the reconstruction process. Many significant confounds are not routinely corrected for and the noise in this data is unnecessarily large and poorly distributed. Moreover, many of the design choices affecting image quality are not well considered, including choices regarding the number and type of filters used when making multispectral measurements and choices regarding the frequency and uniformity of the sampling of both the range and domain of the BLT inverse problem. Finally, progress in BLT image quality is difficult to gauge owing to a lack of realistic gold-standard references that engage the full complexity and uncertainty within a small animal BLT imaging experiment. Within this dissertation, I address all of these issues. I develop a Cerenkov-based gold-standard wherein a Positron Emission Tomography (PET) image can be used to gauge improvements in the accuracy of BLT reconstruction algorithms. In the process of creating this reference, I discover and describe corrections for several confounds that if left uncorrected would introduce artifacts into the BLT images. This includes corrections for the angle of the animal’s skin surface relative to the camera, for the height of each point on the skin surface relative to the focal plane, and for the variation in bioluminescence intensity as a function of luciferin concentration over time. Once applied, I go on to derive equations and algorithms that when employed are able to minimize the noise in the final images under the constraints of a multispectral BLT data acquisition. These equations and algorithms allow for an optimal choice of filters to be made and for the acquisition time to be optimally distributed among those filtered measurements. These optimizations make use of Barrett’s and Moore-Penrose pseudoinverse matrices which also come into play in a paradigm I describe that can be used to guide choices regarding sampling of the domain and range