آموزش آنالیز ژنوم در برنامه نویسی بایوپایتون

4 سال پیش

امتیاز دهید post

آموزش آنالیز ژنوم در برنامه نویسی بایوپایتون

در این درس از آموزش های برنامه نویسی سایت سورس باران، ما در مورد آموزش آنالیز ژنوم در برنامه نویسی بایوپایتون بحث خواهیم کرد.

پیشنهاد ویژه : پکیج آموزش طراحی وب سایت با پایتون

یک ژنوم مجموعه کاملی از DNA می باشد که شامل همه ژن های آن است. تجزیه و تحلیل ژنوم به مطالعه ژن های فردی و نقش آنها در وراثت اشاره دارد.

نمودار ژنوم

نمودار ژنوم اطلاعات ژنتیکی را به عنوان نمودار نشان می دهد. بایوپایتون از ماژول Bio.Graphics.GenomeDiagram برای نشان دادن GenomeDiagram استفاده می کند. ماژول GenomeDiagram نیاز به نصب ReportLab دارد.

مراحل ایجاد نمودار

روند ایجاد نمودار به طور کلی از الگوی ساده زیر پیروی می کند:

برای هر مجموعه جداگانه از ویژگی هایی که می خواهید نمایش دهی ، یک FeatureSet ایجاد کنید و اشیا objects Bio.SeqFeature را به آنها اضافه کنید.
برای هر گرافی که می خواهید نمایش دهید یک GraphSet ایجاد کنید و داده های نمودار را به آنها اضافه کنید.
برای هر مسیری که می خواهید روی نمودار نمودار ایجاد کنید و GraphSets و FeatureSets را به مسیر های مورد نیاز خود اضافه کنید.
یک نمودار ایجاد کنید و مسیر ها را به آن اضافه کنید.
به نمودار دستور ترسیم تصویر را بدهید
تصویر را در یک فایل بنویسید.

اجازه دهید مثالی از فایل GenBank ورودی بیاوریم –

https://raw.githubusercontent.com/biopython/biopython/master/Doc/examples/ls_orchid.gbk و سوابق را از شی SeqRecord بخوانید و سپس در نهایت یک نمودار ژنوم ترسیم کنید. در زیر توضیح داده شده است:

ما ابتدا همه ماژول ها را مانند شکل زیر وارد می کنیم –

&gt;&gt;&gt; from reportlab.lib import colors 
&gt;&gt;&gt; from reportlab.lib.units import cm 
&gt;&gt;&gt; from Bio.Graphics import GenomeDiagram

>>> from reportlab.lib import colors

>>> from reportlab.lib.units import cm

>>> from Bio.Graphics import GenomeDiagram

اکنون، برای خواندن داده ها، ماژول SeqIO را وارد کنید:

&gt;&gt;&gt; from Bio import SeqIO 
record = SeqIO.read("example.gb", "genbank")

1 2	>>> from Bio import SeqIO record = SeqIO.read("example.gb", "genbank")

در اینجا، رکورد توالی فایل genbank را می خواند.

اکنون، یک نمودار خالی برای اضافه کردن مجموعه مسیرها و ویژگی ها ایجاد کنید:

&gt;&gt;&gt; diagram = GenomeDiagram.Diagram(
   "Yersinia pestis biovar Microtus plasmid pPCP1") 
&gt;&gt;&gt; track = diagram.new_track(1, name="Annotated Features") 
&gt;&gt;&gt; feature = track.new_set()

>>> diagram = GenomeDiagram.Diagram(

"Yersinia pestis biovar Microtus plasmid pPCP1")

>>> track = diagram.new_track(1, name="Annotated Features")

>>> feature = track.new_set()

اکنون، می توانیم تغییرات تم رنگ را با استفاده از رنگ های جایگزین از سبز به خاکستری به شرح زیر اعمال کنیم –

&gt;&gt;&gt; for feature in record.features: 
&gt;&gt;&gt;    if feature.type != "gene": 
&gt;&gt;&gt;       continue 
&gt;&gt;&gt;    if len(feature) % 2 == 0: 
&gt;&gt;&gt;       color = colors.blue 
&gt;&gt;&gt;    else: 
&gt;&gt;&gt;       color = colors.red 
&gt;&gt;&gt; 
&gt;&gt;&gt;    feature.add_feature(feature, color=color, label=True)

>>> for feature in record.features:

>>> if feature.type != "gene":

>>> continue

>>> if len(feature) % 2 == 0:

