مطالعه شباهت های ژنتیکی و تجزیه و تحلیل تبار شناختی شترهای تک کوهانه و دو کوهانه ایرانی با هفت گونه اصلی موجود در خانواده شترسانان بر اساس ژنوم میتوکندریایی

نوع مقاله : مقاله پژوهشی

نویسندگان

1 عضو هیات علمی گروه علوم دامی / دانشگاه بوعلی سینا

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

3 گروه علوم دامی، دانشکده کشاورزی، دانشگاه بوعلی سینا، همدان، ایران

4 واحد ژنتیک و اصلاح نژاد حیوانات (AGBU)، دانشگاه نیوانگلند، آرمیدال، استرالیا

5 سازمان تحقیقات آموزش و ترویج کشاورزی، گیلان، ایران

چکیده

توالی‌های ژنوم‌ کامل میتوکندریایی هفت گونه اصلی شترسانان شامل Vicugna vicugna، Lama guanicoe، Vicugna pacos، Lama glama، Camelus ferus، Camelus dromedarius، Camelus bactrianus، 4 شتر تک کوهانه و 6 شتر دو کوهانه ایرانی همراه با توالی‌های ۱۳ ژن رمزگر (PCGs) به ازای هر ژنوم از سایت NCBI استخراج شدند. همچنین، مجموعه‌ای از 3100 چندشکلی تک نوکلئوتیدی (SNP) از توالی‌های کامل گونه‌های مورد بررسی تهیه شد که پس از فیلتراسیون، برای ساخت ماتریس روابط ژنومی (GRM) مورد استفاده قرار گرفتند. نتایج ترسیم درخت‌های تبارشناختی بر اساس توالی‌های کامل، دو خوشه اصلی همراه با چند زیر خوشه را نشان داد به شکلی که گونه‌های Vicugna vicugna، Lama guanicoe، Vicugna pacos، Lama glama و Camelus ferus یک شاخه اصلی و باقی گونه‌های تک کوهانه (Camelus dromedarius) و دو کوهانه (Camelus bactrianus) همراه با Camelus ferus یک زیر شاخه اصلی دیگر با چند زیر شاخه تشکیل دادند. درصد تشابه ژنتیکی بین 7 گونه اصلی و گونه‌های ایرانی بین 4/82 درصد تا 100 درصد متغیر بود. توالی-های نوکلئوتیدی و آمینواسیدی ژن‌های رمزگر، خوشه‌بندی مشابه با توالی‌های کامل نشان دادند. ژن ATP8 بیشترین و ژن COX1 کمترین تنوع ژنتیکی را داشتند. ترسیم نقشه حرارتی (HeatMap) و نمودار PCA مبتنی بر داده‌های SNP و نمودار سطحی (Level Plot) شباهت‌های نوکلئوتیدی، نتایج مشابهی داشت. یافته‌های این مطالعه نشان می‌دهد که مقایسه توالی‌های ژنوم میتوکندریایی می‌تواند گونه‌های اصلی خانواده شترسانان را به‌صورت قوی‌ طبقه‌بندی کند و روابط تکاملی آن‌ها را روشن سازد. این تحقیق بینش‌های ارزشمندی برای مطالعات آینده در زمینه ژنتیک و تکامل شترسانان فراهم می‌کند.

موضوعات


عنوان مقاله [English]

Study of Genetic Similarities and Phylogenetic Analysis of Iranian Dromedary and Bactrian Camels with Seven Major Species of the Camelidae Family Based on Mitochondrial Genome

نویسندگان [English]

  • Ahmad Ahmadi 1
  • Mostafa Hashemi 2
  • Pouya Zamani 3
  • Mohammad Hosein Ferdosi 4
  • Ramin Abdoli 5
1 Member of Animal Sciences Department / Bu Ali Sina University
2 Department of Animal Science, Faculty of Agriculture, Bu-Ali Sina University, Hamedan, Iran
3 Department of Animal Science, Faculty of Agriculture, Bu-Ali Sina University, Hamedan, Iran
4 AGBU, a joint venture of NSW Department of Primary Industries and University of New England, University of New England, Armidale, NSW 2351, Australia
5 Agricultural Research, Education and Extension Organization (AREEO), Guilan, Iran
چکیده [English]

