1.

Introduction

Prostate cancer is the most prevalent cancer in the male population and the second most common cancer-related cause of death.¹ When elevated serum prostate specific antigen (PSA) concentration has raised suspicion, prostate biopsies are the gold standard to confirm the diagnosis of cancer. These biopsies are evaluated by a pathologist and the development of cancer is graded using the Gleason Score, which ranges from 1, unaggressive tumor (well-organized structure with low cellular density), to 5, aggressive tumor (irregular structures with high cellular density).²

The discovery of PSA made the incidence of prostate cancer increase drastically whereas cancer-related mortality remained unchanged.^1,3 Concerns have been raised about detection and overtreatment of these low-risk prostate cancers, since treatment-related side effects (erectile dysfunction and urinary incontinence) may impair quality of life significantly.⁴ For these low- to intermediate-risk patients, focal therapy can be a favorable disease management option.⁵ This form of therapy aims to eradicate known cancer sites, while leaving healthy tissue untouched, which has the potential benefit of diminishing side effects compared to radical treatments and, thus, maintaining quality of life.⁴ However, for focal therapy, accurate identification, localization, demarcation, and grading of a lesion are essential. Imaging modalities such as magnetic resonance imaging (MRI) and transrectal ultrasound (TRUS) are insufficiently accurate to localize prostate tumors with a reported sensitivity of 0.87 and a specificity of 0.30.⁶ Another modality for prostate cancer localization is transperineal three-dimensional (3-D) prostate mapping biopsy. The prostate is sampled every 0.5 cm and the coordinates are correlated to tumor location. The limitations of mapping biopsies are the low volume of tissue sampled in a biopsy and, consequently, the large number of cores needed (average of 50 per prostate), resulting in a larger number of complications and physician-related delay of the results (pathologist has to evaluate all the biopsies), which prevent treatment from occurring in the same session as diagnosis.^7,8 The use of fusion technologies integrating MRI and TRUS to target the biopsy location has been shown to improve on these approaches.⁹ Yet, targeted biopsies still rely largely on standard H&E pathology, which is analyzed at a later stage. If targeted biopsy is combined with a real-time technology which can visualize and analyze tissue structure and architecture similar to H&E histology, simultaneous diagnosis and treatment would be feasible.

Optical coherence tomography (OCT) has the potential to fulfill the premise of a fully digital targeted biopsy. OCT is an in vivo imaging modality, similar to ultrasonography, which allows real-time microscopic imaging and quantitative analysis of the backscattered light of an imaged sample. Using near-infrared light (1300 nm), OCT results in very detailed images with a resolution that can reach $10\mu \mathrm{m}$, although the scattering limits the imaging depth to $\sim 2\mathrm{mm}$.¹⁰ Technological advancements have led to novel applications of OCT in urology,¹¹12.^–13 gastroenterology,¹⁴15.^–16 gynecology,^17,18 dermatology,¹⁹20.^–21 cardiology,²²23.^–24 and neonatology.²⁵ The contrast in an OCT image is based on differences in light scattering between different cellular structures in a sample. The resulting exponential decay of the OCT signal with depth, which is parameterized by the attenuation coefficient [${\mu }_{\mathrm{oct}}({\mathrm{mm}}^{-1})$], directly relates to the scattering properties of the tissue under study.²⁶ During carcinogenesis, a cell is subjected to cellular changes which alter the light scattering properties. Quantification of OCT images using ${\mu }_{\mathrm{oct}}$ measured from lesions in the ureter,¹² kidney,¹¹ vulva,¹⁸ oral tissue,²⁷ and lymph node metastases²⁸ confirmed the ability of OCT to distinguish between tissue types based on ${\mu }_{\mathrm{oct}}$. We hypothesize that the attenuation coefficient also correlates with cellular properties of prostate tissue. To the best of our knowledge, this property was never demonstrated in the prostate before.

This study aims to demonstrate the feasibility of needle-based OCT and functional analysis of OCT data along the full pullback measurement trajectory of the OCT. This also includes validation of the OCT by system and probe calibration. We describe the first in-tissue needle-based OCT imaging and three-dimensional optical attenuation analysis of a single prostate after radical prostatectomy, which we compared with histopathological diagnosis. Our work aligns with stages 1 and 2a of the IDEAL framework of clinical technological innovations.^29,30

2.

Materials and Methods

2.2.

Calibration and Three-Dimensional Attenuation Fitting

Signal attenuation is influenced by the technical properties of the OCT system itself, warranting careful calibration. The rotating imaging probes (like a lighthouse) for this system are provided in sterile packaging so that a priori calibration is not possible.

