Name: high-throughput image analysis of tumor spheroids: a user-friendly software application to measure the size of spheroids automatically and accurately
Uploaded: 2014-07-08T23:30:00
Duration: 8 min 39 s
Description: Watch this Scientific Journal Video about High-throughput Image Analysis of Tumor Spheroids: A User-friendly Software Application to Measure the Size of Spheroids Automatically and Accurately at JoVE.

We would like to thank the Raymond and Beverly Sackler Foundation for their support on our research.

1. 3D tumor spheroids formation, drug treatments and image collection are performed as described in our previous paper15.
2. Software Installation
<ol>
	<li>Install licensed MATLAB&#160;software onto the computer used for image analysis. The following toolboxes from MATLAB&#160;are also required to be installed - Signal Processing toolbox, Image Processing toolbox, and Parallel Computing toolbox* (*needed for parallel computing mode only). 
	NOTE: Many universities purchase and maintain group licenses so that the software is free to be downloaded and used by affiliated scientists.</li>
	<li>Install SpheroidSizer program from SpheroidSizer.zip file (Please see supplemental files):
	<ol>
		<li>Save the zip file in your local file system.</li>
		<li>Unzip the SpheroidSizer file.</li>
		<li>Save the files in a designated directory/folder, which will be subsequently referred to as the &#8220;Installation Directory&#8221;. 
		NOTE: SpheroidSizer has been extensively tested on a Windows 7 operating system. It is expected to work on alternative operating systems with minimal adjustments (not tested).</li>
	</ol>
	</li>
