የሕዝብ ጀነቲክስ ፍሰት አስመሳይ

የሕዝብ ጀነቲክስ ፍሰት አስመሳይ

ይህ የሕዝብ ጀነቲክስ ፍሰት አስመሳይ በዲፕሎይድ Wright–Fisher ሞዴል ላይ በመመስረት ገለልተኛ የጄኔቲክ ፍሰትን (neutral genetic drift) ለመምሰል እና በምስል ለማሳየት የሚያገለግል ነፃ የመስመር ላይ መሣሪያ ነው። ተጠቃሚዎች አስመሳይ የመራቢያ ሕዝብን በማዘጋጀት፣ በበርካታ ገለልተኛ የሕዝብ ስብስቦች ላይ ሊባዙ የሚችሉ የአሌል ድግግሞሽ መንገዶችን መከተል እና ማነጻጸር ይችላሉ። ይህም በጊዜ ሂደት የጄኔቲክ ልዩነት እንዴት እንደሚጠፋ ወይም እንደሚረጋ ለማሳየት ይረዳል። መሣሪያው የዘፈቀደ ናሙና አወሳሰድ (random sampling) በጊዜ ሂደት በአሌል ድግግሞሾች ላይ እንዴት ተጽዕኖ እንደሚያሳድር ለመረዳት የሚያግዝ ሲሆን፣ አስመሳይ ሄትሮዚጎሲቲን፣ መርጋትን እና የመጥፋት መጠኖችን ያሳያል።

የእርስዎ መለኪያዎች እና አስመሳይ የሕዝብ ስብስቦች በዚህ ብሮውዘር ውስጥ ይቆያሉ። ሁሉም መለኪያዎች እና አስመሳይ የሕዝብ ስብስቦች በአካባቢው የሚሰሩ ሲሆን ሙሉ በሙሉ በተጠቃሚው ብሮውዘር ውስጥ ይቆያሉ። ምንም ዓይነት መረጃ ከተጠቃሚው መሣሪያ አይወጣም።

---

የሕዝብ ዝግጅት እና ግብዓቶች

በመሣሪያው በይነገጽ ውስጥ ተጠቃሚዎች የሚከተሉትን መለኪያዎች ማዘጋጀት ይችላሉ፦

• ሁኔታዎች (Presets)፦
  • አነስተኛ የሕዝብ ብዛት
  • ትልቅ የሕዝብ ብዛት
  • አልፎ አልፎ የሚገኝ አሌል
• ውጤታማ የሕዝብ ብዛት፣ Ne፦ በሚራቡ ግለሰቦች ቁጥር ይገለጻል። ከ1 እስከ 100,000 ያለ ሙሉ ቁጥር መሆን አለበት።
• የመጀመሪያ አሌል ድግግሞሽ፣ p₀፦ በድርሻ፣ 0–1 ይገለጻል። ከ0 እስከ 1 ባለው ክልል ውስጥ (ሁለቱንም ወሰኖች ጨምሮ) መሆን አለበት።
• ትውልዶች፦ በትውልዶች ብዛት ይገለጻል። ከ1 እስከ 2,000 ያለ ሙሉ ቁጥር መሆን አለበት።
• ገለልተኛ የሕዝብ ስብስቦች፦ በድግግሞሾች ብዛት ይገለጻል። ከ1 እስከ 200 ያለ ሙሉ ቁጥር መሆን አለበት።
• የዘፈቀደ ሲድ፦ አስመሳይ መንገዶችን መልሶ ለማምረት የሚያገለግል የጽሑፍ መስክ ነው። 100 ቁምፊዎች ወይም ከዚያ በታች መሆን አለበት። ተጠቃሚዎች አዲስ ሲድ የሚለውን ቁልፍ በመጫን ይህንን ግብዓት በዘፈቀደ መቀየር ይችላሉ።

---

አስመሳይ ውጤቶች እና ማሳያዎች

የፍሰት አስመሳይን አሂድ የሚለውን ቁልፍ ከተጫኑ በኋላ መሣሪያው የሚከተሉትን የሁኔታ መልዕክቶች፣ ገበታዎች እና ስታቲስቲክስ ያሳያል፦

የሁኔታ መልዕክቶች

• የሕዝቡን ሁኔታ ይቀይሩ ወይም ይህንን ሲድ የተደረገበትን ምሳሌ እንደገና ያሂዱ።
• የሚቀጥሉትን ትውልዶች ናሙና በመውሰድ ላይ…
• በሲድ “‹seed›” በመጠቀም ‹replicates› የሕዝብ ስብስቦችን በ‹generations› ትውልዶች ውስጥ አስመስሏል።

የአሌል ድግግሞሽ መንገዶች (ገበታ)