>>> color = colors.blue

>>> else:

>>> color = colors.red

>>>

>>> feature.add_feature(feature, color=color, label=True)

اکنون می توانید پاسخ زیر را بر روی صفحه خود مشاهده کنید:

&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d3dc90&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d3dfd0&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x1007627d0&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57290&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57050&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57390&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57590&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57410&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57490&gt; 
&lt;Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d574d0&gt;

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d3dc90>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d3dfd0>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x1007627d0>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57290>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57050>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57390>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57590>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57410>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d57490>

<Bio.Graphics.GenomeDiagram._Feature.Feature object at 0x105d574d0>

اجازه دهید ما یک نمودار برای رکورد ورودی بالا ترسیم کنیم –

&gt;&gt;&gt; diagram.draw(
   format = "linear", orientation = "landscape", pagesize = 'A4', 
   ... fragments = 4, start = 0, end = len(record)) 
&gt;&gt;&gt; diagram.write("orchid.pdf", "PDF") 
&gt;&gt;&gt; diagram.write("orchid.eps", "EPS") 
&gt;&gt;&gt; diagram.write("orchid.svg", "SVG") 
&gt;&gt;&gt; diagram.write("orchid.png", "PNG")

>>> diagram.draw(

format = "linear", orientation = "landscape", pagesize = 'A4',

... fragments = 4, start = 0, end = len(record))

>>> diagram.write("orchid.pdf", "PDF")

>>> diagram.write("orchid.eps", "EPS")

>>> diagram.write("orchid.svg", "SVG")

>>> diagram.write("orchid.png", "PNG")

پس از اجرای دستور بالا، می توانید تصویر زیر را در دایرکتوری بایوپایتون خود ذخیره کنید.

** Result **
genome.png

1 2	Result genome.png

Creating Diagram

با انجام تغییرات زیر می توانید تصویر را به صورت دایره ای ترسیم کنید:

&gt;&gt;&gt; diagram.draw(
   format = "circular", circular = True, pagesize = (20*cm,20*cm), 
   ... start = 0, end = len(record), circle_core = 0.7) 
&gt;&gt;&gt; diagram.write("circular.pdf", "PDF")

>>> diagram.draw(

format = "circular", circular = True, pagesize = (20*cm,20*cm),

... start = 0, end = len(record), circle_core = 0.7)

>>> diagram.write("circular.pdf", "PDF")

بررسی اجمالی کروموزوم ها

مولکول DNA در ساختارهای رشته ای مانند کروموزوم ها بسته بندی می شود. هر کروموزوم از DNA ساخته شده است که بارها و بارها در اطراف پروتئین هایی به نام هیستون هایی که ساختار آن را پشتیبانی می کنند، پیچیده است.

کروموزوم ها در هسته سلول قابل مشاهده نیستند – حتی در زیر میکروسکوپ – هنگامی که سلول تقسیم نمی شود. با این حال، DNA سازنده کروموزوم ها در طول تقسیم سلول محکم بسته می شود و سپس در زیر میکروسکوپ قابل مشاهده است.

در انسان، هر سلول به طور معمول حاوی 23 جفت کروموزوم است که در مجموع 46 عدد است. 22 جفت از این جفت ها که اتوزوم نامیده می شوند، در مردان و زنان یکسان به نظر می رسند. جفت 23 ام، کروموزوم های جنسی، بین زن و مرد متفاوت است. ماده ها دو نسخه از کروموزوم X دارند، در حالی که مردان دارای یک کروموزوم X و یک کروموزوم Y هستند

منبع.

لیست جلسات قبل آموزش برنامه نویسی بایوپایتون

امتیاز دهید post

نویسنده مطلب saber

آموزش آنالیز ژنوم در برنامه نویسی بایوپایتون

آموزش آنالیز ژنوم در برنامه نویسی بایوپایتون

نمودار ژنوم

مراحل ایجاد نمودار

بررسی اجمالی کروموزوم ها

لیست جلسات قبل آموزش برنامه نویسی بایوپایتون

دیدگاه شما

بدون دیدگاه

آموزش آنالیز ژنوم در برنامه نویسی بایوپایتون

نمودار ژنوم

مراحل ایجاد نمودار

بررسی اجمالی کروموزوم ها

لیست جلسات قبل آموزش برنامه نویسی بایوپایتون

دیدگاه شما

برای صرف‌نظر کردن از پاسخ‌گویی اینجا را کلیک نمایید.

بدون دیدگاه