Genetic Diversity of Persian Crocodile Crocodylus Palustris Using Sequencing of D-Loop and Cyt b Regions of Mitochondria

Document Type : Research Paper

Authors

1 PhD Student, Animal Science Department, Faculty of Agriculture, Ferdowsi University of Mashhad, Mashhad, Iran.

2 Professor, Animal Science Department, Faculty of Agriculture, Ferdowsi University of Mashhad. Mashhad, Iran.

3 Associate Professor, School of Molecular and Biomedical Science, The University of Adelaide, Australia.

4 Assistant Professor, Animal Science Department, Faculty of Agriculture, Ferdowsi University of Mashhad, Mashhad, Iran.

Abstract

Identification of the genetic features of Persian mugger crocodile, the Gando, is crucial in management planning for conservation and preservation of this endangered species. Therefore, this study is aimed to investigate the genetic characteristics of Gando in Iran, based on population studies of the conserved and variable regions of the mitochondrial DNA. For this purpose, 12 samples of skin tissue and blood were collected from three regions of Chabahar city, south east of Iran. After DNA extraction, conserved and variable regions of the mitochondrial genome Cyt b and D-loop were amplified using PCR. The PCR products were sequenced and results were analyzed using bioinformatics software based on nucleotide data as well as protein codons. Results revealed that all samples had similar nucleotide sequences and no genetic variation was found among the samples. The closest species to Gando were Crocodylus palustris and then C. porosus. Comparison of sequences with other Crocodile species revealed the presence of 12 polymorphic loci among 1156 loci studied and also 5 haplogroups were introduced. As we didn’t find any nucleotide diversity among studied samples, therefore, to increase the genetic diversity of the Gando, it is essential to include specific conservational plans for this endangered species.

Keywords


بررسی تنوع ژنتیکی گاندو، کروکودیل پوزه کوتاه ایرانی (Crocodylus Palustris)، با استفاده از توالی یابی نواحی D-Loop و Cyt b میتوکندری

 

احمد افشاریان1، محمدرضا نصیری2*، اسماعیل ابراهیمی3، علی جوادمنش4

1 دانشجوی مقطع دکتری، گروه علوم دامی، دانشگاه فردوسی مشهد، مشهد، ایران.

2 استاد، گروه علوم دامی، دانشکده کشاورزی، دانشگاه فردوسی مشهد، مشهد، ایران.

3 دانشیار، گروه بیومدیکال مولکولی، دانشگاه آدلاید استرالیا.

4 استادیار، گروه علوم دامی، دانشکده کشاورزی، دانشگاه فردوسی مشهد، مشهد، ایران.

 

تاریخ دریافت: 19/05/1396، تاریخ پذیرش: 14/07/1396

چکیده

شناسایی ویژگی‌های ژنتیکی تمساح پوزه کوتاه ایرانی (گاندو)، در برنامه‌ریزی مدیریتی به منظور حفاظت و نگهداری آن نقش بسیار مهمی دارد، لذا این پژوهش با هدف، بررسی ژنتیکی کروکودیل پوزه کوتاه ایران بر اساس مطالعه جمعیتی نواحی ثابت و متغیر ژنوم میتوکندری انجام گردید. بدین منظور، 12 نمونه بافت پوست و خون از سه منطقه مختلف شهرستان چابهار جمع‌آوری و پس از استخراج DNA، نواحی ثابت و متغییر ژنوم میتوکندری Cyt b و D-loop با استفاده از روش PCR تکثیر گردید. محصولات PCR توالی یابی و سپس آنالیز نتایج با استفاده از نرم‌افزارهای بیوانفورماتیک و مبتنی بر داده‌های نوکلئوتیدی و کدون‌های پروتئینی انجام گرفت. نتایج این تحقیق نشان داد که تمامی نمونه‌های مورد مطالعه، دارای توالی نوکلوتیدی یکسان هستند و هیچ تنوع ژنتیکی در بین این نمونه‌ها یافت نشد. نزدیکترین گونه به گاندو، گونه Crocodylus palustris و سپس گونه C. porosus بدست آمد. مقایسه توالی‌ها نشان دهنده وجود 12 جایگاه پلی‌مورفیک در بین 1156 جایگاه مورد مطالعه بود و 5 هاپلوگروپ معرفی گردید. لذا، با توجه به عدم مشاهده تنوع ژنتیکی در بین نمونه‌های مورد مطالعه، استفاده از برنامه‌های حفاظتی جهت افزایش تنوع ژنتیکی در گاندو امری بسیار ضروری به نظر می­رسد.

کلمات کلیدی: تمساح گاندو، تنوع ژنتیکی، درخت فیلوژنتیک، ژنوم میتوکندری.

 

 


 

مقدمه

 

جنس کروکودیل شامل تعداد زیادی از گونه‌های کورکودیل می باشد که ساکن مناطق گرمسیری و نیمه گرمسیری نواحی آسیا، استرالیا، آفریقا و نیوزلند می­باشد (Brochu, 2003; McAliley et al, 2006). کروکودیل مردابی یا پوزه کوتاه (گاندو) با نام علمی Crocodylus palustris  اندازه متوسط از جنس کروکودیل است (بیشترین طول قد 4-5 متر) که از ایران به سمت شرق تا آسام هند و به جنوب تا سریلانکا پراکنش دارد و حوزة پراکنش آن در ایران غربی ترین نقطة پراکنش جهانی این گونه به حساب می‌آید (شکل 1). کروکودیل مردابی دارای پهن ترین پوزه در بین تمام اعضای جنس کروکودیلوس (Crocodylus) است و بنابراین بیشترین شباهت ظاهری را به آلیگاتورها دارد (Da Silva & Lenin, 2010).  اگرچه گونه C. palustris ساکن مرداب‌ها می‌باشند، اما بررسی ها نشان داده که این گونه بیشتر ساکن رودخانه‌ها، دریاچه‌ها و برکه‌ها بوده و به خوبی توانسته در برکه‌های مصنوعی و دست ساز، کانالهای آبیاری، دریاچة پشت سد، برکه‌های طبیعی نواحی هند، پاکستان، سریلانکا و ایران سازگار شود (Whitaker & Andrews 2003; Da Silva & Lenin, 2010).

 

 

شکل 1- A: پراکنش کروکودیل گونه C. palustris براساس مطالعه Whitaker and Andrews (2003). B: گاندو در ایستگاه تحقیقات زیست‌-محیطی ریکوکش چابهار، ایران.

Figure 1- A: Distribution of C. palustris crocodile according to Whitaker and Andrews Study (2003). B: Gando in the Environmental Research Station of Rikukash, Chabahar, Iran.

 

 

