Comprendre le modèle de Wright–Fisher et la dérive génétique
Le simulateur de dérive génétique des populations s'appuie sur le modèle diploïde neutre de Wright–Fisher pour illustrer comment les fréquences alléliques fluctuent de manière aléatoire au fil des générations. Dans une population finie, le tirage aléatoire des gamètes lors de la reproduction induit des écarts par rapport aux fréquences théoriques de la génération parente. Ce phénomène purement stochastique, appelé dérive génétique, conduit inévitablement à la perte de diversité génétique ou à la fixation d'allèles, même en l'absence de toute pression sélective.
Le simulateur permet de configurer une population initiale et d'observer en temps réel le comportement de plusieurs réplicats indépendants. On constate ainsi que, bien que chaque population suive une trajectoire unique et imprévisible, l'ensemble des trajectoires obéit à des lois statistiques rigoureuses définies par la taille de la population et la fréquence allélique initiale.
Configuration de la population et paramètres d'entrée
Pour lancer une simulation, l'utilisateur doit configurer plusieurs paramètres clés dans l'interface:
- Scénarios (Presets): Des configurations prédéfinies sont disponibles pour observer rapidement des dynamiques classiques: Petite population, Grande population, ou Allèle rare.
- Taille efficace de la population, Ne: Saisie sous le libellé
taille efficace de la population. Ce paramètre représente le nombre d'individus reproducteurs et doit être un nombre entier compris entre 1 et 100,000. - Fréquence allélique initiale, p₀: Saisie sous le libellé
fréquence allélique initiale. Ce paramètre correspond à la proportion de l'allèle cible au départ, exprimée sous forme deproportion, 0-1. La valeur doit être comprise entre 0 et 1, limites incluses. - Générations: Saisie sous le libellé
générations. Il s'agit du nombre de générations à simuler, qui doit être un nombre entier de 1 à 2,000. - Populations indépendantes: Saisie sous le libellé
populations indépendantes. Ce paramètre définit le nombre de réplicats simulés en parallèle et doit être un nombre entier compris entre 1 et 200. - Graine aléatoire: Un champ textuel de 100 caractères ou moins permettant de reproduire exactement les mêmes trajectoires de simulation. Un bouton
Nouvelle grainepermet de générer une valeur aléatoire.
Le navigateur applique une limite de charge de travail pour garantir la fluidité de l'exécution. Le nombre total de tirages de copies d'allèles prévus, calculé par la formule 2Ne × générations × répétitions, ne doit pas dépasser la limite de 5 000 000 de tirages. Si cette limite est franchie, le message d'erreur suivant s'affiche: Cette configuration prévoit ‹draws› tirages de copies d'allèles, ce qui dépasse la limite du navigateur de ‹limit›. Réduisez la taille de la population, le nombre de générations ou les répétitions.
Analyse des résultats et graphiques
Une fois la simulation exécutée, l'outil affiche un graphique interactif et un ensemble de statistiques résumant le destin des populations simulées.
Le graphique, intitulé Trajectoires des fréquences alléliques (avec la description d'accessibilité Fréquence allélique par génération pour des populations simulées répétées), présente l'évolution des fréquences au fil du temps:
- L'axe des abscisses (X) est intitulé Génération.
- L'axe des ordonnées (Y) est intitulé Fréquence allélique p.
- Une ligne de référence horizontale indique la valeur attendue: p attendu = p₀.
- Une légende précise la proportion de trajectoires affichées: Affichage de
‹shown›trajectoires sur‹total›; chaque répétition est incluse dans le résumé.
Sous le graphique, les statistiques de synthèse permettent d'analyser les résultats globaux à la génération finale T:
- Fixé à p = 1: La proportion de populations où l'allèle a atteint une fréquence de 100 %.
- Perdu à p = 0: La proportion de populations où l'allèle a totalement disparu.
- Non absorbé à T: La proportion de populations toujours en dérive à la génération finale.
- p final moyen: La moyenne des fréquences alléliques finales sur l'ensemble des réplicats.
- H attendu à T: L'hétérozygotie théorique attendue à la génération finale.
- H moyen simulé à T: L'hétérozygotie moyenne observée expérimentalement dans les simulations.
- Génération moyenne d'absorption: La génération moyenne à laquelle les populations ont atteint l'un des deux états d'absorption (0 ou 1). Si aucune population n'a été absorbée, l'outil affiche non atteinte.
Formule, substitution et précision
La section Formule, substitution et précision détaille les fondements mathématiques de la simulation:
À chaque génération, la transition de Wright–Fisher échantillonne le nombre d'allèles suivant et le divise par les copies de gènes disponibles: Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne).
Avant l'échantillonnage, ces moments décrivent la plage des prochaines générations possibles: E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne).
La proportion attendue de paires d'allèles différents diminue sur t générations comme suit: H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ.
Pour la première génération, le simulateur applique ces équations de manière concrète:
- La première génération utilise 2 ×
‹population›=‹copies›tirages de copies d'allèles. - E[X₁] =
‹copies›×‹frequency›=‹mean›copies de A. - Pour la première génération, Var(X₁) =
‹copies›×‹frequency›× (1 −‹frequency›) =‹variance›; SD(p₁) =‹sd›. - La première répétition tracée a échantillonné X₁ =
‹sampled›copies de A, donc p₁ =‹sampled›÷‹copies›=‹frequency›.
Chaque tirage binomial utilise un échantillonnage de Bernoulli exact; aucune approximation normale n'est utilisée. Les calculs conservent la double précision de JavaScript. Les fréquences se situent sur des incréments de 1 ÷ (2Ne); l'affichage s'arrondit au maximum à six décimales.
Méthode, hypothèses et limites
Cette page implémente le modèle diploïde neutre de Wright–Fisher décrit par Nature Education et le manuel Stanford Human Genetics. Chaque nouvelle génération est un échantillon aléatoire de 2Nₑ copies d'allèles de la précédente.
Le modèle repose sur les hypothèses simplificatrices suivantes:
- Taille efficace de la population constante au cours du temps.
- Générations discrètes et sans chevauchement.
- Échantillonnage aléatoire indépendant (panmixie).
- Neutralité sélective totale (pas de sélection naturelle).
- Un unique locus à deux allèles non liés.
- Absence totale de mutation et de migration (pas de flux génique).
En raison de ces hypothèses, ce simulateur ne doit pas être utilisé comme un outil de prévision pour des populations réelles où la sélection, les variations démographiques, un déséquilibre des sexes ou du succès reproducteur, le chevauchement des générations, la consanguinité, la liaison génétique, le flux génique ou la structure de la population jouent un rôle significatif. De plus, un coefficient de sécurité d'ingénierie ne s'applique pas à ce modèle biologique stochastique.
Traitement des données et confidentialité
Le traitement des données s'effectue localement. Vos paramètres et les populations simulées restent dans ce navigateur. Aucune donnée n'est téléversée ou partagée avec des serveurs externes.
FAQ
Dois-je saisir la taille de recensement ou la taille efficace de la population?
Saisissez la taille efficace de la population, Nₑ: la taille d'une population reproductrice idéalisée qui dériverait au même rythme que la population que vous souhaitez étudier. La taille de recensement peut être beaucoup plus grande lorsque les ratios de masculinité, le succès reproducteur, les fluctuations démographiques ou la structure de la population sont inégaux.
Pourquoi la fréquence allélique moyenne peut-elle rester proche de p₀ alors que l'hétérozygotie diminue?
La dérive neutre n'a pas de direction privilégiée sur un grand nombre de populations indépendantes, de sorte que la fréquence moyenne reste p₀ en espérance. Cependant, les populations individuelles se dispersent vers 0 ou 1, et ces deux limites présentent une hétérozygotie nulle. La variation est donc perdue alors même que la moyenne inter-populations reste presque inchangée.
La probabilité neutre de fixation finale est-elle égale à la fréquence initiale?
Oui, dans ce modèle neutre sans mutation ni migration, un allèle commençant à la fréquence p₀ a une probabilité de fixation finale égale à p₀ et une probabilité de perte égale à 1 - p₀. Les taux affichés ci-dessus ne couvrent que le nombre de générations et de répétitions saisi, ils varient donc selon la graine.
Cela peut-il prédire le comportement d'une population réelle?
Non. Elle montre les résultats selon un modèle neutre de Wright–Fisher délibérément idéalisé. Des prévisions réelles nécessitent des estimations défendables de la taille efficace ainsi que de toute sélection, mutation, migration, variation de la taille de la population, chevauchement des générations, liaison génétique et structure de population qui importent pour le cas étudié.