Introduction

One of the most challenging problems in clinical dermatology is the early detection of melanoma [4, 42, 46]. Since clinical diagnosis may be difficult, non-invasive imaging techniques have been developed to enhance early diagnosis in challenging cases [10, 17]. Dermoscopy and reflectance confocal microscopy (RCM) are both able to considerably improve the diagnostic accuracy for melanoma, especially when used synergistically [11, 32, 33, 48]. However, diagnosis strongly depends on the experience of physicians [2, 31].

New diagnostic tools providing automated classification of pigmented skin lesions usable by non-experts have been proposed [23]. Spectral methods fall into this class of emerging new techniques holding the promise to provide quantitative criteria for melanoma diagnosis and to improve early diagnosis [37]. Multispectral information can be assessed both in the spatial domain (multispectral digital dermoscopy [22]) or in the frequency domain (spectroscopic methods such as diffuse-reflectance spectroscopy [21, 41, 54], Raman spectroscopy [36, 52] and fluorescence spectroscopy [15, 34]).

These methods evaluate the colour of a lesion by measuring its reflectance as a function of the wavelength. A comprehensive analysis of pigmented lesions under both ultraviolet and infrared radiations has previously been described in 1976 [37]. It was stated that infrared photographs tends to discriminate the different types of melanocytic lesions, with melanoma showing a relatively high degree of correlation with low infrared reflectance [37].

RCM and optical coherence tomography (OCT) are both technologies operating in the near-infrared (NIR) optical window. Conventional OCT enables the non-invasive imaging of structures with a non-cellular resolution (7.5–15 µm) up to a depth of about 1 mm, while high-definition OCT (HD-OCT) provides cellular (3 µm) resolution up to a depth of 570 µm [3, 5, 6, 8, 18, 19, 25, 44]. Conventional OCT has been used for non-invasive tumour thickness measurement in melanocytic skin lesions. A variable correlation with histopathology could be observed depending on OCT technology used [24, 38]. HD-OCT offers additional structural and cellular information on melanocytic lesions complementary to that of RCM [9]. However, the diagnostic potential of HD-OCT seems to be not high enough for ruling out the diagnosis of melanoma if based on morphology analysis [20].

In a recent study, it was demonstrated that HD-OCT permits to quantify the reflectance of NIR light in a volume of skin as function of depth [7]. Therefore, the aim of this paper was to quantify in vivo optical properties such as light attenuation in melanocytic lesions by HD-OCT. An additional objective was to determine the best critical value of these optical properties for melanoma diagnosis.

Methods

Study sample

The cases for this retrospective pilot study were retrieved from patient files collected at first author’s private practice between 2012 and 2015. The cases contain melanocytic lesions including benign nevi (BN), dysplastic nevi (DN) and melanomas (MM). The inclusion criteria were: (i) selection of clinically suspected melanocytic lesions for full excision based on dermoscopic or combined dermoscopic/RCM evaluation and subsequent histopathologic diagnosis and (ii) availability of relevant (see procedure) HD-OCT images taken prior to excision.

Approval from local ethical committee was obtained (P2015/301). All patients provided informed consent for imaging their lesion by HD-OCT (Skintell®, Agfa Healthcare, Mortsel, Belgium). We conformed to the Helsinki Declaration with respect to human subjects in biomedical research. All international rules governing clinical investigation of human subjects were strictly followed. This study affected neither the routine diagnosis nor the treatment of the lesions presented by the included subjects.

Histopathology

Histopathologic analysis was carried out by two board-certified histopathologists who classified each lesion according to histopathological criteria for melanocytic lesions described elsewhere [12–14].

High-definition optical coherence tomography

Instruments and acquisition methods and technical details have been previously described [5, 9]. For the purpose of this study we will remind what follows.

HD-OCT has four basic design principles: (i) a time domain OCT providing depth information by the position of the reference mirror, (ii) full-field illumination by a light source enabling very low lateral coherence and permitting a high-speed 3-D sharp image acquisition, (iii) high-power tungsten lamp with Gaussian filter and ultrahigh bandwidth (1300 ± 100 nm). This permits high-depth resolution of 3 µm and (iv) full-field domain OCT system with dynamic focus tracking: the focal plane is continuously moved through the skin sample. This ensures a high lateral resolution over the entire depth of 3 µm. HD-OCT offers a constant homogeneous resolution of 3 µm in all three dimensions. The system is capable of capturing a cross-sectional and en face image in real time, as well as of fast real time 3-D acquisition. A sharp image at all depth is guaranteed up to a theoretical depth of 570 µm. The field of view is 1.8 × 1.5 mm. The total light power at the tissue is <3.5 mW.

