Transition from traveling fronts to diffusion-limited growth in expanding populations

Probleemstelling

De groei van niet-bewegende cellen (zoals in bacteriële kolonies, biofilms en tumoren) wordt vaak gemodelleerd met reactie-diffusievergelijkingen. Traditionele modellen, zoals de Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) vergelijking en zijn generalisaties, voorspellen dat kolonies zich voortplanten als reistende fronten met een constante snelheid. Dit impliceert dat de straal van de kolonie lineair toeneemt met de tijd ($R \propto t$).

Echter, experimentele waarnemingen tonen aan dat veel micro-organismen onder diverse omstandigheden een sublineaire groei vertonen, waarbij de straal toeneemt volgens een machtsfunctie met een exponent dicht bij 0,5 ($R \propto t^{1/2}$). Dit betekent dat het oppervlak van de kolonie lineair toeneemt, niet de straal. Bestaande modellen die rekening houden met voedselbeperking (nutrient limitation) voorspellen nog steeds lineaire groei van de straal, maar falen om deze veelvoorkomende $t^{1/2}$-groei te verklaren. De auteurs stellen dat de oorzaak ligt in een over het hoofd gezien factor: de afhankelijkheid van biomassa-herverdeling (motiliteit) van de voedselconcentratie.

Methodologie

De auteurs introduceren een nieuw analytisch en numeriek model dat de standaard reactie-diffusie framework aanpast om de koppeling tussen groei en motiliteit te beschrijven.

1. Het Model:
   Ze modificeren de bestaande vergelijkingen (Eq. 3) door de effectieve diffusiecoëfficiënt van de biomassa ($b$) afhankelijk te maken van de voedselconcentratie ($n$). De nieuwe vergelijkingen (Eq. 4) luiden:
   $$\frac{\partial b}{\partial t} = D_b \nabla \cdot (b n \nabla b) + \gamma b n$$
   $$\frac{\partial n}{\partial t} = D_n \nabla^2 n - \gamma b n$$
   Hierbij is de term $bn$ in de diffusiecruciaal: het impliceert dat biomassa alleen kan worden herverdeeld (bewegen) als er voedsel aanwezig is om groei te ondersteunen. Zonder voedsel is er geen actieve mechanische stress en dus geen beweging.

2. Niet-dimensionalisering:
   Het model wordt niet-gedimensioneerd om de analyse te vereenvoudigen. De enige resterende parameter is een effectieve biomassa-diffusiviteit $D$, gedefinieerd als $D = D_b n_0^2 / D_n$.

3. Analytische en Numerieke Benadering:

   • Numeriek: Ze voeren simulaties uit in 1D (smalle kanalen) en 2D (Petrischalen) om de uitbreiding van de biomassa te volgen.
   • Analytisch: Ze analyseren twee regimes:
     • Reistende fronten: Oplossingen van de vorm $f(x-vt)$.
     • Diffusie-gelimiteerde groei: Oplossingen waarbij de straal schaal met $t^{1/2}$. Hierbij wordt aangenomen dat de biomassa-profiel een bewegende Heaviside-stapfunctie is, terwijl het voedselprofiel diffuus blijft.

Kernbijdragen

1. Identificatie van een Kruispunt: Het artikel toont aan dat er een kritieke drempelwaarde ($D_c \approx 1$) bestaat voor de biomassa-diffusiviteit. Boven deze drempel gedraagt het systeem zich als een klassiek reistend front; eronder treedt een volledig ander dynamisch regime op.
2. Analytische Oplossing voor Sublineaire Groei: Voor $D < D_c$ leiden de auteurs een analytische oplossing af waarbij de straal van de kolonie schaal met $t^{1/2}$. Ze tonen aan dat het biomassa-profiel zijn vorm behoudt (invariant), terwijl het voedselprofiel zich uitbreidt en uitgeput raakt.
3. Koppeling Motiliteit en Groei: Het paper legt de theoretische basis voor waarom biomassa-herverdeling afhankelijk is van voedsel, gebaseerd op mechanische stress binnen de kolonie. Dit verklaart waarom experimenten met niet-bewegende cellen vaak sublineaire groei vertonen.

Belangrijkste Resultaten

• Twee Groeiregimes:
  • Lineair Regime ($D > D_c$): De straal $R$ groeit lineair met de tijd ($R \propto t$). De snelheid $v$ hangt af van de parameters volgens $v \approx \sqrt{D} - 1/\sqrt{D}$. Dit komt overeen met traditionele modellen.
  • Sublineair Regime ($D < D_c$): De straal groeit volgens een wortelwet ($R \propto \sqrt{t}$). De snelheid van de rand neemt af met de tijd, maar de vorm van het biomassa-profiel blijft stabiel.
• Kritieke Overgang: Bij $D_c \approx 1$ divergeert de effectieve diffusieconstante en verandert het gedrag van lineair naar sublineair.
• Invloed van Dimensionaliteit:
  • In 1D (kanalen) is de groeisnelheid onafhankelijk van de voedseldiffusie ($D_n$) bij lage $D$.
  • In 2D (cirkelvormige kolonies) hangt de snelheid wel af van zowel $D_n$ als $D_b$, maar met een andere schaalverhouding ($\kappa \propto n_0$ in 1D versus $\kappa \propto \sqrt{n_0}$ in 2D, tot op logaritmische correcties).
• Kolomhoogte (H): De hoogte van de kolonie ($H$) is omgekeerd evenredig met de wortel van de diffusiviteit ($H \propto 1/\sqrt{D}$). Bij lage diffusie (langzame uitbreiding) stapelt biomassa zich verticaal op (dikke kolonies), terwijl bij hoge diffusie de kolonies dunner zijn.
• Verklaring voor Experimentele Data: Het model biedt een verklaring voor de waarneming dat bij veel microben het oppervlak lineair toeneemt (wat overeenkomt met $R \propto \sqrt{t}$), terwijl andere studies een lineaire straal rapporteren. Dit hangt af van waar het systeem zich bevindt ten opzichte van de kritieke drempel $D_c$.

Significantie

Dit onderzoek is van groot belang voor de biophysica en ecologie omdat het een fundamentele discrepantie tussen theorie en experiment oplost.

• Theoretisch: Het toont aan dat standaard reactie-diffusie modellen, wanneer ze worden aangepast met een realistische afhankelijkheid van voedsel voor motiliteit, een breed scala aan groeigedragingen kunnen verklaren zonder nieuwe complexe mechanismen te hoeven invoeren.
• Experimenteel: Het biedt een voorspelling dat de groei van kolonies sterk afhankelijk is van de initiële voedselconcentratie en de mechanische eigenschappen van de cellen. Het suggereert dat veel natuurlijke microbiële kolonies zich precies bij de overgang tussen lineaire en sublineaire groei bevinden.
• Toepassingen: De inzichten zijn relevant voor het begrijpen van tumorgroei, biofilmvorming, en de verspreiding van invasieve soorten, waar de dynamiek van de uitbreiding cruciaal is voor het beheer en de behandeling.

Kortom, het paper bewijst dat de overgang van reistende fronten naar diffusie-gelimiteerde groei een natuurlijk gevolg is van de koppeling tussen voedselbeschikbaarheid en de mechanische beweging van cellen binnen een dichte populatie.