{"id":"ebe55304-7ebd-451a-bfeb-a592d81308fb","arxiv_id":"2504.15673","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A coarse-grained theory shows that surface-tension-dominated epithelial sheets buckle and wrinkle with thickness scalings and thickness-modulation phases that differ fundamentally from solid plates.","lead":"This paper derives a new elasticity theory for epithelial tissue layers, treating cells as liquid droplets with different surface tensions instead of as solid plates. It predicts thickness scalings for buckling and wrinkling and explains a recently observed inversion of tissue thickness relative to the substrate that classical elasticity cannot capture.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Bulk cell elasticity is omitted but not negligible in the parameter range targeted; the h^0 scalings and phase-inversion threshold are proven only in the ξ≪h limit.","rationale":"I read the paper in good faith as a derivation of a surface-dominated limiting theory. The internal coarse-graining is systematic, the analytic scalings match the vertex-model data shown, and the code is available, so I do not see an internal inconsistency to attack. The load-bearing weakness is exactly the one the reader identified and the paper itself flags: surface and bulk elastic energies have the same order of magnitude over much of the estimated biological parameter range. That makes the biological reach of the central scalings conditional. Since the reader already rendered a CONDITIONAL verdict and my concern does not move it, I recommend UNCHANGED. The proposed bulk-extension test is the minimal check that would either retire the concern or force a more explicit scope restriction.","tokens_in":26432,"tokens_out":11678,"duration_ms":115600,"concrete_test":"Introduce a bulk intracellular modulus into the discrete energy Eq. 1, e.g., add W_bulk=(G/2)Σ_i[(u_{xx}^2+u_{yy}^2)+2ν u_{xx}u_{yy}+(1−ν)u_{xy}^2] built from cell-level displacement gradients, or minimally insert (G h^3/12)(ψ̇)^2 and G h (l−h)^2 into the Lagrangian. Re-derive Eqs. 10, 18, and 26 with G h/Γ_l = 0.1, 1, and 10 (roughly ξ/h = 0.1, 1, 10) and re-run the vertex-model fits in Figs. 2E and 4E. If q0 or Δi changes by more than ~20%, or the h-exponent shifts from 0 toward −1/3, the surface-only scalings do not carry to real epithelia; if the shifts are small, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's novel predictions—μ0=π²/L²∼h^0 (after Eq. 10), q0=(2K)^{1/3}∼h^0 (Eq. 18), and the phase-inversion threshold Δi=2^{-1/3}K^{2/3}/(1−4B/3) (Eq. 26)—are obtained from a Lagrangian in which cell interiors carry no elastic modulus; the only cell-level resistance is the fixed-area constraint. The Discussion's own elastocapillary estimate gives ξ=γ/E≈0.1–10 μm against cell sizes of 1–10 μm and concedes that surface forces 'also may not be dominant.' In the ξ∼h regime, an intracellular bulk modulus E enters the bending channel as ∼Eh^3 and the stretching channel as ∼Eh, i.e., at the same order as the surface and substrate terms. The derived h-exponents and the threshold Eq. 26 are therefore properties of a limiting model, not of epithelia generally. The vertex-model agreement in Figs. 2–4 cannot settle this: the discrete model contains the same omission, and its area penalty κA(A_i−1)^2 constrains area change but does not model bulk shear or compression of the cell interior. Unless ξ/h is shown to be small for a specific tissue, the central claim that surface mechanics governs the observed epithelial elasticity remains an assumption, not an established result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript derives a continuum elasticity theory for a two-dimensional epithelial monolayer in which cell mechanics are dominated by surface tensions on apical, basal, and lateral edges while cell interiors are treated as incompressible fluid. The authors coarse-grain a discrete vertex model, include basement-membrane bending energy and stroma bulk energy, and obtain a harmonic deformation Lagrangian. They analyze linear stability and derive (i) a thickness-independent critical buckling force, (ii) a stroma-induced wrinkling wavelength q0 ~ K^{1/3} independent of thickness, (iii) conditions for buckling-to-wrinkling transitions, and (iv) a phase inversion of groove-to-crest thickness modulation. Analytical predictions are compared with numerical minimizations of the discrete model in Figures 1-4.","tokens_in":26763,"tokens_out":16306,"duration_ms":146633,"significance":"If correct, the results establish a qualitatively different mechanical paradigm for epithelial sheets compared with solid-plate elasticity: critical loads and wavelengths scale differently with layer thickness, and anti-phase thickness modulation emerges without invoking system-spanning fibers. The derivation is careful and largely contained in the SI, and the paper includes concrete, falsifiable predictions (Eqs. 12, 16, 18, 26) and reproducible vertex-model software. The main limitation--neglect of bulk cellular elasticity--is acknowledged in the Discussion, but it is important enough that the paper's scope should be stated more prominently.","major_comments":[{"comment":"Equation (16) as printed gives q[∆,B]_0 = 2^{5/4} sqrt(h^2 - 2/h^2) (|∆| - ∆[B]_c)^{1/2}, i.e., the square-root factor appears in the numerator. The SI derivation (Eqs. S59 and S64) gives the same expression with that factor in the denominator. The denominator version is the physical one: q0 goes to zero at the buckling-to-wrinkling threshold, and it is consistent with the vertex-model data in Fig. 2D. As printed, Eq. (16) would make q0 diverge at threshold and would also change the thickness-modulation prediction in Eq. (24). Please correct Eq. (16) and check all downstream uses.","section":"Wrinkling, Eq. (16)"},{"comment":"The paper's headline scalings--µ0 ~ h^0, q0 ~ h^0, and the phase-inversion threshold Eq. (26)--are derived under the explicit assumption that cell interiors have no bulk elastic modulus. The Discussion's elastocapillary estimate gives ξ ≈ 0.1-10 µm against cell dimensions 1-10 µm, so the assumed hierarchy ξ << h is not guaranteed for the tissues cited in the Introduction. In the ξ ~ h crossover, an intracellular bulk modulus enters the effective bending and stretching channels at the same order as the surface and substrate terms, so the derived exponents and threshold are properties of a limiting model rather than of epithelia generally. The agreement with the vertex model in Figs. 2-4 cannot settle this, because the discrete model embeds the same omission. I ask the authors to state this domain of validity explicitly in the abstract and Introduction and to include a brief quantitative discussion of how a finite bulk modulus would enter the effective Lagrangian and at what ξ/h the scalings cross over to plate-like behavior.","section":"Discussion / domain of validity"}],"minor_comments":[{"comment":"There are several typos: 'embryogensis', 'correspondance', 'addresed' (footnote), 'W rinkling' (section heading), 'cen be' (Section 'Phase inversion'), and 'bucklind' (Fig. 4 caption).","section":"Abstract and footnotes"},{"comment":"The phrase 'An inverse of modulation phase' should read 'An inversion of the modulation phase'.","section":"Section 'Phase inversion'"},{"comment":"The dimensionless energy is denoted e in Eq. (27) but W elsewhere; please align the notation.","section":"Materials and Methods, Eq. (27)"},{"comment":"The vertex model uses a finite area penalty κA = 100 rather than a hard incompressibility constraint; the paper should state explicitly that the continuum derivation assumes the hard-constraint limit and note how the finite penalty is expected to affect the comparison.","section":"Vertex model, Eq. (27)"}],"recommendation":"major_revision","confidential_remarks":"This is a solid theory contribution with a detailed SI and a useful set of falsifiable predictions. The main text contains a load-bearing typographical error in Eq. (16) that contradicts the SI and the figure data, and the scope of the surface-tension-only model should be stated more prominently so that the title and abstract are not over-read. With these revisions the paper would be a good fit for the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about arXiv:2504.15673. First, the paper delivers a clean coarse-grained elasticity theory for epithelial sheets where cell mechanics is dominated by interfacial tensions, with basement membrane bending and stroma compression added. The genuinely new results are the supported-tissue scalings: wrinkling wavenumber q0 ~ K^{1/3} independent of layer thickness h, critical strain ϵc ~ K^{2/3} h^{-2}, critical force µc ~ K^{2/3} h^0, a combined buckling threshold, and a phase-inversion formula for groove-to-crest thickness modulation. These are nontrivial outputs of the Euler–Lagrange equations, not built-in assumptions, and they match vertex-model simulations across the parameter ranges in Figures 2–4. The unsupported cases reduce to the authors' earlier PRL, which is properly cited.\n\nSecond, the paper is honest about its main limitation: cell interiors are modeled as incompressible fluid with no bulk elasticity. The Discussion even gives the elastocapillary length ξ ≈ 0.1–10 μm against cell sizes of 1–10 μm and concedes that surface forces \"also may not be dominant.\" The stress-test note is right: the h^0 scalings and the phase-inversion threshold are properties of the ξ ≪ h limit, and the vertex-model agreement cannot settle the matter because the discrete model shares the same omission. The area penalty in the vertex model constrains area change but does not introduce bulk shear or compression resistance.\n\nThat said, I do not think this is a fatal flaw. The authors explicitly frame the theory as a limiting case and propose tests in systems where bulk elasticity is negligible—adherent lipid vesicle sheets and microfluidic droplet sheets. The analytic predictions are sharp and falsifiable, which is more than many tissue-scale models offer. The phase-inversion prediction is genuinely interesting, though the experimental support is qualitative (blebbistatin organoids, cerebellum sections) rather than quantitative. The continuum derivation is careful, with approximations stated and the resulting regime limits acknowledged. The code is on GitHub, which is a plus for reproducibility.\n\nThe weakest spot is exactly what the stress-test identifies: the elastocapillary length is comparable to cell size for many tissues, so the theory's regime of validity may be narrow. The authors could have done more to delineate where ξ ≪ h holds empirically, and a sentence on how bulk elasticity would enter the Lagrangian (e.g., adding Eh^3 and Eh terms) would sharpen the boundary of the claim. But they do not overclaim; the title says \"governed by interfacial surface mechanics,\" which describes the model, and the Discussion flags the bulk elasticity issue.\n\nThe paper is for soft matter and biophysics readers working on tissue mechanics, buckling of thin sheets, and vertex models. It deserves a serious referee—the derivation is substantial, the scaling laws are counterintuitive and testable, and the phase-inversion mechanism is a real conceptual advance. My own verdict: accept with revisions, conditional on the authors adding a clear discussion of the ξ/h regime and how bulk elasticity would alter the predictions. I'd bring it to a reading group.","headline":"A well-derived surface-tension theory of epithelial sheets gives new scaling laws and a phase-inversion condition, but the h^0 predictions are proven only in the surface-dominated limit, and the paper's own elastocapillary estimate suggests that limit may not hold for many real tissues.","tokens_in":27285,"tokens_out":2373,"would_cite":true,"duration_ms":22660,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Surface-tension-dominated epithelia buckle at a critical force independent of layer thickness, and a differential apico-basal tension can invert the groove-to-crest thickness-modulation phase.","keywords":["epithelial mechanics","surface tension","buckling","wrinkling","coarse-graining","basement membrane","thickness modulation","morphogenesis"],"falsifier":"Measure the critical buckling force or wrinkling wavelength of a well-characterized epithelial monolayer, or a synthetic sheet of adhesive droplets, as a function of layer height $h$ at fixed cell area and aspect ratio. The paper predicts $\\mu_c\\sim h^0$ and $q_0\\sim h^0$ in the surface-tension regime; observing $\\mu_c\\sim h^3$ or $q_0\\sim h^{-1}$, the solid-plate scalings, would falsify the central claim. Alternatively, in an unsupported apico-basally symmetric tissue, any nonzero groove-to-crest thickness modulation contradicts the predicted $\\tau=0$.","tokens_in":26201,"feed_emoji":"🧫","tokens_out":8218,"duration_ms":70167,"temperature":0.7,"pith_summary":"Epithelial tissues are usually interpreted through the elasticity of solid plates, but their cells carry stresses mostly on their surfaces. This paper derives a continuum elasticity theory by coarse-graining a discrete model in which cells are incompressible fluid-filled polygons with constant tensions on their apical, basal, and lateral sides, while the layer is supported by a basement membrane and stroma. It claims that because elasticity originates at cell interfaces, the critical force for buckling is independent of layer thickness and the wrinkling wavelength is independent of thickness, in sharp contrast to solid plates; the same mechanism also makes the groove-to-crest thickness modulation flip phase when apical tension exceeds basal tension. If correct, this gives measurable signatures that tell surface-tension mechanics apart from classical plate elasticity in real epithelia.","feed_headline":"Epithelial buckling force is set by length, not thickness","feed_subtitle":"Cell surface tensions, not bulk elasticity, determine when flat epithelia wrinkle and fold.","key_machinery":"The central object is a discrete vertex model of a tissue as a chain of quadrilateral cell cross-sections with fixed area, then a coarse-graining to a one-dimensional Lagrangian density $L=L_T+L_B+L_K+L_{C_1}+L_{C_2}$ whose terms are the cell-surface energy, basement-membrane bending, stroma bulk elasticity, and two constraints (incompressibility and imposed strain). The load-bearing mathematical object is the dispersion relation $\\mu(q)$ (Eq. 10): minimizing $\\mu$ with respect to $q$ decides between buckling and wrinkling, and the coefficients of $\\mu(q)$ encode every scaling law the conclusions rest on. The physical mechanism is that cell volume is fixed, so deformation energy lives only in the interfacial lengths that change as cells change shape, not in bulk strain of the tissue.","core_discovery":"Within a monolayer of such cells, the coarse-grained deformation energy is harmonic, and the elastic instability of a flat supported sheet is governed by a dispersion relation $\\mu(q)$ for the applied compressive force as a function of wavenumber $q$. Minimizing $\\mu(q)$ selects either buckling ($q_0=0$, or $2\\pi/L$ for a sheet of length $L$) or wrinkling ($q_0>0$). For an unsupported apico-basally symmetric epithelium the critical force is $\\mu_0=\\pi^2/L^2\\sim h^0$, while the critical strain is $\\epsilon_0=\\pi^2/(2h^2L^2)\\sim h^{-2}$; with only the stroma supporting the tissue, $q_0^{(K)}=(2K)^{1/3}\\sim h^0$, $\\mu_c^{(K)}=3\\cdot2^{-4/3}K^{2/3}\\sim h^0$, and $\\epsilon_c^{(K)}=3\\cdot2^{-7/3}h^{-2}K^{2/3}\\sim h^{-2}$. The theory also predicts that differential apico-basal tension $\\Delta=(\\Gamma_a-\\Gamma_b)/\\Gamma_l$ alone can wrinkle an unsupported sheet once $|\\Delta|>\\sqrt{2(1+4B)}$, and that the thickness modulation relative to substrate undulations switches from in-phase to anti-phase when $\\Delta$ exceeds $\\Delta_i=2^{-1/3}K^{2/3}/(1-4B/3)$, so grooves become thicker than crests only for apical tension stronger than basal tension.","pith_inferences":["Editorial inference: because the paper's elastocapillary estimate leaves surface stresses at the edge of dominance in living cells, the $h^0$ scalings are best treated as a limiting law that should cross over to solid-plate scalings when cell bulk elasticity is added; the crossover could be mapped in synthetic droplet sheets.","Editorial inference: the same coarse-graining should carry over to sheets of adherent lipid vesicles or immiscible microfluidic droplets with tunable interfacial tensions, where the thickness-independent buckling and phase inversion can be tested without living-cell variability.","Editorial inference: the model's prediction of zero thickness modulation for an unsupported apico-basally symmetric tissue implies that any observed modulation in such a system would be evidence for a missing ingredient such as bulk elasticity or active contractility."],"forward_implications":["In a surface-tension-dominated epithelium, making the layer thicker does not raise the force needed to buckle it; the critical force depends on the sheet's length, not its height.","When a stroma supports the tissue, the wrinkling wavelength is controlled by the stroma's stiffness $K$, not by epithelial thickness, so a thickness-independent wavelength is a signature of surface-tension mechanics.","Apico-basal differential tension alone can destabilize an unsupported flat sheet into a wrinkle pattern, something classical thin-plate theory does not allow.","The phase of groove-to-crest thickness modulation flips from in-phase to anti-phase when apical tension exceeds basal tension by a threshold fixed by substrate stiffnesses $B$ and $K$; observing this flip identifies the sign of the differential tension.","Measured static quantities, such as the wrinkling wavenumber $q$, the modulation amplitude $\\tau$, and tissue length $L$, can be combined to estimate the ratios $K/\\Delta$ and $B/\\Delta$ from cross-section images."],"supporting_citations":[{"why":"Prior vertex-model result for unsupported sheets showing wrinkling from differential surface tension; supplies the $\\Delta_c=\\sqrt2$ limit and the buckling-force result this continuum theory extends.","marker":"[27]"},{"why":"Supplies the assumption that the stroma's effective bulk modulus scales as $Kq$, the wavenumber-dependent substrate response used in the continuum limit.","marker":"[16]"},{"why":"Gives the solid-plate scalings ($\\epsilon_c\\sim h^0$, $\\mu_c\\sim h$, $q_0\\sim h^{-1}$) against which the paper's scaling claims are compared.","marker":"[37]"},{"why":"Establishes the classical supported-plate result that thickness modulation is in phase with substrate undulations, the baseline the phase-inversion prediction overturns.","marker":"[11]"},{"why":"Documents the observed anti-phase thickness modulation and models it with system-spanning fibers; the paper's surface-tension mechanism is offered as an alternative explanation.","marker":"[31]"},{"why":"Provides the underlying surface-tension-based theory of epithelial elasticity from which the present coarse-graining approach is built.","marker":"[34]"},{"why":"Supplies the liquid-filled cell picture (incompressible cell interiors, surface-localized stresses) on which the discrete cell model rests.","marker":"[21]"}],"fun_headline_variants":["Buckling force for epithelia scales with length, not thickness","Surface tensions dictate when flat epithelia wrinkle and fold","Epithelial buckling is length-limited, not thickness-limited","Differential cell tension alone can wrinkle flat epithelia"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The cells' interiors are taken to be incompressible fluids that carry no bulk elastic stress, so all deformation energy comes from cell surfaces; the paper itself notes that the elastocapillary length $\\xi\\sim\\gamma/E$ is comparable to cell size, so bulk elasticity may contribute and would alter the scalings.","fun_headline_variants_meta":{"raw":{"variants":["Buckling force for epithelia scales with length, not thickness","Surface tensions dictate when flat epithelia wrinkle and fold","Epithelial buckling is length-limited, not thickness-limited","Differential cell tension alone can wrinkle flat epithelia"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000719,"raw_usage":{"total_tokens":3296,"prompt_tokens":1078,"completion_tokens":2218,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":694,"completion_tokens_details":{"reasoning_tokens":2151}},"tokens_in":694,"tokens_out":2218,"duration_ms":13998,"temperature":1.0,"reasoning_tokens":2151,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:20:41.043494+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the critical buckling force or wrinkling wavelength of a well-characterized epithelial monolayer, or a synthetic sheet of adhesive droplets, as a function of layer height $h$ at fixed cell area and aspect ratio. The paper predicts $\\mu_c\\sim h^0$ and $q_0\\sim h^0$ in the surface-tension regime; observing $\\mu_c\\sim h^3$ or $q_0\\sim h^{-1}$, the solid-plate scalings, would falsify the central claim. Alternatively, in an unsupported apico-basally symmetric tissue, any nonzero groove-to-crest thickness modulation contradicts the predicted $\\tau=0$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior vertex-model result for unsupported sheets showing wrinkling from differential surface tension; supplies the $\\Delta_c=\\sqrt2$ limit and the buckling-force result this continuum theory extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the assumption that the stroma's effective bulk modulus scales as $Kq$, the wavenumber-dependent substrate response used in the continuum limit."},{"cited_title":"& Ziherl, P","cited_arxiv_id":null,"evidence_quote":"Gives the solid-plate scalings ($\\epsilon_c\\sim h^0$, $\\mu_c\\sim h$, $q_0\\sim h^{-1}$) against which the paper's scaling claims are compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the classical supported-plate result that thickness modulation is in phase with substrate undulations, the baseline the phase-inversion prediction overturns."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the observed anti-phase thickness modulation and models it with system-spanning fibers; the paper's surface-tension mechanism is offered as an alternative explanation."},{"cited_title":"A., Zhang, T., Lawton, A","cited_arxiv_id":null,"evidence_quote":"Provides the underlying surface-tension-based theory of epithelial elasticity from which the present coarse-graining approach is built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the liquid-filled cell picture (incompressible cell interiors, surface-localized stresses) on which the discrete cell model rests."}],"review_version":1}