Longitudinal Morphological and Physiological Monitoring of Three-dimensional Tumor Spheroids Using Optical Coherence Tomography

Summary

Optical coherence tomography (OCT), a three-dimensional imaging technology, was used to monitor and characterize the growth kinetics of multicellular tumor spheroids. Precise volumetric quantification of tumor spheroids using a voxel counting approach, and label-free dead tissue detection in the spheroids based on intrinsic optical attenuation contrast, were demonstrated.

Abstract

Tumor spheroids have been developed as a three-dimensional (3D) cell culture model in cancer research and anti-cancer drug discovery. However, currently, high-throughput imaging modalities utilizing bright field or fluorescence detection, are unable to resolve the overall 3D structure of the tumor spheroid due to limited light penetration, diffusion of fluorescent dyes and depth-resolvability. Recently, our lab demonstrated the use of optical coherence tomography (OCT), a label-free and non-destructive 3D imaging modality, to perform longitudinal characterization of multicellular tumor spheroids in a 96-well plate. OCT was capable of obtaining 3D morphological and physiological information of tumor spheroids growing up to about 600 µm in height. In this article, we demonstrate a high-throughput OCT (HT-OCT) imaging system that scans the whole multi-well plate and obtains 3D OCT data of tumor spheroids automatically. We describe the details of the HT-OCT system and construction guidelines in the protocol. From the 3D OCT data, one can visualize the overall structure of the spheroid with 3D rendered and orthogonal slices, characterize the longitudinal growth curve of the tumor spheroid based on the morphological information of size and volume, and monitor the growth of the dead-cell regions in the tumor spheroid based on optical intrinsic attenuation contrast. We show that HT-OCT can be used as a high-throughput imaging modality for drug screening as well as characterizing biofabricated samples.

Introduction

Cancer is the second leading cause of death in the world¹. Developing drugs targeting cancer is of crucial importance for patients. However, it is estimated that more than 90% of new anti-cancer drugs fail in the development phase because of a lack of efficacy and unexpected toxicity in clinical trials². Part of the reason can be attributed to the use of simple two-dimensional (2D) cell culture models for compound screening, which provide results with limited predictive values of compound efficacy and toxicity for the following stages of drug discovery^2,3,4. Recently, three-dimensional (3D) tumor spheroid models have been developed to provide clinically relevant physiological and pharmacological data for anti-cancer drug discovery^{3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25}. Since these spheroids can mimic tissue-specific properties of tumors in vivo, such as nutrient and oxygen gradient, hypoxic core as well as drug resistance¹⁹, the use of these models can potentially shorten drug discovery timelines, reduce costs of investment, and bring new medicines to patients more effectively. One critical approach to evaluating compound efficacy in 3D tumor spheroid development is to monitor the spheroid growth and recurrence under treatments^9,26. To do this, quantitative characterizations of the tumor morphology, involving its diameter and volume, with high-resolution imaging modalities, are imperative.

Conventional imaging modalities, such as bright-field, phase contrast^7,9,22,24, and fluorescence microscopy^8,9,16,18,22 can provide a measurement of the spheroid's diameter but cannot resolve the overall structure of the spheroid in 3D space. Many factors contribute to these limitations, including penetration of the probing light in the spheroid; diffusion of the fluorescent dyes into the spheroid; emitting fluorescent signals from excited fluorescent dyes inside or on the opposite surface of the spheroid due to strong absorption and scattering; and depth-resolvability of these imaging modalities. This often leads to an inaccurate volume measurement. Development of the necrotic core in spheroids mimics necrosis in in vivo tumors^{6,10,15,19,25}. This pathological feature is unlikely reproduced in 2D cell cultures^19,25,27,28. With a spheroid size larger than 500 µm in diameter, a three-layer concentric structure, including an outer layer of proliferating cells, a middle layer of quiescent cells, and a necrotic core, can be observed in the spheroid^{6,10,15,19,25}, due to lack of oxygen and nutrients. Live and dead cell fluorescence imaging is the standard approach to label the boundary of the necrotic core. However, again, penetrations of both these fluorescent dyes and visible light hinder the potential to probe into the necrotic core to monitor its development in its actual shape.

