Effect of population structure and stabilizing selection on quantitative genetic variation

Dit onderzoek toont aan dat bij een kritieke migratiedrempel in een gesubdivideerde populatie onder stabiliserende selectie de genetische variantie sterk toeneemt, de vergelijkbaarheid van de genetische basis tussen subpopulaties maximaal is en een analytisch gemiddeld-veldmodel nauwkeurige voorspellingen doet.

De kern

Titel: Hoe de "buurman" en de "stabiliteit" de variatie in een populatie bepalen

Stel je voor dat je een grote groep mensen hebt die allemaal een bepaalde eigenschap hebben, zoals hun lengte. In de natuurkunde en biologie noemen we dit een kwantitatief kenmerk. De vraag die deze auteurs (Juan Li, Joachim Hermisson en Himani Sachdeva) zich stellen, is: Hoeveel variatie blijft er over in deze eigenschap als de populatie verdeeld is in kleine groepjes (eilanden) en er een constante druk is om "gemiddeld" te blijven?

Hier is de uitleg, vertaald naar alledaagse taal met een paar creatieve vergelijkingen.

1. Het Speelveld: Eilanden en de "Gouden Middenweg"

Stel je een archipel voor met honderden eilanden. Op elk eiland woont een kleine groep mensen. Iedereen op elk eiland heeft een bepaalde "lengte".

• Stabiliserende selectie: Er is een perfecte, ideale lengte (bijvoorbeeld 1,80m). Als je te kort of te lang bent, heb je minder kans om kinderen te krijgen. De natuur wil dus dat iedereen rond die 1,80m blijft.
• Migratie: Af en toe reizen mensen van het ene eiland naar het andere. Dit is de "migratie".
• Mutatie: Soms wordt er per ongeluk een baby geboren die net iets anders is dan zijn ouders (een kleine mutatie).

De onderzoekers willen weten: Als we deze eilanden hebben, hoe groot is dan de variatie in lengte? Is er meer variatie dan als iedereen op één groot eiland zou wonen?

2. De Grote Verrassing: De "Kritieke Drempel"

Het meest interessante resultaat is dat er een kritieke drempel is voor migratie. Denk hierbij aan een schuifregelaar.

• Te weinig verkeer (Laag migratie): Als de eilanden bijna gescheiden zijn, gaat iedereen op zijn eigen eiland zijn eigen weg. Sommige eilanden worden heel kort, andere heel lang, maar binnen elk eiland zijn de mensen vrijwel allemaal hetzelfde. De variatie tussen de eilanden is groot, maar binnen een eiland is er weinig variatie.
• Te veel verkeer (Hoog migratie): Als er heel veel mensen tussen de eilanden reizen, mengen de groepen zich volledig. Het wordt één grote, homogene massa. Iedereen is ongeveer even lang. De variatie is laag, omdat iedereen naar hetzelfde gemiddelde wordt getrokken.
• De Gouden Middenweg (De kritieke drempel): Dit is waar het magie gebeurt. Op het moment dat er precies genoeg verkeer is om mensen te mengen, maar niet genoeg om iedereen gelijk te maken, ontstaat er een explosie van variatie.
  • De analogie: Stel je voor dat je een bak met rode en blauwe balletjes hebt. Als je ze niet mengt, heb je een bak rode en een bak blauwe. Als je ze heel hard mengt, krijg je een bak paars (gemengd). Maar als je ze net goed mengt, krijg je een prachtige, chaotische mix waar je alle mogelijke tinten kunt vinden. Op dit punt is de variatie binnen elke kleine groep het grootst.

3. De "Genetische Identiteit" en de Portabiliteit van GWAS

Een ander belangrijk punt is: Zijn de genen die verantwoordelijk zijn voor de variatie op eiland A, ook verantwoordelijk op eiland B?