اطلاع از ساختار ژنتیکی جمعیت­ها می­تواند کمک بزرگی برای برنامه ریزی برای طرح­های اصلاح نژادی و از همه مهمتر، حفظ ذخایر ژنتیکی باشد. روش­های مولکولی و استفاده از نشانگرهای مولکولی در این زمینه یکی از بهترین گزینه‌ها به حساب می­آید، زیرا با توجه به اطلاعات زیادی که به دست می­دهد می­تواند نتایجی که از تجزیه و تحلیل رکوردها با روش­های آماری به دست آمده است را تأیید و تکمیل نموده و حتی ممکن است که آنها را رد کند (Naderi et al., 2007). به علاوه، استفاده از ژنتیک ملکولی فواید زیادی دارد که یکی از این فواید معنی دار تعیین ژنوتیپ افراد برای جایگاه خاصی است (Mousavizadeh et al., 2009) همچنین استفاده از نشانگرهای ژنتیکی در انتخاب و اصلاح نژاد حیوانات ممکن است به طور مهیجی پیشرفت ژنتیکی را تسریع کند (Javanmard et al., 2008) و مطالعه تنوع ژنتیکی نژادهای بومی برای حفاظت از منابع ژنتیکی ذخایر بومی لازم و ضروری است (Mohammadi et al., 2009). حفاظت باید بر اساس دانش عمیقی از منابع ژنتیکی نژادهای خاص باشد، لذا تلاش برای شناسایی و تعیین خصوصیات ژنتیکی نژادهای بومی و محلی بسیار اهمیت دارد (Zamani et al., 2013). مطالعات متنوعی با هدف طبقه بندی ژنتیکی جمعیت‌ها و گونه­های نزدیک به هم، بررسی امکان انشقاق گونه‌های مختلف از یک جد مشترک، مطالعه رابطه فیلوژنی هر موجود با سایر گونه­ها و نژادها و دستیابی به راهکارهایی برای حفظ ذخایر ژنتیکی گونه‌های مختلف کروکودیل، در کشورهای مختلف جهان انجام شده است (Brochu, 2000; Brochu et al., 2010; Brochu and Densmore, 2001; Gatesy et al., 2004;McAliley et al., 2006; Willis et al., 2007; Gatesy and Amato, 2008; Meganathan et al., 2010). این مطالعات به برنامه ریزی صحیح جهت حفظ بقای این گونه که در طبقه بندی کتاب سرخ اتحادیه جهانی حفاظت International Union for the conservation (IUCN) در طبقه آسیب پذیر قرار گرفته است، کمک شایانی می­نماید (Da Silva & Lenin, 2010 ). امروزه استفاده از روش‌های مولکولی مبتنی بر DNA در مطالعات تنوع و ساختار جمعیت­های کروکودیل در اولویت قرار گرفته است. در این بین توالی­یابی ژن­های میتوکندری یکی از کاربردی­ترین روش‌ها برای تعیین رابط فیلوژنتیکی بین جمعیت­ها و گونه­های نزدیک به هم محسوب می­شود (Bruford et al., 2003; Dolati et al., 2016). همچنین DNA میتوکندری نشانگر بسیار مناسبی برای تشخیص بین گونه­ای گروه‌هایی است که برای 100 تا 10000 سال از هم جدا بوده­اند و نیز به دلیل سرعت بالای تکامل DNA میتوکندری در مقایسه با DNA هسته­ای باعث شده که از آن برای مطالعات تکاملی استفاده گردد (Nassiri et al., 2010). DNA میتوکندریایی (DNA میتوکندری) در تشخیص گونه­ها و ترسیم روابط فیلوژنتیکی دارای مزایایی از جمله تعداد زیادی نسخه به ازای هر سلول، اندازه کوچکتر آن از DNA ژنومی، وراثت پذیری مادری، هاپلوئید بودن، عدم وجود نوترکیبی در آن­ها، وجود نواحی حفاظت شده و وجود نواحی حفاظت نشده­ای مانند ناحیه D-Loop برای مطالعات تکاملی گونه‌های وابسته می‌باشد (Bellagamba et al., 2001). شناسایی توالی ژن­های میتوکندری به عنوان مشخصه ژنتیکی ثابت محسوب می­شود و به تهیه شناسنامه ژنتیکی، میزان تنوع موجود در جمعیت، امکان تشخیص گونه­ها و نژاد­ها، نگهداری خلوص نژادهای بومی و به حفظ گونه­های بومی از خطر انقراض و اختلاط ژنتیکی با سایر نژاد­ها کمک می‌کند (Janke et al., 2005; Roos et al., 2007, Feng et al., 2010; Meganathan et al., 2010). بنابراین مطالعه ژنتیکی و فیلوژنتیکی کروکودیل پوزه کوتاه استان سیستان و بلوچستان به منظور دستیابی به اطلاعاتی از قبیل تنوع گونه‌ها و زیرگونه‌ها، اطلاعات در مورد منشأ اهلی شدن مبنای انجام این پژوهش قرار گرفت. در این پژوهش برای اولین بار در ایران، کروکودیل پوزه کوتاه ایران در منطقه چابهار از دیدگاه ژنتیکی، محاسبه پارامترهای ژنتیکی و تعیین نقشه میتوکندریایی آن بر اساس مطالعه جمعیتی نواحی ثابت و متغیر ژنوم میتوکندری که در معرض خطر انقراض می باشد، گزارش می‌شود.

 

مواد و روش ها

با هماهنگی و همکاری اداره محیط زیست منطقه چابهار (مختصات: 25° 17′ 31″ N, 60° 38′ 35″ E)، 12 نمونه بافت از بخش های سطحی پوست و خون از کروکودیل پوزه کوتاه ایرانی، از 4 منطقه این شهرستان جمع آوری و در داخل اتانول مطلق نگهداری شد. از آنجا که نمونه‌های بافت، دارای بافتی به شدت متراکم و دیر هضم بودند، مقدار 1/0 گرم از هر نمونه با استفاده از تیغ اسکاپل جدا شده و سپس با استفاده از نیتروژن مایع هموژن گردید. سپس استخراج DNA با انجام تغییراتی در کیت # K0721  GeneJET Genomic DNA Purification Kit  شرکت Thermo انجام شد. این تغییرات شامل افزایش مدت زمان هضم بافت با استفاده از پروتئیناز K به مدت 16 ساعت بود. مابقی مراحل با توجه به دستورالعمل کیت صورت گرفت. یک ناحیه متغیر و یک ناحیه حفاظت شده از میتوکندری به منظور بررسی پارامترهای ژنتیکی جمعیت مورد بررسی مورد مطالعه قرار گرفت. ناحیه مد نظر برای بخش ثابت کل ناحیه Cyt b و ناحیه مورد نظر برای بخش متغیر، ناحیه‌ای از منطقه کنترل ژنوم میتوکندری حاوی بخش D-Loop (Displacement Loop) بود. بدین منظور توالی میتوکندریایی مربوط به میتوکندری گونه Crocodylus palustris از بانک ژنی (Benson et al., 2013) دریافت شد (Accession Number: HM488007). همچنین به منظور افزایش ضریب دقت توالی مربوط به گونه های دیگر کروکودیل متعلق به جنس Crocodylus نیز دریافت شدند تا بر اساس نواحی حفاظت شده بین گونه ای طراحی پرایمر صورت بگیرد (جدول 1) (شکل 2). طراحی پرایمر بوسیله برنامه Primer Premier v.5 (Lalitha, 2000) صورت پذیرفت (جدول 2).


 

جدول 1- گونه های کروکودیل مورد مطالعه در طراحی پرایمر.

Table 1- Crocodile species that used to design primer pair.

 

گونه   Species

شماره دسترسی   accession number

کروکودیلوس آکیتس

Crocodylus acutus

JF502241

کروکودیلوس اینترمدیس

Crocodylus intermedius

JF502242

کروکودیلوس جانسونی

Crocodylus johnsoni

HM488008

کروکودیلوس مورلتی

Crocodylus moreletii

HQ585889

کروکودیلوس تایلوتیکس

Crocodylus niloticus

AJ810452

کروکودیلوس نوایژنی

Crocodylus novaeguineae

JF502240

کروکودیلوس پالستوریس

Crocodylus palustris

HM488007, FJ173278, FJ173279, FJ173280, FJ173281, FJ173282, FJ173283, FJ173284, FJ173285, FJ173286, FJ173276, FJ173277,

کروکویلوس پایروسس

Crocodylus porosus

AJ810453

کروکودیلوس رومبایفر

Crocodylus rhombifer

JF502247

کروکودیلوس سایمنسس

Crocodylus siamensis

DQ353946

 

 

 

شکل 2- تصویر شماتیک منطقه طراحی پرایمر برای تکثیر ناحیه Cyt b و بخشی از ناحیه کنترل میتوکندری کروکودیلوس پالوستریس.

Figure 2- Schematic showing the primer design of Cyt b and D-loop region of C. palustritis.

 

 

 

جدول 2- مشخصات پرایمرهای استفاده شده برای تکثیر ناحیه Cytbو بخشی از D-loop میتوکندری کروکودیلوس پالوستریس.

Table 2- The primer characteristics used for amplifications of Cyt b and D-loop regions of C. palustritis mitochondria.

نام پرایمر

Primer Name

توالی

Sequence

طول (جفت باز)

Length (bp)

∆G [-kcals/mol]

محتوای گوانین-سیتوزین

GC Content

دمای ذوب

Melting Temperature

اندازه محصول (جفت باز)

Product Size (bp)

A.F

CCACACTAATTATCATTAAGCA