An alternative 3D imaging modality, optical coherence tomography (OCT) is introduced to characterize the tumor spheroids. OCT is a biomedical imaging technique that is capable of acquiring label-free, non-destructive 3D data from up to 1-2 mm depths in biological tissues^{29,30,31,32,33,34}. OCT employs low-coherence interferometry to detect back-scattered signals from different depths of the sample and provides reconstructed depth-resolved images at micron-level spatial resolutions in both lateral and vertical directions. OCT has been widely adopted in ophthalmology^35,36,37 and angiography^38,39. Previous studies have used OCT to observe the morphology of in vitro tumor spheroids in basement membrane matrix (e.g., Matrigel) and evaluate their responses to photodynamic therapy^40,41. Recently, our group established a high-throughput OCT imaging platform to systematically monitor and quantify the growth kinetics of 3D tumor spheroids in multi-well plates⁴². Precise volumetric quantification of 3D tumor spheroids using a voxel counting approach and label-free necrotic tissue detection in the spheroids based on intrinsic optical attenuation contrast were demonstrated. This paper describes the details of how the OCT imaging platform was constructed and employed to obtain high-resolution 3D images of tumor spheroids. The step-by-step quantitative analyses of the growth kinetics of 3D tumor spheroids, including accurate measurements of spheroid diameter and volumes, is described. Also, the method of the non-destructive detection of necrotic tissue regions using OCT, based on the intrinsic optical attenuation contrast is presented.

Protocol

1. Preparation of Cells

1. Obtain cell lines from a qualified supplier.
   NOTE: Verify that cells from the cell lines of interest can form spheroid in the culture media or with the help of a substrate (basement membrane matrix like Matrigel). Look into the literature⁹ or perform one round of a pre-experiment for a check.
2. Thaw the frozen cells following the specific procedure provided by the cell-line supplier. A general procedure can be found elsewhere⁴³.
3. Culture the cells for 1-2 passages in 25 cm² culture flasks. The cells are then ready to use for 3D cell culture.
4. Monitor the health status of the cells every day and maintain them in an incubator under standard conditions (37 °C, 5% CO₂, 95% humidity). Refresh the media as needed.
   NOTE: The culture medium consists of DMEM (4.5 g/L glucose), 1% antibiotic-antimycotic, 10% fetal bovine serum. Subculture cells before they reach confluence in the culture flask. Follow the cell culture guideline provided by the supplier. A general procedure can be found elsehwhere⁴⁴.
5. Perform 3D cell culture in multi-well plates based on the following general protocol⁹.
   1. Remove the culture media from the culture flask and wash it with sterilized phosphate buffered saline (PBS, heated to 37 °C).
   2. Resuspend the cells by adding 1 mL of trypsin ethylenediaminetetraacetic acid (EDTA, 0.5%) into the flask for 3 min. Then, add culture media to dilute the trypsin.
   3. Transfer the cell suspension into a 15 mL centrifuge tube and centrifuge for 5 min at 500 x g and room temperature.
   4. Remove the supernatant and resuspend cells with 4 mL of pre-warmed culture medium. Pipette one drop of sample onto a hemocytometer for cell counting to determine cell concentration. Dilute the cells to appropriate concentration for seeding (e.g., 3,000 cells/mL).
      NOTE: Optimize the initial cell concentration of the spheroid for each cell-line and each type of multi-well plate (96-well, 384-well or 1536-well).
   5. Seed cells into an ultra-low attachment (ULA) round-bottomed multi-well plate. Add 200 µL of cells suspension into each well at the concentration of 3,000 cells/mL so that each well has about 600 cells.
   6. At RT, centrifuge the whole plate using a plate adapter for 7 min, right after seeding, at a speed of 350 x g or the lowest speed available.
      NOTE: The centrifuge helps gather cells to the center of the well to facilitate forming a single, uniform spheroid. The centrifuge step is performed only once at the beginning to form the tumor spheroids. It will not be repeated when the tumor spheroids start growing.
   7. Maintain the multi-well plate at 37 °C and 5% CO₂ in a culture incubator and refresh the culture media every 3 days.
      NOTE: Growth time may vary for different 3D culture conditions. In our study, 3,000 cells/mL is used for both U-87 MG and HCT 116 cell lines in 96-well plates, so that the spheroid can grow to ~500 μm in 4‒7 days for HCT 116 cells. Consider adding media supplements and growth factors for different spheroid models, based on the general 3D culture protocol.
   8. Perform OCT imaging of tumor spheroids every 3‒4 days for a longitudinal study of their growth.
      NOTE: Recommended time points for OCT imaging would be day 4, day 7, day 11, day 14, day 18 and day 21.

