The authors thank Benjamin Krause for the fruitful discussion and his support. The authors further thank Robert G&#252;nther for his assistance with the live cell microscopy.
The work was supported by the DFG financial support NI1167/3-1 (JIMI) to R.N. and Z.C.,&#160;DFG financial support CRC 1278 PolyTarget Project Z01 for Z.C., C01 in TRR130 to R.N. and SFB633, TRR130, Exc257 to A.E.H. and J.B.S. The BfR provided intramural support SFP1322-642 for F.L.K and A.L.

1. Protocol Requirements
<ol>
	<li>Obtain high-resolution deconvolved three-dimensional (3D) microscopy data deconvolved in compliance with the Nyquist criterion with a sampling interval at least twice the highest spatial frequency of the specimen to obtain a high resolution image.</li>
	<li>Use 3D rendering software for the surface reconstruction and export.</li>
	<li>Use 3D animation software capable of running Python scripts (the Python script can be downloaded from the github repository: https://github.com/zcseresn/ShapeAnalysis) to create 2D projections.</li>
	<li>Use Fiji13 to analyze 2D projections and extract the DFT components.
	<ol>
		<li>Use the current Fiji distribution. If there already exists an installed version of Fiji, make sure that the installed version is the latest. This can be easily achieved by running the Help |Update option.</li>
		<li>Use the Active Contour plugin15, which can be downloaded from http://imagejdocu.tudor.lu/doku.php?id=plugin:segmentation:active_contour:start and should be copied into the plugins folder.</li>
		<li>Download the SHADE Fiji plugin from the github repository and copy into the plugins folder.</li>
	</ol>
	</li>
	<li>Use computational mathematics software capable of calculating Self-Organizing Maps.</li>
</ol>
2. Reconstruct the 3D Image.
NOTE: For testing purposes, an example dataset is provided in the github repository (see above).
<ol>
	<li>Start the 3D reconstruction software and open the 3D image data.</li>
	<li>Create a 3D Surface of (all) object(s).
	<ol>
		<li>Select the 3D view option and click on Surfaces. Click on the Next button (blue circle with a white triangle) to proceed with the surface creation wizard.</li>
		<li>Select the image channel for the surface reconstruction.</li>
		<li>Apply a smoothing function to avoid porous surfaces.
		<ol>
			<li>Choose a smoothing value that does not hide the details of the surface but avoids porous surfaces.</li>
		</ol>
		</li>
		<li>Select a thresholding method to find the surfaces.
		<ol>
			<li>Use an absolute intensity threshold when the objects are well-separated from the background and have an approximately uniform brightness level.</li>
			<li>Apply a local contrast threshold when the objects vary in their intensity but can still be separated from the local background and from the other objects surrounding them. Set the local threshold search area according to the value of the expected diameter of the reconstructed objects.</li>
		</ol>
		</li>
		<li>Filter the reconstructed surfaces according to morphological parameters of interest, e.g., volume, sphericity, surface-to-volume ratio, etc., and finish the surface reconstruction.</li>
	</ol>
	</li>
	<li>Save and export the generated surfaces in a format that is compatible with the 3D animation software that will be used in the next step.</li>
</ol>
3. Transform the 3D reconstructed surfaces into 2D projections
<ol>
	<li>Start Blender and go to the output tab in the right-side window. Select the TIFF format from the dropdown menu and set the color depth to 8 bit RGBA.</li>
	<li>Switch into Scripting Mode and open the provided script file &#34;GUI_AutoRotate.py&#34; from the repository provided with this work (https://github.com/zcseresn/ShapeAnalysis).</li>
	<li>Click on Run Script. Choose the folder of the wrl files when prompted for input.</li>
	<li>If needed, create more rotations when working with more complex surfaces: go to the GUI and set the box Rotations to a value above 6. 
	NOTE: A rotation of 6 different angles can be sufficient to distinguish the different cell populations. It is not recommended to create less than six rotations per surface, because of potential information loss.</li>
	<li>Run the script by clicking on the Rotate button in the GUI. Save the projections of the individual surfaces in the same folder that was used as the input folder (step 2.3). By default, the images are saved in an 8 bit Tiff format (see step 2.1), which is the format required by the Fiji plugin SHADE.</li>
</ol>
4. Find the periphery and calculate the Fourier components using Fiji.
<ol>
	<li>Open Fiji and select SHADE in the Plugins menu. Start with the default values and fine-tune the parameters later on. Click OK when ready to run the program.
	<ol>
		<li>Choose a Gradient Threshold value for the thresholding of the input image.</li>
		<li>Choose the Number of Iterations. The higher the Number of Iterations value, the more precise the reconstruction of the periphery. For simpler shapes, a lower number is usually sufficient.</li>
		<li>Use the Number of Dilations parameter to determine how much larger the starting mask is compared to the actual cell. Usually more complex shapes need more dilation steps for proper periphery finding.</li>
		<li>Check the Dark Background checkbox if the projected shapes are brighter than the background.</li>
		<li>Activate the Show Intermediate Results checkbox only when using a small test dataset to determine the performance of SHADE. Activating this option for larger datasets lowers the computational efficiency and could possibly halt a system with low video memory.</li>
		<li>Check the Save Result Tables checkbox to use the results of SHADE as an input for step 5. If the box is checked, all results are saved in individual csv files. A summary of the output data is always generated in a file called &#34;Result_collection_of_all_DFT_calculations.csv&#34;.</li>
	</ol>
	</li>
	<li>Select the input data folder that contains the TIFF files that were created in step 3.</li>
	<li>Provide the output data folder.</li>
	<li>Click OK to start the plugin.</li>
</ol>
5. Self-Organizing Maps
NOTE: SOM networks are only able to classify data when they are trained on a large dataset which contains input from all expected cell types and conditions. For demonstration purposes, such a dataset is provided and can be found in our repository (&#8220;AllCells_summary_normalised.csv&#8221; from https://github.com/zcseresn/ShapeAnalysis
<ol>
	<li>Follow these guidelines if there is no trained SOM available yet for the input data; otherwise proceed to step 5.2.
	<ol>
		<li>Start a computational mathematical software capable of performing neural network classifications.</li>
		<li>Select a data file to be used for training the SOM network. This dataset should contain all experimental conditions in order to train the SOM on the particular cell types and experimental conditions. 
		NOTE: It is also possible to use the provided AllCells_summary_normalised.csv for testing the system.</li>
		<li>Start the training and wait till the training is completed before proceeding. By default, the script is set to run 2000 iterations (&#34;Epochs&#34;). 
		NOTE: The number of iterations depends on the learning rate of the SOM. Depending on the input data it is advisable to test both higher and lower number of epochs and observe the stability of the pattern of the SOM. When using the script provided, the number of iterations can be changed under line 32. The network size can be changed in line 34 (by default it is set to 12 by 12).</li>
		<li>After the training is finished, examine the network&#8217;s topology (neighbor distances, input planes, sample hits, etc.). The network is now trained and can be saved for future use.</li>
	</ol>
	</li>
	<li>Load in the SOM, when using an already trained map (this can come either from Step 5.1 or from other sources) in order to cluster a dataset.
	<ol>
		<li>Import the csv file that is to be tested with the preloaded trained SOM. Select the csv output of the SHADE Plugin from step 4 when using data prepared by the SHADE plugin. 
		NOTE: It is also possible to use the example data files &#8220;InteractingCells_summary_normalised.csv&#8221;, &#8220;MobileCells_summary_normalised.csv&#8221; or &#8220;PhagocytosingCells_summary_normalised.csv&#8221; that are provided via github.</li>
		<li>After the classification is finished, evaluate the results of the SOM as in step 5.1.5.
		<ol>
			<li>Examine the hitmap generated from the csv file. Each cell of the map shows how many times the dataset &#34;hits&#34; that particular cell of the trained SOM. When a group of cells are clustered in a small area of this map, this indicates that the dataset is fairly homogenous. Multiple clusters will indicate that subgroups likely exist in the dataset.</li>
			<li>Examine the neighborhood weight distances. Areas of this map that are well separated correspond to groups of objects that behave very differently from the SOM&#39;s point of view. With DFT components as input data, this means that these cell groups have very dissimilar shapes of the corresponding 3D surfaces.</li>
			<li>Examine the weight planes for information about the contribution by each element of the feature vector. In case of using the 20 DFT components as described earlier, 19 maps will appear here. When using the provided example dataset, the first 5 or 6 weight planes will be different, but the rest of them will appear fairly similar. In this case it can be concluded that it would be enough to use approximately 7 DFT components.</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>

<ol>
	<li>Masedunskas, 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=Intravital+microscopy:+a+practical+guide+on+imaging+intracellular+structures+in+live+animals.">Intravital microscopy: a practical guide on imaging intracellular structures in live animals.</a> Bioarchitecture. 2 (5), 143-157 (2012).</li><li>Niesner, R. A., Hauser, A. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Recent+advances+in+dynamic+intravital+multi-photon+microscopy.">Recent advances in dynamic intravital multi-photon microscopy.</a> Cytometry A. 79 (10), 789-798 (2011).</li><li>Ho, S. Y., 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=NeurphologyJ:+an+automatic+neuronal+morphology+quantification+method+and+its+application+in+pharmacological+discovery.">NeurphologyJ: an automatic neuronal morphology quantification method and its application in pharmacological discovery.</a> BMC Bioinformatics. 12, 230 (2011).</li><li>Yin, Z., 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+screen+for+morphological+complexity+identifies+regulators+of+switch-like+transitions+between+discrete+cell+shapes.">A screen for morphological complexity identifies regulators of switch-like transitions between discrete cell shapes.</a> Nature Cell Biology. 15 (7), 860 (2013).</li><li>Yu, H. Y., Lim, K. P., Xiong, S. J., Tan, L. P., Shim, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Functional+Morphometric+Analysis+in+Cellular+Behaviors:+Shape+and+Size+Matter.">Functional Morphometric Analysis in Cellular Behaviors: Shape and Size Matter.</a> Advanced Healthcare Materials. 2 (9), (2013).</li><li>Johnson, G. R., Buck, T. E., Sullivan, D. P., Rohde, G. K., Murphy, R. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Joint+modeling+of+cell+and+nuclear+shape+variation.">Joint modeling of cell and nuclear shape variation.</a> Molecular Biology of the Cell. 26 (22), 4046-4056 (2015).</li><li>Wang, S. -. H., Cheng, H., Phillips, P., Zhang, Y. -. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Multiple+Sclerosis+Identification+Based+on+Fractional+Fourier+Entropy+and+a+Modified+Jaya+Algorithm.">Multiple Sclerosis Identification Based on Fractional Fourier Entropy and a Modified Jaya Algorithm.</a> Entropy. 20 (4), 254 (2018).</li><li>Kriegel, F. 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=Cell+shape+characterization+and+classification+with+discrete+Fourier+transforms+and+self-organizing+maps.">Cell shape characterization and classification with discrete Fourier transforms and self-organizing maps.</a> Cytometry Part A. 93 (3), 323-333 (2017).</li><li>Styner, M., 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=Framework+for+the+Statistical+Shape+Analysis+of+Brain+Structures+using+SPHARM-PDM.">Framework for the Statistical Shape Analysis of Brain Structures using SPHARM-PDM.</a> Insight Journal. (1071), 242-250 (2006).</li><li>El-Baz, 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=3D+shape+analysis+for+early+diagnosis+of+malignant+lung+nodules.">3D shape analysis for early diagnosis of malignant lung nodules.</a> Medical Image Computing and Computer Assisted Intervention. 14 (Pt 3), 175-182 (2011).</li><li>Williams, E. L., El-Baz, A., Nitzken, M., Switala, A. E., Casanova, M. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Spherical+harmonic+analysis+of+cortical+complexity+in+autism+and+dyslexia.">Spherical harmonic analysis of cortical complexity in autism and dyslexia.</a> Translational Neuroscience. 3 (1), 36-40 (2012).</li><li>Kruggel, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Robust+parametrization+of+brain+surface+meshes.">Robust parametrization of brain surface meshes.</a> Medical Image Analysis. 12 (3), 291-299 (2008).</li><li>Schindelin, J., 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=Fiji:+an+open-source+platform+for+biological-image+analysis.">Fiji: an open-source platform for biological-image analysis.</a> Nature Methods. 9 (7), 676-682 (2012).</li><li>Kohonen, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Essentials+of+the+self-organizing+map.">Essentials of the self-organizing map.</a> Neural Networks. 37, 52-65 (2013).</li><li>Andrey, P., Boudier, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Adaptive+Active+Contours.">Adaptive Active Contours.</a> ImageJ user and developer conference. , (2006).</li><li>Burger, W., Burge, M. 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> Principles of Digital Image Processing. , (2013).</li><li>Lestrel, P. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Fourier Descriptors and their Applications in Biology. , (2008).</li><li>Bayerl, S. H., 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=Time+lapse+in+vivo+microscopy+reveals+distinct+dynamics+of+microglia-tumor+environment+interactions-a+new+role+for+the+tumor+perivascular+space+as+highway+for+trafficking+microglia.">Time lapse in vivo microscopy reveals distinct dynamics of microglia-tumor environment interactions-a new role for the tumor perivascular space as highway for trafficking microglia.</a> Glia. 64 (7), 1210-1226 (2016).</li></ol>

The authors declare that they have no competing financial interests.

The identification of potentially pathological conditions using small, intact tissue samples is of high importance. Such techniques will assure a timely response to infectious diseases and aggressive types of cancer. The kinetic and morphological responses of various immune cells, e.g., microglia and macrophages, are characteristic of the immune response of the body. Although in most cases it is not practical or even possible to monitor the kinetic behavior of these cells, it is fairly straightforward to acquire 3-dimensional images to retrieve their shape. Typically, immune cells assume a complex shape in healthy tissue and a much simpler form under inflamed or cancerous conditions18. Whilst the time-dependent characteristics of such shape change would add to our understanding of the development of the immune response, using only the 3D shape of a representative group of cells can also be sufficient to determine the healthy or pathological nature of the tissue.
Characterizing the 3-dimensional surface of a cell is not a simple task. The application of spherical harmonics is a way to represent a 3D surface with a relatively large number (50-70) of components11,12. In addition, determining the spherical harmonics is computationally expensive; projecting very complex shapes onto the unit sphere is either impossible or very difficult due to the need to apply multiple grids of various fineness on the unit sphere; finally, the meaningful interpretation of the spectra of the spherical harmonics components is far from being trivial.
In our work presented here, we replace the difficult task of direct 3D surface analysis with the much simpler approach of using 2D projections of the original surface to gain sufficient morphological information to identify pathologic conditions. We demonstrated every step of this workflow by using 3D microscopy data from myeloid cells, whilst clearly pointing out that all steps were simple to complete, and the resulting 2-dimensional maps were easy to interpret.
Naturally, a 3D-to-2D projection will lead to information loss about the structure of the surface. In our example dataset of microglia in a mouse cortical tumor model, it was enough to use six angles when creating the 2D projections. However, more complex shapes, or less prominent morphological changes may require that a larger number of projections are created so as to be able to reliably identify cell subgroups with the SOM. For this reason, our approach is designed to be able to generate and analyze any number of projections. Simply by choosing a higher number of projections for more complex shapes, it is possible to scale the information loss to a tolerable minimum. As an example, the interacting cell type in Figure 4a and 4b would require a larger number of projections in order to represent the complex surface properly.
As any approximate method, the hereby proposed workflow had to be tested against the results of a manual classification process of microglia18. The results presented earlier confirmed the reliability of the automated workflow. Furthermore, the workflow is more time efficient compared to conventional analysis. The medical expert who classified the microglia cells manually needed approximately 4 weeks for his analysis of the dataset, whereas our workflow needed only about 1 day. The robustness of our approach was also clearly proven by the reproducibility of the trained SOM to a subset of data that belonged to the same cell type but was not used to train the SOM, as show in Figure 3c.
Even though our approach did not consider kinetic information, we examined the effect of timing on the DFT-based shape analysis. The most typical example for time-dependent behavior was found amongst the mobile cell population, where the contribution from the higher indexed DFT components was clearly observable, as in Figure 4a. This calls attention to the importance of utilizing a high enough number of DFT components when dealing with cell types that are likely to behave in a very time-dependent manner. Due to the automated nature and high execution speed of our software tools, the increased number of DFT components and projections will increase the precision and reliability of the results, whilst they will not appreciably hinder the computational performance.

Timely, simple and precise determination of the pathological status of biological tissue is of the highest interest in biomedical research. Mouse models provide the means to study a range of pathological conditions, such as immune reactions or cancer development, in combination with complex 3D and 4D (3 spatial dimensions and time) microscopy techniques. Microscopy studies can be performed via intravital or excised-tissue 2-photon microscopy, light-sheet microscopy, and -to a limited tissue depth of approximately 100 &#181;m- by confocal microscopy. In order to have time-related information about the cells&#39; behavior under physiological or pathological conditions, it is necessary to monitor the tissue for an extended period of time, which usually requires intravital imaging1,2. Naturally, the applicability of this technique is limited to animal models due to its invasiveness. Non-invasive techniques are also available for human applications, including a variety of tomography methods (MSOT, CT, etc.), but these methods all lack the necessary spatial -and often temporal- resolution to study behavior at the cellular level.
Static information regarding the appearance of cells may be accessible more easily via various 3D imaging techniques executed on excised tissue samples. Here, the kinetic behavior of the cells is not measured, thus it is necessary to adopt novel analysis techniques that are able to determine the pathogenic status of the examined cells based solely on their morphology3. Such an approach was used to link cell shapes and tissue textures to pathological behavior4,5,6.
In the new technique described here, the cells are reconstructed as 3D surfaces and their shapes are characterized via 3D-to-2D projections and successive Fourier-based periphery-shape analysis7,8. By reducing the dimensions from 3 to 2, the problem is simplified. It is also possible to characterize the cell surfaces in 3D by applying spherical harmonics analysis, as it has been done for medical images9. However, spherical harmonics do not handle sharp and rugged shapes well, requiring a multi-scale grid to be established on the unit sphere. In addition, the number of necessary spherical harmonics components can be large (50-70), with the underlying calculations very demanding and the results hard to interpret10,11,12.
With our newly proposed method, the task is reduced to a series of 2D shape descriptions, where the number of the 2D projections is up to the analyst and can be adjusted according to the complexity of the 3D shape. The projections are generated automatically via a Python script that runs inside a 3D animation tool. The 2D projections are described by the discrete Fourier transform (DFT) components of their periphery, calculated by a Fiji13 plugin that is provided here as part of our software package. The DFT is applied here in order to decompose the complex outline of the cell into a series of sin and cos functions. In this way, we can describe the outline with a relatively small number of DFT components, thus reducing the complexity of the problem (for further details see Equations section). The DFT components are put into a trained Self-Organizing Map (SOM14), where the existence of shape clusters can be objectively tested8. SOMs provide a competitive and unsupervised learning tool from the field of artificial intelligence. They consist of a linked array of artificial neurons which communicate with each other via a weighted neighborhood distance function. The neuronal system responds to the first element of the input dataset and the neurons whose response is the strongest are &#34;grouped&#34; nearer to each other. As the neural system receives more and more input, data neurons that repeatedly respond strongly start to form well defined cluster within the system. After proper training on a large dataset that contains 2D shape information in form of a set of DFT components, any individual cell&#39;s DFT components can be put into the trained SOM and reveal whether the cell likely belongs to the healthy or the pathogenic cell group. We expect such tool to become a great addition to the methods of scientific and clinical diagnostics.

<table>
	<tbody>
		<tr>
			<td>Imaris 9.1.2, software</td>
			<td>Bitplane, Z&#252;rich, Switzerland</td>
			<td>v.9.1.2</td>
			<td>3D image reconstruction and surface generation; this was used by us!</td>
		</tr>
		<tr>
			<td>Blender 2.75a, software</td>
			<td>https://www.blender.org/</td>
			<td>v.2.75a</td>
			<td>3D and 4D open source animation software; 2.75a is the required version for this Python</td>
		</tr>
		<tr>
			<td>Fiji /ImageJ, software</td>
			<td>https://fiji.sc/</td>
			<td>ImageJ v.1.52b</td>
			<td>Open source multi-D image analysis toolkit</td>
		</tr>
		<tr>
			<td>MATLAB</td>
			<td>MathWorks, www.mathworks.com</td>
			<td>R2017b</td>
			<td>General computational mathematical software</td>
		</tr>
		<tr>
			<td>MATLAB Machine Learning kit</td>
			<td>MathWorks, www.mathworks.com</td>
			<td>R2017b</td>
			<td>Can only be used together with MATLAB</td>
		</tr>
		<tr>
			<td>Fiji plugins: SHADE</td>
			<td>https://github.com/zcseresn/ShapeAnalysis</td>
			<td>v.1.0</td>
			<td></td>
		</tr>
		<tr>
			<td>Fiji plugins: ActiveContour</td>
			<td>http://imagejdocu.tudor.lu/doku.php?id=plugin:segmentation:active_contour:start</td>
			<td>absnake2</td>
			<td></td>
		</tr>
		<tr>
			<td>Computer</td>
			<td>Any</td>
			<td>NA</td>
			<td>See Imaris instructions for minimum computer requirements</td>
		</tr>
	</tbody>
</table>

morphology-based distinction between healthy and pathological cells utilizing fourier transforms and self-organizing maps

The appearance and the movements of immune cells are driven by their environment. As a reaction to a pathogen invasion, the immune cells are recruited to the site of inflammation and are activated to prevent a further spreading of the invasion. This is also reflected by changes in the behavior and the morphological appearance of the immune cells. In cancerous tissue, similar morphokinetic changes have been observed in the behavior of microglial cells: intra-tumoral microglia have less complex 3-dimensional shapes, having less-branched cellular processes, and move more rapidly than those in healthy tissue. The examination of such morphokinetic properties requires complex 3D microscopy techniques, which can be extremely challenging when executed longitudinally. Therefore, the recording of a static 3D shape of a cell is much simpler, because this does not require intravital measurements and can be performed on excised tissue as well. However, it is essential to possess analysis tools that allow the fast and precise description of the 3D shapes and allows the diagnostic classification of healthy and pathogenic tissue samples based solely on static, shape-related information. Here, we present a toolkit that analyzes the discrete Fourier components of the outline of a set of 2D projections of the 3D cell surfaces via Self-Organizing Maps. The application of artificial intelligence methods allows our framework to learn about various cell shapes as it is applied to more and more tissue samples, whilst the workflow remains simple.

We applied a DFT to calculate the main components of the shape corresponding to the cell projections. The Fourier descriptors were obtained by applying the DFT algorithm to the xy-coordinate pairs of the fitted periphery of the cell projections, obtained as the output of the AbSnake part of our workflow. These xycoordinate pairs can be handled as a complex-valued 2D vector &#34;g&#34;: 
<img alt="Equation 1" src="/files/ftp_upload/58543/58543eq1.jpg" />
From the vector &#34;g&#34;, we use DFT to calculate the complex-valued Fourier spectrum: 
<img alt="Equation 2" src="/files/ftp_upload/58543/58543eq2.jpg" /> 
Based on well-known formulas of the discrete Fourier spectrum, and using the complex-number labeling of &#34;g&#34; as: 
<img alt="Equation 3" src="/files/ftp_upload/58543/58543eq3.jpg" /> 
We get: 
<img alt="Equation 4" src="/files/ftp_upload/58543/58543eq4.jpg" /> (1) 
We can calculate the real (&#34;A&#34;) and imaginary (&#34;B&#34;) components of <img alt="Equation 5" src="/files/ftp_upload/58543/58543eq5.jpg" />: 
<img alt="Equation 6" src="/files/ftp_upload/58543/58543eq6.jpg" /> (2) 
<img alt="Equation 7" src="/files/ftp_upload/58543/58543eq7.jpg" /> (3) 
Here, the first DFT component G0 corresponds to m = 0, which gives: 
<img alt="Equation 8" src="/files/ftp_upload/58543/58543eq8.jpg" /> (4) 
<img alt="Equation 9" src="/files/ftp_upload/58543/58543eq9.jpg" /> (5) 
Consequently, this component describes the geometrical center of the original object. 
The second element of the DFT forward spectrum, G1, corresponds to m = 1: 
<img alt="Equation 10" src="/files/ftp_upload/58543/58543eq10.jpg" /> 
<img alt="Equation 11" src="/files/ftp_upload/58543/58543eq11.jpg" /> (6)
From Eq.6 we conclude that these points form a circle with a radius of <img alt="Equation 12" src="/files/ftp_upload/58543/58543eq12.jpg" /> and starting angle <img alt="Equation 13" src="/files/ftp_upload/58543/58543eq13.jpg" />, where the circle describes one complete revolution while the shape is traced once. The center of the circle is located at the origin (0, 0), the radius is |G1| and the starting point is:
<img alt="Equation 14" src="/files/ftp_upload/58543/58543eq14.jpg" /> <img alt="Equation 15" src="/files/ftp_upload/58543/58543eq15.jpg" /> (7)
In general, for a single Fourier coefficient <img alt="Equation 16" src="/files/ftp_upload/58543/58543eq16.jpg" />, the coordinates are described as:
<img alt="Equation 17" src="/files/ftp_upload/58543/58543eq17.jpg" /> 
<img alt="Equation 18" src="/files/ftp_upload/58543/58543eq18.jpg" /> (8)
Similarly to Eq.6, Eq.8 also describes a circle, but with a radius of Rm=|Gm|, a starting angle <img alt="Equation 19" src="/files/ftp_upload/58543/58543eq19.jpg" /> and a starting point at <img alt="Equation 20" src="/files/ftp_upload/58543/58543eq20.jpg" />, where the contour is traced once whilst the circle runs through &#34;m&#34; full orbits16,17.
Shape parameters as SOM input 
The workflow, as described in Figure 1, was applied to a deconvolved (using a measured Point Spread Function) intravital multi-photon microscopy dataset of microglial cells to characterize their morphological changes in healthy or cancerous cortical tissue18. Twenty DFT components were calculated for each 2D projection of the reconstructed 3D surfaces and the results were used as an input for the SOM training. Under physiological conditions, the microglia presented a rather complex shape with multiple, highly branched processes (Figure 2a). When placed in a cancerous environment (cortical tumor model), the microglia changed to a simpler, more spindle-like shape (Figure 2b).
The trained SOM was tested in order to evaluate its ability to distinguish between healthy and cancerous cells. The healthy cell population was projected onto a single area of the SOM (Figure 2c). The SOM responded to the cancerous microglia dataset with a dumbbell-shaped active region (Figure 2d). A blindly mixed input data set that consisted of DFT shape components from both the healthy and the cancerous group was projected by the SOM into two distinct groups, whilst keeping the shape of their individual contours similar to those of the separated groups (Figure 2e; compare with 2c and 2d). It can be concluded that the mixed dataset was successfully clustered by the SOM.
We tested the performance of the SOM by comparing its projections with the manual analysis of the same data by a medical expert, who classified the dataset based on their spatio-temporal behavior. The expert identified four distinct cell groups (resting cells, phagocytosing cells, interacting cells, and mobile cells18), which were reconstructed and used to train a 12x12 SOM. The trained network (Figure 3a) shows groups of high hit-value artificial neurons, especially in the bottom left and the middle areas of the SOM. The response of the trained network was also tested with four randomly selected subsets (which were not part of the training dataset) of images from the four different groups identified by the expert18. These image subsets resulted in four well-defined responses by the SOM, as shown in Figure 3b. The resting cells exhibit the most complex shape and showed the highest separation level within the neural network (Figure 3b &#34;resting&#34; panel). The other three identified cell types shared a common area of the SOM in the bottom left corner, but were otherwise separated by the SOM. The bottom left corner SOM area thus corresponds to the lower-index DFT values.
The robustness of the SOM approach was tested by using the trained SOM with three random subsets of the same -resting- cell type (not part of the training dataset). The response of the SOM to this input exhibits a very similar response (Figure 3c, subsets 1-3), demonstrating the robustness of our approach.
Time-dependent cell shape changes are precisely characterized by DFT 
In order to examine the effect of time-dependent changes of the cell shape on the DFT components, one to three cells per subgroup (see Figure 3b) were tracked for 13 to 28 time points. Figure 4 shows the first ten DFT components of a mobile cell (Figure 4a) and an interacting cell (Figure 4b), which were plotted as a function of time. The mobile cell exhibits a permanently altering shape (see Supplementary Video 4 in 8), which is reflected by a rougher DFT surface. The bursts of the DFT amplitude in the first third of the time course for the interacting cell coincide with the fast and vast cell shape changes as shown in Supplementary Video 5 in 8.
The time course of all 19 DFT components was also characterized for these two cells at three separate time points during the tracking of a mobile cell (Figure 5a) and of an interacting cell (Figure 5b). The perpendicular axes represent hereby the six rotation angles and indicate that all projections are equally important for the characterization of the shape for both cell types.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/58543/58543fig1.jpg" /> 
Figure 1. Step by step workflow of the data processing to identify cell clustering based on the shape of cells. Surfaces reconstructed in 3D were used as input to Blender for automated 3D-to-2D projections. The periphery of each projection was located and the DFT components were calculated. The components served as an input either to a trained SOM in Matlab, or to train a new SOM. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig1large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/58543/58543fig2.jpg" /> 
Figure 2. Typical appearance of mouse cortical microglia cells under control conditions (a) and in cancerous tissue (b) Screenshots of reconstructed microglia surfaces. SOM projections were created from the three groups of microglia samples from the mouse cortex: control (non-tumorous) cells (c), tumor cells (d), and a mixed population of cells (e). This figure has been modified with permission8. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/58543/58543fig3.jpg" /> 
Figure 3. (a, left) Self-Organizing Map of a mouse microglia dataset consisting of 768 input feature vectors. The dataset was used to train a 12x12 artificial neural network, using hexagonal neighborhood geometry, random initialization and 2000 epochs. (a, right) The corresponding SOM input planes of the first 10 DFT components (b) The responses of the SOM depicted in (a), to one random VRML file subset each from the four cell types &#34;mobile,&#34; &#34;interacting,&#34; &#34;resting,&#34; and &#34;phagocytic&#34; as first described in Figure 5 of Bayerl et al.18. (c) The response of the same SOM as in (a, left) to three random subsets of the entire dataset (which were thus not part of the training dataset) of the &#34;resting cells&#34;-type 3D surfaces. The similarity amongst the three responses is notable. This figure has been modified with permission8. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig3large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/58543/58543fig4.jpg" /> 
Figure 4. (a) Time dependence of the first 10 DFT components during an intravital imaging experiment of mouse microglia. This panel shows data for a cell of the &#34;Mobile Cells&#34; type. The x-axis corresponds to time points of the experiment at 60 s time resolution, the y-axis shows the amplitude of the DFT components in arbitrary units (a.u.), whereas the z-axis corresponds to the DFT component from 1 to 10. (b) As in (a) but for a cell of the &#34;Interacting Cells&#34; type. This figure has been modified with permission8. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig4large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/58543/58543fig5.jpg" /> 
Figure 5. (a) The behavior of all 19 DFT components of a cell of the &#34;Mobile Cells&#34; type at the beginning, at the middle and at the end of the experiment. The numbers on the x-axis correspond to the DFT component ID from 1 to 19. The y-axis shows the DFT component amplitude in arbitrary units (a.u.), whilst the z -axis marks the six random rotation angles. (b) Same as in (a) but for a cell of the &#34;Interacting Cells&#34; type. This figure has been modified with permission8. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig5large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>

Watch this Scientific Journal Video about Morphology-Based Distinction Between Healthy and Pathological Cells Utilizing Fourier Transforms and Self-Organizing Maps at JoVE.com

Morphology-Based Distinction Between Healthy and Pathological Cells Utilizing Fourier Transforms and Self-Organizing Maps

Here, we provide a workflow that allows the identification of healthy and pathological cells based on their 3-dimensional shape. We describe the process of using 2D projection outlines based on the 3D surfaces to train a Self-Organizing Map that will provide objective clustering of the investigated cell populations.

The appearance and the movements of immune cells are driven by their environment. As a reaction to a pathogen invasion, the immune cells are recruited to ...

morphology-based-distinction-between-healthy-pathological-cells

Research

Education

JoVE Core

JoVE Journal

Electrical Engineering

Cancer Research

Unit 1

z-Transform

SCIENTISTS IN ACTION

7.0K Views.  German Federal Institute for Risk Assessment (BfR). Here, we provide a workflow that allows the identification of healthy and pathological cells based on their 3-dimensional shape. We describe the process of using 2D projection outlines based on the 3D surfaces to train a Self-Organizing Map that will provide objective clustering of the investigated cell populations.

Video: Morphology-Based Distinction Between Healthy and Pathological Cells Utilizing Fourier Transforms and Self-Organizing Maps

Analyzing the Communication Between Monocytes and Primary Breast Cancer Cells in an Extracellular Matrix Extract (ECME)-based Three-dimensional System

Embedded in the extracellular matrix (ECM), normal and neoplastic epithelial cells intimately communicate with hematopoietic and non-hematopoietic cells, thus greatly influencing normal tissue homeostasis and disease outcome. In breast cancer, tumor-associated macrophages (TAMs) play a critical role in disease progression, metastasis, and recurrence; therefore, understanding the mechanisms of monocyte chemoattraction to the tumor microenvironment and their interactions with tumor cells is important to control the disease. Here, we provide a detailed description of a three-dimensional (3D) co-culture system of human breast cancer (BrC) cells and human monocytes. BrC cells produced high basal levels of regulated on-activation, normal T-cell expressed and secreted (RANTES), monocyte chemoattractant protein-1 (MCP-1), and granulocyte-macrophage colony-stimulating factor (G-CSF), while in co-culture with monocytes, pro-inflammatory cytokines Interleukin (IL)-1 beta (IL-1&#946;) and IL-8 were enriched together with matrix metalloproteinases (MMP)-1, MMP-2, and MMP-10. This tumor stroma microenvironment promoted resistance to anoikis in MCF-10A 3D acini-like structures, chemoattraction of monocytes, and invasion of aggressive BrC cells. The protocols presented here provide an affordable alternative to study intra-tumor communication and are an example of the great potential that in vitro 3D cell systems provide to interrogate specific features of tumor biology related to tumor aggression.

Analyzing the Communication Between Monocytes and Primary Breast Cancer Cells in an Extracellular Matrix Extract ECME-based Three-dimensional System

Here, we describe a three-dimensional culture method to analyze the morphology of primary breast cancer cells, as well as to study their direct/indirect interactions with monocytes and the outcomes such as collagen degradation, immune cell recruitment, cell invasion, and promotion of cancer-related inflammation.

Embedded in the extracellular matrix (ECM), normal and neoplastic epithelial cells intimately communicate with hematopoietic and non-hematopoietic cells, ...

A Preclinical Mouse Model of Osteosarcoma to Define the Extracellular Vesicle-mediated Communication Between Tumor and Mesenchymal Stem Cells

Within the tumor microenvironment, resident or recruited mesenchymal stem cells (MSCs) contribute to malignant progression in multiple cancer types. Under the influence of specific environmental signals, these adult stem cells can release paracrine mediators leading to accelerated tumor growth and metastasis. Defining the crosstalk between tumor and MSCs is of primary importance to understand the mechanisms underlying cancer progression and identify novel targets for therapeutic intervention. 
Cancer cells produce high amounts of extracellular vesicles (EVs), which can profoundly affect the behavior of target cells in the tumor microenvironment or at distant sites. Tumor EVs enclose functional biomolecules, including inflammatory RNAs and (onco)proteins, that can educate stromal cells to enhance the metastatic behavior of cancer cells or to participate in the pre-metastatic niche formation. In this article, we describe the development of a preclinical cancer mouse model that enables specific evaluation of the EV-mediated crosstalk between tumor and mesenchymal stem cells. First, we describe the purification and characterization of tumor-secreted EVs and the assessment of the EV internalization by MSCs. We then make use of a multiplex bead-based immunoassay to evaluate the alteration of the MSC cytokine expression profile induced by cancer EVs. Finally, we illustrate the generation of a bioluminescent orthotopic xenograft mouse model of osteosarcoma that recapitulates the tumor-MSC interaction, and show the contribution of EV-educated MSCs to tumor growth and metastasis formation. 
Our model provides the opportunity to define how cancer EVs shape a tumor-supporting environment, and to evaluate whether blockade of the EV-mediated communication between tumor and MSCs prevents cancer progression.

Direct injection of cancer-derived extracellular vesicles (EVs) leads to reprogramming of bone marrow supporting tumor progression; however, which cells mediate this effect is unclear. Herein, we describe a step-by-step protocol to investigate EV-mediated tumor-mesenchymal stem cell (MSC) interactions in vivo, revealing a crucial role for EV-educated MSCs in metastasis.

Within the tumor microenvironment, resident or recruited mesenchymal stem cells (MSCs) contribute to malignant progression in multiple cancer types. Under ...

Simple and Rapid Method to Obtain High-quality Tumor DNA from Clinical-pathological Specimens Using Touch Imprint Cytology

It is critical to determine the mutational status in cancer before administration and treatment of specific molecular targeted drugs for cancer patients. In the clinical setting, formalin-fixed paraffin-embedded (FFPE) tissues are widely used for genetic testing. However, FFPE DNA is generally damaged and fragmented during the fixation process with formalin. Therefore, FFPE DNA is sometimes not adequate for genetic testing because of low quality and quantity of DNA. Here we present a method of touch imprint cytology (TIC) to obtain genomic DNA from cancer cells, which can be observed under a microscope. Cell morphology and cancer cell numbers can be evaluated using TIC specimens. Furthermore, the extraction of genomic DNA from TIC samples can be completed within two days. The total amount and quality of TIC DNA obtained using this method was higher than that of FFPE DNA. This rapid and simple method allows researchers to obtain high-quality DNA for genetic testing (e.g., next generation sequencing analysis, digital PCR, and quantitative real time PCR) and to shorten the turnaround time for reporting results.

Obtaining high-quality genomic DNA from tumor tissues is an essential first step for analyzing genetic alterations using next generation sequencing. In this article, we present a simple and rapid method to enrich tumor cells and obtain intact DNA from touch imprint cytology specimens.

It is critical to determine the mutational status in cancer before administration and treatment of specific molecular targeted drugs for cancer patients. ...

Utilizing High Resolution Ultrasound to Monitor Tumor Onset and Growth in Genetically Engineered Pancreatic Cancer Models

The LSL-KrasG12D/+; LSL-Trp53R172H/+; Pdx-1-Cre (KPC) mouse model represents an established and frequently used transgenic model to evaluate novel therapies in pancreatic cancer. Tumor onset is variable in the KPC model between 8 weeks and several months. Therefore, non-invasive imaging tools are required to screen for tumor onset and monitor for response to treatment. To address this issue, different approaches have emerged over the last years. High resolution ultrasound has major advantages such as non-invasiveness, fast session times and a high image resolution without radiation exposure. However, ultrasound in mice is not trivial and sufficient anatomical knowledge and practical skills are required to successfully perform high resolution ultrasound in preclinical pancreatic cancer models. With the following article, a detailed hands-on guide for abdominal ultrasound in murine models with a particular focus on endogenous pancreatic cancer models is presented. Furthermore, a summary of common mistakes and how to avoid them is provided.

This article describes the utilization of high-resolution ultrasound in genetically engineered pancreatic cancer mice. The primary aim is to provide a detailed instruction for detection and evaluation of endogenous pancreatic tumors.

The LSL-KrasG12D/+; LSL-Trp53R172H/+; Pdx-1-Cre (KPC) mouse model represents an established and frequently used transgenic model to evaluate novel ...

Utilizing 18F-FDG PET/CT Imaging and Quantitative Histology to Measure Dynamic Changes in the Glucose Metabolism in Mouse Models of Lung Cancer

A hallmark of advanced tumors is a switch to aerobic glycolysis that is readily measured by [18F]-2-fluoro-2-deoxy-D-glucose positron emission tomography (18F-FDG PET) imaging. Co-mutations in the KRAS proto-oncogene and the LKB1 tumor suppressor gene are frequent events in lung cancer that drive hypermetabolic, glycolytic tumor growth. A critical pathway regulating the growth and metabolism of these tumors is the mechanistic target of the rapamycin (mTOR) pathway, which can be effectively targeted using selective catalytic mTOR kinase inhibitors. The mTOR inhibitor MLN0128 suppresses glycolysis in mice bearing tumors with Kras and Lkb1 co-mutations, referred to as KL mice. The therapy response in KL mice is first measured by 18F-FDG PET and computed tomography (CT) imaging before and after the delivery of MLN0128. By utilizing 18F-FDG PET/CT, researchers are able to measure dynamic changes in the glucose metabolism in genetically engineered mouse models (GEMMs) of lung cancer following a therapeutic intervention with targeted therapies. This is followed by ex vivo autoradiography and a quantitative immunohistochemical (qIHC) analysis using morphometric software. The use of qIHC enables the detection and quantification of distinct changes in the biomarker profiles following treatment as well as the characterization of distinct tumor pathologies. The coupling of PET imaging to quantitative histology is an effective strategy to identify metabolic and therapeutic responses in vivo in mouse models of disease.

Utilizing 18F-FDG PET/CT Imaging and Quantitative Histology to Measure Dynamic Changes in the Glucose Metabolism in Mouse Models of Lung Cancer

In this protocol, we describe how to utilize [18F]-2-fluoro-2-deoxy-D-glucose positron emission tomography and computed tomography (18F-FDG PET/CT) imaging to measure the tumor metabolic response to the targeted therapy MLN0128 in a Kras/Lkb1 mutant mouse model of lung cancer and coupled imaging with high resolution ex vivo autoradiography and quantitative histology.

A hallmark of advanced tumors is a switch to aerobic glycolysis that is readily measured by [18F]-2-fluoro-2-deoxy-D-glucose positron emission tomography ...

Capturing Small Molecule Communication Between Tissues and Cells Using Imaging Mass Spectrometry

Imaging mass spectrometry (IMS) has routinely been applied to three types of samples: tissue sections, spheroids, and microbial colonies. These sample types have been analyzed using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) to visualize the distribution of proteins, lipids, and metabolites across the biological sample of interest. We have developed a novel sample preparation method that combines the strengths of the three previous applications to address an underexplored approach for identifying chemical communication in cancer, by seeding mammalian cell cultures into agarose in coculture with healthy tissues followed by desiccation of the sample. Mammalian tissue and cells are cocultured in close proximity allowing chemical communication via diffusion between the tissue and cells. At specific time points, the agarose-based sample is dried in the same manner as microbial colonies prepared for IMS analysis. Our method was developed to model the communication between high grade serous ovarian cancer derived from the fallopian tube as it interacts with the ovary during metastasis. Optimization of the sample preparation resulted in the identification of norepinephrine as a key chemical component in the ovarian microenvironment. This newly developed method can be applied to other biological systems that require an understanding of chemical communication between adjacent cells or tissues.

A novel method of sample preparation was developed to accommodate cell and tissue coculture to detect small molecule exchange using imaging mass spectrometry.

Imaging mass spectrometry (IMS) has routinely been applied to three types of samples: tissue sections, spheroids, and microbial colonies. These sample ...

Pathological Analysis of Lung Metastasis Following Lateral Tail-Vein Injection of Tumor Cells

Metastasis, the primary cause of morbidity and mortality for most cancer patients, can be challenging to model preclinically in mice. Few spontaneous metastasis models are available. Thus, the experimental metastasis model involving tail-vein injection of suitable cell lines is a mainstay of metastasis research. When cancer cells are injected into the lateral tail-vein, the lung is their preferred site of colonization. A potential limitation of this technique is the accurate quantification of the metastatic lung tumor burden. While some investigators count macrometastases of a pre-defined size and/or include micrometastases following sectioning of tissue, others determine the area of metastatic lesions relative to normal tissue area. Both of these quantification methods can be exceedingly difficult when the metastatic burden is high. Herein, we demonstrate an intravenous injection model of lung metastasis followed by an advanced method for quantifying metastatic tumor burden using image analysis software. This process allows for investigation of multiple end-point parameters, including average metastasis size, total number of metastases, and total metastasis area, to provide a comprehensive analysis. Furthermore, this method has been reviewed by a veterinary pathologist board-certified by the American College of Veterinary Pathologists (SEK) to ensure accuracy.

Intravenous injection of cancer cells is often used in metastasis research, but the metastatic tumor burden can be difficult to analyze. Herein, we demonstrate a tail-vein injection model of metastasis and include a novel approach to analyze the resulting metastatic lung tumor burden.

Metastasis, the primary cause of morbidity and mortality for most cancer patients, can be challenging to model preclinically in mice. Few spontaneous ...

Orthotopic Implantation of Patient-Derived Cancer Cells in Mice Recapitulates Advanced Colorectal Cancer

Over the last decade, more sophisticated preclinical colorectal cancer (CRC) models have been established using patient-derived cancer cells and 3D tumoroids. Since patient derived tumor organoids can retain the characteristics of the original tumor, these reliable preclinical models enable cancer drug screening and the study of drug resistance mechanisms. However, CRC related death in patients is mostly associated with the presence of metastatic disease. It is therefore essential to evaluate the efficacy of anti-cancer therapies in relevant in vivo models that truly recapitulate the key molecular features of human cancer metastasis. We have established an orthotopic model based on the injection of CRC patient-derived cancer cells directly into the cecum wall of mice. These tumor cells develop primary tumors in the cecum that metastasize to the liver and lungs, which is frequently observed in patients with advanced CRC. This CRC mouse model can be used to evaluate drug responses monitored by microcomputed tomography (&#181;CT), a clinically relevant small-scale imaging method that can easily identify primary tumors or metastases in patients. Here, we describe the surgical procedure and the required methodology to implant patient-derived cancer cells in the cecum wall of immunodeficient mice.

This protocol describes the orthotopic implantation of patient-derived cancer cells in the cecum wall of immunodeficient mice. The model recapitulates advanced colorectal cancer metastatic disease and allows for the evaluation of new therapeutic drugs in a clinically relevant scenario of lung and liver metastases.

Over the last decade, more sophisticated preclinical colorectal cancer (CRC) models have been established using patient-derived cancer cells and 3D ...

Risk Prediction Model of Gastric Low-Grade Intraepithelial Neoplasia

Establishment and Evaluation of a Risk Prediction Model for Pathological Escalation of Gastric Low-Grade Intraepithelial Neoplasia

The study aims to explore the risk factors for pathological escalation after endoscopic surgery for gastric low-grade intraepithelial neoplasia (LGIN) and to establish and evaluate a risk prediction model for LGIN. A total of 120 patients diagnosed with gastric LGIN by biopsy and endoscopic submucosal dissection (ESD) between November 2020 and June 2022 were retrospectively analyzed. Gender, age, Helicobacter pylori (HP) infection, lesion size, lesion location, morphology, gastric mucosal congestion, nodules status, surface ulceration and erosion, and ME-observation of all patients were collected and divided into upgraded and non-upgraded groups according to the biopsy and ESD postoperative pathological diagnosis results. Independent risk factors for pathological escalation after ESD surgical treatment were screened by logistic regression analysis, and a risk prediction model was established. Among the 120 patients with gastric LGIN, 49 patients developed postoperative pathological upgrading; the rate of pathological upgrading was 40.83%. Among them, 42 patients were upgraded to high-grade intraepithelial neoplasia (HGIN), 1 case was upgraded to advanced gastric cancer, and 6 cases were upgraded to early gastric carcinoma (EGC). Univariate analysis showed that age, lesion size, gastric mucosal congestion, surface ulcers, and erosion were significantly different between the groups (p &lt; 0.05). Multivariate Logistic regression analysis revealed that age &#8805;60 years, focal length &#8805;2 cm, gastric mucosal congestion, and surface ulceration and erosion were independent risk factors for postoperative pathological escalation in patients with gastric LGIN. Final joint probability forecasting model for P = 1/[1 + e(26.515-0.161 x &#946;1-0.357 x &#946;2+0.039 x &#946;3-0.269 x &#946;4)]. Age, lesion size &#8805;2 cm, gastric mucosal congestion, and lesion surface ulceration and erosion are risk factors for postoperative pathological upgrading in patients with gastric LGIN. The risk prediction model established in this study based on risk factors has predictive value and can provide a scientific reference for the clinical treatment of patients with gastric LGIN.

Author Spotlight: Advancing Early Detection and Treatment of Gastrointestinal Tumors

Here, we performed a systematic evaluation of patients diagnosed with gastric low-grade intraepithelial neoplasia by previous endoscopic biopsy and obtained a pathological diagnosis by complete resection of the lesion by endoscopic submucosal dissection (ESD), analyzing factors that potentially increase the risk of pathological escalation.

The study aims to explore the risk factors for pathological escalation after endoscopic surgery for gastric low-grade intraepithelial neoplasia (LGIN) and ...

Standardizing Ovarian Cancer Organoid Biobanking

Strategy for Biobanking of Ovarian Cancer Organoids: Addressing the Interpatient Heterogeneity across Histological Subtypes and Disease Stages

While the establishment of an ovarian cancer biobank from patient-derived organoids along with their clinical background information promises advances in research and patient care, standardization remains a challenge due to the heterogeneity of this lethal malignancy, combined with the inherent complexity of organoid technology. This adaptable protocol provides a systematic framework to realize the full potential of ovarian cancer organoids considering a patient-specific variability of progenitors. By implementing a structured experimental workflow to select optimal culture conditions and seeding methods, with parallel testing of direct 3D seeding versus a 2D/3D route, we obtain, in most cases, robust long-term expanding lines suitable for a broad range of downstream applications. 
Notably, the protocol has been tested and proven efficient in a great number of cases (N = 120) of highly heterogeneous starting material, including high-grade and low-grade ovarian cancer and stages of the disease with primary debulking, recurrent disease, and post-neoadjuvant surgical specimens. Within a low Wnt, high BMP exogenous signaling environment, we observed progenitors being differently susceptible to the activation of the Heregulin 1 &#223; (HER&#223;-1)-pathway, with HER&#223;-1 promoting organoid formation in some while inhibiting it in others. For a subset of the patient's samples, optimal organoid formation and long-term growth necessitate the addition of fibroblast growth factor 10 and R-Spondin 1 to the medium. 
Further, we highlight the critical steps of tissue digestion and progenitor isolation and point to examples where brief cultivation in 2D on plastic is beneficial for subsequent organoid formation in the Basement Membrane Extract type 2 matrix. Overall, optimal biobanking requires systematic testing of all main conditions in parallel to identify an adequate growth environment for individual lines. The protocol also describes the handling procedure for efficient embedding, sectioning, and staining to obtain high-resolution images of organoids, which is required for comprehensive phenotyping.

Author Spotlight: Advancing Personalized Medicine in Ovarian Cancer 

This protocol offers a systematic framework for the establishment of ovarian cancer organoids from different disease stages and addresses the challenges of patient-specific variability to increase yield and enable robust long-term expansion for subsequent applications. It includes detailed steps for tissue processing, seeding, adjusting media requirements, and immunofluorescence staining.

While the establishment of an ovarian cancer biobank from patient-derived organoids along with their clinical background information promises advances in ...