22

36

32

52

710

A.R

CTGTGAGGATAAACGGGAG

19

35

52

B.F

AGGAGGACCATCCGTCAACAG

21

40

57

61

 

649

B.R

CCGAAGGGGTTTGAATAGTGTG

22

42

50

C.F

GAATCGGAGGACAGCCAGTA

20

38

55

58

846

C.R

CGGGTAGTGGATGAGATCTCC

21

39

57

 


 

واکنش زنجیره ای پلیمراز توسط دستگاه ترموسایکلر Biometra مدل T personal در 35 چرخه حرارتی و با حجم نهایی 50 میکرولیتر انجام شد. از کیت مسترمیکس PCR شرکت آمپلیکون (دانمارک) با غلظت MgCl2، 5/1 میلی مولار جهت تکثیر قطعات هدف استفاده گردید. به منظور تایید تکثیر ناحیه مورد نظر، طی واکنش PCR، الکتروفورز محصولات PCR بر روی ژل آگارز 5/1 درصد با رنگ آمیزی GelRed صورت گرفت. سپس 50 میکرولیتر از محصول PCR با استفاده از کیت GeneJET PCR Purification Kit  شرکت Thermo (امریکا) خالص سازی شد و به همراه 50 میکرولیتر از آغازگرهای رفت و برگشت با غلظت 10 پیکومول به منظور توالی یابی به شرکت ماکروژن کره جنوبی ارسال گردید. این نمونه‌ها با استفاده از دستگاه ABI 3130 به روش اتوماتیک سانگر توالی یابی شدند. داده‌های حاصل از توالی یابی ژنوم میتوکندریایی کروکودیل پوزه کوتاه پس از دریافت به منظور تصحیح خطاهای توالی یابی مورد بازبینی قرار گرفتند. بررسی اولیه توالی‌های بدست آمده با استفاده از برنامه Chromas Lite 2.01 (Mccarthy, 2005) انجام شد. جهت تعیین بالاترین همولوژی کروکودیل پوزه کوتاه از رویه بلاست نوکلئوتیدی تحت پایگاه NCBI (Sayers et al., 2011) استفاده شد. توالی مورد توافق (Consensus) بدست آمده بر اساس معیار اکثریت، با استفاده از نرم‌افزار (Swindell & Plasterer, 1997) SEQMAN تعیین و این توالی توسط ابزار BankIt پایگاه اطلاعاتی NCBI (Sayers et al., 2011) ثبت گردید. بررسی تنوع ژنتیکی در سطح نوکلئوتیدی به کمک ابزار آنلاین DIVEIN)) (Deng et al., 2010) و Arlequin suite ver 3.5 (Excoffier & Lischer, 2010) صورت پذیرفت. به منظور مقایسه داده‌های این مطالعه با دیگر مطالعات از توالی‌های مرجع ثبت شده در بانک داده NCBI موسوم به توالی های مرجع استفاده شد (Pruitt et al., 2007). بدین منظور توالی مایتوژنومیکس 13 گونه مربوط به جنس کروکودیل و همچنین توالی‌های مربوط به جنس Osteolaemus دریافت شد. همچنین از دو زیرخانواده Alligatorinae و Gavialinae  نیز به عنوان کلاستر‌های خواهری استفاده شد. همردیف سازی در سطح نوکلئوتیدی به کمک الگوریتم Clustal نسخه 2 (Larkin et al., 2007) صورت گرفت. رسم درخت فیلوژنتیکی به کمک روش‌های مختلف بنا به شرایط، شامل اتصال نزدیک[1] (Saitou & Nei., 1987)، امساک حداکثر[2] (Farris, 1970; Fitch, 1971)، روش حداکثر درست نمایی (Felsenstein, 1981) و استنباط بیزی بهبود یافته به کمک زنجیره مارکوف مونت کارلو (Yang & Rannala, 1997) بود. ابزار مورد نظر در آنالیز اتصال نزدیک و امساک حداکثر برنامه تجاری PAUP نسخه 4 (Swofford, 2003) و برای دو آنالیز حداکثر درست نمایی و آنالیز بیزی به ترتیب RAxML (Stamatakis, 2006) و PhyloBayes MPI   (Lartillot et al., 2013) انجام گرفت. به منظور ارزیابی و دستیابی به پارامترهای بهینه ، روش بوت استرپ (Stamatakis et al., 2008) و (Felsenstein, 1985) 100 چندگانه مکرر انجام شد. همچنین از شاخص‌های امتیاز حداکثر درست نمایی مقیاس شده و شاخص پایداری مقیاس شده (RC) (Farris, 1989) نیز استفاده شد. فاصله ژنتیکی بر اساس توالی نوکلئوتیدی و همچنین توالی کدکننده پروتئینی مربوطه در ناحیه Cyt b انجام پذیرفت. روش مورد استفاده روش Maximum Composite Likelihood بود و در محاسبه فاصله ژنتیکی، سه جایگاه هر کدون، به شکل مشترک مد نظر قرار گرفت.

 

نتایج و بحث

استخراج DNA از بافت در تمامی نمونه‌ها با موفقیت انجام شد. DNA‌های استخراج شده فاقد شکستگی می‌باشند که نشان دهنده کیفیت مناسب استخراج DNA و حداقل آسیب وارده به آن طی فرایند آماده سازی و استخراج می باشد (شکل 3).

فرآیند تکثیر قطعات مربوط به ناحیه Cyt b و ناحیه کنترل با موفقیت انجام گرفت. کیفیت قطعات تکثیر شده به کمک ژل آگارز 5/1 مورد ارزیابی قرار گرفت. با توجه به طول ژن Cyt b، دو محصول PCR با اندازه‌های 710 و 649 جفت باز (شکل A4 و B4) و برای بخش D-loop، محصول PCR به اندازه 842 جفت باز تکثیر شد (شکل C4).

توالی‌های حاصله از 12 نمونه از لحاظ کیفیت توالی‌یابی برای هر باز توسط نرم افزار Chromas Lite 2.01 و الگوریتم PHRED (Ewing et al., 1998) مورد آنالیز قرار گرفته و کیفیت توالی‌یابی برای هر باز مشخص شد. توالی‌هایی با کیفیت پایین بر اساس انتخاب بلندترین ناحیه با آستانه میانگین خطای کمتر از 1/0 درصد، از دو انتهای قطعات حذف شدند. در نتیجه توالی‌هایی با طول کوتاه‌تر از مقدار مورد انتظار برای نمونه‌ها حاصل شد.

 

 

 

 

 

 

 

 


 

شکل 3- نتایج استخراج DNA بر روی ژل آگارز.

Figure 3- Extracted DNA results on agarose gel.

 

شکل 4- نتایج تکثیر قطعات Cyt b و D-loop روی ژل آگارز. A: قطعه 710 جفت بازی Cyt b، B: قطعه 649 جفت بازی Cyt b، C: قطعه 842 جفت بازی D-loop. چاهک های خالی نشان دهنده کنترل منفی و مارکر استفاده شده 100+ شرکت Thermo می باشد.

Figure 4- Agarose gel electrophoresis of PCR products of Cyt b and D-loop region. A: 710 bp fragments of Cyt b PCR products, B: 649 bp fragments of Cyt b PCR products, C: 842 bp fragments of D-loop PCR products. Empty wells are negative PCR control and 100+ Marker (Thermo, USA) was used in all of gel electrophoresis.


 

 

 

 

به طور میانگین بیش از 97 درصد نوکلئوتیدهای این توالی ها صحت توالی‌یابی معادل 99 درصد را نشان دادند که نشان دهنده ضریب اطمینان بالا در نتیجه گیری بر اساس توالی‌های بدست آمده می‌باشد. از آنجایی که در طراحی پرایمرها نواحی همپوشان بین قطعات در نظر گرفته شده بود، با حذف توالی‌های نامرغوب از دو انتهای توالی‌ها یکپارچگی قطعات مورد مطالعه کاملا حفظ شده بود و در نتیجه یک قطعه به طول 1962 باز از ژنوم میتوکندری بدست آمد که شامل کل ناحیه Cyt b و بخشی از ناحیه کنترل می‌باشد (شکل 5).

 

 

 

