Genetic diversity of Namak chub: Squalius namak (Actinopterygii: Leuciscinae) in the Jajrood River using microsatellite markers

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

The Namak chub (Squalius namak) is an endemic species of the inland waters of Iran. However, no study has been conducted on the structure and genetic diversity of its populations. To address this, 103 Namak chub (S. namak) were sampled from four stations in the Jajrood River, and five microsatellite loci (Sluc4, Sluc5, Sluc7, Sluc13, and SarN2) were used. GenAlEx and CERVUS genetic software were used to analyze the data. The number of alleles ranged from 2 to 37 (average 12.43), and expected heterozygosity ranged from 0 to 0.95 (average 0.74), both of which were higher than the averages typically reported for freshwater fish. The average FST index for genetic differentiation was estimated to be 0.19, and the mean gene flow (Nm) was 1.30. Analysis of molecular variance revealed high genetic diversity within populations and low diversity among them. All loci exhibited a pattern of deviation from Hardy–Weinberg equilibrium (HWE), likely due to the use of non-specific primers, the presence of null alleles, and possibly small sample sizes. According to the findings of the present study, the genetic diversity observed in the populations of the Namak chub at these four sites was high, and considering the relatively low level of gene flow and the high FST, there appears to be significant differentiation among the populations of this species in the studied areas. These results can support efforts in biodiversity conservation and management planning.
Keywords: Population genetics, Genetic differentiation, Microsatellite, Chub, Molecular marker.
 
Introduction
The genetic diversity of a species determines its adaptive capacity and evolutionary potential. Small populations of species with restricted distributions often have low genetic diversity within populations but high genetic differentiation between populations due to genetic drift and limited gene flow (Zhai et al., 2019). Genetic diversity is now a concern in many conservation programs, as low genetic diversity is associated with increased inbreeding, accumulation of deleterious mutations, and reduced adaptive potential (McCusker & Bentzen, 2010). Microsatellite loci are standard genetic markers used for population genetic analysis (Coates et al., 2009). In freshwater ecosystems, fish biodiversity is a good indicator of ecosystem health (Zhou et al., 2024). Freshwater fish species constitute 40% of the world's total fish diversity and are one of the most endangered animal groups (Costa et al., 2021). In terms of diversity and endemism, Iran has a high richness of freshwater fish. According to the latest studies, a total of 300 fish species belonging to 38 families and 110 genera live in the inland waters of Iran (Sayyadzadeh & Esmaeili, 2024). The genus Squalius is one of the important genera in the family Leuciscidae, with five species of this genus identified in the inland waters of Iran (Sayyadzadeh & Esmaeili, 2024), one of which is the Namak chub (Squalius namak Khaefi et al., 2016). This endemic species was recently described and as a result, few studies have been conducted on it (Mouludi-Saleh & Keivany, 2018).
 
Materials & Methods
The study area was the Jajroud River, one of the important aquatic ecosystems in the northeast of Tehran (Naderi et al., 2020). Using a Samus 725 electrofishing device, 103 S. namak were collected. Fish were anesthetized using a clove oil solution. Then, caudal fin clips were obtained from the samples and transferred into the labeled tubes containing 96% ethanol. DNA was extracted from the fin tissue using the phenol-chloroform method (Sambrook & Russel, 2001). PCR reaction was used to amplify five pairs of primers for this work. Since no specific primer was designed for this species, the primers from related species were used. For better separation of PCR products, electrophoresis on 8% polyacrylamide gel was used. Data obtained from AlphaEaseC software were first transferred to Excel. Then genetic indices such as number of alleles per locus (Na), number of effective alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), degree of differentiation (FST) and gene flow (Nm) were calculated using GenAlEx 6.5 software (Peakall & Smouse, 2006). Also, to determine the genetic diversity within and between populations based on FST and RST, the analysis of molecular variance (AMOVA) in the GenAlEx software package was used. Finally, principal component analysis (PCA) was performed to examine the genetic relationship of the four populations and the Hardy-Weinberg equilibrium test for the four loci used with Bonferroni correction in CERVUS 3.0.7 software.
 