• Bij de kritieke drempel: De onderzoekers ontdekten iets tegenintuïtiefs. Op het moment dat de variatie het grootst is (bij de kritieke migratie), zijn de genetische oorzaken van die variatie op de verschillende eilanden het meest op elkaar gelijkend.
• Waarom is dit belangrijk? Dit heeft te maken met GWAS (Genoomwijd Associatie Onderzoek). Dit is een methode waarbij wetenschappers proberen te vinden welke genen bij welke eigenschap horen. Vaak doen ze dit op één populatie en hopen ze dat de resultaten ook gelden voor een andere populatie (bijvoorbeeld: "Dit gen zorgt voor lengte in Nederland, dus zal het dat ook doen in België?").
• De les: Als de populaties te ver uit elkaar liggen (te weinig migratie) of te veel gemengd zijn, zijn de genen die de variatie veroorzaken vaak verschillend. Maar op het "gouden moment" van de kritieke migratie, werken dezelfde genen op alle eilanden. Dit betekent dat voorspellingen over eigenschappen (zoals ziektegevoeligheid of lengte) makkelijker overdraagbaar zijn tussen populaties op dat specifieke moment.

4. De "Effectieve" Krachten

De auteurs hebben ook een slimme wiskundige truc bedacht. Ze laten zien dat je de complexe interacties tussen duizenden genen kunt samenvatten in twee "effectieve" getallen:

1. Effectieve migratie: Migratie voelt alsof het minder sterk is dan het echt is, omdat mensen met een "minder fit" genetische achtergrond (bijvoorbeeld een mix van twee verschillende eilanden) minder succesvol zijn in de voortplanting.
2. Effectieve selectie: De druk om "gemiddeld" te blijven, werkt sterker of zwakker afhankelijk van hoe de genen met elkaar samenwerken.

Met deze twee getallen kunnen ze de hele complexe situatie heel nauwkeurig voorspellen, zonder dat ze elke individuele persoon hoeven te simuleren.

Samenvatting in één zin

Deze studie laat zien dat in een wereld met kleine, gescheiden groepen, de grootste diversiteit en de beste voorspelbaarheid van eigenschappen niet optreedt als de groepen volledig gescheiden zijn of volledig gemengd, maar precies op het moment dat ze net genoeg contact hebben om te "mixen" zonder te "smelten".

De kernboodschap voor de leek:
Net als bij het bakken van een cake: als je de ingrediënten niet mengt, heb je losse bloem en losse suiker. Als je ze te lang mengt, heb je een egaal beslag. Maar op het moment dat je ze net goed hebt gemengd, heb je de perfecte structuur voor variatie. En op dat moment weten we het beste welke ingrediënt (gen) voor welk resultaat zorgt.

Technische samenvatting

Probleemstelling

De auteurs onderzoeken hoe genetische variatie voor polygene kwantitatieve eigenschappen wordt gehandhaafd in een gesubdivideerde populatie (bestaande uit subpopulaties of "eilanden") onder invloed van ruimtelijk homogene stabiliserende selectie. Hoewel eerder werk heeft laten zien dat variatie in tijd of ruimte genetische diversiteit kan vergroten, ontbreekt er een analytisch inzicht in het samenspel van mutatie, stabiliserende selectie, genetische drift en migratie voor eigenschappen die worden bepaald door vele loci met kleine tot gemiddelde effectgroottes.

Specifiek richten de auteurs zich op de volgende vragen:

1. Hoeveel variatie kan worden gehandhaafd op het niveau van individuele loci (genische variatie) en op het niveau van het fenotype (genetische variatie) binnen en tussen subpopulaties?
2. Op welke manier beïnvloedt de populatiestructuur de "draagbaarheid" van Genome-Wide Association Studies (GWAS) tussen subpopulaties? Dat wil zeggen, in welke mate verklaren loci met een groot effect in de ene subpopulatie ook variatie in andere subpopulaties?
3. Kunnen we deze complexe dynamiek nauwkeurig voorspellen met analytische benaderingen die rekening houden met linkagedisequilibrium (LD) tussen loci?

