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

Résumé technique : « Tirée, poussée ou échouée : l'impact démographique d'un forçage génétique peut changer la nature de sa propagation spatiale »

Énoncé du problème
L'article aborde la nécessité critique de comprendre la propagation temporelle et spatiale des allèles de forçage génétique avant de mener des expériences sur le terrain. Alors que les premiers modèles utilisaient des cadres de génétique des populations classiques supposant des générations discrètes et non chevauchantes dans des populations bien mélangées, ces simplifications ne parviennent pas à capturer l'hétérogénéité spatiale et la dynamique démographique continue des scénarios du monde réel. Plus précisément, l'introduction d'individus transgéniques mal adaptés peut réduire la taille de la population, créant un flux démographique des zones denses vers les zones clairsemées qui s'oppose à la propagation de l'allèle de forçage. Le problème central est de clarifier comment les variations de la densité de population (dynamique démographique) affectent la vitesse de propagation et la nature des ondes de forçage génétique à travers l'espace, particulièrement sous des régimes de croissance extrêmes.

Méthodologie
Les auteurs emploient une approche de modélisation par réaction-diffusion déterministe dans un contexte spatial unidimensionnel. L'étude se concentre sur les solutions d'ondes progressives, qui représentent des profils stationnaires se déplaçant à une vitesse constante $v$.

1. Formulation du modèle :

   • Les modèles décrivent la dynamique de trois génotypes : les homozygotes de type sauvage ($WW$), les homozygotes de forçage ($DD$) et les hétérozygotes ($DW$).
   • Le fitness est assigné comme suit : $f_{WW}=1$, $f_{DD}=1-s$ (où $s$ est le coût de fitness) et $f_{DW}=1-sh$ (où $h$ est le coefficient de dominance).
   • La conversion génique (homing) est modélisée avec une probabilité $c \in [0, 1]$, se produisant soit dans le zygote, soit dans la lignée germinale.
   • Le système incorpore un terme de croissance dépendant de la densité, paramétré par le taux de croissance intrinsèque $r$, et un terme de diffusion (normalisé à 1).

2. Approches analytiques et numériques :

   • Régimes extrêmes : Les auteurs étudient analytiquement deux cas limites pour le taux de croissance intrinsèque $r$ :
     • $r \to +\infty$ : Représente un scénario où la population restaure instantanément la capacité de charge, minimisant le flux démographique.
     • $r = 0$ : Représente un scénario où les naissances équilibrent exactement les décès dans une population de type sauvage, maximisant l'impact démographique des coûts de fitness. Ce régime est analysé en établissant un lien mathématique avec les modèles épidémiologiques SI (Susceptible-Infecté).
   • Conversion partielle : Contra�à contrairement aux travaux précédents supposant une conversion de 100 %, cette étude modélise explicitement la conversion partielle ($c < 1$) et la présence résultante d'hétérozygotes.
   • Simulations numériques : Des simulations numériques étendues du problème de Cauchy comblent l'écart entre les régimes $r=0$ et $r=+\infty$, générant des cartes thermiques des vitesses d'invasion à travers les espaces de paramètres ($s, h, c, r$).

Contributions clés et résultats

1. Caractérisation des types d'ondes :
   Le papier distingue entre les ondes tirées (où la vitesse est déterminée par la dynamique à faible densité/fréquence) et les ondes poussées (où la vitesse est strictement supérieure à la vitesse minimale linéarisée, pilotée par le cœur de la population).

   • Dans le régime $r \to +\infty$ : Les résultats s'alignent sur et étendent les études précédentes. Selon le coût de fitness $s$ et le taux de conversion $c$, le forçage peut échouer, envahir en tant qu'onde tirée, envahir en tant qu'onde poussée, ou mener à une invasion du type sauvage.
   • Dans le régime $r = 0$ : Une découverte importante est que les ondes ne peuvent pas être poussées lorsque le taux de croissance intrinsèque est nul. Le système se réduit à une forme analogue aux modèles épidémiologiques SI, où le forçage soit envahit en tant qu'onde tirée (si le coût de fitness est suffisamment bas), soit échoue complètement (élimination du forçage génique).

2. Impact du flux démographique :
   L'étude démontre que les effets démographiques sont les plus prononcés lorsque $r$ est faible. Dans ce régime, le flux démographique (mouvement des zones de haute vers les zones de basse densité) agit comme une barrière à la propagation des allèles de forçage coûteux. Par conséquent, un forçage qui réussit l'invasion dans un régime de croissance élevée ($r \to +\infty$) peut échouer à envahir ou être totalement éliminé dans un régime de faible croissance ($r=0$).

3. La conjecture du « Toujours Tiré » :
   À travers des simulations numériques reliant les deux régimes extrêmes, les auteurs observent que si une onde est tirée dans les deux limites $r=0$ et $r=+\infty$, elle reste tirée pour toutes les valeurs intermédiaires de $r$. Dans ces cas, la vitesse de l'onde semble indépendante du paramètre démographique $r$. Cela se produit lorsque :

   • Le coût de fitness $s$ est suffisamment faible.
   • Il existe une coexistence stable entre l'allèle de forçage et l'allèle de type sauvage après l'invasion.

4. Seuils et stabilité :
   L'article fournit des seuils analytiques pour la transition entre l'invasion par le forçage, l'invasion par le type sauvage et la coexistence. Par exemple, dans la limite $r=0$ avec une conversion zygotique, le forçage envahit si $s < s_{2,z} = \frac{c}{2c + h(1-c)}$. Si $s$ dépasse ce seuil, le forçage est éliminé.

Signification et revendications
L'article prétend délimiter les conditions spécifiques sous lesquelles la dynamique démographique modifie fondamentalement la propagation des forçages génétiques. En connectant la limite $r=0$ aux modèles épidémiologiques SI, les auteurs fournissent des outils analytiques pour prédire l'échec du forçage génique dans les populations à faibles taux de croissance intrinsèque.

L'étude souligne que supposer une densité de population constante (la limite $r \to +\infty$) peut conduire à des prédictions trop optimistes du succès du forçage génique. Plus précisément, la nature de l'onde (tirée vs poussée) et le succès de l'invasion peuvent changer radicalement selon la capacité de la population à se remettre des coûts de fitness imposés par le forçage. Les auteurs concluent que la dynamique démographique n'est pas simplement un facteur de fond, mais peut déterminer si un forçage génique réussit à se propager, échoue ou change la composition génétique d'une population dans un environnement spatialement structuré.