pith. sign in

arxiv: 2603.28147 · v3 · submitted 2026-03-30 · 🧮 math-ph · math.AP· math.MP

A Damage-Driven Model for Duchenne Muscular Dystrophy: Early-Stage Dynamics and Invasion Thresholds

Gaetana Gambino , Francesco Gargano , Alessandra Rizzo , Vincenzo Sciacca This is my paper

Pith reviewed 2026-05-14 02:21 UTC · model grok-4.3

classification 🧮 math-ph math.APmath.MP
keywords Duchenne muscular dystrophydamage-driven modelreaction-diffusion-chemotaxisinvasion thresholdpulled-front mechanismearly-stage dynamicsTuring instability
0
0 comments X

The pith

In Duchenne muscular dystrophy, early disease progression spreads through invading damage fronts once an effective reproduction threshold is exceeded, rather than through diffusion-driven spatial instabilities.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a spatially extended mathematical model for Duchenne muscular dystrophy that treats the condition as driven by tissue damage recruiting immune responses. Linear analysis around the healthy equilibrium demonstrates that diffusion does not trigger Turing instabilities, ruling out pattern formation from diffusion alone. Disease advance instead occurs via invasion, with explicit conditions derived for when this invasion begins, framed as a damage reproduction threshold. The minimal speed of the advancing pathological fronts is characterized, revealing a pulled-front mechanism. Simulations validate the transition from decay to invasion, offering a framework focused on early-stage spatial spreading from localized damage.

Core claim

In the reaction-diffusion-chemotaxis system modeling interactions between healthy tissue, damaged fibers, immune cells, and inflammatory signals, the healthy equilibrium remains stable against diffusion-induced instabilities. Progression instead requires crossing an explicit invasion threshold interpreted as a damage reproduction number, after which pathological fronts propagate at a minimal speed determined by a pulled-front mechanism.

What carries the argument

The effective damage reproduction threshold that governs the onset of invasion and the pulled-front dynamics that set the minimal propagation speed of damage fronts.

If this is right

  • Localized damage expands spatially through traveling fronts when the threshold is crossed.
  • The minimal speed of these fronts is explicitly characterizable from the model parameters.
  • Spatial heterogeneity in the tissue arises from invasion of damage rather than intrinsic pattern-forming instabilities.
  • The early dynamics are governed by pulled-front mechanisms, implying specific asymptotic behaviors for the spreading.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Targeting the damage reproduction threshold could be a strategy to halt progression before widespread invasion occurs.
  • Similar models might apply to other conditions where damage triggers inflammatory spread.
  • Parameter estimation from clinical data could refine the threshold predictions for individual patients.

Load-bearing premise

The chosen functional forms for immune recruitment triggered by damage and the chemotactic responses are sufficiently accurate that the derived invasion threshold remains valid for realistic biological parameter ranges.

What would settle it

An experimental observation of Turing-like spatial patterns emerging in early Duchenne muscular dystrophy without initial localized damage sites, or a mismatch between measured front propagation speeds and the predicted pulled-front minimal speed.

Figures

Figures reproduced from arXiv: 2603.28147 by Alessandra Rizzo, Francesco Gargano, Gaetana Gambino, Vincenzo Sciacca.

