Genotype specificity and spatial arrangement govern the direction and magnitude of selection in variable environments

Résumé technique : La spécificité du génotype et l'agencement spatial régissent la sélection dans des environnements variables

Énoncé du problème
L'hétérogénéité environnementale est une caractéristique omniprésente des systèmes biologiques, allant des gradients de nutriments dans les biofilms microbiens aux gradients d'oxygène et de médicaments dans les tumeurs solides. Bien qu'il soit établi que la structure spatiale influence les résultats évolutifs, la compréhension théorique de la manière dont la variation environnementale interagit avec la structure de la population pour déterminer la probabilité de fixation des mutants bénéfiques reste fragmentée. Des études antérieures ont rapporté des résultats contradictoires : certains modèles suggèrent que l'hétérogénéité spatiale amplifie la sélection, tandis que d'autres indiquent qu'elle la supprime. Les cadres existants échouent souvent à rendre compte systématiquement des interactions spécifiques au génotype — spécifiquement, que la variation environnementale agisse de manière symétrique sur les génotypes en compétition ou de manière asymétrique sur l'un d'eux. Ce manque de principe unificateur empêche la réconciliation de résultats théoriques disparates à travers différentes structures spatiales et distributions de fitness.

Méthodologie
Les auteurs développent un cadre unifié pour la sélection sur des graphes en réseau incorporant l'hétérogénéité environnementale. Le modèle utilise un processus de naissance-mort de Moran en temps discret sur des graphes non orientés (spécifiquement des cycles unidimensionnels et des réseaux carrés bidimensionnels).

• Paysage de fitness : La fitness d'un individu au nœud $i$ est définie par $f_{\alpha,i} = r_\alpha + \sigma_\alpha c_i$, où $r_\alpha$ est la fitness intrinsèque du génotype $\alpha$ (mutant $A$ ou résident $B$), $\sigma_\alpha$ est l'amplitude de l'hétérogénéité spécifique au génotype, et $c_i \in \{+1, -1\}$ représente l'état environnemental local (riche ou pauvre).
• Scénarios d'interaction : L'étude distingue explicitement trois classes d'interactions génotype-environnement :
  1. S1 (Symétrique au génotype) : Les deux génotypes répondent également à la variation environnementale ($\sigma_A = \sigma_B = \sigma$).
  2. S2 (Spécifique au mutant) : Seul le mutant répond à la variation environnementale ($\sigma_A = \sigma, \sigma_B = 0$).
  3. S3 (Spécifique au résident) : Seul le résident répond à la variation environnementale ($\sigma_A = 0, \sigma_B = \sigma$).
• Agencement spatial : La distribution spatiale des états environnementaux est quantifiée à l'aide d'un indice de corrélation spatiale ($M$). Ce paramètre mesure le degré de regroupement par opposition au mélange des sites riches et pauvres, allant d'un mélange maximal (échiquier, $M$ faible) à une ségrégation complète (domaines regroupés, $M$ élevé).
• Analyse : Les auteurs emploient une combinaison de dérivations analytiques pour les configurations limites (échiquier et ségrégué) utilisant des méthodes de martingales et des équations de Kolmogorov en arrière, parallèlement à des simulations numériques extensives et des simulations stochastiques sur 946 configurations environnementales sur des réseaux 1D et 2D.

Contributions clés
La contribution principale de ce travail est l'identification de deux principes directeurs qui déterminent la direction et l'ampleur de la sélection dans des environnements hétérogènes :

1. La spécificité du génotype détermine la direction : Les auteurs démontrent que la direction de l'effet (amplification vs suppression de la sélection) est principalement dictée par le génotype dont la fitness est modulée par l'environnement.

   • Spécifique au résident (S3) : L'hétérogénéité amplifie systématiquement la sélection, augmentant la probabilité de fixation des mutants bénéfiques.
   • Spécifique au mutant (S2) : L'hétérogénéité supprime généralement la sélection, diminuant la probabilité de fixation, en particulier sous une forte hétérogénéité.
   • Symétrique au génotype (S1) : L'hétérogénéité produit un effet d'amplification plus faible que dans le cas S3, résultant souvent d'effets de frontière du second ordre.

2. L'agencement spatial détermine l'ampleur : L'ampleur de l'amplification ou de la suppression est modulée par l'indice de corrélation spatiale ($M$), sans inverser la direction établie par la spécificité du génotype.

   • Dans les scénarios spécifiques au résident, les environnements fortement mélangés (faible $M$) produisent l'amplification la plus forte.
   • Dans les scénarios spécifiques au mutant, les environnements fortement regroupés (fort $M$) contrecarrent la suppression, conduisant à des probabilités de fixation plus élevées que dans les environnements mélangés.
   • Dans les scénarios symétriques au génotype, la relation est non monotone, les probabilités de fixation étant souvent minimisées à des niveaux intermédiaires de mélange spatial.

Résultats

• Limites analytiques : Les solutions analytiques exactes pour les configurations en échiquier (bipartites) confirment que les probabilités de fixation dépendent de la moyenne géométrique effective des rapports de fitness, validant le « théorème isotherme environnemental » pour ces topologies spécifiques. Les approximations pour les configurations ségréguées montrent que la fixation se comporte comme une moyenne de deux environnements homogènes.
• Validation numérique : Les simulations à travers un continuum de corrélations spatiales confirment que les limites analytiques bornent les probabilités de fixation observées dans des configurations aléatoires ou intermédiaires.
• Réconciliation de la littérature : En cartographiant les études antérieures sur ce cadre, les auteurs montrent que les résultats apparemment contradictoires de la littérature sont cohérents lorsque la spécificité du génotype et l'agencement spatial sont explicitement pris en compte. Les études rapportant une amplification correspondent aux modèles ségrégués spécifiques au résident ou faiblement hétérogènes, tandis que celles rapportant une suppression correspondent aux modèles aléatoires spécifiques au mutant ou symétriques.
• Connectivité du graphe : L'étude note que, bien que la spécificité du génotype reste prédictive à travers les types de graphes, l'influence de l'agencement spatial ($M$) diminue à mesure que la connectivité du graphe augmente (approchant la limite du graphe complet), où les corrélations spatiales deviennent sans importance.

Signification
L'article prétend fournir un cadre théorique unifié qui résout les incohérences de longue date en théorie des graphes évolutifs concernant l'hétérogénéité environnementale. En distinguant entre la modulation spécifique au génotype et la modulation symétrique au génotype, les auteurs offrent des critères prédictifs pour l'adaptation dans des paysages hétérogènes. Les résultats suggèrent que l'hétérogénéité environnementale n'amplifie ni ne supprime intrinsèquement la sélection ; plutôt, le résultat émerge de l'interaction entre quel génotype est affecté et comment l'environnement est agencé spatialement. Ce cadre est présenté comme applicable à divers contextes biologiques, y compris les communautés microbiennes, l'évolution somatique et la progression du cancer, où la compréhension de l'interaction spécifique entre le génotype et la niche spatiale est cruciale pour prédire les trajectoires évolutives.