2. High-throughput OCT Imaging Platform

NOTE: See referenced work^{29,30,31,32,33,34} for a thorough review of principles and applications of OCT. See Figure 1 and Huang et al.⁴² for details of the custom OCT imaging system used in this study.

1. Choose an appropriate broadband light source for the OCT system for tumor spheroid imaging.
   NOTE: Here, a superluminescent diode (SLD, Figure 1A,B) with a central wavelength of ~1,320 nm and ~110 nm bandwidth was used as a broadband light source.
2. Construct the reference arm and sample arm of the OCT system following the schematics (See Figure 1A,B for details). See the Table of Materials for a list of optical components to construct the OCT system. Ensure that the optical path length of the reference arm and sample arm are closely matched.
3. Construct the spectrometer, including a collimator, a grating, an F-theta lens and a line-scan camera (See Figure 1C for setup³⁴) for details of spectrometer design of OCT. Alternatively, select a commercial spectrometer that matches the center wavelength of the light source. Make sure that the spectrometer is aligned correctly to cover the entire laser bandwidth, to achieve high photon collection efficiency and to provide slow wash-out of the interference pattern.
4. Characterize the performance of the OCT system, including the following metrics such as sample arm power, total imaging depth, depth-dependent sensitivity, axial resolution, depth of focus and lateral resolution. Place a weak reflector (e.g., a mirror with a neutral density filter) as a sample to measure the depth-dependent sensitivity, axial resolution, and depth of focus. Place a USAF resolution test chart target as the sample to check the lateral resolution.
   NOTE: See references^34,45 for definitions of metrics of OCT performance and protocols to characterize these metrics⁴⁵. See Table 1 for a list of measured parameters for the custom OCT system used in our study.
5. Select a motorized translation stage to provide horizontal movement of the multi-well plate to image tumor spheroids in different wells (See Figure 1B). Use a stage with a travel range larger than 108 mm x 72 mm to ensure a full scanning of all the wells of the multi-well plate. Use a 2D or 3D motorized translation stage with software control to enable precise location of each well and automation of the OCT system for high-throughput imaging.
6. Use a plate adapter or design a plate holder (by 3D printing) to hold the multi-well plate in a fixed position.
7. Correct the tilting and rotation of the multi-well plate using a 2D tilting stage and a rotation stage mounted on the translational stage (See Figure 1D), before conducting any OCT imaging to minimize variation of the focus plane from different wells. Use D2, D11, B6, D6, G6 as the guiding wells when monitoring their relative positions in the OCT images (Figure 1A).
8. Adjust the rotation of the plate to ensure the edges of the plate are parallel with the direction of stage movement so that the wells remain at the same horizontal positions in the OCT images (Figure 1E). Adjust the tilting of the plate to be parallel to the optical table so that the wells remain at the same vertical locations for OCT imaging (Figure 1F).
   NOTE: Adjustment of the tilting angle and focus help optimize the OCT image quality for all the wells. However, variations of the height of culture media in different wells may cause changes in optical paths which may lead to defocusing of the spheroid image. Auto-focus may be implemented to control the focal plane of OCT imaging to achieve optimized image quality. The adjustment step does not resolve poor OCT image quality of the tumor spheroid due to the following issues: the spheroid decentering due to initial seeding location; spheroid elevation when embedded in biofabricated extracellular matrices; poor plate quality with large variations of the height of well bottoms. Additional software control with auto-focus or self-alignment functions can be implemented to optimize the performance of the OCT imaging system.
9. Use a custom computer program to control the OCT image acquisition and the stage movement to collect data from each well sequentially.

3. OCT Scanning and Processing of Tumor Spheroids

1. On the day of the OCT imaging of tumor spheroids, take the multi-well plate from the incubator. Transfer the multi-well plate under the OCT imaging system. Place it on top of the plate adapter.
   NOTE: OCT imaging of tumor spheroids may be performed with the polystyrene plate lid on or off. However, the water condensations on the lid due to evaporation from the wells may affect light transmission and distort the light path, yielding less optimal OCT images from the spheroids.
