Open Access
9 July 2022 In vivo accurate detection of the liver tumor with pharmacokinetic parametric images from dynamic fluorescence molecular tomography
Fei Liu, Peng Zhang, Zeyu Liu, Fan Song, Chenbin Ma, Yangyang Sun, Youdan Feng, Yufang He, Guanglei Zhang
Author Affiliations +
Abstract

Significance: Pharmacokinetic parametric images in dynamic fluorescence molecular tomography (FMT) can describe three-dimensional (3D) physiological and pathological information inside biological tissues, potentially providing quantitative assessment tools for biological research and drug development.

Aim: In vivo imaging of the liver tumor with pharmacokinetic parametric images from dynamic FMT based on the differences in metabolic properties of indocyanine green (ICG) between normal liver cells and tumor liver cells inside biological tissues.

Approach: First, an orthotopic liver tumor mouse model was constructed. Then, with the help of the FMT/computer tomography (CT) dual-modality imaging system and the direct reconstruction algorithm, 3D imaging of liver metabolic parameters in nude mice was achieved to distinguish liver tumors from normal tissues. Finally, pharmacokinetic parametric imaging results were validated against in vitro anatomical results.

Results: This letter demonstrates the ability of dynamic FMT to monitor the pharmacokinetic delivery of the fluorescent dye ICG in vivo, thus, enabling the distinction between normal and tumor tissues based on the pharmacokinetic parametric images derived from dynamic FMT.

Conclusions: Compared with CT structural imaging technology, dynamic FMT combined with compartmental modeling as an analytical method can obtain quantitative images of pharmacokinetic parameters, thus providing a more powerful research tool for organ function assessment, disease diagnosis and new drug development.

1.

Introduction

Dynamic fluorescence molecular tomography (FMT) is a powerful imaging tool that can resolve the three-dimensional (3D) distribution of fluorophores in small animals and has been extensively used for imaging studies of pharmacokinetics parameters.1 Compartment modeling is a classic approach for pharmacokinetic analysis, and the exchange rates of materials between different compartments are defined as pharmacokinetic parameters.2 The combination of dynamic FMT and compartmental modeling can provide 3D images of pharmacokinetic parameters in vivo for tumor detection, drug development, and metabolism studies.3 However, dynamic FMT remains challenging in imaging pharmacokinetic parameters of small animal organs due to the severe ill-posed nature of the inverse problem, making it challenging to determine the location of specific organs and perform the pharmacokinetics analysis.4 Furthermore, when different tissues or organs are assigned with the same average optical parameters, inaccurate FMT forward models can also be generated, further reducing the accuracy of the pharmacokinetic parameter analysis.

Previous studies have shown that the prior information of anatomical structures provided by structural imaging techniques such as computer tomography (CT) or magnetic resonance imaging (MRI) can be used to construct more accurate forward physical models by assigning heterogeneous optical parameters to different tissues or organs after segmenting the structural information; Moreover, the structural information can be used to construct a priori models, which can guide and constrain the solution process of the inverse problem and improve the reconstruction quality.57 Therefore, a hybrid dual-modality FMT/CT system has been employed in this work to image the pharmacokinetic parameters.

The liver is an essential organ for studying metabolism parameters in dynamic FMT.8 Indocyanine green (ICG) is a clinically approved reagent for liver function assessment by the Food and Drug Administration. Meanwhile, ICG is a widely used optical fluorophore in the near-infrared range for FMT. Therefore, liver tumors can be distinguished in vivo according to the different metabolic properties of ICG in normal and liver cancer cells by pharmacokinetic parameters imaging with the help of dynamic FMT.

In this letter, dynamic FMT is conducted to detect liver tumors in vivo. First, the orthotopic liver tumor mouse model was constructed. Then, 3D imaging of ICG metabolic parameters of the mouse liver was achieved with the help of the FMT/CT dual-modality imaging system and a direct reconstruction algorithm, and liver tumors can be discriminated based on the reconstructed parametric images according to the different metabolic properties of ICG in liver cancer cells. Finally, the parametric imaging results were validated against in vitro anatomical results.

