Epithelial cell chirality emerges through the dynamic concentric pattern of actomyosin cytoskeleton

  1. Takaki Yamamoto
  2. Tomoki Ishibashi
  3. Yuko Mimori-Kiyosue
  4. Sylvain Hiver
  5. Naoko Tokushige
  6. Mitsusuke Tarama
  7. Masatoshi Takeichi
  8. Tatsuo Shibata  Is a corresponding author
  1. Laboratory for Physical Biology, RIKEN Center for Biosystems Dynamics Research, Japan
  2. Nonequilibrium Physics of Living Matter RIKEN Hakubi Research Team, RIKEN Center for Biosystems Dynamics Research, Japan
  3. Laboratory for Molecular and Cellular Dynamics, RIKEN Center for Biosystems Dynamics Research, Japan
  4. Laboratory for Cell Adhesion and Tissue Patterning, RIKEN Center for Biosystems Dynamics Research, Japan
  5. Department of Physics, Kyushu University, Japan
10 figures, 7 videos, 1 table and 1 additional file

Figures

Figure 1 with 2 supplements
Chiral nuclear rotation and cytoplasmic flow in singly isolated Caco-2 cells.

(A) Rotating nucleus probed by the rotation of nuclear texture. The endpoints of the red line segments are the positions of tracked landmarks of the nucleus. (B) The cumulative angle of nuclear rotation plotted against time and (C) average angular velocity averaged over the first 10 hr (n=22). Here, positive angle values indicate clockwise rotation. (D) Chiral cytoskeletal structure of F-actin (phalloidin) and microtubule (immunostaining). Scale bar: 20 µm. (E) Schematic diagram of the orientation of actin stress fibers and microtubule.

Figure 1—figure supplement 1
Effect of physical environment on the nuclear rotation.

Angular velocity of nucleus on the different coating applied to the glass substrate: collagen (n = 19), non-coated (n = 21), fibronectin (n = 22) and poly-l-lysine (n = 22). p values were calculated using Mann–Whitney U test (*p < 0.05, **p < 0.01, ***p < 0.001).

Figure 1—figure supplement 2
Actin bundles in the peripheral region of the cell.

In the control cell (DMSO), actin bundles (phalloidin, magenta) in the peripheral region of cell, which are tilted to form chiral pattern, appear to be anchored to vinculin (green), a focal adhesion protein, at their both ends. We call them stress fibers. In the cell treated with SMIFH2, one end of each actin bundle was anchored to vinculin, while the other ends were not anchored. Consequently, the actin bundles extend in the radial direction. We call them radial fibers.

Figure 2 with 14 supplements
Identification of molecular mechanisms responsible for the chiral rotation.

Roles of F-actin, microtubule, Arp2/3, formin-mediated actin polymerization, and Myosin II activity were investigated. Cells were treated with DMSO (0.2%, control), actin polymerization inhibitor latrunculin A (2 µM), actin depolymerization inhibitor Jasplakinolide (10 nM), microtubule inhibitor nocodazole (50 µM), Arp2/3 inhibitor CK666 (200 µM), formin inhibitor SMIFH2 (40 µM), or Myosin II inhibitor blebbistatin (1 µM). (A) Snapshot images from the live image of actin dynamics in cells expressing Lifeact-RFP. Scale bar: 20 µm. (B) The cumulative angle of nuclear rotation averaged over different cells plotted against time for different conditions: DMSO (n=19), nocodazole (n=11), blebbistatin (n=10), CK666 (n=13), and SMIFH2 (n=10). The standard deviation is represented by shaded regions. (C) Angular velocity of cells under different conditions averaged over the first 5 hr of the time-evolution plot in (B). (D) Angular velocity of control and SMIFH2-treated cells averaged over the last 5 hr of the time-evolution plot in (B): DMSO (n=13) and SMIFH2 (n=10). p values were calculated using the Mann–Whitney U test (p<0.05,p<0.01,p<0.001). Here, positive angle values indicate clockwise rotation.

