Comparison of permissible source region and multispectral data using efficient bioluminescence tomography method

ABPO-TVSCAD: alternating Bregman proximity operators approach based on TVSCAD regularization for bioluminescence tomography

Objective.Bioluminescence tomography (BLT) is a promising non-invasive optical medical imaging technique, which can visualize and quantitatively analyze the distribution of tumor cells in living tissues. However, due to the influence of photon scattering effect and ill-conditioned inverse problem, the reconstruction result is unsatisfactory. The purpose of this study is to improve the reconstruction performance of BLT.Approach.An alternating Bregman proximity operators (ABPO) method based on TVSCAD regularization is proposed for BLT reconstruction. TVSCAD combines the anisotropic total variation (TV) regularization constraints and the non-convex smoothly clipped absolute deviation (SCAD) penalty constraints, to make a trade-off between the sparsity and edge preservation of the source. ABPO approach is used to solve the TVSCAD model (ABPO-TVSCAD for short). In addition, to accelerate the convergence speed of the ABPO, we adapt the strategy of shrinking the permission source region, which further improves the performance of ABPO-TVSCAD.Main results.The results of numerical simulations andin vivoxenograft mouse experiment show that our proposed method achieved superior accuracy in spatial localization and morphological reconstruction of bioluminescent source.Significance.ABPO-TVSCAD is an effective and robust reconstruction method for BLT, and we hope that this method can promote the development of optical molecular tomography.

Multispectral Differential Reconstruction Strategy for Bioluminescence Tomography

Bioluminescence tomography (BLT) is a promising in vivo molecular imaging tool that allows non-invasive monitoring of physiological and pathological processes at the cellular and molecular levels. However, the accuracy of the BLT reconstruction is significantly affected by the forward modeling errors in the simplified photon propagation model, the measurement noise in data acquisition, and the inherent ill-posedness of the inverse problem. In this paper, we present a new multispectral differential strategy (MDS) on the basis of analyzing the errors generated from the simplification from radiative transfer equation (RTE) to diffusion approximation and data acquisition of the imaging system. Through rigorous theoretical analysis, we learn that spectral differential not only can eliminate the errors caused by the approximation of RTE and imaging system measurement noise but also can further increase the constraint condition and decrease the condition number of system matrix for reconstruction compared with traditional multispectral (TM) reconstruction strategy. In forward simulations, energy differences and cosine similarity of the measured surface light energy calculated by Monte Carlo (MC) and diffusion equation (DE) showed that MDS can reduce the systematic errors in the process of light transmission. In addition, in inverse simulations and in vivo experiments, the results demonstrated that MDS was able to alleviate the ill-posedness of the inverse problem of BLT. Thus, the MDS method had superior location accuracy, morphology recovery capability, and image contrast capability in the source reconstruction as compared with the TM method and spectral derivative (SD) method. In vivo experiments verified the practicability and effectiveness of the proposed method.

A review of the application of machine learning in molecular imaging

Molecular imaging (MI) is a science that uses imaging methods to reflect the changes of molecular level in living state and conduct qualitative and quantitative studies on its biological behaviors in imaging. Optical molecular imaging (OMI) and nuclear medical imaging are two key research fields of MI. OMI technology refers to the optical information generated by the imaging target (such as tumors) due to drug intervention and other reasons. By collecting the optical information, researchers can track the motion trajectory of the imaging target at the molecular level. Owing to its high specificity and sensitivity, OMI has been widely used in preclinical research and clinical surgery. Nuclear medical imaging mainly detects ionizing radiation emitted by radioactive substances. It can provide molecular information for early diagnosis, effective treatment and basic research of diseases, which has become one of the frontiers and hot topics in the field of medicine in the world today. Both OMI and nuclear medical imaging technology require a lot of data processing and analysis. In recent years, artificial intelligence technology, especially neural network-based machine learning (ML) technology, has been widely used in MI because of its powerful data processing capability. It provides a feasible strategy to deal with large and complex data for the requirement of MI. In this review, we will focus on the applications of ML methods in OMI and nuclear medical imaging.

Development of tomographic reconstruction for three-dimensional optical imaging: From the inversion of light propagation to artificial intelligence

Optical molecular tomography (OMT) is an imaging modality which uses an optical signal, especially near-infrared light, to reconstruct the three-dimensional information of the light source in biological tissue. With the advantages of being low-cost, noninvasive and having high sensitivity, OMT has been applied in preclinical and clinical research. However, due to its serious ill-posedness and ill-condition, the solution of OMT requires heavy data analysis and the reconstruction quality is limited. Recently, the artificial intelligence (commonly known as AI)-based methods have been proposed to provide a different tool to solve the OMT problem. In this paper, we review the progress on OMT algorithms, from conventional methods to AI-based methods, and we also give a prospective towards future developments in this domain.

