Abstract
Inadequate margins in breast-conserving surgery (BCS) are associated with an increased likelihood of local recurrence of breast cancer. Currently, approximately 20% of BCS patients require repeat surgery due to inadequate margins at the initial operation. Implementation of an accurate, intraoperative margin assessment tool may reduce this re-excision rate. This study determined, for the first time, the diagnostic accuracy of quantitative micro-elastography (QME), an optical coherence tomography (OCT)–based elastography technique that produces images of tissue microscale elasticity, for detecting tumor within 1 mm of the margins of BCS specimens. Simultaneous OCT and QME were performed on the margins of intact, freshly excised specimens from 83 patients undergoing BCS and on dissected specimens from 7 patients undergoing mastectomy. The resulting three-dimensional images (45 × 45 × 1 mm) were coregistered with postoperative histology to determine tissue types present in each scan. Data from 12 BCS patients and the 7 mastectomy patients served to build a set of images for reader training. One hundred and fifty-four subimages (10 × 10 × 1 mm) from the remaining 71 BCS patients were included in a blinded reader study, which resulted in 69.0% sensitivity and 79.0% specificity using OCT images, versus 92.9% sensitivity and 96.4% specificity using elasticity images. The quantitative nature of QME also facilitated development of an automated reader, which resulted in 100.0% sensitivity and 97.7% specificity. These results demonstrate high accuracy of QME for detecting tumor within 1 mm of the margin and the potential for this technique to improve outcomes in BCS.
An optical imaging technology probes breast tissue elasticity to provide accurate assessment of tumor margin involvement in breast-conserving surgery.
Introduction
Success of breast-conserving surgery (BCS) is characterized by clear margins and a good cosmetic outcome for the patient. However, approximately 20% of cases require repeat surgery due to inadequate margins (1–3). Re-excision surgery causes substantial physical, psychologic, and financial burdens for patients, with higher risk of complications (4), worse cosmesis, and additional costs, on average, of >$10,000 per patient (5). Intraoperative detection of tumor at the margins would allow more complete resection of malignant tissue in the first operation, provide the surgeon with confidence that no residual cancer remains in the breast, and reduce the number of re-excision surgeries.
Several techniques are currently used for intraoperative margin assessment. Intraoperative pathologic assessment can be performed using frozen section analysis and imprint cytology (6), but these techniques are resource-intensive, sample only a small percentage of the surgical margins, and have limited efficacy, especially for ductal carcinoma in situ (DCIS; ref. 7), and, thus, have not been widely adopted (8). Another available technique, intraoperative specimen radiography (IOSR), provides an X-ray image of the excised specimen. IOSR can show that the main lesion has been removed, which is particularly useful for nonpalpable tumors that are localized using preoperative wire insertion (9) or radioactive seed placement (10). However, IOSR accurately predicts margin status in only 48% of cases (11) and has been shown not to reduce re-excision rates (11, 12). Intraoperative ultrasound guidance of excision has been shown in a small number of studies to reduce re-excision rates by more than half for invasive cancers (13, 14), but ultrasound is operator-dependent and has limited reliability for visualizing in situ or multifocal cancers (15, 16).
To address the shortcomings of existing techniques, a range of margin assessment tools have been proposed, relying on various contrast mechanisms to detect cancer. One commercially available technique uses radiofrequency spectroscopy to measure the dielectric properties of tissue, implemented in a handheld probe (MarginProbe, Dune Medical Devices). The technique shows high accuracy for margin assessment in homogeneous, ex vivo tissues (e.g., if the probe is placed over a large region of tumor, >6 mm diameter) but has lower accuracy (70% sensitivity and 70% specificity) when multiple tissue types are present in the interrogated region (17). In a randomized trial with 596 patients, use of MarginProbe on excised specimens reduced re-excision rate from 25.8% to 19.8% but suffered from low specificity with a 53.6% false positive rate (18).
Fluorescent modalities that utilize molecular contrast to intraoperatively highlight cancer are also under development, potentially enabling surgeons to visualize tumor in the cavity, as well as in the excised lump (19–22). A small number of proof-of-principle studies in humans have been reported (19, 20, 22), but the efficacy of these techniques for reducing re-excision rates in BCS has not been determined. In addition, they typically rely on preoperative, systemic administration of exogenous dyes, requiring extensive dosing and tumor uptake studies, and creating potential barriers to clinical translation (19). Another optical technique, surface-enhanced Raman scattering (23), has recently been proposed, using targeted nanoparticles for multiplexed imaging of cancer biomarkers. A preliminary study on ex vivo tissues dissected from mastectomy and lumpectomy specimens showed promising sensitivity (89.3%) and specificity (92.1%) for cancer detection, and the technique is conducive to personalized biomarker imaging based on tumor-specific molecular profiles. However, it is restricted to surface imaging, and, while many institutions have adopted a “no tumor on ink” criterion for invasive cancers, larger margins are typically desired for in situ cancers (1). Label-free optical techniques have also been proposed, including diffuse reflectance spectroscopy/hyperspectral imaging (24, 25), autofluorescence lifetime imaging, and Raman spectroscopy (26). However, diffuse techniques can suffer from low spatial resolution (∼5 mm; ref. 24), while autofluorescence and Raman techniques typically have low scanning speeds (12–24 minutes/margin; ref. 26), making clinical translation more challenging.
Optical coherence tomography (OCT) is a promising optical technique capable of three-dimensional (3D), high-speed, high-resolution imaging without need for exogenous contrast agents. OCT may be described as an optical analogue to ultrasound. It uses interferometry to effectively measure “time of flight” of light in tissue, creating an image based on the amount of back-scattered light, with microscopic resolution (2–10 μm) up to 1 to 2 mm in depth. These imaging specifications match well with the clinical requirements of margin assessment in BCS. Preliminary studies have reported high sensitivity (80–94%) and specificity (87–93%) for detecting cancer, primarily in mastectomy tissues (27–29). However, these mastectomy samples were typically dissected to expose a bulk of dense, high-grade tumor for imaging. The image contrast in this scenario may not translate to margin assessment in BCS, in which typically lower-grade tumors must be detected at the edge of intact specimens or, for detecting residual tumor, directly within the surgical cavity. The largest study to date that used OCT to assess margins in BCS specimens reported sensitivity of 55% to 65% and specificity of 68% to 70% (30). One reason for this relatively low accuracy may be the limited ability of OCT to distinguish between tumor and surrounding normal stroma (31–33). Stronger contrast between tissue types is expected to aid surgeon decision-making in the intraoperative setting (34, 35).
Beyond its distinct molecular and optical properties, breast cancer also exhibits distinct mechanical properties (36–38). At the cellular scale, atomic force microscopy has revealed unique mechanical signatures within the breast tumor microenvironment caused by the mix of cellular proliferation and desmoplastic stroma (37). At the macro-scale, surgeons rely on mechanical changes as they manually palpate the tissue, feeling for the boundaries of the typically stiff lesion. However, palpation is a subjective tool, and a large proportion of breast lesions are considered “impalpable,” that is, too small or soft to detect through touch (39). Elastography is a technique that creates images of the mechanical properties of tissue, complementing palpation by visualizing mechanical changes in 2D or 3D. Elastography based on ultrasound has been developed for a number of applications, including preoperative diagnosis of breast lesions (38, 40), but has not been applied to intraoperative margin assessment, mainly due to its relatively low spatial resolution (see Supplementary Information, “Note on elastography of the breast.”) In an emerging technique, OCT elastography is used to measure tissue deformation under an applied load, offering 3D maps of mechanical properties with microscale resolution (41). OCT elastography techniques can be classified according to the mechanical loading mechanism used, with compression and shear wave being the most prominent (41). Although early compression OCT elastography studies produced maps of tissue deformation (strain; refs. 31, 42), which is an indirect and qualitative measure of elasticity, introduction of a stress sensing technique to map the local stress in 2D at the tissue surface, has enabled quantitative micro-elastography (QME), providing 3D maps of local elasticity, under the assumption of uniaxial stress (43). Preliminary QME data in mastectomy specimens showed that elasticity images provide additional contrast between tumor and normal tissue compared to OCT and strain, and the technique has been extended to incorporate a wide-field scanning mechanism that enables entire margins of BCS specimens to be imaged within an intraoperative timeframe (44). To build on these promising feasibility studies and to determine the clinical potential of QME for intraoperative margin assessment, it must be established whether clinicians can interpret QME images to accurately identify close or positive tumor margins in BCS specimens. To this end, the goal of this study is to conduct, for the first time, a blinded reader study, with postoperative histology as the gold standard, to determine the diagnostic accuracy (sensitivity and specificity) of QME, compared with OCT alone, for detecting tumor within 1 mm of the margins of freshly excised specimens from patients undergoing BCS.
Patients and Methods
Patient recruitment and imaging
All study procedures were performed after approval by the ethics board of the South Metropolitan Area Health Service in Western Australia. Ninety patients were recruited for this study after informed written consent was obtained: 83 patients (with no prior excision) undergoing BCS for treatment of breast cancer and 7 patients undergoing mastectomy for treatment of breast cancer. Samples from the 7 mastectomy patients and 12 of the BCS patients were used to create a set of pilot data for training readers. Data from the remaining 71 BCS patients were included in the blinded reader study. Table 1 summarizes the disease characteristics of all BCS patients. Supplementary Table S1 summarizes mastectomy patient characteristics.
. | Pilot study (12 patients) . | Reader study (71 patients) . | ||
---|---|---|---|---|
Age, y . | Number . | Percentage . | Number . | Percentage . |
Mean | 60 | 59 | ||
SD | 9 | 11 | ||
Range | 44–74 | 26–76 | ||
≤65 | 8 | 67% | 48 | 68% |
>65 | 4 | 33% | 23 | 32% |
Surgical diagnosisa | ||||
Ductal carcinoma in situ | 6 | 50% | 44 | 62% |
Invasive ductal carcinoma | 8 | 67% | 43 | 61% |
Invasive lobular carcinoma | 2 | 17% | 8 | 11% |
Invasive mucinous carcinoma | 1 | 8% | 3 | 4% |
Mixed invasive ductal lobular carcinoma | 0 | 0% | 2 | 3% |
Invasive micropapillary carcinoma | 0 | 0% | 1 | <1% |
Invasive apocrine carcinoma | 0 | 0% | 1 | <1% |
Invasive solid papillary carcinoma | 0 | 0% | 2 | 3% |
No tumorb | 0 | 0% | 3 | 4% |
Palpability | ||||
Palpable | 7 | 58% | 45 | 63% |
Impalpable | 5 | 42% | 26 | 37% |
Total lesion size (greatest dimension) | ||||
<1 cm | 1 | 8% | 15 | 21% |
1–2 cm | 5 | 42% | 27 | 38% |
>2 cm | 6 | 50% | 26 | 37% |
No tumor | 0 | 0% | 3 | 4% |
Closest tumor marginc | ||||
<1 mmd | 7 | 58% | 32 | 45% |
1–3 mm | 4 | 33% | 25 | 35% |
>3 mm | 1 | 8% | 11 | 15% |
Not applicable (no tumor) | 0 | 0% | 3 | 4% |
. | Pilot study (12 patients) . | Reader study (71 patients) . | ||
---|---|---|---|---|
Age, y . | Number . | Percentage . | Number . | Percentage . |
Mean | 60 | 59 | ||
SD | 9 | 11 | ||
Range | 44–74 | 26–76 | ||
≤65 | 8 | 67% | 48 | 68% |
>65 | 4 | 33% | 23 | 32% |
Surgical diagnosisa | ||||
Ductal carcinoma in situ | 6 | 50% | 44 | 62% |
Invasive ductal carcinoma | 8 | 67% | 43 | 61% |
Invasive lobular carcinoma | 2 | 17% | 8 | 11% |
Invasive mucinous carcinoma | 1 | 8% | 3 | 4% |
Mixed invasive ductal lobular carcinoma | 0 | 0% | 2 | 3% |
Invasive micropapillary carcinoma | 0 | 0% | 1 | <1% |
Invasive apocrine carcinoma | 0 | 0% | 1 | <1% |
Invasive solid papillary carcinoma | 0 | 0% | 2 | 3% |
No tumorb | 0 | 0% | 3 | 4% |
Palpability | ||||
Palpable | 7 | 58% | 45 | 63% |
Impalpable | 5 | 42% | 26 | 37% |
Total lesion size (greatest dimension) | ||||
<1 cm | 1 | 8% | 15 | 21% |
1–2 cm | 5 | 42% | 27 | 38% |
>2 cm | 6 | 50% | 26 | 37% |
No tumor | 0 | 0% | 3 | 4% |
Closest tumor marginc | ||||
<1 mmd | 7 | 58% | 32 | 45% |
1–3 mm | 4 | 33% | 25 | 35% |
>3 mm | 1 | 8% | 11 | 15% |
Not applicable (no tumor) | 0 | 0% | 3 | 4% |
aMultiple tumor types may occur in a given patient, such that diagnosis percentages add to >100%.
bTumor was excised in original core biopsy.
cIncludes deep and superficial margins in addition to radial margins.
dThirteen (18%) of BCS specimens in the reader study had positive margins (i.e., tumor on ink).
For patients undergoing BCS, following surgery, the fresh, intact specimens were transferred to the pathology department at Fiona Stanley Hospital. Pathologists at this institution dictated that tissue should be placed in fixative (formalin) within 1 hour of receiving the specimen from surgery, to avoid any tissue degradation that might influence histological processing. Between excision and imaging, specimens were kept in air at room temperature, and the surface was kept hydrated by applying droplets of saline. Specimen orientation was maintained using clips and sutures, per standards at this institution. In most cases (8/12 pilot cases and 68/71 cases for the reader study), two margins were imaged. In the remaining cases, one margin was imaged due to delays in scanning the specimen within the allotted time. No more than two margins per specimen were scanned as, in addition to imaging, sufficient time was needed to transfer the specimen from pathology, select a margin for scanning and orientate the specimen in the imaging system. Margins for scanning were chosen based on consultation with pathologists, observation of the intraoperative specimen radiography (when available), and surgical notes indicating if extra cavity shavings were performed intraoperatively. The “closest tumor margin” distances reported in Table 1 were determined by postoperative histology and include superficial and deep margins (considering radial margins only, which are most clinically relevant, the rate of tumor margins <1 mm in the reader study was 21%). For patients undergoing mastectomy, nondiagnostic tissue was dissected by a pathologist to create samples of approximately 5 × 5 × 0.5 cm, which were also imaged within approximately 1 hour of excision.
Imaging was performed on a benchtop, wide-field QME system as described previously (Supplementary Fig. S1A; ref. 44). Briefly, the system is based on a Telesto II spectral-domain OCT system (TEL220C1, Thorlabs). It uses a superluminescent diode light source with a central wavelength of 1300 nm and a bandwidth of >170 nm, illuminating the sample with 2.5 mW of power (a low power level that is safe for users). The measured axial and lateral resolutions in air are 5.5 μm and 13 μm, respectively. The system captures one-dimensional axial scans (A-scans) in 14 μs, and the beam is raster scanned to build a 3D image measuring 16 × 16 × 3.5 mm in 55 seconds. Wide-field images are generated by translating the specimen relative to the OCT scan head between 3D acquisitions, as described previously (42, 44). Nine subvolumes are acquired with 1-mm overlap in the lateral plane, resulting in a 45 × 45 × 3.5 mm image captured in under 9 minutes. The partially overlapping subvolumes are stitched to form mosaicked wide-field images, presented in the en face plane. The measured displacement sensitivity of the OCT system is 1.4 nm at an OCT signal-to-noise ratio of 40 dB, acquired under clinical testing conditions in the pathology laboratory (i.e., without a vibration isolation table).
A compliant silicone layer is placed between the specimen and imaging window to ensure contact between the uneven specimen surface and window, as well as to estimate stress for quantification of elasticity (43). After establishing contact, microscale displacement (up to 9.5 μm) was applied to the specimen surface using a piezoelectric actuator. Images were processed to determine stress, strain, and elasticity (Supplementary Fig. S1B), and elasticity was overlaid on the solid regions of the OCT images for visualization, using an algorithm described previously (42). OCT data are displayed in grayscale from 0 to 40 dB, and elasticity data are displayed in color on a logarithmic scale from 3.63 to 363 kPa.
Histopathology and coregistration with imaging data
Following imaging, specimens were fixed in formalin and submitted for standard histopathologic processing. Specimens were inked for orientation and sliced in “bread-loaf” fashion from lateral-to-medial, superficial-to-deep, or superior-to-inferior margins to generate histologic sections at each plane. This sectioning protocol resulted in histology typically being performed in a plane orthogonal to the en face plane in which QME images are displayed. The positions of the histology sections were coregistered with images by cross-referencing the blocking diagram and photographs taken of the specimen during scanning, as illustrated in Supplementary Fig. S2. Pathologists (BFD, BL) determined tissue types present in the histology images. Engineers and pathologists worked together to confirm co-registration by cross-checking tissue features present in the histology images with similar features present in the OCT images, using the depth cross-section (B-scan) view of OCT to match the histology plane where needed, as described in (42).
To focus the histologic analysis, and to facilitate the subsequent reader study, 3D regions of interest (ROI), 10 × 10 × 1 mm, were selected from the wide-field (45 × 45 × 1 mm) scans. At least one ROI was selected on every margin scanned in this study. ROIs were then included in the reader study if they met the following criteria: good physical contact with the specimen; tissue not extensively damaged by thermal effects (due to cauterization during resection) as assessed by histopathology; and availability of a reliable histology match for the ROI. Considering that histology slices are typically taken approximately every 3 mm, 1 to 3 sections were available for a given ROI. ROIs were designated as “positive” for cancer if the pathologist identified any tumor within 1 mm of the margin in the histology sections corresponding to the ROI, similar to designation of “positive” margins in other diagnostic accuracy studies in BCS specimens (26, 30). In each case, the ROI was selected before the histology was analyzed by the pathologist. As such, the researchers selecting the ROI did not know the outcome of the corresponding histology. Note: The clinical standard at the hospital where the study was performed is 1 to 2 mm for invasive ductal carcinoma and 2 mm for DCIS.
Reader study and statistical analysis
The reader study was designed to have three main outcomes: diagnostic accuracy (sensitivity and specificity) of OCT, QME, and OCT plus QME. Seven readers [2 surgeons, 2 engineers, 1 medical sonographer, 1 pathology scientist (equivalent to a pathology assistant) and 1 medical resident], blinded to the histologic results, participated in the study. First, readers were trained to read OCT images in a 1.5-hour session. They then completed a reading of all ROIs (OCT alone, using criteria in Fig. 1A) within 10 days of training and had access to the training images during reading. Two weeks following the OCT training, readers were trained to read QME images in another 1.5-hour session. They then completed a reading of all ROIs once more, this time viewing OCT and QME side-by-side, from which accuracy was calculated for QME (using criterion in Fig. 1B) and OCT plus QME (using criterion in Fig. 1C). To remove potential bias due to reader memory, the order and orientation of ROIs were randomized between readings. In both readings, readers viewed the ROIs (10 × 10 × 1 mm, presented as a stack of en face images displayed every 20 μm up to 1 mm in depth) using the open source software ImageJ (v1.52a; ref. 45), which allowed scrolling through depth to visualize the volume. Readers completed their evaluation using a custom-built interface in which they followed a set of criteria for the presence of cancer in each image type and selected “cancer” or “not cancer” as the endpoint (Fig. 1).
Criteria for cancer in OCT images (Fig. 1A) were determined on the basis of three main observations from prior OCT studies of breast tissues (28, 32, 42, 46, 47) and pilot study data. First, it has been established that adipose tissue is easily distinguished from other tissues in breast based on its “honeycomb” structure in OCT. Therefore, readers were trained to first identify solid features >1 mm in diameter. If the ROI was made up of all adipose tissue or adipose tissue with only small (<1 mm), isolated solid features, it was considered “not cancer”. Second, it was found in preliminary OCT studies on BCS specimens (42) and in this study that cancer near the margin tends to have continuity with depth, uninterrupted by regions of adipose tissue. Thus, if readers found a region of solid tissue, they were to scroll through the 3D image, checking if the feature extended down into the specimen. Finally, a criterion was added to help distinguish benign stroma from regions of cancer. In benign stroma, a striated or banded pattern is observed in OCT, likely due to the organization of the underlying collagen. Cancer tends to disrupt this organization and result in heterogeneous patterns in OCT (31, 48). One special case is mucinous carcinoma (present in four recruited patients), in which tumor cells produce pools of a liquid substance known as mucin. Mucinous pools result in regions of low signal in OCT images, appearing similar to adipose tissue, but can often be identified by the outline of the pool, which has higher OCT signal. Thus, readers evaluated OCT images for the presence of these structures, even if at first it appeared that the tissue was adipose tissue (bottom path, Fig 1A).
In QME images, readers were trained to look for areas of high elasticity, based on studies across spatial scales showing that cancer and its associated stroma are stiffer than benign tissues (36, 37, 44). To avoid mistaking small, isolated patches of stiff stroma for cancer, the region of high elasticity had to cover at least 75% of a 1-mm diameter circle (Fig. 1B). This diameter was chosen empirically based on analysis of how tumor presented in QME images both in previous studies (43, 44) and in the training data. The user interface enabled readers to calculate this precisely by dragging their cursor to any location in the ROI. High elasticity was defined as >26.3 kPa, based on analysis of prior and current study data and is denoted by a black line on the color scale of all presented QME images. To assess OCT + QME in combination, readers were again directed to assess continuity with depth and ensure that the region was solid tissue, in addition to the elasticity criterion (Fig. 1C). Mucinous carcinoma presented a special scenario (second step in Fig. 1C): if an area had high elasticity but appeared in the OCT to lack “solid tissue,” this was likely a mucinous carcinoma, and, as such, results in “cancer” in the decision tree. This was based on our observation that even these liquid-dominant tumors resulted in high elasticity, likely due to interstitial fluid pressure within the tissue.
Finally, the QME criterion, facilitated by its quantitative nature, was implemented into a preliminary automated algorithm. This involved thresholding the elasticity values in the QME images at every depth to generate binary images (with ones representing pixels with an elasticity >26.3 kPa), then convolving these binary images with a 1 mm diameter circular kernel. The values at each pixel in the kernel were normalized such that the result of the convolution equaled the percentage of the circle covered by high elasticity. If, for a given ROI, the convolution produced a value ≥75% at any location in the ROI, that ROI was interpreted as containing cancer. As such, this algorithm was equivalent to the QME reader criteria shown in Fig. 1B. This algorithm was implemented on a standard desktop computer using MATLAB (Mathworks, R2016a), and took approximately one minute to read all ROIs.
Following the readings, the sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and accuracy were calculated for each reader. 95% confidence intervals were estimated using the score interval for a binomial case (“Wald” interval; ref. 49). Aggregates were computed by taking the sum of the individual results of the seven readers (the totals of true positive, false positive, true negative, and false negative counts), and computing the sensitivity, specificity, PPV, NPV, accuracy, and confidence intervals on these summed totals. The performance of QME and OCT plus QME were compared against OCT alone using McNemar test (50) for statistical significance. Interreader agreement was quantified using the Fleiss' Kappa metric (51). Detailed statistical calculations are provided in Supplementary Reader Statistics.
Results
OCT and QME images of malignant and benign breast tissues at the margin
Intact BCS specimens from 71 patients were scanned for the reader study, without any damage to the tissue nor disruption to typical histopathologic protocol. One hundred and seventy-four ROIs from the 71 patients were selected and coregistered with histology. Of these, 154 ROIs were included in the reader study. Table 2 summarizes the tissue types and characteristics of the included ROIs. Reasons for exclusion were extensive thermal damage, as determined by postoperative histology (n = 3); inconclusive registration with histology sections (n = 4); insufficient contact between the specimen, silicone layer, and the imaging window (n = 3); imaging artifacts (e.g., stray reflections or surgical clips precluded image interpretation, n = 4); insufficient elasticity data overlaid on solid tissue (n = 2); a rare form of mucinous DCIS (insufficient data to build decision criteria, n = 2); and the region of tumor within the ROI having a size <1 mm (n = 2). Of the 154 ROIs included in this study, 24 had cancer within 1 mm of the surface, corresponding to a prevalence of 15.6%. Most were made up of a mix of tissue types, including adipose tissue, stroma, and parenchymal tissues. See Supplementary Table S2 and "Note on prevalence" in Supplementary Information for more details on tumor prevalence.
ROI selection . | Number . | . |
---|---|---|
Imaged margins | 139 | |
ROIs coregistered with histology | 174 | |
ROIs included in reader study | 154 | |
Tissue types in included ROIs | Percentage (/154) | |
Cancer within 1 mm (positive ROI) | 24 | 15.6% |
Invasive (solid) | 14 | 9.1% |
Invasive (mucinous) | 3 | 1.9% |
DCIS | 7 | 4.5% |
Clear within 1 mm (negative ROI) | 130 | 84.4% |
Predominantly adipose tissue | 38 | 24.7% |
Adipose tissue with dense stroma | 37 | 24.0% |
Adipose tissue with ducts/vessels | 24 | 15.6% |
Adipose tissue with strands of connective tissue | 19 | 12.3% |
Dense stroma | 12 | 7.8% |
ROI selection . | Number . | . |
---|---|---|
Imaged margins | 139 | |
ROIs coregistered with histology | 174 | |
ROIs included in reader study | 154 | |
Tissue types in included ROIs | Percentage (/154) | |
Cancer within 1 mm (positive ROI) | 24 | 15.6% |
Invasive (solid) | 14 | 9.1% |
Invasive (mucinous) | 3 | 1.9% |
DCIS | 7 | 4.5% |
Clear within 1 mm (negative ROI) | 130 | 84.4% |
Predominantly adipose tissue | 38 | 24.7% |
Adipose tissue with dense stroma | 37 | 24.0% |
Adipose tissue with ducts/vessels | 24 | 15.6% |
Adipose tissue with strands of connective tissue | 19 | 12.3% |
Dense stroma | 12 | 7.8% |
Figures 2 to 4 present ROIs from the reader study, representative of the various tissue types encountered in the study. In each case, en face, 10 × 10 mm OCT and QME images are displayed at a depth of 40 to 60 μm (entire 3D volumes available in Supplementary Videos S1–S6), along with the corresponding hematoxylin and eosin histology. The plane and orientation of the histology is indicated by the red dashed line in each set of images.
Figure 2 shows ROIs containing benign breast tissues: adipose tissue (Fig. 2A–C) and a mix of adipose tissue and stroma containing benign ducts (Fig. 2D–F). The OCT in Fig. 2A (Supplementary Video S1) depicts the honeycomb structure typical of adipose tissue. Thin strands of connective tissue are present throughout this ROI, and QME shows that the elasticity in these regions tends to be low (Fig. 2B). Six of seven readers called this benign using OCT alone, compared to all readers using QME combined with OCT. A larger region of benign stroma is shown at the margin in Fig. 2D, and 3D analysis reveals that this structure lacks continuity, with intermittent adipose tissue present with depth (Supplementary Video S2), similar to the pattern in the histology image (Fig. 2F). However, using OCT alone, three readers incorrectly assessed this as cancer. In QME (Fig. 2E), this stroma exhibits mostly low elasticity. Very small pockets of high elasticity are seen at a few points throughout the stroma, but only features with more widespread elasticity were considered suspicious for cancer. All seven readers correctly called this benign using QME and QME plus OCT.
Figure 3 shows positive ROIs from two BCS cases, each containing invasive ductal carcinoma on the margin. In Fig. 3C, the histology image reveals invasive ductal carcinoma grade 3 with high cellular density, surrounded by adipose tissue. Correspondingly, the OCT (Fig. 3A) captured a region of solid tissue with a nodular appearance (no indication of striated, organized collagen) that matches the location of tumor in the histology image. In QME (Fig. 3B), this region shows elevated elasticity, due to the much higher cellular density in the tumor compared with the surrounding adipose tissue. This detected region of tumor is approximately 2 mm across, making up <0.2% of the total surface area of this particular margin, demonstrating the high resolution of the technique. All seven readers evaluated this ROI correctly using OCT, QME, and the combination. The histology image in Fig. 3F shows invasive ductal carcinoma grade III surrounded by a mix of stroma and adipose tissue. A mix of adipose tissue and solid tissue is also present in the OCT image (Fig. 3D), and QME reveals a portion of the solid tissue to have high elasticity (Fig. 3E). 3D analysis of the images (Supplementary Video 4) shows that the region with high elasticity is also continuous with depth (i.e., adipose tissue is not revealed underneath). On the basis of OCT alone, one of the seven readers called this as a false negative, mistaking it for benign, whereas all seven readers evaluated this as cancer using QME alone or alongside OCT.
Figure 4 shows ROIs containing two malignant tissue types: DCIS (Fig. 4A–C) and mucinous carcinoma (Fig. 4D–F), which were designated “impalpable” preoperatively and required hookwire guidance for lesion excision. Despite being impalpable on a macroscale, changes in mechanical properties on a microscale are detectable by QME for these tissue types (43). Figure 4A–C shows an example of DCIS within 0.15 mm of the margin; the histology image (Fig. 4C) shows the involved duct surrounded by dense, benign stroma. The OCT (Fig. 4A) shows no apparent contrast between this duct and the surrounding stroma, but the QME (Fig. 4B) highlights high elasticity corresponding to the region of DCIS. This elevated elasticity corresponds to the dense, proliferating cells within the duct, as well as a fibrotic stromal response (52) immediately surrounding the duct. Throughout the stroma in the rest of the ROI, there are small, localized areas of elevated elasticity, but most do not meet the size criterion for cancer defined in this study. Six of seven readers evaluated this correctly using OCT alone, as the structure shows continuity with depth (Supplementary Video S5), and all readers evaluated it correctly using QME. In Fig. 4D–F (Supplementary Video S6), a region of benign stroma and adipose tissue neighbors a region of mucinous carcinoma on the margin. In OCT (Fig. 4D), the region corresponding to the mucinous tumor is difficult to distinguish from the surrounding adipose tissue, and one reader designated this as benign. However, in QME (Fig. 4E) the tumor exhibits elevated elasticity, and all readers correctly classified this ROI as cancerous when QME was available.
Blinded reader study results in higher accuracy using QME over OCT alone
Table 3 and Supplementary Fig. S3 summarize the results of the reader study. Readers completed evaluation of the 154 OCT plus QME ROIs in an average of 4 hours, or approximately 90 seconds per ROI. Interreader agreement was nearly perfect for QME and OCT plus QME, and moderate for OCT alone, according to a standard interpretation of the Fleiss' Kappa index for interrater agreement. Reported ranges for sensitivity and specificity in Table 3 indicate 95% confidence intervals. On the basis of the aggregate results, the sensitivity, specificity, and accuracy of OCT for detecting cancer within 1 mm of the margin were 69.0%, 79.0%, and 77.5%, respectively. Using QME, sensitivity, specificity, and accuracy were 92.9%, 96.4%, and 95.8%, respectively. Using combined OCT and QME criteria, sensitivity, specificity, and accuracy were 80.4%, 99.5%, and 96.5%, respectively. Sensitivity was significantly (P < 0.05) improved using QME over OCT alone for 6 of the 7 readers, and specificity was significantly improved using QME over OCT alone for all 7 readers (see Supplementary Reader Statistics for calculations).
. | OCT . | QME . | OCT + QME . | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Reader . | Sensitivity . | Specificity . | PPV . | NPV . | Acc . | Sensitivity . | Specificity . | PPV . | NPV . | Acc . | Sensitivity . | Specificity . | PPV . | NPV . | Acc . |
Eng 1 | 62.5 ± 19.4% | 85.4 ± 6.1% | 44% | 93% | 82% | 100.0 ± 0.0% | 97.7 ± 2.6% | 89% | 100% | 98% | 87.5 ± 13.2% | 100 ± 0.0% | 100% | 98% | 98% |
Eng 2 | 70.8 ± 18.2% | 86.9 ± 5.8% | 50% | 94% | 84% | 100.0 ± 0.0% | 97.7 ± 2.6% | 89% | 100% | 98% | 87.5 ± 13.2% | 100 ± 0.0% | 100% | 98% | 98% |
Surg 1 | 62.5 ± 19.4% | 74.6 ± 7.5% | 31% | 92% | 73% | 75.0 ± 17.3% | 88.5 ± 5.5% | 55% | 95% | 86% | 70.8 ± 18.2% | 96.9 ± 3.0% | 81% | 95% | 93% |
Surg 2 | 66.7 ± 18.9% | 72.3 ± 7.7% | 31% | 92% | 71% | 95.8 ± 8.0% | 97.7 ± 2.6% | 89% | 99% | 97% | 95.8 ± 8.0% | 100 ± 0.0% | 100% | 99% | 99% |
Path | 66.7 ± 18.9% | 83.8 ± 6.3% | 43% | 93% | 81% | 95.8 ± 8.0% | 97.7 ± 2.6% | 89% | 99% | 97% | 70.8 ± 18.2% | 100 ± 0.0% | 100% | 95% | 96% |
Res | 79.2 ± 16.2% | 74.6 ± 7.5% | 37% | 95% | 75% | 87.5 ± 13.2% | 97.7 ± 2.6% | 88% | 98% | 96% | 66.7 ± 18.9% | 99.2 ± 1.5% | 94% | 94% | 94% |
Sonog | 75.0 ± 17.3% | 75.4 ± 7.4% | 36% | 94% | 75% | 95.8 ± 8.0% | 97.7 ± 2.6% | 89% | 99% | 97% | 83.3 ± 14.9% | 100 ± 0.0% | 100% | 97% | 97% |
Agg | 69.0 ± 7.0% | 79.0 ± 2.6% | 38% | 93% | 78% | 92.9 ± 3.9% | 96.4 ± 1.2% | 83% | 99% | 96% | 80.4 ± 6.0% | 99.5 ± 0.5% | 96% | 97% | 97% |
Automated reader | 100.0 ± 0.0% | 97.7 ± 2.6% | 89% | 100% | 98% |
. | OCT . | QME . | OCT + QME . | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Reader . | Sensitivity . | Specificity . | PPV . | NPV . | Acc . | Sensitivity . | Specificity . | PPV . | NPV . | Acc . | Sensitivity . | Specificity . | PPV . | NPV . | Acc . |
Eng 1 | 62.5 ± 19.4% | 85.4 ± 6.1% | 44% | 93% | 82% | 100.0 ± 0.0% | 97.7 ± 2.6% | 89% | 100% | 98% | 87.5 ± 13.2% | 100 ± 0.0% | 100% | 98% | 98% |
Eng 2 | 70.8 ± 18.2% | 86.9 ± 5.8% | 50% | 94% | 84% | 100.0 ± 0.0% | 97.7 ± 2.6% | 89% | 100% | 98% | 87.5 ± 13.2% | 100 ± 0.0% | 100% | 98% | 98% |
Surg 1 | 62.5 ± 19.4% | 74.6 ± 7.5% | 31% | 92% | 73% | 75.0 ± 17.3% | 88.5 ± 5.5% | 55% | 95% | 86% | 70.8 ± 18.2% | 96.9 ± 3.0% | 81% | 95% | 93% |
Surg 2 | 66.7 ± 18.9% | 72.3 ± 7.7% | 31% | 92% | 71% | 95.8 ± 8.0% | 97.7 ± 2.6% | 89% | 99% | 97% | 95.8 ± 8.0% | 100 ± 0.0% | 100% | 99% | 99% |
Path | 66.7 ± 18.9% | 83.8 ± 6.3% | 43% | 93% | 81% | 95.8 ± 8.0% | 97.7 ± 2.6% | 89% | 99% | 97% | 70.8 ± 18.2% | 100 ± 0.0% | 100% | 95% | 96% |
Res | 79.2 ± 16.2% | 74.6 ± 7.5% | 37% | 95% | 75% | 87.5 ± 13.2% | 97.7 ± 2.6% | 88% | 98% | 96% | 66.7 ± 18.9% | 99.2 ± 1.5% | 94% | 94% | 94% |
Sonog | 75.0 ± 17.3% | 75.4 ± 7.4% | 36% | 94% | 75% | 95.8 ± 8.0% | 97.7 ± 2.6% | 89% | 99% | 97% | 83.3 ± 14.9% | 100 ± 0.0% | 100% | 97% | 97% |
Agg | 69.0 ± 7.0% | 79.0 ± 2.6% | 38% | 93% | 78% | 92.9 ± 3.9% | 96.4 ± 1.2% | 83% | 99% | 96% | 80.4 ± 6.0% | 99.5 ± 0.5% | 96% | 97% | 97% |
Automated reader | 100.0 ± 0.0% | 97.7 ± 2.6% | 89% | 100% | 98% |
Abbreviations: Acc, accuracy; Agg, aggregate; Eng, engineer; Surg, surgeon; Path, pathology scientist; Res, medical resident; Sonog, sonographer.
Also reported in Table 3, the automated reader for QME resulted in 100.0%, 97.7%, and 98.1% sensitivity, specificity, and accuracy, respectively.
Discussion
This study is the first to determine the accuracy of QME, or any variant of optical elastography, for assessment of tumor margins in specimens excised during BCS. Prior to this study, recent investigations had demonstrated that QME has the potential to delineate tumor in breast tissue based on elevated elasticity on a micro-scale (43, 44). In these studies, as in this one, the contrast between tumor and healthy tissue based on changes in elasticity was shown to be complementary to the structural contrast provided by the underlying modality, OCT. However, prior studies had not investigated the ability of QME to identify positive margins in specimens from BCS. This study has demonstrated high accuracy of QME (96%) compared with OCT (78%) for detecting cancer within 1 mm of the margin in BCS specimens. The measurements here, on intact BCS specimens immediately following surgery, are also of high clinical relevance, compared with prior OCT studies that relied primarily on dissected mastectomy specimens for testing novel margin assessment techniques (27, 28, 46).
The QME results here indicate that tissue elasticity is an accurate predictor of malignancy. This trend was consistent across multiple tumor types in this study (e.g., solid and mucinous invasive carcinomas, and DCIS), although this must be confirmed with larger sample sizes of each tumor type. Future data collection will also continue to elucidate the sources of elevated elasticity in each tumor type, for example, deposition of stiff desmoplastic stroma, changes in collagen alignment, or interstitial fluid pressure. Better understanding of these mechanisms may foster even greater precision in the use of QME to detect various tumor types intraoperatively. This goal may be facilitated by implementing ultrahigh resolution QME to achieve cellular-scale elasticity resolution (53, 54).
Perhaps counterintuitively, although combining OCT and QME improved overall accuracy compared with OCT alone, the sensitivity using the high elasticity criterion (93%) was greater than that using QME and OCT criteria together (80%). This may be due to variability in reader interpretation of the qualitative OCT criteria (interpretation of structures) versus little variability in interpretation of the quantitative QME criterion (elasticity threshold). In particular, analysis of reader notes revealed variability in whether readers considered features as “continuous” with depth based on OCT (Fig. 1A), which was requisite in this study for an ROI to be considered positive (the engineers who participated as readers had prior experience in reading OCT and elastography images, which may have contributed to their slightly superior reading performance). On the other hand, specificity of OCT and QME together (∼100%) was greater than that of QME (96%) or OCT (79%) alone, indicating that the combination of information reduces false positives (e.g., mistaking benign, fibrous stroma for cancer). To improve intraoperative margin assessment and ultimately reduce the need for costly repeat surgeries, high sensitivity is essential to ensuring that less cancer is missed (i.e., minimizing false negatives). High specificity, on the other hand, ensures that less benign tissue is unnecessarily removed, and, while important to avoid overtreatment of the cancer, could be considered secondary to the need for high sensitivity in margin assessment. QME showed the best overall performance on both sides, but future studies will need to analyze the amount of repeat surgeries avoided versus extra volumes of tissue removed.
Although the accuracy of OCT reported here is similar to that reported in another study in BCS specimens (30), further refinement and quantification of OCT criteria may potentially increase accuracy in future studies. For instance, quantification of the rate of OCT signal attenuation with depth has shown potential for distinguishing benign and malignant regions (28, 32), although it is not clear the extent to which this will improve accuracy. Our preliminary automation of QME reading could also be extended to incorporate both QME and OCT-based information into an algorithm utilizing neural networks.
Although the spatial resolution of OCT and QME images is on the order of micrometers (Supplementary Fig. S4), which is comparable or superior to many existing intraoperative margin assessment techniques, readers were directed to ignore features smaller than approximately 1 mm. This set an effective “diagnostic resolution” of approximately 1 mm, and further study is needed to assess the accuracy of OCT and QME for smaller features, such as very early stage DCIS. However, at a more advanced stage, DCIS tends to result in dilation of ducts, and even small ducts can result in a fibrous stromal response in the surrounding tissue (52); thus, QME is expected to be sensitive to the associated mechanical changes caused by these features. The resolution of the technique may also improve in the future using the inverse method to solve for elasticity (55, 56), rather than rely on assumptions of stress uniformity with depth, as in our current method; however, computation times currently prohibit this from being implemented in clinically feasible timeframes.
Although the results presented here demonstrate the potential of QME for margin assessment in BCS, additional development is required to address the exclusion criteria described in the Results section. In particular, ROIs were excluded because of thermal damage, insufficient contact between the specimen, silicone layer and the imaging window, imaging artifacts and insufficient elasticity data overlaid on solid tissue. In the case of thermal damage, future studies could mitigate its presence in the images by using combined contrast from QME and OCT, including investigation of textural analysis, as well as colorimetric data from photographs simultaneously taken of the specimen (Supplementary Fig. S1), which often show black marks on the tissue surface that may correspond to thermal damage.
Two margins per specimen were imaged in this study, at a rate of approximately 9 minutes per margin. In some specimens, postoperative histology revealed cancer within 1 mm on margins that were not scanned, leading to a discrepancy in the overall prevalence of close margins reported in Table 1 (45%) and the prevalence of ROIs containing cancer (15.6%). Additional details on prevalence are provided in Supplementary Table S2. To facilitate rapid intraoperative QME of all radial margins, which is desirable in the clinical scenario, higher acquisition speed has been achieved using a novel approach (57), and may be further increased by implementation of high-speed OCT systems, an order of magnitude faster than that in this study (58). In addition, reader evaluation of the images took 90 seconds on average, which is conducive to clinical timeframes. However, to enable entire margins to be analyzed intraoperatively will require the development of image processing algorithms that facilitate more rapid image interpretation. Finally, a handheld probe is in development to allow QME assessment of tissues directly within the surgical cavity and to facilitate implementation into the clinical workflow (59).
In this study, analogously to previous studies (27, 30), we selected ROIs (n = 154) on which the reader study was performed rather than performing the study on entire margins. Of utmost importance to our reader study is the ability to correspond reader ratings to a ground truth. Without ground truth histology correspondence to our images, we have no basis to classify regions as tumor or nontumor. As the histology protocol used at our hospital enabled histology matching on only a small subset of each margin, the ROI approach was required. As this study has demonstrated strong correspondence between high elasticity and corresponding co-registered tumor, a next step is to perform a reader study on entire margins.
In summary, QME, an emerging OCT elastography technique that probes tissue elasticity on the microscale, has demonstrated detection of close and positive margins in freshly excised specimens in BCS with 93% sensitivity and 96% specificity. The use of intrinsic tissue contrast without need for exogenous dyes and the optimal trade-off in speed, field of view, and resolution provided by QME make it a promising candidate for improving intraoperative guidance of BCS. More broadly, QME may be applicable in a range of surgical or preoperative biopsy guidance applications, particularly in cancers that are known to exhibit altered mechanical properties.
Disclosure of Potential Conflicts of Interest
K.M. Kennedy is a consultant at OncoRes Medical. L. Chin has ownership interest (including patents) in OncoRes Medical. A. Curatolo is a consultant at OncoRes Medical Pty and has ownership interest (including patents) in OncoRes Medical Pty. S.L. Chin is an advisor at OncoRes Medical. K. Giles is the CEO at OncoRes Medical and has ownership interest (including patents) in OncoRes Medical. B. Latham has ownership interest (including patents) in OncoRes. B.F. Kennedy is the chief scientific officer at, reports receiving a commercial research grant from, and has ownership interest (including patents) in OncoRes Medical. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: K.M. Kennedy, W.M. Allen, B. Latham, C.M. Saunders, B.F. Kennedy
Development of methodology: K.M. Kennedy, R. Zilkens, W.M. Allen, L. Chin, H. DeJong, C.M. Saunders, B.F. Kennedy
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): K.M. Kennedy, R. Zilkens, W.M. Allen, Q. Fang, L. Chin, R.W. Sanderson, H.E.I. Tan, B. Kunjuraman, C. Yeomans, S.L. Chin, H. DeJong, B.F. Dessauvagie, B. Latham, C.M. Saunders, B.F. Kennedy
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): K.M. Kennedy, R. Zilkens, W.M. Allen, K.Y. Foo, L. Chin, P. Wijesinghe, A. Curatolo, N. Morin, C. Yeomans, S.L. Chin, K. Giles, B.F. Dessauvagie, B. Latham, C.M. Saunders, B.F. Kennedy
Writing, review, and/or revision of the manuscript: K.M. Kennedy, R. Zilkens, W.M. Allen, K.Y. Foo, Q. Fang, L. Chin, J. Anstie, P. Wijesinghe, A. Curatolo, H.E.I. Tan, B. Kunjuraman, S.L. Chin, H. DeJong, K. Giles, B.F. Dessauvagie, B. Latham, C.M. Saunders, B.F. Kennedy
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): R. Zilkens, W.M. Allen, L. Chin, H.E.I. Tan, C. Yeomans
Study supervision: C.M. Saunders, B.F. Kennedy
Acknowledgments
We gratefully acknowledge Dr. Max Bulsara from University of Notre Dame Australia for discussion of biostatistics. We acknowledge funding from OncoRes Medical and the Department of Health Western Australia.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.