Introduction: Mitochondrial DNA (mtDNA), also referred to as the mitogenome, has emerged as a pivotal tool in evolutionary biology and phylogenetics due to its maternal inheritance, high mutation rate, and lack of recombination. These characteristics make it an ideal molecular marker for tracing lineage divergences and understanding genetic relationships among species. This study leverages complete mitochondrial genome sequences, protein-coding genes (PCG), and single nucleotide polymorphism (SNP) data to elucidate the genetic similarities and phylogenetic relationships between Iranian dromedary (Camelus dromedarius) and Bactrian (Camelus bactrianus) camels and seven major species of the Camelidae family, including wild and domesticated members (including Vicugna vicugna, Lama guanicoe, Vicugna pacos, Lama glama, Camelus ferus, Camelus dromedarius and Camelus bactrianus). The Iranian dromedary and Bactrian camels play a vital role in desert ecosystems and local livelihoods, yet their populations face challenges due to habitat loss and declining traditional pastoralism The primary objectives of this research were to: 1) compare the complete mitochondrial genomes of Iranian camels with those of other Camelidae species to identify evolutionary divergences and conserved genetic regions, 2) analyze nucleotide and amino acid sequences of 13 protein-coding mitochondrial genes per each mitogenomes to assess functional conservation and divergence, 3) utilize SNP-based genomic data generated based on complete mitogenomes to construct a Genomic Relationship Matrix (GRM) and validate phylogenetic inferences, and 4) resolve the taxonomic classification of Camelidae species and clarify the evolutionary position of Iranian camels within the family.
Materials and Methods: Complete mitochondrial genomes of seven Camelidae species (Vicugna vicugna, Lama guanicoe, Vicugna pacos, Lama glama, Camelus ferus, Camelus dromedarius and Camelus bactrianus) and 10 Iranian camel haplotypes (four Camelus dromedary i.e., Haplotypes 1, 2, 3, and 4, and six Camelus Bactrian i.e., Haplotypes 1, 2, 3, 4, 5, and 6) were retrieved from NCBI. In following nucleotide and amino acid sequences of 13 PCGs (ATP6, ATP8, COX1, COX2, COX3, ND1, ND2, ND3, ND4, ND5, ND6, ND4L and CYTB) were extracted for comparative analysis. Also, a total of 3100 SNPs were derived from whole mitogenome alignments, filtered for quality, and used for population genomic analyses. Maximum Likelihood (ML) method were applied to whole mitogenome sequences and PCG alignments to reconstruct evolutionary relationships. Nucleotide diversity were calculated for PCGs and mitogenome regions. GRM was generated to quantify pairwise genetic distances. Principal Component Analysis (PCA), Heatmap clustering and Level Plot of common nucleotides were performed on the SNP dataset to visualize genetic structuring. A heatmap is a data visualization method based on color coding that is widely used to uncover hidden patterns in genomic data. A Level Plot is a graphical tool for displaying nucleotide similarity or differences between various samples based on SNP data, and is depicted as a color-coded or spectral matrix
Results and discussion: Phylogenetic trees revealed two primary clades: Clade/Claster A: Comprised Vicugna and Lama species (Vicugna vicugna, Lama guanicoe, Vicugna pacos and Lama glama), forming a distinct group. Clade/Claster B: Included Camelus ferus, Camelus dromedarius and Camelus bactrianus with Iranian dromedaries and bactrians forming subclusters closely related to Camelus ferus. Sequence identity between Iranian camels and other species ranged from 82.4% to 100%, with the highest similarity observed among C. dromedarius haplotypes. Analysis of nucleotide and amino acid sequences of all protein-coding genes (PCGs) showed a clustering pattern similar to that of the complete mitochondrial genomes. However, ATP8 exhibited the highest variability, while COX1 was the most conserved, reflecting differential evolutionary pressures. Heatmap, PCA, and Level Plot analyses of SNP-based nucleotide similarities and GRM data confirmed the results from the complete mtDNA analysis (the separation of Vicugna/Lama and Camelus camelids, with Iranian camels grouping distinctly within the latter), demonstrating that SNPs are an effective tool for phylogenetic and evolutionary studies. SNP-based analyses proved highly congruent with traditional mtDNA methods, validating their utility for high-resolution phylogenetic studies.
Conclusion: This study provides a comprehensive genomic framework for understanding the evolutionary history of Camelidae, demonstrating the power of mitochondrial genomics and SNP-based approaches in phylogenetics as comparative genomic analysis of mtDNA can robustly classify the major species of the Camelidae family and clarify their evolutionary relationships. The results of the present study not only refine the classification of Iranian camels but also offer a foundation for future research on camelid adaptation, domestication, and conservation.