2. Adjust the height of the plate by moving along the z-direction of the translation stage. Maintain the focal plane position at ~100–200 μm below the top surface of each spheroid, to minimize the effect of the non-uniform depth-wise focal profile.
3. Set a proper OCT scanning range (e.g., 1 mm x 1 mm) in the custom software to cover the whole tumor spheroid according to its development stages. Click Save Parameters to save the setting.
4. Use the custom software to acquire 3D OCT images of tumor spheroids one by one for all the wells of the plate containing spheroids. Click the Preview button to view the preview image and click the Acquire button to acquire the OCT image.
   NOTE: Ensure that the OCT spheroid data are collected when the stage is not in motion. The spheroid is usually located at the center of the U-bottom well. However, the spheroid may be shifted in the culture media when the stage is accelerating or decelerating due to the inertia of the spheroid in the culture media.
5. Process 3D OCT datasets of tumor spheroids to generate OCT structural images with a custom C++ processing code. See Figure 2A for a flowchart of post-processing of OCT data.
   NOTE: See Figure 4A for the generated 3D OCT structural images.
   1. See Chapter 5 of Drexler and Fujimoto³⁴ and Jian et al.⁴⁶ for detail descriptions of post-processing steps of OCT data. Calibrate the pixel size in all three dimensions. Re-scale the OCT structural images on corrected scales.
      NOTE: The distance in the axial direction (z-direction) of OCT images is a measure of the optical path difference between the reference arm and sample arm. Thus, the refractive index of the sample (n) needs to be taken into consideration when calibrating the pixel size in the axial direction for rescaling. In our study, we use n = 1.37 as the refractive index of the tumor spheroid⁴².
6. Generate the collage of spheroid images using 2D OCT Images in three cross-sectional XY, XZ, and YZ planes across the centroid of the spheroid. See Figure 3C–E for the representative output of collages of spheroid images. Perform image registration for all the spheroids, using the MATLAB function dftregistration⁴⁷, to ensure that the centroids of all the spheroid are located approximately at the same location.
7. Obtain 3D rendering of the spheroid using a commercial or custom software.
   NOTE: The following steps show how to obtain the 3D rendering of tumor spheroids using a commercial software.
   1. Load the 3D OCT data into the software.
   2. Click the Surpass panel. Then, click Add New Volume. Choose the Blend mode to use for 3D rendering.
   3. Adjust the viewing angle by dragging the image using the mouse pointer.

4. Morphological Quantification of 3D Tumor Spheroids

NOTE: A custom written code in MATLAB processes this quantification. Click the Run button to initiate the process. See Figure 2B for the flowchart of the steps of morphological quantification of spheroids.

1. Quantify spheroid diameter, height, and diameter-based volume.
   1. Select 2D OCT Images in three cross-sectional XY, XZ, and YZ planes that cross the centroid of the spheroid.
   2. Measure the diameter and height of the spheroid in XY and XZ planes, respectively.
   3. Calculate diameter-based spheroid volume using: , with a presumption of the spherical shape of tumor.
2. Quantify voxel-based spheroid volume.
   1. Apply a 3D averaging filter on the OCT structural data of spheroid to remove speckles.
   2. Segment tumor spheroids using the Canny edge detection⁴⁸ filter, frame by frame, with a proper threshold separating the tumor spheroid region from the well bottom.
   3. Group connective voxels for 3D data (see built-in function: bwconncomp).
   4. Calculate the mean distance between each connective voxel in the group and the spheroid centroid (manually chosen), for each group. Identify the spheroid region as the group with the minimum mean distance.
   5. Count the number of voxels within the spheroid region and then multiply by the actual volume of an individual voxel (volume/voxel), yielding the total volume of the spheroid.

5. Dead-Cell Region Detection of 3D Tumor Spheroids

NOTE: In a homogeneous medium, OCT back-scattered intensity detected as a function of depth (I(z)) can be described by the Beer-Lambert Law⁴⁹: , where z represents the depth, μ is the optical attenuation coefficient, and I₀ is the incident intensity to the sample. Hence the derived optical attenuation coefficient can be expressed as: . Since OCT images are often plotted on a logarithmic scale, the slope of the OCT intensity profile can be retrieved to derive the optical attenuation coefficient. See Figure 2C for a flowchart of the generation of the optical attenuation maps.

1. Perform segmentation to remove unwanted regions outside the spheroid. Perform 3D average filter to suppress the speckle noise that is inherent in OCT images.
2. Obtain pixel-wise optical attenuation coefficients by linear fitting the log-scale OCT intensity profile over a certain depth range (moving window), extract its slope, and multiply the slope by -1/2.
   NOTE: The attenuation coefficient at each voxel within the segmented spheroid region is calculated based on the slope of OCT intensity profile in a 10-voxel depth window (~40 μm in depth), with the voxel located in the middle of the window.
