Travelling waves of invasion in microbial communities with phenotypic switching

Enunciado do Problema
Comunidades microbianas frequentemente exibem estruturas espaciais complexas onde populações bacterianas clonais se dividem em fenótipos diferentes. Embora a importância ecológica dessa estrutura de subpopulação seja reconhecida, particularmente no que tange à coexistência de espécies em ambientes bem misturados, seu efeito sobre a dinâmica da competição espacial permanece menos compreendido. Especificamente, este artigo investiga como a mudança fenotípica — tanto estocástica quanto responsiva a competidores — afeta as ondas viajantes de invasão que ocorrem quando uma espécie bacteriana invade outra em um habitat espacialmente estendido. O estudo foca em um modelo mínimo envolvendo duas espécies competidoras: Espécie A (um único fenótipo) e Espécie B, que pode alternar entre um fenótipo em crescimento (B) e um fenótipo persistente (P) que é resiliente à competição.

Metodologia
Os autores desenvolvem um modelo mínimo de reação-difusão para descrever a competição espacial entre as duas espécies. O modelo consiste em três equações diferenciais parciais acopladas que rastreiam as abundâncias do fenótipo em crescimento ($B$), do fenótipo persistente ($P$) e do competidor ($A$). As dinâmicas incluem:

• Termos de crescimento logístico e competição de Lotka-Volterra.
• Taxas de mudança estocástica entre $B$ e $P$.
• Mudança responsiva de $B$ para $P$ induzida pela presença do competidor $A$.
• Termos difusivos representando a motilidade bacteriana.

A análise prossegue por meio de duas abordagens complementares:

1. Derivação Analítica: Os autores transformam o sistema em um conjunto de equações diferenciais ordinárias usando uma coordenada de onda viajante $z = x - ct$. Eles analisam as dinâmicas linearizadas próximas aos estados estacionários homogêneos (onde apenas uma espécie existe) para derivar condições algébricas necessárias para a existência de conexões heteroclínicas (ondas viajantes). Isso envolve calcular os autovalores das matrizes jacobianas nas bordas de liderança e de cauda da onda para determinar a velocidade mínima da onda ($c_{min}$) requerida para soluções não negativas.
2. Simulação Numérica: O sistema completo de reação-difusão e um modelo "médio" correspondente (que colapsa as subpopulações fenotípicas em uma única população efetiva) são integrados numericamente no MATLAB. Os autores amostram aleatoriamente espaços de parâmetros para comparar as velocidades de onda reais ($c^*$) selecionadas pelas dinâmicas contra as velocidades mínimas derivadas analiticamente e contra as velocidades previstas pelo modelo médio.

Principais Contribuições e Resultados
O artigo produz duas descobertas primárias, um tanto contra-intuitivas, sobre o efeito da mudança fenotípica na velocidade de invasão:

1. Invasão por Competidores ($c > 0$): Quando a espécie competidora (A) invade a espécie com mudança (B), a presença de mudança fenotípica não tem efeito sobre a velocidade da onda viajante. A velocidade mínima da onda para o modelo completo com mudança é idêntica à do modelo médio sem mudança. Este resultado vale independentemente da taxa de mudança responsiva ($\beta$).

2. Invasão de Competidores ($c < 0$): Quando a espécie com mudança (B) invade o competidor (A), a mudança fenotípica pode aumentar a velocidade da onda de invasão. Especificamente, a velocidade da onda no modelo completo pode exceder a velocidade prevista pelo modelo médio, desde que a taxa de mudança responsiva não seja muito alta. Isso sugere que a capacidade de alternar para um fenótipo resiliente pode atuar como um mecanismo ofensivo, acelerando a tomada do território de um competidor.

Os resultados analíticos estabelecem que a velocidade mínima da onda é determinada pelas dinâmicas linearizadas na borda de liderança da onda. Para $c > 0$, os termos de mudança não aparecem na expressão de ordem principal para a velocidade mínima. Para $c < 0$, os termos de mudança modificam os parâmetros efetivos de crescimento e difusão, permitindo propagação mais rápida sob condições específicas.

Significado e Afirmações
Os autores concluem que a persistência bacteriana, frequentemente vista como uma estratégia defensiva para sobreviver ao estresse ou antibióticos, pode funcionar como uma estratégia ecológica ofensiva na competição espacial. Ao alternar para um fenótipo persistente, uma população pode paradoxalmente acelerar sua invasão de um competidor. Esta descoberta desafia a expectativa intuitiva de que um "pico" de um fenótipo resiliente retardaria uma onda invasora ao atuar como um amortecedor.

O artigo situa este resultado no contexto mais amplo da variação fenotípica estabilizando a coexistência em populações bem misturadas, sugerindo que o mesmo mecanismo (alternar para um fenótipo raro) também pode impulsionar a expansão espacial. Os autores observam que, embora seu modelo seja unidimensional, ele é relevante para o movimento bacteriano em meios porosos como o solo. Eles identificam o mecanismo de seleção de velocidade da onda em sistemas com mudança fenotípica como uma questão matemática aberta para pesquisas futuras, particularmente no que diz respeito a como a imigração proveniente da mudança afeta a seleção de velocidades de onda em modelos clássicos de dinâmica populacional onde $\dot{X}=0$ se $X=0$.