The attenuation coefficient (${\mu }_{\mathrm{oct}}$) was determined as described before.^18,26 Briefly, the decay of light intensity with depth (${\mathrm{mm}}^{-1}$) is quantified by fitting the OCT data to a single exponential decay using

(1)

$$i(z)\propto \sqrt{[{\mathrm{e}}^{-2{\mu }_{\mathrm{oct}}z}]}+y_0,$$

where $2z$ is the axial round-trip path length of the light in the sample and ${\mu }_{\mathrm{oct}}$ is the attenuation coefficient; $y_0$ is a term accounting for the fraction of noise offset. The square root accounts for the fact that the detector current $i(z)$ is proportional to the field returning from the sample, rather than intensity.²⁶ Prior to fitting, the signal is corrected for system-induced attenuation, i.e., reduction in amplitude with increasing distance to the focal point (quantified through the confocal point spread function³¹ and reduction in amplitude with increasing distance from the zero-delay point, the so-called sensitivity roll-off due to finite spectral line width of the swept source and finite integration time of the detector³²).

The ability of the St. Jude OCT system to measure correct values of ${\mu }_{\mathrm{oct}}$ was validated by measuring ${\mu }_{\mathrm{oct}}$ for increasing concentrations of a scattering medium (Intralipid®) as our group described previously.^33,34 The OCT system induced signal attenuation (due to the combined effect of the confocal point spread function and the sensitivity roll-off, see below) was calibrated on a highly diluted Intralipid sample (0.0005 volume%). After correction for water absorption ($\sim 0.2{\mathrm{mm}}^{-1}$), an attenuation of $1.1{\mathrm{mm}}^{-1}$ remained, which is subtracted from subsequently measured ${\mu }_{\mathrm{oct}}$ values. Validation of ${\mu }_{\mathrm{oct}}$ measurements was performed by determining the attenuation coefficient from Intralipid dilutions with increasing scattering coefficient. To ensure that each single use sterile imaging probe is comparable to the measurements of the other imaging probes, the probe-to-probe variations were analyzed by measuring attenuation coefficients of 0.03, 6.5, 11.4, 17, and 22.7 volume% Intralipid®. Experiments were repeated five times per probe for four probes.

In previous studies, we concluded that attenuation coefficients in tissue can reliably be determined on depth segments down to $50\mu \mathrm{m}$ length (provided 50 to 100 A-scans can be averaged to suppress signal variation due to speckle).³⁵ Recently, algorithms have been presented that allow determination of the attenuation coefficient down to the pixel level.³⁶ Since our goal is to localize tumor tissue along the trajectory of the biopsy needle, we devised the procedure outlined below to reduce the amount of data.

Visual inspection of the OCT images reveals that the data appear largely homogeneous along the pullback, so that fitting regions can be set at 50 data points in the automated analysis. To practically implement the attenuation analysis on a full 3-D pullback dataset acquired with the St. Jude OCT console, a custom-made plugin was developed for ImageJ.³⁷ The ${\mu }_{\mathrm{oct}}$ analysis is performed on the raw data (Fig. 1). First, to optimize the fitting procedure, the data were arranged in 31 discrete radial wedges (pie-slices) per circular segment. Each wedge consists of 16 radial A-scans per wedge and 5 axial B-scans [Fig. 1(b)]. The resulting average wedge is further smoothed using a Savitsky-Golay filter with a width of 21 data points [Fig. 1(c)].³⁸ Second, the first 90 data points in an A-scan were removed since they contain only inner reflections of the probe itself [Fig. 1(d)]. Third, the decay of the OCT signal is analyzed starting from the first maximum data point until 50 points ($500\mu \mathrm{m}$) in depth, using Eq. (1) [Fig. 1(e)]. Fourth, subtraction of the ${\mu }_{\mathrm{oct}}$ of the calibration measurement ($1.1{\mathrm{mm}}^{-1}$) from the signal decay yields the attenuation coefficient. Finally, the OCT data were visually inspected for extremely high or low scattering structures to account for extraordinary high or low attenuation coefficients created by structures like calcifications, cysts or air gaps between tissue and probe. These values were excluded from the correlation analysis with pathology. For visualization, the data were grouped with a corresponding color representing a specific ${\mu }_{\mathrm{oct}}$ range according to their attenuation coefficient; ${\mu }_{\mathrm{oct}}=0$ to $5.5{\mathrm{mm}}^{-1}=\log$ java hue, 0.65 to 1.00, ${\mu }_{\mathrm{oct}}>5.5{\mathrm{mm}}^{-1}$ and $<0.5=\text{gray}$ (Fig. 2).

Fig. 1