2.

Materials and Methods

2.1.

FMT/CT Dual-Modality Imaging System

The dual-modality system consists of two imaging modalities, FMT and CT, which are placed vertically at a 90-deg angle.9 The schematic diagram of the system is shown in Fig. 1. A laser (i) (Asahi Spectra, Torrance, California) is employed to excite the imaged small animal (viii). Full-angle measurements are implemented by a rotation stage working under a step-by-step mode. Data acquisition is performed by an electron-multiplying charge-coupled device camera (iii) (iXon DU-897, Andor Technologies, Belfast, Northern Ireland) coupled with a 35-mm f/1.6 lens (iv) (C3514-M, Pentax, Japan). Corresponding excitation filter (ii) and emission filter (v) pairs are employed for data acquisition. The x-ray source (vi) is an XTF5011 X-ray tube (Oxford Instrument, California), and the x-ray detector (vii) is a 1056×1056  pixel CMOS array detector (Hamamatsu, Shizuoka, Japan) for acquiring CT projection data.

Fig. 1

Schematic diagram of the FMT/CT dual-modality imaging system.

JBO_27_7_070501_f001.png

In the dual-modality system, the laser was used as the excitation light source, and the shape of the excitation light source was a linear excitation light with a length of 2 cm and a width of 2 mm. The field of view of the charged-coupled device (CCD) camera was 8.2  cm×8.2  cm. The CCD pixels used for fluorescence projection data acquisition were 2×2, the exposure time was set as 0.5 s, and the excitation light power was 1.5 mW.

2.2.

Experiments

2.2.1.

Cell culture

HepG2 cells of a human hepatocellular carcinoma cell line were cultured in Dulbecco’s modified Eagle’s Medium (Invitrogen, Carlsbad, California) with 10% fetal bovine serum (Gibco) in a humidified incubator at 37°C in a 5% atmosphere. Cells were harvested by trypsinization and resuspended in phosphate buffer solution (PBS) before injection.

2.2.2.

Animal studies

Animal experiments were conducted under the protocol approved by the Institutional Animal Care and Use Committee of Beihang University. An 8-week BALB/c mice (female) mouse was anesthetized by intraperitoneal injection of 2% sodium pentobarbital solution at a dose of 0.3 ml/100 g body weight. The mouse was placed in a supine position and fixed on a mouse plate. The abdominal skin was wiped with iodine and alcohol for routine disinfection. An incision was made with surgical scissors along the left upper abdominal pararectal line and peritoneum to expose the liver, and 0.02 ml of 4×105 HepG2 cells in PBS was injected into the right lobe of the liver. After inoculation, the abdominal incision was closed with surgical sutures, and gentamycin eye ointment was applied to the suture to prevent infection of the incision.

2.2.3.

Data acquisition

The data acquisition process for FMT/CT dual-modality imaging system is as follows. During the acquisition of fluorescence projection data for dynamic fluorescence molecular imaging, the imaged mouse was rotated continuously for multiple frames on the rotating platform, and a total of 110 frames were collected in this study. The fluorescence data from each frame was used for FMT reconstruction and to obtain the 3D image. After the fluorescence projection data acquisition was completed, the data acquisition of the CT system was performed to obtain the anatomical structure of the imaged mouse. For the acquisition of CT projection data, only one rotation of the imaged mouse was required.

Additionally, to achieve accurate registration of the imaging results of the FMT and CT systems, the metal marker dot was attached to the body surface of the nude mouse, and a white light image and a CT projection image were acquired separately to provide the height information required for the registration of the two imaging systems.

2.3.

Pharmacokinetic Parameter Reconstruction Method

The compartment modeling can be expressed by the differential equation as follows:

Eq. (1)