Figure 2—figure supplement 1
Microtubule inhibition by nocodazole.

Typical snapshots of microtubule of a cell 0.5 and 3.5 hr after the addition of nocodazole (50 µM).The scale bar is 20 µm.

Figure 2—figure supplement 2
Effect of formin depletion on the nuclear rotation.

(A) RNA-sequencing (RNA-seq) analysis showing gene expression levels (transcripts per million; TPM) of major mammalian formin family members (DIAPHs and DAAMs) in Caco-2 cells. Circles represent individual samples (circles of the same color indicate replicates), and bars indicate mean values (n = 3). Western blot showing protein levels of DIAPH2 (B) and DAAM1 (C) in Caco-2 cell treated with siRNAs.GAPDH (bottom row) was used as an internal control. (D) Angular velocity of nucleus of Caco-2 cells that were treated with siRNA for DIAPH2 or DAAM1. The rotation of the nucleus was tracked for 5 hr following the start of the live imaging. Negative control (N.C.) (n = 10), DIAPH2 (n = 18), and DAAM1 (n = 17). p values were calculated using Mann–Whitney U test (*p < 0.05, **p < 0.01, ***p < 0.001). Snapshot images from the live image of actin dynamics in cells expressing Lifeact-RFP in DIAPH2 (E, E’) and DAAM1 (F, F’) depleted cells. Scale bar: 20 µm.

Figure 2—figure supplement 2—source data 1

PDF file containing original western blots for Figure 2—figure supplement 2B, C, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/102296/elife-102296-fig2-figsupp2-data1-v2.pdf
Figure 2—figure supplement 2—source data 2

Original files for western blot analysis displayed in Figure 2—figure supplement 2B, C.

https://cdn.elifesciences.org/articles/102296/elife-102296-fig2-figsupp2-data2-v2.zip
Figure 2—figure supplement 3
Effect of Myosin II depletion on the nuclear rotation.

(A) Angular velocity of nucleus of Caco-2 cells that were treated with siRNA for Myosin II A and/or B heavy chains. The rotation of the nucleus was tracked for 5 hr following the start of the live imaging. Negative control (N.C.) (n = 10), Myosin II A (n = 19), Myosin II B (n = 15), and Myosin IIA and B (n = 20). p values were calculated using Mann–Whitney U test (*p < 0.05, **p < 0.01, ***p < 0.001). (B) Western blot showing protein levels of Myosin IIA (top row) and IIB (middle row) heavy chains in Caco-2 cells treated with siRNAs. GAPDH (bottom row) was used as an internal control. Quantification of fold change relative to the negative control (N.C.), normalized to GAPDH protein level, is shown in bar graphs (bottom panel). Images of Myosin IIA (C) or Myosin IIB (D) (immunostaining, green) and actin filaments (phalloidin, magenta) in cells treated with siRNAs.

Figure 2—figure supplement 3—source data 1

PDF file containing original western blots for Figure 2—figure supplement 3B, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/102296/elife-102296-fig2-figsupp3-data1-v2.pdf
Figure 2—figure supplement 3—source data 2

Original files for western blot analysis displayed in Figure 2—figure supplement 3B.

https://cdn.elifesciences.org/articles/102296/elife-102296-fig2-figsupp3-data2-v2.zip
Figure 2—figure supplement 4
Effect of vinculin depletion on the nuclear rotation.

(A) Angular velocity of nucleus of Caco-2 cells that were treated with siRNA for vinculin. The rotation of the nucleus was tracked for 5 hr following the start of the live imaging. Negative control (N.C.) (n = 10) and vinculin (n = 17). p values were calculated using Mann–Whitney U test (*p < 0.05, **p < 0.01, ***p < 0.001). (B) Western blot showing protein levels of vinculin in Caco-2 cell treated with siRNAs. GAPDH (bottom row) was used as an internal control. Quantification of fold change relative to the negative control (N.C.), normalized to GAPDH protein level, is shown in bar graphs (bottom panel). (C) Images of actin filaments (phalloidin, magenta) and vinculin (immunostaining, green) in cells treated with siRNAs. Scale bar: 20 µm.