شکل 5- توالی های بدست آمده بعد از پردازش دارای همپوشانی کافی به منظور یکپارچه سازی می­باشند.

Figure 5- Sequences obtained after processing, had sufficient overlap in order to assemble the consensus.

 

 

توالی‌های به دست آمده با نرم افزار DIVEIN و Arlequin suite v.3.5  بررسی شدند. این ابزار تعداد جهش مشاهده شده، جایگاه‌های پلی­مورفیک، سطح تنوع نوکلئوتیدی و ساختار و تنوع هاپلوتایپی را در جمعیت مورد مطالعه مشخص می نماید. نتایج نشان داد که در هر سه ناحیه مورد بررسی بین 12 نمونه مورد مطالعه هیچ اختلاف نوکلئوتیدی وجود ندارد و جمعیت در این بخش از ژنوم میتوکندری فاقد تنوع ژنتیکی می‌باشد. این پدیده می­تواند نشان دهنده شدت بالای همخونی و روابط شدید خانوادگی در جمعیت کروکودیل پوزه کوتاه مورد مطالعه باشد و به طور آشکاری نشان‌دهنده وضعیت بحرانی تنوع زیستی این جمعیت است. تنوع زیستی لازمه بقای یک گونه در محیط می باشد و سطح پایین تنوع، کل جمعیت را در خطر انقراض در مواجه با شرایط خاص مانند تغییرات اقلیمی و یا پیدایش بیماری‌های خاص قرار می دهد. تنوع ژنتیکی درون جمعیتی از دو لحاظ برای بقای بلند مدت جمعیت ها مهم می باشد. اول اینکه تنوع فنوتیپی وابسته به ژنوتیپ که تعیین کننده میزان سازگاری یک جمعیت می باشد محصول تنوع ژنتیکی است (Blows & Hoffmann, 2005؛Bell & Collins 2008; Hoffmann & Sgrò, 2011). ثانیا تنوع ژنتیکی خنثی (Neutral Genetic Diversity) در یک جمعیت طبیعی انعکاس دهنده هم‌خونی و رانش ژنتیکی است که هر دو مورد از عوامل مشخص کاهش توان بقای جمعیت­ها هستند (Reed & Frankham 2003; Puurtinen et al., 2004; O’grady et al., 2006). با این وجود نمی‌توان انتظار داشت که تنوع ژنتیکی خنثی مستقیما مترادف با تنوع در هر ویژگی فنوتیپیک خاصی باشد (Reed & Frankham, 2001)، اما بخاطر ارتباط آن با اندازه جمعیت موثر می تواند انعکاس دهنده پتانسیل سازگاری جمعیت ها در بلند مدت باشد (Frankham, 2005; Willi et al., 2006; Lanfear et al., 2014). در نتیجه کاهش تنوع ژنتیکی درون جمعیتی (چه تنوع ژنتیکی انطباقی و در نتیجه فشار انتخاب باشد و چه تنوع ژنتیکی خنثی و نهفته در ژنوم باشد) با افزایش خطر انقراض جمعیت‌های طبیعی مرتبط است.DNA میتوکندریایی چون تحت تاثیر نوترکیبی قرار نمی‌گیرد (زیرا هاپلوئید است)، تنوع خود را تنها مدیون انتخاب طبیعی است و در نتیجه تنوع ژنتیکی آن وابسته به اندازه جمعیت نیست. همچنین انتخاب شدیدی بر ضد جهش‌های مضر ژنوم میتوکندری در لایه زاینده جنس ماده صورت می‌گیرد که منجر به همگن شدن بیشتر سطح تنوع آن می‌گردد (انتخاب پس زمینه ای) (Stewart et al., 2008). با این حال چندین مطالعه در سالهای اخیر DNA میتوکندریایی به عنوان مارکر اندازه جمعیت را زیر سوال برده­اند (Ballard & Whitlock, 2004; Hurst & Jiggins, 2005; Bazin et al., 2006). مهمترین این مطالعات شواهدی مبنی بر اینکه تنوع ژنتیکی در DNA میتوکندریایی با اندازه جمعیت همبستگی ندارد را ارائه می‌دهد (Bazin et al., 2006). این نتایج با مطالعه و مقایسه گروه های تاکسونومیک بزرگ جانوری بدست آمده است. به عنوان مثال، سطح پلی مورفیسم بین مهره داران و بی مهرگان مشابه یکدیگر بوده‌اند، در حالیکه بی مهره گان اندازه جمعیت بسیار بزرگتری در طبیعت دارند.

نتایج بررسی فیلوژنتیکی بر اساس توالی کامل Cyt b کروکودیل پوزه کوتاه ایران با دیگر اعضای خانواده کروکودیلینه (taxid:8493) نمایش داده شده است (شکل 6). این مقایسه با استفاده از توالی‌های ژنوم مرجع هر گونه انجام شده است. به طور کلی 1149 جایگاه و 12 گونه مختلف خانواده کروکودیلینه در این بررسی مورد استفاده قرار گرفته است و نتایج بدست آمده مبتنی بر داده های نوکلئوتیدی و کدون­های پروتئینی به طور مشترک می‌باشد.

درخت فیلوژنتیکی بدست آمده براساس فاصله ژنتیکی و به کمک روش حداکثر درست‌نمایی ترکیبی بدست آمده است (شکل 7). برای محاسبه این مقادیر علاوه بر داده‌های نوکلئوتیدی از داده‌های پروتئینی و نرخ جهش در سه جایگاه کدون پروتئینی استفاده شده است. در نتیجه علاوه بر مدل تکاملی خنثی، مدل تکامل انطباقی نیز در محاسبات مد نظر قرار گرفته است. معیار فاصله ژنتیکی تعداد جایگزینی‌های نوکلئوتید در هر جایگاه بین توالی‌ها می­باشد.

همانطور که در شکل 7 مشاهده می شود، نزدیکترین گونه به گونه گاندو Crocodylus palustris یا ماگر است که مطابق با یافته های مورفولوژیکی می باشد که قبلا ارائه شده است. این نکته بدین مفهوم است که جایگاه رده بندی کروکودیل پوزه کوتاه به عنوان گونه  palustris از دیدگاه ژنتیک مولکولی قابل تایید است. از کلاستر خواهری بین این دو در ادامه با عنوان کلاستر گاندو یاد خواهیم کرد (شکل 6).

پس از ماگر نزدیکترین گونه به گاندو Crocodylus porosus یا کروکودیل آب شور می‌باشد. این در حالی می باشد که این کروکودیل از لحاظ مورفولوژیکی تفاوت های زیادی با گاندو دارد. به عنوان مثال کروکودیل آب شور بزرگترین خزنده دنیا محسوب می شود و طول جنس نر آن به حدود 3/6 بالغ می شود در حالیکه جنس نر گاندو حدود 4 الی 5 متر رشد می‌کند. این کروکودیل عمدتا در سواحل شرقی هند زندگی می‌کند و در نواحی دیگر هند به ندرت یافت می‌شود (Devos, 1984). تشابه بین این دو گونه ریشه در تاریخ تکامل آن دو دارد و اگر به خاستگاه هر دو گونه توجه شود، به کشور هند خواهیم رسید.

در حال حاظر ماگر و کروکودیل آب شور دو گونه از سه گونه کروکودیل می باشد که در هند زندگی می‌کنند. مطالعات متعدد دیگری نیز وجود دارند که خویشاوندی بین این دو گونه را مورد تایید قرار می‌دهند. با توجه به اینکه کروکودیل آب شور و کروکودیل سیامنسیس تشکیل یک کلاستر خواهری را می دهند. این کلاستر در ادامه با عنوان کلاستر پروسوس و کلاستر خواهری گاندو به عنوان یک کلاد گاندو معرفی می شوند (شکل 6).

 


 

 

 