</ol>
3. Preparation for the Image Analysis by SpheroidSizer
<ol>
	<li>Determine the image scale/resolution of the imaging system (absolute scale of the image in microns per pixel (&#956;m/pix)). 
	NOTE: If the size of each pixel on the camera chip is known, the image scale can be computed as objective magnification x pixel size (&#956;m/pix). This value can be obtained from the imaging software equipped with the microscope as embedded metadata or with help from the imaging system vendor. This value will be required in Step 4.6.</li>
	<li>Convert any proprietary image file formats to the accepted file formats - TIFF, JPEG, and other common image file formats.</li>
	<li>Name image files and arrange directories (Figure 5A). 
	NOTE: The software relies on the proper layout of directory structure and filenames to reformat the results into plate format:
	<ol>
		<li>Name the image files in the following format: [plate name]_[row][column].[extension] or [plate name][space][row][column].[extension]. [row] follows the alphabetical order and [column] follows the numerical order.&#160; 
		NOTE: Freely available automatic batch renaming software can be found to aid users in this step.</li>
		<li>Arrange the directory/folder based on the experiment in the following way: Each experiment should have one directory. Under each experiment&#8217;s directory, there should be subdirectories for each time point. Under each time point subdirectory, there should be all the images from all plates. 
		NOTE: In order for the analysis results to be optimally sorted in the formatted results, we suggest keeping the same number of digits for each identifier by filling 0&#8217;s on the left, e.g., time points are named as 000H, 072H, and 144H.</li>
	</ol>
	</li>
</ol>
4. Image Analysis of the Spheroids by SpheroidSizer
<ol>
	<li>Open MATLAB, then open the &#8220;Command Window&#8221;, type cd &#8216;[Installation Directory]&#8217; and press [RETURN].</li>
	<li>Type &#8220;SpheroidSizer1_0&#8221; in the &#8220;Command Window&#8221; and press [RETURN] to launch the SpheroidSizer program.</li>
	<li>Click the &#8220;Browse&#8221; button in the SpheroidSizer1.0 window to select the experiment directory that contains all the images.</li>
	<li>Select &#8220;Include Subfolders&#8221; toggle beneath the &#8220;Folder&#8221; text field to process multiple nested image folders under the designated directory. 
	NOTE: If &#8220;Include Subfolders&#8221; toggle is not selected, only the images directly under the directory are processed, and all of the subfolders are ignored.</li>
	<li>Select the &#8220;On-the-fly display&#8221; option to display each segmented image on top of its source image for quality control as the computation is being executed. 
	NOTE: Computation speed is faster if &#8220;On-the-fly display&#8221; option is not chosen.</li>
	<li>Specify &#8220;Resolution&#8221; (image scale/resolution in &#956;m/pix) of the analyzed images in the box, for the program to correctly convert the spheroid&#8217;s measurements from pixel to &#956;m. 
	NOTE: All the images in the same folder or analyzed together should be taken under the same microscope with the same objective so that the image scale/resolution remains constant for each experiment.</li>
	<li>(Optional) users can follow Step 5 for Advanced User-defined Settings.</li>
	<li>Click &#8220;Compute&#8221; to start the computation. 
	NOTE: The software conducts an automatic filename check before proceeding to the computation. If a dialog box shows up indicating - &#8220;Error exists in filename&#8221;, click &#8220;Exit and show list of errors&#8221; and fix the errors in listed filenames (See Step 3.3). Then, click &#8220;Compute&#8221; to start the computation again.</li>
	<li>Click &#8220;Pause&#8221; button to pause the computation; and computation can be resumed by clicking the same button showing &#8220;Continue&#8221;. 
	NOTE: The &#8220;Results Table&#8221; displays the &#8220;Folder&#8221;, &#8220;File&#8221;, &#8220;Volume&#8221; (in mm3), &#8220;Length&#8221; (in &#956;m), &#8220;Width&#8221; (in &#956;m), and &#8220;Valid&#8221; (check box) for all of the analyzed spheroids (Figure 5C). Volume is calculated based on the measured major axis (length) and minor axis (width) (V=0.5 * Length * Width2). The &#8220;Valid&#8221; check box is an option for the user to choose if the analysis of the image is valid or invalid after quality control, see Step 6.</li>
</ol>
5. Advanced User-defined Settings
<ol>
	<li>Click &#8220;Advanced&#8221; button in the SpheroidSizer1.0 window to bring up the Advanced Configurations window in order to adjust user-defined settings (Figure 5B).</li>
	<li>Enter the filenames of interest in the &#8220;Format Output&#8221; and &#8220;List Output&#8221; boxes in the Advanced Configurations window under &#8220;Input&#8221;.</li>
	<li>Enter a number from &#8220;2-10&#8221; in the &#8220;Reduce&#8221; box. This is a coefficient for the software to reduce the image size in computation in order to improve the computation speed. The larger the number is, the faster the computation speed. The default &#8220;Reduce&#8221; is set to 10.</li>
	<li>Enter the to-be-processed image file extension in the &#8220;Include type&#8221; box.</li>
	<li>Enter the image file extensions or endings that are not going to be processed by the program in the &#8220;Exclude type&#8221; box as follows: &#8220;_crude.jpg&#8221;.</li>
	<li>Choose &#8220;None&#8221; for &#8220;Special Color&#8221; to process 8 bit and 16 bit color images properly; choose &#8220;12 bit&#8221; for &#8220;Special Color&#8221; to process 12-bit color images properly</li>
	<li>Check &#8220;Use Parallel Computing&#8221; if the computer used for image analysis is equipped with multiple CPUs and/or multi-core CPUs. If this is true, then go to Step 5.7.1; if not, then skip Steps 5.7.1 and 5.7.2. 
	NOTE: An error will occur if the computer being used does not support the chosen configuration.
	<ol>
		<li>Check the &#8220;Use Parallel Computing&#8221; option in the Advanced Configurations window. 
		NOTE: Only use the parallel computing mode when 4 or more cores are available for the computer.</li>
		<li>Enter a number from 4-12 in the &#8220;#Workers&#8221; (computing cores) box. 
		NOTE: This number has to be equal to or less than the number of computing cores in user&#8217;s computer. A maximum of 12 is imposed by the MATLAB&#160;parallel computing toolbox supports the maximum of 12 cores. When parallel computing is being executed, a small dialogue box displays asking the user to wait for the parallel computation to finish; the computation cannot be paused, and nor is the &#8220;On-the-fly display&#8221; feature being executed in parallel computing mode.</li>
	</ol>
	</li>
</ol>
6. Quality Control
<ol>
	<li>Click the corresponding cell in the &#8220;Results Table&#8221; to confirm the accurate boundary contour of a spheroid in the analyzed images, 
	NOTE: The original and the quality control images will show up on the right side for review. The user can examine all the images sequentially using the down-arrow on the keyboard.</li>
	<li>Refine the boundary of the spheroid on the selected image using the following two tools, if needed:
	<ol>
		<li>Click the &#8220;Manual Initialize&#8221; button to display the original image. Then click-and-hold the mouse right outside the spheroid and drag the ellipse tool to cover the spheroid on the original image. 
		NOTE: The active-contour algorithm initiates using the user-submitted contour and executes to converge on the desired spheroid outline. The &#8220;Results Table&#8221; will be automatically updated with the new results. The &#8220;Manual Initialize&#8221; tool allows the user to manually provide initialization for active-contour.</li>
		<li>Click the &#8220;Hand Draw&#8221; button to display the original image. Then use the mouse or a touch enabled screen to precisely draw the boundary of the spheroid. 
		NOTE: This outline is directly measured to generate major and minor axes, which are updated in the &#8220;Results Table&#8221;. The &#8220;Hand Draw&#8221; tool is used only when the &#8220;Manual Initialize&#8221; tool fails to converge on the desired boundary of the spheroid.</li>
	</ol>
	</li>
	<li>Uncheck the checkbox in the &#8220;Valid&#8221; column in the corresponding row of the &#8220;Results Table&#8221;, when an image does not contain any valid spheroid upon inspection. An &#8220;Invalid&#8221; label appears at the top-left corner of the quality-control image. If &#8220;Valid&#8221; is unchecked, the values of all the measurements are empty for the spheroid in the exported formatted and output results files. 
	NOTE: the following keyboard shortcuts are available to use in the &#8220;Result Table&#8221;: &#8220;down arrow&#8221; for next image; &#8220;v&#8221; for valid/invalid; &#8220;m&#8221; for &#8220;Manual Initialize&#8221; tool and &#8220;h&#8221; for &#8220;Hand Draw&#8221; tool.</li>
</ol>
7. Saving and Exporting Data
<ol>
	<li>Click the &#8220;Export Study&#8221; button in the SpheroidSizer1.0 window to export the intermediate state of the analysis, if the user needs to exit the software before finishing the project. Specify the name and directory of the file to be saved.</li>
	<li>Click &#8220;Import Study&#8221; button to bring back the above intermediate state result from &#8220;Export Study&#8221; and continue to work on it. 
	NOTE: The intermediate state files are in a native MATLAB&#160;format (.mat) and are not directly readable by any other software programs. A safety feature built into the software makes automatic exports of the open project in case the program exits unintentionally. When needed, the user can find this file, whose name starts with &#8220;~tmp&#8221; that contains the corresponding time stamp in the [Installation Directory].</li>
	<li>Click &#8220;Format Results&#8221; in the SpheroidSizer1.0 window to save the results. 
	NOTE: Two forms of results are saved in the experiment&#8217;s directory. The exported filenames can be configured in the Advanced Configurations window (see Protocol Step 5.2). The format output file is a tab-delineated table that organizes the volume value into the original plate format in the order of an ascending plate number for each time point; and all the time points are organized in an ascending order (Figure 5D). The list output file is a tab-delineated table that contains all of the measurements in the form of ordered lists (Figure 5E).</li>
</ol>

<ol>
	<li>Hamilton, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Multicellular+spheroids+as+an+in+vitro+tumor+model.">Multicellular spheroids as an in vitro tumor model.</a> Cancer Lett. 131, 29-34 (1998).</li><li>Sutherland, R. M., McCredie, J. A., Inch, W. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Growth+of+multicell+spheroids+in+tissue+culture+as+a+model+of+nodular+carcinomas.">Growth of multicell spheroids in tissue culture as a model of nodular carcinomas.</a> J Natl Cancer Inst. 46, 113-120 (1971).</li><li>Inch, W. R., McCredie, J. A., Sutherland, R. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Growth+of+nodular+carcinomas+in+rodents+compared+with+multi-cell+spheroids+in+tissue+culture.">Growth of nodular carcinomas in rodents compared with multi-cell spheroids in tissue culture.</a> Growth. 34, 271-282 (1970).</li><li>Hirschhaeuser, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Multicellular+tumor+spheroids:+an+underestimated+tool+is+catching+up+again.">Multicellular tumor spheroids: an underestimated tool is catching up again.</a> J Biotechnol. 148, 3-15 (2010).</li><li>Friedrich, J., Seidel, C., Ebner, R., Kunz-Schughart, L. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Spheroid-based+drug+screen:+considerations+and+practical+approach.">Spheroid-based drug screen: considerations and practical approach.</a> Nat Protoc. 4, 309-324 (2009).</li><li>Yuhas, J. M., Li, A. P., Martinez, A. O., Ladman, A. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+simplified+method+for+production+and+growth+of+multicellular+tumor+spheroids.">A simplified method for production and growth of multicellular tumor spheroids.</a> Cancer Res. 37, 3639-3643 (1977).</li><li>Ayers, G. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Volume+of+preclinical+xenograft+tumors+is+more+accurately+assessed+by+ultrasound+imaging+than+manual+caliper+measurements.">Volume of preclinical xenograft tumors is more accurately assessed by ultrasound imaging than manual caliper measurements.</a> J Ultrasound Med. 29, 891-901 (2010).</li><li>Carpenter, A. E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=CellProfiler:+image+analysis+software+for+identifying+and+quantifying+cell+phenotypes.">CellProfiler: image analysis software for identifying and quantifying cell phenotypes.</a> Genome Biol. 7, (2006).</li><li>Kamentsky, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Improved+structure,+function+and+compatibility+for+CellProfiler:+modular+high-throughput+image+analysis+software.">Improved structure, function and compatibility for CellProfiler: modular high-throughput image analysis software.</a> Bioinformatics. 27, 1179-1180 (2011).</li><li>Monazzam, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+new,+fast+and+semi-automated+size+determination+method+(SASDM)+for+studying+multicellular+tumor+spheroids.">A new, fast and semi-automated size determination method (SASDM) for studying multicellular tumor spheroids.</a> Cancer Cell Int. 5, 32 (2005).</li><li>Schneider, C. A., Rasband, W. S., Eliceiri, K. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=NIH+Image+to+ImageJ:+25+years+of+image+analysis.">NIH Image to ImageJ: 25 years of image analysis.</a> Nat Methods. 9, 671-675 (2012).</li><li>Chan, V. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Active+contours+without+edges.">Active contours without edges.</a> IEEE Trans Image Process. 10, 266-277 (2001).</li><li>Bernard, O., Friboulet, D., Thevenaz, P., Unser, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Variational+B-spline+level-set:+a+linear+filtering+approach+for+fast+deformable+model+evolution.">Variational B-spline level-set: a linear filtering approach for fast deformable model evolution.</a> IEEE Trans Image Process. 18, 1179-1191 (2009).</li><li>Kass, M. W. A., Terzopoulos, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Snakes:+Active+contour+models.">Snakes: Active contour models.</a> International Journal of Computer Vision. 1, 321-331 (1987).</li><li>Wong, C., Vosburgh, E., Levine, A. J., Cong, L., Xu, E. Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Human+neuroendocrine+tumor+cell+lines+as+a+three-dimensional+model+for+the+study+of+human+neuroendocrine+tumor+therapy.">Human neuroendocrine tumor cell lines as a three-dimensional model for the study of human neuroendocrine tumor therapy.</a> J Vis Exp. (66), (2012).</li></ol>

No conflicts of interest declared.

This study presents a fast, flexible, effective and automated program - SpheroidSizer for accurate determination of the size of 3D tumor spheroids. SpheroidSizer is easy to use and requires minimal user input. The most critical steps for accurate, smooth and successful running of SpheroidSizer include: that the spheroids are imaged in the center of the field without touching the edge of the well; all the files to be analyzed together as one project should be imaged under the same microscope with the same objective; all the files to be analyzed are correctly named and arranged as indicated in the protocol; and correct user-defined settings are entered before computation.
The advantages of SpheroidSizer include its ability to tolerate gradual background change in the image as well as to generate smooth contours which correspond to the general spherical shapes of the spheroids utilizing the active contour algorithm. Performance of active contour can be compromised in two situations: poor initialization, or presence of other local edges distracting from the desired contour. Specifically in our tested cases, the second situation sometimes happen when the necrotic core of a large spheroid attracts the active contour resulting in smaller contour being reported. It is worth noting that other automated threshold-based methods also suffer in this situation unless the threshold is specifically set by hand. The software therefore puts forward effort to help users detect and remedy compromised segmentation by providing easy quality-control features. If segmentation error happens from poor initialization, uses can use &#8220;Manual Initialize&#8221; tool to override the automated initialization. When image quality is too poor for active contour, users can easily &#8220;Hand Draw&#8221; the contour that feeds into quantification. Existing software such as CellProfiler can be adapted for this application in a semi-automated fashion. The workflow can be cumbersome when large amounts of images with different imaging conditions are presented or when a subset of images needs more human intervention to measure correctly. SpheroidSizer provides an all-in-one suite for computation and quality control to manage high-throughput image analysis workflow.
SpheroidSizer is currently limited to detecting one spheroid per image and only measures the axial length of the spheroid. The program can be extended to support further quantification needed by researchers such as quantification on spheroids with necrotic core, detecting multiple spheroids in one image or monitoring the shape of the spheroids. Furthermore, the program can be modified to detect and measure the size of the excised tumors from animal or human, which will certainly be beneficial for researchers when conducting in vivo pre-clinical or clinical research. Post-processing of the detected spheroids can also be investigated aiming at reducing human effort needed for quality control and further improving throughput. SpheroidSizer is a generalized image analysis application for 3D tumor spheroids that are produced from any cell types, and therefore can be used by a wide cancer research community.

Three-dimensional (3D) tumor spheroids are &#8220;spherically symmetric aggregates of tumor cells analogous to tissues, with no artificial substrate for cell attachment&#8221;1-3. The cytology and morphology of the tumor spheroids better mimics the in vivo tumor tissue organization and microenvironments than monolayer two-dimensional (2D) cells. 3D tumor spheroids have become a practical in vitro model for high-throughput screenings of anti-cancer therapeutic drugs or examining the efficacy of candidate drugs before in vivo animal or clinical testing4. Clinically, the efficacy of any anti-cancer drug treatment is evaluated based on reduced tumor growth. Analogously, spheroid volume can be used as a measure of the efficacy for in vitro cancer drug studies. Spheroid volume (V = 0.5 * Length * Width2) is determined based on the major and minor axial length (more commonly known as length and width) of the spheroids6,7. Most researchers have to manually draw the length and width on each spheroid, often using the software offered by microscopy companies and sold together with the imaging instruments. This technique becomes problematic when high-throughput drug screens are performed and more than hundreds of images are produced. Some recent studies reported the use of open source image analysis software toolboxes such as CellProfiler8-10 and ImageJ11 to develop rudimentary segmentation routines/macros which involved illumination correction and simple thresholding. These routines often need to be re-adjusted for different batches of images according to the illumination condition and the image contrast change; therefore, these software packages cannot meet the robustness requirement of the high-throughput image analysis. Friedrich and collaborators (2009) used proprietary software to measure the spheroid&#8217;s volume semi-automatically5. The method described in Monazzam and his colleagues&#8217; paper10 was a semi-automated method to measure spheroid&#8217;s size only for a small numbers of images. Therefore, there exists a clear need for robust, flexible, automated and ready-to-use image analysis tools for the 3D tumor spheroids.
In this study, we describe SpheroidSizer &#8211; a MATLAB-based and open-source software application to measure the size of tumor spheroids automatically and accurately. SpheroidSizer is designed to process many different batches of images of 3D spheroids in the same session. Utilizing the active contour algorithm12-14, SpheroidSizer can tolerate image contrast change, robustly ignore gradual change in background illumination and recognize spheroids in the image. It can also tolerate many usual artifacts, e.g., debris, originated from the specimen. The workflow is designed so that users can perform quality control during or after the computation. Manual overwrite of the analysis result can easily be performed on the spot. Taking advantage of the parallel computing toolbox, the analysis speed can be further boosted by coordinating multiple computing cores to work on the computation simultaneously on a user&#8217;s computer. Moreover, SpheroidSizer outputs the results in two different forms to allow easy interfacing with downstream analysis tools.

<table border="0" cellpadding="0" cellspacing="0" width="569">
	<colgroup>
		<col />
		<col />
		<col />
		<col />
	</colgroup>
	<tbody>
		<tr height="51">
			<td height="51" style="height:51px;width:123px;">Axiovert 200M inverted microscope</td>
			<td style="width:176px;">Carl Zeiss Microscopy, LLC</td>
			<td style="width:88px;"></td>
			<td style="width:183px;">microscope for imaging</td>
		</tr>
		<tr height="17">
			<td height="17" style="height:17px;width:123px;">Vostro 1720</td>
			<td style="width:176px;">Dell Inc.</td>
			<td style="width:88px;"></td>
			<td style="width:183px;">single-core regular laptop</td>
		</tr>
		<tr height="35">
			<td height="35" style="height:35px;">HP Z820</td>
			<td style="width:176px;">HP Inc.</td>
			<td style="width:88px;"></td>
			<td style="width:183px;">multi-core performance workstation</td>
		</tr>
		<tr height="35">
			<td height="35" style="height:35px;width:123px;">Matlab and Simulink R2013a</td>
			<td>Mathworks, Inc, Natick, MA</td>
			<td style="width:88px;"></td>
			<td style="width:183px;">Matlab software&#160;</td>
		</tr>
	</tbody>
</table>

high-throughput image analysis of tumor spheroids: a user-friendly software application to measure the size of spheroids automatically and accurately

The increasing number of applications of three-dimensional (3D) tumor spheroids as an in vitro model for drug discovery requires their adaptation to large-scale screening formats in every step of a drug screen, including large-scale image analysis. Currently there is no ready-to-use and free image analysis software to meet this large-scale format. Most existing methods involve manually drawing the length and width of the imaged 3D spheroids, which is a tedious and time-consuming process. This study presents a high-throughput image analysis software application &#8211; SpheroidSizer, which measures the major and minor axial length of the imaged 3D tumor spheroids automatically and accurately; calculates the volume of each individual 3D tumor spheroid; then outputs the results in two different forms in spreadsheets for easy manipulations in the subsequent data analysis. The main advantage of this software is its powerful image analysis application that is adapted for large numbers of images. It provides high-throughput computation and quality-control workflow. The estimated time to process 1,000 images is about 15 min&#160;on a minimally configured laptop, or around 1 min&#160;on a multi-core performance workstation. The graphical user interface (GUI) is also designed for easy quality control, and users can manually override the computer results. The key method used in this software is adapted from the active contour algorithm, also known as Snakes, which is especially suitable for images with uneven illumination and noisy background that often plagues automated imaging processing in high-throughput screens. The complimentary &#8220;Manual Initialize&#8221; and &#8220;Hand Draw&#8221; tools provide the flexibility to SpheroidSizer in dealing with various types of spheroids and diverse quality images. This high-throughput image analysis software remarkably reduces labor and speeds up the analysis process. Implementing this software is beneficial for 3D tumor spheroids to become a routine in vitro model for drug screens in industry and academia.

SpheroidSizer is designed to produce automated detection, delineation and measurement of the 3D spheroids, with remarkably reduced labor and acutely increased efficiency for large amounts of images. Figure 1A shows the workflow of SpheroidSizer. The core computation steps include automated initialization, active contour algorithm and contour quantification. After automated computation, the quality-control feature uses a combination of &#8220;Manual Initialize&#8221; and &#8220;Hand Draw&#8221; tools to salvage any imperfect segmentation. Figure 1B illustrates the detailed automated active contour algorithm. The initialization step (0th iteration) utilizes basic image processing steps to generate approximate size and location of the spheroid and generate a spherical initiation contour with an estimated size. The initiation contour feeds into the active contour algorithm. In turn it iterates to adjust according to the local image gradient and shape curvature. The active contour algorithm finishes when the contour stabilizes (converges), i.e. 477 iterations for this image, or when the pre-defined maximum number of iterations is executed. In this example, the initialization contour is purposely enlarged to better showcase the algorithm. In reality, the initialization is usually very close to the actual boundary and much fewer iterations are needed for the algorithm to converge. Subsequently, the algorithm takes morphometric measurements of the detected spheroid boundary. The major and minor axes of the spheroid are measured using MATLAB&#160;image processing toolbox (Figure 1C). The major axis is defined as the line segment connecting a single pair of the farthest points on the contour, which is referred to length (L). The minor axis is defined as the longest line perpendicular to the major axis, which is referred to width (W). In this case, the values of L and W are very close since the spheroid is spherical. The volume of the spheroid is calculated as V = 0.5 *L * W2.
One of the features of SpheroidSizer is its automated detection of the boundary of the spheroids even on images with uneven or noisy background utilizing the active contour algorithm (Figures 2B-D). Computational processing of bright field images is often plagued by uneven background, which misleads adaptive thresholding-based methods to produce unwanted thresholding results. The issue is especially apparent when multi-well plates are used and walls of the wells may create shadowing effects on the images. However, because the active contour algorithm is not sensitive to the gradual shading change in background, it is able to identify spheroid segmentation in these bright field images with proper initialization. Figure 2 shows a few examples of images with uneven or noisy background, like uneven illumination (Figure 2B), debris (Figure 2C) or necrotic core (Figure 2D). With automated active contour algorithm, SpheroidSizer delineates these spheroids accurately in all these images as shown in the red contour at the lower panel of each figure.
The quality-control feature of SpheroidSizer is key to a high-throughput workflow. The &#8220;Manual Initialize&#8221; and &#8220;Hand Draw&#8221; tools are the valuable complimentary tools for this application. Among hundreds or thousands of images, it is unavoidable that the automated algorithm is not able to correctly detect the spheroids in some images. As illustrated in Figure 3A, when improper detection of the spheroid is caused due to the initialization step, i.e.&#160;improper size or location of the initiation contour in the image (top panel), the &#8220;Manual Initialize&#8221; tool works by allowing the user to properly define the location and size of the spheroid manually (lower panel). It triggers the active contour algorithm to initiate with manually defined contour and execute to converge on the desired outline. For these tough images like the original image in Figure 3B, the spheroid is located in a distracting and noisy background. SpheroidSizer is not able to properly identify the spheroid by the automated method (top panel) or by the &#8220;Manual Initialize&#8221; tool with proper initialization (middle panel). In this case, the &#8220;Hand Draw&#8221; tool can be used to manually draw the outline of the spheroid as illustrated in the lower panel. The program uses the user-defined boundary to measure the major and minor axes of the spheroid and calculate the volume. All the corrected results are immediately incorporated into the &#8220;Results Table&#8221; and can be accordingly exported.
To determine the performance of SpheroidSizer in larger data sets, we first compare operation time by analyzing the same set of 288 images using 1) manual measurements with microscope vendor supplied software; 2) SpheroidSizer with a single-core regular laptop; and 3) SpheroidSizer with a multi-core parallel computing performance workstation. The manual measurements follow our typical protocol prior to developing the software: the length and width of each spheroid are drawn by hand and measured using the vendor program (as seen the red lines in the top panel of Figure 4A); then the user copies down the values of the measurements. SpheroidSizer processes each image by generating the spheroid boundary (as shown the red outline in the lower panel of Figure 4A), measuring the axial major and minor length, and exporting the results in spreadsheets. As seen in Table 1, based on the calculation from 288 images, it takes an average of 31.67 sec&#160;to measure one spheroid per image manually; while it only takes SpheroidSizer less than 2 sec&#160;when running on a single-core regular laptop; and less than 1 sec&#160;when running on a 12-core performance workstation. Therefore, image analysis is over 18x&#160;faster per image using SpheroidSizer than manual measurements. It dramatically reduces labor when more than thousands of images are analyzed. Next, we compare the variability in measurements of the 24 spheroids shown in Figure 4A between manual measurements and SpheroidSizer. The 24 spheroids are measured three times by both methods; and the standard deviation of each individual spheroid is calculated. As seen in Figure 4B, the standard deviation from SpheroidSizer (green line and dots) is close to zero except for the three spheroids that are corrected at the quality-control step, which still show smaller standard deviation than those from the manual measurements method. All these indicate that SpheroidSizer performs the image analysis more efficiently and accurately.
We conducted a drug screen using human BON-1 3D tumor spheroids to find out which compounds in combination with an hsp90 inhibitor are the potential candidates for testing the anti-tumor effects in vivo. Human BON-1 3D tumor spheroids were grown on agarose-coated 96-well plates as described in the previous paper15. Eight different compounds with six serial dilutions plus media and vehicle were screened for their single and combinatorial effects with 10 nM and 20 nM hsp90 inhibitor in duplicates respectively. Two spheroids were used for each concentration of the individual compound or the combined compounds. Four 96-well plates with total 384 spheroids were used. All the spheroids were imaged at 0, 72, 144, 168, and 192 hr. A total of 1,920 images were produced from this experiment. It took SpheroidSizer only 30 min&#160;to complete the computational analysis of the 1,920 images with an additional 50 min&#160;for the quality control and data export. SpheroidSizer speeds up the image analysis process enormously. Figure 5A shows a screen shot of the folder arrangements and filenames for this experiment as an example for Protocol Step 3.3. Figures 5B-E shows the screen shots of the image analysis pop-up windows and results using SpheroidSizer as illustrations for Protocol Steps 4, 5, and 7. Taking the volumes of the 3D spheroids from the formatted results table exported from SpheroidSizer, we made graphs - growth of 3D tumor spheroids upon compound treatments versus time of treatments. Two representative graphs from this experiment are shown in Figure 5F and 5G. Figure 5F shows that the combined treatments of hsp90 inhibitor and cladribine (green line) inhibit the growth of 3D spheroids more than the single treatment of hsp90 inhibitor (purple line) or cladribine (orange line), suggesting that the combined treatments of hsp90 inhibitor and cladribine may have anti-tumor effects in vivo. Figure 5G shows that the combined treatments of hsp90 inhibitor and adriamycin (green line) does not inhibit the growth of 3D spheroids more than the single treatment of adriamycin (orange line) or hsp90 inhibitor (purple line), suggesting that the combined treatments of hsp90 inhibitor and adriamycin may not have anti-tumor effects in vivo. This experiment helped us better select the compounds to test their anti-tumor effects in vivo and SpheroidSizer is the key to the speedy experimental data analysis.
<img alt="Table 1" fo:content-width="6in" src="/files/ftp_upload/51639/51639table1.jpg" width="600" /> 
Table 1. Comparison of operation time on image analysis between manual measurements and SpheroidSizer when analyzing the same set of 288 images. <a href="https://www.jove.com/files/ftp_upload/51639/51639table1.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Figure 1" fo:content-width="5in" fo:src="/files/ftp_upload/51639/51639fig1highres.jpg" src="/files/ftp_upload/51639/51639fig1.jpg" /> 
Figure 1. SpheroidSizer - an open-source software application for measuring the spheroid&#8217;s size. A)&#160;The core workflow of the application. B)&#160;Illustration of the active contour algorithm at different stages of iteration. Please note that the initialization contour (iteration 0) was purposely enlarged in order to showcase the algorithm. C)&#160;The major and minor axial length measurements and calculation of the volume by SpheroidSizer. L - major axis: the line segment connecting a single pair of farthest points on the contour (referred to length); W - minor axis: the longest line perpendicular to the major axis (referred to width).
<img alt="Figure 2" fo:content-width="5in" fo:src="/files/ftp_upload/51639/51639fig2highres.jpg" src="/files/ftp_upload/51639/51639fig2.jpg" /> 
Figure 2. Representative results from the automated segmentation of SpheroidSizer, showing robustness against various image conditions. A)&#160;Typical good quality images. B)&#160;Images with different brightness and contrast. C)&#160;Images with distracting debris. D)&#160;Images of spheroids with necrotic core. Images at the top panel of each figure are the source/original images; images at the lower panel of each figure are the quality-control images; and the red outline is the spheroid segmentation drawn by automated computation.&#173;&#173;
<img alt="Figure 3" fo:content-width="5in" fo:src="/files/ftp_upload/51639/51639fig3highres.jpg" src="/files/ftp_upload/51639/51639fig3.jpg" /> 
Figure 3. Illustration of the &#8220;Manual Initialize&#8221; and &#8220;Hand Draw&#8221; tools. A)&#160;The &#8220;Manual Initialize&#8221; tool allows the drawing of a fitting ellipse shape across the spheroid for initialization, when inaccurate spheroid segmentation occurs after the automated initialization. B)&#160;The &#8220;Hand Draw&#8221; tool allows accurate hand drawing of the spheroid boundary, when inaccurate spheroid segmentations occur with both automated and manual initialization. The blue line around the spheroid shows the initialization contour; the red outline is the identified spheroid boundary. Please note that the spheroid in &#8220;Manual Initialize&#8221; in A) and the spheroid in &#8220;Hand Draw&#8221; in B) is purposely enlarged to better showcase the tools.
<img alt="Figure 4" fo:content-width="5in" fo:src="/files/ftp_upload/51639/51639fig4highres.jpg" src="/files/ftp_upload/51639/51639fig4.jpg" /> 
Figure 4. Comparison of image analysis performance between SpheroidSizer and manual measurements when analyzing the same set of 24 images. A)&#160;Representative spheroids to show how the length and width of spheroids are determined by manual measurements and SpheroidSizer. Top 24 images contain hand-drawn length/width of each spheroid in red lines using manual measurements; lower 24 images (the same 24 images) contain computer-drawn spheroid boundary in red outline using SpheroidSizer. B)&#160;Standard deviation of length or width from three measurements on each individual spheroid.
<img alt="Figure 5" fo:content-width="6in" src="/files/ftp_upload/51639/51639fig5highres.jpg" width="600" /> 
Figure 5. A representative example of utilization of SpheroidSizer in drug screen - image analysis on the spheroids&#8217; images that were collected from a drug screen using BON-1 3D tumor spheroids. A)&#160;A screen shot of the folder arrangements and filenames for this project. B)&#160;A screen shot of the advanced configurations window in SpheroidSizer. C)&#160;A screen shot of SpheroidSizer1.0 window with displayed Results Table. D)&#160;A screen shot of the format output file exported from SpheroidSizer. E)&#160;A screen shot of the list output file exported from SpheroidSizer. F)&#160;Growth of 3D tumor spheroids upon the treatments with hsp90 inhibitor and cladribine. G)&#160;Growth of 3D tumor spheroids upon the treatments with hsp90 inhibitor and adriamycin. <a href="https://www.jove.com/files/ftp_upload/51639/51639fig5highres.jpg" target="_blank">Please click here to view a larger version of this figure.</a>

