Travelling waves of invasion in microbial communities with phenotypic switching

Probleemstelling
Microbiële gemeenschappen vertonen vaak complexe ruimtelijke structuren waarbij klonale bacteriële populaties zich splitsen in verschillende fenotypes. Hoewel het ecologische belang van deze subpopulatiestructuur wordt erkend, met name met betrekking tot het samenleven van soorten in goed gemengde omgevingen, blijft het effect daarvan op ruimtelijke competitiedynamiek minder goed begrepen. Specifiek onderzoekt dit artikel hoe fenotypische schakeling – zowel stochastisch als responsief op concurrenten – de reizende golven van invasie beïnvloedt die optreden wanneer één bacteriesoort een andere binnendringt in een ruimtelijk uitgebreide habitat. De studie richt zich op een minimaal model dat twee concurrerende soorten omvat: Soort A (een enkel fenotype) en Soort B, die kan schakelen tussen een groeiend fenotype (B) en een persister-fenotype (P) dat veerkrachtig is tegen competitie.

Methodologie
De auteurs ontwikkelen een minimaal reactie-diffusiemodel om de ruimtelijke competitie tussen de twee soorten te beschrijven. Het model bestaat uit drie gekoppelde partiële differentiaalvergelijkingen die de abundantie volgen van het groeiende fenotype ($B$), het persister-fenotype ($P$) en de concurrent ($A$). De dynamiek omvat:

• Logistische groei en Lotka-Volterra competitie-termen.
• Stochastische schakelsnelheden tussen $B$ en $P$.
• Responsieve schakeling van $B$ naar $P$ veroorzaakt door de aanwezigheid van concurrent $A$.
• Diffusieve termen die bacteriële motiliteit voorstellen.

De analyse verloopt via twee complementaire benaderingen:

1. Analytische Afleiding: De auteurs transformeren het systeem naar een set gewone differentiaalvergelijkingen met behulp van een reizende-golfcoördinaat $z = x - ct$. Ze analyseren de linearisatie van de dynamiek nabij de homogene stationaire toestanden (waar slechts één soort bestaat) om noodzakelijke algebraïsche voorwaarden af te leiden voor het bestaan van heterokliene connecties (reizende golven). Dit omvat het berekenen van de eigenwaarden van de Jacobiaanse matrices aan de voor- en achterrand van de golf om de minimale golfsnelheid ($c_{min}$) te bepalen die vereist is voor niet-negatieve oplossingen.
2. Numerieke Simulatie: Het volledige reactie-diffusiesysteem en een corresponderend "gegemiddeld" model (dat de fenotypische subpopulaties samenvoegt tot één effectieve populatie) worden numeriek geïntegreerd in MATLAB. De auteurs willekeurig parameter ruimtes om de werkelijke golfsnelheden ($c^*$) die door de dynamiek worden geselecteerd te vergelijken met de analytisch afgeleide minimale snelheden en met de snelheden voorspeld door het gegemiddelde model.

Belangrijkste Bijdragen en Resultaten
Het artikel levert twee primaire, enigszins contra-intuïtieve bevindingen op met betrekking tot het effect van fenotypische schakeling op invasiesnelheid:

1. Invasie door Concurrenten ($c > 0$): Wanneer de concurrerende soort (A) de schakelende soort (B) binnendringt, heeft de aanwezigheid van fenotypische schakeling geen effect op de snelheid van de reizende golf. De minimale golfsnelheid voor het volledige model met schakeling is identiek aan die van het gegemiddelde model zonder schakeling. Dit resultaat geldt ongeacht de snelheid van responsieve schakeling ($\beta$).

2. Invasie van Concurrenten ($c < 0$): Wanneer de schakelende soort (B) de concurrent (A) binnendringt, kan fenotypische schakeling de snelheid van de invasiegolf verhogen. Specifiek kan de golfsnelheid in het volledige model de snelheid die door het gegemiddelde model wordt voorspeld, overschrijden, mits de snelheid van responsieve schakeling niet te hoog is. Dit suggereert dat het vermogen om naar een veerkrachtig fenotype te schakelen kan fungeren als een offensief mechanisme, waardoor de overname van het territorium van een concurrent wordt versneld.

De analytische resultaten stellen vast dat de minimale golfsnelheid wordt bepaald door de linearisatie van de dynamiek aan de voorrand van de golf. Voor $c > 0$ komen de schakeltermen niet voor in de uitdrukking van de leidende orde voor de minimale snelheid. Voor $c < 0$ modificeren de schakeltermen de effectieve groeiparameters en diffusieparameters, waardoor snellere voortplanting onder specifieke omstandigheden mogelijk wordt.

Betekenis en Beweringen
De auteurs concluderen dat bacteriële persistentie, vaak gezien als een defensieve strategie om stress of antibiotica te overleven, kan fungeren als een offensieve ecologische strategie in ruimtelijke competitie. Door te schakelen naar een persister-fenotype kan een populatie paradoxaal genoeg de invasie van een concurrent versnellen. Deze bevinding daagt de intuïtieve verwachting uit dat een "bult" van een veerkrachtig fenotype een binnenvallende golf zou vertragen door te fungeren als buffer.

Het artikel plaatst dit resultaat in de bredere context van fenotypische variatie die het samenleven stabiliseert in goed gemengde populaties, en suggereert dat hetzelfde mechanisme (schakelen naar een zeldzaam fenotype) ook ruimtelijke expansie kan aandrijven. De auteurs merken op dat hoewel hun model eendimensionaal is, het relevant is voor bacteriële beweging in poreuze media zoals bodem. Zij identificeren het mechanisme van golfsnelheidsselectie in systemen met fenotypische schakeling als een open wiskundige vraag voor toekomstig onderzoek, met name met betrekking tot hoe immigratie vanuit schakeling de selectie van golfsnelheden beïnvloedt in klassieke populatiedynamische modellen waar $\dot{X}=0$ als $X=0$.