Improved Block Sparse Bayesian Learning Method Using K-Nearest Neighbor Strategy for Accurate Tumor Morphology Reconstruction in Bioluminescence Tomography

OBJECTIVE

Bioluminescence tomography (BLT) is a non-invasive technique designed to enable three-dimensional (3D) visualization and quantification of viable tumor cells in living organisms. However, despite the excellent sensitivity and specificity of bioluminescence imaging (BLI), BLT is limited by the photon scattering effect and ill-posed inverse problem. If the complete structural information of a light source is considered when solving the inverse problem, reconstruction accuracy will be improved.

METHODS

This article proposed a block sparse Bayesian learning method based on K-nearest neighbor strategy (KNN-BSBL), which incorporated several types of a priori information including sparsity, spatial correlations among neighboring points, and anatomical information to balance over-sparsity and morphology preservation in BLT. Furthermore, we considered the Gaussian weighted distance prior in a light source and proposed a KNN-GBSBL method to further improve the performance of KNN-BSBL.

RESULTS

The results of numerical simulations and in vivo glioma-bearing mouse experiments demonstrated that KNN-BSBL and KNN-GBSBL achieved superior accuracy for tumor spatial positioning and morphology reconstruction.

CONCLUSION

The proposed method KNN-BSBL incorporated several types of a priori information is an efficient and robust reconstruction method for BLT.

Filtered maximum likelihood expectation maximization based global reconstruction for bioluminescence tomography

The reconstruction of bioluminescence tomography (BLT) is severely ill-posed due to the insufficient measurements and diffuses nature of the light propagation. Predefined permissible source region (PSR) combined with regularization terms is one common strategy to reduce such ill-posedness. However, the region of PSR is usually hard to determine and can be easily affected by subjective consciousness. Hence, we theoretically developed a filtered maximum likelihood expectation maximization (fMLEM) method for BLT. Our method can avoid predefining the PSR and provide a robust and accurate result for global reconstruction. In the method, the simplified spherical harmonics approximation (SP_N) was applied to characterize diffuse light propagation in medium, and the statistical estimation-based MLEM algorithm combined with a filter function was used to solve the inverse problem. We systematically demonstrated the performance of our method by the regular geometry- and digital mouse-based simulations and a liver cancer-based in vivo experiment. Graphical abstract The filtered MLEM-based global reconstruction method for BLT.

A hybrid clustering algorithm for multiple-source resolving in bioluminescence tomography

Bioluminescence tomography is a preclinical imaging modality to locate and quantify internal bioluminescent sources from surface measurements, which experienced rapid growth in the last 10 years. However, multiple-source resolving remains a challenging issue in BLT. In this study, it is treated as an unsupervised pattern recognition problem based on the reconstruction result, and a novel hybrid clustering algorithm combining the advantages of affinity propagation (AP) and K-means is developed to identify multiple sources automatically. Moreover, we incorporate the clustering analysis into a general multiple-source reconstruction framework, which can provide stable reconstruction and accurate resolving result without providing the number of targets. Numerical simulations and in vivo experiments on 4T1-luc2 mouse model were conducted to assess the performance of the proposed method in multiple-source resolving. The encouraging results demonstrate significant effectiveness and potential of our method in preclinical BLT applications.

Bioluminescence Tomography Based on Gaussian Weighted Laplace Prior Regularization for In Vivo Morphological Imaging of Glioma

Bioluminescence tomography (BLT) is a powerful non-invasive molecular imaging tool for in vivo studies of glioma in mice. However, because of the light scattering and resulted ill-posed problems, it is challenging to develop a sufficient reconstruction method, which can accurately locate the tumor and define the tumor morphology in three-dimension. In this paper, we proposed a novel Gaussian weighted Laplace prior (GWLP) regularization method. It considered the variance of the bioluminescence energy between any two voxels inside an organ had a non-linear inverse relationship with their Gaussian distance to solve the over-smoothed tumor morphology in BLT reconstruction. We compared the GWLP with conventional Tikhonov and Laplace regularization methods through various numerical simulations and in vivo orthotopic glioma mouse model experiments. The in vivo magnetic resonance imaging and ex vivo green fluorescent protein images and hematoxylin-eosin stained images of whole head cryoslicing specimens were utilized as gold standards. The results demonstrated that GWLP achieved the highest accuracy in tumor localization and tumor morphology preservation. To the best of our knowledge, this is the first study that achieved such accurate BLT morphological reconstruction of orthotopic glioma without using any segmented tumor structure from any other structural imaging modalities as the prior for reconstruction guidance. This enabled BLT more suitable and practical for in vivo imaging of orthotopic glioma mouse models.