(a) Raw optical coherence tomography (OCT) data consist of 1024 radially directed A-scans per B-scan and 541 B-scans in the $z$ direction, covering $541*100\mu \mathrm{m}=54.1\mathrm{mm}$. (b) 16 A-scans and 5 B-scans were averaged to smoothen the data, resulting in 31 radial averaged A-scans and 108 B-scans over the 54.1 mm length. (c) Savitzky-Golay filter was applied to smooth the A-lines. (d) The first 90 data points were removed from each averaged A-scan to remove scattering properties from the probe itself. (e) Attenuation analysis was performed over first 50 points after the maximum of the resulting curve (red part).

Fig. 2

The raw data were plotted (upper part). With the help of a customized Image J plug-in, the attenuation coefficient was analyzed (middle part). Each range of attenuation coefficients was given a color for interpretation corresponding to the scale bar on the left side (lower part).

2.3.

Measurement of Prostate Tissue

Directly after radical resection, the prostate was transported to the pathology department for ex vivo measurements. First, a silicon catheter was placed to indicate the urethra. Six intravenous (IV) catheters [Terumo Surflo® $18\mathrm{G}\times 2(1/2)$] were introduced in the prostate, two in the transitional zone and four in the peripheral zone. The needle was removed from the catheter and the OCT probe was inserted in the IV catheter. During measurements, the IV catheter was retracted to ensure direct contact of the C7 Dragonfly™ OCT imaging probe with prostate tissue. Measurements (OCT pullbacks) started at the apex and extended to the base of the prostate. After data collection, the catheters were reintroduced over the OCT imaging probe and left behind in the tissue to mark the imaging trajectory for correlation with the pathological diagnosis (Fig. 3).

Fig. 3

Ex vivo OCT measurement directly after radical prostatectomy. (a) Performance of the measurement with the 0.9 mm C7 Dragonfly™ OCT imaging probe positioned in the tissue. Measurement pullbacks start at the apex (left on the picture) and end at the base of the prostate (right in the picture). (b) Marking the trajectory by replacing the intravenous (IV) catheters in the tissue after measurements for optimal correlation with whole mount histopathology.

3.

Results

The calibration measurements in various Intralipid concentrations are depicted in Fig. 4, demonstrating a nonlinear increase of attenuation coefficient with increasing concentration. These values and nonlinear behavior, which is due to dependent and multiple scattering, are in concordance with our earlier published results.^33,34 As shown in Fig. 5 and Table 1, the smallest difference in ${\mu }_{\mathrm{oct}}$ per probe is $0.03{\mathrm{mm}}^{-1}$; the largest difference in ${\mu }_{\mathrm{oct}}$ per probe is $0.33{\mathrm{mm}}^{-1}$. The smallest difference in ${\mu }_{\mathrm{oct}}$ between probes is $0.15{\mathrm{mm}}^{-1}$; the largest difference in ${\mu }_{\mathrm{oct}}$ between probes is $0.65{\mathrm{mm}}^{-1}$. For Intralipid concentrations within the scattering range of prostate tissue ($3-6{\mathrm{mm}}^{-1}$) two times the standard deviation of the measurements, an indication of the precision (reproducibility) of our measurement technique, was $<8\%$ of the mean of the measured values.

Table 1

Mean optical attenuation coefficients per probe with range (visual representation in Figs. 1 and 2). Note that per probe the smallest and largest ranges in μoct are 0.06 and 0.33  mm−1, respectively. Between probes, the smallest and largest differences in measured μoct are 0.15 and 0.65  mm−1, respectively.

Fig. 4

Attenuation measurements were validated by measurement of ${\mu }_{\mathrm{OCT}}$ of samples with increasing concentrations of Intralipid®. The boxplots represent mean and range for all probes combined (see also Table 1).

Fig. 5

Interprobe variability was tested by measuring four probes five times in increasing concentrations of Intralipid. The boxplots represent the mean and range (see also Table 1).

Analysis of the OCT data obtained in the prostate resulted in mean attenuation coefficients of $3.8{\mathrm{mm}}^{-1}$ for parts of the tissue that were classified as benign (SD: $0.8{\mathrm{mm}}^{-1}$, minimum: $2.2{\mathrm{mm}}^{-1}$, maximum: $8.9{\mathrm{mm}}^{-1}$) and $4.1{\mathrm{mm}}^{-1}$ for parts of tissue that were classified as malignant (SD: $1.2{\mathrm{mm}}^{-1}$, minimum: $2.5{\mathrm{mm}}^{-1}$, maximum: $9.0{\mathrm{mm}}^{-1}$).

