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

Technische Zusammenfassung: Genotypspezifität und räumliche Anordnung bestimmen die Selektion in variablen Umgebungen

Problemstellung
Umweltheterogenität ist ein allgegenwärtiges Merkmal biologischer Systeme, das von Nährstoffgradienten in mikrobiellen Biofilmen bis hin zu Sauerstoff- und Wirkstoffgradienten in soliden Tumoren reicht. Zwar ist bekannt, dass die räumliche Struktur evolutionäre Ergebnisse beeinflusst, doch das theoretische Verständnis dafür, wie Umweltvariation mit der Populationsstruktur interagiert, um die Fixierungswahrscheinlichkeit vorteilhafter Mutanten zu bestimmen, bleibt fragmentiert. Frühere Studien haben widersprüchliche Ergebnisse berichtet: Einige Modelle deuten darauf hin, dass räumliche Heterogenität die Selektion verstärkt, während andere anzeigen, dass sie sie unterdrückt. Bestehende Rahmenwerke versagen oft darin, genotypspezifische Interaktionen systematisch zu berücksichtigen – insbesondere, ob Umweltvariation symmetrisch auf konkurrierende Genotypen oder asymmetrisch auf einen einzigen wirkt. Das Fehlen eines vereinheitlichenden Prinzips verhindert die Versöhnung unterschiedlicher theoretischer Ergebnisse über verschiedene räumliche Strukturen und Fitnessverteilungen hinweg.

Methodik
Die Autoren entwickeln ein vereinheitlichtes Rahmenwerk für die Selektion auf Gittergraphen unter Einbeziehung von Umweltheterogenität. Das Modell verwendet einen diskretzeitlichen Moran-Geburt-Todes-Prozess auf ungerichteten Graphen (speziell eindimensionale Zyklen und zweidimensionale quadratische Gitter).

• Fitnesslandschaft: Die Fitness eines Individuums am Knoten $i$ ist definiert als $f_{\alpha,i} = r_\alpha + \sigma_\alpha c_i$, wobei $r_\alpha$ die inhärente Fitness des Genotyps $\alpha$ (Mutante $A$ oder Resident $B$) ist, $\sigma_\alpha$ die genotypspezifische Heterogenitätsamplitude darstellt und $c_i \in \{+1, -1\}$ den lokalen Umweltzustand (reich oder arm) repräsentiert.
• Interaktionsszenarien: Die Studie unterscheidet explizit drei Klassen von Genotyp-Umwelt-Interaktionen:
  1. S1 (Genotyp-symmetrisch): Beide Genotypen reagieren gleich auf Umweltvariation ($\sigma_A = \sigma_B = \sigma$).
  2. S2 (Mutantenspezifisch): Nur die Mutante reagiert auf Umweltvariation ($\sigma_A = \sigma, \sigma_B = 0$).
  3. S3 (Residentenspezifisch): Nur der Resident reagiert auf Umweltvariation ($\sigma_A = 0, \sigma_B = \sigma$).
• Räumliche Anordnung: Die räumliche Verteilung der Umweltzustände wird mittels eines räumlichen Korrelationsindex ($M$) quantifiziert. Dieser Parameter misst den Grad der Clusterbildung versus Durchmischung reicher und armer Standorte und reicht von maximal durchmischt (Schachbrettmuster, niedriger $M$) bis vollständig segregiert (geclusterte Domänen, hoher $M$).
• Analyse: Die Autoren verwenden eine Kombination aus analytischen Herleitungen für Grenzkonfigurationen (Schachbrettmuster und segregiert) unter Verwendung von Martingal-Methoden und rückwärtigen Kolmogorov-Gleichungen sowie umfangreichen numerischen und stochastischen Simulationen über 946 Umweltkonfigurationen auf 1D- und 2D-Gittern.

Hauptbeiträge
Der primäre Beitrag dieser Arbeit ist die Identifizierung zweier leitender Prinzipien, die Richtung und Ausmaß der Selektion in heterogenen Umgebungen bestimmen:

