{"id":"c4d804df-ffc1-4ed2-a1a7-6475dbf6dd15","arxiv_id":"2501.11177","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In simulations of bidisperse growing bacterial colonies, cells with larger aspect ratio spontaneously sort to the colony periphery through mechanical interactions alone.","lead":"This paper simulates growing colonies of rod-shaped bacteria of two different shapes, and finds that longer cells spontaneously end up at the colony's outer edge. The result suggests that purely physical pushing and packing, not genes or nutrients, could create an evolutionary advantage for elongated cell shapes in real bacterial colonies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Periphery metric may be biased by rod length: a cell counts as peripheral if any part of its rod lies within 5d of the boundary, so longer rods are more likely to be counted even absent sorting; SI Fig. 8's thickness dependence matches this null geometric effect.","rationale":"The central claim is that long cells progressively enrich the periphery in equal-division-time colonies. The weakest link is not the dynamics but the observable: the periphery is defined by whether any part of the rod lies within a fixed distance of the alpha complex (Methods, Fig. 3). Because a rod's reach toward the boundary scales with its length, the counting rule itself can generate an apparent overrepresentation of long cells in a spatially uniform mixture. The SI's thickness scan (Fig. 8) shows the effect weakens with larger wp, exactly the behavior expected from this geometric bias. The reader flagged the initialization bias as the main weakness; the present concern is broader because it applies even to unbiased initial configurations and to the time-series. A center-based or label-shuffle control would settle this within the authors' existing computational framework. Since such a control is absent and the conclusion is an evolutionary claim, the verdict remains CONDITIONAL pending this test.","tokens_in":15040,"tokens_out":5519,"duration_ms":63658,"concrete_test":"Recompute the periphery fraction ϕA and ϕB using only the cell center (or a fixed point on the rod, e.g., the midpoint) within wp of the alpha complex, for the same saved trajectories used in Figs. 3 and 8. If the positive slope in ϕB vs aB and the time-increasing enrichment of long cells persist under center-based membership, the sorting is real. If the signal vanishes or reverses, the reported effect is dominated by the length-dependent rod-body criterion. As a second, independent null check, take the final configurations and randomly relabel cell identities (or randomly assign lengths to the existing positions/orientations, preserving the length distribution) and recompute the same periphery statistics over many label shuffles; the null distribution should be compared directly with Fig. 3B.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the definition of \"periphery\" in Methods: a cell belongs to the periphery if any point on its rod lies within wp=5d of the alpha complex (Fig. 3 and Methods). For a rod of length L and random orientation, the distance from the center to the boundary at which the rod still intersects a strip of width w is w + (L/2)|cosθ|. Therefore the probability that a uniformly placed cell is counted is an increasing function of L. In the equal-division-time protocol, A and B have equal cell numbers but different L, so even a perfectly mixed colony would show a positive correlation between aspect ratio and periphery fraction by this rod-body criterion. The SI's own Fig. 8 shows the slope of ϕB vs aB decreases as wp grows from 1d to 20d; that is precisely the signature of this geometric cross-section effect (the L-dependent term becomes relatively less important as w grows), rather than evidence for a real boundary-localized population. The paper's observation that sorting develops over time lessens but does not eliminate the concern, since the boundary layer of a growing colony changes thickness and orientational order over time, and the initial condition is explicitly acknowledged to place longer cells in a radially biased and short cells centrally localized configuration (Methods). No control simulation or null calculation is provided that separates length-dependent detection probability from genuine spatial sorting.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses a two-dimensional Brownian dynamics model of growing, dividing circo-rectangle \"cells\" to study colonies containing two populations with different heritable aspect ratios. Under an equal-division-time protocol, the authors report that the longer-aspect-ratio population becomes progressively overrepresented at the colony periphery, and they propose two mechanisms: radial nematic microdomains that drive longer cells outward, and tangential anchoring plus active nematic mixing that retain them at the boundary. They also report that bidispersity suppresses local genetic mixing when one population is below the nematic-ordering threshold, and that under an equal-elongation-rate protocol longer cells are still overrepresented at the periphery relative to their abundance in the bulk. The manuscript includes statistical averaging over 5-10 independent runs, standard errors, robustness checks on the periphery-thickness parameter, lineage tracking, and tracer-particle control simulations.","tokens_in":15307,"tokens_out":4891,"duration_ms":52253,"significance":"If the central sorting claim holds, the result is interesting and potentially important: it would show that purely mechanical interactions in a growing colony can sort microbial populations by cell shape, with possible evolutionary consequences under nutrient-limited conditions where the periphery is the growth zone. The paper's strengths are its careful statistical reporting (multiple independent runs, standard errors), the use of two different growth-rate protocols, the explicit robustness scan over periphery thickness in SI Fig. 8, and the lineage-history analysis that attempts to connect sorting to radial order. However, the central claim currently rests on a periphery definition that is biased by rod length: longer rods are more likely to be counted as peripheral even in a perfectly mixed colony. Because this confound directly affects the paper's main quantitative evidence, the result is not yet established to the standard required for publication.","major_comments":[{"comment":"The periphery definition counts a cell if any of five points along its rod lies within wp = 5d of the alpha complex. For a rod of length L at angle θ to the boundary normal, the center can be up to (L/2)|cos θ| farther from the boundary and still be counted, so longer rods have a higher detection probability in a perfectly mixed colony. The positive slope of ϕB versus aB in Fig. 3B and the decrease of that slope with increasing wp in SI Fig. 8 are exactly the signatures of this length-dependent cross-section effect. The caption of Fig. 3A says the periphery comprises cells whose centers lie within 5d, which conflicts with the Methods description of using five points along each rod; if centers were used the bias would largely disappear, but the reported analysis appears to be body-based. Please provide a null-model correction (e.g., random label shuffling or a center-based periphery metric) or otherwise demonstrate that the sorting signal survives removal of this geometric detection bias.","section":"Methods, 'Periphery Definition'; Fig. 3; SI Fig. 8"},{"comment":"The initialization protocol is acknowledged to introduce a radial orientation bias and to place shorter cells slightly more centrally on average. Because the paper's key claim is that the periphery enrichment emerges over time, the analysis must show that this initial bias does not seed or amplify the observed effect. No control simulation with an alternative initial condition, nor a quantitative bound on the initial bias, is provided. The abstract's statement that the findings are robust across initial conditions is not supported by the presented experiments, which vary the growth-rate protocol but not the initialization protocol.","section":"Methods, 'Simulation Initialization and Stopping Conditions'; Fig. 3A"},{"comment":"The equal-elongation-rate comparison uses the same body-based periphery metric as the equal-division-time analysis. In this protocol the longer cells are also the rarer population, so the overrepresentation measure ϕA/NA − ϕB/NB is inflated by the same length-dependent detection probability: a small number of long rods contributes a disproportionately large periphery count. A center-based or otherwise length-unbiased periphery measure should be applied before concluding that longer cells are overrepresented at the periphery relative to the bulk in this protocol.","section":"Fig. 7D-E (equal elongation rate)"}],"minor_comments":[{"comment":"The Fig. 3A caption states that periphery cells are those whose centers lie within 5d of the alpha complex, while Methods states that five equally spaced points along each rod are used. Please reconcile this discrepancy, as it is directly relevant to the geometric-bias concern.","section":"Fig. 3A caption vs. Methods"},{"comment":"The paragraph after Fig. 4 says 'Replacing the circle of passive tracers with a circle of actively growing Population B in this scenario further confirms...', but Fig. 4A-B show passive tracers and Fig. 4C-D show actively growing rings; the cross-references to panels are confusing and should be corrected.","section":"Fig. 4 and surrounding text"},{"comment":"The Introduction cites Refs. 40-41 for the Brazil nut effect, but those references are about phase transitions in mixtures of rods and spheres, not granular convection. The actual Brazil nut references (48-49) are cited later; please fix the citation placement.","section":"References 40-41"},{"comment":"The term 'circo-rectangles' is used without a definition or standard reference; please define it at first use or use a more common term such as spherocylinders or rounded rectangles.","section":"Model description"},{"comment":"The text states that at early times the more radially aligned cell type is that with the higher aspect ratio, and at late times larger-aspect-ratio cells have higher tangential alignment; it would help to state explicitly that the crossover time is not a fitted parameter and to quantify its uncertainty if possible.","section":"Fig. 5D discussion"}],"recommendation":"major_revision","confidential_remarks":"The geometric confound in the periphery definition is serious enough that I cannot recommend acceptance in the current form. However, the underlying phenomenon may well be real; the authors have the tools to resolve this with a null-model calculation or a center-based periphery metric, plus an initialization control. If those controls confirm the sorting, the paper would be a solid contribution to the active-matter and microbial-evolution literature. The scope fits the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short take: the bidisperse active nematic colony is a genuinely new setup and the periphery enrichment of long rods is visually striking, but the length-dependent periphery metric means the headline effect may be partly a detection artifact. The paper deserves review, but the authors need to add a null control.\n\nWhat's new: previous work (their own Ref. 37) only treated monodisperse colonies; here they mix two aspect ratios and find suppressed mixing when one type is short, and spatial sorting of longer cells to the edge. They also link this to radial nematic microdomains and tangential anchoring, and check robustness across growth protocols and periphery thickness definitions. That's a solid simulation study with clear hypotheses.\n\nSoft spots: the biggest is the periphery definition. A cell is counted as peripheral if any point on its rod lies within 5d of the alpha complex. For a rod of length L, the probability of intersecting a boundary strip of width w scales like w + (L/2)|cosθ|, so longer rods are more likely to be counted even with no spatial sorting. The SI's Fig. 8 shows the slope of φB vs aB decreases as w grows from 1d to 20d — exactly what this geometric cross-section effect predicts. The paper interprets this as evidence the effect is peripheral, but it's also the signature of the detection bias. The initial condition is acknowledged to bias longer cells radially and shorter cells centrally, and no control simulation separates length-dependent detection from genuine sorting. A simple fix: recompute periphery fractions using cell centers only, or run a null simulation with identical dynamics but swapped labels. Until then, the quantitative claim is not clean.\n\nMinor: the proposed mechanisms (radial microdomains, tangential anchoring fence) are plausible but speculative; the granular convection analogy is not tested. No code released, which limits checking.\n\nVerdict: this is a useful paper for the soft matter / microbial evolution crowd, and the authors are thinking carefully. But the load-bearing metric needs a control. I'd send it to review with a request for that null test.","headline":"Periphery enrichment of long rods is new and plausible, but the length-dependent periphery metric means the headline effect may be partly a detection artifact.","tokens_in":15858,"tokens_out":1539,"would_cite":false,"duration_ms":15580,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In growing colonies of rod-shaped bacteria, cells with longer aspect ratio increasingly dominate the periphery through purely mechanical interactions, even with equal division times.","keywords":["bacterial colonies","cell shape","aspect ratio","spatial sorting","active nematics","nematic order","genetic demixing","Brownian dynamics"],"falsifier":"A decisive test would be to repeat the growth protocol from an initially well-mixed but radially unbiased, or deliberately tangentially biased, dense droplet and check whether the longer cells still accumulate at the periphery; if the enrichment disappears or reverses, the effect is an artifact of the initial radial bias rather than an emergent mechanical sorting.","tokens_in":14827,"feed_emoji":"🦠","tokens_out":4877,"duration_ms":44639,"temperature":0.7,"pith_summary":"This paper argues that in a growing two-dimensional colony of rod-shaped bacteria that differ only in how elongated they are, the longer cells end up concentrated at the colony's expanding edge even when both types divide at the same rate. Using overdamped Brownian dynamics simulations of lengthening and dividing circo-rectangles, it shows that this sorting emerges over time and is driven entirely by mechanical interaction. The proposed mechanisms are that longer cells are pushed outward inside radially aligned nematic domains, and once at the boundary they anchor tangentially and tend to stay, while shorter cells get convected back into the bulk. If true, this gives a purely physical, nutrient-independent route by which natural selection could favor elongated cell shapes whenever the periphery is the advantageous zone.","feed_headline":"Longer cells take over a growing colony's rim","feed_subtitle":"A purely mechanical edge: longer cells take the periphery in growing bacterial colonies, hinting at selection for elongation.","key_machinery":"The central object is the bidisperse colony of rod-like agents modeled as lengthening, dividing circo-rectangles (rectangles with semicircular caps) of fixed width and two heritable maximum aspect ratios. Dynamics are overdamped Brownian motion with Hertzian repulsion between overlapping rods, so all ordering arises from steric forces during growth. The quantitative measures that carry the argument are the periphery fraction (cells within 5 cell widths of an $\\alpha$-complex boundary), the local radial alignment parameter $\\bar{s}_r = \\frac{1}{N}\\sum_i \\cos(2(\\theta_i - \\phi_i))$ that distinguishes radially versus tangentially oriented cells, and local heterozygosity for mixing. Together they link nematic order in radial microdomains to transport of long cells outward and tangential anchoring to their retention at the edge.","core_discovery":"The central discovery is a mechanical shape-sorting effect: in a bidisperse colony with equal division times and no growth-rate or nutrient differences, the fraction of the periphery occupied by the longer-aspect-ratio population rises steadily over time, while shorter cells cluster into monoallelic pockets in the bulk. The enrichment at the boundary scales linearly with aspect-ratio difference and survives changes in the width used to define the periphery. The paper traces the effect to two mechanisms: radial expansion pressure organizes longer cells into radial nematic microdomains that help them break through to the boundary, and tangential active anchoring at the periphery makes long cells harder to pull back into the bulk, while active mixing flows carry shorter cells inward. Even under an equal-elongation-rate protocol where short cells dominate by number, long cells remain overrepresented at the periphery relative to their abundance.","pith_inferences":["A direct experimental test would be a microfluidic 2D colony of two isogenic strains engineered to have different lengths but identical growth rates; tracking edge composition over roughly nine doublings should reproduce the linear periphery-fraction slope if the mechanism is purely mechanical.","The same sorting logic may apply to mixtures differing in stiffness or adhesion rather than aspect ratio, since the mechanism only requires shape-dependent nematic order and boundary anchoring.","Because the initial radial bias is a possible confound, the authors' claim would be sharpened by showing that a tangentially biased or fully random well-mixed initial state still produces long-cell enrichment; absent that control, the early-time sorting signal could partly stem from initialization."],"forward_implications":["Periphery composition is a function of cell shape alone: with equal division times, the longer cell type's periphery fraction rises over time and the effect grows linearly with the aspect-ratio gap.","Bidispersity suppresses intermixing: when one cell type is below the nematic-ordering threshold, the colony develops large monoallelic clusters and lower local heterozygosity than a monodisperse colony of the same average shape.","A growing colony acts as an active analogue of granular convection: longer rods rise to the periphery as shorter rods flow inward, analogous to the Brazil nut effect.","Under nutrient-limited conditions, this mechanical sorting implies a selective pressure favoring high aspect ratio, because the periphery is where growth would occur.","Even when shorter cells divide faster, longer cells are still overrepresented at the edge relative to their share of the population."],"supporting_citations":[{"why":"Provides the prior model of active nematic self-mixing dynamics that this work extends, and supplies the baseline for how mixing maintains local genetic diversity.","marker":"[37]"},{"why":"Supplies the framework for nematic microdomains in growing bacterial colonies, which the paper invokes for the radial transport of longer cells.","marker":"[32]"},{"why":"Establishes defect-mediated morphologies and active anchoring at colony boundaries, used here to explain tangential retention of longer cells.","marker":"[30]"},{"why":"Originates the mechanical simulation approach for quasi-two-dimensional microbial colony growth used in this paper.","marker":"[5]"},{"why":"Shows that mechanical interactions in bacterial colonies affect surfing probability of beneficial mutations, motivating the evolutionary relevance of the sorting effect.","marker":"[7]"},{"why":"Provides prior evidence that collective alignment gives longer cells a competitive advantage at the colony periphery during range expansion.","marker":"[16]"},{"why":"Experimental evidence that cell shape affects colony growth under physical confinement, supporting the idea that aspect ratio can be a fitness-relevant trait.","marker":"[24]"},{"why":"Classifies a growing bacterial colony as an active nematic, grounding the interpretation of the observed mixing and buckling dynamics.","marker":"[29]"}],"fun_headline_variants":["Cell shape sorts bacteria: longer cells migrate outward","In mixed colonies, elongated cells win the rim","Bacterial colony demixes by shape alone: long cells edge out","Longer rod-shaped bacteria naturally drift to colony edge","Shape sorting in bacteria: long cells take the periphery"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The main load-bearing premise is that the sorting is emergent and not inherited from how the colony is started: the initialization protocol packs a small isotropic droplet, which already gives a slight radial orientation bias and puts shorter cells slightly closer to the center, and the paper does not test a control with the opposite or zero initial bias.","fun_headline_variants_meta":{"raw":{"variants":["Cell shape sorts bacteria: longer cells migrate outward","In mixed colonies, elongated cells win the rim","Bacterial colony demixes by shape alone: long cells edge out","Longer rod-shaped bacteria naturally drift to colony edge","Shape sorting in bacteria: long cells take the periphery"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000187,"raw_usage":{"total_tokens":1309,"prompt_tokens":908,"completion_tokens":401,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":324}},"tokens_in":524,"tokens_out":401,"duration_ms":3972,"temperature":1.0,"reasoning_tokens":324,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T18:33:12.691208+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be to repeat the growth protocol from an initially well-mixed but radially unbiased, or deliberately tangentially biased, dense droplet and check whether the longer cells still accumulate at the periphery; if the enrichment disappears or reverses, the effect is an artifact of the initial radial bias rather than an emergent mechanical sorting.","supporting_citations":[],"review_version":1}