{"id":"087e8aa5-bc18-4630-ad9c-60787d0e27c5","arxiv_id":"2510.27262","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Motile and non-motile Pseudomonas putida produce similar total biomass in a porous chip, but reduce permeability by 78% vs 94%, showing biomass spatial organization controls flow resistance more than biomass amount.","lead":"Using a microfluidic chip with pillars, the authors grew motile and non-motile bacteria under constant pressure and found the same amount of biomass blocked flow very differently: 78% permeability loss for swimmers versus 94% for non-swimmers. The study concludes that where biomass sits inside pores, not just how much there is, controls clogging.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim depends on uncalibrated GFP intensity as a proxy for occluding biomass; if WT and ΔfliC differ in fluorescence per unit volume, the 'same total biomass' comparison collapses and the model's quantitative predictions lose their basis.","rationale":"The reader already identified the lack of fluorescence calibration as the weakest assumption, and my reading agrees. The central claim is an empirical comparison of total biomass in two strains, and that comparison is made entirely in uncalibrated GFP intensity units. If the fluorescence-to-occluding-volume ratio differs between strains, the conclusion that 'not total amount' is controlling is not supported. This is a testable missing measurement rather than an incurable flaw: the authors can add an orthogonal biomass/volume calibration and, if needed, adjust the model. The direct flow measurements and triplicate experiments are real evidence, and the model could still work after recalibration, so the appropriate disposition remains conditional on that added validation—i.e., the reader's original CONDITIONAL verdict is appropriate and no verdict adjustment is required. A secondary gap is the placeholder data/code availability statement, which also supports the conditional framing but is not the central concern.","tokens_in":17355,"tokens_out":4457,"duration_ms":69488,"concrete_test":"Repeat the colonization experiments with an orthogonal, volume-sensitive measurement of biofilm: for at least three WT and three ΔfliC chips at t=44 h, acquire confocal z-stacks or optical coherence tomography scans in the same pores used for GFP imaging, and measure actual biofilm thickness/volume; in parallel, quantify total cell biomass (e.g., protein content, dry weight, or cell counts) per chip. Compute strain-specific fluorescence-to-volume conversion factors and re-run Eqs. (1)-(2) with those factors. If the strains have equal true biomass/volume but different permeability, the central claim is confirmed; if their true volumes differ or their fluorescence per unit volume differs, the headline comparison must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline conclusion—spatial organization, not total biomass, controls permeability—rests on the observation that WT and ΔfliC reach nearly identical total biomass (K≈4.6×10^4 vs 4.7×10^4 in arbitrary GFP units) while reducing permeability by 78±7% and 94±4%. However, all biomass quantities are defined through GFP pixel intensities: B_exp is normalized integrated intensity, and the pore-scale biofilm thickness d2i(t) used in Eq. (2) is set equal to ρ_pi^i D_i^MIC/I_max with ρ_pi^i = ΣI_mic/(I_max A_mic) (Methods F, Eq. 8). No calibration relates fluorescence to biofilm volume, thickness, dry mass, or cell number, and no check establishes equal GFP expression per unit biomass across WT and ΔfliC. If the non-motile mutant produces more EPS per cell, has altered GFP expression in the stationary phase, or exhibits different fluorescence quantum yield when densely packed, then the supposedly matched 'total biomass' may conceal a real difference in occluding volume, and the permeability gap could be explained by total biomass after all. This is the load-bearing assumption of the empirical test, not a peripheral modeling detail.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper combines microfluidic porous-media experiments, constant-pressure flow control, and time-lapse fluorescence microscopy to compare biofilm formation by motile wild-type (WT) and non-motile (ΔfliC) Pseudomonas putida. It reports that the two strains reach nearly identical integrated GFP biomass while reducing permeability by different amounts (78±7% for WT vs 94±4% for ΔfliC), and interprets this as evidence that spatial organization of biomass, not total biomass, controls permeability. A pore-scale model is developed in which fluorescence intensity is converted into an effective biofilm thickness within each pore, and the permeability of the medium is computed by treating each pore as two parallel flow pathways. The model, with one fitted parameter k_bf, reproduces the measured permeability curves. The paper also presents idealized nutrient-limited vs space-limited growth scenarios to argue that growth regime interacts with pore-size heterogeneity to determine permeability decline.","tokens_in":17673,"tokens_out":3849,"duration_ms":41931,"significance":"If the central claim is correct, the paper provides a valuable experimental demonstration that biomass spatial distribution can matter more than total biomass for bioclogging, with direct implications for filtration, bioremediation, and subsurface engineering. The experimental setup is a strength: triplicate experiments, direct flow-rate measurements under controlled pressure, spatial maps of biomass, and a pore-scale modeling framework. The authors also state that data and code are available. However, the quantitative conclusions rest on an uncalibrated GFP-to-biomass thickness conversion and on a model whose single free parameter is fitted to the same data it is compared with. These issues must be addressed before the central claim is fully supported, but they are not disqualifying if calibration or independent validation can be supplied.","major_comments":[{"comment":"The conversion from fluorescence intensity to biofilm thickness is uncalibrated and dimensionally suspect. The text defines ρ_pi^i = ΣI_mic^i/(I_max A_mic^i) and then sets d2i(t) = ρ_pi^i(t) D_i^MIC / I_max, which appears to use I_max twice and yields an inconsistent dimension. More importantly, no independent calibration links GFP pixel intensity to biofilm volume, thickness, dry mass, or cell number, and no check establishes that WT and ΔfliC fluoresce identically per unit occluding biomass. Since d2i is the key input to Eq. (2), the model's quantitative agreement in Fig. 5 depends on this unverified proportionality. The authors should either calibrate the fluorescence-to-thickness relationship (e.g., confocal thickness measurements, cell counts, or EPS quantification) or clearly reframe the model as a qualitative demonstration. This issue also affects the global 'same total biomass' c","section":"§II.A and Methods F, Eq. (8)"},{"comment":"The model's agreement with experiments is not a prediction in the current form, because k_bf is estimated as the best fit to the same k_exp curves shown in Fig. 5d,e. The statement that the model 'fully consistent' with measurements is therefore circular at the parameter level. The authors should obtain k_bf independently (e.g., from the cited literature range, from a separate set of experiments, or via cross-validation) and then show the resulting permeability curves. This would not change the core empirical observation, but it would strengthen the claim that the pore-scale biomass distribution quantitatively explains the permeability dynamics.","section":"§II.A, Fig. 5d,e; Eq. (2)"},{"comment":"The central empirical claim—'nearly identical total biomass' causing different permeability reductions—relies entirely on GFP fluorescence as a proxy for biomass amount. If the ΔfliC mutant produces more extracellular polymeric substance per cell, has different GFP expression in stationary phase, or packs more densely in the pore space, then the two strains may not actually have the same occluding biomass. The paper provides no independent biomass measurement (e.g., cell counts, dry weight, confocal volume, or staining of EPS). Without such a control, the possibility remains that the permeability difference is due to a difference in total occluding volume rather than spatial organization. This is a load-bearing assumption for the headline conclusion and should be explicitly tested or discussed as a limitation.","section":"§I.C and §I.D"}],"minor_comments":[{"comment":"The formula for d2i(t) appears to have a typo: d2i(t) = ρ_pi^i(t) D_i^MIC / I_max is dimensionally inconsistent because ρ_pi^i already includes 1/I_max. This should be corrected.","section":"§II.A, text around Eq. (1)"},{"comment":"The mutant is referred to as 'ΔfilC' in one place and 'ΔfliC' elsewhere; the spelling should be consistent (the gene is fliC).","section":"§Methods B"},{"comment":"The Logistic Growth (LG) model is cited to reference [55], which is an ISME paper on multispecies biofilms and not a standard source for the logistic equation. A methods or modeling reference would be more appropriate.","section":"§Methods F"},{"comment":"The contributions list 'N.W. developed the feedback loop', but the author list contains Nolwenn Delouche (N.D.), not N.W. This appears to be a typo.","section":"Author Contributions"},{"comment":"The homogeneous-structure control in Fig. 6d is described in the text but not clearly distinguished in the figure legend; adding an explicit label or color key would improve readability.","section":"§II.B, Fig. 6"}],"recommendation":"major_revision","confidential_remarks":"The paper has a strong experimental core and a clear, interesting message. My main concern is that the quantitative equivalence of biomass between the two strains is established only through GFP intensity, with no calibration to occluding biomass. This is fixable with additional experiments or an explicit, well-argued limitation, so I do not recommend rejection. The modeling claim of 'prediction' should be softened or supported by independent estimation of k_bf. I would be comfortable with acceptance after these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the motile/non-motile comparison is the real contribution, and it is probably right in the qualitative direction. But the paper labels a fit as a prediction and leans on uncalibrated fluorescence as a measure of occluding biomass, so the headline claim needs qualification before it is taken at face value.\n\nWhat is genuinely new: a controlled comparison of WT and ΔfliC under constant pressure drop, with triplicates and direct flow measurements. Both strains reach nearly the same integrated GFP signal (K about 4.6 vs 4.7 in arbitrary units) yet permeability drops 78% and 94% respectively. The spatial maps show the non-motile strain keeps filling downstream pores while the motile one does not. That asymmetry is the interesting result.\n\nWhat the paper does well: the model is a sensible extension of the authors' earlier pore-network permeability framework. Representing each pore as a clean channel in parallel with a biofilm annulus of permeability k_bf, and using measured biomass maps, is reasonable. Getting two very different permeability curves with a single fitted k_bf is not nothing.\n\nSoft spots: the stress-test note lands. All biomass quantities are GFP pixel intensities (Eq. 8). There is no calibration linking intensity to biofilm thickness, volume, cell count, or EPS content, and no check that WT and ΔfliC produce the same fluorescence per unit occluding biomass. If the mutant packs more tightly or makes more EPS, the 'same total biomass' premise collapses. That is load-bearing, not peripheral. Second, k_bf is fit to the same k_exp curves in Fig. 5, so using 'predicts' in the abstract is an overstatement. The model is consistent, not independently predictive. Third, the data availability statement is a placeholder.\n\nFor whom: biofilm and porous-media experimentalists, particularly people working on bioclogging and microfluidic models. It deserves a serious referee: the experimental system is well-built and the comparison is worth publishing. But referees should require either an independent biomass/thickness calibration or substantial rewriting of the claims so they do not rest on uncalibrated GFP and on a fit presented as a prediction.\n\nRecommendation: send it to peer review, not desk reject. Expect revision.","headline":"Worth a referee: the motile/non-motile permeability contrast is a real result, but 'same biomass' is uncalibrated GFP and the model's quantitative match is a fit wearing a prediction's clothes.","tokens_in":18135,"tokens_out":3780,"would_cite":true,"duration_ms":41734,"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":"Spatial biomass layout, not total amount, controls how much biofilms clog porous media.","keywords":["biofilm","porous media","permeability","bioclogging","bacterial motility","microfluidics","pore-scale model","Pseudomonas putida"],"falsifier":"Measure the actual biofilm thickness and volume in the same microfluidic pores using confocal microscopy or optical coherence tomography and compare these to the fluorescence-derived thicknesses used in the model; if fluorescence-to-thickness conversion is off by more than the stated uncertainties, the quantitative permeability predictions would fail for both strains.","tokens_in":17259,"feed_emoji":"🦠","tokens_out":1677,"duration_ms":17379,"temperature":0.7,"pith_summary":"This paper claims that the spatial organization of biomass, not its total quantity, is the primary control on permeability in biofilm-colonized porous media. The evidence comes from microfluidic experiments comparing motile wild-type bacteria with non-motile mutants: both strains reach nearly identical total biomass, yet the motile strain reduces permeability by 78±7% while the non-motile strain reduces it by 94±4%. The difference arises because motile bacteria stop accumulating downstream where nutrients are scarce, whereas non-motile cells continue growing and progressively clog the entire pore space. The authors build a pore-scale mechanistic model that predicts measured permeability dynamics using only one fitted parameter—biofilm permeability—by treating each pore as a central open channel surrounded by a biomass annulus. If correct, this reframes how bio-clogging should be predicted: measurements of total biomass alone are insufficient; where that biomass sits in the pore network matters more.","feed_headline":"Where biomass sits, not how much, controls clogging","feed_subtitle":"Motile and non-motile bacteria reach identical biomass, yet cut permeability 78% vs 94%: placement decides the clog.","key_machinery":"The central object is a two-pathway pore-scale flow model: each pore is idealized as a central biofilm-free pipe of diameter d1i with Hagen-Poiseuille permeability d1i²/32, surrounded by an annular biofilm layer of thickness d2i/2 with an effective biofilm permeability kbf. The equivalent pore permeability is the thickness-weighted average of these two parallel flow paths. Biomass thickness in each pore is inferred from fluorescence intensity converted to a biomass density ρ_pi, which shrinks the effective open-pore diameter. The macroscopic permeability is then computed as the harmonic mean over ~896 serial pore elements, with kbf as the only fitted parameter. This machinery links pore-scal","core_discovery":"The central claim is that the spatial organization of biomass, not its total amount, is the primary factor controlling intrinsic permeability of porous systems. This is supported by experiments where wild-type (motile) and ΔfliC (non-motile) Pseudomonas putida reach nearly identical total biomass carrying capacity (K≈4.6–4.7×10^4 in the logistic growth model), yet cause drastically different permeability reductions: 78±7% for motile versus 94±4% for non-motile. Motility limits downstream biomass accumulation—motile cells escape resource-depleted zones and are advected away—whereas non-motile cells continue slow growth even under resource limitation, filling and clogging the entire system. A","pith_inferences":["A testable extension: the same experimental setup applied to a fixed flow-rate boundary (instead of constant pressure) should show different permeability dynamics because local shear would promote detachment; the model could be adapted with a detachment term.","The framework suggests that engineered bio-barriers could be made more effective by suppressing bacterial motility, while antifouling strategies in filtration should target spatial coverage rather than total biofilm mass.","If fluorescence does not linearly report occluding biomass volume, the quantitative comparison between WT and ΔfliC may be biased; independent thickness measurements would strengthen or correct the central claim.","The contrast between nutrient-limited and space-limited regimes hints that the pore-size distribution of a medium could be used to predict whether biomass will clog it uniformly or leave flow pathways open."],"forward_implications":["Predicting bioclogging in soils, filters, and aquifers requires spatial biomass distribution data, not just total biomass or bulk porosity measurements.","Bacterial motility is a key control: motile strains that escape nutrient-depleted regions leave downstream pores open, preserving higher overall permeability.","The two-regime growth model predicts that space-limited growth clogs pores far more effectively per unit biomass than nutrient-limited growth, and that heterogeneous pore networks clog faster than homogeneous ones at the same biomass.","The single fitted biofilm permeability value (kbf=2.5 darcy) is consistent with independent literature estimates, suggesting the model captures a physically meaningful biofilm property rather than an arbitrary fit.","Under constant-pressure conditions, permeability can continue dropping long after total biomass plateaus, as seen for the non-motile strain, implying post-plateau biomass redistribution is hydraulically significant."],"fun_headline_variants":["Biomass layout, not load, dictates clogging","Motility steers clogging: 78% vs 94% permeability loss","Spatial pattern of bacteria, not total, controls flow","Where microbes settle decides pore clogging"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The model assumes that GFP fluorescence intensity is a quantitative, linearly comparable measure of flow-occluding biomass volume, converting fluorescence into biofilm thickness without independent calibration against actual biomass density or thickness.","fun_headline_variants_meta":{"raw":{"variants":["Biomass layout, not load, dictates clogging","Motility steers clogging: 78% vs 94% permeability loss","Spatial pattern of bacteria, not total, controls flow","Where microbes settle decides pore clogging"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000166,"raw_usage":{"total_tokens":1054,"prompt_tokens":674,"completion_tokens":380,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":418,"completion_tokens_details":{"reasoning_tokens":311}},"tokens_in":418,"tokens_out":380,"duration_ms":33260,"temperature":1.0,"reasoning_tokens":311,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T07:00:03.030440+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the actual biofilm thickness and volume in the same microfluidic pores using confocal microscopy or optical coherence tomography and compare these to the fluorescence-derived thicknesses used in the model; if fluorescence-to-thickness conversion is off by more than the stated uncertainties, the quantitative permeability predictions would fail for both strains.","supporting_citations":[],"review_version":1}