CC BY-NC-ND 4.0 · Endosc Int Open 2017; 05(08): E775-E783
DOI: 10.1055/s-0043-111790
Original article
Eigentümer und Copyright ©Georg Thieme Verlag KG 2017

A feasibility study of photoacoustic imaging of ex vivo endoscopic mucosal resection tissues from Barrett’s esophagus patients

Liang Lim
1   Princess Margaret Cancer Centre, Toronto, Ontario, Canada
,
Catherine J. Streutker
2   St. Michael’s Hospital, Toronto, Ontario, Canada
,
Norman Marcon
2   St. Michael’s Hospital, Toronto, Ontario, Canada
,
Maria Cirocco
2   St. Michael’s Hospital, Toronto, Ontario, Canada
,
Alexandra Lao
2   St. Michael’s Hospital, Toronto, Ontario, Canada
,
Vladimir V. Iakovlev
2   St. Michael’s Hospital, Toronto, Ontario, Canada
,
Ralph DaCosta
1   Princess Margaret Cancer Centre, Toronto, Ontario, Canada
,
Brian C. Wilson
1   Princess Margaret Cancer Centre, Toronto, Ontario, Canada
› Author Affiliations
Further Information

Corresponding author

Liang Lim, PhD
University Health Network – Princess Margaret Cancer Centre
101 College Street
PMCRT #15-301V Toronto
Ontario M5G 1L7
Canada    

Publication History

submitted 20 December 2016

accepted after revision 02 May 2017

Publication Date:
07 August 2017 (online)

 

Abstract

Background and study aims Accurate endoscopic detection of dysplasia in patients with Barrett’s esophagus (BE) remains a major clinical challenge. The current standard is to take multiple biopsies under endoscopic image guidance, but this leaves the majority of the tissue unsampled, leading to significant risk of missing dysplasia. Furthermore, determining whether there is submucosal invasion is essential for proper staging. Hence, there is a clinical need for a rapid in vivo wide-field imaging method to identify dysplasia in BE, with the capability of imaging beyond the mucosal layer. We conducted an ex vivo feasibility study using photoacoustic imaging (PAI) in patients undergoing endoscopic mucosal resection (EMR) for known dysplasia. The objective was to characterize the esophageal microvascular pattern, with the long-term goal of performing in vivo endoscopic PAI for dysplasia detection and therapeutic guidance.

Materials and methods EMR tissues were mounted luminal side up. The tissues were scanned over a field of view of 14 mm (width) by 15 mm (depth) at 680, 750, and 850 nm (40 MHz acoustic central frequency). Ultrasound and photoacoustic images were simultaneously acquired. Tissues were then sliced and fixed in formalin for histopathology with hematoxylin and eosin staining. A total of 13 EMR specimens from eight patients were included in the analysis, which consisted of co-registration of the photoacoustic images with corresponding pathologist-classified histological images. We conducted mean difference test of the total hemoglobin distribution between tissue classes.

Results Dysplastic and nondysplastic BE can be distinguished from squamous tissue in 84 % of region-of-interest comparisons (42/50). However, the ability of intrinsic PAI to distinguish dysplasia from NDBE, which is the clinically important challenge, was only about 33 % (10/30).

Conclusion We demonstrated the technical feasibility of this approach. Based on our ex vivo data, changes in total hemoglobin content from intrinsic PAI (i. e. without exogenous contrast) can differentiate BE from squamous esophageal mucosa. However, most likely intrinsic PAI is unable to differentiate dysplastic from nondysplastic BE with adequate sensitivity for clinical translation.


#

Introduction

In North America approximately 19 000 new cases of esophageal cancer are diagnosed annually. The 5-year survival rate is only about 15 %, with 18 000 cancer-related deaths annually [1] [2]. The most important risk factor for esophageal adenocarcinoma is Barrett’s esophagus (BE), which affects some 3.5 million adults in North America. This is a condition in which the outermost layer of esophageal squamous mucosa (whitish in appearance on standard endoscopy) transforms to gastric/intestinal-type mucosa (pink-red appearance) due to chronic gastroesophageal reflux. Patients with BE have a > 30-fold increased risk of developing esophageal adenocarcinoma [3].

The current clinical standard for esophageal cancer diagnosis in patients with BE is targeted biopsies of focal visible abnormalities followed by 4-quadrant random biopsy every 1 – 2 cm along the length of the Barrett’s segment under white-light endoscopic guidance. However, the majority of the tissue remains unsampled, leading to significant risk of missing dysplasia, especially in areas with no visible abnormalities. Biopsied BE tissue is then staged by histopathology: in order of cancer progression, from nondysplastic BE (NDBE), to low grade dysplasia (LGD), high grade dysplasia (HGD) and, ultimately, cancer (e. g. adenocarcinoma). Whereas less than 10 % of BE patients progress to HGD or cancer, patients with HGD are at considerably increased risk of progressing to cancer [4], with yearly rates of 4 % – 11.8 % [5] [6].

Early detection of dysplasia, especially HGD, is crucial in managing these patients. Newer methodologies such as cytology from brushes [7] and sponges [8] are being developed but are not yet in common usage. Several clinical trials with a variety of chromoendoscopy have shown promising results [9]. However, chromoendoscopy only enhances the surface features of the tissue and does not provide submucosal information. As well as detecting HGD and intramucosal adenocarcinoma, determining whether there is submucosal invasion is essential for proper staging. Thus, the absence of submucosal invasion indicates that minimally invasive endoscopic mucosal resection (EMR) is potentially curative, whereas its presence indicates the need for esophagectomy, which results in significant mortality as well as detriment to quality of life. Hence, there is an unmet need for a rapid, wide-field endoscopic imaging method to identify dysplasia in BE, with the capability of imaging beyond the mucosal layer.