Figure 2—figure supplement 4—source data 1

PDF file containing original western blots for Figure 2—figure supplement 4B, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/102296/elife-102296-fig2-figsupp4-data1-v2.pdf
Figure 2—figure supplement 4—source data 2

Original files for western blot analysis displayed in Figure 2—figure supplement 4B.

https://cdn.elifesciences.org/articles/102296/elife-102296-fig2-figsupp4-data2-v2.zip
Figure 2—video 1
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with latrunculin A.

Scale bar: 20 µm.

Figure 2—video 2
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with jasplakinolide.

Scale bar: 20 µm.

Figure 2—video 3
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with nocodazole.

Scale bar: 20 µm.

Figure 2—video 4
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with CK666.

Scale bar: 20 µm.

Figure 2—video 5
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with SMIFH2.

Scale bar: 20 µm.

Figure 2—video 6
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with siRNA for DIAPH2.

Scale bar: 20 µm.

Figure 2—video 7
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with siRNA for DAAM1.

Scale bar: 20 µm.

Figure 2—video 8
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with blebbistatin.

Scale bar: 20 µm.

Figure 2—video 9
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with siRNA for Myosins IIA and B.

Scale bar: 20 µm.

Figure 2—video 10
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with siRNA for vinculin.

Scale bar: 20 µm.

Organization of F-actin and Myosin IIA.

(A) Control cells treated with DMSO show a chiral tilted pattern of F-actin and Myosin II visualized by phalloidin and immunofluorescence with an antibody against Myosin IIA, respectively. (B) SMIFH2 (40 μM) treated cells show a concentric pattern of F-actin and Myosin II. (C) The chiral tilted pattern of F-actin and Myosin II is suppressed in cells treated with blebbistatin (1 μM). The bottom panel shows vertical cross-sections. Scale bars: 20 μm (horizontal) and 5 μm (vertical).

Figure 4 with 1 supplement
Expansion microscopy (ExM) imaging of F-actin and Myosin IIA.

Maximum intensity projection (MIP) images of F-actin (A, C) and Myosin IIA (B, D) in DMSO (A, B) and SMIFH2 (C, D) treated cells. The color indicates the height along the z-axis, where the height was measured after the samples were swollen (color bar, right). Magnified views of the white boxes are shown in the right top panels, and corresponding outlines of F-actin are shown in the right bottom panels, where the bold and dotted lines indicate stress fibers and dorsal actomyosin fibers, respectively. The vertical cross-sections (xz) are shown in the bottom panels, where the bold and dotted lines indicate the peripheral and dorsal inner regions, respectively. Scale bars: 20 µm (horizontal) and 10 µm (vertical). (A’, C’) Composite F-actin images of the ventral (red) and dorsal (green) sides. (A’) In the DMSO-treated cell, the thickness of the ventral and dorsal sides is 2.7 and 6.5 µm, respectively. (C’) In the SMIFH2-treated cells, the thickness of the ventral and dorsal sides is 4.2 and 5.7 µm, respectively.

Figure 4—figure supplement 1
Expansion microscopy (ExM) imaging of actin filaments and Myosin II in cells treated with blebbistatin.

Maximum intensity projection (MIP) images of F-actin (A) and Myosin IIA (B) in blebbistatin-treated cells. The color indicates the height along the z-axis, where the height was measured after the sample was swollen (colorbar, the most right). Magnified views of the white boxes are shown in the right top panels, and corresponding outlines of F-actin are shown in the right bottom panels, where the bold and dotted lines indicate thick and thin fibers, respectively. The vertical section (xz) images are shown in the bottom panels, where the bold and dotted lines indicate peripheral and dorsal inner region, respectively. Scale bars: 20 µm (horizontal), 10 µm (vertical).

