Travelling waves of invasion in microbial communities with phenotypic switching

Enunciato del Problema
Le comunità microbiche spesso esibiscono strutture spaziali complesse in cui popolazioni batteriche clonali si dividono in diversi fenotipi. Sebbene l'importanza ecologica di questa struttura di sottopopolazioni sia riconosciuta, in particolare per quanto riguarda la coesistenza di specie in ambienti ben miscelati, il suo effetto sulla dinamica della competizione spaziale rimane meno compreso. Nello specifico, questo articolo indaga come il passaggio fenotipico, sia stocastico che responsivo ai competitori, influenzi le onde di invasione che si verificano quando una specie batterica invade un'altra in un habitat spazialmente esteso. Lo studio si concentra su un modello minimale che coinvolge due specie in competizione: la Specie A (un singolo fenotipo) e la Specie B, che può passare da un fenotipo in crescita (B) a un fenotipo persistente (P) resiliente alla competizione.

Metodologia
Gli autori sviluppano un modello minimale di reazione-diffusione per descrivere la competizione spaziale tra le due specie. Il modello consiste in tre equazioni differenziali alle derivate parziali accoppiate che tracciano le abbondanze del fenotipo in crescita ($B$), del fenotipo persistente ($P$) e del competitore ($A$). La dinamica include:

• Termini di crescita logistica e di competizione di Lotka-Volterra.
• Tassi di passaggio stocastico tra $B$ e $P$.
• Passaggio responsivo da $B$ a $P$ indotto dalla presenza del competitore $A$.
• Termini diffusivi che rappresentano la motilità batterica.

L'analisi procede attraverso due approcci complementari:

1. Derivazione Analitica: Gli autori trasformano il sistema in un insieme di equazioni differenziali ordinarie utilizzando una coordinata d'onda viaggiante $z = x - ct$. Analizzano la dinamica linearizzata vicino agli stati stazionari omogenei (dove esiste solo una specie) per derivare le condizioni algebriche necessarie per l'esistenza di connessioni eterocliniche (onde viaggianti). Ciò comporta il calcolo degli autovalori delle matrici Jacobiane ai bordi anteriore e posteriore dell'onda per determinare la velocità minima dell'onda ($c_{min}$) richiesta per soluzioni non negative.
2. Simulazione Numerica: Il sistema completo di reazione-diffusione e un corrispondente modello "mediato" (che collassa le sottopopolazioni fenotipiche in un'unica popolazione efficace) sono integrati numericamente in MATLAB. Gli autori campionano casualmente gli spazi dei parametri per confrontare le velocità d'onda effettive ($c^*$) selezionate dalla dinamica con le velocità minime derivate analiticamente e con le velocità previste dal modello mediato.

Contributi e Risultati Chiave
L'articolo produce due principali risultati, in qualche modo controintuitivi, riguardanti l'effetto del passaggio fenotipico sulla velocità di invasione:

1. Invasione da parte dei Competitori ($c > 0$): Quando la specie competitiva (A) invade la specie con passaggio fenotipico (B), la presenza del passaggio fenotipico non ha alcun effetto sulla velocità dell'onda viaggiante. La velocità minima dell'onda per il modello completo con passaggio è identica a quella del modello mediato senza passaggio. Questo risultato vale indipendentemente dal tasso di passaggio responsivo ($\beta$).

2. Invasione dei Competitori ($c < 0$): Quando la specie con passaggio fenotipico (B) invade il competitore (A), il passaggio fenotipico può aumentare la velocità dell'onda di invasione. Nello specifico, la velocità dell'onda nel modello completo può superare la velocità prevista dal modello mediato, a condizione che il tasso di passaggio responsivo non sia troppo elevato. Ciò suggerisce che la capacità di passare a un fenotipo resiliente può agire come un meccanismo offensivo, accelerando la conquista del territorio di un competitore.

I risultati analitici stabiliscono che la velocità minima dell'onda è determinata dalla dinamica linearizzata al bordo anteriore dell'onda. Per $c > 0$, i termini di passaggio non compaiono nell'espressione di ordine principale per la velocità minima. Per $c < 0$, i termini di passaggio modificano i parametri di crescita e diffusione efficaci, permettendo una propagazione più rapida in condizioni specifiche.

Significato e Affermazioni
Gli autori concludono che la persistenza batterica, spesso considerata una strategia difensiva per sopravvivere a stress o antibiotici, può funzionare come una strategia ecologica offensiva nella competizione spaziale. Passando a un fenotipo persistente, una popolazione può paradossalmente accelerare la propria invasione di un competitore. Questa scoperta sfida l'aspettativa intuitiva secondo cui un "rigonfiamento" di un fenotipo resiliente rallenterebbe un'onda invasiva agendo come un cuscinetto.

L'articolo colloca questo risultato nel contesto più ampio della variazione fenotipica che stabilizza la coesistenza in popolazioni ben miscelate, suggerendo che lo stesso meccanismo (passaggio a un fenotipo raro) può anche guidare l'espansione spaziale. Gli autori notano che, sebbene il loro modello sia unidimensionale, è rilevante per il movimento batterico in mezzi porosi come il suolo. Identificano il meccanismo di selezione della velocità dell'onda in sistemi con passaggio fenotipico come una questione matematica aperta per la ricerca futura, in particolare per quanto riguarda come l'immigrazione dal passaggio influisca sulla selezione delle velocità d'onda nei modelli classici di dinamica di popolazione dove $\dot{X}=0$ se $X=0$.