Imaging system
We have converted an inverted microscope to a DUV-FSM to image the surfaces of fresh tissues from breast surgical specimens. A schematic of the DUV-FSM system is shown in Figure 1. A 285 nm LED (M285L4, Thorlabs, Newton, NJ) is mounted on the right side of an inverted fluorescence microscope (EXI-310, Accu-scope, Commack, NY) for oblique back-illumination. A 325 nm short-pass filter (XUV0325, Asahi Spectra, Torrance, CA) is placed in front of the LED to block emission spectrum tails in the visible range, avoiding possible overlap with fluorescence signals. A fused silica ball lens (model #67-388, Edmund Optics, Barrington, NJ) is used as a condenser to converge the LED radiation for a smaller illumination field and improved power density. The LED, short pass filter and ball lens are mounted inside a lens tube. A 3D printed arm holds the lens tube and is mounted on an optical post to allow for easy adjustment of the LED height and illumination angle so that the illumination area is slightly larger than the field-of-view (FOV) of the 4x microscope objective. Once the position of the LED was optimized, the entire system was fixed on an optical breadboard. To image a lumpectomy specimen, the specimen is placed on one of its six margins in a 70 mm diameter quartz dish to minimize autofluorescence of the glass. The quartz dish is mounted on a robotic, stepper-motor controlled XY stage custom designed for fast mosaic imaging (ABĒMIS LLC, Cleveland, OH). The excitation/emission filter block of the microscope is switched to the empty position so that the fluorescent signals of multiple fluorophores can be captured by a color camera without having to switch emission filters during the imaging process. A cooled, USB3.0 camera (MTR3CCD06000KPA, Hangzhou ToupTek Photonics Co., Ltd, Hangzhou, China) was selected for its large image sensor and pixel size, very low dark noise, and high image transfer speed, which are very important for fast image acquisition in intraoperative margin assessment. The camera has 2748 × 2200 pixels, pixel size of 4.54 µm and active area of 14.6 × 12.8 mm2. A 4x apochromatic long working distance objective lens with a numerical aperture of 0.13 was selected as a compromise between good lateral resolution (2~3 µm) and a large effective imaging area of 3.48 x 2.78 mm. The FOV of the objective lens is slightly larger than the imaging area of the camera to avoid distortions at the edge of the FOV. The microscope is housed inside a dark enclosure to prevent personnel exposure to DUV light and to eliminate background from room light.
Figure 1 The DUV-FSM margin imaging system: (A) Principle schematic of the system and (B) The simplified imaging light path. A 45 mW, 285 nm LED with a short-pass filter and a fused silica ball lens has been added to the inverted fluorescence microscope for fluorescence excitation. Fluorescence emission is collected by a Plan Fluor 4x (numerical aperture = 0.13) objective. The excitation/emission filter block of the microscope is switched to the empty position. Large specimen mosaic scanning is achieved by a motorized XY stage.
Breast tissue sample preparation
Thirty-seven fresh human breast tissues obtained from both lumpectomy and mastectomy specimens were acquired from the Medical College of Wisconsin (MCW) Tissue Bank. Information about the tissue subtype, number of samples, and surface area are provided in Table 1. Specimens were grossly examined and procured by Tissue Bank staff (MP), then placed in phosphate-buffered saline (PBS) solution, transported on ice to the research lab immediately, and stored in a refrigerator. Samples were imaged during the same day of excision. Propidium iodide (PI, P21493, Thermo Fisher Scientific, Waltham, MA) was used for nuclear staining and eosin Y (EY, 230251-25G, Sigma-Aldrich, St. Louis, MO) for staining of cytoplasm and connective tissues. Both PI and EY can be effectively excited at 285 nm. PI has a fluorescence emission in the yellow-red spectral range and EY has an emission in the green-yellow spectral range. For staining, PI and EY were dissolved in PBS (pH 7.2) to a concentration of 100 µg/ml and 1.0 mg/ml, respectively. Each specimen was stained in the PI solution for 1 minute, then in the EY solution for 20 seconds, and finally rinsed in PBS for 10 seconds. Once staining was completed, the specimen was placed onto the quartz plate of the specimen holder. A wide pallet knife was used to gently flatten the tissue against the quartz plate to remove air bubbles between the tissue and plate. Once the tissue was in the correct position, excess liquid was removed from the edges using a Kimwipe.
Table 1 Tissue types, number of samples, surface area, and number of patches (2 mm × 2 mm)
Tissue type
|
No. of Samples
|
Surface Area (cm2) (Median)
|
Patches/Sample (Median)
|
Total Patches
|
IDC
|
18
|
0.44 - 5.5 (1.3)
|
3 - 48 (19)
|
358
|
ILC
|
4
|
2 - 3.9 (2.8)
|
4 - 55 (15.5)
|
90
|
Adipose-rich Normal
|
3
|
3.9 - 5.29 (4.8)
|
48 - 70 (58)
|
176
|
Fibrous/Glandular-rich Normal
|
12
|
2.4 - 9 (4.4)
|
36 - 141 (84)
|
1,005
|
Imaging protocol
The specimen holder loaded with a tissue specimen was immobilized on the motorized XY stage after a specimen size measurement by a caliper. The focal plane was set at the bottom surface of the specimen. Mosaic images were collected with conservative overlapping regions of 0.75 mm in the X direction and 0.60 mm in the Y direction for a tradeoff between speed and stitching accuracy. The number of scanning steps are decided by the specimen size, effective imaging area and overlapping region dimension. The temperature of the camera was set to -18 ºC to reduce the electronic noise level. All images of each specimen were captured with a constant exposure time ranging from 50 to 100 ms, depending on the sample tissue type. The image acquisition and motor movements were controlled by a customized software developed in Microsoft Visual C# .NET. Image files were saved in TIFF format with 2748 x 2200 pixels. After imaging, the raw images were transformed to hue-saturation-lightness (HSV) color space images. The open source image processing package Fiji (fiji.sc/) was used to process the tissue images. A Fiji plugin named BaSiC [45] was applied to the saturation and lightness channels to correct for background and shadings caused by uneven and tilted illumination. The color space transform is necessary to preserve the original color information during illumination correction. After transforming back to red-green-blue (RGB) color space, image stitching was performed using a Fiji plugin developed by Preibisch et al.[46] Lastly, histogram equalization was applied to the R and G color channels of the stitched image to enhance the visual contrast.
Histopathology evaluation
Routine histopathology was used for final diagnosis of the tissue samples. Fereidouni et al. has previously shown PI and EY staining does not interfere with downstream histopathology processes.[42] Following DUV-FSM imaging, tissue specimens were returned to MCW Tissue Bank for formalin-fixed paraffin-embedded (FFPE) tissue processing. In order to obtain full face sections for histologic evaluation, an average cut depth of ~200 µm into the embedded tissue block was used during microtomy. The tissue sections were transferred to glass slides and stained with H&E. All slides were digitalized by a Panoramic 250 Flash II slide scanner (3DHistech Ltd., Budapest, Hungary). An unblinded qualitative side by side comparison of the H&E and DUV-FSM images was performed by an experienced breast pathologist (JMJ).
Visual inspection of DUV images
Visual inspection of DUV images was performed by three trained non-pathologists to evaluate the accuracy of non-pathologists to differentiate cancer from non-cancer tissue. The 37 breast tissue samples were divided into two groups: a training and a test group. The training group included 3 invasive carcinomas (2 IDC and 1 ILC) and 2 normal tissues (1 fibrotic and 1 adipose-rich breast sample), while the test set included 3 ILC, 16 IDC, 2 adipose-rich and 11 non-adipose-rich normal samples. Three non-medical inspectors (TGS, DHY and AE) who were blinded to pathological diagnosis were trained by the pathologist (JMJ) and imaging engineer (TL) during a one-hour session to visually identify the diagnostically useful features (such as adipose, ducts, cell density, infiltration, etc.) in the training DUV images using the associated H&E images. After training, each inspector was provided DUV images of samples in the test group without access to correlative H&E images. The inspectors interpreted DUV images and provided a diagnosis (invasive carcinoma vs. normal) for each of the test samples.
Quantitative image analysis
Quantitative analysis was applied to DUV tissue images to extract diagnostically useful parameters that may be useful for detecting positive tumor margins of lumpectomy specimens during BCS. Previous studies have shown that breast cancer cells have irregular cell size and shape, enlarged nuclei, and increased N/C.[47, 48] In this study, we investigated the feasibility of using N/C as a biomarker to differentiate invasive carcinoma from normal breast parenchyma at the surface of the tissue samples. Tumor region(s) on the stitched DUV image was outlined based on the corresponding H&E image. Since PI-stained cell nuclei primarily emit lights in the yellow-red wavelength range, only the R-channel of the stitched images in RGB color space was extracted and used to calculate the N/C.
The process for N/C calculation is illustrated in Figure 2. First, the color image (A&E) was converted to a R-channel image (not shown). Segmentation of the R channel image was implemented by combining edge detection and intensity thresholding.[49] The edge detection Sobel operator with adaptive threshold detects the edge information, while intensity thresholding eliminates textures caused by other features. The intensity threshold was set to 70% of the 1% brightest pixels in the image based on our experience. The segmented full R image (B&F) was divided into a set of small patches of 250 µm x 250 µm (or 198 x 198 pixels) in size, and N/C was calculated for each small patch (C&G). Then, N/Cs of neighboring 8 x 8 small patches were averaged and merged to form a large patch of 2 mm x 2 mm in size (D&H). Larger patches were used to reduce the sensitivity to small features, such as lobules, ducts and blood vessels. In general, larger patch size results in lower spatial resolution and sensitivity for cancer cell detection, but also lower false positive rate. To minimize the effect of tissue patches at the boundaries, large patches with more than half (or 32) small patches that have no cells (i.e., N/C = 0) were excluded in further analysis. A window size of 2 mm x 2 mm was selected to match the spatial resolution of standard breast pathology which samples at a step of 2 mm.
Figure 2 Demonstration of N/C calculation process with an IDC (A-D) and a normal sample (E-H). (A&E) A 4 mm x 4 mm area from the fluorescence images of an IDC and a normal sample, respectively. (B&F) The binarized image after nuclei segmentation from the same region as in (A&E). (C&G) The N/C image calculated with small patches (250 µm x 250 µm). (D&H) 8 x 8 small patches in (C&G) are binned to form large patches (2 mm x 2 mm).
The large patches were manually classified into adipose-rich, non-adipose-rich normal, ILC and IDC in accordance with H&E images. The number of large patches per tissue sample and total number of patches for each tissue subtype are presented in Table 1. This resulted in a total of 1,629 large patches of N/C images, including 358 patches from 18 IDC, 90 patches from 4 ILC, 176 patches from 3 adipose-rich and 1,005 patches from 12 non-adipose-rich normal tissues. All patches were used for the following comparison and classification studies. The mean N/C was compared between the 4 tissue subtypes (IDC, ILC, adipose-rich, non-adipose-rich normal groups) and between invasive (IDC, ILC) and normal tissue using Generalized Estimating Equations (GEE) to account for repeated observations per sample. Tukey’s adjustment was used for multiple pairwise comparisons.[50] For classifications, ROC curves were constructed using patch-level N/C to predict invasive versus normal tissue, IDC versus ILC among invasive samples, and adipose-rich versus non-adipose-rich tissue among normal samples. The Youden Index, which weighs false positive and false negative errors equally, was used to determine the cutoff point for the calculation of patch-level sensitivity and specificity in differentiating invasive and normal tissue. The analysis was done using SAS 9.4 (SAS Institute, Cary, NC).