3. Apply the methods above (steps 5.1 and 5.2) to each axial scan in a frame and each frame in a 3D dataset containing the segmented spheroid region until optical attenuation coefficients for all voxels of the segmented spheroid region are calculated.
4. Perform the binary thresholding to highlight the high-attenuation region.
   NOTE: See Huang et al.⁴² for the determination of threshold of the high-attenuation region using histogram analysis.
5. Highlight the binarized optical attenuation map on the original image to label the dead-cell region (blending). Generate the 3D-rendered image of the blended attenuation map to visualize the 3D distribution of the dead-cell region.

6. Histology and Immunohistochemistry

NOTE: Histology and immunohistochemistry (IHC) stained images of tumor spheroids are obtained to correlate with the corresponding OCT results.

1. At selected time points, select 1-2 tumor spheroids from the multi-well plate for histology and IHC staining. Use a pipette with the 1 mL pipette tips to transfer the spheroid from the well to a 1.5 mL centrifuge tube.
   NOTE: Cut the 1 mL pipette tip before transfer to ensure that the opening of the tip is larger than the size of the tumor spheroid to avoid damaging the structure of the spheroid.
2. Collect each tumor spheroid in a single 1.5 mL microcentrifuge tube filled with 10% formaldehyde and fix for 48 h.
3. Perform the histology and IHC processes for each spheroid, using standard paraffin embedding techniques.
   NOTE: Stain 5 μm thick sections of tumor spheroids for hematoxylin and eosin (H&E) and terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) apoptosis detection. A counterstaining of hematoxylin is applied to TUNEL. A digital slide scanner was used to scan the stained sample and obtain high-resolution histological and IHC images.

Results

High Throughput Optical Coherence Tomography Imaging of Spheroids in a 96-well Plate

Figure 3 exhibits the result of HT-OCT scanning of a 96-well plate with HCT 116 tumor spheroids on Day 3. The sequential scan of the whole plate starts from the bottom-right well (H12). Figure 3B shows the flow chart of the software implementation of the HT-OCT system. After one spheroid data were collected and processed, the plate would move to next well, wait for ~2 s to allow the spheroid to rest, and collect the next spheroid data. Each OCT data consist of 400 x 400 x 1024 voxels, which corresponded to an actual volume of 1.0 x 0.84 x 2.3 mm³. Figure 3C shows a collage of en-face OCT images of HCT 116 spheroids generated from the processed data. The result is comparable with images from other 2D high-throughput imaging system²². Given the 3D imaging capability of the OCT, we could also generate the collage of 2D cross-sectional spheroid images from 96 wells (Figure 3D) to monitor spheroid heights and visualize spheroid inhomogeneity in the vertical direction. A collage of 3D-rendered spheroid images is also feasible from any predefined angle (Figure 3E) to visualize the overall 3D shape and evaluate the roundness of the spheroid. Note that the overall OCT imaging and process time for the whole 96-well plate would be ~21 min and ~25 min when the line-scan camera is running at a speed of 92 kHz and 47 kHz, respectively. See Video 1 for an example.

Longitudinal Morphological and Physiological Monitoring of the Tumor Spheroid

After we obtained the OCT structural images of tumor spheroids from the plate for multiple time points, we could further analyze these data by quantifying the morphological and physiological information of the tumor spheroids. Figure 4 shows the different approaches to characterize tumor spheroids and obtain longitudinal morphological and physiological information from them.

Figure 4B shows different ways to visualize the tumor spheroid. With the aid of either commercial or free software, we could load the 3D data into the software and create a "volume" of the tumor spheroid (3D rendering), which shows the overall structure of tumor spheroid in 3D space. With proper thresholding, we could generate a surface plot of tumor spheroid (Figure 4B), which could be used to segment the spheroid and measure the volume. We could also generate the orthogonal slides (ortho slides) from different cross-section planes in different orientations (Figure 4B, XZ, YZ, and XY) and measure the diameter and height of the tumor spheroid from these ortho slides.