A wide variety of optical imaging and spectroscopic techniques have been investigated to address this need, including diffuse reflectance spectroscopy [10], tissue autofluorescence [11], contrast-based fluorescence [12], Raman spectroscopy [13], optical coherence tomography [14] [15], and narrow-band imaging [16], as well as various combinations of these techniques [17]. However, the clinical need is still largely unmet: for example, optical coherence tomography looks promising in providing large-area scanning at high spatial resolution, but is marginal in terms of depth (~1.5 – 2 mm) to reliably assess submucosal invasion.

The mechanistic rationale for using photoacoustic imaging (PAI) in BE is that studies using narrow-band imaging [16] or diffuse reflectance spectroscopy [10] have reported microvasculature changes between NDBE and dysplasia, and this is consistent with angiogenesis in tumor progression [18] [19]. PAI, based on the detection of acoustic waves generated by local absorption of short (~ns) pulses of laser light, is specifically sensitive to detecting blood vessels in tissue as a result of light absorption by haemoglobin, and is also able to image much deeper into tissue (up to 3.5 cm) [20] than other optical endoscopic imaging techniques [21]. Hence, PAI has the potential to provide imaging of microvasculature from the tissue surface to well below the mucosal layer.

We hypothesize that PAI can detect microvasculature changes between nondysplastic and dysplastic BE. Here, we report a feasibility study using PAI in ex vivo tissues from BE patients with known dysplasia who were undergoing EMR. We have previously reported the technical feasibility of this technique [22], and here we provided analysis from our first group of patients. This analysis consisted of acoustic image co-registration with pathologist-annotated histological images, followed by mean difference comparison between the intrinsic PAI intensity images of the EMR specimen and the pathological classification. Initial results show clearly that intrinsic PAI can differentiate BE from squamous tissue (42/50 = 84 %), but is likely not sufficiently sensitive to diagnose dysplasia in BE.


#

Materials and methods

Photoacoustic system

We used a commercial photoacoustic tomography imaging system (Vevo LAZR Visualsonics, Toronto, Ontario, Canada), which has been described in detail elsewhere [23]. This included an integrated fiberoptic linear-array transducer (LZ550, Visualsonics, 40 MHz central acoustic frequency, 256 elements) and a crossed-laser-beam geometry for optical excitation [23]. Photoacoustic and ultrasound images were collected at a frame rate of 5 Hz. The linear transducer generated images with a field of view of approximately 15 mm (depth) by 14 mm (width). The PAI system was located adjacent to a clinical endoscopy suite, so that tissue specimens could be imaged immediately after endoscopic resection.


#

Patient enrollment

All patients were scheduled for routine screening, diagnostic, and/or therapeutic endoscopy at the Gastroenterology Clinic, St. Michael’s Hospital, Toronto, Canada. Patients were 18 years or above, able to sign the informed consent, and had known pathology indicative for EMR. We excluded patients with any underlying medical condition contraindicating EMR or who had gastrointestinal (GI) malignancy treated with chemotherapy or radiotherapy within the past 3 years. All patients provided written informed consent prior to their endoscopy procedure. The study was carried out with institutional Research Ethics Board approval.

Eight patients were enrolled in this initial feasibility study, from whom PAI and ultrasound scans were obtained in 13 ex vivo EMR tissue samples. Patients’ pre-EMR diagnosis ranged from inflamed mucosa to intramucosal carcinoma. Patients’ post-EMR diagnosis ranged from NDBE to intramucosal carcinoma. EMRs were directed at areas with field defect most suspicious for disease. Detailed information of pre-EMR diagnosis, endoscopic descriptions, and resected specimen diagnosis are included in the supplementary materials (see S1. Patient List).


#

EMR tissue preparation

The study neither changed standard endoscopic procedure nor increased risk to the patients, as the ultrasound and PAI scans were performed on excised tissues. We also confirmed that the imaging procedure had no detectable effect on the tissues and did not compromise subsequent clinical histopathology. When BE with suspected dysplasia is identified during endoscopy, patients may undergo an EMR procedure to excise BE tissues for diagnosis. There are variations of EMR procedures [24]. In the current study, all EMRs were performed using a Duette suction device (Duette; Cook Medical, Bloomington, Indiana, USA), which excises the tissue by suctioning to draw the tissue into the cap, followed by ligation (rubber banding) to a polyp-like shape and resection using an electrosurgical snare.


#

Data acquisition