شکل 6- درخت فیلوژنتیک بین گونه های مختلف خانواده کروکودیل و جایگاه گاندو در این خانواده. طول شاخه ها مبین فاصله ژنتیکی و رنگ هر یک نشان دهنده سطح اطمینان محاسبه بر اساس 100 تکرار آزمون بوت استرپ می باشد. خط چین قرمز نشان دهنده کلاستر خواهری گاندو، خط چین آبی نشان دهنده  کلاستر پروسوس و خط چین سبز نشان دهنده کلاد گاندو می باشد.

Figure 6- Phylogenetic tree of different species of crocodylinae family and Gando position in this tree. The length of the branches indicates the genetic distance and the color of each represents the confidence level of the calculation based on 100 repetitions of the Bootstrap test. The red dashed line represents the Gando's sister cluster, the blue dashed line represents the Prosus cluster and the green dashed line indicates the Gando's clad.

 

 

شکل 7- فاصله ژنتیکی بین گاندو و دیگر گونه های خانواده کروکودیلینه که بر اساس روش حداکثر درست نمایی بدست آمده است. مقادیر فاصله ژنتیکی در بخش پایینی جدول و مقادیر خطای استاندارد محاسبه در بخش فوقانی جدول گزارش شده است.

Figure 7- The genetic distance between Gando and other species of the crocodylinae family obtained based on the maximum likelihood method. The genetic distance values in the lower part of the table and the values of the standard error calculation are reported in the upper part of the table.

 

 

وجود رابطه خواهری بین کروکودیل آب شور و کروکودیل سیامنسیس و تشکیل کلاد با ماگر در چند مطالعه دیگر نیز گزارش شده است (Mcaliley et al., 2006; Meganathan et al., 2010). البته در مطالعات دیگری رابطه نزدیک بین کروکودیل آب شور و ماگر و تشکیل کلاد با کروکودیل سیامنسیس گزارش شده است (Densmore & Owen, 1989; Gatesy & Amato, 2008). مطالعاتی دیگری نیز با بررسی 12 ژن کد کننده پروتئین در میتوکندری وجود رابطه نزدیک بین کروکودیل سیامیس و ماگر و تشکیل کلاد با کروکودیل آب شور نشان داده شده است (Man et al., 2011). همین نتیجه در یک بررسی دیگر نیز تایید شده است (Oaks, 2011). در مطالعه ای مشابه قرار گیری ماگر، سیامیس و کروکودیل آب شور در یک کلاد بر اساس مطالعه کل DNA میتوکندری مورد تایید قرار گرفته است (Meredith et al., 2011).

باید توجه داشت که فاصله ژنتیکی بین گاندو و دیگر گونه‌ها، مانند کروکودیل گینه نو کمتر از فاصله آنها با کلاستر پروسوس می باشد، اما روند تکاملی بین آنها متفاوت از یکدیگر می باشد و در نتیجه گاندو با پروسوس تشکیل کلاد داده است و کروکودیل گینه نو در این کلاد قرار نگرفته است. با بررسی بانک ژنی، 16 توالی برای Cyt b در گونه Crocodylus palustris شناسایی شد. مقایسه این توالی‌ها با توالی Cyt b گاندو که در این مطالعه بدست آمد نشان دهنده وجود 12 جایگاه پلی مورفیک در بین 1156 جایگاه مورد مطالعه بود. 3 مورد از این جایگاه‌ها بصورت سینگلتون با دو واریانت می باشند و 9 جایگاه دیگر از نوع جایگاه‌های حاوی اطلاعات پارسیمونی هستند که در تعیین جایگاه هر نمونه در ساختار فیلوژنتیکی موثر می باشند (جدول 3 ).

 

 

جدول 3- جایگاه های پلی مورفیک در توالی های Cyt b گونه C. palustris (جایگاه های پارسیمونیک بصورت پر رنگ نمایش داده شده است).

Table 3- Polymorphic positions in the Cyt b sequence of C. palustris (parsimony positions were shown in bold).

جایگاه/نمونه

139

270

289

306

475

511

570

603

712

750

876

1127

JF315283

G

A

C

A

G

G

C

C

G

G

T

T

HQ595025

.

.

.

.

.

A

.

.

.

.

.

.

JF315254

.

.

.

.

.

A

.

.

.

.

.

.

GU144286

.

.

.

.

.

A

.

.

.

.

.

.

FJ173286

A

G

.

G

A

.

.

.

.

A

C

C

HM488007

.

.

A

G

A

.

.

G

A

A

C

C

FJ173285

.

.

A

G

A

.

.

G

A

A

C

C

FJ173284

.

.

A

G

A

.

.

G

A

A

C

C

FJ173283

.

.

A

G

A

.

.

G

A

A

C

C

FJ173282

.

.

A

G

A

.

.

G

A

A

C

C

FJ173281

.

.

A

G

A

.

.

G

A

A

C

C

FJ173280

.

.

A

G

A

.

.

G

A

A

C

C

FJ173279

.

.

A

G

A

.

.

G

A

A

C

C

FJ173278

.

.

A

G

A

.

.

G

A

A

C

C

FJ173277

.

.

A

G

A

.

.

G

A

A

C

C

FJ173276

.

.

A

G

A

.

.

G

A

A

C

C

Gando  

.

.

.

.

.

.

T

.

.

.

.

.

 

 

تنوع هاپلوتایپی بدست آمده برابر 574/0 بود، که نشان دهنده میزان پراکنش بالایی در ساختار ژنتیکی کروکودیل ماگر می باشد. با آنالیز مجدد فیلوژنتیکی بر اساس مجموعه توالی های Cyt b مربوط به مطالعات مختلف در کنار نتایج این مطالعه به درخت فیلوژنتیکی جدیدی منجر گردید که نشان‌دهنده اختلاف درون گونه‌ای برای کروکودیل ماگر می‌باشد (شکل 8). حتی با استفاده از چندین نمونه مختلف بازهم تمامی اعضای گونه ماگر در یک کلاستر قرار می گیرند و با فاصله زیادی تشکیل یک کلاد با دو کروکودیل اب شور و پروسوس را می دهند که نشان دهنده صحت نتیجه گیری اولیه در مورد رابطه بین کلاستر پروسوس و گاندو می باشد. این مسئله بدین مفهوم است که تنوع درون گونه ای ماگر منجر به تغییر جایگاه آن در ساختار فیلوژنتیکی خانواده کروکودلینه نمی شود.

نتایج این مطالعه نشان داد که مطابق با یافته های مورفولوژیکی، نزدیکترین گونه از لحاظ ژنتیکی به گاندو، گونه Crocodylus palustris است. با این حال، با توجه به منطقه جغرافیایی نمونه‌های مورد بررسی، فاصله ژنتیکی بین این دو وجود دارد. نکته حائز اهمیت در این پژوهش، عدم تنوع ژنتیکی در میان گاندوهای سه منطقه مورد بررسی در شهرستان چابهار است. فقدان تنوع ژنتیکی در داخل یک گونه می تواند توانایی آن گونه در سازگاری با تغییرات شرایط محیطی را تحت تاثیر قرار داده و در نتیجه خطر انقراض آنها را افزایش بدهد. از آنجاییکه هدف نهایی برنامه‌های حفاظت از گونه‌ها، بقای بلند مدت آنهاست، حفاظت از تنوع ژنتیکی درون جمعیتی باید از جایگاه ویژه ای برخوردار باشد. لذا جلوگیری از هم خونی بیشتر در میان گاندوها با استفاده از برنامه‌های تلاقی مناسب پیشنهاد می‌گردد.

 

تشکر و قدردانی

این مقاله مستخرج شده از پروژه پژوهشی به شماره 100807 مصوب معاونت پژوهش و فناوری دانشگاه فردوسی مشهد می باشد. لذا، بدین وسیله از آن معاونت محترم تقدیر و تشکر می گردد. تصویب این پروژه از محل حمایت‌های مالی اداره کل حفاظت محیط زیست سیستان و بلوچستان می‌باشد که مراتب تشکر و قدردانی را از آن اداره کل محترم داریم. پروژه در آزمایشگاه دانشکده کشاورزی انجام گرفته و بدین ترتیب از مسئولان مربوطه کمال تشکر را داریم.

 

 

 

 

 

 