• የገበታው ARIA መለያ፦ ለተደጋገሙ አስመሳይ የሕዝብ ስብስቦች የአሌል ድግግሞሽ በትውልድ
• የአግድም ዘንግ (X-axis) መለያ፦ ትውልድ
• የአቀባዊ ዘንግ (Y-axis) መለያ፦ የአሌል ድግግሞሽ p
• የማጣቀሻ መስመር መለያ፦ የሚጠበቅ p = p₀
• የመንገዶች ብዛት አመልካች፦ ከ‹total› መንገዶች ውስጥ ‹shown› እያሳየ ነው፤ እያንዳንዱ ድግግሞሽ በማጠቃለያው ውስጥ ተካቷል።

የማጠቃለያ ስታቲስቲክስ

• በp = 1 ላይ የረጋ
• በp = 0 ላይ የጠፋ
• በT ያልተመጠጠ
• አማካይ የመጨረሻ p
• የሚጠበቅ H በT ላይ
• አማካይ አስመሳይ H በT ላይ
• አማካይ የመምጠጫ (absorption) ትውልድ (ውጤቱ ካልተገኘ አልተደረሰበትም የሚል ይበራል)

---

ፎርሙላ፣ መተካት እና ትክክለኛነት

መሣሪያው ስሌቶችን ለማከናወን እና ውጤቶችን ለማሳየት የሚከተሉትን የሂሳብ ቀመሮች ይጠቀማል፦

• በእያንዳንዱ ትውልድ፣ የ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)።
• የማይመሳሰሉ የአሌል ጥንዶች የሚጠበቀው ድርሻ በ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›።

---

ዘዴ፣ ግምቶች እና ገደቦች

ይህ ገጽ በNature Education እና በStanford Human Genetics ጽሑፍ የተገለጸውን ገለልተኛ ዲፕሎይድ Wright–Fisher ሞዴልን ተግባራዊ ያደርጋል። እያንዳንዱ አዲስ ትውልድ ከቀደመው የተወሰደ የ2Ne አሌል ቅጂዎች የዘፈቀደ ናሙና ነው። ምንጮች የተረጋገጡት ኦገስት 6፣ 2026 ነው።

የናሙና አወሳሰድ እና ትክክለኛነት

እያንዳንዱ ባይኖሚያል መውጣት ትክክለኛ የBernoulli ናሙናን ይጠቀማል፤ ምንም ዓይነት መደበኛ ግምት (normal approximation) ጥቅም ላይ አይውልም። ስሌቶች የJavaScript ባለ ሁለት ትክክለኛነትን (double precision) ይጠብቃሉ። ድግግሞሾች በ1 ÷ (2Ne) ጭማሪዎች ላይ ይገኛሉ፤ ማሳያው ቢበዛ እስከ ስድስት አስርዮሽ ቦታዎች ያጠጋጋል።

የሥራ ጫና ገደብ

ብሮውዘሩ ቢበዛ 5,000,000 የታቀዱ የአሌል ቅጂ መውጣቶችን ብቻ ይፈቅዳል፣ ይህም በ2Ne × ትውልዶች × ድግግሞሾች ይሰላል። ይህ ምርት ከገደቡ ካለፈ ስራው ይታገዳል።

ግምቶች እና ተፈጻሚነት የሌላቸው ሁኔታዎች

ሞዴሉ የሚከተሉትን ግምቶች ይወስዳል፦ ቋሚ ውጤታማ የሕዝብ ብዛት፣ ተደራራቢ ያልሆኑ የተነጠሉ ትውልዶች፣ የዘፈቀደ ገለልተኛ ናሙና፣ መራጭ ገለልተኝነት፣ አንድ ያልተገናኘ ባለ ሁለት አሌል ሎከስ፣ እና ምንም ሚውቴሽን ወይም ፍልሰት የለም። ምርጫ፣ የሕዝብ ለውጥ፣ የጾታ ወይም የመራባት ስኬት አለመመጣጠን፣ ተደራራቢ ትውልዶች፣ የቅርብ ዝምድና መራባት፣ ትስስር፣ የጂን ፍሰት ወይም የሕዝብ መዋቅር አስፈላጊ በሚሆኑበት ጊዜ ይህንን እንደ ትንበያ አይጠቀሙበት። የኢንጂነሪንግ ደህንነት ሁኔታ ለዚህ ስቶካስቲክ ባዮሎጂያዊ ሞዴል ተፈጻሚ አይሆንም።

---

የስህተት መልዕክቶች

በግብዓት ወቅት ስህተቶች ሲከሰቱ የሚከተሉት መልዕክቶች ይወጣሉ፦