Figure 6(a) shows a 3-D representation of the OCT measurement with visualization of three separate B-scans. In Fig. 6(b), which was characterized as malignant tissue on pathology, glandular appearing structures can be observed. In Fig. 6(c), which was characterized as benign tissue on pathology, homogeneous prostate tissue can be seen. In Fig. 6(d), a cyst filled with clear fluid is shown, which the automated analysis recognizes as a high attenuation coefficient. The ring around the B-scans visually represents the attenuation coefficient (the same color scale is used as in Fig. 2). It is clear that the attenuation coefficient differs in prostate tissue.

Fig. 6

Structural properties of prostate tissue. (a) Three-dimensional (3-D) representation of an OCT measurement of prostate tissue with a visualization of three separate B-scans. (b) B-scan from malignant area (Gleason $3+4$) of the prostate. Glandular appearing structures are observed in these malignant areas (arrow). Optical attenuation coefficients in this area are high, indicated by the red ring around the raw data (visualizing the output of the automated attenuation analysis). (c) B-scan from a benign region of the prostate. In these benign areas, homogeneous tissue was seen, resulting in a low attenuation coefficient, indicated by the blue ring (visualizing the output of the automated attenuation analysis). (d) B-scan from an area in which a cyst was seen. In these irregular areas, attenuation coefficient matching is not always accurate, resulting in a high variation of attenuation coefficients, indicated by the multiple colors on the ring (visualizing the output of the automated attenuation analysis). The scale bar is 5 mm.

Figure 7 shows a prostate after histopathological evaluation. After fixation in formalin, the prostate was painted orange on the left side and blue on the right side. Subsequently, the prostate was sliced from base (slice 1) to apex (slice 8), and all slices were positioned with the same side up. The OCT measurement trajectories, indicated by the IV catheters, are very well visible. Four H&E stained microscopy slides per prostate slice were obtained and the pathologist indicated areas of tumor on these slides, as marked in red. These microscopic slides were overlaid on the overview picture of the prostate slices. Using AMIRA™, we reconstructed the prostate contour and lesion contours by stacking the slices in 3-D [Fig. 8(a)]. By plotting the St. Jude OCT measurements [Fig. 8(c)] and attenuation maps [Fig. 8(d)] in these reconstructed contours by using the holes of the IV catheters as a guide, we could estimate the distance between the slices (average slice thickness). Also, we were able to estimate which part of the OCT measurements went through tumor and which parts went through benign tissue. Finally, the information on coregistration was used to determine which B-scans corresponded to the sampled areas in the prostate slices. The attenuation information from these matched B-scans was also plotted in the picture of the prostate slices (Fig. 7, 1 to 8). Not all areas in which a high attenuation coefficient became apparent corresponded to areas of prostate cancer.

Fig. 7

Histopathological evaluation. The prostate was sliced from base (slice 1), to apex (slice 8), and all slices were positioned with the same side up. The pathologist indicated areas of tumor on the slides, as marked in red. The attenuation information from these matched B-scans were also plotted in the picture of the prostate slices (1 to 8). # corresponds to Fig. 6(b), * to Fig. 6(c), and ^ to Fig. 6(d).

Fig. 8

3-D correlation to histopathology. (a) Using AMIRA™, we reconstructed the prostate contour (green) by stacking the slices from Fig. 6 in 3-D. (b) Because the pathologist outlined the tumor contours on the slides, we were able to reconstruct the tumors in 3-D. By plotting the (c) St. Jude OCT measurements and (d) attenuation maps in these reconstructed contours, based on the IV catheters that were visible in the slices, we could estimate which B-scans corresponded to the sampled areas in the prostate slices. The attenuation maps of these corresponding B-scans are overlaid in Fig. 7.

4.

Discussion

We demonstrate a method that allows for 3-D quantitative analysis of OCT optical attenuation coefficient datasets that may indicate a correlation with tissue disease status. This study is the first to show the feasibility of quantitative needle-based OCT in the prostate. This study is the first step toward real-time objective diagnosis of prostate cancer, a largely emphasized challenge in urology.

Full clinical translation of OCT should build on three fundamental pillars. First, qualitative properties should be objectified with quantitative information, which can be spatial measurements from images (e.g., layer thicknesses), attenuation coefficients such as in our study, or more complex statistics related to tissue organization. A major challenge is relating measured optical properties to gold standard pathology diagnosis. This gap cannot be bridged without fundamental studies—far beyond the scope of our present contribution—that include both advanced OCT signal modeling and leveraging the potential of quantification of digital pathology.