1. Genotypspezifität bestimmt die Richtung: Die Autoren zeigen, dass die Richtung des Effekts (Verstärkung versus Unterdrückung der Selektion) primär davon bestimmt wird, welcher Genotyp durch die Umwelt in seiner Fitness moduliert wird.

   • Residentenspezifisch (S3): Heterogenität verstärkt die Selektion konsistent und erhöht die Fixierungswahrscheinlichkeit vorteilhafter Mutanten.
   • Mutantenspezifisch (S2): Heterogenität unterdrückt die Selektion im Allgemeinen und verringert die Fixierungswahrscheinlichkeit, insbesondere unter starker Heterogenität.
   • Genotyp-symmetrisch (S1): Heterogenität erzeugt einen schwächeren Verstärkungseffekt im Vergleich zu S3, der oft aus Sekundär-Effekten an den Grenzen resultiert.

2. Räumliche Anordnung bestimmt das Ausmaß: Das Ausmaß der Verstärkung oder Unterdrückung wird durch den räumlichen Korrelationsindex ($M$) moduliert, ohne die durch die Genotypspezifität festgelegte Richtung umzukehren.

   • In residentenspezifischen Szenarien liefern stark durchmischte (niedriger $M$) Umgebungen die stärkste Verstärkung.
   • In mutantenspezifischen Szenarien wirken stark geclusterte (hoher $M$) Umgebungen der Unterdrückung entgegen und führen zu höheren Fixierungswahrscheinlichkeiten als in durchmischten Umgebungen.
   • In genotyp-symmetrischen Szenarien ist die Beziehung nicht-monoton, wobei die Fixierungswahrscheinlichkeiten oft bei intermediären Niveaus der räumlichen Durchmischung minimiert werden.

Ergebnisse

• Analytische Grenzen: Exakte analytische Lösungen für Schachbrett-(bipartite) Konfigurationen bestätigen, dass Fixierungswahrscheinlichkeiten vom effektiven geometrischen Mittel der Fitnessverhältnisse abhängen und damit das „Umwelt-isotherme Theorem" für diese spezifischen Topologien validieren. Näherungen für segregierte Konfigurationen zeigen, dass sich die Fixierung als Durchschnitt zweier homogener Umgebungen verhält.
• Numerische Validierung: Simulationen über ein Kontinuum räumlicher Korrelationen bestätigen, dass die analytischen Grenzen die in zufälligen oder intermediären Konfigurationen beobachteten Fixierungswahrscheinlichkeiten begrenzen.
• Versöhnung der Literatur: Durch die Abbildung früherer Studien auf dieses Rahmenwerk zeigen die Autoren, dass scheinbar widersprüchliche Ergebnisse in der Literatur konsistent sind, wenn Genotypspezifität und räumliche Anordnung explizit berücksichtigt werden. Studien, die eine Verstärkung berichten, stimmen mit residentenspezifischen oder schwach heterogenen segregierten Modellen überein, während jene, die eine Unterdrückung berichten, mit mutantenspezifischen oder symmetrischen zufälligen Modellen übereinstimmen.
• Graph-Konnektivität: Die Studie stellt fest, dass die Genotypspezifität zwar über Graphentypen hinweg vorhersagbar bleibt, der Einfluss der räumlichen Anordnung ($M$) jedoch mit zunehmender Graph-Konnektivität abnimmt (Annäherung an den vollständigen Graphen-Grenzwert), wo räumliche Korrelationen irrelevant werden.

Bedeutung
Die Arbeit behauptet, ein vereinheitlichtes theoretisches Rahmenwerk zu liefern, das langjährige Inkonsistenzen in der evolutionären Graphentheorie bezüglich Umweltheterogenität auflöst. Durch die Unterscheidung zwischen genotypspezifischer und genotyp-symmetrischer Modulation bieten die Autoren prädiktive Kriterien für die Anpassung in heterogenen Landschaften. Die Ergebnisse legen nahe, dass Umweltheterogenität die Selektion nicht inhärent verstärkt oder unterdrückt; vielmehr ergibt sich das Ergebnis aus dem Zusammenspiel davon, welcher Genotyp betroffen ist und wie die Umwelt räumlich angeordnet ist. Dieses Rahmenwerk wird als anwendbar auf diverse biologische Kontexte präsentiert, einschließlich mikrobieller Gemeinschaften, somatischer Evolution und Krebsprogression, wo das Verständnis der spezifischen Interaktion zwischen Genotyp und räumlicher Nische entscheidend für die Vorhersage evolutionärer Trajektorien ist.