HD-OCT focus tracking

Skin is a highly turbid medium. Light which propagates through skin is scattered and absorbed resulting in attenuation of this light. In contradiction to a “normal” light detector HD-OCT is only detecting photons which reach the detector on an almost straight path—the ballistic portion of the photons or photons which were only forward scattered a few times [49].

The resulting attenuation coefficient was estimated by Jacques et al. [27–29]. Since absorption is much smaller than scattering for tissues in the near-infrared spectrum, absorption can be neglected [29]. HD-OCT operates in the second diagnostic window of the NIR spectrum [47] and hence the main attenuation mechanism for HD-OCT is scattering. Measuring tissue optical properties in vivo is only applicable to OCT when operating focus tracking mode which is the case for HD-OCT [27]. Table 1 provides definitions and abbreviations of optical properties adapted from the literature [27, 29, 30, 35, 49, 50].

Table 1 Definition of optical properties adapted from the literature [33, 35, 36, 38–40]
Full size table

Procedure (Figs. 1, 2, 3)

The selection of spots in melanocytic lesions imaged by HD-OCT was based on dermoscopic or combined dermoscopic/RCM evaluation. Subsequently the selected spots were marked with a plastic ring of 2 mm diameter. A 3-D HD-OCT DICOM (digital imaging and communication in medicine) image taken at the centre of this ring was selected. This DICOM file was then opened using ImageJ® software. A single square region of interest (ROI) in the en face view (green square) was selected in each melanocytic lesion based on the presence of relevant morphologic features as described in previous papers [9, 20]. To exclude obvious “artefacts” such as air bubbles or skin furrows, a square ROI was carefully chosen. The size of the ROI was set at least to 300 × 300 µm.

Fig. 1
figure 1