{dc1(t)dt=f10+i=2n[k1ici(t)ki1ci(t)]k01c1(t)dci(t)dt=ki1ci(t)k1ici(t),
where c1 denotes the concentration of the fluorescent probe in the blood, ci (i=2n) denotes the concentration of the fluorescent probe in different tissues and organs, ki1 denotes the rate coefficient of fluorescent probe inflow from blood to organ i (i=2n), k1i denotes the rate coefficient of fluorescent probe inflow into the blood from different organs i (i=2n), k01 denotes the rate coefficient of the fluorescent probe excreted from the blood, and f10 is the rate at which the fluorescent probe is injected into the blood.

By solving Eq. (1), the metabolism curves of ICG in different tissues and organs can be obtained:

Eq. (2)

Ci(t)=Aieαit+Bieβit,
where Ci(t) is the ICG concentration at the time of t; A, B, α, and β are the metabolic parameters that can characterize the pharmacokinetic properties of each tissue and organ. α denotes the uptake rate of the fluorescent probe by the tissue, β denotes the excretion rate of the fluorescent probe by the tissue, A and B denote the concentration gain of the model, which is related to the concentration of the fluorescent probe.

To obtain the metabolic parameters of dynamic FMT reconstruction, the internal anatomy of nude mice was first obtained by CT as a priori information for the dynamic FMT reconstruction algorithm. Then, the dynamic FMT imaging region of the mouse was divided into four parts: lung, liver, kidney, and other tissues based on the structural prior information provided by the CT results, and finally, the corresponding optical parameters were assigned to the different organs of the mice based on Table 1 to construct the heterogeneous forward model.

Table 1

Heterogeneous optical parameters in mice.10

Organsμa (cm−1)μs′ (cm−1)
Lung0.2530
Liver0.513
Kidney0.17520.9
Other tissues0.310

Then, the double exponential curve model obtained from Eq. (2) is substituted into the heterogeneous forward model, and the results of the forward prediction function are obtained as follows:

Eq. (3)

f(X)=j=1NWx(tk)=j=1NW[Aexp(αtk)+Bexp(βtk)].
Finally, the inverse problem of Eq. (3) is solved by the full-direct FMT method9 to reconstruct the parametric images. This method can make full use of the temporal correlations of boundary measurements to model the high temporal variation of fluorophore concentration. Exactly, the objective function of the direct reconstruction algorithm can be obtained by substituting Eq. (3) into the Laplace regularization matrix:

Eq. (4)

Ψ(X)=yf(X)22+λ1LA22+λ2LB22+λ3Lα22+λ4Lβ22,=yf(X)22+i=14λiLxi22,
where X denotes the fluorophore concentration distribution, and L denotes the Laplace matrix, and A, B, α, and β are the pharmacokinetic parametric images to be reconstructed. The pharmacokinetic parameters A and B have an arbitrary unit (a.u.), while α (min1) and β (min1) are the uptake and excretion rates of fluorophores, which have physiological significance for quantitative evaluation of organ function.11 In this letter, the regularization parameters were set as [λ1,λ2,λ3,λ4]=[2×106,5×106,5×103,1×104], and the number of iterations was 200. The CT data were reconstructed using the Feldkamp–Davis–Kress algorithm,12 and the size of the reconstructed voxels was 0.1×0.1×0.1  mm3.

3.

Results and Discussion

The anatomical imaging results of the mouse are shown in Fig. 2. Figure 2(a) shows a representative CT projection image of the tumor-bearing mouse. Figure 2(b) shows three crosssectional images of the CT reconstruction results of the tumor-bearing mouse with the heights corresponding to the positions of the red, green, and blue dashed lines in Fig. 2(a), respectively. The whole imaging area of the mouse was manually segmented into lung, liver, kidney, and other tissues [Fig. 2(c)] according to the gray values of the CT images, and these segmentation results will be used as structural a priori information for the metabolic parameters of the FMT reconstruction.

Fig. 2

CT imaging results of the tumor-bearing nude mouse. (a) A representative CT projection image of the tumor-bearing mouse. (b) Three cross-sectional images of the CT results. (c) CT segmentation results of lung, liver, kidney, and other tissues.

JBO_27_7_070501_f002.png

Figure 3 shows the dynamic FMT reconstruction results of liver metabolic parameters. Figure 3(a) shows a CT projection, and the green dotted lines indicate the corresponding heights of the reconstructed images in Figs. 3(b)3(d). Figure 3(b) shows an image of a frozen mouse section obtained after imaging, the height of which was approximately matched to the height of the tomographic results by manual marking. The actual location of the tumor in the liver is shown by the green arrow in Fig. 3(b). The liver tumor region can be clearly identified from the imaging results of metabolic parameters shown in Fig. 3(d). Moreover, based on the reconstructed metabolic parameter images, four metabolic parameter images A, B, α, and β can be obtained. By substituting the parameter values into Eq. (2), we can quantitatively recover the ICG metabolic curve of the liver. As shown in Fig. 3(e), the blue curve represents the ICG metabolic curve at the location of the blue point in the normal liver region in Fig. 3(b), and the red curve represents the ICG metabolic curve at the location of the red point in the tumor region in Fig. 3(b). It can be found that compared with the normal liver tissue, the uptake rate of ICG is increased and the excretion rate of ICG is decreased in the tumor area, which is consistent with the theoretical expectation.

Fig. 3

Imaging results of metabolic parameters in liver tumors. (a) CT projection image. (b) Cryosection image of mouse with a green arrow indicating the tumor location. (c) Cross-sectional view of CT imaging results. (d) The metabolic parameter images A (a.u.), B (a.u.), α (min1), and β (min1) of the mouse. (b)–(d) correspond to the heights represented by the green-dashed lines in (a). (e) ICG metabolic curves of the liver.

JBO_27_7_070501_f003.png

4.

Conclusion

Pharmacokinetic parameters can describe the dynamics of drug absorption, metabolism and excretion processes and enable quantitative monitoring of physiological properties such as permeability and perfusion rates of tissues and organs. Metabolic parameter imaging technology can visualize the pharmacokinetic parameters of each voxel in the imaging area, allowing a more intuitive description of the physiological processes of drug absorption and excretion in different tissues and organs, which is of great significance for physiological and pathological studies.

In this letter, the pharmacokinetic parameters of ICG have been monitored through a hybrid FMT/CT system and the direct reconstruction algorithm. Although liver tumors can also be directly observed from CT imaging results, as a structural imaging technique, the functional state of their tissues and organs cannot be reflected in the CT images. As a functional imaging technique, quantitative assessment of liver function can be obtained by pharmacokinetic parametric images derived from dynamic FMT. Furthermore, ICG metabolic curves at any location in the liver can be achieved according to the reconstructed pharmacokinetic parameters. Compared with CT structural imaging, the imaging of metabolic parameters with dynamic FMT can obtain quantitative metabolic parameter images, providing a more powerful research tool for organ function assessment, disease diagnosis, and new drug development. Future work will focus on imaging pharmacokinetic parameters of pathological livers and other organs.

Disclosures

The authors have no relevant financial interests in this article and no potential conflicts of interest to disclose. This work involved human subjects or animals in its research. Approval of all ethical and experimental procedures and protocols was granted by the Biological and Medical Ethics Committee of Beihang University (Grant No. BM20200162).

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 61871022), the Beijing Natural Science Foundation (Grant No. 7202102), and the Academic Excellence Foundation of BUAA for PhD Students.

