The effects of coating culture dishes with collagen on fibroblast cell shape and swirling pattern formation

Motile human-skin fibroblasts form macroscopic swirling patterns when grown to confluence on a culture dish. In this paper, we investigate the effect of coating the culture-dish surface with collagen on the resulting pattern, using human-skin fibroblast NB1RGB cells as the model system. The presence of the collagen coating is expected to enhance the adherence of the fibroblasts to the dish surface, and thereby also enhance the traction that the fibroblasts have as they move. We find that, contrary to our initial expectation, the coating does not significantly affect the motility of the fibroblasts. Their eventual number density at confluence is also unaffected. However, the coherence length of cell orientation in the swirling pattern is diminished. We also find that the fibroblasts cultured in collagen-coated dishes are rounder in shape and shorter in perimeter, compared with those cultured in uncoated polystyrene or glass culture dishes. We hypothesise that the rounder cell-shape which weakens the cell–cell nematic contact interaction is responsible for the change in coherence length. A simple mathematical model of the migrating fibroblasts is constructed, which demonstrates that constant motility with weaker nematic interaction strength does indeed lead to the shortening of the coherence length. Electronic supplementary material The online version of this article (10.1007/s10867-020-09556-3) contains supplementary material, which is available to authorized users.

observed in single-cell organisms. A well-known example is aggregation by which Dictyostelium cells form a slug-like structure when starved [3]. Deciphering the mechanisms that drive robust and precise collective migration of a large number of cells is of vital importance in understanding development, differentiation, and evolution, with many possible applications in cancer therapies, regenerative medicine, and tissue engineering [2,[4][5][6].
In vitro cultivation studies can provide important insights into collective cell migration. When cultivated densely, complex alignment patterns are known to spontaneously form in several types of cells. The characteristic coherence length of the resulting alignment pattern is of special interest since it provides a measure of the number of cells that can migrate collectively. One of the key factors determining the coherence length is the strength of the cell-cell contact interaction. Although a migrating cell has vectorial polarity, moving toward a specific heading, the alignment of these cells is often nematic, that is, neighbouring cells tend to migrate either in parallel or anti-parallel directions. This nematic alignment can be observed in several cellular-scale objects including gliding microtubules, actin filaments, bacteria, and cultured cells [7][8][9][10][11][12][13][14]. Nematic alignment is also observed in purely mechanical systems whose members interact via the excluded-volume effect, such as in a population of rod-like objects [15]. In this case, the strength of the nematic interaction is determined by the shape of the rods; the interaction being stronger between longer rods.
In the present paper, we investigate the alignment pattern of human-skin fibroblast NB1RGB cells grown to confluence in a culture dish. Skin fibroblasts are cells that provide structural support for the skin and are easily cultivated in vitro. They represent a convenient model system for studying the collective migration of cells. Here, we examine how the alignment pattern of the fibroblasts depends on whether or not the dish surface is coated with collagen. To check for possible dependence of the result on the substrate material, the experiment is performed on two types of dishes: polystyrene and glass.
Since fibroblasts adhere to collagen via integrins, a natural expectation would be that the collagen-coated surface would provide enhanced traction to the fibroblasts leading to their enhanced motility, which in turn would lead to a change in the alignment pattern. In vitro, the collagenous extracellular matrix (ECM) is known to stimulate skin-fibroblast motility [16], which adds support to this expectation.
As reported in [12], when fibroblasts are seeded onto a culture dish at low density, the cells adhere to the dish individually and then migrate randomly into cell-free areas while only occasionally coming into contact with other cells. As they proliferate, their density increases and confluence is eventually achieved. Confluent fibroblasts align locally along their elongated axes and form macroscopic swirling patterns (see Supplementary Videos 1 and 2). A similar swirling pattern, the storiform pattern, is often observed in fibrohistiocytic lesions in vivo [17,18].
We find that for both polystyrene and glass dishes, the characteristic coherence length of the swirling pattern decreases as the density of the collagen coated onto the dish surface is increased. Moreover, we observe that the cells become rounder in shape with the increase in coated-collagen density, whereas cell-number density and, unexpectedly, cell motility remain unchanged. From these experimental results, we hypothesise that the difference in the coherence length mainly follows from the difference in the strength of the nematic contact interaction between the cells; i.e., rounder cells experience weaker nematic interactions, and thereby the coherence length becomes shorter. To test the feasibility of this hypothesis, we construct a simple mathematical model of migrating cells in which the strength of the nematic interaction can be controlled by a single parameter. Numerical simulations of this model demonstrate that the coherence length does indeed correlate positively with the nematic interaction strength. We thus propose that the collagen coating first leads to the change in the fibroblast cell shape, which in turn shortens the coherence length. Possible mechanisms by which the collagen coating leads to the rounding of the fibroblasts are also discussed.

Human-skin fibroblast NB1RGB cells form macroscopic swirling patterns
Human-skin fibroblast NB1RGB cells form macroscopic swirling patterns when cultivated in a culture dish due to their elongated shape, motility, and proliferation within the confined twodimensional surface. Our objective is to quantify the difference in the swirling patterns when the fibroblasts are cultured on uncoated and collagen-type-I-coated dishes. As discussed above, we perform the experiment on both polystyrene-and glass-bottom dishes. The collagen-coating procedure is detailed in Sect. 4. The resulting collagen density on the surface of the culture dish depends on the concentration of collagen type I in the initial collagen solution used in the coating process. We use the coating obtained from a 10.0-μg/mL collagen type-I solution as the standard reference.
Typical results for uncoated and collagen-coated dishes are shown in Fig. 1. See also Supplementary Videos 1 (uncoated) and 2 (collagen-coated). In the first row of Fig. 1a, still images of the confluent fibroblast swirling patterns on uncoated (left) and collagen-coated (right) polystyrene dishes are shown. The black scale bar in the lower-right corner of the right image is 1 mm long. For ease of comparison, the orientations of the individual fibroblasts in these images are read using Orientation J software [19] and shown colour coded in the second row of Fig. 1a, with light-blue and red respectively indicating horizontal and vertical orientations (see Sect. 4.4 in Sect. 4 for details). Other colours indicate orientations in between, as shown on the scale on the right margin of Fig. 1b.
Comparing the colour-coded images by inspection, one discerns that the fibroblasts form into patches of cells with similar orientation, and that these patches are slightly larger for the uncoated dish compared with the collagen-coated dish.
We quantify this observation by extracting the correlations of cell orientations from the images following the procedure detailed in Sect. 4. First, the 1600 × 1200 pixel image is divided into a 50 × 38 grid, each subdivision being 32 × 32 pixels in size. Then, the blockaveraged orientation is calculated for each subdivision, the results of which are shown in Fig.  1b. The correlations between the block-averaged orientations for each separation of the blocks is then calculated.
The resulting correlation functions are plotted in Fig. 1c (uncoated polystyrene) and d (collagen-coated polystyrene) for the images shown in Fig. 1a. A blow-up of the graphs between 0 and 1 mm is shown in the subframe inside Fig. 1d (c' and d'). We can see that the orientation correlation falls off more quickly for the collagen-coated dish compared with the uncoated dish. This difference is consistently reproduced for multiple dishes, over multiple repetitions of the experiment, for both dish materials.
To characterise the fall-off of the correlation function with distance, we define the length d 50 to be the distance in which the correlation is reduced to one half (50%) of the maximum value (100%). How this length is determined is illustrated in Fig. 1c

c d
A clear trend can be seen in Fig. 1e in which the orientation coherence length decreases with the initial collagen concentration, indicating that the cell orientations vary more rapidly with distance when the collagen coatings on the culture dishes are denser, resulting in the cells forming more swirling patterns overall. Figure S3B, though having only two data points, confirms this trend.
The culture dishes in this study are coated with collagen molecules dissolved in acidic conditions and not with collagen fibrils. To test whether the structure of the collagen coating affects the coherence length, we culture the cells in dishes coated with various types of collagen, all starting from an acidic (pH 3.0) solution of concentration 10.0 μg/mL. These are: 1. Collagen type I (the reference) 2. Heat-denatured collagen type I 3. Gelatin, which comprise thermal denatured collagen fibrils, 4. Collagen type IV, which is a non-fibrillar form of collagen.
In all of these cases, the values of d 50 are decreased compared with those of the control (uncoated). See Fig. S2A-F in the Supplementary information. This indicates that the reduction of orientation coherence occurs independently of the structure of collagen.

Collagen coating affects cell shape, but not single-cell motility or cell-number density
To clarify the effects of the collagen coating on individual fibroblast cells, we analyse cell motility and cell-number density, both of which are likely to affect the alignment pattern [20,21]. To investigate cell motility, we track the movements of individual cells cultured at low and high densities (Fig. 2a). This measurement is possible for glass dishes only due to limitations of our equipment. At low density, the instantaneous velocity is slightly smaller for the collagen-coated dish compared with the uncoated dish but not by a significant amount. At high cell density, the instantaneous velocities are indistinguishable between the uncoated and collagen-coated cases. A similar tendency is observed in the 10 haverage velocity (Fig. S3D, E), suggesting that directional persistence in the cell motility has little dependence, if any, on the dish coating. Moreover, we measure the number of cells in the high-density cultures, for both dish materials, and find no significant difference in the numbers of cells ( Fig. 2b; Fig. and S3C).
We also analyse the effects of the collagen coating on the fibroblast cell shape. We measure the area S and perimeter L of each cell on uncoated and 0.1, 1.0, and 10.0 μg/mL collagen Fig. 1 Patterns formed by human-skin NB1RGB fibroblast cells cultured in a polystyrene dish. a Typical patterns of confluent fibroblasts cultured for 144 h on the uncoated control (left) and the dish coated from the 10.0 μg/mL collagen type-I solution (right). The black scale bar in the lower-right corner of the rightmonochrome image is 1 mm long. The smaller images in the right margin are dishes coated from 0.1 μg/mL (top) and 1.0 μg/mL (bottom) collagen type-I solutions. For each pair of images, the lower colour-map encodes the orientations of the fibroblasts in the upper bright-field image in accordance with the colour scale shown to the right of (b). b The block-averaged orientations of the fibroblasts on uncoated (left) and 10.0 μg/mL collagen-coated (right) dishes shown in (a). c, d Correlation functions of cell orientation for the uncoated (c) and collagen-coated (d) cases. The graphs c' and d' in the sub-window show the behaviour of the two functions near the origin on the same axes for the ease of comparison. e d 50 values averaged over 40 images. Data represent the mean ± standard error of mean (SEM) from four independent cultures. Ten images were captured from each culture dish. ***p < 0.001 under Student's one-tailed t test when compared with the control type-I-coated polystyrene dishes at low density (Fig. 2c, f), and quantify the cell roundness in terms of the circularity 4πS/L 2 (Fig. 2e). The comparison between uncoated and 10.0-μg/mL collagen type-I-coated glass dishes is shown in Fig. S3F, G in the Supplementary information.
These measurements indicate that the NB1RGB fibroblasts at low cell density become rounder and smaller as the density of the collagen on the coating is increased, for both polystyrene and glass dishes. The same tendency is observed for heat-denatured collagen type I, gelatin, and collagen type IV on polystyrene (Fig. S2G, H), which suggests that interactions Samples are compared with Student's t test and labelled n.s. (no significance); *p < 0.05; **p < 0.01; ***p < 0.001. a Migration velocities of the fibroblasts cultured at low and high densities on uncoated and 10.0μg/mL collagen-type-I-coated glass dishes. Each data point represents the mean ± SEM of 20 cells chosen randomly from two cell cultures. b Number of fibroblasts cultured for 24-144 h on uncoated and 10-μg/mL collagen-type-I-coated polystyrene dishes. Data points represent the mean ± SEM of four dishes. c-f Analysis of fibroblast cell morphology. c, d Fibroblasts cultured on uncoated and 10-μg/mL collagen-type-I-coated polystyrene dishes at 24 h (c; low density) and 72 h after culture start (d; high density). Scale bar, 100 μm. The fibroblasts were observed and photographed with a Nikon ECLIPSE TS 100 phase-contrast microscope (NIKON Corp., Tokyo, Japan). e, f Circularity and area of fibroblasts cultured on uncoated and 0.1-, 1.0-, and 10.0-μg/mL collagen-type-I-coated polystyrene dishes at low density. Each data point represents the mean ± SEM of 20 cells chosen randomly from two cell cultures with the collagen coating make the shape of NB1RGB fibroblasts rounder and smaller. At high cell density, it is difficult to quantify the cell shape since the cell boundaries become difficult to discern (Fig. 2d). Regarding cell area, there should be no significant difference between the uncoated and 10.0 μg/mL collagen type-I-coated dishes at high cell density given that the final cell densities are indistinguishable.

Simulation: interaction strength increases coherence length
Among our observations, the only significant difference between the cells cultured in coated and uncoated dishes is in their shapes at low density. Although it is possible that the difference in the cell shape is not maintained at high density, the shape of isolated single cells is likely to affect the orientation dynamics in highly dense cell populations.
Motivated by previous studies on the collective dynamics of rod-like objects [15], we consider the hypothesis that the nematic interactions between the cells in collagen-coated dishes are weaker than those in the uncoated dishes due to the cells becoming rounder in the presence of collagen, and it is this weaker nematic interaction that leads to the shorter coherence length.
The viability of this hypothesis depends on whether varying the strength of nematic interactions alone is sufficient in changing the coherence length. To this end, we construct a simple mathematical model with just a few tunable parameters, as described in Sect. 4. The model includes effects of cell motility, cell proliferation, excluded volume effect, and nematic interactions. We run numerical simulations of our model to assess the effect of different nematic interaction strengths on the cell alignment patterns.
Despite its simplicity, our model reproduces the experimentally observed dynamical behaviour well (see Supplementary Videos 3 and 4). Figure 3 shows some typical numerical results, where the parameter K (in units of hour −1 ) quantifies the strength of the nematic interaction. Figure 3a shows the alignment pattern obtained for K = 0.020/h, after running the simulation from a random initial condition for an equivalent of 144 h. Figure 3c shows its block-averaged cell orientations, and Fig. 3e the corresponding distance dependence of the orientation correlation. Figure 3b, d, and f shows those for K = 0.015/h. As shown in Fig. 3g, an increase in the nematic interaction strength K, with all other model parameters kept fixed, results in an increase in the coherence length. Figure 3g also shows how the K-dependence of d 50 changes with the cell migration velocity v. The three values of v chosen are the experimentally observed values at low cell density with and without collagen coating (Fig. 2a) and their average. There is little variation, demonstrating that our conclusion is independent of v. We also note that the observed values of d 50 develop rather late in the simulation as shown in Fig. 3h.
To gain deeper understanding of the pattern formation process, we perform further numerical simulations under various conditions. The results are shown in Fig. 4 and Supplementary Videos 5-7. In all of the cases except for K = 0, we confirm the formation of swirling patterns and the increase of d 50 values with the strength of nematic interaction.
Recall that in Fig. 3h, d 50 values grow rapidly at late times. In contrast, in the "from high density" case in Fig. 4b, d 50 grows from earlier times. These observations suggest that the growth of the coherence length starts after the cell density becomes sufficiently high.
To check the robustness of our results, we further perform additional numerical simulations in which noise is included in the angular dynamics (Fig. 4c, d and Supplemental movies 6 and 7). We introduce white noise of the strength μ (with units hour −1 ) in the angular dynamics of each cell. The correlation time of angular dynamics of single isolated cells is approximately given by 1/μ, i.e., the direction of a cell at time t is approximately independent of that at time t − 1/μ. We choose a reference value of μ to be 0.1/h in our simulations. As seen  Fig. 3 with v = 0.1864 units/h, (ii) "zero speed", where v = 0, (iii) "random walk", where the migration direction of each cell at each time step is set to be a random value irrespective of its orientation, (iv) "from high density", where the cell proliferation is absent and the initial number of cells is 14,080, (v) "random orientation", where the orientation of the newborn cell is assigned a random value, (vi) "zero speed, from high density, no excluded volume effect", where we employ (ii) and (iv) and further assume K C = 0. c d 50 vs. K at t = 144 h with noise strength μ = 0.1/h. d d 50 vs. μ at t = 144 h for K = 0.2/h in the movies, swirling patterns develop even in the presence of noise though the strength of nematic interaction required is larger compared with the noiseless case. In Fig. 4c, it is observed that substantial coherence emerges when the coupling strength K is comparable with or larger than the noise strength μ = 0.1/h, and the coherence length increases with the coupling strength K. Conversely, Fig. 4d shows the dependence of d 50 on the noise strength μ with K fixed to 0.20/h, and the value of d 50 can be seen to increase as the noise strength μ is lowered. As there is no discernible experimental difference in the directional persistency between coated and uncoated dishes (Fig. S3d, e), the difference in actual noise level must be small, and it is unlikely such a difference would account for the observed difference between the d 50 values in the fibroblast cultures.
Our simple mathematical model thus demonstrates that the effect of collagen coating on the orientation patterns can be accounted for solely by collagen's effect on the strength of the nematic interactions between the fibroblast cells.

Discussion
Human-skin fibroblasts cultured at high density form macroscopic swirling patterns in the culture dishes. This study reveals that the cells become rounder and form less coherent patterns as the density of the collagen coated on the culture dishes is increased.
There are clear differences in the morphological properties between cells cultured in the uncoated and coated dishes; cells cultured in the coated dishes are shorter in perimeter and rounder (Fig. 2c, e, and f). The following molecular mechanism could underlie this change. It is known that the elongated cells have developed stress fibres; conversely, lamellipodia, which are formed by actin projections on the edge of the cell, rather than a stress fibre, are activated to make the cells round [22]. Two members of the Rho GTPase family, RhoA and Rac1, are necessary for the regulation of various cellular behaviours, including microfilament network organisation. The formation of stress fibre and lamellipodia are induced by RhoA and Rac1, respectively [23][24][25][26][27]. RhoA and Rac1 inhibit each other's activation, and this competition between the two is one of the factors that regulate cell morphogenesis [23,28]. It has been reported that collagen type I increases the activity of Rac1 in human platelets [29]. In fibroblasts, it is unknown whether collagen controls the Rac1 activity, whereas the inhibition of Rac1 activation promotes the expression of collagen protein, suggesting that collagen interacts with Rac1 in fibroblasts [30]. Therefore, it is possible that the collagen coating induces the development of lamellipodia and decreases that of stress fibres via Rac1 activation, thus changing the NB1RGB cell roundness.
Collagen is also known to control fibroblast motility. Li et al. [16] report that collagen coating promotes cell migration into cell-free space in scratch assays. This result seemingly differs from our results shown in Fig. 2a, wherein no significant difference was found between the migration speeds in the uncoated and coated dishes. However, they do not necessarily contradict each other since Li et al. [16] studied a different scenario: it analysed the migration of the cell assembly from densely populated to scratched regions. Moreover, the migration speed of individual cells may be sensitive to its measurement method, such as the discretisation of cell trajectories and the observation time interval. Depending on the method, a significant difference between the uncoated and coated dishes could arise. However, our numerical simulations reveal that the coherence length is insensitive to small variations in migration velocity (Fig. 3g), suggesting that the difference in migration velocities at low density, if any actually exists, does not have a substantial effect on pattern formation at high density.
As the collagen density is increased, the cells become rounder, whereas cell motility remains unchanged, as shown in Figs. 1e and 2e. The same tendency is found for different types of coating materials (collagen type I, heat-denatured collagen type I, gelatin, and collagen type IV), as shown in Fig. S2G. We thus hypothesise that the change in coherence length results from a change in the cell-cell contact interactions mediated by a change in cell shape. However, we do not have any direct evidence of this, and it therefore remains an open issue. Use of genetic manipulation or inhibitors that exclusively control the cell shape would be required for such a study.
Concerning the mathematical models of fibroblast orientation proposed in the literature, some focus on the interplay between the cells and the extra-cellular matrix (ECM) [31,32], known as dynamic reciprocity. Fibroblast orientation has also been modelled in terms of individual cell migration [33]. Systems of reaction-diffusion and integro-partial differential equations have also been used to model fibroblast orientation; however, these require large numbers of parameters and computational complexity that would unnecessarily complicate the isolation of the key target parameter. Instead, our goal in using mathematical modelling is to test the hypothesis that changes in the coherence length of the orientation patterns can be driven solely by changes in the strength of cell-cell interactions. We therefore consider a model in which the strength of the nematic interaction can be controlled by a single parameter. In addition, to keep the model as simple as possible, we regarded each cell as a self-driven particle with constant speed, following various theoretical studies on collective migration [34,35]. Despite the simplicity of our model, it qualitatively reproduces our experimental findings (Fig. 3) and represented rather realistic dynamics of collective migration ( Supplementary  Videos 3 and 4). Our model thus demonstrates that the coherence length increases with the strength of the cell-cell interaction, lending support to our hypothesis. Using our mathematical model, we perform various in silico experiments ( Fig. 4; Supplementary Videos 5 and 6). Similar swirling patterns are observed in all the conditions we employ except for the case of K = 0. In particular, qualitatively the same results are obtained even when spontaneous mobility, reproduction, and excluded volume effect are all turned off. Therefore, we propose that the nematic interaction is the primal factor for the alignment process and the formation of swirling patterns in cell cultures.
When only the nematic interaction is considered, the system is essentially the same as a population of identical oscillators distributed in a two-dimensional space with a synchronisation interaction between close neighbours. The synchronisation process of such a system can be described by the model dθi where ω denotes the natural frequency of each oscillator, b K denotes the interaction strength, and j = 〈i〉 indicates a sum over the nearest neighbours of the i-th oscillator. The parameter m is introduced for convenience to toggle between the two cases of ferromagnetic (m = 1) and nematic (m = 2) interactions. The latter corresponds to our situation, whereas the former is usually considered for a system of interacting oscillators. We may set ω = 0 without loss of generality as it corresponds to the change θ i → θ i − ωt. The parameter m may also be set to unity without loss of generality as it corresponds to the changes mθ i →θ i ; mω→ω; m b K→ b K: Thus, let us assume ω = 0 and m = 1.
Then, it is clear that b K determines only the time scale of the process; i.e., the process becomes faster for larger b K without any other changes. Moreover, the system can be given a variational . Thus, the system evolves with time toward a minimum of G. The global minimum of G corresponds to θ j = θ i for all coupled pairs, suggesting that perfect alignment should eventually be achieved. When a random initial condition for the θ i values is employed, many topological defects typically arise. The coherence length increases with time as the number of defects decrease by annihilation, as observed in various systems (see ref. [37] and references therein). Thus, with a given observation time, a larger coherence length is obtained for larger b K. We suppose that this speed-up effect is a primal mechanism underlying the increase of coherence length with the nematic coupling strength, denoted by K in our model of motile cells. These observations are in line with ref. [12], where the author observes pattern formation in cultures of fibroblasts with elongated cell shapes and finds that the number of swirling patches progressively decreases after confluence and a single parallel array is eventually obtained.
Correlation analysis reveals that the patterns formed under the conditions of this experiment are on the order of several cell-lengths long; however, larger cell assemblies may migrate collectively in vivo [38]. Possible factors that hamper larger-scale coherence in vitro include noisy single-cell migration, cell division, and topological defects [10]. To realise the largescale collective migration in vivo, several other factors may play substantial or complementary roles, such as cell adhesion, spontaneous role assignment among cells, and polarity alignments among adjacent cells [38][39][40]. Mechanical tension applied to the tissue would also contribute to collective migration [41,42].
Our experimental and numerical results support our hypothesis that cell shape affects largescale coherence by mediating the strength of cell-cell contact interactions. Although our finding is based on an in vitro system, such a mechanism may be at work in vivo as well. Thus, our study suggests that cell shape may play an essential role in cell-cell communication in single and multi-cellular organisms.

Coating of culture dish
Collagen type I, collagen type IV, and gelatin were purchased from Nitta Gelatin Inc. (Osaka, Japan. Product names: Cellmatrix Type IA, Cellmatrix Type IV, and GLS250). Collagen types I and IV are provided as acidic solutions of pH 3 and concentration 3 mg/mL. These were diluted with a 1 mM HCl solution (pH 3.0) to the desired concentrations. Gelatin was dissolved in the 1 mM HCl solution to the desired concentration. Maintaining the acidic condition prevents the formation of collagen fibrils from the collagen molecules. The heatdenatured collagen type-I solution was prepared by heating the 10.0 μg/ml Cellmatrix Type-IA solution for 30 min at 60°C.
The polystyrene culture dishes used in this experiment were 100 mm in diameter and obtained from AS ONE (Osaka, Japan). The glass culture dishes were 35 mm in diameter and obtained from Iwaki Co., Ltd. (Tokyo, Japan).
The culture dishes were incubated with the various types of collagen solutions, or just the vehicle solution (1 mM HCl) for the uncoated control, overnight at 4°C. Subsequently, the dishes were air dried and washed four times in phosphate-buffered saline (PBS).
The amount of collagen that was attached to the dish surface was estimated as follows. Following the same procedure as the culture dishes, each well of a 96-well polystyrene plate was incubated with 50 μL of the collagen or vehicle solution overnight at 4°C, then air dried and washed four times in PBS. The amount of collagen attached to the well surface was then measured with the Collagen Quantitation Kit (Cosmo Bio, Tokyo, Japan) following the manufacturer's protocol. By comparing the amount of coated collagen prepared with 1 and 10 μg/mL collagen solutions, we confirmed that the collagen coated on the wells increased approximately in proportion to the dosage (Fig. S1). Note that the detection limit of the Quantitation Kit was 0.4 μg/mL, hence the amount of collagen in the well prepared with the 0.1 μg/mL solution was below the detection limit.

Cell culture
Normal human-skin fibroblasts, RIKEN original (NB1RGB), were provided by the RIKEN BioResource Research Center through the National BioResource Project of MEXT, Japan. The cells were cultured in minimum essential medium alpha (MEMα; Life technologies, Carlsbad, CA) supplemented with 10% fetal bovine serum (FBS; Biowest, Nuaille, France) at 37°C in a humidified incubator with a 5% CO 2 atmosphere.
The 100-mm diameter polystyrene dishes were seeded with 5.0 × 10 5 NB1RGB cells, and incubated for up to 144 h (6 days). The exception was the experiment reported in Fig. S2, which started out with 1.0 × 10 6 cells and incubated for 72 h (3 days), the larger initial cellcount leading to an earlier attainment of confluence. The 35-mm diameter glass dishes were seeded with 6.0 × 10 4 NB1RGB cells so that the initial cell-density will be the same as the 100mm diameter dish with 5.0 × 10 5 cells. These were incubated for up to 90 h. After incubation, the dishes were washed in ice-cold PBS and the cells fixed in ice-cold 100% methanol.

Cell number measurement
For the cell-number density measurements reported in Fig. 2b and Fig. S3C, the cells were dissociated from the dishes with 0.025% Trypsin-EDTA, and the number of cells counted using an automated cell counter (BACKMAN COULTER, Brea, CA).

Orientation analysis
For cell-orientation analyses reported in Fig. 1 and Figs. S2 and S3, images of the swirling patterns were captured at × 50 magnification with a digital microscope (VH-Z20W, KEYENCE, Osaka, Japan). Each image had 1600 × 1200 pixels, corresponding to an area of 6960 × 5220 μm 2 . Since this is quite small compared with the total area of the 100-mm diameter polystyrene dish, images of ten different non-overlapping fields of the dish were collected from each. For the 35-mm diameter glass dishes, with a 12% area compared with the 100-mm dish, one image was taken from each. These images were analysed with the ImageJ plugin OrientationJ [19] to generate the colour-coded maps shown in Fig. 1a and Figs. S2A-E and S3A.
OrientationJ determines the local orientation θ of an image as follows. The 2D monochrome image is essentially an intensity function f(x, y) defined for every pixel (x, y) of the frame. OrientationJ overlays a Gaussian window w(x − x 0 , y − y 0 ) on the field and computes the structure tensor matrix for every (x 0 , y 0 ). Here, w(x − x 0 , y − y 0 ) is a Gaussian centered at (x 0 , y 0 ) with user-specified width, but with its tail truncated outside the local region of interest. The dominant orientation u ! ¼ cosθ; sinθ ð Þat (x 0 , y 0 ) is a vector of norm 1 which maximises i.e. the norm of the directional derivative of f(x 0 , y 0 ) in the direction of u ! . It is the normalised eigenvector of the largest eigenvalue of the structure tensor matrix J ij (x 0 , y 0 ). The value of the orientation θ is then colour coded according to the scale shown on the right margin of Fig. 1b.
To obtain the average-orientation images, each 1600 × 1200-pixel image was subdivided into a 50 × 38 grid, each subdivision being 32 × 32 pixels (139.2 × 139.2 μm 2 ) in size. We label the subdivisions with a pair of integers (k, ℓ), where 1 ≤ k ≤ l x = 50, and 1 ≤ ℓ ≤ l y = 38. The index k labels the columns of the grid from left to right, while ℓ labels the rows of the grid from top to bottom. The orientation θ kℓ of the subdivision (k, ℓ) was then obtained by setting a Gaussian window of the size of 32 × 32 pixels in OrientationJ.

Correlation functions of average orientation
The correlation functions of average orientation are computed as follows. The correlation between region (i, j) and region (k, ℓ) is defined as C i, j, k, ℓ = cos 2(θ ij − θ kℓ ). The total correlation between region (i, j) and all other regions at distance d from region (i, j) is calculated as follows: The correlation function C(d) is computed as the average of this total correlation over all regions: The distance value d 50 is defined as C(d 50 ) = 0.5. Note that the distance d in these expressions is given in units of 32 pixels (139.2 μm) so it must be multiplied by 139.2 μm to convert to physical units.

Cell movies
To obtain high-resolution images, we used the 35-mm diameter glass-bottom dishes (Iwaki Co., Ltd., Tokyo, Japan) along with a high numerical aperture objective lens. The smaller size of the glass dishes allowed us to place up to three dishes simultaneously inside a temperatureand humidity-controlled microscope (BZ-X700, KEYENCE, Osaka, Japan), enabling continued observation of multiple cell cultures incubating under identical conditions. The glass-bottom culture dishes were collagen-coated following the procedure described above from the 10.0 μg/mL collagen type-I solution. Uncoated controls were prepared using only the vehicle solution (1 mM HCl) in the same procedure. NB1RGB cells (6.0 × 10 4 cells/ 35-mm diameter dish) were seeded in the glass-bottom dishes and cultured for 90 h while taking time-laps images at 15-min intervals using a microscope (BZ-X700, KEYENCE, Osaka, Japan) (see Supplementary Videos 1 and 2).
The VW-H2MA motion analyser (KEYENCE), which performs cell tracking, was used to measure the cell migration velocity. The displacement of each isolated cell during each 15-min interval was measured, from which the average velocity of each cell during that time-interval (which is essentially the instantaneous velocity due to the cells moving slowly) was determined (Fig. 2a). This velocity was averaged over 20 cells. The instantaneous velocities for low-and high-density conditions were respectively based on data from 0 to 30 h and 60 to 90 h after culture start. The 10-h-average velocity, shown in Fig. S3D, is based on the linear displacement during each 10 h interval. This velocity was averaged over 30 cells. The 10-h average velocities for low-and high-density conditions were respectively based on data from 0 to 30 h and 60 to 90 h after culture start (Fig. S3E). The area S and perimeter L of the NB1RGB cells were measured using ImageJ after culturing for 24 h, and circularity was calculated as 4πS/L 2 (Fig. 2e, f).

Mathematical model
To obtain insight into the role of collagen in the 2D patterns formed by fibroblasts, we introduced a simple mathematical model of cell collective motion. The model equation is given as: where the variables (x i (t); y i (t)) and θ i (t) are the position and the orientation of cell i at time t, respectively. In this model, cell i migrates spontaneously in the direction θ i (t) with constant speed v. The terms R x and R y denote the repulsive force due to cell-cell excluded-volume interactions described as a Gaussian soft-core potential H r ð Þ ¼ σ 2 K C exp − r 2 2σ 2 with interaction strength K c , distance r between cells, and repulsion length σ. Explicitly, the terms R x and R y are given as: where r ij stands for the distance between cell i and cell j. Equation (3) describes the nematic interaction with strength K between the cells that leads to the nematic alignment of cell orientations. A similar interaction was considered by Sumino Y. et al. [7]. The effective interaction strength in Eq. (3) is Kexp − r 2 ij 2λ 2 with characteristic interaction length λ. An additive noise is also introduced in Eq. (3), where μ is the noise strength and ξ i (t) is white noise with zero mean and unit variance, i.e., 〈ξ i (t)ξ j (s)〉 = δ ij δ(t − s).

Fixed parameters
The interaction lengths σ and λ are both fixed to σ = λ = 0.696 units (70 μm) in all simulations. We take 140 μm to be the typical cell length, cf. Figure 2c, d, and have chosen σ and λ to be ½ of this value. Except when the excluded volume effect is turned off for Fig. 4, case (vi) (K C = 0), the soft-core interaction strength is taken to be K C = 0.89/h. This value was found via trial and error to reproduce the experimental characteristics well.

Simulations
The continuous-time model of Eqs. (1)-(3) was discretised for each cell using Euler's method with a step size of 30 min, for a simulation length of 144 h (6 days), resulting in 144 × 2 + 1 = 289 frames.
To match the experimental data, each frame's dimensions were adjusted to 69.6 × 52.2 units, where 1 unit is 100 μm. Instead of considering periodic boundaries, which would not match the experimental conditions, we added 15-unit wide margins to the simulation frames. Therefore, the actual frame dimensions for the simulation were 99.6 × 82.2 units, within which cells were initially placed. When creating the simulation video and while taking quantitative measurements for comparison with the experimental results, these margins were ignored.
The initial number of cells is N 0 = 3434, which corresponds to the initial cell density in the experiments,~8000 cells/cm 2 . All the simulated cells were initialised with random positions (x i (0), y i (0)) in the window and with a random orientation θ i (0) in the range [0, 2π).
Based on Fig. 2b, we assume that the increase in the cell number induced by cell divisions starts at t = 24 h and is set to 1.005 cells/0.5 h so that approximately 11,000 cells result after 144 h of cultivation. Each new cell is divided from a randomly chosen existing cell, and positioned 0.35 units (35 μm) away from the parent cell, either in front or behind on the parent's line of motion with equal probability, and with the same orientation and direction of motion. The distance of 35 μm (=σ/2) was chosen so that the position of the offspring will initially lie within the length of the parent, cf. Fig. 2c, d. We do not take into account cell cycles; i.e. each cell proliferates irrespective of its proliferation history. We expect that the detail of the proliferation rule does not considerably affect our results since the cell number is large enough for random and cyclic proliferation rules to be statistically almost equivalent.
GIF animations from the 289 frames were generated using Gnuplot. The block-average orientation θ kℓ of the cells in region (k, ℓ) is defined as the exponent appearing in the equation below: where i is the square root of − 1, and N kℓ is the number of cells in region (k, ℓ). Correlation data were then measured at t = 6, 24, 48, 72, 96, 120, and 144 h, following the protocol described above for the experimental observations. To correct for the effects of randomness in the simulation, ten simulations with the same parameters were run and averaged data were reported as the simulation results.

Dependence of coherence length on migration speed
The simulations were conducted for the following three cases with regard to the constant speed v: 1) The measured migration velocity at low density on an uncoated dish: v = 0.1864 unit/h (1 unit =100 μm, cf. Fig. 2a), 2) The measured migration velocity at low density on a coated dish: v = 0.1598 unit/h (cf. Fig. 2a), 3) The average of the above two: v = 0.1731 unit/h. The dependence of the coherence length on these selections is shown in Fig. 3g.
The solid, dashed, and dotted lines are the average d 50 values over ten simulations for cases (1), (2), and (3), respectively. The error bars are the standard deviation of ten runs for case (1).

Statistical analyses
The data were analysed with the one-tailed Student's t test. The values were expressed as the mean ± mean standard error. Changes were considered to be significant if the p value from the Student's t test was less than 0.05.