Figure 1
Figure 1. Figure 1: While multiple equilibria may arise in the general case, the previous result shows that, [PITH_FULL_IMAGE:figures/full_fig_p012_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Temporal and spatial dynamics of the system for different values of [PITH_FULL_IMAGE:figures/full_fig_p039_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Left: Comparison between the analytical dispersion relation s(γ) and numerical estimates of the front speed obtained from exponential initial conditions of the form e −γx . The minimal speed s ∗ corresponds to the minimum of the dispersion curve and is selected for γ ≥ γ ∗ . Right: Comparison between analytical and numerical front speeds for Gaussian initial conditions, for different domain sizes L. The ag… view at source ↗
Figure 4
Figure 4. Figure 4: Comparison between the analytical wave speed predicted by the linearized system [PITH_FULL_IMAGE:figures/full_fig_p041_4.png] view at source ↗
read the original abstract

We introduce a spatially extended mathematical model for Duchenne muscular dystrophy based on a damage-driven paradigm, in which immune recruitment is triggered by tissue injury. The model is formulated as a reaction--diffusion--chemotaxis system describing the interaction between healthy tissue, damaged fibers, immune cells and inflammatory signals. We establish the global well-posedness of the system and investigate the early-stage dynamics through linearization around the healthy equilibrium. Our analysis shows that diffusion does not induce Turing instabilities, so that spatial heterogeneity cannot arise from diffusion-driven mechanisms. Instead, disease progression occurs through invasion processes. We derive explicit conditions for the onset of invasion, interpreted as an effective damage reproduction threshold and characterize the minimal propagation speed of pathological fronts, showing that the dynamics is governed by a pulled-front mechanism. Numerical simulations support the analytical results and confirm the transition between decay and invasion. These results provide a mathematical framework for early-stage disease progression and indicate that spatial spreading arise from the expansion of localized damage rather than from intrinsic pattern-forming mechanisms.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The paper introduces a reaction-diffusion-chemotaxis PDE model for early-stage Duchenne muscular dystrophy in which tissue damage recruits immune cells and inflammatory signals. It establishes global well-posedness of the system, linearizes around the healthy equilibrium to show that diffusion does not produce Turing instabilities, derives explicit invasion thresholds interpreted as an effective damage-reproduction number, and characterizes the minimal speed of pathological fronts via traveling-wave analysis, concluding that progression is governed by a pulled-front mechanism. Numerical simulations are presented to illustrate the transition between decay and invasion regimes.

Significance. If the derivations hold, the work supplies a mathematically consistent framework that cleanly separates invasion-driven spatial spreading from diffusion-induced pattern formation in DMD modeling. The explicit threshold and pulled-front speed results, obtained from standard linear marginal-stability analysis, constitute falsifiable predictions that could be tested against experimental propagation data. The absence of hidden parameter fitting in the threshold derivation is a methodological strength.

major comments (2)
  1. §3 (linearization and invasion threshold): the explicit condition for invasion onset is derived under the assumption that the chemotactic and recruitment functions keep the spatially homogeneous healthy state linearly stable; the manuscript should state the precise inequalities on the recruitment parameters that guarantee this stability, because violation would invalidate the threshold interpretation.
  2. §4 (traveling-wave analysis): the claim that the front is pulled is supported by the linear marginal-stability calculation, but the paper does not compare the analytically predicted minimal speed against the numerically observed speed for the full nonlinear system; a quantitative table or plot of this comparison is needed to confirm the pulled-front regime holds beyond the leading edge.
minor comments (3)
  1. The abstract and introduction should include a brief statement of the precise functional forms chosen for immune recruitment and chemotaxis (e.g., linear or saturating) so that readers can immediately assess the biological scope of the derived threshold.
  2. Figure captions for the numerical simulations should report the exact parameter values used and the spatial domain size, allowing direct reproduction of the decay-to-invasion transition.
  3. A short discussion of possible extensions (e.g., stochastic damage or fiber heterogeneity) would clarify the model's limitations without altering the core claims.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and the constructive comments on our manuscript. We address each major comment below and will revise the paper to incorporate the suggested clarifications and additions.

read point-by-point responses

  1. Referee: §3 (linearization and invasion threshold): the explicit condition for invasion onset is derived under the assumption that the chemotactic and recruitment functions keep the spatially homogeneous healthy state linearly stable; the manuscript should state the precise inequalities on the recruitment parameters that guarantee this stability, because violation would invalidate the threshold interpretation.

    Authors: We agree that the domain of validity for the invasion threshold should be stated explicitly. In the revised manuscript we will add the precise inequalities on the recruitment parameters (derived from the requirement that all eigenvalues of the Jacobian at the healthy equilibrium have negative real parts) that guarantee linear stability of the spatially homogeneous healthy state. This will clarify the parameter regime in which the threshold interpretation holds. revision: yes

  2. Referee: §4 (traveling-wave analysis): the claim that the front is pulled is supported by the linear marginal-stability calculation, but the paper does not compare the analytically predicted minimal speed against the numerically observed speed for the full nonlinear system; a quantitative table or plot of this comparison is needed to confirm the pulled-front regime holds beyond the leading edge.

    Authors: We acknowledge that a direct quantitative comparison would strengthen the pulled-front claim. In the revised manuscript we will include a new table (or figure) that reports the analytically predicted minimal speed from linear marginal-stability analysis alongside the front speeds measured from numerical simulations of the full nonlinear system across a range of parameter values, thereby confirming that the pulled-front regime persists beyond the leading edge. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The paper derives absence of Turing instability and explicit invasion thresholds by linearizing the reaction-diffusion-chemotaxis system around the healthy equilibrium and applying standard traveling-wave marginal-stability analysis to the leading edge. These steps follow directly from the PDE structure and equilibrium stability assumption without reducing to fitted parameters, self-citations, or ansatzes that presuppose the target result. The invasion threshold is obtained as a condition on the linearized growth rate, and the pulled-front speed is the minimal speed from the dispersion relation; both are independent mathematical consequences rather than tautological re-statements of inputs. No load-bearing self-citation chains or renamings of known empirical patterns are required for the central claims.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Model rests on standard existence theory for reaction-diffusion-chemotaxis systems and on domain-specific assumptions about damage-triggered immune recruitment; no new entities are postulated.

axioms (2)
  • standard math The system admits a global classical solution for suitable initial data
    Invoked to justify linearization around the healthy equilibrium
  • domain assumption Immune recruitment is triggered solely by tissue injury with no other independent activation pathways
    Central modeling choice stated in the abstract

pith-pipeline@v0.9.0 · 5489 in / 1226 out tokens · 55094 ms · 2026-05-14T02:21:47.835811+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    On the foundations of cancer modelling: Selected topics, speculations, and perspectives.Mathematical models and methods in applied sciences, 18(4):593–646, 2008

    N Bellomo, N Li, and P Maini. On the foundations of cancer modelling: Selected topics, speculations, and perspectives.Mathematical models and methods in applied sciences, 18(4):593–646, 2008

    work page 2008

  2. [2]

    Clinical analysis of 155 patients with Duchenne muscular dystrophy.Chinese Journal of Contemporary Neurology and Neurosurgery, 15(5):380 – 386, 2015

    Qi Bing, Jing Hu, Zhe Zhao, Hong-Rui Shen, and Na Li. Clinical analysis of 155 patients with Duchenne muscular dystrophy.Chinese Journal of Contemporary Neurology and Neurosurgery, 15(5):380 – 386, 2015

    work page 2015

  3. [3]

    Ca˜ nizo, Laurent Desvillettes, and Klemens Fellner

    Jos´ e A. Ca˜ nizo, Laurent Desvillettes, and Klemens Fellner. Improved duality estimates and applications to reaction-diffusion equations.Communications in Partial Differential Equations, 39(6):1185–1204, 2014. 43

    work page 2014

  4. [4]

    K. N. Chueh, C. C. Conley, and J. A. Smoller. Positively invariant regions for systems of nonlinear diffusion equations.Indiana University Mathematics Journal, 26:373–392, 1977

    work page 1977

  5. [5]

    Stability and phase transitions in a mathemat- ical model of Duchenne muscular dystrophy.Journal of Theoretical Biology, 260(2):283 – 289, 2009

    Guido Dell’Acqua and Filippo Castiglione. Stability and phase transitions in a mathemat- ical model of Duchenne muscular dystrophy.Journal of Theoretical Biology, 260(2):283 – 289, 2009

    work page 2009

  6. [6]

    Desvillettes, V

    L. Desvillettes, V. Giunta, J. Morgan, and B.Q. Tang. Global well-posedness and non- linear stability of a chemotaxis system modelling multiple sclerosis.Proceedings of the Royal Society of Edinburgh Section A: Mathematics, 152(4):826–856, 2022

    work page 2022

  7. [7]

    Desvillettes, T

    L. Desvillettes, T. Lepoutre, A. Moussa, and A. Trescases. On the entropic structure of reaction-cross diffusion systems.Communications in Partial Differential Equations, 40(9):1705–1747, 2015

    work page 2015

  8. [8]

    Desvillettes and A

    L. Desvillettes and A. Trescases. New results for triangular reaction cross diffusion system.Journal of Mathematical Analysis and Applications, 430(1):32–59, 2015

    work page 2015

  9. [9]

    Paul Dowling, Dieter Swandulla, and Kay Ohlendieck. Cellular pathogenesis of Duchenne muscular dystrophy: progressive myofibre degeneration, chronic inflammation, reactive myofibrosis and satellite cell dysfunction.European Journal of Translational Myology, 33(4), 2023

    work page 2023

  10. [10]

    Duchenne muscular dystrophy.Nature Reviews Disease Primers, 7(1), 2021

    Dongsheng Duan, Nathalie Goemans, Shin’ichi Takeda, Eugenio Mercuri, and Annemieke Aartsma-Rus. Duchenne muscular dystrophy.Nature Reviews Disease Primers, 7(1), 2021

    work page 2021

  11. [11]

    Abdelbaset Mohamed Elasbali, Waleed Abu Al-Soud, Saleha Anwar, et al. A review on mechanistic insights into structure and function of dystrophin protein in pathophysiology and therapeutic targeting of Duchenne muscular dystrophy.International Journal of Biological Macromolecules, 264:130544, 2024

    work page 2024

  12. [12]

    R. A. Fisher. The wave of advance of advantageous genes.Annals of Eugenics, 7(4):355– 369, 1937

    work page 1937

  13. [13]

    Gargano, M.C

    F. Gargano, M.C. Lombardo, R. Rizzo, M. Sammartino, and V. Sciacca. Cytokine- induced instabilities in a reaction-diffusion-chemotaxis model of multiple sclerosis: Bi- forcation analysis and well-posedness.International Journal of Non-Linear Mechanics, 161:104672, 2024

    work page 2024

  14. [14]

    Birkh¨ auser, Boston, 2010

    Mi-Ho Giga, Yoshikazu Giga, and J¨ urgen Saal.Nonlinear Partial Differential Equa- tions, volume 79 ofProgress in Nonlinear Differential Equations and Their Applications. Birkh¨ auser, Boston, 2010

    work page 2010

  15. [15]

    Gilbarg and N.S

    D. Gilbarg and N.S. Trudinger.Elliptic Partial Differential Equations of Second Order. Classics in Mathematics. Springer, Berlin Heidelberg, 2001. 44

    work page 2001

  16. [16]

    Judith Haslett, Despina Sanoudou, Alvin Kho, Richard Bennett, Steven Greenberg, Isaac Kohane, Alan Beggs, and Louis Kunkel. Gene expression comparison of biopsies from Duchenne muscular dystrophy (dmd) and normal skeletal muscle.Proceedings of the National Academy of Sciences of the United States of America, 99:15000–5, 12 2002

    work page 2002

  17. [17]

    Hillen and K.J

    T. Hillen and K.J. Painter. A user’s guide to pde models for chemotaxis.Journal of Mathematical Biology, 58(1-2):183 – 217, 2009

    work page 2009

  18. [18]

    Houston, Amanda N

    Matthew T. Houston, Amanda N. Cameron, and Juan B. Gutierrez. A review of math- ematical models for muscular dystrophy: A systems biology approach.PLOS Currents Muscular Dystrophy, 2018

    work page 2018

  19. [19]

    Houston and Juan B

    Matthew T. Houston and Juan B. Gutierrez. The frind model: A mathematical model for representing macrophage plasticity in muscular dystrophy pathogenesis.Bulletin of Mathematical Biology, 81(10):3976 – 3997, 2019

    work page 2019

  20. [20]

    Evans, Robert W

    Abdul Salam Jarrah, Filippo Castiglione, Nicholas P. Evans, Robert W. Grange, and Reinhard Laubenbacher. A mathematical model of skeletal muscle disease and immune response in the mdx mouse.BioMed Research International, 2014:871810, 2014

    work page 2014

  21. [21]

    A. N. Kolmogorov, I. G. Petrovskii, and N. S. Piskunov. A study of the diffusion equation with increase in the amount of substance, and its application to a biological problem. In V. M. Tikhomirov, editor,Selected Works of A. N. Kolmogorov, Volume I, pages 248–

    work page

  22. [22]

    Translated by V

    Kluwer Academic Publishers, Dordrecht, 1991. Translated by V. M. Volosov from Bull. Moscow Univ. Math. Mech. 1 (1937), 1–25

    work page 1991

  23. [23]

    Ladyzenskaija, V

    O.A. Ladyzenskaija, V. A. Solonnikov, and N.N. Ural’ceva.Linear and quasi-linear equa- tions of parabolic type, volume 23 ofTranslations of Mathematical Monographs. American Mathematical Society, Providence, RI, 1988

    work page 1988

  24. [24]

    Demyelination patterns in a mathematical model of multiple sclerosis.Journal of Math- ematical Biology, 75:373–417, 2017

    Maria Carmela Lombardo, Rachele Barresi, Francesco Gargano, and Marco Sammartino. Demyelination patterns in a mathematical model of multiple sclerosis.Journal of Math- ematical Biology, 75:373–417, 2017

    work page 2017

  25. [25]

    On interfaces between cell populations with different mobilities.Kinetic and Related Models, 10(1):299 – 311, 2017

    Tommaso Lorenzi, Alexander Lorz, and Benoˆ ıt Perthame. On interfaces between cell populations with different mobilities.Kinetic and Related Models, 10(1):299 – 311, 2017

    work page 2017

  26. [26]

    Boundedness for reaction-diffusion systems with lya- punov functions and intermediate sum conditions.Nonlinearity, 33(7):3105–3133, 2020

    Jeff Morgan and Bao Quoc Tang. Boundedness for reaction-diffusion systems with lya- punov functions and intermediate sum conditions.Nonlinearity, 33(7):3105–3133, 2020

    work page 2020

  27. [27]

    Zhang, Mansi V

    Michael Ogundele, Jesslyn S. Zhang, Mansi V. Goswami, Marissa L. Barbieri, Utkarsh J. Dang, James S. Novak, Eric P. Hoffman, Kanneboyina Nagaraju, and Yetrib Hathout. Validation of chemokine biomarkers in Duchenne muscular dystrophy.Life, 11(8), 2021

    work page 2021

  28. [28]

    Frontiers in Mathematics

    Benoˆ ıt Perthame.Transport Equations in Biology. Frontiers in Mathematics. Birkh¨ auser, Basel, 2006. 45

    work page 2006

  29. [29]

    Global existence in reaction-diffusion systems with control of mass: A survey.Milan Journal of Mathematics, 78(2):417–455, 2010

    Michel Pierre. Global existence in reaction-diffusion systems with control of mass: A survey.Milan Journal of Mathematics, 78(2):417–455, 2010

    work page 2010

  30. [30]

    Porter, Wei Guo, Anita P

    John D. Porter, Wei Guo, Anita P. Merriam, Sangeeta Khanna, Georgiana Cheng, Xiao- hua Zhou, Francisco H. Andrade, Chellah Richmonds, and Henry J. Kaminski. Persistent over-expression of specific cc class chemokines correlates with macrophage and t-cell re- cruitment in mdx skeletal muscle.Neuromuscular Disorders, 13(3):223 – 235, 2003

    work page 2003

  31. [31]

    Porter, Sangeeta Khanna, Henry J

    John D. Porter, Sangeeta Khanna, Henry J. Kaminski, J. Sunil Rao, Anita P. Merriam, Chelliah R. Richmonds, Patrick Leahy, Li Jingjin, Guo Wei, and Francisco H. Andrade. A chronic inflammatory response dominates the skeletal muscle molecular signature in dystrophin-deficient mdx mice.Human Molecular Genetics, 11(3):263 – 272, 2002

    work page 2002

  32. [32]

    Senesac, Alison M

    Claudia R. Senesac, Alison M. Barnard, Donovan J. Lott, Kavya S. Nair, Ann T. Har- rington, Rebecca J. Willcocks, Kirsten L. Zilke, William D. Rooney, Glenn A. Walter, and Krista Vandenborne. Magnetic resonance imaging studies in Duchenne muscular dys- trophy: Linking findings to the physical therapy clinic.Physical Therapy, 100(11):2035 – 2048, 2020

    work page 2035

  33. [33]

    Oxford Series in Ecology and Evolution

    Nanako Shigesada and Kohkichi Kawasaki.Biological Invasions: Theory and Practice. Oxford Series in Ecology and Evolution. Oxford University Press, UK, 1997

    work page 1997

  34. [34]

    J. Simon. Compact sets in the spaceL p(0, T;B).Annali di Matematica Pura e Applicata, 146(1):65–96, 1986

    work page 1986

  35. [35]

    Smoller.Shock Waves and Reaction-Diffusion Equations, volume 258

    J. Smoller.Shock Waves and Reaction-Diffusion Equations, volume 258. Springer-Verlag New York, NY, 2 edition, 1994

    work page 1994

  36. [36]

    Yu Song, Yingkun Guo, Rong Xu, Ziqi Zhou, Ting Xu, Hang Fu, Ke Xu, Xuesheng Li, Xiaoyu Niu, Ying Ren, Yilin Zhang, Yunan Zhang, Huayan Xu, and Xiaotang Cai. Quantitative magnetic resonance imaging of gluteal muscle groups detects stage-specific progression and early-stage damage in Duchenne muscular dystrophy: A 12-month lon- gitudinal study.Muscle & Nerv...

    work page 2026

  37. [37]

    Spencer and James G

    Melissa J. Spencer and James G. Tidball. Do immune cells promote the pathology of dystrophin-deficient myopathies?Neuromuscular Disorders, 11(6-7):556 – 564, 2001

    work page 2001

  38. [38]

    Changes of serum creatine kinase levels in children with Duchenne muscular dystrophy.Chinese Journal of Contemporary Pediatrics, 10(1):35 – 37, 2008

    Shun-Chang Sun, Yun-Sheng Peng, and Jing-Bo He. Changes of serum creatine kinase levels in children with Duchenne muscular dystrophy.Chinese Journal of Contemporary Pediatrics, 10(1):35 – 37, 2008

    work page 2008

  39. [39]

    James G. Tidball. Inflammatory processes in muscle injury and repair.American Journal of Physiology-Regulatory, Integrative and Comparative Physiology, 288(2):R345–R353,

    work page

  40. [40]

    Front propagation into unstable states.Physics Reports, 386(2-6):29 – 222, 2003

    Wim van Saarloos. Front propagation into unstable states.Physics Reports, 386(2-6):29 – 222, 2003. 46

    work page 2003

  41. [41]

    Histopathology of Duchenne muscular dystrophy in cor- relation with changes in proteomic biomarkers.Histology and Histopathology, 37(2):101– 116, 2022

    Margit Zweyer, Hemmen Sabir, Paul Dowling, Stephen Gargan, Sandra Murphy, Dieter Swandulla, and Kay Ohlendieck. Histopathology of Duchenne muscular dystrophy in cor- relation with changes in proteomic biomarkers.Histology and Histopathology, 37(2):101– 116, 2022. 47

    work page 2022