References

1. 

G. Zhang et al., “Bayesian framework based direct reconstruction of fluorescence parametric images,” IEEE Trans. Med. Imaging, 34 (6), 1378 –1391 (2015). https://doi.org/10.1109/TMI.2015.2394476 ITMID4 0278-0062 Google Scholar

2. 

B. Alacam and B. Yazici, “Direct reconstruction of pharmacokineticrate images of optical fluorophores from NIR measurements,” IEEE Trans. Med. Imaging, 28 (9), 1337 –1353 (2009). https://doi.org/10.1109/TMI.2009.2015294 ITMID4 0278-0062 Google Scholar

3. 

G. Zhang et al., “A direct method with structural priors for imaging pharmacokinetic parameters in dynamic fluorescence molecular tomography,” IEEE Trans. Biomed. Eng., 61 (3), 986 –990 (2014). https://doi.org/10.1109/TBME.2013.2292714 IEBEAX 0018-9294 Google Scholar

4. 

X. Liu et al., “Unmixing dynamic fluorescence diffuse optical tomography images with independent component analysis,” IEEE Trans. Med. Imaging, 30 (9), 1591 –1604 (2011). https://doi.org/10.1109/TMI.2011.2134865 ITMID4 0278-0062 Google Scholar

5. 

X. Guo et al., “A combined fluorescence and micro-computed tomography system for small animal imaging,” IEEE Trans. Biomed. Eng., 57 (12), 2876 –2883 (2010). https://doi.org/10.1109/TBME.2010.2073468 IEBEAX 0018-9294 Google Scholar

