Travelling waves of invasion in microbial communities with phenotypic switching

Enunciado del Problema
Las comunidades microbianas a menudo exhiben estructuras espaciales complejas donde poblaciones bacterianas clonales se dividen en diferentes fenotipos. Si bien se reconoce la importancia ecológica de esta estructura de subpoblaciones, particularmente en lo que respecta a la coexistencia de especies en entornos bien mezclados, su efecto sobre la dinámica de competencia espacial permanece menos comprendido. Específicamente, este artículo investiga cómo el cambio fenotípico —tanto estocástico como responsivo a los competidores— afecta las ondas viajeras de invasión que ocurren cuando una especie bacteriana invade a otra en un hábitat espacialmente extendido. El estudio se centra en un modelo mínimo que involucra dos especies competidoras: la Especie A (un solo fenotipo) y la Especie B, que puede cambiar entre un fenotipo en crecimiento (B) y un fenotipo persistente (P) que es resiliente a la competencia.

Metodología
Los autores desarrollan un modelo mínimo de reacción-difusión para describir la competencia espacial entre las dos especies. El modelo consiste en tres ecuaciones diferenciales parciales acopladas que rastrean las abundancias del fenotipo en crecimiento ($B$), el fenotipo persistente ($P$) y el competidor ($A$). La dinámica incluye:

• Términos de crecimiento logístico y competencia de Lotka-Volterra.
• Tasas de cambio estocástico entre $B$ y $P$.
• Cambio responsivo de $B$ a $P$ inducido por la presencia del competidor $A$.
• Términos difusivos que representan la motilidad bacteriana.

El análisis procede mediante dos enfoques complementarios:

1. Derivación Analítica: Los autores transforman el sistema en un conjunto de ecuaciones diferenciales ordinarias utilizando una coordenada de onda viajera $z = x - ct$. Analizan la dinámica linealizada cerca de los estados estacionarios homogéneos (donde solo existe una especie) para derivar condiciones algebraicas necesarias para la existencia de conexiones heteroclínicas (ondas viajeras). Esto implica calcular los valores propios de las matrices jacobianas en los bordes delantero y trasero de la onda para determinar la velocidad mínima de la onda ($c_{min}$) requerida para soluciones no negativas.
2. Simulación Numérica: El sistema completo de reacción-difusión y un modelo "promediado" correspondiente (que colapsa las subpoblaciones fenotípicas en una sola población efectiva) se integran numéricamente en MATLAB. Los autores muestrean aleatoriamente los espacios de parámetros para comparar las velocidades de onda reales ($c^*$) seleccionadas por la dinámica contra las velocidades mínimas derivadas analíticamente y contra las velocidades predichas por el modelo promediado.

Contribuciones y Resultados Clave
El artículo arroja dos hallazgos principales, algo contraintuitivos, sobre el efecto del cambio fenotípico en la velocidad de invasión:

1. Invasión por Competidores ($c > 0$): Cuando la especie competidora (A) invade a la especie con cambio fenotípico (B), la presencia de cambio fenotípico no tiene efecto sobre la velocidad de la onda viajera. La velocidad mínima de la onda para el modelo completo con cambio es idéntica a la del modelo promediado sin cambio. Este resultado se mantiene independientemente de la tasa de cambio responsivo ($\beta$).

2. Invasión de Competidores ($c < 0$): Cuando la especie con cambio fenotípico (B) invade al competidor (A), el cambio fenotípico puede aumentar la velocidad de la onda de invasión. Específicamente, la velocidad de la onda en el modelo completo puede exceder la velocidad predicha por el modelo promediado, siempre que la tasa de cambio responsivo no sea demasiado alta. Esto sugiere que la capacidad de cambiar a un fenotipo resiliente puede actuar como un mecanismo ofensivo, acelerando la toma del territorio de un competidor.

Los resultados analíticos establecen que la velocidad mínima de la onda está determinada por la dinámica linealizada en el borde delantero de la onda. Para $c > 0$, los términos de cambio no aparecen en la expresión de orden principal para la velocidad mínima. Para $c < 0$, los términos de cambio modifican los parámetros efectivos de crecimiento y difusión, permitiendo una propagación más rápida bajo condiciones específicas.

Significado y Afirmaciones
Los autores concluyen que la persistencia bacteriana, a menudo vista como una estrategia defensiva para sobrevivir al estrés o a los antibióticos, puede funcionar como una estrategia ecológica ofensiva en la competencia espacial. Al cambiar a un fenotipo persistente, una población puede paradójicamente acelerar su invasión de un competidor. Este hallazgo desafía la expectativa intuitiva de que un "bulto" de un fenotipo resiliente ralentizaría una onda invasora al actuar como un amortiguador.

El artículo sitúa este resultado dentro del contexto más amplio de la variación fenotípica estabilizando la coexistencia en poblaciones bien mezcladas, sugiriendo que el mismo mecanismo (cambiar a un fenotipo raro) también puede impulsar la expansión espacial. Los autores señalan que, aunque su modelo es unidimensional, es relevante para el movimiento bacteriano en medios porosos como el suelo. Identifican el mecanismo de selección de velocidad de onda en sistemas con cambio fenotípico como una pregunta matemática abierta para futuras investigaciones, particularmente en lo que respecta a cómo la inmigración desde el cambio afecta la selección de velocidades de onda en modelos clásicos de dinámica poblacional donde $\dot{X}=0$ si $X=0$.