Dynamics of F-actin in Caco-2 cells live-imaged by lattice light-sheet microscopy (LLSM).

(A) Maximum intensity projection (MIP) image of Caco-2 treated with DMSO. (B) Snapshot images in the green rectangle in (A), obtained from the z-slice at z=0 is defined as the plane closest to the substrate. The green line is the same as the green circle in (A). Arrowheads indicate the position of filaments. (C) Snapshot images in the red rectangle in (A) obtained from the z-slice at z=0.5μm. The red line is the same as the red circle in (A). Arrowheads indicate the position of filaments. (D) Kymograph along the green circle in (A), obtained from a slice at z=0. (E) Kymograph along the yellow line in (A), obtained from the z-slice at z=0. (F) Kymograph along the red circle in (A), obtained from the z-slice at z=0.5μm. inset: schematic diagram of F-actin (black lines) passing through the circle. (G) Kymograph along the yellow line in (A), obtained from the z-slice z=0.5μm. (H) MIP image of Caco-2 treated with SMIFH2 (40μm). (I) Kymograph along the red circle in (H), obtained from the MIP image. (J) Kymograph along the yellow line in (H), obtained from the MIP image. (K) Schematic diagram of F-actin structure in control and SMIFH2 treated cells. Stress fibers (dark red) were immobile, while the dorsal actin fibers formed an ‘actomyosin ring’ (light blue), moved in centripetal and clockwise directions. Scale bar: 10 µm.

Figure 6 with 4 supplements
Angular velocity obtained by the particle image velocimetry (PIV) analysis.

Spatial profile of angular velocity (color code) obtained from the time average of the PIV vector field in a control cell (A: DMSO) or in a cell treated with SMIFH2 (B) superimposed on a snapshot F-actin image. Scale bar: 20 µm. (C) Average angular velocity as a function of the distance from the center. (D) Average angular velocity as a function of a distance scaled by the inner radius of the actomyosin ring of individual cells. Here, positive angular velocity indicates clockwise rotation. Sample averages for two conditions are indicated by the solid lines. Error bars and shaded areas represent standard errors of the means (SEM).

Figure 6—figure supplement 1
Particle image velocimetry (PIV) analysis of cell treated with DMSO and SMIFH2.

PIV analysis of cell treated with DMSO and SMIFH2. Spatial distribution of azimuthal velocity (top), radial velocity (middle), and angular velocity (bottom) for cell treated with DMSO (control) (A–D) and with SMIFH2 (E–H). (I) Angular averaged azimuthal velocity plotted against the distance from the cell center. (J) Angular averaged radial velocity plotted against the distance from the cell center. Scale bar: 20 µm. Positive azimuthal and angular velocities indicate clockwise rotation.

Figure 6—video 1
Particle image velocimetry (PIV) analysis performed at each time step on the cell shown in Video F-actin ring circulates along the dorsal membrane.

Angular, azimuthal, and radial velocities around the cell center are indicated by color map.

Figure 6—video 2
MIP fluorescent images of Caco-2 expressing Lifeact-mEmelard treated with SMIFH2 taken by LSM880, Zeiss.

Fluorescent debris and nucleus are marked by yellow and white lines.

Figure 6—video 3
MIP fluorescent images of Caco-2 expressing Lifeact-mEmelard treated with SMIFH2 taken by lattice light-sheet microscopy (LLSM).

Two fluorescent debris are marked by red and yellow circles. Scale bar: 10 µm.

Figure 7 with 2 supplements
Theoretical model for the chiral cytoplasmic flow.