Gathering the OCT data of the same spheroid from multiple time point, we could quantify the morphological information and generate the growth curve of the spheroid to show its longitudinal changes. Figure 4C shows representative data of an HCT 116 tumor spheroid being monitored for 21 days. From the segmented data and ortho slides, we measured the diameter, height and voxel-based volume of the spheroid for all the time point, which were listed in the table. We also calculated the diameter-based volume for a comparison. The growth curves in size and volume were plotted, respectively. From the growth curves, we could see that this HCT 116 tumor spheroid followed a linear growth pattern in volume before day 11. Before this time point, the spheroid kept growing and maintained a relatively uniform shape. However, after day 11, the spheroid became disrupted, flattened and fully collapsed on day 21. The growth curve of voxel-based volumes clearly shows the trend, with a gradually decreased volumes after day 11.

Based on the OCT data, we can also obtain the physiological information of the distribution of dead-cells within the tumor spheroids by analyzing the pixel-by-pixel optical attenuation from 2D cross-sectional images. Following the methods illustrated in Figure 2 and Protocol 5, we could quantitatively determine the dead-cell regions and monitor the growth of these regions as a function of time. Figure 4D shows a representative result of longitudinal tracking of the increase of dead-cell areas in the tumor spheroid. The areas highlighted in red, which had high optical attenuation, show the labeled necrotic areas. From the 3D rendered optical attenuation map during the 14 day development, we could see the red sector expanding, indicating the increase of the necrotic regions. As the percentage of the necrotic areas increased, the tumor spheroid could not maintain its perfect shape. Therefore, they would tend to disrupt and collapse, which were seen in the longitudinal monitoring of tumor morphology in Figure 4C.

The proposed nondestructive dead tissue region detection technique was verified by comparing OCT optical attenuation map of HCT 116 tumor spheroid with corresponding images obtained by histology and IHC. Figure 4D presents such a comparison with a Day 14 HCT 116 spheroid. A good match between the OCT attenuation map and corresponding H&E and TUNEL slices were found, which was indicated by analyzing the features within the regions in H&E and TUNEL slices marked by dash lines derived from the contour of OCT high attenuation regions. In H&E slices, the dead tissue regions were indicated by less dense and aggregated structure located within the dashed line region. In TUNEL slices, a good match was observed between high attenuation region and TUNEL-labeled apoptotic cellular region.

Figure 1: Construction of a high-throughput optical coherence tomography (HT-OCT) system for tumor spheroid imaging. (A) Schematics of the HT-OCT system. A diagram of the 96-well plate is plotted next to the OCT system. Five wells (D2, D11, B6, D6, G6) labeled in yellow are used for the fine adjustment of the stages in (D). (B) The actual configuration of HT-OCT system. See Table of Materials for optical components used for each part of the system. (C) Spectrometer design for the HT-OCT system. (D) Stage setup for the HT-OCT system. Proper alignment of the 6-axis stage and synchronization between the OCT acquisition and the stage movement are required for high-throughput imaging. (E) and (F) show the effects of rotation and tilting on final image of different wells. Rotation causes the OCT images of different wells to shift horizontally while tilting will lead to vertical shifting of different wells. Please click here to view a larger version of this figure.

Figure 2: Data Processing for OCT images of tumor spheroids. (A) Flowchart of general post-processing steps for OCT data. (B) Flowchart of morphological quantification of the tumor spheroid. (C) Flowchart of dead cell region detection of the tumor spheroid. Scale bar: 100 µm for all the subfigures. Please click here to view a larger version of this figure.

Figure 3: High-throughput OCT scanning of a 96 well plate containing U-87 MG tumor spheroids. (A) The actual setup with the 96-well plate under the objective. (B) Flow Chart of the software implementation of HT-OCT system. Collages of 96 en face (C), cross-sectional (D) and 3D rendered maximum intensity projection (MIP) (E) OCT images of Day 3 HCT 116 spheroids were generated from the processed data. Scale Bar: 200 µm for all the subfigures. Please click here to view a larger version of this figure.

Figure 4: Longitudinal Morphological and Physiological Quantification of Tumor Spheroids with 3D OCT data. (A) Obtained 3D OCT structural images of a tumor spheroid after general OCT post-processing. From the OCT data, we can generate a 3D surface plot and XZ, YZ and XY orthogonal slices to visualize the structure of the tumor spheroid in any direction (B). We can perform longitudinal monitoring of a single tumor spheroid (C), characterizing its diameter, height and voxel-based volume (listed in the Table of Materials) and plotting the growth curves in size and volume during the 21-day development. In the example, as the spheroid developed, it became disrupted on day 11 and fully collapsed on day 21. We can further monitor the physiological status of a tumor spheroid longitudinally based on the optical intrinsic attenuation contrast (D). 3D rendered images of a tumor spheroid showed the appearance and growth of dead-cell regions from day 7 to day 14. The high-attenuation-labeled dead-cell areas in red were matched with histological and immunohistochemical (IHC) results. OCT attenuation map, H&E, and TUNEL result in Figure 4D are modified from Ref. 42. Scale bars: 100 µm for all the subfigures. Please click here to view a larger version of this figure.