Methodologie

De studie combineert geavanceerde theoretische modellering met individuen-gebaseerde simulaties.

1. Het Model:

• Populatiestructuur: Een oneindig-eilandmodel (infinite-island model) met $D$ demes (subpopulaties), elk met $N$ diploïde individuen. Migratie vindt plaats via een migrantenpool met een migratiekans $m$.
• Eigenschap: Een additieve kwantitatieve eigenschap wordt bepaald door $L$ ongekoppelde, biallelische loci. Allelen hebben effectgroottes $\alpha_i$.
• Selectie: Stabiliserende selectie met een ruimtelijk uniform optimum ($z_{opt}$) en een fitnessfunctie met breedte $V_s$.
• Dynamiek: Het model omvat mutatie, selectie, migratie en genetische drift.

2. Analytische Benadering:

• Linkage Equilibrium (LE) Benadering: De auteurs behandelen eerst de eigenschap als een verzameling onafhankelijk evoluerende loci. Stabiliserende selectie op een polygene eigenschap manifesteert zich lokaal als onderdominante selectie (heterozygote nadeel) op individuele loci. Ze gebruiken de diffusie-benadering (diffusion approximation) om de evenwichtsverdeling van allelfrequenties af te leiden.
• Effectieve Parameters (LD Benadering): Omdat stabiliserende selectie en migratie negatief LD creëren (associaties tussen loci met tegengestelde effecten), passen de auteurs een "mean-field" aanpak toe. Ze introduceren effectieve parameters:
  • Een effectieve selectiecoëfficiënt ($s_{eff}$) die de vermindering van heterozygote nadeel door negatief LD binnen een deme compenseert.
  • Een effectieve migratiesnelheid ($m_{eff}$) die de selectie tegen immigranten (vanwege hun minder fitte genetische achtergrond) kwantificeert.
• Deze effectieve parameters worden iteratief opgelost om voorspellingen te doen voor genische en genetische variantie, $F_{ST}$ (differentiatie op locusniveau) en $Q_{ST}$ (differentiatie op eigenschapsniveau).

3. Simulaties:

• Discrete-generatie simulaties worden uitgevoerd met $D=100$ demes en $N=200$ individuen per deme.
• De simulaties omvatten zowel loci met gelijke effectgroottes als loci met effectgroottes getrokken uit een exponentiële verdeling.
• De resultaten worden vergeleken met de analytische LE en LD voorspellingen.

Belangrijkste Resultaten

1. Kritieke Migratiedrempel en Genische Variatie:

• Er bestaat een kritieke migratiedrempel ($N_{m_{crit}}$), afhankelijk van de effectgrootte van de loci ($s = \alpha^2/2V_s$) en de mutatiegraad.
• Onder de drempel ($Nm < N_{m_{crit}}$): Subpopulaties kunnen verschillende adaptieve pieken bezetten (verschillende allelen zijn dominant in verschillende demes). Dit leidt tot een sterke inflatie van de totale genische variatie ($V_T^{(gen)}$) in de gehele populatie, vaak veel hoger dan in een panmictische populatie.
• Boven de drempel ($Nm > N_{m_{crit}}$): De evolutie van de demes is sterk gekoppeld; hetzelfde allel wordt dominant in de meeste demes, wat de totale variatie doet dalen naar het niveau van een panmictische populatie.
• Maximale variatie binnen demes: De variatie binnen een individuele subpopulatie ($V_W^{(gen)}$) is maximaal dicht bij de kritieke migratiedrempel. Bij zeer lage migratie is de variatie laag (fixatie), en bij zeer hoge migratie is deze ook lager door homogenisatie.

2. Differentiatie ($F_{ST}$ vs. $Q_{ST}$):