Reconstruction algorithm for fluorescence molecular tomography using sorted L-one penalized estimation

Fluorescence molecular tomography (FMT) has been a promising imaging tool that provides convenience for accurate localization and quantitative analysis of the fluorescent probe. In this study, we present a reconstruction method combining sorted L-one penalized estimation with an iterative-shrinking permissible region strategy to reconstruct fluorescence targets. Both numerical simulation experiments on a three-dimensional digital mouse model and physical experiments on a cubic phantom were carried out to validate the accuracy, effectiveness, and robustness of the proposed method. The results indicate that the proposed method can produce better location and satisfactory fluorescent yield with computational efficiency, which makes it a practical and promising reconstruction method for FMT.

A Multi-Camera System for Bioluminescence Tomography in Preclinical Oncology Research

Bioluminescent imaging (BLI) of cells expressing luciferase is a valuable noninvasive technique for investigating molecular events and tumor dynamics in the living animal. Current usage is often limited to planar imaging, but tomographic imaging can enhance the usefulness of this technique in quantitative biomedical studies by allowing accurate determination of tumor size and attribution of the emitted light to a specific organ or tissue. Bioluminescence tomography based on a single camera with source rotation or mirrors to provide additional views has previously been reported. We report here in vivo studies using a novel approach with multiple rotating cameras that, when combined with image reconstruction software, provides the desired representation of point source metastases and other small lesions. Comparison with MRI validated the ability to detect lung tumor colonization in mouse lung.

Hybrid multilevel sparse reconstruction for a whole domain bioluminescence tomography using adaptive finite element

Quantitative reconstruction of bioluminescent sources from boundary measurements is a challenging ill-posed inverse problem owing to the high degree of absorption and scattering of light through tissue. We present a hybrid multilevel reconstruction scheme by combining the ability of sparse regularization with the advantage of adaptive finite element method. In view of the characteristics of different discretization levels, two different inversion algorithms are employed on the initial coarse mesh and the succeeding ones to strike a balance between stability and efficiency. Numerical experiment results with a digital mouse model demonstrate that the proposed scheme can accurately localize and quantify source distribution while maintaining reconstruction stability and computational economy. The effectiveness of this hybrid reconstruction scheme is further confirmed with in vivo experiments.

Three-dimensional noninvasive monitoring iodine-131 uptake in the thyroid using a modified Cerenkov luminescence tomography approach

Background

Cerenkov luminescence tomography (CLT) provides the three-dimensional (3D) radiopharmaceutical biodistribution in small living animals, which is vital to biomedical imaging. However, existing single-spectral and multispectral methods are not very efficient and effective at reconstructing the distribution of the radionuclide tracer. In this paper, we present a semi-quantitative Cerenkov radiation spectral characteristic-based source reconstruction method named the hybrid spectral CLT, to efficiently reconstruct the radionuclide tracer with both encouraging reconstruction results and less acquisition and image reconstruction time.

Methodology/Principal Findings

We constructed the implantation mouse model implanted with a 400 µCi Na¹³¹I radioactive source and the physiological mouse model received an intravenous tail injection of 400 µCi radiopharmaceutical Iodine-131 (I-131) to validate the performance of the hybrid spectral CLT and compared the reconstruction results, acquisition, and image reconstruction time with that of single-spectral and multispectral CLT. Furthermore, we performed 3D noninvasive monitoring of I-131 uptake in the thyroid and quantified I-131 uptake in vivo using hybrid spectral CLT. Results showed that the reconstruction based on the hybrid spectral CLT was more accurate in localization and quantification than using single-spectral CLT, and was more efficient in the in vivo experiment compared with multispectral CLT. Additionally, 3D visualization of longitudinal observations suggested that the reconstructed energy of I-131 uptake in the thyroid increased with acquisition time and there was a robust correlation between the reconstructed energy versus the gamma ray counts of I-131 (). The ex vivo biodistribution experiment further confirmed the I-131 uptake in the thyroid for hybrid spectral CLT.

Conclusions/Significance

Results indicated that hybrid spectral CLT could be potentially used for thyroid imaging to evaluate its function and monitor its treatment for thyroid cancer.