Video 1: High-throughput OCT imaging of tumor spheroids. A workflow of 3D OCT imaging, basic OCT processing and stage movement was presented in the video with a 5x speed. Previews of processed OCT structural images of spheroids were also presented. Please click here to view this video. (Right-click to download.)

Discussion

Tumor activity is highly relevant to its morphological structure. Similar to monitoring characteristic growth curve for 2D cell cultures, tracking the growth curve for 3D tumor spheroids is also a conventional approach to characterize the long-term spheroid growth behavior for different cell lines. Notably, we can characterize the drug response by analyzing tumor degradation or tumor regrowth directly reflected in the growth curve. Therefore, quantitative assessment of 3D tumor spheroids, including the size and volume, to derive the growth curve, is of great importance for the characterization of tumor spheroids and the evaluation of compound effect. Currently, imaging platforms based on bright field, phase contrast or fluorescent imaging have been established for routine imaging and analysis of morphology or functions of the 3D tumor spheroids^8,9,18,22. However, they are unable to resolve the entire, large tumor structure due to limited depth penetration as well as low-resolution depth-resolvability. In the representative results, we have demonstrated OCT to visualize the entire 3D structure of the tumor spheroid developing over time. 3D OCT imaging could provide the view of the spheroid in any orientation and any cross-section with high-resolvability, which was not available in conventional imaging modalities that lack the resolution along the depth. Furthermore, voxel-based volume quantification based on 3D OCT data yielded an accurate quantification of spheroid volumes without assuming their original shapes. Therefore, we have demonstrated that OCT is a robust imaging modality for 3D morphology characterization of tumor spheroids, which ensures accurate measurements of characteristic growth patterns for different cell lines and potentially can serve as an alternative for drug response evaluation.

Viability tests using fluorescent staining remain a popular approach for functional analyses of tumor spheroids, especially for drug screening¹⁸. However, the disruptive nature of fluorescent dyes indicates that these tests are only suitable for end-point studies. In our representative results (Figure 4D), we demonstrated an alternative method that can characterize cell viability within the entire spheroid. Our results have shown that OCT could distinguish the dead-cell region from the viable region in the spheroid based on intrinsic optical attenuation contrast. In addition, with 3D imaging capability and non-destructive nature of the OCT system, quantitative evaluation of the dead-cell distributions and in situ monitoring of the progression of dead-cell regions within the spheroid are feasible, which potentially provide more valuable information of the spheroid growth pattern. However, we should note that, in our representative results, we are not able to differentiate different types of cell death modes, such as apoptosis and necrosis, in the binary OCT attenuation map.

