Video tracking of single cells to identify clustering behavior

Cancer cell clustering is a critical factor in metastasis, with cells often believed to migrate in groups as they establish themselves in new environments. This study presents preliminary findings from an in vitro experiment, suggesting that co-culturing cells provides an effective method for observing this phenomenon, even though the cells are grown as monolayers. We introduce a novel single-cell tracking approach based on graph theory to identify clusters in PC3 cells cultivated in both monoculture and co-culture with PC12 cells, using 66-h time-lapse recordings. The initial step consists of defining “linked” pairs of PC3 cells, laying the foundation for the application of graph theory. We propose two alternative definitions for cell pairings. The first method, Method 1, defines cells as “linked” at a given time t if they are close together within a defined time period before and after t. A second potential alternative method, Method 2, pairs cells if there is an overlap between the convex hulls of their respective tracks during this time period. Pairing cells enables the application of graph theory for subsequent analysis. This framework represents a cell as a vertex (node) and a relation between two cells as an edge. An interconnected set of high-degree nodes (nodes with many connections or edges) forms a subgraph, or backbone, that defines a patch (cluster) of cells. All nodes connected to this backbone are part of the subgraph. The backbone of high-degree nodes functions as a partition (or cut) of the initial graph. Two consecutive clusters in the video are considered to share the same identity if the following cluster contains at least p = 75 % of the cells from the preceding cluster, and the mean positions of their cells are within △r = 75μm. PC3 cells grown in co-culture appear to form persistent clusters exceeding 10 cells after 40–50 h incubation following seeding. In contrast, PC3 cells cultured alone (mono-culture) did not exhibit this behavior. This approach is experimental and requires further validation with a broader dataset.

1 Introduction

This work aims to improve access to data on clustering within cancer cell populations. A cluster is defined as a small group of cells that maintain spatial proximity and exhibit coordinated behavior over time (Gopinathan and Gov, 2019). This phenomenon is thought to play a critical role in metastasis, and new observations could enhance our understanding and lead to novel approaches to cancer treatment. Understanding the characteristics and dynamics of these clusters is essential for uncovering the mechanisms behind cancer progression, and single-cell clustering assays offer a powerful approach to studying how such cell groups form, interact, and evolve over time.

Single-cell clustering assays can reveal how cells communicate and influence each other, particularly when they cluster and respond collectively. Analyzing clustering can help identify early patterns of abnormal cell behavior, potentially leading to earlier disease diagnosis. Several authors have contributed to advancing this field (Liu et al., 2023; Maity et al., 2024; Shannon et al., 2024).

The present preliminary in vitro experiments suggest that co-culturing cells may promote increased clustering. This finding stems from a clustering analysis approach inspired by graph theory (Ioannides et al., 2022). Our method utilizes cell movement data to identify clusters, as demonstrated by analyzing 66-h time-lapse recordings of PC3 prostate cancer cells in both monoculture and co-culture with PC12 nerve cells. These data suggest that PC3 cells interact in ways that influence their movements, particularly when co-cultured with PC12 cells. Lasting clusters begin to form after ~35 h of recording when PC3 cells are co-cultured with PC12 cells.

The current focus on using monolayer (2D) cells can, in principle, be generalized to apply to analyze recordings of cells grown in 3D, which offer an environment more similar to in vivo situations. However, it is advantageous to first test the 2D approach, as it can more easily provide tracking data over several generations. Recent advances in deep learning-based cell tracking may change this situation (Mosier et al., 2021; Wen and Kimura, 2022; Merino-Casallo et al., 2022; Freckmann et al., 2022; Wiggins et al., 2023).

The rationale for this showcase of methods is the general idea that unicellular organisms can exhibit collective behavior, such as flocking (Ling et al., 2019). Bacterial biofilms, where bacteria work together to form complex structures, are well-established examples of how unicellular organisms can enhance their survival and proliferation. Cells are in general competent to produce a quorum signal (Niu and Wang, 2012). Cooperation is useful for individuals to reach a collective benefit, share information or neutralize threats (Wrenn et al., 2021).

Cancer cells are often associated with selfish behavior. However, they can cooperate and maintain physical contacts by forming clusters. This may facilitate the metastatic cascade and promote disease progression (Archetti and Pienta, 2019). Such clusters of cells appear to include distinct cellular states/phenotypes in which “leaders” and “followers” can affect the migratory pattern of clusters, as observed in melanoma and breast cancer Haeger et al. (2020); Khalil et al. (2017). Such “leaders” and “followers” need to sense the attractant through the extracellular matrix integrin signaling and adhere to each other to coordinate their movements with robust directionality (Colak-Champollion et al., 2019).

A leader-follower organization among cells is important for successful invasion and metastasis (Wrenn et al., 2021). Disrupting this organization can therefore greatly suppress the collective migration and metastatic potential (Cheung et al., 2013; Gao et al., 2017; Zhang et al., 2019; Yang et al., 2019; Khalil et al., 2020). A recent study by Gómez-de Mariscal et al. (2024) demonstrated collective cell behavior over a 14-h period, with images captured at 10-min intervals. They showed, for example, that cells closer to the leading edge exhibited more directional movement compared to those farther away. While the distance of cells to the leading edge remained constant from the beginning to the end of the tracking period, the distribution of the data was broad.

Co-culture experimental systems to identify cooperative phenotypes between “leaders” and “followers” have shown that leader cells can maintain their invasive phenotype and that the traditional EMT signature alone cannot be utilized to identify them (Konen et al., 2017). Differences in metabolism, epigenetic modifications and gene mutations might also be important for identification of cooperative behavior (Zoeller et al., 2019; Commander et al., 2020; Summerbell et al., 2020).

The outline of this paper is as follows. Sections 2.2–2.4 detail the production of test video of PC3 cells in mono-culture as well as in co-culture with PC12 cells. Sections 2.5.1–2.5.2 describes refinement of cell positional data obtained from singe-cell tracking, as well as demonstrating how to distinguish between PC3 and PC12 cells in video based on a nearly simultaneous fluorescence image at start. Section 2.5.3 provides two alternative methods to identify patches among PC3 cells based on track data from them. Section 3.1 brings a comparison between these two methods, suggesting that the choice of method may not be critical for identification of patches. Section 3.2 shows an example where PC3 cells in co-culture seem to form more stable clusters as compared to when they live in mono-culture. Section 4 discusses the limitations of the present method showcase, which is based on limited test data.

2 Materials and methods

2.1 General work flow

Our approach for cluster identification involves the following steps:

1. Cultivation of PC3 and PC12 cells (Section 2.2).