شکل 8- جایگاه کروکودیل­های گونه ماگر در قیاس با کلاستر پروسوس. تنوع درون گونه­ای ماگر باعث تغییر موقعیت آن نمی­شود.

Figure 8- Position of Mugger species of crocodylus in comparison to Crocodylus porosus cluster. Inter-species variation of the Mugger does not change its position.


منابع

 

Ballard JWO, Whitlock MC (2004). The incomplete natural history of mitochondria. Molecular Ecology 13: 729-744.

Bazin E, Glemin S, Galtier N (2006). Population size does not influence mitochondrial genetic diversity in animals Science. 31: 570-572.

Bell G, Collins S (2008). Adaptation, extinction and global change. Evolutionary Applications 1: 3-16.

Bellagamba F, Moretti VM, Comincini S, Valfare F (2001). Identification of species in animal feedstuffs by polymerase chain reaction restriction fragment length polymorphism analysis of mitochondrial DNA. Journal of Agricultural and Food Chemistry 49: 75-81.

Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ (2013). GenBank. Nucleic Acids Research 41: 36-42.

Blows MW, Hoffmann AA (2005). A reassessment of genetic limits to evolutionary change. Ecology 86: 1371-1384.

Brochu C (2000). Phylogenetic relationships and divergence timing of Crocodylus based on morphology and the fossil record. Copeia 3: 657–673.

Brochu C, Densmore L (2001). Crocodile phylogenetics: a summary of current progress. In: Grigg, G., Seebacher, F., Franklin, C. (Eds.), Crocodilian Biology and Evolution. Surrey Beatty and Sons, Chipping Norton, NSW, pp. 3–8.

Brochu C, Njau J, Blumenschine R, Densmore L (2010). A new horned crocodile from the Plio-Pleistocene hominid sites at Olduvai Gorge, Tanzania. PloS ONE 5: 1-13.

Bruford MW, Bradley DG, Luikart G (2003). DNA markers reveal the complexity of livestock domestication. Nature Reviews Genetics 4: 900-910.

Da Silva A, Lenin J (2010). Mugger crocodile Crocodylus palustris. Crocodile: Status Survey and Conservation Action Plan (Third ed.), Eds. Manolis SC, Stevenson 94-98.

de Vos A. 1984. Crocodile conservation in India. Biological conservation 29: 183-189.

Deng W, Maust BS, Nickle DC, Learn GH, Liu Y (2010). DIVEIN: a web server to analyze phylogenies, sequence divergence, diversity, and informative sites. Biotechniques 48: 405-417.

Densmore LD, Owen RD (1989). Molecular systematics of the order Crocodilia. American Zoologist 29: 831-841.

Dolati V, Evrigh N, Nikbin S, Behmaram R (2016). Phylogenetic analyzing of Khalkhali goats using cytochrome b gene. Journal of Agriculture Biotechnology 9: 75-86 (In Farsi).

Ewing B, Green P (1998). Base-Calling of Automated Sequencer Traces Using Phred. II. Error Probabilities. Genome Research 8: 186-194.