Excised tissues were mounted using ultrasound gel onto a plastic dish (#430167; Corning Inc., New York, USA) filled with clear gelatin (5 % weight of 40 mL volume) to minimize acoustic signals from the dish. Digital photographs were taken and included a scale bar. Phosphate buffered saline was then added to the dish to provide good acoustic contact with the transducer and to prevent tissue dehydration. A transducer cover (CIV-Flex 610 – 637; Civco, Iowa, USA) was mounted onto the transducer probe to protect the tissue sample and avoid cross contamination. The transducer was lowered in contact with the phosphate buffered saline and the sample was centered on the transducer field of view. Because of the crossed-laser-beam geometry of the PAI fiber bundle attached to the transducer, there is an optimal photoacoustic zone located between 9 and 11 mm in front of the probe face. Hence, we positioned the top surface of the tissue at 9 mm for consistency between samples.

PAI scans were made at each of 680, 750, and 850 nm. Additional measurements were made at 824 and 970 nm in anticipation of future work in which these data will serve as controls in studies of topical biomarker-targeted porphyrin-lipid nanoparticle contrast [25] [26]. The transducer was mounted onto a motorized translational stage and scanned over a 3 D volume at ~0.06 mm step size. The time for complete data acquisition over a 17 mm scan distance was approximately 5 minutes (~1 minute per wavelength). The PAI laser power was below the ANSI maximum permissible exposure for skin (20 – 70 mJ/cm2 for 680 – 970 nm) [27], avoiding any tissue heating that might damage the tissue, as confirmed on subsequent histopathology.


#

Histopathology procedure

Once data acquisition was completed, a pathologist marked the orientation of the specimen for subsequent histopathology and sectioning along the ultrasound/photoacoustic imaging plane. Larger samples were grossly cut into 2 mm thick slabs prior to fine sectioning. Sectioned tissues were assigned a unique identifier and fixed in 10 % formalin. The tissues were kept moist with phosphate buffered saline from the time of excision until fixation in formalin. The whole procedure, from excision to placing the tissue in formalin, was completed within 60 minutes to avoid compromising subsequent clinical histopathology.


#

Data analysis

PAI and ultrasound data were analyzed in Matlab (Mathworks; Natick, Massachusetts, USA). The PAI intensities at 750 nm and 850 nm were used as input to estimate the relative concentrations of oxy-hemoglobin (HbO2) and deoxy-hemoglobin (Hb) using the corresponding known molar extinction coefficients (εHbO2 and εHb) [23] [28] [29]. The relative total hemoglobin (HbT) was then given by the sum of the relative HbO2 and Hb values.

A digitized whole-section hematoxylin and eosin-stained slide from each tissue was segmented and color coded by an expert GI pathologist (C.S. and A.L.) according to the following tissue classes: squamous mucosa (blue), inflammation (gray), intestinal metaplasia (IM; green), dysplasia (yellow), and adenocarcinoma (orange). [Fig. 1] shows two typical histological images ([Fig. 1a,c]), overlaid with the color-coded segmentation. The histopathology slices were then co-registered with the closest corresponding ultrasound and photoacoustic image ([Fig. 1b,d]), and regions of interest (ROIs) for each tissue class were mapped manually onto these acoustic images. PAI with the above probe-tissue geometry was able to detect signal from blood located up to 2 mm beneath the tissue surface.

Zoom Image
Fig. 1 Two examples of co-registered histopathology (hematoxylin and eosin staining) and ultrasound/photoacoustic images. a, c Histopathology overlaid with color-coded tissue classification (green: intestinal metaplasia [IM]; blue: squamous [SQ]; yellow: dysplasia [Dys]; gray: inflammation). The black lines are 2 mm scale bars. b, d Corresponding ultrasound (grayscale) and photoacoustic images (blue: 750 nm; red 850 nm). Regions of interest for different tissue classes are indicated. Scales are in mm.

The ROI comparisons were limited to regions within 9 – 11 mm from the transducer surface in order to minimize the effect of nonuniform distribution of the excitation light, and regions within the same field of view to minimize the effects of intra- and inter-patient and site heterogeneity in BE, which is known to be high. To make PAI comparison across patients and specimens without accounting for this heterogeneity would be misleading; therefore, histological mapping was done within each histological slice ([Fig. 1a,c]) and comparisons were made only within the same field of view. Comparisons were made for all possible ROI permutations within each image: for example, for the case shown in  [Fig. 1a,b], these were ROI-IM vs. ROI-SQ1, ROI-IM vs. ROI-SQ2, and ROI-SQ1 vs. ROI-SQ2, but the third comparison was not included in the final analysis as both ROI-SQ1 and ROI-SQ2 belonged to the same diagnostic class (squamous).

[Fig. 2] shows an example of comparison between two ROIs from [Fig. 1b]. Here, ROI-IM (green) had a higher mean relative HbT value than ROI-SQ1 (blue), the difference being statistically significant (Wilcoxon Rank Sum test, P < 0.05). We considered ROI-IM to have higher relative HbT than that of ROI-SQ1 when the mean HbT was higher and the difference in their HbT distributions was statistically significant (P < 0.05). If these conditions did not hold, then the relative HbT comparison for that ROI pair was considered to be negative.

Zoom Image
Fig. 2 Example of frequency histograms of photoacoustic imaging intensity for two of the regions of interest (ROIs) shown in [Fig. 1b], corresponding to the middle green ROI (intestinal metaplasia [IM]) and left blue ROI (squamous mucosa [SQ1). a 750 nm. b 850 nm. c Relative total hemoglobin.

#
#

Results

 [Fig. 3] shows 3 D reconstruction of a complete EMR tissue volume using the acoustic images from a typical EMR tissue. Ultrasound ( [Fig. 3a])provided structural information, while PAI ([Fig. 3b,c]) provided complementary information about the distribution of absorbers (i. e. blood) within the tissue.

Zoom Image
Fig. 3 3D reconstruction of the photoacoustic imaging (PAI) signal. a Ultrasound B-mode (grayscale). b PAI at 750 nm (blue), mainly dominated by hemoglobin absorption. c PAI at 850 nm (red), mainly dominated by oxy-hemoglobin absorption. Axis scale in mm.

The ROI comparisons of relative HbT are summarized in  [Table 1]. In general, the relative HbT differences between the subclasses of BE (dysplasia, NDBE with inflammation, and IM) were not consistent. However, relative HbT was higher in BE than in squamous mucosa in the majority of comparisons.

Table 1

Summary of region-of-interest comparisons of relative total hemoglobin between the subclasses of Barrett’s esophagus.

HbT

Dysplasia vs. NDBE

NDBE

BE vs. Squamous

Dysplasia > Inflammation[*]

Dysplasia > IM

Inflammation[*] > IM

Inflammation[*] > Squamous

Dysplasia > Squamous

IM > Squamous

Positive

 6

 4

 8

20

 9

13

Negative

12

 8

19

 4

 1

 3

Total

18

12

27

24

10

16

HbT, total hemoglobin; BE, Barrett’s esophagus; NDBE, nondysplastic BE; IM, intestinal metaplasia.

* NDBE with inflammation.


Based on the relative HbT content, NDBE and BE can be distinguished from squamous tissue in 84 % of ROI comparisons (42/50). This is not surprising, as the chronic inflammation and resulting mucosal transformation in BE is associated with increased vascularity, as evidenced by the marked color difference between BE and squamous mucosa visible on white-light endoscopy. However, the ability of intrinsic PAI to distinguish dysplasia from NDBE, which is the clinically important challenge, was only about 33 % (10/30), even though the wavelengths used were hemoglobin sensitive so that PAI was expected to pick up alterations in the tissue microvascular density.


#

Discussion and conclusions

In terms of technical issues in performing these ex vivo PAI scans, as previously mentioned the crossed-laser-beam geometry of the PAI probe resulted in an optimal photoacoustic zone within the field of view located between 9 and 11 mm from the probe surface. This raised challenges in both data acquisition and analysis. For data acquisition, accurate placement of the tissue sample within the field of view is required to acquire consistent PAI signals from the same tissue volume. However, this is difficult to achieve because the EMR tissue ligation results in a domed or polyp-like specimen shape. Furthermore, the sizes of the EMR tissue samples were highly variable, from < 3 × 3 × 2 mm3 up to 10 × 12 × 6 mm3. Given these constraints, the simplest and most reliable solution was to position the top of the tissue at 9 mm. For data analysis, the crossed-beam geometry has higher spatial variation in intensity than a flat-top beam. Even with simpler beam geometry, the relationship of photoacoustic signal to physiologically relevant parameters is complex, and estimating the hemoglobin-related concentrations is not strictly quantitative, as documented in other studies [30] [31] [32]. Hence, as an initial analysis without requiring complex inverse modeling, we included only ROIs lying within the optimal (9 – 11 mm) PAI region. We also limited the comparisons to ROIs within the same image, so as to mitigate the effects of inter-site and inter-patient variations [33]. Our sample pool did not include tissue with submucosal invasion. However, we consistently acquired PAI data beyond 2 mm from the tissue surface, which supports the imaging depth potential of PAI. Future studies will include samples with submucosal invasion. We note, however, that we have successfully imaged PAI-active nanoparticles lying below full-thickness ex vivo porcine intestine that was 2 mm thick (data not shown).

As regards the main clinical findings, the successful identification of BE relative to squamous mucosa indicated that the experimental protocol and image analysis used were not technically flawed. Furthermore, we repeated the comparison analysis using only the intrinsic PAI intensity (detailed in Supplementary Material S2) and this did not substantially change the results. The observations remained consistent at the ROI, slice, sample, and patient levels, indicating that the findings are not patient- or sample-biased (detailed in Supplementary Material S3). We realized our sample size was less than ideal, but it did serve the purposes of this first feasibility study (power and sample size analysis in Supplementary Material S4).

One might have expected the performance of PAI in identifying BE vs. squamous mucosa to be even better, given how distinct BE appears visually on endoscopy. However, there is very high heterogeneity in BE, and not all regions are equally inflamed and hypervascularized. The more important point was the inability, at least in this ex vivo study, to reliably distinguish dysplasia from NDBE, which was hypothesized to be possible based on altered microvascularity. However, it should be noted that there is a relatively poor concordance between even expert GI pathologists in identifying dysplasia, particularly low grade, in BE owing to considerable histological overlap between inflammatory/reactive changes and premalignant transformation [34].

A further limitation in interpreting the present results is that the ex vivo PAI data may not translate directly to the in vivo setting, primarily because of blood loss from the EMR tissue sample. Although most of the blood content is preserved by the ligation process, the blood does not necessarily fully fill the full vascular volume and there may be redistribution within the sample. The oxygen saturation (SO2) value was omitted from the analysis to account for probable hypoxia post resection. Douplik et al. reported an algorithm using diffuse reflectance spectroscopy that involved a combination of the HbT and SO2 values [10], so that accessing the SO2 information with in vivo PAI at multiple wavelengths may improve the diagnostic accuracy and, in the best-case scenario, result in clinically useful detection of dysplasia. To this end, we are developing an endoscopic PAI system and will use it to conduct an in vivo trial in the near future to confirm or alter the negative conclusion of the present study. A prototype endoscopic PAI instrument has been reported by Yang et al. [30], but to our knowledge the device has been used only in preclinical in vivo models in the esophagus. Even if intrinsic PAI is not adequate for detecting dysplasia, it should still have value in endoscopic staging of disease, because of its imaging depth capability.

Finally, if the intrinsic differences in microvascularity of dysplasia in BE does indeed prevent adequate endoscopic detectability using intrinsic PAI, then photoacoustic contrast agents may be a useful addition. Hence, in parallel, we are investigating the utility of PAI in combination with nontargeted and biomarker-targeted porphyrin-lipid nanoparticles (porphysomes) [25] [26], which may be administered systemically or by topical application to the tissue to minimize cost and potential toxicity.


#

Supplementary Materials

S1. Patient List

Table S1

Detailed list of pre-EMR diagnosis, endoscopic descriptions, and resected specimens’ diagnosis.

Patient ID#

Pre-EMR diagnosis

Endoscopic descriptions

Paris Classification

EMR ID#

EMR diagnosis

(Staging – Stolte Classification)

 2

Suspicious for intramucosal carcinoma

Nodular lesion

0-Is

 1

Adenocarcinoma

(T1M2)

 2

IM

 3

Inflamed mucosa, no IM or dysplasia identified

Nodular GEJ

0-Is

 6

NDBE

12

LGD

 4

HGD

Nodular BE

0-Is

 1

HGD

11

Intramucosal carcinoma

Irregular esophageal mucosa

0-IIa

 4

LGD

 5

LGD

15

LGD

Post-ablation, area nonepithelialized, slightly nodular

0-IIa

 1

NDBE

16

Intramucosal carcinoma

Irregular slightly raised mucosa

0-IIa

 1

Adenocarcinoma (T1M2) + HGD

 3

HGD

17

HGD

Slight nodularity and erosion at GEJ

0-IIa; 0-IIc

 1

NDBE

 2

NDBE

18

HGD

Slightly raised mucosa

0-IIa

 1

Adenocarcinoma

(T1M4) + HGD

EMR, endoscopic mucosal resection; HGD, high grade dysplasia; LGD, low grade dysplasia; GEJ, gastroesophageal junction; BE, Barrett’s esophagus; IM, intestinal metaplasia; NDBE, nondysplastic Barrett’s esophagus;


#

S2. PAI Comparison Analysis using PAI intensity

We repeated the region-of-interest (ROI) analysis using the PAI intensity at 680 nm (maximum laser output) as well as at 750 and 850 nm. In the case shown in  [Fig. 1a,b] and [Fig. 2], we considered the photoacoustic intensity of ROI-1 to be higher than that of ROI-2 (i. e. to have a “positive” score) when 2 or 3 variables have higher mean value and the difference in their frequency distributions is statistically significant (P  < 0.05).  [Table S2] summarizes the corresponding comparisons using all 3 wavelengths. Overall, the findings results did not change significantly from the use of the relative HbT values, so that HbT adequately represents the collected data.

Table S2

Summary of ROI comparison analysis using majority poll of PA intensity.

Score

Dys vs. NDBE

NDBE

BE vs. Squ

Dys > Infl

Dys > IM

Infl > IM

Infl > Squ

Dys > Squ

IM > Squ

Positive

 6

 4

13

20

 9

14

Negative

12

 8

14

 4

 1

 2

Total

18

12

27

24

10

16

Dys-dysplasia, NDBE, nondysplastic Barrett’s esophagus; Squ, squamous; Infl, NDBE with inflammation; IM, intestinal metaplasia.


#

S3. Comparison Level: ROI vs. Slice vs. Sample vs. Patient

We acquired data from at least one tissue sample from each patient. Each sample can have multiple slices, depending on the size. Within each slice, multiple ROI comparisons are possible, depending on the number of permutations between ROIs. One potential problem is then that the results could be biased by a subgroup of tissue samples or patients if these contributed a major fraction of the data. To check this, we carried out the comparison analysis at each of the ROI, tissue slice, tissue sample, and patient levels. At the ROI level, each comparison was counted as an observation. At the tissue slice level, ROI comparisons from the same slice were combined, with the majority of the observations determining the outcome for that slice. At the tissue sample level, observations from all slices from the sample were combined, and at the patient level observations from all the tissue sample(s) were combined. The most important comparisons are summarized in  [Table S3].

Table S3

Comparison analysis at different levels.

ROI

Slice

Sample

Patient

Dys vs. NDBE

Positive

10

 3

1

0

Negative

20

 5

4

3

Same

 0

 3

0

1

Total

30

11

5

4

BE vs. Squ

Positive

42

12

6

5

Negative

 8

 3

2

1

Same

 0

 1

0

0

Total

50

16

8

6

Dys, dysplasia; NDBE, nondysplastic Barrett’s esophagus; BE, Barrett’s esophagus; Squ, squamous; ROI, region of interest.

Again, the relative HbT in BE was higher than in squamous esophagus at least 75 % of the time at every level, whereas it was not consistent at any level between dysplasia and NDBE.


#

S4. Power and Sample Size Analysis

Notwithstanding the small number of patients in this feasibility study, we conducted a power analysis to compare the mean difference between ROIs. With expected 50 observations (at least 50 data points/pixels for each ROI) in each ROI, we have at least 90 % power to detect significant difference between the two ROIs. The calculation is based on the assumed effect size (mean difference) at 1 SD. In order to adjust for multiple comparisons, the significance level is set at 0.0005 (0.05/100, as we conducted 107 tests in total).

We also conducted sample size calculation for a future pilot study. If the power for a single sample test is 0.90 as we conducted above, a total sample size of 35 EMR samples can produce a two-sided 95 % confidence interval with a width equal to 0.2 (0.80 – 1.00).


#
#
#

Competing interests

None

Acknowledgments

This work was supported by an Emerging Frontiers grant (#1022) from the Terry Fox Research Institute. The authors also thank Visualsonics Inc. for equipment support.

  • References

  • 1 Canada Cancer Society. Canadian Cancer Statistics 2014. Available from: http://www.cancer.ca/en/support-and-services/resources/publications/ Accessed 24 Feb 2015
  • 2 American Cancer Society. Cancer Facts & Figures 2015. Available from: https://www.cancer.org/research/cancer-facts-statistics/all-cancer-facts-figures/ Accessed 24 Feb 2015
  • 3 Solaymani-Dodaran M, Logan RFA, West J. et al. Risk of oesophageal cancer in Barrett’s oesophagus and gastro-oesophageal reflux. Gut 2004; 53: 1070-1074
  • 4 Shaheen NJ, Richter JE. Barrett’s oesophagus. Lancet 2009; 373: 850-861
  • 5 Schnell TG, Sontag SJ, Chejfec G. et al. Long-term nonsurgical management of Barrett’s esophagus with high-grade dysplasia. Gastroenterology 2001; 120: 1607-1619
  • 6 Reid BJ, Levine DS, Longton G. et al. Predictors of progression to cancer in Barrett’s esophagus: baseline histology and flow cytometry identify low- and high-risk patient subsets. Am J Gastroenterol 2000; 95: 1669-1676
  • 7 Vennalaganti PR, Naag KanakadandiV, Gross SA. et al. Inter-observer agreement among pathologists using wide-area transepithelial sampling with computer-assisted analysis in patients with Barrett’s esophagus. Am J Gastroenterol 2015; 110: 1257-1260
  • 8 Haisley KR, Dolan JP, Olson SB. et al. Sponge sampling with fluorescent in situ hybridization as a screening tool for the early detection of esophageal cancer. J Gastrointest Surg 2017; 21: 215-221
  • 9 Mannath J, Ragunath K. Role of endoscopy in early oesophageal cancer. Nat Publ Gr 2016; 13: 720-730
  • 10 Douplik A, Zanati S, Saiko G. et al. Diffuse reflectance spectroscopy in Barrett’s Esophagus: developing a large field-of-view screening method discriminating dysplasia from metaplasia. J Biophotonics 2014; 7: 304-311
  • 11 Kara M, DaCosta RS, Wilson BC. et al. Autofluorescence-based detection of early neoplasia in patients with Barrett’s esophagus. Dig Dis 2004; 22: 134-141
  • 12 Messmann H, Endlicher E, Freunek G. et al. Fluorescence endoscopy for the detection of low and high grade dysplasia in ulcerative colitis using systemic or local 5-aminolaevulinic acid sensitisation. Gut 2003; 52: 1003-1007
  • 13 Bergholt MS, Zheng W, Lin K. et al. In vivo diagnosis of esophageal cancer using image-guided Raman endoscopy and biomolecular modeling. Technol Cancer Res Treat 2011; 10: 103-112
  • 14 Yang VXD, Tang SJ, Gordon ML. et al. Endoscopic Doppler optical coherence tomography in the human GI tract: Initial experience. Gastrointest Endosc 2005; 61: 879-890
  • 15 Bouma B, Tearney G. High-resolution imaging of the human esophagus and stomach in vivo using optical coherence tomography. Gastrointest Endosc 2000; 51: 467-474
  • 16 Singh R, Shahzad MA, Tam W. et al. Preliminary feasibility study using a novel narrow-band imaging system with dual focus magnification capability in Barrett’s esophagus: Is the time ripe to abandon random biopsies?. Dig Endosc 2013; 25: 151-156
  • 17 Georgakoudi I, Jacobson BC, Van Dam J. et al. Fluorescence, reflectance, and light-scattering spectroscopy for evaluating dysplasia in patients with Barrett’s esophagus. Gastroenterology 2001; 120: 1620-1629
  • 18 Raica M, Cimpean AM, Ribatti D. Angiogenesis in pre-malignant conditions. Eur J Cancer 2009; 45: 1924-1934
  • 19 Hanahan D, Folkman J. Patterns and emerging mechanisms of the angiogenic switch during tumorigenesis. Cell 1996; 86: 353-364
  • 20 Kim C, Erpelding TN, Jankovic L. et al. Deeply penetrating in vivo photoacoustic imaging using a clinical ultrasound array system. Biomed Opt Express 2010; 1: 278-284
  • 21 Xu M, Wang LV. Photoacoustic imaging in biomedicine. Rev Sci Instrum 2006; 77: 41101
  • 22 Lim L, Streutker CJ, Marcon N. et al. Clinical study of ex vivo photoacoustic imaging in endoscopic mucosal resection tissues. Proc SPIE 2015; 9323: 932307
  • 23 Needles A, Heinmiller A, Sun J. et al. Development and initial application of a fully integrated photoacoustic micro-ultrasound system. IEEE Trans Ultrason Ferroelectr Freq Control 2013; 60: 888-897
  • 24 Conio M, Ponchon T, Blanchi S. et al. Endoscopic mucosal resection. Am J Gastroenterol 2006; 101: 653-663
  • 25 Lovell JF, Jin CS, Huynh E. et al. Porphysome nanovesicles generated by porphyrin bilayers for use as multimodal biophotonic contrast agents. Nat Mater 2011; 10: 324-332
  • 26 Ng KK, Shakiba M, Huynh E. et al. Stimuli-responsive photoacoustic nanoswitch for in vivo sensing applications. ACS Nano 2014; 8: 8363-8373
  • 27 Laser Institute of America. American national standard for safe use of lasers ANSI Z136.1—2014. Orlando: Laser Inst Am; 2014
  • 28 Laufer J, Elwell C, Delpy D. et al. In vitro measurements of absolute blood oxygen saturation using pulsed near-infrared photoacoustic spectroscopy: accuracy and resolution. Phys Med Biol 2005; 50: 4409-4428
  • 29 Wang X, Xie X, Ku G. et al. Noninvasive imaging of hemoglobin concentration and oxygenation in the rat brain using high-resolution photoacoustic tomography. J Biomed Opt 2006; 11: 24015
  • 30 Maslov K, Sivaramakrishnan M, Zhang HF. et al. Technical considerations in quantitative blood oxygenation measurement using photoacoustic microscopy in vivo. Proc SPIE, Photons Plus Ultrasound Imaging Sens 2006; 6086: 60860R1
  • 31 Cox BT, Arridge SR, Beard PC. Estimating chromophore distributions from multiwavelength photoacoustic images. J Opt Soc Am A Opt Image Sci Vis 2009; 26: 443-455
  • 32 Laufer J, Delpy D, Elwell C. et al. Quantitative spatially resolved measurement of tissue chromophore concentrations using photoacoustic spectroscopy: application to the measurement of blood oxygenation and haemoglobin concentration. Phys Med Biol 2007; 52: 141-168
  • 33 Lim L, Nichols B, Migden MR. et al. Clinical study of noninvasive in vivo melanoma and nonmelanoma skin cancers using multimodal spectral diagnosis. J Biomed Opt 2014; 19: 117003
  • 34 Kerkhof M, Van Dekken H, Steyerberg EW. et al. Grading of dysplasia in Barrett’s oesophagus: substantial interobserver variation between general and gastrointestinal pathologists. Histopathology 2007; 50: 920-927
  • 35 Yang J-M, Li C, Chen R. et al. Catheter-based photoacoustic endoscope. J Biomed Opt 2014; 19: 66001

Corresponding author

Liang Lim, PhD
University Health Network – Princess Margaret Cancer Centre
101 College Street
PMCRT #15-301V Toronto
Ontario M5G 1L7
Canada    

  • References

  • 1 Canada Cancer Society. Canadian Cancer Statistics 2014. Available from: http://www.cancer.ca/en/support-and-services/resources/publications/ Accessed 24 Feb 2015
  • 2 American Cancer Society. Cancer Facts & Figures 2015. Available from: https://www.cancer.org/research/cancer-facts-statistics/all-cancer-facts-figures/ Accessed 24 Feb 2015
  • 3 Solaymani-Dodaran M, Logan RFA, West J. et al. Risk of oesophageal cancer in Barrett’s oesophagus and gastro-oesophageal reflux. Gut 2004; 53: 1070-1074
  • 4 Shaheen NJ, Richter JE. Barrett’s oesophagus. Lancet 2009; 373: 850-861
  • 5 Schnell TG, Sontag SJ, Chejfec G. et al. Long-term nonsurgical management of Barrett’s esophagus with high-grade dysplasia. Gastroenterology 2001; 120: 1607-1619
  • 6 Reid BJ, Levine DS, Longton G. et al. Predictors of progression to cancer in Barrett’s esophagus: baseline histology and flow cytometry identify low- and high-risk patient subsets. Am J Gastroenterol 2000; 95: 1669-1676
  • 7 Vennalaganti PR, Naag KanakadandiV, Gross SA. et al. Inter-observer agreement among pathologists using wide-area transepithelial sampling with computer-assisted analysis in patients with Barrett’s esophagus. Am J Gastroenterol 2015; 110: 1257-1260
  • 8 Haisley KR, Dolan JP, Olson SB. et al. Sponge sampling with fluorescent in situ hybridization as a screening tool for the early detection of esophageal cancer. J Gastrointest Surg 2017; 21: 215-221
  • 9 Mannath J, Ragunath K. Role of endoscopy in early oesophageal cancer. Nat Publ Gr 2016; 13: 720-730
  • 10 Douplik A, Zanati S, Saiko G. et al. Diffuse reflectance spectroscopy in Barrett’s Esophagus: developing a large field-of-view screening method discriminating dysplasia from metaplasia. J Biophotonics 2014; 7: 304-311
  • 11 Kara M, DaCosta RS, Wilson BC. et al. Autofluorescence-based detection of early neoplasia in patients with Barrett’s esophagus. Dig Dis 2004; 22: 134-141
  • 12 Messmann H, Endlicher E, Freunek G. et al. Fluorescence endoscopy for the detection of low and high grade dysplasia in ulcerative colitis using systemic or local 5-aminolaevulinic acid sensitisation. Gut 2003; 52: 1003-1007
  • 13 Bergholt MS, Zheng W, Lin K. et al. In vivo diagnosis of esophageal cancer using image-guided Raman endoscopy and biomolecular modeling. Technol Cancer Res Treat 2011; 10: 103-112
  • 14 Yang VXD, Tang SJ, Gordon ML. et al. Endoscopic Doppler optical coherence tomography in the human GI tract: Initial experience. Gastrointest Endosc 2005; 61: 879-890
  • 15 Bouma B, Tearney G. High-resolution imaging of the human esophagus and stomach in vivo using optical coherence tomography. Gastrointest Endosc 2000; 51: 467-474
  • 16 Singh R, Shahzad MA, Tam W. et al. Preliminary feasibility study using a novel narrow-band imaging system with dual focus magnification capability in Barrett’s esophagus: Is the time ripe to abandon random biopsies?. Dig Endosc 2013; 25: 151-156
  • 17 Georgakoudi I, Jacobson BC, Van Dam J. et al. Fluorescence, reflectance, and light-scattering spectroscopy for evaluating dysplasia in patients with Barrett’s esophagus. Gastroenterology 2001; 120: 1620-1629
  • 18 Raica M, Cimpean AM, Ribatti D. Angiogenesis in pre-malignant conditions. Eur J Cancer 2009; 45: 1924-1934
  • 19 Hanahan D, Folkman J. Patterns and emerging mechanisms of the angiogenic switch during tumorigenesis. Cell 1996; 86: 353-364
  • 20 Kim C, Erpelding TN, Jankovic L. et al. Deeply penetrating in vivo photoacoustic imaging using a clinical ultrasound array system. Biomed Opt Express 2010; 1: 278-284
  • 21 Xu M, Wang LV. Photoacoustic imaging in biomedicine. Rev Sci Instrum 2006; 77: 41101
  • 22 Lim L, Streutker CJ, Marcon N. et al. Clinical study of ex vivo photoacoustic imaging in endoscopic mucosal resection tissues. Proc SPIE 2015; 9323: 932307
  • 23 Needles A, Heinmiller A, Sun J. et al. Development and initial application of a fully integrated photoacoustic micro-ultrasound system. IEEE Trans Ultrason Ferroelectr Freq Control 2013; 60: 888-897
  • 24 Conio M, Ponchon T, Blanchi S. et al. Endoscopic mucosal resection. Am J Gastroenterol 2006; 101: 653-663
  • 25 Lovell JF, Jin CS, Huynh E. et al. Porphysome nanovesicles generated by porphyrin bilayers for use as multimodal biophotonic contrast agents. Nat Mater 2011; 10: 324-332
  • 26 Ng KK, Shakiba M, Huynh E. et al. Stimuli-responsive photoacoustic nanoswitch for in vivo sensing applications. ACS Nano 2014; 8: 8363-8373
  • 27 Laser Institute of America. American national standard for safe use of lasers ANSI Z136.1—2014. Orlando: Laser Inst Am; 2014
  • 28 Laufer J, Elwell C, Delpy D. et al. In vitro measurements of absolute blood oxygen saturation using pulsed near-infrared photoacoustic spectroscopy: accuracy and resolution. Phys Med Biol 2005; 50: 4409-4428
  • 29 Wang X, Xie X, Ku G. et al. Noninvasive imaging of hemoglobin concentration and oxygenation in the rat brain using high-resolution photoacoustic tomography. J Biomed Opt 2006; 11: 24015
  • 30 Maslov K, Sivaramakrishnan M, Zhang HF. et al. Technical considerations in quantitative blood oxygenation measurement using photoacoustic microscopy in vivo. Proc SPIE, Photons Plus Ultrasound Imaging Sens 2006; 6086: 60860R1
  • 31 Cox BT, Arridge SR, Beard PC. Estimating chromophore distributions from multiwavelength photoacoustic images. J Opt Soc Am A Opt Image Sci Vis 2009; 26: 443-455
  • 32 Laufer J, Delpy D, Elwell C. et al. Quantitative spatially resolved measurement of tissue chromophore concentrations using photoacoustic spectroscopy: application to the measurement of blood oxygenation and haemoglobin concentration. Phys Med Biol 2007; 52: 141-168
  • 33 Lim L, Nichols B, Migden MR. et al. Clinical study of noninvasive in vivo melanoma and nonmelanoma skin cancers using multimodal spectral diagnosis. J Biomed Opt 2014; 19: 117003
  • 34 Kerkhof M, Van Dekken H, Steyerberg EW. et al. Grading of dysplasia in Barrett’s oesophagus: substantial interobserver variation between general and gastrointestinal pathologists. Histopathology 2007; 50: 920-927
  • 35 Yang J-M, Li C, Chen R. et al. Catheter-based photoacoustic endoscope. J Biomed Opt 2014; 19: 66001

Zoom Image
Fig. 1 Two examples of co-registered histopathology (hematoxylin and eosin staining) and ultrasound/photoacoustic images. a, c Histopathology overlaid with color-coded tissue classification (green: intestinal metaplasia [IM]; blue: squamous [SQ]; yellow: dysplasia [Dys]; gray: inflammation). The black lines are 2 mm scale bars. b, d Corresponding ultrasound (grayscale) and photoacoustic images (blue: 750 nm; red 850 nm). Regions of interest for different tissue classes are indicated. Scales are in mm.
Zoom Image
Fig. 2 Example of frequency histograms of photoacoustic imaging intensity for two of the regions of interest (ROIs) shown in [Fig. 1b], corresponding to the middle green ROI (intestinal metaplasia [IM]) and left blue ROI (squamous mucosa [SQ1). a 750 nm. b 850 nm. c Relative total hemoglobin.
Zoom Image
Fig. 3 3D reconstruction of the photoacoustic imaging (PAI) signal. a Ultrasound B-mode (grayscale). b PAI at 750 nm (blue), mainly dominated by hemoglobin absorption. c PAI at 850 nm (red), mainly dominated by oxy-hemoglobin absorption. Axis scale in mm.