Watch this Scientific Journal Video about High-throughput Image Analysis of Tumor Spheroids: A User-friendly Software Application to Measure the Size of Spheroids Automatically and Accurately at JoVE.

High-throughput Image Analysis of Tumor Spheroids: A User-friendly Software Application to Measure the Size of Spheroids Automatically and Accurately

We present a high-throughput image analysis software application to measure the size of three-dimensional tumor spheroids imaged with bright-field microscopy. This application provides a fast and effective way to examine the effects of therapeutic drugs on spheroids, which is beneficial for researchers who wish to use spheroids in drug screens.

The increasing number of applications of three-dimensional (3D) tumor spheroids as an in vitro model for drug discovery requires their adaptation to ...

Research

JoVE Journal

Biology

Tomato Analyzer: A Useful Software Application to Collect Accurate and Detailed Morphological and Colorimetric Data from Two-dimensional Objects

Measuring fruit morphology and color traits of vegetable and fruit crops in an objective and reproducible way is important for detailed phenotypic ...

High Throughput Microinjections of Sea Urchin Zygotes

Microinjection into cells and embryos is a common technique that is used to study a wide range of biological processes. In this method a small amount of ...

A High Throughput MHC II Binding Assay for Quantitative Analysis of Peptide Epitopes

Biochemical assays with recombinant human MHC II molecules can provide rapid, quantitative insights into immunogenic epitope identification, deletion, or ...

