Wright–Fisher मॉडल और जेनेटिक ड्रिफ्ट
पॉपुलेशन जेनेटिक ड्रिफ्ट सिम्युलेटर एक स्वतंत्र ऑनलाइन उपकरण है जो उपयोगकर्ताओं को द्विगुणित (diploid) Wright–Fisher मॉडल का उपयोग करके तटस्थ जेनेटिक ड्रिफ्ट (neutral genetic drift) का अनुकरण और विज़ुअलाइज़ करने की अनुमति देता है। एक सिम्युलेटेड प्रजनन समष्टि स्थापित करके, उपयोगकर्ता कई स्वतंत्र समष्टियों में पुनरुत्पादित होने वाले एलील-आवृत्ति पथों का अनुसरण और तुलना कर सकते हैं ताकि यह देखा जा सके कि पीढ़ियों के दौरान आनुवंशिक विविधता कैसे समाप्त या स्थिर (fixed) हो जाती है। यह उपकरण उपयोगकर्ताओं को यह समझने में मदद करता है कि रैंडम नमूनाकरण समय के साथ एलील आवृत्तियों को कैसे प्रभावित करता है, और यह वास्तविक दुनिया के विकासवादी परिणामों की भविष्यवाणी करने का दावा किए बिना सिम्युलेटेड विषमलुब्धता (heterozygosity), फिक्सेशन और विलोपन (loss) दरों को प्रदर्शित करता है।
यह मॉडल मानता है कि प्रत्येक पीढ़ी पिछली पीढ़ी से 2Ne एलील प्रतियों का एक रैंडम नमूना होती है। इस प्रक्रिया में, प्रत्येक द्विपद (binomial) ड्रॉ सटीक बर्नौली नमूनाकरण का उपयोग करता है; कोई सामान्य सन्निकटन (normal approximation) उपयोग नहीं किया जाता है। सभी गणनाएँ JavaScript डबल प्रिसिजन रखती हैं, जिससे आवृत्तियाँ 1 ÷ (2Ne) के गुणकों पर होती हैं और प्रदर्शन अधिकतम छह दशमलव स्थानों तक राउंड किया जाता है।
समष्टि सेटअप और इनपुट पैरामीटर
सिमुलेशन शुरू करने के लिए, उपयोगकर्ता "प्रभावी समष्टि आकार, Ne" और अन्य मापदंडों को कॉन्फ़िगर कर सकते हैं। इंटरफ़ेस में निम्नलिखित इनपुट विकल्प उपलब्ध हैं:
- परिदृश्य (Presets): उपयोगकर्ता त्वरित सेटअप के लिए "छोटी समष्टि", "बड़ी समष्टि", या "दुर्लभ एलील" जैसे पूर्व-निर्धारित परिदृश्यों का चयन कर सकते हैं।
- प्रभावी समष्टि आकार, Ne: इसे "प्रजनन करने वाले जीव" के रूप में दर्ज किया जाता है। यह 1 से 100,000 तक का एक पूर्णांक होना चाहिए।
- शुरुआती एलील आवृत्ति, p₀: इसे "अनुपात, 0–1" के रूप में दर्ज किया जाता है। यह मान दोनों सीमाओं सहित 0 और 1 के बीच होना चाहिए।
- पीढ़ियाँ: इसे "पीढ़ियाँ" के रूप में दर्ज किया जाता है। यह 1 से 2,000 तक का एक पूर्णांक होना चाहिए।
- स्वतंत्र समष्टियाँ: इसे "रेप्लिकेट्स" के रूप में दर्ज किया जाता है। यह 1 से 200 तक का एक पूर्णांक होना चाहिए।
- रैंडम सीड: सिम्युलेटेड पथों को दोबारा उत्पन्न करने के लिए 100 वर्णों या उससे कम का टेक्स्ट सीड दर्ज किया जा सकता है। उपयोगकर्ता रैंडम इनपुट के लिए "नया सीड" बटन पर भी क्लिक कर सकते हैं।
ब्राउज़र सुरक्षा और प्रदर्शन के लिए, सिस्टम अधिकतम 5,000,000 नियोजित एलील-कॉपी ड्रॉ की सीमा लागू करता है, जिसकी गणना 2Ne × पीढ़ियाँ × रेप्लिकेट्स के रूप में की जाती है। यदि यह उत्पाद सीमा से अधिक हो जाता है, तो सिमुलेशन ब्लॉक हो जाता है।
सिमुलेशन परिणाम और चार्ट विश्लेषण
सिमुलेशन चलाने के बाद, इंटरफ़ेस "एलील-आवृत्ति पथ" शीर्षक के तहत एक चार्ट प्रदर्शित करता है। इस चार्ट का ARIA लेबल "बार-बार सिम्युलेट की गई समष्टियों के लिए पीढ़ी दर पीढ़ी एलील आवृत्ति" है।
- X-अक्ष: "पीढ़ी" को दर्शाता है।
- Y-अक्ष: "एलील आवृत्ति p" को दर्शाता है।
- संदर्भ रेखा: "अपेक्षित p = p₀" पर एक स्थिर रेखा प्रदर्शित होती है।
- पथ संकेतक: चार्ट के नीचे "
‹total›में से‹shown›पथ दिखाए जा रहे हैं; प्रत्येक रेप्लिकेट को सारांश में शामिल किया गया है।" संदेश दिखाई देता है।
सिमुलेशन के दौरान निम्नलिखित स्थिति संदेश प्रदर्शित होते हैं:
- प्रारंभिक अवस्था: "समष्टि परिदृश्य बदलें या इस सीडेड उदाहरण को फिर से चलाएँ।"
- प्रसंस्करण के दौरान: "अगली पीढ़ियों का नमूना लिया जा रहा है…"
- पूरा होने पर: "सीड “
‹seed›” का उपयोग करके‹generations›पीढ़ियों के माध्यम से‹replicates›समष्टियों का अनुकरण किया गया।"
सांख्यिकीय सारांश और गणितीय सूत्र
सिम्युलेटर "सूत्र, प्रतिस्थापन और सटीकता" अनुभाग के तहत सैद्धांतिक और सिम्युलेटेड मूल्यों की तुलना करता है।
उपयोग किए जाने वाले मुख्य सूत्र:
- संक्रमण मॉडल: प्रत्येक पीढ़ी पर, Wright–Fisher संक्रमण अगली एलील संख्या का नमूना लेता है और इसे उपलब्ध जीन प्रतियों से विभाजित करता है: Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne)
- मोमेंट्स (Moments): नमूनाकरण से पहले, ये मोमेंट्स अगली पीढ़ियों की संभावित सीमा का वर्णन करते हैं: E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne)
- विषमलुब्धता (Heterozygosity) का ह्रास: असमान एलील जोड़ों के अपेक्षित हिस्से में t पीढ़ियों के दौरान इस प्रकार गिरावट आती है: H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ
पहली पीढ़ी के लिए प्रतिस्थापन उदाहरण:
- "पहली पीढ़ी 2 ×
‹population›=‹copies›एलील-कॉपी ड्रॉ का उपयोग करती है।" - "E[X₁] =
‹copies›×‹frequency›=‹mean›A प्रतियाँ।" - "पहली पीढ़ी के लिए, Var(X₁) =
‹copies›×‹frequency›× (1 −‹frequency›) =‹variance›; SD(p₁) =‹sd›।" - "पहले प्लॉट किए गए रेप्लिकेट ने X₁ =
‹sampled›A प्रतियों का नमूना लिया, इसलिए p₁ =‹sampled›÷‹copies›=‹frequency›।"
आउटपुट मेट्रिक्स तालिका:
| मीट्रिक लेबल | विवरण |
|---|---|
| p = 1 पर स्थिर (Fixed) | उन समष्टियों का अनुपात जहाँ एलील पूरी तरह से स्थिर हो गया है। |
| p = 0 पर विलुप्त (Lost) | उन समष्टियों का अनुपात जहाँ एलील पूरी तरह से लुप्त हो गया है। |
| T तक अवशोषित नहीं हुआ | उन समष्टियों का अनुपात जो अभी भी अंतिम पीढ़ी तक अपवाहित हो रही हैं। |
| माध्य अंतिम p | सभी रेप्लिकेट्स में अंतिम पीढ़ी पर औसत एलील आवृत्ति। |
| T पर अपेक्षित H | पीढ़ी T पर सैद्धांतिक रूप से बची हुई विषमलुब्धता। |
| T पर माध्य सिम्युलेटेड H | सिम्युलेटेड समष्टियों में पाई गई वास्तविक औसत विषमलुब्धता। |
| माध्य अवशोषण पीढ़ी | वह औसत पीढ़ी जिस पर एलील स्थिर या विलुप्त हो गया (यदि कोई भी समष्टि अवशोषित नहीं होती है, तो यह "नहीं पहुँचे" प्रदर्शित करता है)। |
विधि, मान्यताएँ और सीमाएँ
यह उपकरण Nature Education और Stanford Human Genetics पाठ्यपुस्तक द्वारा वर्णित सिद्धांतों पर आधारित है। इस मॉडल की कुछ महत्वपूर्ण सीमाएँ और मान्यताएँ हैं जिन्हें उपयोगकर्ताओं को ध्यान में रखना चाहिए:
- मान्यताएँ: निरंतर प्रभावी समष्टि आकार, असतत गैर-ओवरलैपिंग पीढ़ियाँ, रैंडम स्वतंत्र नमूनाकरण, चयनात्मक तटस्थता, एक अनलिंक्ड दो-एलील लोकस, और कोई उत्परिवर्तन या प्रवास नहीं।
- गैर-अनुप्रयुक्तता: जब चयन, समष्टि परिवर्तन, लिंग या प्रजनन-सफलता असंतुलन, ओवरलैपिंग पीढ़ियाँ, इनब्रीडिंग, लिंकेज, जीन प्रवाह, या समष्टि संरचना मायने रखती हो, तो इसे पूर्वानुमान के रूप में उपयोग न करें।
- सुरक्षा कारक: इस स्टोकेस्टिक जैविक मॉडल पर कोई इंजीनियरिंग सुरक्षा कारक लागू नहीं होता है।
- डेटा गोपनीयता: आपके पैरामीटर और सिम्युलेट की गई समष्टियाँ इसी ब्राउज़र में रहती हैं। सभी प्रसंस्करण स्थानीय रूप से उपयोगकर्ता के डिवाइस पर होते हैं और कोई भी डेटा बाहर नहीं भेजा जाता है।
अक्सर पूछे जाने वाले प्रश्न (FAQ)
प्रश्न: क्या मुझे वास्तविक गणना आकार दर्ज करना चाहिए या प्रभावी समष्टि आकार? उत्तर: प्रभावी समष्टि आकार, Ne दर्ज करें: एक आदर्श प्रजनन समष्टि का आकार जो उसी दर से अपवाहित (drift) होगा जिस दर से वह समष्टि जिसका आप अध्ययन करना चाहते हैं। जब लिंग अनुपात, प्रजनन सफलता, समष्टि में उतार-चढ़ाव, या समष्टि संरचना असमान हो, तो वास्तविक गणना आकार (census size) बहुत बड़ा हो सकता है।
प्रश्न: विषमलुब्धता गिरने के दौरान भी माध्य एलील आवृत्ति p₀ के करीब क्यों रह सकती है? उत्तर: कई स्वतंत्र समष्टियों में तटस्थ अपवाह (neutral drift) की कोई पसंदीदा दिशा नहीं होती है, इसलिए औसत आवृत्ति की अपेक्षा p₀ बनी रहती है। हालाँकि, व्यक्तिगत समष्टियाँ 0 या 1 की ओर फैलती हैं, और दोनों सीमाओं पर विषमलुब्धता (heterozygosity) शून्य होती है। इसलिए, समष्टि-व्यापी माध्य लगभग अपरिवर्तित रहने पर भी विविधता समाप्त हो जाती है।
प्रश्न: क्या अंतिम फिक्सेशन की तटस्थ प्रायिकता शुरुआती आवृत्ति के बराबर होती है? उत्तर: हाँ, बिना किसी उत्परिवर्तन (mutation) या प्रवास (migration) वाले इस तटस्थ मॉडल में, आवृत्ति p₀ से शुरू होने वाले एलील की अंतिम फिक्सेशन प्रायिकता p₀ और विलोपन (loss) प्रायिकता 1 − p₀ होती है। ऊपर दिखाई गई दरें केवल दर्ज की गई पीढ़ियों और दोहराव की संख्या को कवर करती हैं, इसलिए वे सीड के साथ बदलती रहेंगी।
प्रश्न: क्या यह अनुमान लगा सकता है कि एक वास्तविक समष्टि क्या करेगी? उत्तर: नहीं। यह जानबूझकर आदर्श बनाए गए एक तटस्थ Wright–Fisher मॉडल के तहत परिणाम दिखाता है। वास्तविक पूर्वानुमानों के लिए प्रभावी आकार और किसी भी चयन, उत्परिवर्तन, प्रवास, बदलते समष्टि आकार, ओवरलैपिंग पीढ़ियों, लिंकेज और समष्टि संरचना के विश्वसनीय अनुमानों की आवश्यकता होती है जो उस मामले के लिए मायने रखते हैं।