(A) Actomyosin generates a force dipole and a torque dipole. (B–D) Numerical simulation of Equation 1 assuming that the cell shape is axisymmetric around the cell center. (B) Actomyosin is distributed along the dorsal membrane with a concentric orientation. (C) Azimuthal velocity vφ showing negative values indicating that the flow is generated in a clockwise direction. (D) Velocity in the radial ρ- and z- directions, (vρ, vz), indicated by vectors. Circulating flow is generated in the ρ-z plane. (E) Top: A concentric orientational field on a ring generates active torque in the center direction (magenta arrow). Middle: Active torque (magenta clockwise arrow) generated by a concentric orientational field on a ring. Bottom: A side view of a cell from the outside toward the cell center. The concentration of actomyosin increases in z (gray color), leading to a gradient of active torque (magenta clockwise arrow) in the z direction, resulting in a rotational flow clockwise (black arrows). (F) Flow profile along the dorsal side showing an inward sinistral swirling pattern. (G) Flow profile at the ventral side showing an outward dextral swirling pattern. (H) Angular velocity averaged in the z direction plotted along the radial direction, showing a peak at around ρ/ρa=1. Here, ρa is the leftmost position where S0.8.

Figure 7—figure supplement 1
Cell shape used in the numerical simulation.

Cell shape used in the numerical simulation. Here, cell shape is assumed to be axisymmetric with respect to the cell center. Z0 is the cell height, R0 is the cell radius, and r1, ρ2 and α are the parameters that determine the cell shape.

Figure 7—figure supplement 2
Comparison of radial and azimuthal velocities to determine the model parameters.

(A–H) Radial and azimuthal velocities shown in Figure 6—figure supplement 1. (I) Radial and azimuthal velocities in the numerical simulation. In the seplots, the radial velocity is positive in the centripetal direction.

Figure 8 with 1 supplement
Depletion of dorsal actin and myosin by Rho Activator II stops the rotating motion of the nucleus.

(A) Actin (magenta) and Myosin II (yellow) showing localization with the dorsal marker Ezrin (cyan) in the DMSO-treated cell. (B) SMIFH2-treated cell showing an increase in dorsal actomyosin and a decrease in ventral actomyosin. (C) Rho Activator II (CN03) treated cell showing a decrease in dorsal actomyosin and an increase in ventral actomyosin. Control (DMSO) and SMIFH2-treated cells showed clockwise (CW) rotation, while CN03-treated cells did show rotation. Scale bars: 20 µm (horizontal) and 5 µm (vertical). Ratio of dorsal F-actin (D) and MyoII (E) to the ventral ones. p values were calculated using the Mann–Whitney U test. p<0.05, n.s.: p0.05.

Figure 8—video 1
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP treated with CN03.

Scale bar: 20 µm.

Author response image 1
Author response image 2

Videos

Video 1
Differential interference contrast (DIC) and fluorescent images of Caco-2 expressing Lifeact-RFP.

Scale bar: 20 µm.

Video 2
Differential interference contrast (DIC) images of Caco-2 with beads attached to the dorsal membrane.

Scale bar: 40 µm.

Video 3
Fluorescent images of Caco-2 expressing Lifeact-RFP and EMTB-3XGFP.

Scale bar: 20 µm.

Video 4
Fluorescent images of Caco-2 expressing Lifeact-RFP treated and MRLC-EGFP.

Scale bar: 20 µm.

Video 5
Fluorescent images of Caco-2 expressing Lifeact-mEmelard taken by lattice light-sheet microscopy (LLSM) at three different z-planes (z = 0, 0.5, and 1.0 µm).
Video 6
MIP fluorescent images of Caco-2 expressing Lifeact-mEmelard treated with SMIFH2 taken by lattice light-sheet microscopy (LLSM).

Scale bar: 10 µm.