Research Findings
The total number of alleles at the loci level ranged from 2 to 37, with the Sluc5 locus showing the highest number of alleles (37) and the Sluc4 locus showing the lowest number (2). The average number of alleles per locus was 10.25 in the Khojir 1 population, 26.25 in the Khojir 2 population, 11.75 in the Jajrood 1 population, and 1.5 in the Jajrood 2 population. The average observed heterozygosity in the Khojir 1 population was 0.65, in the Khojir 2 population was 0.68, in the Jajrood 1 population was 0.51, and in the Jajrood 2 population was 0.75. Expected heterozygosity values (He) ranged from 0.50 to 0.96 (average across all four populations: 0.74). At the population level, the highest and lowest average values were 0.94 and 0.37, observed in the Khojir 2 and Jajrood 2 populations, respectively. At the locus level, the degree of differentiation (FST) and gene flow (Nm) were calculated, and the average gene flow among the populations was 1.30, with the highest value observed at the Sluc5 locus (1.93) and the lowest at the SarN2 locus (0.49). The average genetic differentiation was 0.19, with the highest differentiation observed at the SarN2 locus (0.33). AMOVA based on RST indicated that 93% of the observed genetic variation was related to differences among individuals within populations, and 7% was due to differences among the four studied populations. The chi-square test with Bonferroni correction, used to examine Hardy–Weinberg equilibrium, showed that all loci deviated from equilibrium.
 
Discussion of Results & Conclusion
In this study, the average heterozygosity observed in the four populations was 0.64, which is higher than the values observed in freshwater fish (0.46) (Hedayati et al., 2017). If heterozygosity values are higher than the average, the level of genetic diversity is considered high; if they are lower, genetic diversity in that population is reduced (Farasati et al., 2020). The average total number of alleles for the four populations was 12.43, which was higher than the value calculated for freshwater fish (7.5) (Hedayati et al., 2017). This value is strongly affected by the number of samples (Gorjipoor & Nazari, 2014). In the present study, the average gene flow rate was 1.30, and the reason for the high genetic diversity observed in this species can be attributed to the existence of gene flow between populations. When Nm < 1, gene flow is the dominant factor causing differentiation, and when Nm > 1, genetic drift is the main factor causing genetic differentiation (Karami Nasab et al., 2014). However, high gene flow values do not necessarily guarantee genetic diversity, because a low number of parents and populations, and consequently interbreeding among small populations, can reduce diversity (Farasati et al., 2020). However, the number obtained in this study does not represent all real populations. Another point is that the FST rate does not reach unity in most cases because the effect of polymorphism (caused by mutations) effectively reduces the FST rate (Silavi et al., 2024). Based on the available analyses, it appears that the S. namak species exhibits favorable genetic diversity in the studied areas. Given the ecological importance of this endemic species in rivers, preserving its genetic diversity is deemed necessary and important. 

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

(Gigliarelli et al., 2012)

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

(Mesquita et al., 2003)

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 دارای تنوع ژنتیکی مطلوبی در مناطق بررسی‌شده است و با توجه به اهمیت اکولوژیک این گونه بومزاد در رودخانه‌ها حفظ تنوع ژنتیکی آن لازم و ضروری به نظر می‌رسد. همچنین نتایج نشان دادند نشانگر ریزماهواره از توانایی بالایی برای نشان‌دادن میزان تنوع ژنتیکی در این ماهی برخوردار است. پیشنهاد می‌شود به‌منظور دستیابی به اطلاعات دقیق‌تر در ارتباط با این گونه با ارزش بومزاد ایران، در آینده مطالعاتی روی ژنوم میتوکندریایی منابع این گونه صورت گیرد.

 

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