Second, the technology should be compatible with existing procedures and protocols. We presently use the clinically proven OCT equipment C7-XR™ OCT Intravascular Imaging System interfaced to a C7 Dragonfly™ Intravascular Imaging Probe (St. Jude Medical, St. Paul, Minnesota). Since the sterile probes are unpackaged just before measurements, prior calibration of the probes is not possible. We, therefore, devised an efficient, cost-effective procedure based on an easily obtained Intralipid™ suspension. Moreover, this study shows that intra- and interprobe variation is minimal, which suggests high data fidelity even if post hoc calibration is not possible (for example, because blood has entered the imaging catheter). However, the study shows that postmeasurement calibration of the probes improves measurement precision, since there is some small variation between the probes.

Third, the future of prostate cancer will be image guided targeted diagnosis, likely by a combination of imaging technologies.³⁹ Using TRUS, a location estimation of the OCT probe can be obtained. In addition to this, OCT visualizes the position of a lesion along the optical biopsy axis. It has been shown that when MRI data were fused with data from a conventional TRUS, the sensitivity increased drastically. More biopsies were found positive in the fused group and more malignant tissue was found per biopsy,⁴⁰ even accomplishing results similar to transperineal 3-D prostate mapping biopsy.⁴¹ Our hypothesis is that integrating OCT in the combined results of MRI/TRUS fusion will further improve the diagnostic accuracy. Moreover, results will become objective and real time. For the diagnosis of kidney tumors, similar developments are ongoing. Projects are running to test the ability of OCT as a means of optical biopsy for kidney cancer (NCT02073110, Ref. 42) using the same OCT device as is used in our study. The few studies performed regarding OCT in the prostate focus on the qualitative interpretation of optical findings to identify surgical margins and neurovascular bundles.⁴³44.45.^–46 One manuscript described a difference between malignant and benign prostate tissue structures visualized with OCT; however, these results were not quantitatively analyzed.⁴⁷ When the OCT technology further improves (smaller, higher-resolution, faster machines) and analysis software evolves (e.g., structure recognition, automatic learning, etc.), the technology will most likely become faster, more objective and more accurate.

OCT is a new diagnostic modality in the prostate. We acknowledge the limitations that this study entails. First, the St. Jude system evaluates tissue every 0.1 mm, whereas whole mount pathology assesses the tissue in theory every 4 mm. Due to free hand slicing, this slice thickness varies in practice between 3 and 5 mm, creating $\sim 20$ OCT images uncertainty in the first slice. This matching-uncertainty is of non-negligible proportion, since every prostate slice has this variable slice thickness of 3 to 5 mm. Furthermore, we assumed in this study that the slice thickness is constant throughout the prostate slice. However, it is very well possible that a slice might be slightly wedge shaped in reality. Also, when a tumor is not present throughout the whole slice, it can cause matching issues (Fig. 9). All these aspects create uncertainty in the 3-D OCT histopathology matching and have to be overcome for further validation studies. Solving this issue is currently in progress by using a customized tool for dedicated pathology matching and slicing. Second, the St. Jude imaging probe does not have a marker that indicates angular probe position in tissue, so a method should be designed to register this. Further studies are in progress that provide larger numbers of patients and address the slicing and pathology-matching issues described above. Finally, in ex vivo measurements, blood flow and tissue perfusion are not present. In vivo measurements are needed to determine whether or not the results are reproducible.

Fig. 9

Schematic representation of histopathology matching challenges. Each blue block (numbered 1 to 8) represents a slice of prostate. The bold red lines on each left side of a blue block represent the pathology slide (part of the prostate that is actually visualized under the microscope). The striped bar in the middle represents an OCT scan, and each stripe represents a B-scan. The letters on the left side represent the possible scenarios. A, slice of prostate that looks benign but actually is half tumor (false negative); B, prostate slice that looks malignant and is malignant (true positive); C, prostate slice that looks malignant but is half benign (false positive); and D, prostate slice that looks negative and is negative (true negative). It is clear that because of the matching issues, the difference between benign and malignant attenuation coefficients becomes smaller.

5.

Conclusion

We demonstrated the feasibility of needle-based 3-D quantitative 1300 nm OCT in prostate tissue as a first step toward objective and real-time digital diagnosis of prostate cancer. Fully automated attenuation coefficient analysis was performed over the full pullback. Optimal correlation with pathology was achieved by coregistration of 3-D OCT attenuation maps with 3-D pathology of the prostate. This approach may contribute to the current challenge of prostate imaging and the rising interest in focal therapy for reduction of side effects occurring with current therapies. However, in further research, the challenges of exact histopathology correlation as well as ${\mu }_{\mathrm{oct}}$ analysis of different cell types of the prostate need to be addressed.