Since a drug compound library can be extensive (>10,000), a high-throughput and robust system to characterize tumor spheroids in multi-well plates during drug screening is imperative. The current high-throughput imaging system can achieve a screening of the whole 96-well plate in <5 min in 2D scan mode²². OCT can be adapted for high-throughput screening purpose, with the aid of a motorized stage. One can also obtain a commercially available OCT system (See Table of Materials for a list of commercial OCT systems) with a similar performance to our custom OCT system, and incorporate a motorized stage into the system. However, efforts must be taken to modify the commercial OCT system to integrate the motorized stage. Also, custom software implementation to realize the synchronization between the OCT acquisition trigger and stage movement trigger is required. For our prototype HT-OCT system, it took 2-18 seconds to acquire one 3D OCT data from a single tumor spheroid, depending on the choice of camera speed. Thus, the total acquisition time can be as short as ~3.2 min for a 96-well plate using HT-OCT system. However, the intermediate steps for the current HT-OCT system, including data processing, reading and writing data on hard drives, and stage movements, remained time consuming. Additional ~18 min would be needed on top of the ~3.2 min minimum data acquisition time. The total imaging time can be further reduced in several aspects: use state-of-the-art OCT systems equipped with a high-speed tunable laser source^50,51; optimized the workflow by arranging critical steps (data acquisition, data processing, writing, stage movement) working in parallel; employ a parallel OCT imaging with a space-division multiplexing setup⁵². With system optimization, the high-throughput OCT system can be utilized in cancer drug discovery as well as characterization of other 3D bio-fabricated samples (e.g., 3D tissue organoids) for various biomedical applications.

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Custom Spectral Domain OCT imaging system			Developed in our lab
Superluminescent Diode (SLD)	Thorlabs	SLD1325	light source
2×2 single mode fused fiber coupler, 50:50 splitting ratio	AC Photonics	WP13500202B201
Reference Arm
Lens Tube	Thorlabs
Adapter	Thorlabs
Collimating Lens	Thorlabs	AC080-020-C
Focusing Lens	Thorlabs
Kinematic Mirror Mount	Thorlabs
Mirror	Thorlabs
1D Translational Stage	Thorlabs
Continuous neutral density filter	Thorlabs
Pedestrial Post	Thorlabs
Clamping Fork	Thorlabs
Sample Arm
Lens Tube	Thorlabs
Adapter	Thorlabs
Collimating Lens	Thorlabs	AC080-020-C
Galvanometer	Thorlabs
Relay Lens	Thorlabs	AC254-100-C	two Relay lens to make a telescope setup
Triangle Mirror Mount	Thorlabs
Mirror	Thorlabs
Objective	Mitutoyo
Pedestrial Post	Thorlabs
Clamping Fork	Thorlabs
Polarization Controller	Thorlabs
30mm Cage Mount	Thorlabs
Cage Rod	Thorlabs
Stage
3D motorized translation stage	Beijing Mao Feng Optoelectronics Technology Co., Ltd.	JTH360XY
2D Tilting Stage
Rotation Stage
Plate Holder			3D printed
Spectrometer
Lens Tube	Thorlabs
Adapter	Thorlabs
Collimating Lens	Thorlabs	AC080-020-C
Grating	Wasatch		G = 1145 lpmm
F-theta Lens	Thorlabs	FTH-1064-100
InGaAs Line-scan Camera	Sensor Unlimited	SU1024-LDH2
Name	Company	Catalog Number	Comments
Cell Culture Component
HCT 116 Cell line	ATCC	CCL-247
Cell Culture Flask	SPL Life Sciences	70025
Pipette	Fisherbrand	14388100
Pipette tips	Sorenson Bioscience	10340
Gibco GlutaMax DMEM	Thermo Fisher Scientific	10569044
Fetal Bovine Serum, certified, US origin	Thermo Fisher Scientific	16000044
Antibiotic-Antimycotic (100X)	Thermo Fisher Scientific	15240062
Corning 96-well Clear Round Bottom Ultra-Low Attachment Microplate	Corning	7007
Gibco PBS, pH 7.4	Thermo Fisher Scientific	10010023
Gibco Trypsin-EDTA (0.5%)	Thermo Fisher Scientific	15400054
Forma Series II 3110 Water-Jacketed CO2 Incubators	Thermo Fisher Scientific	3120
Gloves	VWR	89428-750
Parafilm	Sigma-Aldrich	P7793
Transfer pipets	Globe Scientific	138080
Centrifuge	Eppendorf	5702 R	To centrifuge the 15 mL tube
Centrifuge	NUAIRE	AWEL CF 48-R	To centrifuge the 96-well plate
Microscope	Olympus
Name	Company	Catalog Number	Comments
Histology & IHC
Digital slide scanner	Leica	Aperio AT2	Obtain high-resolution histological images
Histology Service	Histowiz		Request service for histological and immunohistological staining of tumor spheroid
Name	Company	Catalog Number	Comments
List of Commerical OCTs
SD-OCT system	Thorlabs	Telesto Series
SD-OCT system	Wasatch Photonics	WP OCT 1300 nm
Name	Company	Catalog Number	Comments
Software for Data Analyses
Basic Image Analysis	NIH	ImageJ	Fiji also works.
3D Rendering	Thermo Fisher Scientific	Amira	Commercial software. Option 1
3D Rendering	Bitplane	Imaris	Commercial software. Option 2. Used in the protocol
OCT acquisition software			custom developed in C++.
Stage Control	Beijing Mao Feng Optoelectronics Technology Co., Ltd.	MRC_3	Incorporated into the custom OCT acquisition code
OCT processing software			custom developed in C++. Utilize GPU. Incorporated into the custom OCT acquisition code.
Morphological and Physiological Analysis			custom developed in MATLAB