Document Type : Original Article
Authors
1 Ph. D. Graduated, Department of Fisheries, Faculty of Natural Resources, University College of Agriculture & Natural Resources, University of Tehran, Karaj, Iran
2 Professor, Department of Fisheries, Faculty of Natural Resources, University College of Agriculture & Natural Resources, University of Tehran, Karaj, Iran
3 Associate Professor, Department of Fisheries, Faculty of Natural Resources, University College of Agriculture & Natural Resources, University of Tehran, Karaj, Iran
Abstract
Keywords
Main Subjects
مقدمه
بیشتر مطالعات جهانی، تنوع زیستی را در سطح گونهها (تنوع گونهای) بررسی میکنند و مطالعات کمی تنوع ژنها در موجودات، یعنی تنوع ژنتیکی و تنوع اکوسیستمها (تنوع اکولوژیکی) را بررسی میکنند Ellegren & Galtier, 2016)). درواقع، هزینه نمونهبرداری و تعیین ژنوتیپ تعداد کافی از افراد در گونهها، درک ما را از عوامل تعیینکنندة تنوع ژنتیکی درونگونهای، بهویژه در مقیاس بزرگ، محدود کرده است ( Manel et al., 2020). تنوع زیستی را میتوان بهعنوان تعداد، تنوع و تغییرپذیری موجودات زنده در مقیاس زمانی و مکانی تعریف کرد که شامل آرایههای بالاتر از گونه (جنسها، خانوادهها و غیره) تا گونهها، جمعیتها، زیرجمعیتها، افراد و ژنها میشود ( Bickham et al., 2000). برهمخوردن تعادل ژنتیکی در هریک از این سطوح، ارتباط مستقیمی با کاهش تنوع و بهتبع افزایش آسیبپذیری در برابر استرسهای محیطی و انقراض گونهها دارد. هم انقراض و هم فرگشت گونههای جدید، بخشی جداییناپذیر از تاریخ حیات روی زمین بودهاند؛ اما امروزه بهدلیل تخریب یا تغییر سریع زیستگاههای طبیعی، افزایش فعالیتهای کشاورزی و آلودگیهای شیمیایی، میزان انقراض کنونی 10 تا 100 برابر پیشینه تاریخی تخمین زده میشود ( Pimm et al., 1995).
تنوع ژنتیکی، بیانکنندة تفاوت در تعداد و نوع آللهای موجود در لوکوسهای کروموزومی است و طی چند سال اخیر علاقهمندان زیادی را بهدلیل نقش کلیدی در ساختار ذخایر طبیعی به خود جلب کرده است ( Shabani et al., 2012). درواقع تنوع ژنتیکی یک گونه، ظرفیت سازگاری و پتانسیل تکاملی آن را تعیین میکند. جمعیتهای کوچک متعلق به گونههایی با توزیع محدود، اغلب دارای تنوع ژنتیکی کم در بین جمعیتها هستند؛ اما تمایز ژنتیکی بالایی در میان جمعیتها بهدلیل رانش ژنتیکی و محدودیت جریان ژنی دارند ( Zhai et al., 2019). امروزه تنوع ژنتیکی در بسیاری از برنامههای حفاظتی درخور توجه قرار گرفته است؛ زیرا تنوع ژنتیکی کم با بالارفتن میزان همخونی، تجمع جهشهای مضر و کاهش پتانسیل سازگاری همراه است ( McCusker & Bentzen, 2010)؛ بنابراین، تنوع یا گوناگونی ژنتیکی یک ویژگی اساسی هر جمعیت برای شایستگی بقای افراد آن جمعیت و همچنین بقای کل جمعیت است که امکان سازگاری با شرایط متغیر محیطی و استرس را فراهم میکند. ازاینرو، تضعیف تنوع ژنتیکی یک گونه، قابلیت سازگاری آن را کاهش و خطر انقراض آن را افزایش میدهد ( Mukhopadhyay & Bhattacharjee, 2014).
طی دهههای گذشته، ژنتیک حفاظت از یک حوزه عمدتاً مبتنی بر نظریات بیولوژی جمعیت به یک حوزه تجربی کامل ارتقا یافته است و پیشرفتهای تکنولوژیکی در ژنتیک مولکولی منجر به استفاده گسترده از نشانگرهای مولکولی در زیستشناسی حفاظت شده است ( Ouborget al., 2010). ریزماهوارهها بهعنوان رایجترین و چندکارهترین نوع نشانگر برای کاربردهای زیستمحیطی ظهور کردهاند (Selkoe & Toonen, 2006). لوکوسهای ریزماهواره نشانگرهای ژنتیکی استانداردی برای تجزیه و تحلیل ژنتیکی جمعیت هستند (Coates et al., 2009). ریزماهوارهها که با نامهای توالیهای کوتاه تکراری (STR)، توالیهای ساده تکراری (SSRs) یا پلیمورفیسمهای طولی با توالی ساده (SSLP) شناخته میشوند، موتیفهای تکرارشوندة پشت سر هم از 1 تا 6 جفتباز هستند که در تمام ژنومهای پروکاریوتی و یوکاریوتی که تا به امروز تجزیه و تحلیل شدهاند، به وفور حضور دارند. با توجه به تعداد نوکلئوتیدها در واحد تکرار، ریزماهوارهها بهعنوان مونو، دی، تری، تترا، پنتا یا هگزانوکلئوتید طبقهبندی میشوند. بیشتر ریزماهوارههای ژنومی ریزماهوارههای هستهای هستند؛ با این حال، ریزماهوارهها در میتوکندریها و کلروپلاستها نیز توزیع شدهاند (Kalia et al., 2011).
اکوسیستمهای آب شیرین ازجمله در معرض خطرترین اکوسیستمهای جهان هستند و اگرچه تنها 8/0 درصد از سطح زمین را تشکیل میدهند، میزبان 6 درصد از تنوع گونهها هستند (Dudgeon et al., 2006). تهدیدهای کلیدی در اکوسیستمهای آب شیرین شامل گونههای مهاجم، تخریب و تکهتکهشدن زیستگاه، برداشت بیرویه آب، آلودگیهای صنعتی، کشاورزی، نوری، صوتی و فعالیتهای شیلاتی است (Faulks et al., 2017). در اکوسیستمهای آب شیرین، تنوع زیستی ماهیها بهعنوان شاخص خوبی برای سلامت اکوسیستم در نظر گرفته میشود (Zhou et al., 2024). ماهیان آب شیرین 40 درصد از تنوع کل ماهیهای جهان را تشکیل میدهند و یکی از گروههای جانوری در معرض خطر هستند (Costa et al., 2021). در دهههای اخیر، جمعیت ماهیان آب شیرین در نتیجة تهدیدات مختلف، ازجمله نابودی و تخریب زیستگاه، معرفی گونههای مهاجم، صید بیرویه، آلودگی و تغییرات اقلیمی، دچار افت شدید و کاهش دامنه شده است (Dudgeon, 2010). ازنظر تنوع و بومزایی، کشور ایران دارای غنای بالایی از ماهیان آب شیرین است (Sayyadzadeh & Esmaeili, 2024)؛ بهطوریکه ماهیان آبهای شیرین ایران حتی بدون در نظر گرفتن ماهیان آبهای لب شور دریای خزر بسیار متنوع و جالب توجه هستند (Bagherian nejad, 2014).
طبق آخرین مطالعات، در مجموع 300 گونه ماهی متعلق به 38 خانواده و 110 جنس در آبهای داخلی ایران زیست میکنند (Sayyadzadeh & Esmaeili, 2024). جنس ماهی سفید رودخانهای Squalius یکی از جنسهای مهم در خانواده سفیدماهیان (Leuciscidae) به شمار میرود. پنج گونه از این جنس در آبهای داخلی ایران شناسایی شدهاند (Sayyadzadeh & Esmaeili, 2024) که یکی از آنها سفید ماهی رودخانهای نمک (Squalius namak Khaefi et al., 2016) است. این گونه بومزاد اخیراً شناسایی شده و در نتیجه مطالعات کمی روی آن صورت گرفته است Mouludi-Saleh & Keivany, 2018)). گونه S. namak دارای قدرت باروری بالا و سرعت رشد سریع است .(CABI, 2023) ازنظر تنوع ژنتیکی ماهی سفید رودخانهای حوضه نمک دارای چهار نوکلئوتید متمایز در جایگاه COI در mtDNA است (Mouludi Saleh et al., 2018).
نظر به اینکه ماهی سفید رودخانهای نمک یک گونه بومزاد و یکی از ماهیان مهم اکولوژیکی آبهای داخلی ایران است (Khaefi et al., 2016) و تاکنون هیچگونه مطالعهای روی ساختار و تنوع ژنتیکی جمعیتهای این گونه صورت نگرفته است، این پژوهش میتواند راهگشای تعدادی از سؤالات درخصوص ویژگیهای جمعیتی این ماهی باشد. همچنین این پژوهش به حفظ تنوع زیستی و نیز حفاظت پایدار از این گونه کمک خواهد کرد.
مواد و روشها
منطقه مطالعهشده: منطقه مطالعهشده رودخانه جاجرود، یکی از رودخانههای با جریان دائمی و از مهمترین بومسازگانهای آبی است که در شمالشرق تهران جریان دارد (Naderi et al., 2020). حوضه آبخیز رودخانه جاجرود به مساحت 674 کیلومتر مربع و طول جغرافیایی 51 درجه و 51 دقیقه تا 51 درجه و 22 دقیقه و عرض جغرافیایی 35 درجه و 45 دقیقه تا 36 درجه قرار دارد. طول رودخانه اصلی حوضه جاجرود 40 کیلومتر است (Najimi et al., 2023). این رودخانه در حوضه آبریز نمک واقع است و ابتدا وارد سد لتیان و سپس سد ماملو در پایین دست میشود (Shabanloo et al., 2022).
نمونهبرداری ماهیان: با استفاده از دستگاه الکتروشوکر 725 SAMUS تعداد 103 ماهی سفید رودخانهای نمک برای نمونهبرداری صید شدند (شکل 1). طول هر ایستگاه نمونهبرداری در حدود 90 متر بهصورت سه تکرار 30 متری بود (Shabanloo et al., 2021). نمونهها با استفاده از محلول میخک (10 گرم در هر لیتر آب مقطر) بیهوش شدند. سپس حدود 1 سانتیمتر از باله دمی بهعنوان نمونه جدا شد و به داخل تیوپ حاوی اتانول 96 درصد منتقل شد؛ زیرا هدف کمترین آسیب به سلامت ماهیان بود و در این زمینه عنوان شده است که نمونهبرداری از باله دمی نسبت به بقیه بالهها آسیب کمتری بر ماهیان دارد (Ferguson, 2020). در پایان، بعد از اطمینان از بازیابی قدرت شنای نمونهها، ماهیان صیدشده به زیستگاه خود بازگردانده شدند و تنها تعداد اندکی به آزمایشگاه منتقل شدند.
شکل 1- تصویر از ماهی سفید رودخانهای نمک (Squalius namak) (شعاع بالگان: سفید ماهیان) از رودخانه جاجرود.
Fig. 1. Namak chub (Squalius namak) (Actinopterygii: Leuciscinae) from the Jajrood River
نمونهبرداری در فصل پاییز سال 1400 در طول مسیر رودخانه جاجرود از پاییندست به سمت بالادست رودخانه در 4 ایستگاه انجام شد (جدول 1). ایستگاههای نمونهبرداری بهنحوی انتخاب شدند که تمام تنوع زیستگاهی در دسترس را شامل شوند. نمونهگیری در زیستگاههای مختلف موجود در رودخانه، در سه مسیر مختلف در خلاف جهت جریان انجام شد. همچنین بهمنظور جلوگیری از فرار ماهیان، در تمام ایستگاهها از تورهای پشتیبان در بالادست و پاییندست استفاده شد (Eagderi et al., 2024).
جدول 1- مختصات جغرافیایی ایستگاههای نمونهبرداری و تعداد نمونههای ماهی سفید رودخانهای نمک Squalius namak از رودخانه جاجرود
Table 1. Geographic coordinates of sampling stations and number of samples of Namak chub (Squalius namak) along the Jajrood River
استخراجDNA : با روش فنل - کلروفرم - ایزوآمیل الکل (Sambrook & Russel, 2001) استخراج DNA از بافت بالههای ماهیان انجام گرفت. DNAهای استخراجی پس از افزودن 30 میکرولیتر آب مقطر استریل برای استفاده در مراحل بعدی، یک شب در یخچال نگه داشته و سپس به فریزر 20- درجه سانتیگراد منتقل شدند. بررسی کمیّت و کیفیت DNAهای استخراجی با استفاده از الکتروفورز ژل آگارز 1 درصد و نانو دراپ انجام شد. غلظت DNA تمام نمونهها 100 تا 180 نانوگرم بود و با توجه به نتایج الکتروفورز ژل آگارز، از کیفیت مناسبی برخوردار بودند. واکنش PCR برای تکثیر پنج جفت پرایمر برای این آزمایش استفاده شد. بهدلیل اینکه پرایمر اختصاصی برای این گونه طراحی نشده بود، از پرایمرهای گونههای نزدیک استفاده شد (جدول 2).
جدول 2- مشخصات پرایمرهای (آغازگرهای) استفادهشده برای بررسی ساختار ژنتیکی ماهی سفید رودخانهای نمک Squalius namak
Table 2. The list of PCR primers and sequences used to investigate the genetic structure of the Namak chub (Squalius namak)
|
لوکوس (نام پرایمر) Primer name |
شماره number |
(5'-3') توالی پرایمر (آغازگر) sequences |
دامنه توالی تکرار (Domain) |
اندازه آلل (bp) |
طول (Length)
|
منابع References |
|
Sluc4-F |
1 |
AAGCATTACCCATGCAGAGC |
(AG)26 |
136-198 |
20 |
|
|
Sluc4-R |
AGCTGCAACACAACCTCCAT |
20 |
||||
|
Sluc5-F |
2 |
GAGAAAGAGAGACCAATCCATAGTT |
(TTC) 10 (TTA)5 |
242-358 |
25 |
|
|
Sluc5-R |
CAAAGCAAGCATCAAACCTG |
20 |
||||
|
Sluc7-F |
3 |
GGAAACTGACACATCGCTTG |
(CTTT)16 |
200-298 |
20 |
|
|
Sluc7- R |
CGAAGGACTGGACTGGAAAG |
20 |
||||
|
Sluc13-F |
4 |
CACCCAGGCAATAAACAAGG |
(CA)35 |
228-236 |
20 |
|
|
Sluc13-R |
GGGTTAAGGGTCGGTTTAGG |
20 |
||||
|
SarN2-F |
5 |
GAACAAACATCACTGAAGCACTCT |
(AC)10 |
110-120 |
24 |
|
|
SarN2-R |
ACGTCAGACTTCAGGCATCC |
20 |
واکنش زنجیرهای پلیمراز (PCR) : تکثیر لوکوسهای ریزماهواره با استفاده از واکنشهای زنجیرهای پلیمراز در حجم 15 میکرولیتر شامل دو میکرولیتر DNA، یک میکرولیتر از هر پرایمر، 5/7 میکرولیتر مسترمیکسPCR (Ampliqon, Denmark) و آب مقطر تا رسیدن به حجم انجام شد. چرخه دمایی برای هر لوکوس ژنی عبارت بود از یک چرخة سه دقیقهای در دمای 94 درجه سانتیگراد (واسرشتسازی اولیه) 35 چرخه در 94 درجه به مدت 30 ثانیه (واسرشتسازی)، 72 درجه به مدت یک دقیقه (بسط) و یک چرخه 72 درجهای به مدت سه دقیقه که بهعنوان مرحلة بسط نهایی صورت گرفت.
الکتروفورز عمودی: برای جداسازی بهتر محصولاتPCR ، از الکتروفورز روی ژل پلیاکریل آمید 8 درصد استفاده شد. الکتروفورز عمودی برای هر ژل پلیاکریل آمید با جریان 150 ولت به مدت چهار ساعت انجام شد. درنهایت، ژل از صفحات شیشهای جدا شد و رنگآمیزی توسط نیترات نقره انجام شد (شکل 2) (Mesquita et al., 2003, O'connell & Wright, 1997) و پس از عکسبرداری از ژلها، از نرمافزار AlphaEaseC برای امتیازدهی به باندها و محاسبه وزن مولکولی آنها استفاده شد.
شکل 2 - آرایش باندهای تشکیلشده محصول PCR نشانگرها روی ژل پلیاکریل آمید 8 درصد
Fig. 2. Bands in the 8% polyacrylamide gel produced via PCR using primer 1 for amplification of Sluc4 gene
آنالیز آماری: دادههای حاصل از AlphaEaseC ابتدا به اکسل منتقل شدند و سپس شاخصهای ژنتیکی نظیر تعداد آلل در هر لوکوس (Na)، تعداد آلل مؤثر (Ne)، هتروزیگوسیتی مشاهدهشده (Ho)، هتروزیگوسیتی مورد انتظار (He)، میزان تمایز (Fst) و جریان ژنی (Nm) با استفاده از نرمافزار GeneAlex 6.5 (Peakall & Smouse, 2006) محاسبه شدند. همچنین بهمنظور تعیین تنوع ژنتیکی درون و بین جمعیتی براساس مدل اللی بینهایت (Fst) و مدل جهش پلهای (RST) از آنالیز واریانس مولکولی (AMOVA) بسته نرمافزاری GeneAlex استفاده شد. درنهایت، تجزیه براساس مؤلفههای اصلی (PCA) برای بررسی ارتباط ژنتیکی 4 جمعیت و آزمون تعادل هاردی واینبرگ برای چهار لوکوس بهکاررفته با تصحیح بونفرونی در نرمافزار Cervus 3.0.7اجرا شد.
نتایج
در این مطالعه پنج جفت پرایمر استفاده شد که چهار جفت از آنها پلیمورف (چندشکل) بودند و برای ارزیابی تنوع ژنتیکی، نتایج قابل قبولی را نشان دادند. درمجموع 372 باند قابل ارزشدهی در میان جمعیتهای 4 نقطه بررسیشده مشخص شد. تعداد کل آلل در سطح لوکوس در دامنه 2 تا 37 به دست آمد، لوکوس Sluc5 بیشترین تعداد آلل (37) و لوکوس Sluc4 کمترین آلل (2) را نشان دادند (جدول 3). تعداد متوسط آلل بهازای هر لوکوس در جمعیت یک خجیر 25/10، جمعیت دو خجیر 25/26، جمعیت یک جاجرود 75/11 و جمعیت دو جاجرود 5/1 به دست آمد. میانگین هتروزیگوسیتی مشاهدهشده در جمعیت یک خجیر 65/0، جمعیت دو خجیر 68/0، جمعیت یک جاجرود 51/0 و جمعیت دو جاجرود 75/0 محاسبه شد. مقادیر هتروزیگوسیتی مورد انتظار (He)، در دامنه 50/0 تا 96/0 (میانگین کل چهار جمعیت 74/0) و در سطح مناطق نیز بیشترین و کمترین مقدار متوسط بهترتیب 94/0 و 37/0 متعلق به جمعیتهای 2 خجیر و 2 جاجرود بود.
در سطح لوکوسها میزان تمایز (Fst) و جریان ژنی (Nm) محاسبه شد (جدول 4) و متوسط جریان ژنی (Nm) در بین مناطق 30/1 و بیشترین جریان ژنی در لوکوس Sluc5 (93/1) و کمترین مقدار آن در لوکوس SarN2 (49/0) مشاهده شد. میانگین تمایز ژنتیکی نیز 19/0 بود و بیشترین مقدار تمایز در لوکوس SarN2 (33/0) مشاهده شد.
جدول 3 - تنوع ژنتیکی چهار لوکوس مطالعهشده در چهار جمعیت از ماهی سفید رودخانهای نمک Squalius namak در رودخانه جاجرود.
Table 3. Genetic diversity of four studied loci in four populations of Namak chub (Squalius namak) in the Jajrood River based on microsatellite markers
|
جمعیت Population |
لوکوس Locus |
Na |
Ne |
Ho |
He |
|
جمعیت 1 خجیر |
Sluc4 |
7 |
12/6 |
28/0 |
83/0 |
|
Sluc5 |
16 |
33/13 |
80/0 |
92/0 |
|
|
Sluc13 |
12 |
66/10 |
75/0 |
90/0 |
|
|
SarN2 |
6 |
60/3 |
77/0 |
72/0 |
|
|
|
میانگین |
25/10 |
43/8 |
65/0 |
85/0 |
|
جمعیت 2 خجیر |
Sluc4 |
15 |
50/9 |
42/0 |
89/0 |
|
Sluc5 |
37 |
54/23 |
90/0 |
95/0 |
|
|
Sluc13 |
28 |
78/19 |
60/0 |
94/0 |
|
|
SarN2 |
25 |
75/17 |
77/0 |
94/0 |
|
|
|
میانگین |
25/26 |
64/17 |
68/0 |
94/0 |
|
جمعیت 1 جاجرود |
Sluc4 |
6 |
20/2 |
07/0 |
57/0 |
|
Sluc5 |
17 |
65/11 |
92/0 |
91/0 |
|
|
Sluc13 |
15 |
22/10 |
40/0 |
90/0 |
|
|
SarN2 |
9 |
11/7 |
62/0 |
85/0 |
|
|
|
میانگین |
75/11 |
80/7 |
51/0 |
80/0 |
|
جمعیت 2 جاجرود |
Sluc4 |
2 |
2 |
1 |
50/0 |
|
Sluc5 |
2 |
2 |
1 |
50/0 |
|
|
Sluc13 |
2 |
2 |
1 |
50/0 |
|
|
SarN2 |
0 |
0 |
0 |
0 |
|
|
|
میانگین |
5/1 |
5/1 |
75/0 |
37/0 |
Na: تعداد آلل، Ne: تعداد آلل مؤثر، HO: هتروزیگوسیتی مشاهدهشده، He: هتروزیگوسیتی مورد انتظار
جدول 4- میزان جریان ژنی (Nm) و تمایز (Fst) در سطح چهار لوکوس (جایگاه ژنی) استفادهشده
Table 4. Gene flow (Nm) and genetic differentiation (Fst) rates at the four studied loci
|
لوکوس (جایگاه ژنی) |
Sluc4 |
Sluc5 |
Sluc13 |
SarN2 |
میانگین |
|
Nm |
10/1 |
93/1 |
68/1 |
49/0 |
30/1 |
|
Fst |
18/0 |
11/0 |
12/0 |
33/0 |
19/0 |
آنالیز واریانس مولکولی (AMOVA) در Rst بیانکنندة این بود که 93 درصد تمایز ژنتیکی مشاهدهشده به تفاوت در میان افراد جمعیتهای مختلف مربوط است و 7 درصد نیز بهدلیل تفاوت بین 4 جمعیت مطالعهشده است (جدول 5). همچنین نتایج Fst حاصل از AMOVA در سطح 99 درصد نشان داد در گونه ماهی سفید رودخانهای نمک S.namak)) 42 درصد از تنوع مشاهدهشده مربوط به افراد درون جمعیتهای (Within Indiv) هر ایستگاه بوده و 54 درصد از تنوع به افراد بین جمعیتهای مختلف (Among Indiv) مربوط است و سهم تنوع در میان دو منطقه (خجیر و جاجرود) و چهار جمعیت (Among Pops) تنها 2 درصد است (شکل 3).
جدول 5- آنالیز واریانس مولکولی AMOVA)) برحسب Rst
Table 5. Analysis of molecular variance (AMOVA) in terms of Rst
|
df |
P |
درصد |
|
|
Among Pops |
2 |
01/0 |
7 |
|
Within Indiv |
103 |
01/0 |
93 |
|
Total |
205 |
01/0 |
100 |
شکل 3 - چگونگی توزیع تنوع ژنتیکی مشاهدهشده در ماهی سفید رودخانهای نمک Squalius namak تحت معیار Fst
Fig. 3. genetic diversity distribution among and within subgroups in the Namak chub (Squalius namak) based on estimated Fst values
نمودار دوبعدی تحلیل مؤلفههای اصلی (PCA) برای بررسی روابط بین نمونههای ماهی سفید رودخانهای نمک S.namak))، نشان داد ماهیان ایستگاههای 1 خجیر و 2 خجیر و همچنین ماهیان ایستگاههای 1 جاجرود و 2 جاجرود همپوشانی بالایی دارند؛ اما ایستگاههای 1خجیر و 1 جاجرود تفکیکپذیر هستند (شکل 4).
شکل 4- رابطه ژنتیکی میان چهار گروه مطالعهشده از ماهیان سفید رودخانهای نمک با استفاده از تجزیه مؤلفههای اصلی (PCA)
Fig. 4. Genetic relationship between the four studied groups of Namak chub using principal component analysis (PCA)
آزمون مربع کا همراه با اعمال ضریب تصحیح بونفرونی برای بررسی تعادل هاردی - واینبرگ نشان داد تمام لوکوسها انحراف از تعادل هاردی - واینبرگ داشتند (جدول 6).
جدول 6- بررسی تعادل هاردی - واینبرگ در 4 لوکوس ریزماهواره برای جمعیتهای نمونهبرداریشده ماهی Squalius namak در رودخانه جاجرود.
Table 6. Hardy-Weinberg equilibrium analysis for microsatellite loci applied to studied groups of the Namak chub (Squalius namak) in the Jajrood River
|
نام لوکوس |
مربع کا (خی2) |
درجه آزادی |
مقدار P |
|
Sluc4 |
1979 |
666 |
00/0 |
|
Sluc5 |
3977 |
3655 |
00/0 |
|
Sluc13 |
3572 |
2145 |
00/0 |
|
SarN2 |
2394 |
1711 |
00/0 |
بحث
هتروزیگوسیتی در مطالعه ساختار جمعیت گونهها ارزش بسیار دارد؛ زیرا هر هتروزیگوت ناقل آللهای متفاوتی بوده که نشاندهندة گوناگونی است (Ghodsi et al., 2011). در این بررسی میانگین هتروزیگوسیتی مشاهدهشده در چهار جمعیت 64/0 بود که نسبت به مقادیر مشاهدهشده در ماهیان آب شیرین (46/0) بالاتر است (Hedayati et al., 2017)؛ اما براساس نتایج تحلیل واریانس مولکولی بین جمعیتهای دو ناحیه بررسیشده (خجیر و جاجرود) از این نظر تفاوت معنیداری وجود نداشت (P<0.05). در صورتی که مقادیر بهدستآمده برای هتروزیگوسیتی بالاتر از حد متوسط باشد، میزان تنوع ژنتیکی نیز بالا است و اگر پایینتر باشد، میزان تنوع ژنتیکی در آن جمعیت کاهش پیدا کرده است (Farasati et al., 2020). ریزماهوارهها بهعلت بالابودن تعداد آللهایشان، در بین تمام نشانگرها، بالاترین میزان هتروزیگوسیتی را نشان میدهند (Karami nasab et al., 2014). ممکن است در جمعیتهای طبیعی، هتروزیگوسیتی از طریق همخونی، نمونهبرداری غیرتصادفی، ساختار درون جمعیت، رانش ژنتیکی، آلل نول، فشار صید بیرویه یا ترکیبی از عوامل فوق رخ دهد (Shabani & Kolangi Miandare, 2014). در این مطالعه، مقدار متوسط هتروزیگوسیتی مشاهدهشدهHo) ) در سه جمعیت از مقدار متوسط هتروزیگوسیتی مورد انتظار (He) پایینتر بود که این امر نشاندهندة روند کاهشی تنوع ژنتیکی و وجود احتمالی درونآمیزی در افراد جمعیتهای ماهیان S.namak است. بهطور کلی دلایل زیستشناختی کاهش هتروزیگوسیتی بهخوبی شناخته نشده است و آللهای تکثیرنشده یا نول، از مهمترین عوامل ایجادکنندة کاهش در جایگاههای ژنی ریزماهواره هستند (Shabani et al., 2012). بهعنوان یک اصل، در صورتی که مقادیر هتروزیگوسیتی در محدوده 3/ تا 8/0 باشد، نشانگرهای ریزماهواره میتوانند در نشاندادن تنوع ژنتیکی مفید واقع شوند (Farasati et al., 2020).
در بررسیهای تنوع ژنتیکی، غنای آللی نسبت به هتروزیگوسیتی دارای ارزش بالاتری است. درواقع، بالابودن غنای آللی، نشاندهندة بالابودن اندازة مؤثر جمعیت بوده و برای ارزیابی تنوع ژنتیکی در جمعیتهایی مناسبتر است که برای برنامههای بهگزینی یا حفاظت انتخاب شدهاند.(Karami nasab et al., 2014) در مطالعه حاضر مجموع آللهای موجود برای تعداد چهار لوکوس در جمعیت یک خجیر 41، در جمعیت دو خجیر 105، در جمعیت یک جاجرود 47 و در جمعیت دو جاجرود 6 آلل بود. میانگین کل تعداد آللها برای چهار جمعیت برابر 43/12 بود که این مقدار بیشتر از مقدار محاسبهشده برای ماهیان آب شیرین (5/7) بود (Hedayati et al., 2017). این مقدار بهشدت تأثیرگرفته از تعداد نمونهها است؛ بر همین اساس، این امکان وجود دارد که در آزمایشهای گوناگون با تعداد نمونههای متفاوت، تعداد آللهای واقعی مختلفی برای یک لوکوس معین به دست آید. پتانسیل بیولوژیکی هر جمعیت، در حقیقت به آللهای موجود در آن جمعیت بستگی دارد (Gorjipoor & Nazari, 2014). تغییراتی که تعداد آللها میتوانند از نسلی به نسل دیگر داشته باشند، در جمعیت کوچک سریعتر از جمعیت بزرگ رخ میدهد؛ بنابراین، چنین احتمالی وجود دارد که در نسلهای بعدی غنای آللی کاملاً متفاوت با وضعیت کنونی باشد (Hosseini et al., 2016).
یکی از ویژگیهای مشخص در مطالعات ژنتیکی جمعیت، جریان ژنی در میان زیرجمعیتها است. با سطوح بالای مهاجرت و جریان ژنی بین جمعیت، میزان شباهت جمعیتها افزایش مییابد (Shabani et al., 2012). درواقع میزان جریان ژنی (Nm) به تعداد مولدین مهاجر از یک منطقه به منطقه دیگر گفته میشود؛ هرچه این میزان بین دو منطقه بیشتر باشد، به این معنی است که مهاجرت بین دو منطقه، بیشتر و اختلاف ژنتیکی، کمتر و میزان تنوع ژنتیکی در دو منطقه، بیشتر میشود (Silavi et al., 2024). در این مطالعه میزان متوسط جریان ژنی 30/1 بود و دلیل مشاهده تنوع ژنتیکی بالا در این گونه را میتوان وجود جریان ژنی بین جمعیتها دانست. هرگاه 1< Nm باشد، جریان ژنی اصلیترین عامل ایجاد تمایز ژنتیکی است و هرگاه 1> Nm باشد، رانش ژنی عامل اصلی ایجاد تمایز ژنتیکی میشود (Karami nasab et al., 2014). مقدار متوسط جریان ژنی میتواند ناشی از مهاجرت طبیعی تغذیهای و تولیدمثلی بین مناطق باشد (Shabani et al., 2012). در عین حال، بالابودن مقادیر جریان ژنی نمیتواند تنوع ژنتیکی را تضمین کند؛ زیرا پایینبودن تعداد مولدین و جمعیتها و در نتیجه آمیزش بین جمعیتهای کوچک میتواند تنوع را کاهش دهد (Farasati et al., 2020).
فاکتور Fst توصیفکنندة تمایز جمعیتها در سطوح مختلف ساختار ژنتیکی است (Ghasemi & Pourjam, 2021). درواقع Fst تفاوت فراوانی آللی لوکوسهای مطالعهشده در بین جمعیتهای مختلف را ارزیابی میکند و مقادیر عددی این شاخص از صفر تا یک متغیر است. Fst برابر با صفر نشاندهندة عدم وجود تفاوت در بین جمعیتها، آللهای یکسان و فراوانی یکسان است (Hosseini et al., 2016). طبق یافتههای مطالعه حاضر مقدار متوسط Fst برای چهار لوکوس برابر 19/0 بود که نشاندهندة تمایز نسبتاً بالا است. بهطور کلی هرگاه میزان Fst کمتر از 05/0 باشد، نشاندهندة وجود تمایز کم، مقدار بین 05/0 تا 15/0 نشاندهندة تمایز متوسط، مقدار بین 15/0 تا 25/0 نشاندهندة تمایز بالا و مقدار بالای 25/0 نیز نشاندهندة تمایز ژنتیکی خیلی زیاد در بین جمعیتها و جدایی کامل جمعیتهاست (Karami nasab et al., 2015)؛ با این حال، عدد بهدستآمده در این مطالعه نمایانگر همه جمعیت حقیقی نیست. نکته دیگر اینکه میزان Fst در اکثر موارد به یک نمیرسد؛ زیرا اثر پلیمورفیسم (ناشی از جهش) بهطور مؤثری میزان Fstرا کاهش میدهد (Silavi et al., 2024). با تبادل افراد، تبادل ژنها نیز پیش میآید و تبادل بیشتر منجر به کمشدن تمایز ژنتیکی بین جمعیتها میشود. مهاجرت زیاد، از جدایی ژنتیکی جمعیتها جلوگیری میکند و در ماهیان بین مقدار Fst و قابلیت پراکنش همبستگی منفی وجود دارد (Norouzi et al., 2013).
آنالیز واریانس مولکولی (AMOVA) بهعنوان روشی مناسب در تعیین ساختار و تمایز ژنتیکی بین جمعیتها مطرح است (Shabani et al., 2016). براساس نتایج این مطالعه، آنالیز واریانس مولکولی تحت معیار Rst و Fst بیانکنندة این بود که بخش عمده تفاوت ژنتیکی مربوط به افراد جمعیتهاست (Rst=93% , Fst=96% ) و تفاوت بین جمعیتها و دو ناحیه اندک است. برای ریزماهوارهها میزان Rst میتواند بیشتر از Fst باشد. با توجه به اینکه Rstاز دانستههای مربوط به اندازه آللی استفاده میکند و وابسته به جهش نیست، میتواند دادههای بیولوژیک بهتری نسبت به معیار Fstفراهم کند (Karami nasab et al., 2015). اختلاف درون جمعیتی ممکن است بهدلیل تفاوت ژنتیکی بین افراد باشد. در الگوی کلی ساختار ژنتیک جمعیت، اختلاف ژنتیکی درون جمعیتی در ماهیان، بیشتر از اختلاف ژنتیکی بین جمعیتهاست که این امر ممکن است ناشی از موانع اکولوژیکی باشد که باعث جلوگیری از پراکندگی میشوند (Ghasemi & Pourjam, 2021). بالاتربودن تنوع درون جمعیتی نسبت به بین جمعیتی نشان میدهد در بین جمعیتهای مختلف ساختار ژنتیکی بارزی وجود ندارد (Ghodsi et al., 2011). با تشدید توسعه شهرنشینی و نواحی صنعتی، بعضی از زیستگاههای تخمریزی، مکانهای نوزادگاهی و تغذیهای ماهیان از بین میروند و این امر میتواند جمعیتهای متفاوتی را در بعضی از زیستگاههای جداشده درون یک ناحیه با فاصله جغرافیایی کوتاه به وجود آورد (Archangi et al., 2015).
نتایج PCA نشان دادند جمعیتهای دو خجیر و دو جاجرود کمترین همپوشانی را با هم دارند و بنابراین براساس شاخصههای ژنتیکی گفته میشود ماهی سفید رودخانهای نمک در رودخانه جاجرود دارای جمعیتهای مجزا است و به نظر میرسد جریان ژنی بین جمعیتها متأثر از فاصله جغرافیایی، تمایز جمعیتها را افزایش دهد. انحراف از تعادل هاردی - واینبرگ که نشاندهندة بهگزینی، اختلاط جمعیتها یا جفتگیری غیرتصادفی است، ممکن است در جمعیتهای بسیاری از ماهیان وحشی رخ دهد؛ بنابراین، بررسی این موضوع جزء اولین مراحل در مطالعات مربوط بـه ساختار جمعیتها است (Shabani et al., 2016). در بررسیهای تعادل هاردی - واینبرگ در این مطالعه تمامی لوکوسها از تعادل خارج بودند. بهطور کلی در جمعیتهای ماهیان، انحراف از تعادل هاردی - واینبرگ زیاد دیده میشود که ضریب تصحیح انحراف از تعادل هاردی - واینبرگ میتواند علل متعددی داشته باشد (Karami nasab et al., 2015). این انحراف از تعادل میتواند به کاهش اندازه جمیعت، وجود رابطههای خویشاوندی، وجود آللهای خنثی و در نتیجه افزایش کاذب هموزیگوسیتی مرتبط باشد (Farasati et al., 2020). انحراف از تعادل هاردی - واینبرگ را همچنین میتوان به تعداد کم نمونهها و خطای نمونهبرداری (نمونهبرداری از افراد خویشاوند) نسبت داد (Askari et al., 2014). در اینجا علاوه بر دلایل ذکرشده، انحراف از تعادل را میتوان به غیراختصاصیبودن پرایمرها یا تعداد کم نمونه نیز نسبت داد. با توجه به نبود ثبت شجره ماهیان، باید مطالعات بیشتری را در آینده در دستور کار قرار داد (Gorjipoor and Nazari, 2014).
براساس آنالیزهای موجود، به نظر میرسد گونه S.namak دارای تنوع ژنتیکی مطلوبی در مناطق بررسیشده است و با توجه به اهمیت اکولوژیک این گونه بومزاد در رودخانهها حفظ تنوع ژنتیکی آن لازم و ضروری به نظر میرسد. همچنین نتایج نشان دادند نشانگر ریزماهواره از توانایی بالایی برای نشاندادن میزان تنوع ژنتیکی در این ماهی برخوردار است. پیشنهاد میشود بهمنظور دستیابی به اطلاعات دقیقتر در ارتباط با این گونه با ارزش بومزاد ایران، در آینده مطالعاتی روی ژنوم میتوکندریایی منابع این گونه صورت گیرد.
تشکر و قدردانی
این مطالعه با حمایت مالی دانشگاه تهران انجام شده است.