Το μοντέλο Wright–Fisher και η γενετική παρέκκλιση
Η γενετική παρέκκλιση αποτελεί μία από τις θεμελιώδεις δυνάμεις της εξέλιξης των ειδών. Ο «Προσομοιωτής Γενετικής Παρέκκλισης Πληθυσμού» επιτρέπει τη μελέτη και την οπτικοποίηση αυτής της διαδικασίας μέσω του ουδέτερου διπλοειδούς μοντέλου Wright–Fisher. Το εργαλείο αυτό προσομοιώνει πώς η τυχαία δειγματοληψία των γαμετών από γενιά σε γενιά οδηγεί στη σταδιακή απώλεια ή την καθήλωση των αλληλομόρφων, χωρίς την επίδραση της φυσικής επιλογής, της μετάλλαξης ή της μετανάστευσης [i, j].
Σε κάθε γενιά, η μετάβαση Wright–Fisher δειγματίζει τον επόμενο αριθμό αλληλομόρφων και τον διαιρεί με τα διαθέσιμα αντίγραφα γονιδίων: Xₜ₊₁ ~ Binomial(2Ne, pₜ)· pₜ₊₁ = Xₜ₊₁ ÷ (2Ne).
Πριν από τη δειγματοληψία, αυτές οι ροπές περιγράφουν το εύρος των πιθανών επόμενων γενεών: E[Xₜ₊₁] = 2Ne·pₜ· Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ)· Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne).
Με αυτόν τον τρόπο, ο χρήστης μπορεί να παρακολουθήσει πώς το δραστικό μέγεθος πληθυσμού, Ne, επηρεάζει άμεσα τη διακύμανση των συχνοτήτων και την ταχύτητα με την οποία οι πληθυσμοί οδηγούνται στην ομοζυγωτία [i, j].
Ρύθμιση πληθυσμού και παράμετροι εισαγωγής
Για την εκτέλεση μιας προσομοίωσης, ο χρήστης ορίζει τις παρακάτω παραμέτρους στην ενότητα «Ρύθμιση πληθυσμού»:
- Σενάρια: Έτοιμες προκαθορισμένες ρυθμίσεις που περιλαμβάνουν τις επιλογές «Μικρός πληθυσμός», «Μεγάλος πληθυσμός» και «Σπάνιο αλληλόμορφο».
- Δραστικό μέγεθος πληθυσμού, Ne: Εισάγεται σε «αναπαραγωγικά άτομα». Πρέπει να είναι ακέραιος αριθμός από 1 έως 100.000.
- Αρχική συχνότητα αλληλομόρφων, p₀: Εισάγεται ως «αναλογία, 0–1». Η τιμή πρέπει να βρίσκεται μεταξύ 0 και 1, συμπεριλαμβανομένων και των δύο ορίων.
- Γενιές: Εισάγεται σε «γενιές». Πρέπει να είναι ακέραιος αριθμός από 1 έως 2.000.
- Ανεξάρτητοι πληθυσμοί: Εισάγεται σε «επαναλήψεις». Πρέπει να είναι ακέραιος αριθμός από 1 έως 200.
- Τυχαίος σπόρος (seed): Ένα πεδίο κειμένου έως 100 χαρακτήρων που επιτρέπει την ακριβή αναπαραγωγή των προσομοιωμένων διαδρομών. Ο χρήστης μπορεί να πατήσει το κουμπί «Νέος σπόρος» για να δημιουργήσει μια νέα τυχαία τιμή.
Μέθοδος, υποθέσεις και όρια
Αυτή η σελίδα υλοποιεί το ουδέτερο διπλοειδές μοντέλο Wright–Fisher όπως περιγράφεται από το Nature Education και το εγχειρίδιο Stanford Human Genetics [j]. Κάθε νέα γενιά είναι ένα τυχαίο δείγμα 2Ne αντιγράφων αλληλομόρφων από την προηγούμενη [i, j]. Έλεγχος πηγών στις 6 Αυγούστου 2026 [i, j].
Μαθηματικές υποθέσεις του μοντέλου
Η λειτουργία του προσομοιωτή βασίζεται αυστηρά στις ακόλουθες παραδοχές:
- Σταθερό δραστικό μέγεθος πληθυσμού.
- Διακριτές μη επικαλυπτόμενες γενιές.
- Τυχαία ανεξάρτητη δειγματοληψία (τυχαία σύζευξη).
- Επιλεκτική ουδετερότητα (όλα τα γονίδια έχουν την ίδια προσαρμοστικότητα).
- Ένας μη συνδεδεμένος τόπος δύο αλληλομόρφων.
- Απουσία μετάλλαξης ή μετανάστευσης.
Μην το χρησιμοποιείτε ως πρόβλεψη όταν παίζει ρόλο η επιλογή, η μεταβολή του πληθυσμού, η ανισορροπία των φύλων ή της αναπαραγωγικής επιτυχίας, οι επικαλυπτόμενες γενιές, η αιμομιξία, η σύνδεση, η γονιδιακή ροή ή η δομή του πληθυσμού [i, j]. Ένας συντελεστής ασφαλείας μηχανικής δεν εφαρμόζεται σε αυτό το στοχαστικό βιολογικό μοντέλο.
Υπολογιστική ακρίβεια και όρια φόρτου
Κάθε διωνυμική κλήρωση χρησιμοποιεί ακριβή δειγματοληψία Bernoulli· δεν χρησιμοποιείται κανονική προσέγγιση. Οι υπολογισμοί διατηρούν διπλή ακρίβεια JavaScript. Οι συχνότητες βρίσκονται σε βήματα του 1 ÷ (2Ne)· η οθόνη στρογγυλοποιεί σε το πολύ έξι δεκαδικά ψηφία.
Για την αποφυγή υπερφόρτωσης του συστήματος, το πρόγραμμα περιήγησης επιβάλλει ένα ανώτατο όριο 5.000.000 προγραμματισμένων κληρώσεων αντιγράφων αλληλομόρφων, το οποίο υπολογίζεται ως: Κληρώσεις = 2Ne × γενιές × επαναλήψεις
Αν η ρύθμιση υπερβαίνει αυτό το όριο, η εκτέλεση εμποδίζεται και εμφανίζεται σχετικό σφάλμα.
Ανάλυση αποτελεσμάτων και στατιστικά στοιχεία
Μετά την εκτέλεση της προσομοίωσης, εμφανίζεται το διάγραμμα με τίτλο «Διαδρομές συχνότητας αλληλομόρφων» (με ARIA label «Συχνότητα αλληλομόρφων ανά γενιά για επαναλαμβανόμενους προσομοιωμένους πληθυσμούς»). Ο οριζόντιος άξονας ορίζεται ως «Γενιά» και ο κατακόρυφος ως «Συχνότητα αλληλομόρφων p». Μια οριζόντια γραμμή αναφοράς δείχνει το «Αναμενόμενο p = p₀».
Κάτω από το διάγραμμα, το μήνυμα «Εμφάνιση ‹shown› από ‹total› διαδρομές· κάθε επανάληψη περιλαμβάνεται στη σύνοψη.» ενημερώνει για το πλήθος των γραμμών που σχεδιάστηκαν.
Τα συγκεντρωτικά αποτελέσματα παρουσιάζουν τις εξής μετρικές:
- Καθηλώθηκε στο p = 1: Το ποσοστό των πληθυσμών όπου το αλληλόμορφο επικράτησε πλήρως.
- Χάθηκε στο p = 0: Το ποσοστό των πληθυσμών όπου το αλληλόμορφο εξαφανίστηκε.
- Δεν απορροφήθηκε έως το T: Οι πληθυσμοί που παρέμειναν πολυμορφικοί στο τέλος της προσομοίωσης.
- Μέσο τελικό p: Η μέση συχνότητα του αλληλομόρφου σε όλες τις επαναλήψεις κατά την τελευταία γενιά.
- Αναμενόμενο H στο T: Η θεωρητική εναπομένουσα ετεροζυγωτία.
- Μέσο προσομοιωμένο H στο T: Η μέση ετεροζυγωτία που παρατηρήθηκε πειραματικά στις επαναλήψεις.
- Μέση γενιά απορρόφησης: Ο μέσος όρος των γενεών που χρειάστηκαν για να φτάσουν οι πληθυσμοί στην καθήλωση ή την απώλεια (εμφανίζει «δεν επιτεύχθηκε» αν κανένας πληθυσμός δεν έφτασε σε αυτά τα όρια).
Τύπος, αντικατάσταση και ακρίβεια
Για την κατανόηση των μαθηματικών υπολογισμών του πρώτου βήματος, ο προσομοιωτής αναλύει τις εξισώσεις χρησιμοποιώντας τις τρέχουσες τιμές εισαγωγής:
- Το αναμενόμενο μερίδιο των ανόμοιων ζευγών αλληλομόρφων μειώνεται κατά τη διάρκεια t γενεών ως εξής: H₀ = 2p₀(1 − p₀)· Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ [i, j].
- Η πρώτη γενιά χρησιμοποιεί 2 ×
‹population›=‹copies›κληρώσεις αντιγράφων αλληλομόρφων. - E[X₁] =
‹copies›×‹frequency›=‹mean›αντίγραφα Α. - Για την πρώτη γενιά, Var(X₁) =
‹copies›×‹frequency›× (1 −‹frequency›) =‹variance›· SD(p₁) =‹sd›. - Η πρώτη σχεδιασμένη επανάληψη δειγματίστηκε X₁ =
‹sampled›αντίγραφα Α, επομένως p₁ =‹sampled›÷‹copies›=‹frequency›.
Διαχείριση σφαλμάτων
Σε περίπτωση λανθασμένης εισαγωγής δεδομένων, η εφαρμογή εμφανίζει συγκεκριμένα μηνύματα σφάλματος:
- Αν μια τιμή δεν είναι αριθμητική:
‹field›: το «‹token›» δεν είναι αριθμός. - Αν εισαχθεί δεκαδικός αριθμός εκεί που απαιτείται ακέραιος:
‹field›: εισαγάγετε έναν ακέραιο αριθμό. - Αν το μέγεθος πληθυσμού είναι εκτός ορίων:
Το ενεργό μέγεθος πληθυσμού πρέπει να είναι ακέραιος αριθμός από 1 έως 100.000. - Αν η αρχική συχνότητα είναι εκτός ορίων:
Η αρχική συχνότητα αλληλομόρφων πρέπει να είναι μεταξύ 0 και 1, συμπεριλαμβανομένων και των δύο ορίων. - Αν οι γενιές είναι εκτός ορίων:
Οι γενιές πρέπει να είναι ακέραιος αριθμός από 1 έως 2,000. - Αν οι επαναλήψεις είναι εκτός ορίων:
Οι ανεξάρτητοι πληθυσμοί πρέπει να είναι ακέραιος αριθμός από 1 έως 200. - Αν λείπει ο σπόρος:
Εισαγάγετε έναν τυχαίο σπόρο (seed) ώστε να είναι δυνατή η αναπαραγωγή της προσομοίωσης. - Αν ο σπόρος είναι πολύ μεγάλος:
Κρατήστε τον τυχαίο σπόρο (seed) στους 100 χαρακτήρες ή λιγότερους. - Αν ξεπεραστεί το υπολογιστικό όριο:
Αυτή η ρύθμιση προβλέπει ‹draws› κληρώσεις αντιγράφων αλληλομόρφων, υπερβαίνοντας το όριο ‹limit› του προγράμματος περιήγησης. Μειώστε το μέγεθος πληθυσμού, τις γενιές ή τις επαναλήψεις.
Προστασία προσωπικών δεδομένων
Όλες οι παράμετροι εισαγωγής και τα δεδομένα των προσομοιωμένων πληθυσμών υποβάλλονται σε επεξεργασία τοπικά. Οι υπολογισμοί εκτελούνται εξ ολοκλήρου μέσα στο πρόγραμμα περιήγησης του χρήστη και καμία πληροφορία δεν μεταφέρεται ή αποθηκεύεται σε εξωτερικούς διακομιστές.
Συχνές Ερωτήσεις (FAQ)
Πρέπει να εισαγάγω το απογραφικό μέγεθος ή το δραστικό μέγεθος πληθυσμού;
Εισαγάγετε το δραστικό μέγεθος πληθυσμού, Ne: το μέγεθος ενός εξιδανικευμένου αναπαραγωγικού πληθυσμού που θα παρουσίαζε παρέκκλιση με τον ίδιο ρυθμό όπως ο πληθυσμός που θέλετε να μελετήσετε. Το απογραφικό μέγεθος μπορεί να είναι πολύ μεγαλύτερο όταν οι αναλογίες των φύλων, η αναπαραγωγική επιτυχία, οι διακυμάνσεις του πληθυσμού ή η δομή του πληθυσμού δεν είναι ομοιόμορφες.
Γιατί η μέση συχνότητα αλληλομόρφων μπορεί να παραμένει κοντά στο p₀ ενώ η ετεροζυγωτία μειώνεται;
Η ουδέτερη γενετική παρέκκλιση δεν έχει προτιμώμενη κατεύθυνση σε πολλούς ανεξάρτητους πληθυσμούς, επομένως η μέση συχνότητα παραμένει p₀ ως προσδοκία [i, j]. Ωστόσο, οι μεμονωμένοι πληθυσμοί διασπείρονται προς το 0 ή το 1, και τα δύο όρια έχουν μηδενική ετεροζυγωτία [i, j]. Συνεπώς, η ποικιλομορφία χάνεται ακόμη και όταν ο μέσος όρος μεταξύ των πληθυσμών παραμένει σχεδόν αμετάβλητος [i, j].
Είναι η ουδέτερη πιθανότητα τελικής καθήλωσης ίση με την αρχική συχνότητα;
Ναι, σε αυτό το ουδέτερο μοντέλο χωρίς μετάλλαξη ή μετανάστευση, ένα αλληλόμορφο με αρχική συχνότητα p₀ έχει τελική πιθανότητα καθήλωσης p₀ και πιθανότητα απώλειας 1 − p₀ [i, j]. Τα ποσοστά που εμφανίζονται παραπάνω καλύπτουν μόνο τον αριθμό γενεών και επαναλήψεων που έχουν εισαχθεί, επομένως θα ποικίλλουν ανάλογα με τον σπόρο (seed).
Μπορεί αυτό να προβλέψει τι θα κάνει ένας πραγματικός πληθυσμός;
Όχι. Δείχνει αποτελέσματα υπό ένα εσκεμμένα εξιδανικευμένο ουδέτερο μοντέλο Wright–Fisher [i, j]. Οι πραγματικές προβλέψεις απαιτούν βάσιμες εκτιμήσεις του δραστικού μεγέθους και οποιασδήποτε επιλογής, μετάλλαξης, μετανάστευσης, μεταβαλλόμενου μεγέθους πληθυσμού, επικαλυπτόμενων γενεών, σύνδεσης και δομής πληθυσμού που επηρεάζουν την περίπτωση [i, j].