2. Capture a still image of the cells (with PC3 cells labeled with GFP) and record a video of them in monoculture and co-culture with PC12 cells (Section 2.4).

3. Identify PC3 and PC12 cells in the co-culture by comparing the still image to the video image nearest in time (see Section 2.5.2 and Figure 1).

4. Track individual cells, ensuring that PC3 cells maintain their GFP label throughout their lineage. Refine and manually correct the tracks as needed (Section 2.5.1).

5. Identify clusters based on the refined tracks (Section 2.5.3).

Figure 1

Figure 1. Manual classification of cells in a label-free video image [(right), at the start of recording] using nearly simultaneous GFP fluorescence imagery (left), where GFP labels PC3 cells with green fluorescence. Green dots in the label-free video image correspond to GFP-labeled PC3 cells, while red dots indicate classified PC12 cells. Classifying a single cell in the video image based on its GFP fluorescence (green) allows for the classification of all members of its lineage.

2.2 Cell culture

Prostate cancer PC-3 luc2-GFP cells were purchased from Caliper (# 133416) and neuronal PC-12 Adh cell lines were purchased from ATCC (# CRL-1721.1). The PC-3 cells were maintained in RPMI-1640 (# R8758; Sigma) with 10 % FBS (fetal bovine serum)(androgen-proficient medium; F7524; Sigma) and the PC-12 Adh cells were cultivated in F-12K medium ((# 21127022; Thermo Fisher Scientific) with 2.5 % FBS and 15 % horse serum (# 16050122); Thermo Fisher Scientific). Cells were maintained at 37 °C in a humidified 5 % CO₂ atmosphere. During the cell tracking experiments, both cell lines were maintained in phenol red-free RPMI1640 (# 32404014; Thermo Fisher Scientific) with 10 % FBS. Cells were routinely checked for mycoplasma contamination.

2.3 Co-culture establishment

Prostate cancer PC-3 luc2-GFP and PC-12 Adh cells were added in co-culture in a 48 wells plates from Costar (# 3548). Both cell lines were cultured at a density of 1,500 cells per well in RPMI-1640 (# R8758; Sigma) medium without phenol red and 10 % FBS. Cells were then incubated at 37 °C in a humidified 5 % CO₂ atmosphere for 48 h. Cells in mono-culture were also seeded at a density of 3,000 cells per well and incubated using the same conditions as the co-culture. All wells were visually inspected to confirm that there were no abnormalities in the wells and that cells look viable.

2.4 Time-lapse video microscopy

PC-12 Adh and PC-3 GFP cells were cultured in mono- or co-culture for time-lapse imaging. Cells were imaged with Cytation5 (BioTek) with temperature and gas control set to 37 °C and 5 % CO₂ atmosphere, respectively. Sequential imaging of each well was taken using 10 × objective and with an interval of 5 min for 66 h incubation period. Four images with. 5 % overlap from each well were used to stitches together the images used for tracking (Gen5 software; BioTek).

2.5 Video-tracking of cells

2.5.1 Cell positional data

Single-cell tracking was performed using the in-house experimental software Kobio_Celltrack¹ for convenience. Alternatively, Fiji (Schindelin et al., 2012) and TrackMate (Ershov et al., 2022) could also generate the required tracking data. Another viable option is Btrack (Ulicna et al., 2021). The in-house system allows users to define a rectangular region in the center of the video, sized to contain a specified minimum number of cells of each type at the start of recording (Korsnes and Korsnes, 2015, 2018; Quinsgaard et al., 2024). In the co-culture experiment, a frame was set to include 40 PC3 cells at the start, which automatically included 80 PC12 cells. In the monoculture experiment, the frame was defined to contain 50 PC3 cells. These cells and their descendants were tracked throughout the experiments. Videos demonstrating the performance of the tracking system are available for validation through sample tests (Korsnes et al., 2024b,c).

The initial step of the tracking provides cell positional data, which is susceptible to errors that can influence estimates of track length and cell speed. Irregularities in cell movements can inherently make track length scale-dependent or dependent on both temporal and spatial resolution (Korsnes and Korsnes, 2023). Also note that the current sampling interval of 5 min results in missing intermediate positions. Additionally, the precise definition of a cell's position may not be straightforward. While the center of the nucleus might be considered a natural reference point, it's not always discernible in various scenarios, and some cells exhibit multiple nuclei.

The present work takes a pragmatic approach to the above problematic issues and employs a Gaussian filtering of the track positions as follows (Quinsgaard et al., 2024). Assume that the positional vector r_i represents the location of a cell at time t_i. Define the (Gaussian) weighted sum for all these n observed positions of the cell:

$\begin{array}{ll}{s}_i=\frac{1}{M_i}{\displaystyle \sum _{k=1}^n}{r}_{k}exp(-\frac{{(t_i-t_k)}^2}{T^2}).& (1)\end{array}$

for i = 1, 2, …, n and where $M_i=\sum _{k=1}^{n}exp(-\frac{{(t_i-t_k)}^2}{T^2})$. This renders s_i as a smoothed (filtered) version of the time series r_i with bandwidth T. This work applies T = 15min, considering the sample period is 150 s and that cells tend to move a distance less than their diameter during 15 min. Tests show that final results seem not to be sensitive to perturbations of T.

2.5.2 Cell type identification

Video tracking enables the classification of cells into distinct cell types based on single, near-simultaneous GFP fluorescence images, typically acquired at the commencement of the recording. Once a cell's classification is determined, its entire lineage inherits the same designation. Figure 1 exemplifies this approach.

2.5.3 Automatic patch identification

This work demonstrates two types of algorithms for identifying clusters of PC3 cells that move together as a unit, even as individual cells within the cluster rearrange positions. The approaches rely on analyzing cell movement data over short time periods, limiting the temporal resolution of the methods. These methods, exemplified by Figure 2, differ in how they define which cells are considered “linked” or “related” based on their movement patterns. Graph theory then serves as the basis for the subsequent data treatment, where cells are represented by vertices (also known as nodes. When two cells are defined to be related or “linked”, this relation is called an edge. The number of edges (links) of a node x is normally written G(x).

Figure 2

Figure 2. Illustration of two methods, based on call trajectories (black curves), to identify a patch of PC3 cells ($\text{green}$ labels) among PC12 cells ($\text{red}$ labels).

• (Upper left) Subset of video image showing a patch of PC3 cells among PC12 cells at 50 h after the start of recording. The patch is identified as “id4” (cf. Figure 4 below).

• (Upper right) Trajectories of cells in the patch within during their lifetime, starting at 35 h from the beginning of the recording until its end, suggesting a stable cluster.

• (Lower left) Illustration of Method 1 for patch identification.

• The black curves here represent trajectories of cells 2.5 h before and after the image time.

• $\text{Red}$ and $\text{blue}$ lines depict the backbone and peripheral parts of the network, respectively (see Section 2.5.3).

• Backbone nodes (cells) have at least five connections (G(x) = G_min≥5), in contrast to others.

• (Lower right) Illustration of Method 2 where $\text{aquamarine}$ lines show the convex hull of the tracks of individual cells during 2.5 h before and after the image time (total 5 h). Overlapping convex hulls of two cells' tracks define a link between them, creating a network ($\text{red}$ lines) that defines a patch.

2.5.3.1 Method 1: Contact Network Formation

This method defines two PC3 cells as “linked” at time t if they come within △d = 30μm of each other within a time frame of △t = 2.5h before and after t. A patch (of PC3 cells) is then defined as a collection of nodes (cells) connected via a “backbone” of nodes x with degree G(x)≥G_min = 5. Cells directly linked to this backbone are also considered part of the patch. The backbone network serves to avoid the influence of stray cells on the identification of patches. It also prevents groups of stationary cells from forming patches, as it is unlikely for a cell to have at least five close neighbors simultaneously. Note that the use of a backbone can reduce the ability to observe small clusters of size less than five cells.

2.5.3.2 Method 2: Convex Hull Overlap

This method relies on computing the convex hulls of PC3 cell tracks within △t = 2.5h before and after t. Two cells are considered linked at time t if their respective convex hulls overlap. In this case, the illustrations below do not apply a backbone as in Method,1, and a patch is simply defined as an interconnected subset of the linked PC3 cells. This method serves only as an example to illustrate approaches alternative to Method,1. It may be of interest as an alternative or complement to Method,1 since it enables the identification of small clusters.

2.5.3.3 Identity of patches during time

The following outline applies for both methods above. Let C₁ and C₂ represent the set of cells in patches at subsequent image times t₁ and t₂, respectively. They are given the same unique identity if

$\begin{array}{ll}\frac{|C_1\cap C_2|}{\left|C_1\right|}\ge p& (2)\end{array}$

where p = 0.75 and |·| represents the size (cardinality) of a set (i.e., the number of cells). Additionally, the mean positions must be within △r = 75μm. In other words, patches retain their identity across subsequent image times if most of the cells at one time step remain in the next, and the patches do not move too much. ID conservation enables determination of duration of patches. This is important to distinguish between behaviors of cells under different treatments.

3 Results

3.1 Two methods for patch identification

Figure 2 illustrates the above two approaches (Method 1 and Method 2) to identify patches of PC3 cells in co-culture with PC12 cells.

These methods differ in that Method 2 relies more on cell movements compared to its counterpart, Method 1, which links two cells at a specific time point, t, if they are close together (within a distance △d = 30μm) for a moment within a time window △t = 2.5h before and after t. This means cells don't necessarily need to move much to be linked. However, a cell must link to at least five other cells to become part of the actual network's backbone, which prevents groups of stagnant cells from forming a patch (see Section 2.5.3). Hence, if cell density makes cells to slow down their movements, Method 1 is not likely to find patches among them.

Method 2 uses cell movements to define links between cells. Two cells are considered “linked” at time t if the convex hulls of their respective paths intersect during the same time window, as in Method 1. This means that the cells must cross each other's paths to be linked. Therefore, the two methods for patch identification can yield different results (Figure 3).

Figure 3

Figure 3. Comparison of Method 1 and Method 2 for identifying clustering in PC3 cells ($\text{green}$) in co-culture with PC12 cells ($\text{red}$). The images were captured 50 h after the start of recording. Method,2 relies more heavily on cell movement compared to Method,1, leading to differences in results, particularly for short-lasting patches. The labels provide the patch ID, duration of the patch, and the current number of cells within the patch.

3.2 Comparison of cells in mono- and co-culture

The following data treatment uses Method 1 for patch identification, leaving Method 2 as a potential alternative for further work.

Patch formation among PC3 cells seems to depend on their environment, with patches appearing to dissolve in mono-culture (Figure 4).

Figure 4

Figure 4. Development of PC3 cell patch sizes in mono-culture and co-culture with PC12 cells, derived from the present experiment data (Korsnes et al., 2024c,b). Both cultures included 12 patches, identified as id1, id2, ..., id12. The thickness of the elongated green spots represents the number of cells in a patch at a given time after recording began. The numbers to the right of each green spot tells the number of cells at the end of its observation. The blue dotted line indicates the end of the video recording. This figure includes only patches lasting longer than 5 h, which is presumably well above the temporal resolution of the current method for identifying patches (cf. 2.5.3). Notably, patches in the mono-culture appear to dissolve shortly after 10 h. In contrast, many patches in the co-culture persist to the end of recording.

Several of these short-lived patches of PC3 cells in mono-culture may be artifacts of the current identification method.

Patches observed among PC3 cells in co-culture appear to persist longer than those seen in monoculture. Notably, the patch labeled id4 in the co-culture appears around 35 h and remains throughout the entire recording (Figures 4, 5). Figure 2 shows this same patch at 50 h, along with the tracks of the cells contained within it. These tracks indicate cell-cell cohesion. If the cells were moving independently, the observed track lengths suggest that they would disperse across a much larger area. Hence, the patch cannot be merely an artifact of the patch identification method. Figure 5 further supports this conclusion by depicting the tracks of the initial cells in id4 and their descendants.

Figure 5

Figure 5. Trajectories of the initial PC3 cells in patch id4 of the co-culture example, including their descendants. $\text{Blue}$ and $\text{Green}$ represent the first and second generations from the start of the patch, respectively. There is no third generation. The patch appears at 35 h from the start of video recording and persists throughout the recording, which ends at 66 h.

It also indicates movement of the whole patch.

The observed data, visualized in Figure 6, reveals segregation among PC3 cells in mono-culture. Here, groups of cells exhibit distinct movement speeds. This segregation of cells suggests intriguing cellular interactions, deserving further investigation through concepts like flocking or synchronization.

Figure 6

Figure 6. PC3 cells in monoculture: Illustration of differences in cell velocity among groups of cells 40 h after start of recording. The $\text{blue}$ and $\text{yellow}$ curves represent the trajectories of cells within 30 min before and after the image time, respectively. Therefore, the length of the tracks indicates individual cell velocities. The cells encircled by the $\text{aquamarine}$ curve (in the lower left corner) are significantly more stagnant than those in the patches encircled by $\text{magenta}$ lines. As in Figure 3, the labels provide the patch ID, duration of the patch, and the current number of cells within the patch.

Figure 7 further supports the concept of segregation among the PC3 cells.

Figure 7

Figure 7. Probability density function (PDF) for 6-h average cell speed during the period from 10 to 55 h after the start of video recording for all tracked PC3 cells in mono-culture (left) and co-culture with PC12 cells (right). The scale bar indicates the color coding corresponding to different levels of probability density. The cells in mono-culture show a broader variation compared to the cells in co-culture.

It shows distributions of 6-h average cell speed of PC3 cells in mono-culture and co-culture. Compared to co-culture, PC3 cells in mono-culture exhibit a broader distribution of speeds, indicating a greater heterogeneity in their movement patterns.

A linear increase in the mean squared displacement (MSD) suggests diffusive motion (Qian et al., 1991; Huda et al., 2018; Reynolds, 2018). Figure 8 builds on this concept and offers a preliminary test for potential cell clustering. It illustrates the temporal evolution of the MSD of PC3 cells in both mono- and co-culture over the period from 40 to 60 h after recording began. The movement of PC3 cells in mono-culture follows a pattern characteristic of diffusive motion, whereas in co-culture, the cells appear to be more confined to clusters. Figure 9 further indicates that cohesion between spatially close PC3 cells affects their movements, leading them to form clusters.

Figure 8

Figure 8. Mean square displacement (MSD) for all complete PC3 cell trajectories between 40 and 60 h from the start of video recording. The $\text{orange}$ curve represents PC3 cells in mono-culture (35 trajectories), while the $\text{blue}$ curve represents PC3 cells co-cultured with PC12 cells (30 trajectories). Please note that these descriptive statistics are based on a single experiment, where cells in the co-culture appear to follow a particular configuration. Therefore, this data summary does not provide an average across multiple experiments.

Figure 9

Figure 9. Cell trajectories reveal dynamic interactions between PC3 cells. (Left panel) Phase contrast image of PC3 cells in mono-culture at the start of recording. Colored arrows mark four cells. (Right panel) Trajectories of the marked cells over 15 h, highlighting intercellular interactions. These observations suggest that dynamic cell-cell interactions contribute to the initiation of clustering. See supplementary video (Korsnes et al., 2024a).

4 Discussion

The current data suggest that PC3 cells in co-culture with PC12 cells tend to form lasting clusters. A simple intuition is that cohesion between PC3 cells causes them to expel PC12 cells, leading to the formation of clusters. Subsequently, they may synchronize internal processes and become a distinct and prevailing entity. This presumably depends on the size of the patches. Contact with an environment of similar (PC3) cells may lead to the dissolution of a patch, as indicated by Figure 4. This could prevent further synchronization to initiate. A treatment that affects cell cohesion could help test this hypothesis.

The tendency to form persistent clusters in co-culture provides an opportunity for more in-depth study of clustering. The approach used to identify clusters, inspired by graph clustering theory (Schaeffer, 2007), seems to offer a valuable tool for such investigations. However, developing a generic cluster identification algorithm that can be applied across a wide range of experimental conditions and cell types would require more extensive testing. For example, the current algorithm, Method 1, relies on specific values for the parameters p, △t, △d, G_min, and △r. A potential generalization could involve calculating these parameters automatically to identify the most persistent clusters below a certain size, as long-lived small clusters may be particularly relevant to many applications related to cancer spread.

The following equation can replace Equation 2 to define the conservation of patch identity over time:

$\begin{array}{ll}\frac{\left|C_1\Delta C_2\right|}{|C_1\cup C_2|}\le q& (3)\end{array}$

for a value of q well below 1, where Δ represents the symmetric set difference (i.e., |C₁ΔC₂| is the number of cells that are in either C₁ or C₂ but not in both). This ensures that identified patches do not conserve identity if they suddenly increase in size.

The current video recordings are too short to capture situations where cells slow down due to high confluence, whereas Method 1 and Method 2 depend on cell movement. If the cells in a patch or their descendants do not (slowly) disperse, then the patch may be assumed to exist. Notably, if cells within certain clusters have entered a migratory state, one might expect them to be the last to slow down due to increased confluence. This could provide an opportunity for video-based observation of such changes in states.

Theories on clustering graph ² are relevant for a wide range of applications and provide inspiration for generalizing the present approach. Note that a graph consists of vertices (nodes) that are pairwise connected by edges, with associated weights reflecting “similarities” between two vertices (cells) or their spatial distance. A possible modification of Method 1 is to apply as weight a Gaussian kernel:

$\begin{array}{ll}K(r_i,r_j)=exp(-\gamma ||r_i-r_j|{|}^2)& (4)\end{array}$

where γ = 1/(2σ²) and σ determines the width of the kernel (Von Luxburg, 2007). Here, r_i and r_j represent the spatial position of the nodes i and j.

Several measures of “similarities” between two cells can be combined to reflect membership in a cluster or partition of a graph through cutting wedges based on “cost” of cuts (Nascimento and De Carvalho, 2011). This approach provides the opportunity to include cellular phenotypes in addition to movements in the process of identifying cellular clusters. The present test data indicate that distributions of cell speed, displacement and clustering parameters differ for PC3 cells in mono- and co-culture (cf. Figures 4, 7, 8). Such distributions may therefore serve as simple general indicators of clustering tendencies and help check for consistency of results from experiments. A more comprehensive data analysis is necessary to definitively incorporate cell speed and morphology into an algorithmic method for cluster identification.

The framework of graph theory allows for the application of theories on percolation and diffusion on graphs (Stauffer and Aharony, 2018; Chung, 1997) to further refine cluster identification from basic tracking data. These theories provide measures of how well a graph or cluster (subgraph) is connected.

An intriguing speculation is that small clusters of closely interacting cells may serve as potential “sites of ignition” for collective changes in behavior among cancer cells through synchronization, a well-documented phenomenon observed in interacting organisms that influences fundamental cellular processes (Hannon and Ruth, 2014). The present identification of clustering may therefore represent only an initial step in exploring collective changes among subsets of cells. The current co-culture experimental setup has the potential to facilitate the formation of small clusters, supporting such investigations of collective behavior. Qualitative studies are a natural initial step in exploring this theory of possible sporadic “sites of ignition.”

A simple, intuitive analogy for this work is that juggling multiple challenges at once is akin to being caught in a “fire overlap” situation, where things become significantly more complex compared to handling each challenge individually. This creates opportunities for cancer cells to migrate as a group. There may be links between cell cluster formation, poor disease prognosis, and resistance to drug treatments (Au et al., 2016; Lee et al., 2017; Bithi and Vanapalli, 2017; Andrei et al., 2020). Their greater colonization efficiency may stem from factors such as protection against anchorage-dependent apoptosis, cooperation among heterogeneous cell phenotypes within the clusters, and protection from assaults by immune cells (Hou et al., 2012; Yu et al., 2013; Hong and Zu, 2013).

Knowledge of cell cluster biogenesis and organization can contribute to a better understanding of cancer dissemination. Clusters serve as useful monitoring tools and prognostic biomarkers, especially in metastasis, likely due to cell jamming or collective migration mechanisms that remain largely unknown (Amintas et al., 2020). A promising therapeutic approach suggests that inhibiting or disrupting cluster formation could reduce cancer dissemination (Au et al., 2016; Gkountela et al., 2019).

The present data analysis shows that PC3 cells can exhibit coordinated movement within stable clusters when grown in co-culture with PC12 cells. Cluster sizes from these data align with previous findings of Au et al. (2016), who identified clusters of 20 cells in malignant melanoma. However, clusters of 2–6 cells are more commonly reported in the literature (Molnar et al., 2001).

One aspect worthy of investigation is whether cancer cells may collaborate to achieve perineural invasion (PNI), a complex phenomenon characterized by reciprocal interactions between cancer cells and the surrounding nerve micro-environment (Magnon, 2015; Grigore et al., 2015; Bakst and Wong, 2016). PNI can occur in various malignant tumors (Demir et al., 2012; Schmitd et al., 2018) and it is believed to play a role in the development of neuroendocrine prostate cancer (NEPC), an aggressive variant of prostate cancer where prostate adenocarcinoma cells (PCA) trans-differentiate into a neuroendocrine (NE) cell phenotype to escape anti-androgen therapies (Braadland et al., 2015; Kaarijärvi et al., 2021). The infiltration of cancer cells into this space can actively stimulate cancer progression and metastasis (Chen et al., 2019; March et al., 2020).

Schwann nerve cells can act as “leaders”, reorganizing cancer cell clusters to facilitate PNI when in contact with them (Deborde et al., 2016). They can also collectively function as tumor-activated Schwann cell tracks (TAST), promoting cancer cell invasion and migration (Deborde et al., 2022). Physical contact among cells is required to enhance cancer invasion. Cluster detection often correlates with higher rates of disease progression and poorer treatment responses (Wrenn et al., 2021).

It is becoming evident that single-cell profiling allow researchers to address the co-occurrence of molecular events in individual cells, however, it does not need to be limited to gene expression (Kumar et al., 2017). The incorporation of visualization data on a per-cell basis will complement the knowledge gained from single-cell molecular profiling, including epigenetic modifications, which often contributes to “stochastic” expression (Angermueller et al., 2016).

There have been some limited contributions using in vitro co-cultures to study cell clustering behavior (Wrenn et al., 2021). Recent research, however, emphasizes the role of specific molecular mechanisms in promoting cluster formation, which is linked to an increased risk of metastasis. Tumor cells that cluster within the tumor microenvironment (TME) tend to gain an evolutionary advantage. The TME is a complex ecosystem comprising tumor cells, immune cells, fibroblasts, blood vessels, and extracellular matrix (ECM). Clustering facilitates cooperation among tumor cells, enhancing their ability to survive, migrate collectively, spread, and evade the immune system. The presence of cell clustering is often associated with poor prognosis and a heightened likelihood of metastasis in various cancer types (Rozenberg et al., 2023). Disrupting these mechanisms could offer a promising avenue for cancer treatment, making efforts to uncover them of significant interest.

5 Conclusions

PC3 cells cultured in monolayers in test wells, alongside PC12 cells, tend to form stable clusters. This behavior may be attributed to cohesion between PC3 cells, providing an opportunity to generate data on cell clustering that could be of interest for further studies and may promote new therapeutic ideas. Graph theory offers a framework for developing tools to identify clusters within cell populations. A key step in this process is defining relationships (“similarities”) between pairs of cells, which creates an initial large graph. This pairing can be based on tracking individual cells over time. The next step involves partitioning the initial large graph into disjoint subgraphs (clusters). A third step is to find similarities between consecutive clusters. This study demonstrates an example of applying this approach.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

MK: Conceptualization, Data curation, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. HR: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – review & editing. KT: Conceptualization, Funding acquisition, Investigation, Writing – review & editing. RK: Conceptualization, Data curation, Formal analysis, Methodology, Software, Visualization, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by Olav Raagholt and Gerd Meidel Raagholts legacy, Gidske and Peter Jacob Sorensen Research Foundation, Astri and Birger Torsteds legacy, Oslo University Hospital, and University of Oslo.

Conflict of interest

MK is the owner of the upstart firm Korsnes Biocomputing (KoBio) aimed to participate in research and development of methods for single-cell analysis.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

1. ^https://www.korsnesbiocomputing.no/

2. ^https://www.sciencedirect.com/topics/computer-science/clustering-graph

References

Amintas, S., Bedel, A., Moreau-Gaudry, F., Boutin, J., Buscail, L., Merlio, J.-P., et al. (2020). Circulating tumor cell clusters: united we stand divided we fall. Int. J. Mol. Sci. 21:2653. doi: 10.3390/ijms21072653

PubMed Abstract | Crossref Full Text | Google Scholar

Andrei, L., Kasas, S., Garrido, I. O., Stanković, T., Korsnes, M. S., Vaclavikova, R., et al. (2020). Advanced technological tools to study multidrug resistance in cancer. Drug Resist. Updat. 48:100658. doi: 10.1016/j.drup.2019.100658

PubMed Abstract | Crossref Full Text | Google Scholar

Angermueller, C., Clark, S. J., Lee, H. J., Macaulay, I. C., Teng, M. J., Hu, T. X., et al. (2016). Parallel single-cell sequencing links transcriptional and epigenetic heterogeneity. Nat. Methods 13, 229–232. doi: 10.1038/nmeth.3728

PubMed Abstract | Crossref Full Text | Google Scholar

Archetti, M., and Pienta, K. J. (2019). Cooperation among cancer cells: applying game theory to cancer. Nat. Rev. Cancer 19, 110–117. doi: 10.1038/s41568-018-0083-7

PubMed Abstract | Crossref Full Text | Google Scholar

Au, S. H., Storey, B. D., Moore, J. C., Tang, Q., Chen, Y.-L., Javaid, S., et al. (2016). Clusters of circulating tumor cells traverse capillary-sized vessels. Proc. Natl. Acad. Sci. U. S. A. 113, 4947–4952. doi: 10.1073/pnas.1524448113

PubMed Abstract | Crossref Full Text | Google Scholar

Bakst, R. L., and Wong, R. J. (2016). Mechanisms of perineural invasion. J. Neurol. Surg. B Skull Base 77, 96–106. doi: 10.1055/s-0036-1571835

PubMed Abstract | Crossref Full Text | Google Scholar

Bithi, S. S., and Vanapalli, S. A. (2017). Microfluidic cell isolation technology for drug testing of single tumor cells and their clusters. Sci. Rep. 7:41707. doi: 10.1038/srep41707

PubMed Abstract | Crossref Full Text | Google Scholar

Braadland, P. R., Ramberg, H., Grytli, H. H., and Taskén, K. A. (2015). β-adrenergic receptor signaling in prostate cancer. Front. Oncol. 4:375. doi: 10.3389/fonc.2014.00375

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, S.-H., Zhang, B.-Y., Zhou, B., Zhu, C.-Z., Sun, L.-Q., and Feng, Y.-J. (2019). Perineural invasion of cancer: a complex crosstalk between cells and molecules in the perineural niche. Am. J. Cancer Res. 9:1.

PubMed Abstract | Google Scholar

Cheung, K. J., Gabrielson, E., Werb, Z., and Ewald, A. J. (2013). Collective invasion in breast cancer requires a conserved basal epithelial program. Cell 155, 1639–1651. doi: 10.1016/j.cell.2013.11.029

PubMed Abstract | Crossref Full Text | Google Scholar

Chung, F. R. (1997). Spectral Graph Theory, Vol. 92. Providence, RI: American Mathematical Soc.

Google Scholar

Colak-Champollion, T., Lan, L., Jadhav, A. R., Yamaguchi, N., Venkiteswaran, G., Patel, H., et al. (2019). Cadherin-mediated cell coupling coordinates chemokine sensing across collectively migrating cells. Curr. Biol. 29, 2570–2579. doi: 10.1016/j.cub.2019.06.061

PubMed Abstract | Crossref Full Text | Google Scholar

Commander, R., Wei, C., Sharma, A., Mouw, J., Burton, L., Summerbell, E., et al. (2020). Subpopulation targeting of pyruvate dehydrogenase and glut1 decouples metabolic heterogeneity during collective cancer cell invasion. Nat. Commun. 11:1533. doi: 10.1038/s41467-020-15219-7

PubMed Abstract | Crossref Full Text | Google Scholar

Deborde, S., Gusain, L., Powers, A., Marcadis, A., Yu, Y., Chen, C.-H., et al. (2022). Reprogrammed schwann cells organize into dynamic tracks that promote pancreatic cancer invasion. Cancer Discov. 12, 2454–2473. doi: 10.1158/2159-8290.CD-21-1690

PubMed Abstract | Crossref Full Text | Google Scholar

Deborde, S., Omelchenko, T., Lyubchik, A., Zhou, Y., He, S., McNamara, W. F., et al. (2016). Schwann cells induce cancer cell dispersion and invasion. J. Clin. Invest. 126, 1538–1554. doi: 10.1172/JCI82658

PubMed Abstract | Crossref Full Text | Google Scholar

Demir, I. E., Friess, H., and Ceyhan, G. O. (2012). Nerve-cancer interactions in the stromal biology of pancreatic cancer. Front. Physiol. 3:97. doi: 10.3389/fphys.2012.00097

PubMed Abstract | Crossref Full Text | Google Scholar

Ershov, D., Phan, M.-S., Pylvänäinen, J. W., Rigaud, S. U., Le Blanc, L., Charles-Orszag, A., et al. (2022). Trackmate 7: integrating state-of-the-art segmentation algorithms into tracking pipelines. Nat. Methods 19, 829–832. doi: 10.1038/s41592-022-01507-1

PubMed Abstract | Crossref Full Text | Google Scholar

Freckmann, E. C., Sandilands, E., Cumming, E., Neilson, M., Román-Fernández, A., Nikolatou, K., et al. (2022). Traject3d allows label-free identification of distinct co-occurring phenotypes within 3D culture by live imaging. Nat. Commun. 13:5317. doi: 10.1038/s41467-022-32958-x

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, X.-l., Wu, J.-s., Cao, M.-x., Gao, S.-y., Cen, X., Jiang, Y.-p., et al. (2017). Cytokeratin-14 contributes to collective invasion of salivary adenoid cystic carcinoma. PLoS ONE 12:e0171341. doi: 10.1371/journal.pone.0171341

PubMed Abstract | Crossref Full Text | Google Scholar

Gkountela, S., Castro-Giner, F., Szczerba, B. M., Vetter, M., Landin, J., Scherrer, R., et al. (2019). Circulating tumor cell clustering shapes dna methylation to enable metastasis seeding. Cell 176, 98–112. doi: 10.1016/j.cell.2018.11.046

PubMed Abstract | Crossref Full Text | Google Scholar

Gómez-de Mariscal, E., Grobe, H., Pylvänäinen, J. W., Xénard, L., Henriques, R., Tinevez, J.-Y., et al. (2024). CellTracksColab is a platform that enables compilation, analysis, and exploration of cell tracking data. PLoS Biol. 22:e3002740. doi: 10.1371/journal.pbio.3002740

PubMed Abstract | Crossref Full Text | Google Scholar

Gopinathan, A., and Gov, N. S. (2019). Cell cluster migration: Connecting experiments with physical models. Semin. Cell Dev. Biol. 93, 77–86. doi: 10.1016/j.semcdb.2018.09.009

PubMed Abstract | Crossref Full Text | Google Scholar

Grigore, A. D., Ben-Jacob, E., and Farach-Carson, M. C. (2015). Prostate cancer and neuroendocrine differentiation: more neuronal, less endocrine? Front. Oncol. 5:37. doi: 10.3389/fonc.2015.00037

PubMed Abstract | Crossref Full Text | Google Scholar

Haeger, A., Alexander, S., Vullings, M., Kaiser, F. M., Veelken, C., Flucke, U., et al. (2020). Collective cancer invasion forms an integrin-dependent radioresistant niche. J. Exp. Med. 217:1184. doi: 10.1084/jem.20181184

PubMed Abstract | Crossref Full Text | Google Scholar

Hannon, B., and Ruth, M. (2014). Biosynchronicity (Cham: Springer International Publishing), 241–245.

Google Scholar

Hong, B., and Zu, Y. (2013). Detecting circulating tumor cells: current challenges and new trends. Theranostics 3:377. doi: 10.7150/thno.5195

PubMed Abstract | Crossref Full Text | Google Scholar

Hou, J.-M., Krebs, M. G., Lancashire, L., Sloane, R., Backen, A., Swain, R. K., et al. (2012). Clinical significance and molecular characteristics of circulating tumor cells and circulating tumor microemboli in patients with small-cell lung cancer. J. Clin. Oncol. 30, 525–532. doi: 10.1200/JCO.2010.33.3716

PubMed Abstract | Crossref Full Text | Google Scholar

Huda, S., Weigelin, B., Wolf, K., Tretiakov, K. V., Polev, K., Wilk, G., et al. (2018). Lévy-like movement patterns of metastatic cancer cells revealed in microfabricated systems and implicated in vivo. Nat. Commun. 9:4539. doi: 10.1038/s41467-018-06563-w

PubMed Abstract | Crossref Full Text | Google Scholar

Ioannides, A. A., Kourouyiannis, C., Karittevlis, C., Liu, L., Michos, I., Papadopoulos, M., et al. (2022). A unified framework for analyzing complex systems: Juxtaposing the (kernel) pca method and graph theory. Front. Appl. Math. Stat. 8:947053. doi: 10.3389/fams.2022.947053

Crossref Full Text | Google Scholar

Kaarijärvi, R., Kaljunen, H., and Ketola, K. (2021). Molecular and functional links between neurodevelopmental processes and treatment-induced neuroendocrine plasticity in prostate cancer progression. Cancers 13:692. doi: 10.3390/cancers13040692

PubMed Abstract | Crossref Full Text | Google Scholar

Khalil, A. A., Ilina, O., Gritsenko, P. G., Bult, P., Span, P. N., and Friedl, P. (2017). Collective invasion in ductal and lobular breast cancer associates with distant metastasis. Clin. Exp. Metastasis 34, 421–429. doi: 10.1007/s10585-017-9858-6

PubMed Abstract | Crossref Full Text | Google Scholar

Khalil, A. A., Ilina, O., Vasaturo, A., Venhuizen, J.-H., Vullings, M., Venhuizen, V., et al. (2020). Collective invasion induced by an autocrine purinergic loop through connexin-43 hemichannels. J. Cell Biol. 219:11120. doi: 10.1083/jcb.201911120

PubMed Abstract | Crossref Full Text | Google Scholar

Konen, J., Summerbell, E., Dwivedi, B., Galior, K., Hou, Y., Rusnak, L., et al. (2017). Image-guided genomics of phenotypically heterogeneous populations reveals vascular signalling during symbiotic collective cancer invasion. Nat. Commun. 8:15078. doi: 10.1038/ncomms15078

PubMed Abstract | Crossref Full Text | Google Scholar

Korsnes, M. S., and Korsnes, R. (2015). Lifetime distributions from tracking individual BC3H1 cells subjected to yessotoxin. Front. Bioeng. Biotechnol. 3:166. doi: 10.3389/fbioe.2015.00166

PubMed Abstract | Crossref Full Text | Google Scholar

Korsnes, M. S., and Korsnes, R. (2018). Single-cell tracking of a549 lung cancer cells exposed to a marine toxin reveals correlations in pedigree tree profiles. Front. Oncol. 8:260. doi: 10.3389/fonc.2018.00260

PubMed Abstract | Crossref Full Text | Google Scholar

Korsnes, M. S., and Korsnes, R. (2023). Initial refinement of data from video-based single-cell tracking. Cancer Innov. 2, 416–432. doi: 10.1002/cai2.88

PubMed Abstract | Crossref Full Text | Google Scholar

Korsnes, M. S., Ramberg, H., Taskén, K. A., and Korsnes, R. (2024a). PC3 cell dance video. Mendeley Data. 2024:1. doi: 10.17632/k93k4z48rk.1

Crossref Full Text | Google Scholar

Korsnes, M. S., Ramberg, H., Taskén, K. A., and Korsnes, R. (2024b). PC3 cell in co-culture video. Mendeley Data. 2024:1. doi: 10.17632/fkzmhznghs.1

Crossref Full Text | Google Scholar

Korsnes, M. S., Ramberg, H., Taskén, K. A., and Korsnes, R. (2024c). PC3 cell in mono-culture video. Mendeley Data. 2024:1. doi: 10.17632/fb97xkb6vx.1

Crossref Full Text | Google Scholar

Kumar, P., Tan, Y., and Cahan, P. (2017). Understanding development and stem cells using single cell-based analyses of gene expression. Development 144, 17–32. doi: 10.1242/dev.133058

PubMed Abstract | Crossref Full Text | Google Scholar

Lee, M., Kim, E. J., Cho, Y., Kim, S., Chung, H. H., Park, N. H., et al. (2017). Predictive value of circulating tumor cells (CTCS) captured by microfluidic device in patients with epithelial Ovarian cancer. Gynecol. Oncol. 145, 361–365. doi: 10.1016/j.ygyno.2017.02.042

PubMed Abstract | Crossref Full Text | Google Scholar

Ling, H., Mclvor, G. E., Westley, J., van der Vaart, K., Vaughan, R. T., Thornton, A., et al. (2019). Behavioural plasticity and the transition to order in jackdaw flocks. Nat. Commun. 10:5174. doi: 10.1038/s41467-019-13281-4

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, C. C., Greenwald, N. F., Kong, A., McCaffrey, E. F., Leow, K. X., Mrdjen, D., et al. (2023). Robust phenotyping of highly multiplexed tissue imaging data using pixel-level clustering. Nat. Commun. 14:4618. doi: 10.1101/2022.08.16.504171v1

PubMed Abstract | Crossref Full Text | Google Scholar

Magnon, C. (2015). Role of the autonomic nervous system in tumorigenesis and metastasis. Mol. Cell. Oncol. 2:e975643. doi: 10.4161/23723556.2014.975643

PubMed Abstract | Crossref Full Text | Google Scholar

Maity, D., Sivakumar, N., Kamat, P., Zamponi, N., Min, C., Du, W., et al. (2024). Profiling dynamic patterns of single-cell motility. bioRxiv. doi: 10.1101/2022.09.21.508955v2

Crossref Full Text | Google Scholar

March, B., Faulkner, S., Jobling, P., Steigler, A., Blatt, A., Denham, J., et al. (2020). Tumour innervation and neurosignalling in prostate cancer. Nat. Rev. Urol. 17, 119–130. doi: 10.1038/s41585-019-0274-3

PubMed Abstract | Crossref Full Text | Google Scholar

Merino-Casallo, F., Gomez-Benito, M. J., Hervas-Raluy, S., and Garcia-Aznar, J. M. (2022). Unravelling cell migration: defining movement from the cell surface. Cell Adhes. Migrat. 16, 25–64. doi: 10.1080/19336918.2022.2055520

PubMed Abstract | Crossref Full Text | Google Scholar

Molnar, B., Ladanyi, A., Tanko, L., Sréter, L., and Tulassay, Z. (2001). Circulating tumor cell clusters in the peripheral blood of colorectal cancer patients. Clin. Cancer Res. 7, 4080–4085.

Google Scholar

Mosier, J. A., Wu, Y., and Reinhart-King, C. A. (2021). Recent advances in understanding the role of metabolic heterogeneities in cell migration. Facul. Rev. 10:8. doi: 10.12703/r/10-8.eCollection

Crossref Full Text | Google Scholar

Nascimento, M. C., and De Carvalho, A. C. (2011). Spectral methods for graph clustering–a survey. Eur. J. Operat. Res. 211, 221–231. doi: 10.1016/j.ejor.2010.08.012

Crossref Full Text | Google Scholar

Niu, B., and Wang, H. (2012). Bacterial colony optimization. Discr. Dyn. Nat. Soc. 2012:698057. doi: 10.1155/2012/698057

Crossref Full Text | Google Scholar

Qian, H., Sheetz, M. P., and Elson, E. L. (1991). Single particle tracking. Analysis of diffusion and flow in two-dimensional systems. Biophys. J. 60, 910–921.

Google Scholar

Quinsgaard, E. M. B., Korsnes, M. S., Korsnes, R., and Moestue, S. A. (2024). Single-cell tracking as a tool for studying emt-phenotypes. Exp. Cell Res. 437:113993. doi: 10.1016/j.yexcr.2024.113993

PubMed Abstract | Crossref Full Text | Google Scholar

Reynolds, A. M. (2018). Current status and future directions of lévy walk research. Biol. Open 7:bio030106. doi: 10.1242/bio.030106

PubMed Abstract | Crossref Full Text | Google Scholar

Rozenberg, J. M., Buzdin, A. A., Mohammad, T., Rakitina, O. A., Didych, D. A., Pleshkan, V. V., et al. (2023). Molecules promoting circulating clusters of cancer cells suggest novel therapeutic targets for treatment of metastatic cancers. Front. Immunol. 14:1099921. doi: 10.3389/fimmu.2023.1099921

PubMed Abstract | Crossref Full Text | Google Scholar

Schaeffer, S. E. (2007). Graph clustering. Comput. Sci. Rev. 1, 27–64. doi: 10.1016/j.cosrev.2007.05.001

Crossref Full Text | Google Scholar

Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682. doi: 10.1038/nmeth.2019

PubMed Abstract | Crossref Full Text | Google Scholar

Schmitd, L., Scanlon, C., and D'silva, N. (2018). Perineural invasion in head and neck cancer. J. Dent. Res. 97, 742–750. doi: 10.1177/0022034518756297

PubMed Abstract | Crossref Full Text | Google Scholar

Shannon, M. J., Eisman, S. E., Lowe, A. R., Sloan, T. F., and Mace, E. M. (2024). cellPLATO–an unsupervised method for identifying cell behaviour in heterogeneous cell trajectory data. J. Cell Sci. 137:564355. doi: 10.1101/2023.10.28.564355

PubMed Abstract | Crossref Full Text | Google Scholar

Stauffer, D., and Aharony, A. (2018). Introduction to Percolation Theory. New York, NY: Taylor & Francis.

Google Scholar

Summerbell, E. R., Mouw, J. K., Bell, J. S., Knippler, C. M., Pedro, B., Arnst, J. L., et al. (2020). Epigenetically heterogeneous tumor cells direct collective invasion through filopodia-driven fibronectin micropatterning. Sci. Adv. 6:eaaz6197. doi: 10.1126/sciadv.aaz6197

PubMed Abstract | Crossref Full Text | Google Scholar

Ulicna, K., Vallardi, G., Charras, G., and Lowe, A. R. (2021). Automated deep lineage tree analysis using a bayesian single cell tracking approach. Front. Comput. Sci. 3:734559. doi: 10.3389/fcomp.2021.734559

Crossref Full Text | Google Scholar

Von Luxburg, U. (2007). A tutorial on spectral clustering. Stat. Comput. 17, 395–416. doi: 10.48550/arXiv.0711.0189

Crossref Full Text | Google Scholar

Wen, C., and Kimura, K. D. (2022). Tracking moving cells in 3D time lapse images using 3DeeCellTracker. Bio-protocol 12:e4319. doi: 10.21769/BioProtoc.4319

PubMed Abstract | Crossref Full Text | Google Scholar

Wiggins, L., Lord, A., Murphy, K. L., Lacy, S. E., O'Toole, P. J., Brackenbury, W. J., et al. (2023). The cellPhe toolkit for cell phenotyping using time-lapse imaging and pattern recognition. Nat. Commun. 14:1854. doi: 10.1038/s41467-023-37447-3

PubMed Abstract | Crossref Full Text | Google Scholar

Wrenn, E., Huang, Y., and Cheung, K. (2021). Collective metastasis: coordinating the multicellular voyage. Clin. Exp. Metastasis 38, 373–399. doi: 10.1007/s10585-021-10111-0

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, C., Cao, M., Liu, Y., He, Y., Du, Y., Zhang, G., et al. (2019). Inducible formation of leader cells driven by CD44 switching gives rise to collective invasion and metastases in luminal breast carcinomas. Oncogene 38, 7113–7132. doi: 10.1038/s41388-019-0899-y

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, M., Bardia, A., Wittner, B. S., Stott, S. L., Smas, M. E., Ting, D. T., et al. (2013). Circulating breast tumor cells exhibit dynamic changes in epithelial and mesenchymal composition. Science 339, 580–584. doi: 10.1126/science.1228522

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, J., Goliwas, K. F., Wang, W., Taufalele, P. V., Bordeleau, F., and Reinhart-King, C. A. (2019). Energetic regulation of coordinated leader–follower dynamics during collective invasion of breast cancer cells. Proc. Natl. Acad. Sci. U. S. A. 116, 7867–7872. doi: 10.1073/pnas.1809964116

PubMed Abstract | Crossref Full Text | Google Scholar

Zoeller, E. L., Pedro, B., Konen, J., Dwivedi, B., Rupji, M., Sundararaman, N., et al. (2019). Genetic heterogeneity within collective invasion packs drives leader and follower cell phenotypes. J. Cell Sci. 132:jcs231514. doi: 10.1242/jcs.231514

PubMed Abstract | Crossref Full Text | Google Scholar