6. 

A. Ale et al., “FMT-XCT: in vivo animal studies with hybrid fluorescence molecular tomography–X-ray computed tomography,” Nat. Methods, 9 (6), 615 –620 (2012). https://doi.org/10.1038/nmeth.2014 1548-7091 Google Scholar

7. 

P. Zhang et al., “A review of advances in imaging methodology in fluorescence molecular tomography,” Phys. Med. Biol., 67 (10), 10TR01 (2022). https://doi.org/10.1088/1361-6560/ac5ce7 PHMBA7 0031-9155 Google Scholar

8. 

G. Zhang et al., “Imaging of pharmacokinetic rates of indocyanine green in mouse liver with a hybrid fluorescence molecular tomography/x-ray computed tomography system,” J. Biomed. Opt., 18 (4), 040505 (2013). https://doi.org/10.1117/1.JBO.18.4.040505 JBOPFO 1083-3668 Google Scholar

9. 

P. Zhang et al., “UHR-DeepFMT: ultra-high spatial resolution reconstruction of fluorescence molecular tomography based on 3-D fusion dual-sampling deep neural network,” IEEE Trans. Med. Imaging, 40 (11), 3217 –3228 (2021). https://doi.org/10.1109/TMI.2021.3071556 ITMID4 0278-0062 Google Scholar

10. 

Y. An et al., “A novel region reconstruction method for fluorescence molecular tomography,” IEEE Trans. Biomed. Eng., 62 (7), 1818 –1826 (2015). https://doi.org/10.1109/TBME.2015.2404915 IEBEAX 0018-9294 Google Scholar

11. 

G. Zhang et al., “Full-direct method for imaging pharmacokinetic parameters in dynamic fluorescence molecular tomography,” Appl. Phys. Lett., 106 (8), 081110 (2015). https://doi.org/10.1063/1.4913690 APPLAB 0003-6951 Google Scholar

12. 

T. Rodet et al., “The cone‐beam algorithm of Feldkamp, Davis, and Kress preserves oblique line integrals,” Med. Phys., 31 (7), 1972 –1975 (2004). https://doi.org/10.1118/1.1759828 Google Scholar
CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Fei Liu, Peng Zhang, Zeyu Liu, Fan Song, Chenbin Ma, Yangyang Sun, Youdan Feng, Yufang He, and Guanglei Zhang "In vivo accurate detection of the liver tumor with pharmacokinetic parametric images from dynamic fluorescence molecular tomography," Journal of Biomedical Optics 27(7), 070501 (9 July 2022). https://doi.org/10.1117/1.JBO.27.7.070501
Received: 23 May 2022; Accepted: 28 June 2022; Published: 9 July 2022
Lens.org Logo
CITATIONS
Cited by 3 scholarly publications.
Advertisement
Advertisement
KEYWORDS
Liver

Tumors

Computed tomography

Tissues

Imaging systems

Luminescence

Tomography

Back to Top