1.IntroductionNonalcoholic fatty liver disease (NAFLD) is characterized by excessive hepatic fat accumulation associated with an increasing rate of obesity, diabetes, and metabolic syndrome. The prevalence of NAFLD is globally, afflicting about 1 billion individuals worldwide, which is becoming the predominant cause of chronic liver disease and a substantial clinical and economic burden around the world.1–3 In the context of NAFLD, the first diagnostic requirement is to show the presence of fat in the liver accurately. The presence of fat to 10% is considered abnormal.1,4 Liver biopsy is the gold standard for the diagnosis and classification of NAFLD. However, it suffers from being invasive, expensive, associated with potential complications and sampling error, and plagued with interobserver variability.5–7 Magnetic resonance spectroscopy or magnetic resonance imaging (MRI) quantifies the proton density fat fraction (PDFF) to measure hepatic triglyceride amount. However, both techniques are costly with restricted accessibility and high requirements for expertise.2,8–10 Compared to MRI-PDFF, the transient elastography-based controlled attenuation parameter is less accurate for the diagnosis of hepatic steatosis.11 Optical method, such as near-infrared hyperspectral imaging, can be an alternative that is a noninvasive, noncontact, nonionizing, and label-free technique. However, it suffers from penetration depth limitation due to the scattering of light.12 Although a variety of techniques are available for hepatic steatosis evaluation, little is known about the liver function reserve (LFR) with the progression of NAFLD.13 LFR is the remnant functional capacity of the liver after liver injury, which is considered an important factor for treatment and prognosis in chronic liver diseases. Indocyanine green (ICG) clearance test is the acknowledged approach for LFR assessment, the nature of which is to accurately detect ICG concentration change over time in blood. ICG can be selectively up taken by liver cells, eliminated to bile and excreted with feces through the bowel in the original form. ICG was not uptake by other organs or tissue and is an ideal drug for LFR reflection. In hepatic surgery, ICG retention rate at 15 min is an essential condition for safe hepatectomy. is the contraindication for hepatectomy because posthepatectomy liver failure can occur with a large probability.14–16 Photoacoustic imaging (PAI) is an emerging noninvasive modality focusing on the measurement of optical absorption of endogenous or exogenous absorbers, which can afford high-resolution images in superficial or deep tissue. Its clinical applicability has been demonstrated in breast cancer, thyroid cancer, melanoma, arthritis, bowel disease,17–21 and in the detection of endogenous molecular chromophores, such as melanin, hemoglobin, lipids, and collagen.22–25 Steatosis in fatty liver mice diagnosed with PAI indicated a prominent difference from normal mice.26–28 However, previous applications were focused on small animals. Here we want to test the feasibility in a larger animal model with similar imaging depth to check its potential for clinical translations. In addition, combining the LFR assessment by PAI with general contrast agents29 and the multispectral method has not been demonstrated on large animal models. Based on this prior study,29 the goal of the current study is to show that PAI is capable of not only assessing LFR but also visualizing lipid content in fatty liver rabbits and to demonstrate PAI as a new noninvasive imaging biomarker for the assessment of NAFLD. 2.Methods2.1.Experimental SetupBoth in vivo and ex vivo PAI were performed in this study. Before photoacoustic (PA) scanning, the anatomical regions of interest are scanned using an ultrasound system with a linear array working at 8.5 MHz, 100 Hz frame rate (iNSIGHT 23R, Saset Healthcare, Inc. Chengdu, China) by a single professional operator. The setup for in vivo PA imaging [Fig. 1(a)] is based on a -switched Nd:YAG-pumped optical parameter oscillator (Surelite, Continuum, California), a 128-element concave ultrasound transducer array (Japan Probe Corporation, Yokohama, Japan), and a 64-channel DAQ system (PXIe5105, National Instrument, USA). The diameter of the array is 100 mm. The pitch size of each element is 2 mm. The laser provides tunable laser pulses from 700 to 960 nm at a repetition rate of up to 20 Hz. The output laser light is coupled into a custom-made fiber bundle. The output shape of the fiber bundle is a rectangular region with a area. The maximum light fluence from the fiber bundle at 800 nm is estimated to be (within the safety limit of the American National Standards Institute regulations). The concave ultrasound transducer array has a center frequency of 5 MHz ( two-way relative bandwidth) with about 2.1 mm in radius one-way Nyquist zone30,31 and a spatial resolution of for cross-sectional imaging. The array is enclosed with a resin shell covered by a membrane ( thickness), and water is filled into the space between the membrane and transducer. In the animal study, the PAI probe was put on the animal body (ear or skin surface), and the fiber bundle output was adjusted so that the illumination area was right at the center of the transducer [about 25 deg illumination angle, 15 to 20 mm height as shown in Figs. 1(b) and 1(c)]. Detected PA signals are amplified through a custom-built low-noise preamplifier (55 dB) and transferred to a 64-channel analog-to-digital system (PXIe5105, National Instrument, Texas) after 2:1 multiplexing with 12-bit digital resolution and sampling rate. Thus the frame rate of the system is 10 Hz. Signals are saved in the onboard computer (PXIe8840, National Instruments, Texas, USA) which also works as a control panel for the PAI system. The setup for the ex vivo PA imaging is adapted from the in vivo system, as illustrated in Fig. 2(a). The output laser light is coupled into a custom-made fiber bundle to illuminate the whole area of a sample from the above, achieving an averaged light fluence of on the sample surface at 800 nm (again within the safety limit of the American National Standards Institute regulations). A sample piece is cut from an intact lobe of a rabbit liver with a thickness of 1 to 2 mm. The same concave ultrasound transducer array is used to perform 2D tomographic imaging of an ex vivo liver. Both the tissue and transducer are put in a small water tank. The liver tissue is fixed in a cylindrical sample holder made from resin and put in the center of the transducer to ensure maximum ultrasound sensitivity and the laser beam from the fiber was 15 mm above the tissue to produce a rounded area of illumination [Fig. 2(b), water tank and sample holder is not shown here]. After amplification by our custom-built preamplifier, the PA signal is recorded with the DAQ system. 2.2.Animal Model PreparationNew Zealand rabbits (, male, weight: 2.5 to 2.7 kg) were used in this study. The NAFLD group () was fed with a standard diet for the first 8 weeks, followed by a high-fat diet (75% standard feed with 10% butter, 2% cholesterol, 5% sugar, and 8% yolk powder) for the 12 weeks thereafter. The control group () was fed with the standard diet (100% standard feed) for 20 weeks.32 These diets were provided by Dashuo Laboratory Animal Co., LTD, Chengdu, China. 2.3.In Vivo LFR AssessmentFor the in vivo LFR assessment, 11 rabbits in two groups underwent PAI scanning of the liver region and ICG injection. All the animals were first anesthetized by isoflurane gas inhalation and were put on a warming pad and fixed to a custom-made animal pad to set in the right lateral position. We first used traditional ultrasonic imaging (US) to find the left lateral lobe of the liver for reference. The PA imaging probe was then placed at the same position as the US probe to acquire PA images, as shown in Fig. 1(c). Posture and imaging probe placements (both US and PA) were standardized and marked to ensure scanning reproducibility. For in vivo imaging, multispectral PA signals were acquired at 760, 800, 850, and 930 nm to obtain oxygenated hemoglobin (), deoxygenated hemoglobin (HbR), and lipid according to previous studies.21 Twenty images were continuously acquired at one wavelength without averaging. After the multispectral imaging, the PAT imaging probe was placed above the right central auricular artery (CAA) as shown in Fig. 1(b) and the ICG was injected through the left auricular vein. Corresponding US images were acquired at the same position. Before the PAI scan, ICG was dissolved in distilled water to obtain a concentration of . PAI scanning at 805 nm was performed for 20 min over the right CAA. The acquisition rate was 1 frame per second (time-averaged for 10 consecutive frames). One minute after a baseline PAI scan, a dose of ICG solution was injected via the left auricular vein along with a 5-ml saline flush for 10 s. Then the PA signal was continuously acquired for additional 19 min. 2.4.Ex Vivo and Histological ExaminationAfter the in vivo imaging examination, representative tissue specimens were taken from previously imaged anatomical regions. After harvesting, tissue specimens from the same lobe were divided into two parts, one for ex vivo PA imaging, and the other one was fixed in a 4% formaldehyde/PBS solution and then embedded in paraffin. For the ex vivo experiments, multispectral PA signals were acquired at 760, 800, 850, and 930 nm. For histological examination, liver tissue sections were stained with haematoxylin–eosin (HE) and Masson. The histologic features were evaluated with the NAFLD activity score (NAS), which was calculated as the unweighted sum of the scores for steatosis (0 to 3), lobular inflammation (0 to 3), and ballooning (0 to 2). Based on the NAS score, NAFLD was classified as “NAFLD but not-NASH” (), “borderline-NASH” ( to 4), and “definite-NASH” ( to 8). The classification process was performed by a pathologist with 8 years of experience in liver histology, who was blinded to the PA results. To estimate the lipid area in H&E staining images, we used a machine learning toolbox (WeKa3)33 in Fiji (National Institutes of Health). One expert labeled H&E staining images at which have a large area of lipid cells to train the network. Then the trained network was applied to classify all slices then the area of lipid cells was segmented by setting a threshold. The results were shown using the lipid cells area divided by the total area of the image. 2.5.Image Reconstruction and Spectral UnmixingPA images were first reconstructed in Labview (National Instruments, USA) in real time using a back-projection algorithm34 with real-time display. The raw PA data were also stored for offline processing by MATLAB (R2016b, Mathworks, Inc., MA, USA). In the offline processing, 0.2 to 10 MHz filtering was applied to remove high-frequency noise, and then PA images were reconstructed offline using back-projection schemes.35 To remove the breathing-related artifact, we built a MATLAB code as described previously.36 Briefly, after we got 20 frames (each frame obtained without averaging) from one wavelength, a correlation coefficient (CC) was calculated among these frames. Frames with movement showing significantly low CC (slight movement: 0.9 to 1, breathing: 0.5 to 0.7) were rejected. In most cases, about 12 to 14 frames were effective and averaged as one wavelength image. We then repeated the same process above for different wavelengths and registered all averaged images to that for the first wavelength for spectral unmixing. A linear unmixing was applied to obtain the signal from three endogenous absorbers following previous studies,21,37 and then values of HbR, , and lipid were obtained. Lipid unmixing was based on all wavelengths (760, 800, 850, and 930 nm), whereas the HbR and signals were calculated from a subrange (760, 800, and 850 nm), which is more accurate in unmixing due to increased water absorptivity at higher wavelengths.38 The total hemoglobin () and lipid values were calculated from (long × width) regions of interest close to the surface of the skin to minimize error due to the linear unmixing.39 The averaged HbT signal, , lipid signal, and were calculated for statistical analysis. To minimize the error from unknown light fluence, we first measured the laser energy below the imaging probe and normalized the energy among different wavelengths. Then we divided our image into two parts, skin and liver based on the clear structural differences in PA images and referred to corresponding US images. A depth compensation method based on recently published research was applied to each wavelength40 for in vivo data. The relevant optical parameters were obtained from two other papers.41,42 For ex vivo data, we measured the laser energy on the tissue surface and normalize the energy. Since only the liver tissues were imaged in 2D form. No depth compensation was applied. 2.6.Statistical AnalysisStatistical analysis and graphical display of data were performed using GraphPad software (version 7.00; GraphPad Software, San Diego, CA, USA). In the ICG results, continuous variables were given as means and SEM in the manually selected ROI region at the CAA. Continuous PA images were registered to the first frame. Similar to the previous study,29 averaged PA signals in the CAA were normalized as the relative PA signal intensity () to reflect the ICG concentration change as follows: is the base PSI before ICG injection which is generated from hemoglobin. is the PSI at time . To quantify the kinetics of change after the ICG injection, an exponential decay model43 was used as follows: where the elimination time () indicates the elapsed time from to the baseline. The rate constant () gives the decay rate.A mean total hemoglobin () and lipid signal () were calculated over anatomical regions. In vivo and ex vivo statistical analyses between two groups (NAFLD versus control) of PA parameters and histological results were conducted using independent samples -tests. A value of was considered statistically significant. The correlation between in vivo and ex vivo lipid signals, between lipid signal and lipid/unit area, were calculated using linear regression. 3.Results3.1.In Vivo ICG Clearance StudiesICG clearance studies were conducted first by real-time PAI for both groups (NAFLD and control group). PAI and US sections of the CAA in vivo were displayed as the region of interest (ROI) in Fig. 3(a) (top and bottom panels, respectively) (yellow dashed circle and arrow). Figure 3(c) shows the in vivo ICG measurements from the CAA. The CCA can be found easily in both the PA and US images. After the injection of ICG at the 60-s time point, a high-intensity peak was observed and gradually decreased over time [Fig. 3(c)]. Statistical results for display a high value in the control group and a lower value in the lipid group. The time to recovery also showed significant differences between the two groups. The control group showed a rapid decrease after the peak point, whereas the elimination of ICG for the NAFLD group was slower [Fig. 3(c)]. Two major parameters were analyzed here. We calculated to show a difference in elimination rate ( versus , versus , ). In addition, the elimination time (), which indicates the elapsed time after the summit of the injection to the baseline, was also analyzed ( versus , versus , ). 3.2.Lipid Detection in a NAFLD ModelICG clearance can serve as a biomarker for liver function but cannot provide lipid content, critical for accessing NAFLD. Hence, following the ICG clearance studies, multispectral PA studies were conducted. US and PA imaging was performed in NAFLD rabbits and the control group. Spectral unmixing was applied to obtain and lipid images [Figs. 4(a) and 4(b)]. From the PA images, distinguishable features are present, indicating different tissue types around the liver. High value regions reveal vascular-rich tissue areas, the dermal layer, and liver vessels. The high lipid signal regions were the subcutaneous lipid and liver. Due to the spectral unmixing, we can find significant differences in lipid concentration across the liver. To select the ROI in the liver region for analysis, ultrasound images were used as a guide to ensure the liver’s location. Then the images showed distinguishable features which helped us to find the liver region below the skin layer and superficial vessels. Moreover, due to the limitation of our linear unmixing method,39 a rectangle region was selected close to the top of the liver to calculate the mean lipid values () and lipid/HbT ratio () in this region. A significant difference in the value was observed ( arbitrary units (a.u.) control versus NAFLD a.u., ). Second, when comparing total hemoglobin between control and NAFLD groups ( a.u. versus a.u., ), no significant difference was found. The lipid to hemoglobin ratio per independent rabbit between groups was also computed ( a.u. versus a.u., ) and a showed significant difference. 3.3.Ex Vivo Fatty Liver ImagingTo verify the in vivo findings, corresponding tissue specimens were taken from both groups. Representative PA images of and lipid were shown in Fig. 5(a) and corresponding histological results were shown in Fig. 5(b). Similar to in vivo analysis, mean lipid values () and lipid/HbT ratio () were compared between two groups. First, a significant difference in was observed ( versus a.u., ). Second, a significant difference was also seen in total hemoglobin between the two groups ( versus a.u., ). Finally, lipid/HbT ratio showed a significant difference between the two groups ( versus a.u., ). 3.4.Relationship Between PA Signal and HistologyWe also compared the PA signal (i.e., ) with histopathological classification. The details of the histopathological scores on steatosis (0 to 3), lobular inflammation (0 to 3), and ballooning (0 to 2) of each rabbit are shown in Table 1. Based on the NAS score, NAFLD group rabbits were classified as “NAFLD but not-NASH” (), “borderline-NASH” ( to 4), and “definite-NASH” ( to 8). Table 1Detailed histopathological scores and PA values NAFLD rabbits and control group (C1 to C3).
As shown in Fig. 6(a), there was a progressive elevation in the in vivo with increased severity from not-NASH to borderline-NASH. There was a significant difference in between control and not-NASH ( versus a.u. ). The increased significantly in borderline-NASH groups and showed significant differences between not-NASH and borderline-NASH ( versus a.u. ). Significant differences were also shown between control and borderline-NASH ( versus a.u., ). In the ex vivo results, significant differences in were shown between control and not-NASH ( versus a.u. ). As the increased, significant differences were shown between control and borderline-NASH ( versus a.u. ). The between not-NASH and borderline-NASH also showed significant differences ( versus a.u. ). According to data in Table 1, the linear regression between ex vivo and in vivo lipid data is shown in Fig. 6(c) (, ). To evaluate the correlation between PA lipid estimation and the histology estimates, we calculated the relative area of lipid/unit area in each staining slice. Figures 6(d) and 6(e) showed that lipid value from PA images has a relatively high correlation with the lipid area. In vivo mean lipid values () had a good correlation with histology estimates (, ), and ex vivo were well correlated with histology estimates (, ). The H&E staining images showed normal hepatic tissue and steatosis pathology. 4.Discussion and ConclusionIn our previous work, we demonstrated the capability of PAI-based ICG clearance for LFR assessment in rabbits,29 the accuracy of which was confirmed by the gold standard method spectrophotometry. In this work, we applied PAI/ICG to assess LFR in rabbits with NAFLD. We have shown that time-resolved PA signal can differentiate the kinetics of ICG uptake and excretion in control and NAFLD groups in vivo, indicating LFR impairment occurred during the development process of NAFLD. Although hepatic ICG clearance serves as a biomarker for liver function in the assessment of NAFLD rabbit, NAFLD cannot be quantified by it. The pathology of NAFLD is a spectrum of lesions ranging from pure steatosis to a complex pattern with significant necroinflammatory injury, fibrosis, or even cirrhosis.44 Therefore, LFR assessment was necessary but not enough for the diagnosis and treatment of the pathological change of NAFLD. Multispectral PAI is a target-specific, quantitative, noninvasive, imaging modality38 when compared with established imaging modalities, including ultrasound imaging and MRI. To obtain more specific information about NAFLD liver, in vivo and ex vivo multispectral PAI were performed and the results from our study demonstrated a significant increase in lipid content in the liver of the NAFLD group compared to the control group, which agreed well with the pathologic results, indicating lipid metabolism malfunction in the NAFLD liver tissue.26,45 Both in vivo and ex vivo mean lipid values displayed a high correlation with histology estimates, which indicated the accuracy of using PAT for steatosis assessment of NAFLD. This was in line with previous studies in mice of NAFLD.27,46 In addition, the linear correlation between lipid in vivo and ex vivo mean lipid is also relatively high. The divergent outcomes for the in vivo and ex vivo HbT results may come from several reasons. First, the liver has a strong signal as it is a vascular-rich tissue and stands for strong endogenous absorbers in vivo. In an ex vivo setting, signal decreased drastically due to blood loss in the tissue. In general, pathological changes associated with steatosis could be unequivocally detected using multispectral PAI both in vivo and ex vivo, consistent with similar findings from recent studies in mice.46 In addition, the pathological changes in those two situations (in vivo and ex vivo), we note that there are some limitations in our current study. The light absorption of lipid in the 760- to 930-nm range is much weaker than that at 1200 and 1720 nm which might be a better wavelength to detect lipid regions.47 As the limitation of our laser system, output energy at 1200 and 1720 nm is too low for PA imaging. As an alternative method, we refer to previous research21,38 to use 930 nm and other wavelengths (e.g., 760, 800, 850, and 930 nm). Those wavelengths with some water absorption may induce inconsistent lipid areas as shown in Fig. 4(b) and show some “patchy” in the NAFLD liver. Although we have a good correlation between PA lipid value and histology lipid area, we should notice that the is below 0.9 and more studies are needed. For example, we can briefly calculate the lipid area in H&E staining and compare it with our PA image results. Red oil O staining will be a better choice to estimate the lipid area. In addition, we need to optimize our system settings including the illumination geometry and fiber bundle parameters in both in vivo and ex vivo setups. Except for the hardware, since we applied zero-thresholding to remove the negative signals in figures, some artifacts in reconstructed images may affect the quantitative results.48 A recently developed method based on Hilbert transformation can be applied in the future to remove such noise and give more strict quantitative measurements.31,48 In addition, the linear unmixing method may be another reason which affects our results. Advanced spectral unmixing methods21,49 need to be developed to achieve the capability of refined quantification of PA signal. In sum, this study showed significantly injured LFR in the NAFLD group compared to the control group. A future study with a larger sample size is needed to investigate the LFR change during various NAFLD pathologies. Due to the limitation of the sample size used in this study, refined PA classification of the NAFLD is not achieved. We are developing an advanced NAFLD animal model with various pathology states and upgrading our PAI system for this challenging study. AcknowledgmentThis study was supported by the National Natural Science Foundation of China (NSFC) (Grant No. 82001833) and the Post-Doctor Research Project, West China Hospital, Sichuan University (Grant No. 2019HXBH014). Code, Data, and Materials AvailabilityData used in the research can be made available to obtain the results reported in the article; please contact the corresponding author. ReferencesQ. M. Anstee, D. Seth and C. P. Day,
“Genetic factors that affect risk of alcoholic and nonalcoholic fatty liver disease,”
Gastroenterology, 150
(8), 1728
–1744.e7 https://doi.org/10.1053/j.gastro.2016.01.037 GASTAB 0016-5085
(2016).
Google Scholar
P. S. Dulai, C. B. Sirlin and R. Loomba,
“MRI and MRE for non-invasive quantitative assessment of hepatic steatosis and fibrosis in NAFLD and NASH: clinical trials to clinical practice,”
J. Hepatol., 65
(5), 1006
–1016 https://doi.org/10.1016/j.jhep.2016.06.005
(2016).
Google Scholar
Z. M. Younossi et al.,
“Global epidemiology of nonalcoholic fatty liver disease—meta-analytic assessment of prevalence, incidence, and outcomes,”
Hepatology, 64
(1), 73
–84 https://doi.org/10.1002/hep.28431 HEPADF 0161-0538
(2016).
Google Scholar
EASL-EASD-EASO,
“Clinical practice guidelines for the management of non-alcoholic fatty liver disease,”
Obesity Facts, 9
(2), 65
–90 https://doi.org/10.1159/000443344
(2016).
Google Scholar
V. Ratziu et al.,
“Sampling variability of liver biopsy in nonalcoholic fatty liver disease,”
Gastroenterology, 128
(7), 1898
–1906 https://doi.org/10.1053/j.gastro.2005.03.084 GASTAB 0016-5085
(2005).
Google Scholar
Z. M. Younossi et al.,
“Nonalcoholic fatty liver disease: assessment of variability in pathologic interpretations,”
Mod. Pathol., 11
(6), 560
–565 MODPEO 0893-3952
(1998).
Google Scholar
Z. M. Younossi,
“Long-term outcomes of nonalcoholic fatty liver disease: from nonalcoholic steatohepatitis to nonalcoholic steatofibrosis,”
Clin. Gastroenterol. Hepatol., 15
(8), 1144
–1147 https://doi.org/10.1016/j.cgh.2017.05.029
(2017).
Google Scholar
L. Castera, V. Vilgrain and P. Angulo,
“Noninvasive evaluation of NAFLD,”
Nat. Rev. Gastroenterol. Hepatol., 10
(11), 666
–675 https://doi.org/10.1038/nrgastro.2013.175
(2013).
Google Scholar
Z. Permutt et al.,
“Correlation between liver histology and novel magnetic resonance imaging in adult patients with non-alcoholic fatty liver disease – MRI accurately quantifies hepatic steatosis in NAFLD,”
Alimentary Pharmacol. Ther., 36
(1), 22
–29 https://doi.org/10.1111/j.1365-2036.2012.05121.x APTHEN 1365-2036
(2012).
Google Scholar
M. Noureddin et al.,
“Utility of magnetic resonance imaging versus histology for quantifying changes in liver fat in nonalcoholic fatty liver disease trials,”
Hepatology, 58
(6), 1930
–1940 https://doi.org/10.1002/hep.26455 HEPADF 0161-0538
(2013).
Google Scholar
C. C. Park et al.,
“Magnetic resonance elastography vs transient elastography in detection of fibrosis and noninvasive measurement of steatosis in patients with biopsy-proven nonalcoholic fatty liver disease,”
Gastroenterology, 152
(3), 598
–607.e2 https://doi.org/10.1053/j.gastro.2016.10.026 GASTAB 0016-5085
(2017).
Google Scholar
K. Okubo et al.,
“Visualization of quantitative lipid distribution in mouse liver through near-infrared hyperspectral imaging,”
Biomed. Opt. Express, 12
(2), 823
–835 https://doi.org/10.1364/BOE.413712 BOEICL 2156-7085
(2021).
Google Scholar
L. Castera, M. Friedrich-Rust and R. Loomba,
“Noninvasive assessment of liver disease in patients with nonalcoholic fatty liver disease,”
Gastroenterology, 156
(5), 1264
–1281.e4 https://doi.org/10.1053/j.gastro.2018.12.036 GASTAB 0016-5085
(2019).
Google Scholar
H. Imamura et al.,
“Assessment of hepatic reserve for indication of hepatic resection: decision tree incorporating indocyanine green test,”
J. Hepato-Biliary-Pancreatic Surg., 12
(1), 16
–22 https://doi.org/10.1007/s00534-004-0965-9
(2005).
Google Scholar
Y.-Y. Wang et al.,
“Comparison of the ability of Child-Pugh score, MELD score, and ICG-R15 to assess preoperative hepatic functional reserve in patients with hepatocellular carcinoma,”
J. Surg. Oncol., 118
(3), 440
–445 https://doi.org/10.1002/jso.25184 JSONAU 0022-4790
(2018).
Google Scholar
Q. Liu et al.,
“Lipid metabolism and functional assessment of discarded human livers with steatosis undergoing 24 hours of normothermic machine perfusion,”
Liver Transplant., 24
(2), 233
–245 https://doi.org/10.1002/lt.24972
(2018).
Google Scholar
A. Becker et al.,
“Multispectral optoacoustic tomography of the human breast: characterisation of healthy tissue and malignant lesions using a hybrid ultrasound-optoacoustic approach,”
Eur. Radiol., 28
(2), 602
–609 https://doi.org/10.1007/s00330-017-5002-x
(2018).
Google Scholar
M. Yang et al.,
“Photoacoustic/ultrasound dual imaging of human thyroid cancers: an initial clinical study,”
Biomed. Opt. Express, 8
(7), 3449
–3457 https://doi.org/10.1364/BOE.8.003449 BOEICL 2156-7085
(2017).
Google Scholar
J. R. Rajian et al.,
“Characterization and treatment monitoring of inflammatory arthritis by photoacoustic imaging: a study on adjuvant-induced arthritis rat model,”
Biomed. Opt. Express, 4
(6), 900
–908 https://doi.org/10.1364/BOE.4.000900 BOEICL 2156-7085
(2013).
Google Scholar
F. Knieling et al.,
“Multispectral optoacoustic tomography for assessment of Crohn’s disease activity,”
N. Engl. J. Med., 376
(13), 1292
–1294 https://doi.org/10.1056/NEJMc1612455 NEJMAG 0028-4793
(2017).
Google Scholar
G. Diot et al.,
“Multispectral optoacoustic tomography (MSOT) of human breast cancer,”
Clin. Cancer Res., 23
(22), 6912 https://doi.org/10.1158/1078-0432.CCR-16-3200
(2017).
Google Scholar
L. V. Wang and S. Hu,
“Photoacoustic tomography: in vivo imaging from organelles to organs,”
Science, 335
(6075), 1458 https://doi.org/10.1126/science.1216210 SCIEAS 0036-8075
(2012).
Google Scholar
X. Shu et al.,
“Quantifying melanin concentration in retinal pigment epithelium using broadband photoacoustic microscopy,”
Biomed. Opt. Express, 8
(6), 2851
–2865 https://doi.org/10.1364/BOE.8.002851 BOEICL 2156-7085
(2017).
Google Scholar
G. S. Sangha, E. H. Phillips and C. J. Goergen,
“In vivo photoacoustic lipid imaging in mice using the second near-infrared window,”
Biomed. Opt. Express, 8
(2), 736
–742 https://doi.org/10.1364/BOE.8.000736 BOEICL 2156-7085
(2017).
Google Scholar
Y. Yan et al.,
“Photoacoustic imaging of the uterine cervix to assess collagen and water content changes in murine pregnancy,”
Biomed. Opt. Express, 10
(9), 4643
–4655 https://doi.org/10.1364/BOE.10.004643 BOEICL 2156-7085
(2019).
Google Scholar
G. Xu et al.,
“The functional pitch of an organ: quantification of tissue texture with photoacoustic spectrum analysis,”
Radiology, 271
(1), 248
–254 https://doi.org/10.1148/radiol.13130777 RADLAX 0033-8419
(2014).
Google Scholar
O. Rom et al.,
“Nitro-fatty acids protect against steatosis and fibrosis during development of nonalcoholic fatty liver disease in mice,”
EBioMedicine, 41 62
–72 https://doi.org/10.1016/j.ebiom.2019.02.019
(2019).
Google Scholar
C. Tian et al.,
“Imaging and sensing based on dual-pulse nonlinear photoacoustic contrast: a preliminary study on fatty liver,”
Opt. Lett., 40
(10), 2253
–2256 https://doi.org/10.1364/OL.40.002253 OPLEDP 0146-9592
(2015).
Google Scholar
T. Qiu et al.,
“Assessment of liver function reserve by photoacoustic tomography: a feasibility study,”
Biomed. Opt. Express, 11
(7), 3985
–3995 https://doi.org/10.1364/BOE.394344 BOEICL 2156-7085
(2020).
Google Scholar
P. Hu et al.,
“Spatiotemporal antialiasing in photoacoustic computed tomography,”
IEEE Trans. Med. Imaging, 39
(11), 3535
–3547 https://doi.org/10.1109/TMI.2020.2998509 ITMID4 0278-0062
(2020).
Google Scholar
L. Li et al.,
“Single-impulse panoramic photoacoustic computed tomography of small-animal whole-body dynamics at high spatiotemporal resolution,”
Nat. Biomed. Eng., 1
(5), 0071 https://doi.org/10.1038/s41551-017-0071
(2017).
Google Scholar
M. A. Van Herck, L. Vonghia and S. M. Francque,
“Animal models of nonalcoholic fatty liver disease—a starter’s guide,”
Nutrients, 9
(10), 1072 https://doi.org/10.3390/nu9101072
(2017).
Google Scholar
M. Hall et al.,
“The WEKA data mining software: an update,”
ACM SIGKDD Explor. Newslett., 11
(1), 10
–18 https://doi.org/10.1145/1656274.1656278
(2009).
Google Scholar
M. Xu and L. V. Wang,
“Universal back-projection algorithm for photoacoustic computed tomography,”
Phys. Rev. E, 71
(1), 016706 https://doi.org/10.1103/PhysRevE.71.016706
(2005).
Google Scholar
J. Yang et al.,
“Photoacoustic assessment of hemodynamic changes in foot vessels,”
J. Biophotonics, 12
(6), e201900004 https://doi.org/10.1002/jbio.201900004
(2019).
Google Scholar
T. F. Fehm et al.,
“In vivo whole-body optoacoustic scanner with real-time volumetric imaging capacity,”
Optica, 3
(11), 1153
–1159 https://doi.org/10.1364/OPTICA.3.001153
(2016).
Google Scholar
V. Ntziachristos and D. Razansky,
“Molecular imaging by means of multispectral optoacoustic tomography (MSOT),”
Chem. Rev., 110
(5), 2783
–2794 https://doi.org/10.1021/cr9002566 CHREAY 0009-2665
(2010).
Google Scholar
A. P. Regensburger et al.,
“Detection of collagens by multispectral optoacoustic tomography as an imaging biomarker for Duchenne muscular dystrophy,”
Nat. Med., 25
(12), 1905
–1915 https://doi.org/10.1038/s41591-019-0669-y 1078-8956
(2019).
Google Scholar
H. Roman et al.,
“Estimating blood oxygenation from photoacoustic images: can a simple linear spectroscopic inversion ever work?,”
J. Biomed. Opt., 24
(12), 121914 https://doi.org/10.1117/1.JBO.24.12.121914 JBOPFO 1083-3668
(2019).
Google Scholar
S. Park et al.,
“Normalization of optical fluence distribution for three-dimensional functional optoacoustic tomography of the breast,”
J. Biomed. Opt., 27
(3), 036001 https://doi.org/10.1117/1.JBO.27.3.036001 JBOPFO 1083-3668
(2022).
Google Scholar
P. Parsa, S. L. Jacques and N. S. Nishioka,
“Optical properties of rat liver between 350 and 2200 nm,”
Appl. Opt., 28
(12), 2325
–2330 https://doi.org/10.1364/AO.28.002325 APOPAI 0003-6935
(1989).
Google Scholar
S. Kumari and A. K. Nirala,
“Study of light propagation in human, rabbit and rat liver tissue by Monte Carlo simulation,”
Optik, 122
(9), 807
–810 https://doi.org/10.1016/j.ijleo.2010.06.006 OTIKAJ 0030-4026
(2011).
Google Scholar
N. Brillant et al.,
“Dynamic and accurate assessment of acetaminophen-induced hepatotoxicity by integrated photoacoustic imaging and mechanistic biomarkers in vivo,”
Toxicol. Appl. Pharmacol., 332 64
–74 https://doi.org/10.1016/j.taap.2017.07.019 TXAPA9 0041-008X
(2017).
Google Scholar
S. L. Friedman et al.,
“Mechanisms of NAFLD development and therapeutic strategies,”
Nat. Med., 24
(7), 908
–922 https://doi.org/10.1038/s41591-018-0104-9 1078-8956
(2018).
Google Scholar
K. Pei et al.,
“An overview of lipid metabolism and nonalcoholic fatty liver disease,”
BioMed. Res. Int., 2020 4020249 https://doi.org/10.1155/2020/4020249
(2020).
Google Scholar
S. Huang et al.,
“Functional multispectral optoacoustic tomography imaging of hepatic steatosis development in mice,”
EMBO Mol. Med., 13 e13490 https://doi.org/10.15252/emmm.202013490
(2021).
Google Scholar
T. Allen et al.,
“Spectroscopic photoacoustic imaging of lipid-rich plaques in the human aorta in the 740 to 1400 nm wavelength range,”
J. Biomed. Opt., 17
(6), 061209 https://doi.org/10.1117/1.JBO.17.6.061209 JBOPFO 1083-3668
(2012).
Google Scholar
G. Li et al.,
“Multiview Hilbert transformation for full-view photoacoustic computed tomography using a linear array,”
J. Biomed. Opt., 20
(6), 066010 https://doi.org/10.1117/1.JBO.20.6.066010 JBOPFO 1083-3668
(2015).
Google Scholar
S. Tzoumas et al.,
“Effects of multispectral excitation on the sensitivity of molecular optoacoustic imaging,”
J. Biophotonics, 8
(8), 629
–637 https://doi.org/10.1002/jbio.201400056
(2015).
Google Scholar
|
CITATIONS
Cited by 1 scholarly publication.
Liver
In vivo imaging
Photoacoustic imaging
Tissues
Animal model studies
Imaging systems
Animals