Travelling waves of invasion in microbial communities with phenotypic switching

Problemstellung
Mikrobielle Gemeinschaften weisen häufig komplexe räumliche Strukturen auf, bei denen klonale Bakterienpopulationen in verschiedene Phänotypen aufspalten. Obwohl die ökologische Bedeutung dieser Subpopulationstruktur anerkannt ist, insbesondere im Hinblick auf das Koexistieren von Arten in gut durchmischten Umgebungen, bleibt ihr Einfluss auf die Dynamik der räumlichen Konkurrenz weniger verstanden. Insbesondere untersucht dieser Artikel, wie phänotypisches Switching – sowohl stochastisch als auch als Reaktion auf Konkurrenten – die Wanderwellen der Invasion beeinflusst, die auftreten, wenn eine Bakterienart in einer räumlich ausgedehnten Lebensraum eine andere Art invadiert. Die Studie konzentriert sich auf ein Minimodell, das zwei konkurrierende Arten umfasst: Art A (ein einzelner Phänotyp) und Art B, die zwischen einem wachsenden Phänotyp (B) und einem persistierenden Phänotyp (P) wechseln kann, der gegenüber Konkurrenz resistent ist.

Methodik
Die Autoren entwickeln ein Minimodell der Reaktions-Diffusion, um die räumliche Konkurrenz zwischen den beiden Arten zu beschreiben. Das Modell besteht aus drei gekoppelten partiellen Differentialgleichungen, die die Häufigkeiten des wachsenden Phänotyps ($B$), des persistierenden Phänotyps ($P$) und des Konkurrenten ($A$) verfolgen. Die Dynamik umfasst:

• Logistisches Wachstum und Lotka-Volterra-Konkurrenzterme.
• Stochastische Switching-Raten zwischen $B$ und $P$.
• Responsives Switching von $B$ zu $P$, ausgelöst durch die Anwesenheit des Konkurrenten $A$.
• Diffusionsterme, die die bakterielle Motilität darstellen.

Die Analyse erfolgt über zwei komplementäre Ansätze:

1. Analytische Herleitung: Die Autoren transformieren das System unter Verwendung einer Wanderwellenkoordinate $z = x - ct$ in ein System gewöhnlicher Differentialgleichungen. Sie analysieren die linearisierte Dynamik in der Nähe homogener stationärer Zustände (wo nur eine Art existiert), um notwendige algebraische Bedingungen für das Vorhandensein heterokliner Verbindungen (Wanderwellen) abzuleiten. Dies beinhaltet die Berechnung der Eigenwerte der Jacobi-Matrizen an der Vorder- und Hinterkante der Welle, um die minimale Wellengeschwindigkeit ($c_{min}$) zu bestimmen, die für nicht-negative Lösungen erforderlich ist.
2. Numerische Simulation: Das vollständige Reaktions-Diffusions-System und ein entsprechendes „gemitteltes" Modell (das die phänotypischen Subpopulationen zu einer einzigen effektiven Population zusammenfasst) werden numerisch in MATLAB integriert. Die Autoren probieren Parameterräume zufällig ab, um die tatsächlichen Wellengeschwindigkeiten ($c^*$), die von der Dynamik ausgewählt werden, mit den analytisch abgeleiteten Mindestgeschwindigkeiten und mit den vom gemittelten Modell vorhergesagten Geschwindigkeiten zu vergleichen.

Hauptbeiträge und Ergebnisse
Der Artikel liefert zwei primäre, etwas kontraintuitive Erkenntnisse bezüglich des Einflusses des phänotypischen Switchings auf die Invasionsgeschwindigkeit:

1. Invasion durch Konkurrenten ($c > 0$): Wenn die Konkurrentenart (A) die switchende Art (B) invadiert, hat das Vorhandensein von phänotypischem Switching keinen Einfluss auf die Geschwindigkeit der Wanderwelle. Die minimale Wellengeschwindigkeit für das vollständige Modell mit Switching ist identisch mit der des gemittelten Modells ohne Switching. Dieses Ergebnis gilt unabhängig von der Rate des responsiven Switchings ($\beta$).

2. Invasion von Konkurrenten ($c < 0$): Wenn die switchende Art (B) den Konkurrenten (A) invadiert, kann phänotypisches Switching die Geschwindigkeit der Invasionswelle erhöhen. Insbesondere kann die Wellengeschwindigkeit im vollständigen Modell die vom gemittelten Modell vorhergesagte Geschwindigkeit überschreiten, sofern die Rate des responsiven Switchings nicht zu hoch ist. Dies deutet darauf hin, dass die Fähigkeit, zu einem resistenten Phänotyp zu wechseln, als offensives Mechanismus wirken und die Übernahme des Territoriums eines Konkurrenten beschleunigen kann.

Die analytischen Ergebnisse zeigen, dass die minimale Wellengeschwindigkeit durch die linearisierte Dynamik an der Vorderkante der Welle bestimmt wird. Für $c > 0$ erscheinen die Switching-Terme nicht im führenden Ausdruck für die Mindestgeschwindigkeit. Für $c < 0$ modifizieren die Switching-Terme die effektiven Wachstums- und Diffusionsparameter und ermöglichen unter spezifischen Bedingungen eine schnellere Ausbreitung.

Bedeutung und Behauptungen
Die Autoren schließen daraus, dass bakterielle Persistenz, die oft als defensive Strategie zum Überleben von Stress oder Antibiotika betrachtet wird, in der räumlichen Konkurrenz als offensive ökologische Strategie fungieren kann. Durch den Wechsel zu einem persistierenden Phänotyp kann eine Population paradoxerweise ihre Invasion eines Konkurrenten beschleunigen. Diese Erkenntnis stellt die intuitive Erwartung in Frage, dass ein „Buckel" eines resistenten Phänotyps eine invadierende Welle verlangsamen würde, indem er als Puffer wirkt.

Der Artikel verortet dieses Ergebnis im weiteren Kontext der phänotypischen Variation, die das Koexistieren in gut durchmischten Populationen stabilisiert, und legt nahe, dass derselbe Mechanismus (Switching zu einem seltenen Phänotyp) auch die räumliche Expansion antreiben kann. Die Autoren stellen fest, dass ihr Modell zwar eindimensional ist, aber für die bakterielle Bewegung in porösen Medien wie Boden relevant ist. Sie identifizieren den Mechanismus der Wellengeschwindigkeitsauswahl in Systemen mit phänotypischem Switching als eine offene mathematische Frage für zukünftige Forschung, insbesondere hinsichtlich der Frage, wie die Immigration aus dem Switching die Auswahl der Wellengeschwindigkeiten in klassischen Populationsdynamik-Modellen beeinflusst, in denen $\dot{X}=0$ gilt, wenn $X=0$.