Pulled, pushed or failed: the demographic impact of a gene drive can change the nature of its spatial spread

Resumen Técnico: "Pulled, pushed or failed: the demographic impact of a gene drive can change the nature of its spatial spread"

Planteamiento del problema
El artículo aborda la necesidad crítica de comprender la propagación temporal y espacial de los alelos de la impulsión genética (gene drive) antes de realizar experimentos de campo. Mientras que los modelos tempranos utilizaron marcos de genética de poblaciones clásicos que asumían generaciones discretas y no solapadas en poblaciones bien mezcladas, estas simplificaciones no logran capturar la heterogeneidad espacial y la dinámica demográfica continua de los escenarios del mundo real. Específicamente, la introducción de individuos transgénicos mal adaptados puede reducir el tamaño de la población, creando un flujo demográfico desde áreas densas hacia áreas dispersas que se opone a la propagación del alelo de la impulsión. El problema central es esclarecer cómo las variaciones en la densidad de población (dinámicas demográficas) afectan la velocidad de propagación y la naturaleza de las ondas de impulsión genética a través del espacio, particularmente bajo regímenes de crecimiento extremos.

Metodología
Los autores emplean un enfoque de modelado de reacción-difusión determinista en un contexto espacial unidimensional. El estudio se centra en soluciones de ondas viajeras, que representan perfiles estacionarios que se mueven a una velocidad constante $v$.

1. Formulación del modelo:

   • Los modelos describen la dinámica de tres genotipos: homocigotos silvestres ($WW$), homocigotos de la impulsión ($DD$) y heterocigotos ($DW$).
   • La aptitud (fitness) se asigna como $f_{WW}=1$, $f_{DD}=1-s$ (donde $s$ es el costo de aptitud) y $f_{DW}=1-sh$ (donde $h$ es el coeficiente de dominancia).
   • La conversión génica (homing) se modela con una probabilidad $c \in [0, 1]$, que ocurre ya sea en el zigoto o en la línea germinal.
   • El sistema incorpora un término de crecimiento dependiente de la densidad parametrizado por la tasa de crecimiento intrínseco $r$ y un término de difusión (normalizado a 1).

2. Enfoques analíticos y numéricos:

   • Regímenes extremos: Los autores investigan analíticamente dos casos límite para la tasa de crecimiento intrínseco $r$:
     • $r \to +\infty$: Representa un escenario donde la población restaura instantáneamente la capacidad de carga, minimizando el flujo demográfico.
     • $r = 0$: Representa un escenario donde los nacimientos equilibran exactamente las muertes en una población silvestre, maximizando el impacto demográfico de los costos de aptitud. Este régimen se analiza estableciendo una conexión matemática con los modelos epidemiológicos SI (Susceptibles-Infectados).
   • Conversión parcial: A diferencia de trabajos previos que asumen una conversión del 100%, este estudio modela explícitamente la conversión parcial ($c < 1$) y la presencia resultante de heterocigotos.
   • Simulaciones numéricas: Extensas simulaciones numéricas del problema de Cauchy cierran la brecha entre los regímenes de $r=0$ y $r=+\infty$, generando mapas de calor de las velocidades de invasión a través de los espacios de parámetros ($s, h, c, r$).

Contribuciones clave y resultados

1. Caracterización de los tipos de ondas:
   El artículo distingue entre ondas tiradas (pulled waves) (donde la velocidad está determinada por la dinámica en baja densidad/frecuencia) y ondas empujadas (pushed waves) (donde la velocidad es estrictamente mayor que la velocidad mínima linealizada, impulsada por la masa de la población).

   • En el régimen $r \to +\infty$: Los resultados se alinean con estudios previos y los extienden. Dependiendo del costo de aptitud $s$ y la tasa de conversión $c$, la impulsión puede fallar, invadir como una onda tirada, invadir como una onda empujada o permitir la invasión del tipo silvestre.
   • En el régimen $r = 0$: Un hallazgo significativo es que las ondas no pueden ser empujadas cuando la tasa de crecimiento intrínseco se anula. El sistema se reduce a una forma análoga a los modelos SI epidemiológicos, donde la impulsión o bien invade como una onda tirada (si el costo de aptitud es lo suficientemente bajo) o falla completamente (depuración de la impulsión genética).

2. Impacto del flujo demográfico:
   El estudio demuestra que los efectos demográficos son más pronunciados cuando $r$ es bajo. En este régimen, el flujo demográfico (movimiento de alta a baja densidad) actúa como una barrera para la propagación de los alelos de la impulsión que conllevan un costo. Consecuentemente, una impulsión que logra invadir con éxito en un régimen de alto crecimiento ($r \to +\infty$) puede fallar en invadir o depurarse por completo en un régimen de bajo crecimiento ($r=0$).

3. La conjetura de "Siempre Tirada" (Always Pulled):
   A través de simulaciones numéricas que conectan los dos regímenes extremos, los autores observan que si una onda es tirada tanto en el límite $r=0$ como en el $r=+\infty$, permanece siendo tirada para todos los valores intermedios de $r$. En estos casos, la velocidad de la onda parece ser independiente del parámetro demográfico $r$. Esto ocurre cuando:

   • El costo de aptitud $s$ es suficientemente pequeño.
   • Existe una coexistencia estable de los alelos de la impulsión y el tipo silvestre tras la invasión.

4. Umbrales y estabilidad:
   El artículo proporciona umbrales analíticos para la transición entre la invasión de la impulsión, la invasión del tipo silvestre y la coexistencia. Por ejemplo, en el límite $r=0$ con conversión zigótica, la impulsión invade si $s < s_{2,z} = \frac{c}{2c + h(1-c)}$. Si $s$ excede este umbral, la impulsión se depura.

Significancia y afirmaciones
El artículo afirma delimitar las condiciones específicas bajo las cuales la dinámica demográfica altera fundamentalmente la propagación de las impulsiones genéticas. Al conectar el límite $r=0$ con los modelos epidemiológicos SI, los autores proporcionan herramientas analíticas para predecir el fallo de la impulsión en poblaciones con bajas tasas de crecimiento intrínseco.

El estudio enfatiza que asumir una densidad de población constante (el límite $r \to +\infty$) puede conducir a predicciones excesivamente optimistas sobre el éxito de la impulsión genética. Específicamente, la naturaleza de la onda (tirada vs. empujada) y el éxito de la invasión pueden cambiar drásticamente dependiendo de la capacidad de la población para recuperarse de los costos de aptitud impuestos por la impulsión. Los autores concluyen que la dinámica demográfica no es meramente un factor de fondo, sino que puede determinar si una impulsión genética se propaga con éxito, falla o cambia la composición genética de la población en un entorno espacialmente estructurado.