این مطالعه با حمایت مالی دانشگاه تهران انجام شده است.

Archangi, B., Maturian, H., Salari, M. A., & Ronagh, M. T. (2015). Genetic Stock Assessment of Tenualosa ilisha of Persian Gulf Waters and Western Coast of Malaysia Using MtDNA Genome Sequencing. Journal of Oceanography, 6(23), 19-25.‏ http://joc.inio.ac.ir/article-1-846-fa.html [In Persian]
Askari, G., Shabani, A., & Rezaei, H. (2014). Study of genetic diversity of Garra rufa (Heckel, 1843) in Fars Province using microsatellite markers. Modern Genetics Journal, 8(4), 341-348.‏ https://mg.genetics.ir/article-1-1222-fa.html [In Persian]
Bagherian nejad, S. (2014). Fonestical study of freshwater Fishes from Poldoukhtar county [Unpublished master thesis]. Arak University. [In Persian]
Bickham, J. W., Sandhu, S., Hebert, P. D., Chikhi, L., & Athwal, R. (2000). Effects of chemical contaminants on genetic diversity in natural populations: implications for biomonitoring and ecotoxicology. Mutation research/Reviews in Mutation research, 463(1), 33-51.‏ https://doi.org/10.1016/S1383-5742(00)00004-1
CABI (Centre for Agriculture and Biosciences International), (2023). Squalius (European chub). CABI press.
Coates, B. S., Sumerford, D. V., Miller, N. J., Kim, K. S., Sappington, T. W., Siegfried, B. D., & Lewis, L. C. (2009). Comparative performance of single nucleotide polymorphism and microsatellite markers for population genetic analysis. Journal of Heredity, 100(5), 556-564.‏ https://doi.org/10.1093/jhered/esp028
Costa, M.J., Duarte, G., Segurado, P. and Branco, P. (2021). Major threats to European freshwater fish species. Science of the Total Environment, 797, 149105. https://doi.org/10.1016/j.scitotenv.2021.149105 .
Dudgeon, D. (2010). Prospects for sustaining freshwater biodiversity in the 21st century: linking ecosystem structure and function. Current Opinion in Environmental Sustainability, 2(5-6), 422-430. https://doi.org/10.1016/j.cosust.2010.09.001
Dudgeon, D., Arthington, A. H., Gessner, M. O., Kawabata, Z. I., Knowler, D. J., Lévêque, C., ... & Sullivan, C. A. (2006). Freshwater biodiversity: importance, threats, status and conservation challenges. Biological reviews, 81(2), 163-182.‏ https://doi.org/10.1017/S1464793105006950
Eagderi, S., Zamani-Faradonbeh, M., Mouludi-Saleh, A., & Poorbagher, H. (2024). Comparison of habitat suitability of Namak scraper, Capoeta buhsei Kessler, 1877 in autumn and winter seasons in the Jajroud River, Namak Lake basin. Environmental Researches, 14(28), 41-54.‏ https://doi.org/10.22034/eiap.2024.191702 [In Persian]
Ellegren, H., & Galtier, N. (2016). Determinants of genetic diversity. Nature Reviews Genetics, 17(7), 422-433.‏ https://doi.org/10.1038/nrg.2016.58
Farasati, S., Khoshkholgh, M., & Yarmohammadi, M. (2020). Genetic diversity of cultured beluga sturgeon (Huso huso) brood stocks by using Microsatellite method. Aquatic Physiology and Biotechnology, 8(8), 55-72.‏ https://doi.org/10.22124/japb.2020.12704.1319 [In Persian]
Faulks, L. K., Kerezsy, A., Unmack, P. J., Johnson, J. B., & Hughes, J. M. (2017). Going, going, gone? Loss of genetic diversity in two critically endangered Australian freshwater fishes, Scaturiginichthys vermeilipinnis and Chlamydogobius squamigenus, from Great Artesian Basin springs at Edgbaston, Queensland, Australia. Aquatic Conservation: Marine and Freshwater Ecosystems, 27(1), 39-50.‏ https://doi.org/10.1002/aqc.2684
Ferguson, R. (2020). A comparative analysis of lethal and non-lethal fish sampling techniques‏ [Master thesis, North Carolina State University]. https://B2n.ir/nz1408
Ghasemi, S. A., & Pourjam, F. (2021). Population Genetic Structure and Genetic Diversity of Cobia (Rachycentron canadum) in the Persian Gulf and Makran Sea. Taxonomy and Biosystematics, 13(48), 33-46.‏ https://doi.org/10.22108/tbj.2021.128640.1160 [In Persian]
Ghodsi, Z., Shabani, A., & Shabanpour, B. (2011). Genetic diversity of Liza aurata (Risso, 1810) in the coastal regions of Golstan province, using microsatellite marker. Taxonomy and Biosystematics, 3(6), 35-46.‏ https://tbj.ui.ac.ir/article_17400.html [In Persian]
Gigliarelli, L., Puletti, M., Giannetto, D., Franchi, E., Lanfaloni, L., Panara, F., Lorenzoni, M. and Lucentini, L. (2012) Isolation of microsatellite markers in Squalius lucumonis (Bianco, 1983) and cross-species amplification within the family Cyprinidae and other freshwater fish species. Italian journal of zoology, 79(2), 169-174. https://doi.org/10.1080/11250003.2011.642900
Gorjipoor, E., & Nazari, S. (2014). Genetic variation in farmed population of Rainbow trout (Onchorhynchus mykiss) based on microsatellite markers in Iran. Iranian Scientific Fisheries Journal, 22(4), 93-106.‏ http://isfj.ir/article-1-1174-fa.html [In Persian]
Hedayati, S. A. A., Askari, G., Shabani, A., & Soltanifar, Z. (2017). Genetic comparison of Garra rufa (Heckel, 1843) in Shahpour (Kazeron) and Berim (Gachsaran) rivers by using microsatellite markers. Cellular and Molecular Research (Iranian Journal of Biology), 30(1), 106-116.‏ https://cell.ijbio.ir/article_639.html [In Persian]
Hosseini, S. A, Shabani, A., & Rezaei, H.R. (2016). Genetic diversity of zagros zebrafish (Aphanius vladykovi) in Madar-o-dokhtar spring and Chehelgazi spring in the Chaharmahall-o-Bakhtiari Province by using microsatellite markers. MGj, 11(1), 91-99. https://mg.genetics.ir/article-1-1407-fa.html [In Persian]
Kalia, R. K., Rai, M. K., Kalia, S., Singh, R., & Dhawan, A. K. (2011). Microsatellite markers: an overview of the recent progress in plants. Euphytica, 177(3), 309-334.‏ https://doi.org/10.1007/s10681-010-0286-9
Karami Nasab, M, Hosseinnia Z, Kolangi Miandare H, Shabany A. (2014). Genetic diversity of Barbus grypus in the Karkheh River in Khuzestan Province and cultured fish studied using a microsatellite marker. Genetic Engineering and Biosafety Journal, 3(2), 114-123.‏ https://gebsj.ir/article-1-171-fa.html [In Persian]
Karami Nasab, M., Shabani, A., Kolangi Miandare, H., & Sharbaty, S. (2015). Genetic diversity of Barbus grypus (Heckel, 1843) in Karoun and Dez rivers in Khuzestan province by using microsatellite marker. JAIR, 2(1), 63-74. http://jair.gonbad.ac.ir/article-1-177-en.html [In Persian]
Khaefi, R., Esmaeili, H. R., Sayyadzadeh, G., Geiger, M. F., & Freyhof, J. (2016). Squalius namak, a new chub from Lake Namak basin in Iran (Teleostei: Cyprinidae). Zootaxa, 4169(1), 145-159.‏ https://doi.org/10.11646/zootaxa.4169.1.7
Manel, S., Guerin, P. E., Mouillot, D., Blanchet, S., Velez, L., Albouy, C., & Pellissier, L. (2020). Global determinants of freshwater and marine fish genetic diversity. Nature communications11(1), 692.‏ https://doi.org/10.1038/s41467-020-14409-7
McCusker, M. R., & Bentzen, P. (2010). Positive relationships between genetic diversity and abundance in fishes. Molecular Ecology, 19(22), 4852-4862.‏ https://doi.org/10.1111/j.1365-294X.2010.04822.x
Mesquita, N., Cunha, C., Hänfling, B., Carvalho, G., Zé‐Zé, L., Tenreiro, R. and Coelho, M. (2003) Isolation and characterization of polymorphic microsatellite loci in the endangered Portuguese freshwater fish Squalius aradensis (Cyprinidae). Molecular Ecology Notes, 3(4), 572-574. https://doi.org/10.1046/j.1471-8286.2003.00515.x
Mouludi Saleh, A., Keivany, Y., and Jalali, S. A. H. (2018). Biometry of Chub (Squalius namak Khaefi et al., 2016) in rivers of Namak Basin. Experimental animal Biology, 7(1), 107-118.‏ https://eab.journals.pnu.ac.ir/article_4927.html [In Persian]
Mouludi-Saleh, A., & Keivany, Y. (2018). Morphometric analysis of Squalius namak Khaefi et al. 2016 in Khaznagh and Ghare-Chai rivers. Sri Lanka Jouranl of Aquatic Science, 23(2), 173-178.‏ http://doi.org/10.4038/sljas.v23i2.7558
Mukhopadhyay, T., & Bhattacharjee, S. (2014). Study of the Genetic Diversity of the Ornamental Fish Badis badis (HamiltonBuchanan, 1822) in the Terai Region of Sub‐Himalayan West Bengal, India. International Journal of Biodiversity, 2014(1), 791364.‏ https://doi.org/10.1155/2014/791364
Naderi, M., Alioghli, S., Jahandideh, O., Rajabizadeh, Y., & Salarijazi, M. (2020). Determination of Optimal and Desirable Environmental Flow Release from Latian Dam reservoir with Consideration of Ecohydraulic, Hydrological and Hydromorphological Characteristics to Protect the Habitat of the Jajrood River. Iranian Journal of Irrigation & Drainage, 14(4), 1277-1300.‏ https://idj.iaid.ir/article_114898.html [In Persian]
Najimi, F., Aminnejad, B., & Nourani, V. (2023). Assessment of Climate Change’s Impact on Flow Quantity of the Mountainous Watershed of the Jajrood River in Iran Using Hydroclimatic Models. Sustainability, 15(22), 15875.‏ https://doi.org/10.3390/su152215875
Norouzi, M., Nazemi, A., Pourkazemi, M., Samiei, M. H., Daneshvar, F., & Amirjanati, A. (2013). Genetic variability and differentiation of common Kilka fish (Clupeonella cultriventris Nordmann, 1840) in the southern coasts of Caspian Sea. Iranian Scientific Fisheries Journal, 22(1), 139-148.‏ http://isfj.ir/article-1-879-fa.html [In Persian]
O'connell, M. and Wright, J.M. (1997) Microsatellite DNA in fishes. Reviews in fish biology and fisheries, 7(3), 331-363. https://doi.org/10.1023/A:1018443912945.
Ouborg, N. J., Pertoldi, C., Loeschcke, V., Bijlsma, R. K., & Hedrick, P. W. (2010). Conservation genetics in transition to conservation genomics. Trends in genetics, 26(4), 177-187.‏ https://doi.org/10.1016/j.tig.2010.01.001
Peakall, R. O. D., & Smouse, P. E. (2006). GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular ecology notes, 6(1), 288-295.‏ https://doi.org/10.1111/j.1471-8286.2005.01155.x
Pimm, S. L., Russell, G. J., Gittleman, J. L., & Brooks, T. M. (1995). The future of biodiversity. Science, 269(5222), 347-350.‏ https://doi.org/10.1126/science.269.5222.347
Sambrook, J., & Russel DW. (2001). Rapid isolation of yeast DNA. In Molecular cloning, a laboratory manual (pp. 631-632). Cold Spring Harbor Laboratory.
Sayyadzadeh, G., & Esmaeili, H. R. (2024). Freshwater lamprey and fishes of Iran: Reappraisal and updated checklist with a note on Eagderi et al. (2022). Zootaxa, 5402(1), 1-99.‏ https://doi.org/10.11646/zootaxa.5402.1.1
Selkoe, K. A., & Toonen, R. J. (2006). Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers. Ecology letters, 9(5), 615-629.‏ https://doi.org/10.1111/j.1461-0248.2006.00889.x
Shabani, A., & Kolangi Miandare, H. (2014). Genetic structure comparison of Paraschistura nielseni (Nalbant and Bianco, 1998 (of Shapour, Dalaki (Fars province) and Mand (Busher province) Rivers using microsatellite markers. Journal of Aquatic Ecology, 4(2), 71-79.‏ https://jae.hormozgan.ac.ir/article-1-31-fa.html [In Persian]
Shabani, A., Ghodsi, Z., & Naderi, L. (2012). Comparison of genetic diversity in two species of Golden mullet (Risso., 1810) (Liza aurata) and small mullet (Liza saliens) (Risso., 1810) in the southern Caspian Sea, using microsatellite markers. Journal of Utilization and Cultivation of Aquatics, 1(4), 81-93. https://japu.gau.ac.ir/article_1345.html?lang=en [In Persian]
Shabani, A., Kashiri, H., & Ghodsi, Z. (2016). Genetic diversity of Golden Mullet Liza aurata (Risso, 1810) in Gharesou and Gomishan regions and its comparison with dead samples using microsatellite markers. Journal of Applied Ichthyological Research, 4(3), 15-28.‏ http://jair.gonbad.ac.ir/article-1-260-en.html [In Persian]
Shabanloo, H., Poorbagher, H., & Eagderi, S. (2021). Effects of environmental parameters on morphological traits of Squalius namak in the Jajrood River. Aquaculture Sciences, 9(1), 172-181.‏ https://www.aquaculturesciences.ir/article_137435.html [In Persian]
Shabanloo, H., Poorbagher, H., & Eagderi, S. (2022). Effect of size on the relationship between habitat suitability index and niche overlap in Namak chub (Squalius namak Khaefi et al., 2016) in Jajrood River, Namak Lake basin. Iranian Journal of Ecohydrology, 9(4), 751-759.‏ https://doi.org/10.22059/ije.2022.343012.1640 [In Persian]
Silavi, M., Rajabzadeh Ghatrami, E., & Nazemroaya, S. (2024). Determining the genetic diversity and inbreeding coefficient of cultured Anzeh (Luciobarbus esocinus) population using fluorescently labeled microsatellite markers (In press). Journal of Marine Science and Technology. https://doi.org/10.22113/jmst.2024.447140.2587 [In Persian]
Zhai, D. D., Li, W. J., Liu, H. Z., Cao, W. X., & Gao, X. (2019). Genetic diversity and temporal changes of an endemic cyprinid fish species, Ancherythroculter nigrocauda, from the upper reaches of Yangtze River. Zoological research, 40(5), 427.‏ https://dx.doi.org/10.24272/j.issn.2095-8137.2019.027
Zhou, J., Mogollón, J.M., van Bodegom, P.M., Beusen, A.H. and Scherer, L. (2024). Global regionalized characterization factors for phosphorus and nitrogen impacts on freshwater fish biodiversity. Science of the Total Environment, 912, 169108. https://doi.org/10.1016/j.scitotenv.2023.169108