Author response video 1

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Cell line (Homo sapiens)Caco-2ATCCOzawa et al., 2020
Transfected construct (Homo sapiens)Lifeact-RFPOzawa et al., 2020
Transfected construct (Homo sapiens)EMTB-3XGFPAddgene, Miller and Bement, 2009Plasmid \# 26741
Transfected construct (Homo sapiens)Lifeact-mEmeraldNakamura et al., 2012pLVSIN-EF1a-Lifeact-mEmerald-IRES-pur
Sequence-based reagentsiRNA: Myosin IIAInvitrogenMYH9HSS106871
Sequence-based reagentsiRNA: Myosin IIBInvitrogenMYH10HSS106875
Sequence-based reagentsiRNA: DAAM1InvitrogenDAAM1HSS177085
Sequence-based reagentsiRNA: DIAPH2InvitrogenDIAPH2HSS102773
Sequence-based reagentsiRNA: VinculinInvitrogenVCL s14764
Antibodyanti-Myosin IIA (rabbit monoclonal)Sigma-AldrichM8064IF (1:1000)
WB (1:1000)
Antibodyanti-Myosin IIB (rabbit monoclonal)Cell Signaling Technology8824WB (1:1000)
Antibodyanti-GAPDH (mouse monoclonal)Santa Cruz Biotechnologysc-166574WB (1:1000)
AntibodyHRP-conjugated anti-rabbit IgG (goat monoclonal)InvitrogenT20926WB (1:1000–1:10,000)
AntibodyHRP-conjugated anti-mouse IgG (goat monoclonal)InvitrogenT20912WB (1:1000–1:10,000)
Antibodyanti-Vinculin (mouse monoclonal)SigmaV9131IF (1:200)
WB (1:1000)
Antibodyanti-Ezrin (mouse monoclonal)Abcamab4069IF (1:1000)
Antibodyanti-DAAM1 (mouse monoclonal)Santa Cruz Biotechnologysc-100942WB (1:1000)
Antibodyanti-Diaph2 (mouse monoclonal)Santa Cruz Biotechnologysc-55540WB (1:500)
AntibodyAlexa Fluor 488 anti-rabbit IgG (goat monoclonal)Sigma-Aldrich11034IF (1:1000)
AntibodyAlexa Fluor 488 anti-mouse IgG (goat monoclonal)InvitrogenA11029IF (1:1000)
AntibodyAlexa Fluor 647 anti-mouse IgG (donkey monoclonal)Sigma-AldrichAP192SA6IF (1:1000)
Antibodyanti-Alexa Fluor 488 (rabbit monoclonal)abcamab150077IF (1:500)
AntibodyAlexa Fluor 488-conjugated anti-rabbit IgG (goat monoclonal)Park et al., 2020A11034IF (1:1000)
Chemical compound, drugLatranculin ASigmaL5163-100UG
Chemical compound, drugJasplakinolideToronto Research Chemicals IncJ210700
Chemical compound, drugNocodazoleSigma-AldrichM1404
Chemical compound, drugCK666Sigma-AldrichSML0006-5MG
Chemical compound, drugSMIFH2Wako4401/10
Chemical compound, drugBlebbistatinSigmaB0560-1MG
Chemical compound, drugRho Activator IICytoskeleton, IncCat. #CN03
Chemical compound, drugAlexa Fluor 568 phalloidinInvitrogenA123801:400
Chemical compound, drugAlexa Fluor 488 phalloidinInvitrogenA123791:500
Software, algorithmImageJ FijiSchindelin et al., 2012Version 2.14.0
Software, algorithmTrim GaloreKrueger et al., 2023Version 0.6.10
Software, algorithmHISAT2Kim et al., 2019Version 2.2.1
Software, algorithmSamtoolsDanecek et al., 2021Version 1.9
Software, algorithmStringtieKovaka et al., 2019Version 2.2.1
Software, algorithmPIVlabThielicke and Sonntag, 2021Version 2.56
Software, algorithmMATLABMathworksR2024a
Software, algorithmFreeFEM++Hecht, 2012Version 4.15

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Takaki Yamamoto
  2. Tomoki Ishibashi
  3. Yuko Mimori-Kiyosue
  4. Sylvain Hiver
  5. Naoko Tokushige
  6. Mitsusuke Tarama
  7. Masatoshi Takeichi
  8. Tatsuo Shibata
(2025)
Epithelial cell chirality emerges through the dynamic concentric pattern of actomyosin cytoskeleton
eLife 14:e102296.
https://doi.org/10.7554/eLife.102296