Benign nevus. a, b A 3-D HD-OCT DICOM (digital imaging and communication in medicine) image of each melanocytic lesion is selected. To exclude artefacts such as skin furrows or air bubbles, a square region of interest is chosen in the en face image based on the presence of relevant morphologic features of melanocytic lesions (green square). A plot z-axis profile of the scanned volume is performed. c, d The obtained graph displays the reflectance of the ballistic photons [OCT signal: measured on the y-axis with arbitral units (AU)] versus imaging depth which is indicated on the x-axis (slice numbers from 0 to 200 axial position of depth). The slice numbers are multiplied by factor 3 (slice thickness equals 3 µm) to correspond to the real depth in micron (from 0 to 600 µm: see x-axis in (e). e Semi-log plot: an exponential function becomes a straight line given by y = ax + b. Four successive layers (epidermis, upper papillary dermis, deeper papillary dermis and superficial reticular dermis) with clear exponential decay are identified and plotted. A straight line is fitted in each of the four layers (i = 1–4) resulting in equation of the type y i  = ax i  + b whereby a is proportional to the attenuation coefficient for each of the four layers

Full size image
Fig. 2
figure 2

Dysplastic nevus. a Cross-sectional image, b en face image, c section of the cross-sectional image, d graph displaying reflectance versus depth of focus and e semi-log plot

Full size image
Fig. 3
figure 3

Malignant melanoma. a Cross-sectional image, b en face image, c section of the cross-sectional image, d graph displaying reflectance versus depth of focus and e semi-log plot

Full size image

The obtained graph displayed the reflectance and attenuation of the ballistic photons [OCT signal: measured on the y-axis with arbitral units (AU)] versus imaging depth which was indicated on the x-axis (from 0 to 200 axial position of depth). In a next step, the offset corresponding to the mean signal within gel was removed from each OCT signal. The first peak corresponded to the skin entrance signal (SES). All the HD-OCT values will be divided by SES to normalize them to the SES. The natural logarithm ln(x) of those values will be taken. The slice numbers needed to be multiplied by factor 3 to correspond to the real depth in micron (from 0 to 600 µm).

The decay with imaging depth of reflectance (backscattered intensity) of ballistic photons is a process that can be well modelled by an exponential relationship [26, 29].

When operating in the NIR diagnostic window, absorption becomes negligible and the slope of the exponential attenuation is proportional to the reduced scattering coefficient µ′s. Analysis by semi-log plot (an exponential function becomes a straight line given by y = ax + b) of the exponential profile of light attenuation, can provide information on tissue scattering properties [26, 29].

Four successive layers (epidermis, upper papillary dermis, deeper papillary dermis and superficial reticular dermis) with clear exponential decay were identified and plotted. A straight line was fitted in each of the four layers (i = 1–4) resulting in equation of the type y i  = ax i  + b whereby a was proportional to the attenuation coefficient for each of the four layers given by µ = a/2.

Three optical properties were measured:

  • Relative attenuation factor normalized to skin entrance signal for each of the four layers (µ raf1-4; µm−1)

  • Skin entrance signal (SES; arbitrary unit)

  • z value under the skin surface at which the OCT signal equals half of the SES, half value layer (z 1/2; µm)

Statistical analysis

One-way analysis of variance (ANOVA) was used to compare means of three samples (BN,DN and MM) using the F distribution. Prior to the Anova test, Levene’s Test for “Equality of Variances” was performed. If the Levene’s test was positive (p < 0.05) then the variances in the different groups were different (the groups are not homogeneous) and a logarithmic transformation to the data has been performed. Moreover, Scheffé test was used for all pairwise comparisons.

The best critical value of all HD-OCT assessed optical properties was determined by applying the receiver operating characteristic curves (ROC curves). This is a graph displaying the relationship between the true-positive rate (on the vertical axis) and the false-positive rate (on the horizontal axis) (Table 2).

Table 2 In vivo HD-OCT measured values of optical properties of melanocytic lesions according to thee subgroups
Full size table

Based on these critical values absolute and relative frequencies were calculated for the three groups. Chi-squared (χ2) test was employed to compare each melanocytic group versus the other melanocytic groups. The phi (φ) coefficient, employed to weight diagnostic power of each significant parameter, is a measure of association of two binary variables and is related to the Chi-squared (χ2) statistic by the formula: φ 2 = χ2/n, where n equals the total number of observations.

Calculations were made using MedCalc® statistical software version 14.12.0.

Results

Subjects

We included in this retrospective study 45 cases. These cases were retrieved from 2249 patient files collected between 2012 and 2015. Based on dermoscopic or dermoscopic/RCM assessment, 321 cases underwent full excision for histopathologic diagnosis. Forty-nine of these cases have been imaged by HD-OCT whereof 45 cases were retained based on the presence of relevant morphologic features of melanocytic lesions as described previously [9, 20]. The cases comprise 15 BN, 15 DN and 15 melanomas. All melanomas were histopathologically classified as superficial spreading melanomas with a Breslow-index varying between 0.31 and 0.72. The lesions belonged to 25 females and 20 males with skin type I–III and ages ranging from 25 to 70 years (median 51 years).

Quantitative in vivo evaluation (Table 3)

Table 3 Quantitative evaluation of optical properties of melanocytic lesions
Full size table
  1. 1.

    µ raf1: a significant difference (p < 0.001) could be observed between the three subgroups in the upper layer. µ raf1 decreased progressively from benign to malignant lesions. No significant differences could be demonstrated between the three groups for the other layers µ raf2-4. (Fig. 4).

    Fig. 4
    figure 4

    Comparison of trend lines of relative attenuation factor at different skin layers in benign nevus (upper part), dysplastic nevus (middle part) and melanoma (lower part)

    Full size image
  2. 2.

    SES: a significant (p < 0.001) reduction in skin entrance signal could be observed in the melanoma group compared to the two other groups.

  3. 3.

    z 1/2: a significant difference (p < 0.001) could be observed between the three subgroups. This value increased progressively with malignancy.

Best critical value of optical properties

Critical values permitting differentiation by HD-OCT of the three subgroups could be defined by applying the receiver operating characteristic curves (ROC curves) and are displayed in Table 4.

Table 4 Selection of best critical value of optical properties of melanocytic lesions by applying receiver operating characteristic curves
Full size table

Absolute and relative frequencies of optical properties (Table 5)

Table 5 Absolute and relative frequencies of optical properties of melanocytic lesions in relation to critical values as assessed in vivo by high-definition optical coherence tomography
Full size table
  1. 1.

    Differentiation of MM from non-malignant lesions:

    • µ raf1 < 0.0067/µm was a high sensitive (SS) and specific (SP) feature of MM group (93.3 and 96.7 %, respectively; ϕ = 0.90 and χ 2 = 36.45, p < 0.001). Moreover, the negative predictive value (NPV) was 96.7 % and the positive predictive value (PPV) was 93.3 %. The diagnostic accuracy was 95.6 %.

    • SES < 475 a.u. was a sensitive and specific feature of MM group (80.0 and 83.3 %, respectively, NPV = 89.3 %, PPV = 70.6 %, ϕ = 0.62 and χ 2 = 17.06, p < 0.001).

    • z 1/2 > 8.18 µm was a high sensitive and specific feature of MM group (93.3 and 86.7 %, respectively, NPV = 96.3 %, PPV = 77.8 %, ϕ = 0.77 and χ 2 = 26.67, p < 0.001).

  2. 2.

    Differentiation of benign nevi from non-benign nevi (DN and MM)

    • µ raf1 > 0.0248/µm was a sensitive and high specific feature of the BN group (80.0 and 90.0 %, respectively, NPV = 93.0 %, PPV = 81.0 %, ϕ = 0.60, χ 2 = 16.25, p < 0.001). Values higher than this cutoff value were absent in all lesions of the MM group.

    • SES > 678 a.u. was moderate sensitive and specific for the BN group (80.0 and 76.7 %, respectively, NPV = 88.0 % and PPV = 63.0 %, ϕ = 0.54, χ 2 = 13.16, p < 0.01). In MM no value higher than 678 a.u. could be observed.

    • z 1/2 < 6.64 µm was a moderate sensitive and specific feature of the BN group (both 73.3 %, NPV = 85.0 % and PPV = 58.0 %, ϕ = 0.45 and χ 2 = 8.93, p < 0.01). Values lower than this cutoff value were absent in all lesions of the MMgroup.

Discussion

Non-invasive imaging techniques have been introduced to improve the early detection of melanoma which can be often challenging with the naked eye alone [10, 17]. Dermoscopy and RCM were shown to improve diagnostic accuracy in this field [11, 32, 33, 48]. In a recent study, almost no melanomas were misclassified and consequently undertreated when both techniques were used in combination: sensitivity 97.82 %, specificity 92.44 %, PPV 87.37 % and NPV 98.75 % [1]. Moreover, a significant reduction in numbers needed to excise could be determined by combining both devices [1, 43]. The problem with dermoscopy and RCM, however, is that their diagnostic performances varies with the user’s experience and would be consequently poor for non-expert physicians [2, 31]. A user-independent, automated classification of pigmented skin lesions based on their optical properties could be achieved by other techniques, such as spectrophotometric technologies [23]. The diagnostic potential of HD-OCT is traditionally thought to be insufficient for ruling out the diagnosis of melanoma, due to limitations of its cellular resolution, as compared to RCM [20]. However, it was shown that HD-OCT allows quantifying the skin reflectance as function of depth [7]. The present study is—to our knowledge—the first one describing the optical properties of benign and malignant melanocytic skin lesions assessed in vivo by means of HD-OCT. We assessed the scattering properties of 45 melanocytic lesions, including 15 benign nevi, 15 dysplastic nevi and 15 melanomas.

Only ballistic photons backscattered to the detection system of the HD-OCT contribute to the image; the decay with depth of focus of reflectance of ballistic photons is a process that can be well modelled by an exponential relationship [26, 29]. The technique of semi-log plot whereby an exponential function becomes a straight line given by y = ax + b has been implemented on HD-OCT signals coming from four successive skin layers (epidermis, upper papillary dermis, deeper papillary dermis and superficial reticular dermis). This permitted the HD-OCT in vivo measurement of skin entrance signal (SES), relative attenuation factor normalized for the skin entrance signal (µ raf1) and half value layer (z 1/2).

A highly significant difference between each melanocytic group could be observed with regard to the measured optical properties. A reduction in µ raf1 and SES and increase of z 1/2 could be noticed with increasing malignancy. Spectrophotometric studies provided evidence of decreased reflectance with increased malignancy of melanocytic lesions [37]. From a clinical point of view, it is more disastrous if MM is misdiagnosed as BN then if a DN is mistaken for a MM. Hence, decision criteria minimizing the false-negative rate deserve priority. Therefore, critical values could be determined for these optical properties permitting differentiation of malignant melanoma from non-malignant melanocytic lesions with high area under the roc curve (AUC) values (0.98, 0.86 and 0.96, respectively) and high positive likelihood ratio’s (28, 4.8 and 7, respectively).

The diagnostic performance of HD-OCT in discriminating MM from non-melanoma based on 3-D cellular and micro-architectural morphological features has been found to be moderate with an NPV (Negative Predictive Value) of 89.7 % [20]. High false-negative rates were observed in very thin MM and high false-positive rates in DN [20]. The present study dealt with superficial spreading melanomas and provided evidence that the diagnostic accuracy of HD-OCT based on the optical properties, µ raf1 , SES and z 1/2 scored much better (95.6, 82.2 and 88.9 %, respectively). High NPV could be found for these optical properties (96.7, 89.3 and 96.3 %, respectively) reducing the risk of mistreating a malignant lesion to a more acceptable level (3.3–7.3 % instead of 11.3 %). The role of spectrophotometry (SPT) in early diagnosis of melanoma in 10 studies has been compared [2]. The authors concluded that low diagnostic accuracy with moderate NPV represented the main hamper for the introduction of SPT technology in clinical practice. HD-OCT enables the rather unique combination of in vivo morphological analysis of cellular and micro-architectural structures with in vivo analysis of optical properties of tissue scatterers.

Optical property measurements sensitive for scattering and anisotropy factor g offer the best means of characterizing the micro-architecture of cells and tissues in general and melanocytic lesions more specifically [27, 45, 47, 49, 55]. Although a lot of variation in data exists, in general there appears to be a trend toward increasing g and decreasing reduced scattering coefficient µ′s as the wavelength increases. In normal skin with photo type I-III, at 1300 nm (second NIR optical window) the mean value for reduced scattering coefficient µ′s is 0.001 (±0.0005) [53]. This value is proportional to the measured µ raf measured by HD-OCT. In the upper layer of normal skin at inner site of upper arm significant differences in µ raf1 according to intrinsic ageing could be measured by HD-OCT (unpublished data). In pre-menopausal and post-menopausal females the relative attenuation factor was found to be 0.0114/µm (±0.0012/µm) and 0.0375/µm (±0.0037/µm), respectively. It also appeared that µ raf4 at superficial reticular dermis was significantly higher in pre-menopausal women compared to post-menopausal women; 0.016/µm (±0.0012/µm) and 0.0094/µm (0.0005/µm), respectively (unpublished data). This suggested an intrinsic ageing-related decrease in anisotropy of scattering in upper layers and increase in anisotropy of scattering in lower layers.

In the epidermis, a significant difference (p < 0.001) in relative attenuation factor normalized to skin entrance signal could be observed among the three melanocytic subgroups. This factor is proportional to the µ′s when a light source of 1300 nm is used [26, 29]. A strong reduction of g and consequently a robust increase in µ′s could be noticed in BN compared to normal values [29]. Hence, a reduction of factor µ raf1 with increasing malignancy which could be observed in our study implied an important difference in anisotropy factor g among the three melanocytic subgroups. Lower g values (less than 0.8) corresponding with more isotropic scattering seemed to be characteristic for BN. In BN the µ′s of the upper layer was approximately 3 times higher compared to normal skin. High g values (higher than 0.9) corresponding with strong anisotropic scattering seemed to be characteristic for melanoma. Forward scattering became more likely. In melanoma the upper layer µ′s seemed to be 10 times lower than normal skin. Hence, the epidermis became much more translucent compared to BN and DN. As a consequence NIR light penetrated much deeper in melanoma than in DN or BN. This probably explains the observed differences in dermoscopic imaging.

Cellular and micro-architectural risk parameters used for diagnosing MM have been determined for HD-OCT such as large roundish pagetoid cells, atypical cell clusters at dermo-epidermal junction, totally disarranged epidermal/dermal pattern and large vertical icicle-shaped structures [9, 20]. With regard to the impact of scattering, alterations of distribution, size, density and orientation of scatterers with malignancy probably are of higher importance than the content of pigment itself because scattering is dependent on the refractive index mismatch between cellular components, extracellular matrix fibres and the extracellular fluid [51]. Dysplastic nevi are characterized by nuclear enlargement, slight irregularity and hyperchromasia, with clumping of chromatin and sometimes with prominent nucleoli. In melanoma besides nuclear atypia, nests and single melanocytes of variable sizes, shapes are present in the epidermis in a pagetoid pattern [16]. A decline of the relative attenuation factor demonstrated in our study suggesting an increase of g value is in line with former observations stating that the average scatterer size in melanoma is significantly larger than those of DN and BN [21, 39, 40].

Some important limitations need to be addressed. This study is a retrospective pilot study. However, this type of study is essential to gather information to design an appropriate (e.g. adequate sample size) prospective study. A second issue is the selected region of interest for optical analysis. This selection needs some experience with morphology analysis of HD-OCT images of melanocytic lesions. This fact questions to some extent the accuracy of HD-OCT optical analysis in the hands of non-experts. Last but not least, the effect of age and gender on optical properties was not taken into account.

In conclusion, HD-OCT seems to enable the rather unique combination of in vivo morphological analysis of cellular and 3-D micro-architectural structures with in vivo analysis of optical properties of tissue scatterers in melanocytic lesions. In vivo HD-OCT analysis of optical properties permits melanoma diagnosis with higher accuracy than in vivo HD-OCT analysis of morphology alone. The diagnostic performance of HD-OCT in MM should be further assessed in other clinical settings combining both types of analysis.