Excoffier L, Lischer HE (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular ecology resources 10: 564-567.

Farris JS (1970). Methods for computing Wagner trees. Systematic Biology. 19: 83-92.

Farris JS (1989). The retention index and the rescaled consistency index. Cladistics 5: 417-419.

Felsenstein J (1981). Evolutionary trees from DNA sequences: a maximum likelihood approach. Journal of molecular evolution 17: 368-376.

Felsenstein J (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution 783-791.

Feng G, Wu X, Yan P, Li X (2010). Two complete mitochondrial genomes of Crocodylus and implications for crocodilians phylogeny. Amphibia – Reptilia 31: 299–309.

Fitch WM (1971). Toward defining the course of evolution: minimum change for a specific tree topology. Systematic Biology 20: 406-416.

Gatesy J, Amato G (2008). The rapid accumulation of consistent molecular support for intergeneric crocodylian relationships. Molecular Phylogenetics and Evolution 48: 1232– 1237.

Gatesy J, Baker R, Hayashi C (2004). Inconsistencies in arguments for the supertree approach: supermatrices versus supertrees of Crocodylia. systematic biology journal. 53: 342–355.

Hoffmann AA, Sgrò CM (2011). Climate change and evolutionary adaptation. Nature 470: 479-485.

Hurst GD, Jiggins FM (2005). Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts. Proceedings of the Royal Society of London B: Biological Sciences 272: 1525-1534.

Janke A, Gullberg A, Hughes S, Aggarwal R, Arnason U (2005). Mitogenomic analyses place the gharial (Gavialis gangeticus) on the crocodile tree and provide pre-K/T divergence times for most crocodilians. Journal of molecular evolution 61: 620–626.

Javanmard A, Mohammadabadi MR, Zarrigabayi GE, Gharahedaghi AA, Nassiry MR, Javadmansh A, Asadzadeh N (2008). Polymorphism within the intron region of the bovine leptin gene in Iranian Sarabi cattle (Iranian Bos taurus). Russian Journal of Genetics 44: 495-497.

Lalitha S (2000). Primer premier 5. Biotech Software & Internet Report: The Computer Software Journal for Science 1: 270-272.

Lanfear R, Kokko H, Eyre-Walker A (2014). Population size and the rate of evolution. Trends in Ecology & Evolution 9: 33-44.

Larkin MA, Blackshields G, Brown N, Chenna R, McGettigan PA. (2007). Clustal W and Clustal X version 2.0. Bioinformatics 23: 2947-2948.

Lartillot N, Rodrigue N, Stubbs D, Richer J (2013). PhyloBayes MPI. Phylogenetic reconstruction with infinite mixtures of profiles in a parallel environment. Systematic Biology 62: 611-616.

Man Z, Yishu W, Peng Y, Xiaobing W (2011). Crocodilian phylogeny inferred from twelve mitochondrial protein-coding genes, with new complete mitochondrial genomic sequences for Crocodylus acutus and Crocodylus novaeguineae. Molecular Phylogenetics and Evolution 60: 62-67.

McAliley LR, Willis RE, Ray DA, White PS, Brochu CA (2006). Are crocodiles really monophyletic? Evidence for subdivisions from sequence and morphological data. Molecular Phylogenetics and Evolution 39: 16-32.

McCarthy C (2005). Chromas, pp. Technelysium Pty Ltd.

Meganathan P, Dubey B, Batzer M, Ray D, Haque I (2010). Molecular phylogenetic analyses of genus Crocodylus (Eusuchia, Crocodylia, Crocodylidae) and the taxonomic position of Crocodylus porosus. Molecular Phylogenetics and Evolution 57: 393–402.

Meredith RW, Hekkala ER, Amato G, Gatesy J (2011). A phylogenetic hypothesis for Crocodylus (Crocodylia) based on mitochondrial DNA: evidence for a trans-Atlantic voyage from Africa to the New World. Molecular Phylogenetics and Evolution 60: 183-191.

Mohammadi A, Nassiry MR, Mosafer J, Mohemmadabadi MR, Sulimova GE (2009). Distribution of BoLA-DRB3 allelic frequencies and identification of a new allele in the Iranian cattle breed Sistani (Bos indicus). Russian Journal of Genetics 45: 198-202.

Mousavizadeh A, Mohammad Abadi MR, Torabi A, Nassiry MR, Ghiasi H, Esmailizadeh AK (2009). Genetic polymorphism at the growth hormone locus in Iranian Talli goats by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Iranian Journal of Biotechnology 7: 51-53.

Naderi S, Rezaei HR, Taberlet P, Zundel S, Rafat SA, Naghash HR, El-Barody MAA, Ertugrul O, Pompanon F (2007) Large-Scale Mitochondrial DNA Analysis of the Domestic Goat Reveals Six Haplogroups with High Diversity. PLoS ONE  2(10): 1-12.

Nassiri MR, Parizadeh SA, Mahdavi M, Ariannejad H (2010). Genetic and phylogenetic analysis of D-Loop in Persian Leopard. Journal of Agriculture Biotechnology 3: 82-95 (In Farsi).

Nassiry M, Tahmorespour M, Javadmanesh A, Soltani M, Foroutani Far S (2006). Calpastatin polymorphism and its association with daily gain in Kurdi sheep. Iranian Journal of Biotechnology 4: 188-192.

O’Grady JJ, Brook BW, Reed DH, Ballou JD, Tonkyn DW (2006). Realistic levels of inbreeding depression strongly affect extinction risk in wild populations. Biological Conservation 133: 42-51.

Oaks JR (2011). A time‐calibrated species tree of crocodylia reveals a recent radiation of the true crocodiles. Evolution 65: 3285-3297.

Pruitt KD, Tatusova T, Maglott DR (2007). NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic acids research 35: 61-65.

Puurtinen M, Knott KE, Suonpää S, Ooik Tv, Kaitala V (2004). Genetic variability and drift load in populations of an aquatic snail. Evolution 58: 749-756.

Reed DH, Frankham R (2003). Correlation between fitness and genetic diversity. Conservation Biology 17: 230-237.

Roos J, Aggarwal R, Janke A (2007). Extended mitogenomic phylogenetic analyses yield new insight into crocodilian evolution and their survival of the Cretaceous– Tertiary boundary. Molecular Phylogenetics and Evolution 45: 663–673.

Saitou N, Nei M (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Biology and Evolution 4: 406-425.

Sayers EW, Barrett T, Benson DA, Bolton E, Bryant SH (2011). Database resources of the national center for biotechnology information. Nucleic Acids Research 39: 38-51.

Stamatakis A (2006). RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688-2690.

Stamatakis A, Hoover P, Rougemont J (2008). A rapid bootstrap algorithm for the RAxML web servers. Systematic Biology 57: 758-771.

Stewart JB, Freyer C, Elson JL, Larsson N (2008). Purifying selection of mtDNA and its implications for understanding evolution and mitochondrial disease. Nature Reviews Genetics 9: 657-662.

Swindell S, Plasterer T (1997). SEQMAN, in Sequence Data Analysis Guidebook, edited by S. Swindell. Springer New York. pp. 75-89

Swofford DL (2003). PAUP. Phylogenetic analysis using parsimony (and other methods). Version 4.

Whitaker R, Andrews H (2003). Crocodile conservation, Western Asia Region: an update. Journal of the Bombay Natural History Society 100: 432-445.

Willi Y, Van Buskirk J, Hoffmann AA (2006). Limits to the adaptive potential of small populations. Annual Review of Ecology, Evolution and Systematics 37: 433-458.

Willis R, McAliley L, Neeley E, Densmore L (2007). Evidence for placing the false gharial (Tomistoma schlegelii) into the family Gavialidae: inferences from nuclear gene sequences. Molecular Phylogenetics and Evolution 43: 787–794.

Yang Z, Rannala B (1997). Bayesian phylogenetic inference using DNA sequences: A Markov chain Monte Carlo method. Molecular Biology and Evolution 14: 717-724.

Zamani P, Akhondi M, Mohamadabadi M, Saki AA, Ershadi A, Banabazi MH, Abdolmohammadi AR (2013). Genetic variation of Mehraban sheep using two inter simple sequence repeat (ISSR) markers. African Journal of Biotechnology 10: 1812-1817.

 

 

 

 


 

Genetic Diversity of Persian Crocodile Crocodylus Palustris Using Sequencing of D-Loop and Cyt b Regions of Mitochondria

 

Afsharian A.1, Nassiri M.R.*2, Ebrahimi E.3, Javadmanesh A.4

 

1PhD Student, Animal Science Department, Faculty of Agriculture, Ferdowsi University of Mashhad, Mashhad, Iran.

2Professor, Animal Science Department, Faculty of Agriculture, Ferdowsi University of Mashhad. Mashhad, Iran.

3 Associate Professor, School of Molecular and Biomedical Science, The University of Adelaide, Australia.

 4 Assistant Professor, Animal Science Department, Faculty of Agriculture, Ferdowsi University of Mashhad, Mashhad, Iran.

 

Abstract

Identification of the genetic features of Persian mugger crocodile, the Gando, is crucial in management planning for conservation and preservation of this endangered species. Therefore, this study is aimed to investigate the genetic characteristics of Gando in Iran, based on population studies of the conserved and variable regions of the mitochondrial DNA. For this purpose, 12 samples of skin tissue and blood were collected from three regions of Chabahar city, south east of Iran. After DNA extraction, conserved and variable regions of the mitochondrial genome Cyt b and D-loop were amplified using PCR. The PCR products were sequenced and results were analyzed using bioinformatics software based on nucleotide data as well as protein codons. Results revealed that all samples had similar nucleotide sequences and no genetic variation was found among the samples. The closest species to Gando were Crocodylus palustris and then C. porosus. Comparison of sequences with other Crocodile species revealed the presence of 12 polymorphic loci among 1156 loci studied and also 5 haplogroups were introduced. As we didn’t find any nucleotide diversity among studied samples, therefore, to increase the genetic diversity of the Gando, it is essential to include specific conservational plans for this endangered species.

Keywords: Gando Crocodile, Genetic Diversity, Phylogenetic Tree, Mitochondria Genome.



*  نویسنده مسئول: احمد افشاریان                               تلفن: 09153114119 Email:  nassiryr@um.ac.ir                                      

[1] Neighbor Joining

[2] Maximum Parsimony

* Corresponding author: Nassiri M.R.           Tel: +98 51 38795618                              Email: nassiryr@um.ac.ir

 
Ballard JWO, Whitlock MC (2004). The incomplete natural history of mitochondria. Molecular Ecology 13: 729-744.
Bazin E, Glemin S, Galtier N (2006). Population size does not influence mitochondrial genetic diversity in animals Science. 31: 570-572.
Bell G, Collins S (2008). Adaptation, extinction and global change. Evolutionary Applications 1: 3-16.
Bellagamba F, Moretti VM, Comincini S, Valfare F (2001). Identification of species in animal feedstuffs by polymerase chain reaction restriction fragment length polymorphism analysis of mitochondrial DNA. Journal of Agricultural and Food Chemistry 49: 75-81.
Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ (2013). GenBank. Nucleic Acids Research 41: 36-42.
Blows MW, Hoffmann AA (2005). A reassessment of genetic limits to evolutionary change. Ecology 86: 1371-1384.
Brochu C (2000). Phylogenetic relationships and divergence timing of Crocodylus based on morphology and the fossil record. Copeia 3: 657–673.
Brochu C, Densmore L (2001). Crocodile phylogenetics: a summary of current progress. In: Grigg, G., Seebacher, F., Franklin, C. (Eds.), Crocodilian Biology and Evolution. Surrey Beatty and Sons, Chipping Norton, NSW, pp. 3–8.
Brochu C, Njau J, Blumenschine R, Densmore L (2010). A new horned crocodile from the Plio-Pleistocene hominid sites at Olduvai Gorge, Tanzania. PloS ONE 5: 1-13.
Bruford MW, Bradley DG, Luikart G (2003). DNA markers reveal the complexity of livestock domestication. Nature Reviews Genetics 4: 900-910.
Da Silva A, Lenin J (2010). Mugger crocodile Crocodylus palustris. Crocodile: Status Survey and Conservation Action Plan (Third ed.), Eds. Manolis SC, Stevenson 94-98.
de Vos A. 1984. Crocodile conservation in India. Biological conservation 29: 183-189.
Deng W, Maust BS, Nickle DC, Learn GH, Liu Y (2010). DIVEIN: a web server to analyze phylogenies, sequence divergence, diversity, and informative sites. Biotechniques 48: 405-417.
Densmore LD, Owen RD (1989). Molecular systematics of the order Crocodilia. American Zoologist 29: 831-841.
Dolati V, Evrigh N, Nikbin S, Behmaram R (2016). Phylogenetic analyzing of Khalkhali goats using cytochrome b gene. Journal of Agriculture Biotechnology 9: 75-86 (In Farsi).
Ewing B, Green P (1998). Base-Calling of Automated Sequencer Traces Using Phred. II. Error Probabilities. Genome Research 8: 186-194.
Excoffier L, Lischer HE (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular ecology resources 10: 564-567.
Farris JS (1970). Methods for computing Wagner trees. Systematic Biology. 19: 83-92.
Farris JS (1989). The retention index and the rescaled consistency index. Cladistics 5: 417-419.
Felsenstein J (1981). Evolutionary trees from DNA sequences: a maximum likelihood approach. Journal of molecular evolution 17: 368-376.
Felsenstein J (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution 783-791.
Feng G, Wu X, Yan P, Li X (2010). Two complete mitochondrial genomes of Crocodylus and implications for crocodilians phylogeny. Amphibia – Reptilia 31: 299–309.
Fitch WM (1971). Toward defining the course of evolution: minimum change for a specific tree topology. Systematic Biology 20: 406-416.
Gatesy J, Amato G (2008). The rapid accumulation of consistent molecular support for intergeneric crocodylian relationships. Molecular Phylogenetics and Evolution 48: 1232– 1237.
Gatesy J, Baker R, Hayashi C (2004). Inconsistencies in arguments for the supertree approach: supermatrices versus supertrees of Crocodylia. systematic biology journal. 53: 342–355.
Hoffmann AA, Sgrò CM (2011). Climate change and evolutionary adaptation. Nature 470: 479-485.
Hurst GD, Jiggins FM (2005). Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts. Proceedings of the Royal Society of London B: Biological Sciences 272: 1525-1534.
Janke A, Gullberg A, Hughes S, Aggarwal R, Arnason U (2005). Mitogenomic analyses place the gharial (Gavialis gangeticus) on the crocodile tree and provide pre-K/T divergence times for most crocodilians. Journal of molecular evolution 61: 620–626.
Javanmard A, Mohammadabadi MR, Zarrigabayi GE, Gharahedaghi AA, Nassiry MR, Javadmansh A, Asadzadeh N (2008). Polymorphism within the intron region of the bovine leptin gene in Iranian Sarabi cattle (Iranian Bos taurus). Russian Journal of Genetics 44: 495-497.
Lalitha S (2000). Primer premier 5. Biotech Software & Internet Report: The Computer Software Journal for Science 1: 270-272.
Lanfear R, Kokko H, Eyre-Walker A (2014). Population size and the rate of evolution. Trends in Ecology & Evolution 9: 33-44.
Larkin MA, Blackshields G, Brown N, Chenna R, McGettigan PA. (2007). Clustal W and Clustal X version 2.0. Bioinformatics 23: 2947-2948.
Lartillot N, Rodrigue N, Stubbs D, Richer J (2013). PhyloBayes MPI. Phylogenetic reconstruction with infinite mixtures of profiles in a parallel environment. Systematic Biology 62: 611-616.
Man Z, Yishu W, Peng Y, Xiaobing W (2011). Crocodilian phylogeny inferred from twelve mitochondrial protein-coding genes, with new complete mitochondrial genomic sequences for Crocodylus acutus and Crocodylus novaeguineae. Molecular Phylogenetics and Evolution 60: 62-67.
McAliley LR, Willis RE, Ray DA, White PS, Brochu CA (2006). Are crocodiles really monophyletic? Evidence for subdivisions from sequence and morphological data. Molecular Phylogenetics and Evolution 39: 16-32.
McCarthy C (2005). Chromas, pp. Technelysium Pty Ltd.
Meganathan P, Dubey B, Batzer M, Ray D, Haque I (2010). Molecular phylogenetic analyses of genus Crocodylus (Eusuchia, Crocodylia, Crocodylidae) and the taxonomic position of Crocodylus porosus. Molecular Phylogenetics and Evolution 57: 393–402.
Meredith RW, Hekkala ER, Amato G, Gatesy J (2011). A phylogenetic hypothesis for Crocodylus (Crocodylia) based on mitochondrial DNA: evidence for a trans-Atlantic voyage from Africa to the New World. Molecular Phylogenetics and Evolution 60: 183-191.
Mohammadi A, Nassiry MR, Mosafer J, Mohemmadabadi MR, Sulimova GE (2009). Distribution of BoLA-DRB3 allelic frequencies and identification of a new allele in the Iranian cattle breed Sistani (Bos indicus). Russian Journal of Genetics 45: 198-202.
Mousavizadeh A, Mohammad Abadi MR, Torabi A, Nassiry MR, Ghiasi H, Esmailizadeh AK (2009). Genetic polymorphism at the growth hormone locus in Iranian Talli goats by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Iranian Journal of Biotechnology 7: 51-53.
Naderi S, Rezaei HR, Taberlet P, Zundel S, Rafat SA, Naghash HR, El-Barody MAA, Ertugrul O, Pompanon F (2007) Large-Scale Mitochondrial DNA Analysis of the Domestic Goat Reveals Six Haplogroups with High Diversity. PLoS ONE  2(10): 1-12.
Nassiri MR, Parizadeh SA, Mahdavi M, Ariannejad H (2010). Genetic and phylogenetic analysis of D-Loop in Persian Leopard. Journal of Agriculture Biotechnology 3: 82-95 (In Farsi).
Nassiry M, Tahmorespour M, Javadmanesh A, Soltani M, Foroutani Far S (2006). Calpastatin polymorphism and its association with daily gain in Kurdi sheep. Iranian Journal of Biotechnology 4: 188-192.
O’Grady JJ, Brook BW, Reed DH, Ballou JD, Tonkyn DW (2006). Realistic levels of inbreeding depression strongly affect extinction risk in wild populations. Biological Conservation 133: 42-51.
Oaks JR (2011). A time‐calibrated species tree of crocodylia reveals a recent radiation of the true crocodiles. Evolution 65: 3285-3297.
Pruitt KD, Tatusova T, Maglott DR (2007). NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic acids research 35: 61-65.
Puurtinen M, Knott KE, Suonpää S, Ooik Tv, Kaitala V (2004). Genetic variability and drift load in populations of an aquatic snail. Evolution 58: 749-756.
Reed DH, Frankham R (2003). Correlation between fitness and genetic diversity. Conservation Biology 17: 230-237.
Roos J, Aggarwal R, Janke A (2007). Extended mitogenomic phylogenetic analyses yield new insight into crocodilian evolution and their survival of the Cretaceous– Tertiary boundary. Molecular Phylogenetics and Evolution 45: 663–673.
Saitou N, Nei M (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Biology and Evolution 4: 406-425.
Sayers EW, Barrett T, Benson DA, Bolton E, Bryant SH (2011). Database resources of the national center for biotechnology information. Nucleic Acids Research 39: 38-51.
Stamatakis A (2006). RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688-2690.
Stamatakis A, Hoover P, Rougemont J (2008). A rapid bootstrap algorithm for the RAxML web servers. Systematic Biology 57: 758-771.
Stewart JB, Freyer C, Elson JL, Larsson N (2008). Purifying selection of mtDNA and its implications for understanding evolution and mitochondrial disease. Nature Reviews Genetics 9: 657-662.
Swindell S, Plasterer T (1997). SEQMAN, in Sequence Data Analysis Guidebook, edited by S. Swindell. Springer New York. pp. 75-89
Swofford DL (2003). PAUP. Phylogenetic analysis using parsimony (and other methods). Version 4.
Whitaker R, Andrews H (2003). Crocodile conservation, Western Asia Region: an update. Journal of the Bombay Natural History Society 100: 432-445.
Willi Y, Van Buskirk J, Hoffmann AA (2006). Limits to the adaptive potential of small populations. Annual Review of Ecology, Evolution and Systematics 37: 433-458.
Willis R, McAliley L, Neeley E, Densmore L (2007). Evidence for placing the false gharial (Tomistoma schlegelii) into the family Gavialidae: inferences from nuclear gene sequences. Molecular Phylogenetics and Evolution 43: 787–794.
Yang Z, Rannala B (1997). Bayesian phylogenetic inference using DNA sequences: A Markov chain Monte Carlo method. Molecular Biology and Evolution 14: 717-724.
Zamani P, Akhondi M, Mohamadabadi M, Saki AA, Ershadi A, Banabazi MH, Abdolmohammadi AR (2013). Genetic variation of Mehraban sheep using two inter simple sequence repeat (ISSR) markers. African Journal of Biotechnology 10: 1812-1817.