• $F_{ST}$: Bij lage migratie is de differentiatie op locusniveau ($F_{ST}$) hoog en vaak hoger dan bij neutrale loci. Bij hogere migratie daalt $F_{ST}$ en kan deze zelfs onder de neutrale verwachting zakken.
• $Q_{ST}$: De differentiatie in de eigenschapswaarde ($Q_{ST}$) is door stabiliserende selectie aanzienlijk lager dan $F_{ST}$. Dit komt omdat stabiliserende selectie de fenotypische variatie tussen demes onderdrukt, zelfs als de onderliggende genetische variatie hoog is.
• De auteurs leiden een analytische formule af: $Q_{ST} \approx \frac{1}{1 + 4Nm + 4N(V_W/V_s)}$.

3. Invloed op GWAS en "Portability":

• Overeenkomst in genetische basis: De genetische basis van variatie tussen subpopulaties is het meest vergelijkbaar (dezelfde loci dragen bij) dicht bij de kritieke migratiedrempel.
• Paradoxaal effect: Bij hogere migratie (waar men zou verwachten dat populaties meer op elkaar lijken) neemt de "portability" van GWAS-resultaten af. Hoewel hoge migratie leidt tot dezelfde allelfrequenties, zorgt dit ervoor dat individuele loci minder variatie bijdragen (ze fixeren), waardoor het moeilijker wordt om significante associaties te vinden die over populaties heen gelden.
• Outliers: Loci met een atypisch grote bijdrage aan de variatie (outliers) zijn het meest frequent en gedeeld tussen demes bij intermediaire migratie.

4. Nauwkeurigheid van de Benadering:

• De analytische benadering met effectieve parameters (LD-benadering) voorspelt de simulatieresultaten zeer nauwkeurig, zelfs voor complexe scenario's met verschillende effectgroottes.
• De LE-benadering (zonder LD) mist de kwantitatieve details, vooral bij sterke selectie, maar vangen wel het kwalitatieve gedrag.

Bijdrage en Relevantie

Theoretische Vooruitgang:
Het artikel vult een belangrijke leemte in de evolutionaire genetica door een analytisch raamwerk te bieden dat drift, selectie, migratie en LD combineert voor polygene eigenschappen. Het toont aan dat de complexe effecten van LD kunnen worden samengevat in effectieve parameters, wat de voorspelbaarheid van het model aanzienlijk verbetert.

Implicaties voor Empirisch Onderzoek:

• GWAS en Voorspelling: De resultaten waarschuwen dat de overdraagbaarheid van polygene scores tussen populaties niet alleen afhangt van de genetische afstand, maar ook van de migratiegeschiedenis en de sterkte van de selectie. Intermediaire migratie kan de gedeelde variatie maximaliseren, terwijl zeer hoge migratie deze juist kan verminderen door het "verlies" van grote bijdragen van individuele loci.
• Interpretatie van $Q_{ST}$: Het paper biedt een betere interpretatie van $Q_{ST}$ in vergelijking met $F_{ST}$. Een lage $Q_{ST}$ ten opzichte van $F_{ST}$ is een sterk signaal van ruimtelijk homogene stabiliserende selectie, maar de grootte van het verschil hangt af van de mutatiegraad en de populatiestructuur.
• Cryptische Variatie: Populatiestructuur kan grote hoeveelheden "cryptische" genetische variatie opslaan (variatie die niet zichtbaar is in het fenotype binnen een deme, maar wel in de totale populatie). Deze variatie kan worden "ontgrendeld" door kruisingen tussen subpopulaties.

Conclusie:
De studie demonstreert dat ruimtelijke populatiestructuur, zelfs onder uniforme stabiliserende selectie, een krachtige drijver is voor het handhaven van genetische diversiteit. De interactie tussen migratie en selectie creëert een niet-monotoon patroon van variatie, waarbij intermediaire migratie niveaus optimaal zijn voor het behoud van variatie binnen subpopulaties en voor de gedeelde genetische architectuur van eigenschappen. Dit heeft directe gevolgen voor hoe we evolutionaire geschiedenissen reconstrueren en hoe we genetische voorspellingen in menselijke en niet-menselijke populaties interpreteren.