High-throughput Screening for Chemical Modulators of Post-transcriptionally Regulated Genes

Both transcriptional and post-transcriptional regulation have a profound impact on genes expression. However, commonly adopted cell-based screening assays ...

qPCRTag Analysis - A High Throughput, Real Time PCR Assay for Sc2.0 Genotyping

The Synthetic Yeast Genome Project (Sc2.0) aims to build 16 designer yeast chromosomes and combine them into a single yeast cell. To date one synthetic ...

High-Throughput, Multi-Image Cryohistology of Mineralized Tissues

There is an increasing need for efficient phenotyping and histopathology of a variety of tissues. This phenotyping need is evident with the ambitious ...

Applications of pHluorin for Quantitative, Kinetic and High-throughput Analysis of Endocytosis in Budding Yeast

Green fluorescent protein (GFP) and its variants are widely used tools for studying protein localization and dynamics of events such as cytoskeletal ...

Precise, High-throughput Analysis of Bacterial Growth

Bacterial growth is a central concept in the development of modern microbial physiology, as well as in the investigation of cellular dynamics at the ...

Diagonal Method to Measure Synergy Among Any Number of Drugs

A synergistic drug combination has a higher efficacy compared to the effects of individual drugs. Checkerboard assays, where drugs are combined in many ...

A Sectioning, Coring, and Image Processing Guide for High-Throughput Cortical Bone Sample Procurement and Analysis for Synchrotron Micro-CT

Bone is a dynamic and mechanically active tissue that changes in structure over the human lifespan. The products of the bone remodeling process have been ...