• ግብዓቱ ቁጥር ካልሆነ፦ ‹field›፦ “‹token›” ቁጥር አይደለም።
• ሙሉ ቁጥር በሚያስፈልግበት ቦታ አስርዮሽ ከገባ፦ ‹field›፦ ሙሉ ቁጥር ያስገቡ።
• የሕዝብ ብዛት ከወሰን ውጭ ከሆነ፦ ውጤታማ የሕዝብ ብዛት ከ1 እስከ 100,000 ያለ ሙሉ ቁጥር መሆን አለበት።
• የመጀመሪያው ድግግሞሽ ከወሰን ውጭ ከሆነ፦ የመጀመሪያው የአሌል ድግግሞሽ በ0 እና በ1 መካከል መሆን አለበት (ሁለቱንም ወሰኖች ጨምሮ)።
• ትውልዶች ከወሰን ውጭ ከሆኑ፦ ትውልዶች ከ1 እስከ 2,000 ያለ ሙሉ ቁጥር መሆን አካባቸው።
• ገለልተኛ የሕዝብ ስብስቦች ከወሰን ውጭ ከሆኑ፦ ገለልተኛ የሕዝብ ስብስቦች ከ1 እስከ 200 ያለ ሙሉ ቁጥር መሆን አለባቸው።
• ሲድ ከጠፋ፦ አስመሳይነቱን መልሶ ማባዛት እንዲቻል የዘፈቀደ ሲድ (seed) ያስገቡ።
• ሲድ በጣም ረጅም ከሆነ፦ የዘፈቀደ ሲዱን በ100 ቁምፊዎች ወይም ከዚያ በታች ያድርጉት።
• የሥራ ጫና ገደብ ካለፈ፦ ይህ ዝግጅት ‹draws› የአሌል ቅጂ መውጣቶችን ያቅዳል፣ ይህም ከብሮውዘሩ ገደብ ‹limit› በላይ ነው። የሕዝቡን ብዛት፣ ትውልዶችን ወይም ድግግሞሾችን ይቀንሱ።

---

ተደጋግመው የሚነሱ ጥያቄዎች (FAQ)

ጥያቄ፦ ማስገባት ያለብኝ የቆጠራ መጠንን ነው ወይስ ውጤታማ የሕዝብ ብዛትን?
መልስ፦ ውጤታማ የሕዝብ ብዛት Ne ያስገቡ፦ ይህ ልታጠኑት ከፈለጉት የሕዝብ ስብስብ ጋር በእኩል መጠን ሊፈስ የሚችል ሃሳባዊ የመራቢያ ሕዝብ መጠን ነው። የጾታ ጥምርታ፣ የመራባት ስኬት፣ የሕዝብ መዋዠቅ ወይም የሕዝብ መዋቅር እኩል ባልሆኑበት ጊዜ የቆጠራው መጠን በጣም የላቀ ሊሆን ይችላል።

ጥያቄ፦ ሄትሮዚጎሲቲ በሚቀንስበት ጊዜ አማካይ የአሌል ድግግሞሽ በp₀ አቅራቢያ ሊቆይ የሚችለው ለምንድን ነው?
መልስ፦ ገለልተኛ ፍሰት (neutral drift) በብዙ ገለልተኛ የሕዝብ ስብስቦች ላይ የተለየ አቅጣጫ የለውም፣ ስለዚህ አማካይ ድግግሞሹ በሚጠበቀው መሠረት p₀ ሆኖ ይቆያል። ሆኖም ግን፣ ነጠላ የሕዝብ ስብስቦች ወደ 0 ወይም 1 ይሰራጫሉ፣ እና ሁለቱም ወሰኖች ዜሮ ሄትሮዚጎሲቲ አላቸው። ስለዚህ በሕዝቦች መካከል ያለው አማካይ ሳይለወጥ በሚቆይበት ጊዜም እንኳ ልዩነቱ ይጠፋል።

ጥያቄ፦ የገለልተኛ የመጨረሻ መርጋት ዕድል ከመጀመሪያው ድግግሞሽ ጋር እኩል ነው?
መልስ፦ አዎ፣ ምንም ዓይነት ሚውቴሽን ወይም ፍልሰት በሌለበት በዚህ ገለልተኛ ሞዴል ውስጥ፣ በድግግሞሽ p₀ የሚጀምር አሌል የመጨረሻ የመርጋት (fixation) ዕድሉ p₀ ሲሆን የመጥፋት ዕድሉ ደግሞ 1 − p₀ ነው። ከላይ የሚታዩት መጠኖች ግን የገቡትን የትውልዶች እና የድግግሞሾች ብዛት ብቻ የሚሸፍኑ በመሆናቸው እንደ ሲዱ ይለያያሉ።

ጥያቄ፦ ይህ እውነተኛ የሕዝብ ስብስብ የሚያደርገውን መተንበይ ይችላል?
መልስ፦ አይደለም። ይህ ሆን ተብሎ ሃሳባዊ በተደረገ ገለልተኛ የWright–Fisher ሞዴል ስር ያሉ ውጤቶችን ያሳያል። እውነተኛ ትንበያዎች ለጉዳዩ አስፈላጊ የሆኑትን የውጤታማ መጠን እና የማንኛውንም ምርጫ፣ ሚውቴሽን፣ ፍልሰት፣ ተለዋዋጭ የሕዝብ ብዛት፣ ተደራራቢ ትውልዶች፣ ትስስር እና የሕዝብ መዋቅር ተገቢ ግምቶችን ይፈልጋሉ።

---

ማጣቀሻዎች

• Nature Education, Genetic Drift and Effective Population Size
• Stanford Human Genetics, Genetic Drift: What happens to alleles over time