Authors
1 Department of Biology, Faculty of Sciences, Shiraz University, Shiraz, Iran
2 Department of Limnology, Faculty of Sciences, University of Gent, Gent, Belgium
Abstract
Keywords
Main Subjects
مقدمه.
سنجاقک با نام علمی Calopteryx splendens Harris, 1782 پراکنش وسیعی در اوراسیا دارد، اما جمعیتهای مختلف آن، به طور قابل توجهی از لحاظ ریختشناسی و انتخاب زیستگاه یک شکل هستند (Dumont et al., 2005). آنها یک کمپلکس گونهای هولارکتیک (Holarctic) مشتمل بر حدود 20 گونه هستند (Misof et al., 2000). اعضای این گونه را میتوان با توجه به اندازه، رنگ بدن و بال و همچنین رفتارشان و به ندرت از روی مورفولوژی آنها تشخیص داد (Asahina, 1976). تاکنون به طور معمول، اندازه و تراکم رنگدانههای بالها یا لکه بالی برای مطالعه تاکسونومی و دامنه پراکنش گونههای جنس Calopteryx Leach, 1815 استفاده شده است Dumont et al., 1993)؛ Siva-Jothy, 1999؛ Rantala et al., 2000). اما چون صفات متریک و مریستیک معمولاً ممکن است تنها بخشی از اطلاعات برگرفته از موقعیت نسبی لندمارکهایی را ارایه کند که این اندازهگیریها بر مبنای آنها هستند، ارتباطات فضایی بین متغیرهای اندازهگیری شده نادیده گرفته میشوند. مورفومتری هندسی، مقایسه را بر اساس مختصات لندمارکهای هم خاستگاه (homologous) اجرا میکند و در نتیجه برخی برتریها را نسبت با روش مورفومتری سنتی ارایه میدهد. دلایل تنوع در اندازه لکه بالی هنوز مشخص نیست اما عقیده عمومی بر این است که رنگدانههای بال تحت تأثیر عواملی نظیر: دما، انگلهای داخلی Siva-Jothy, 2000)؛ (Koskimäki et al., 2004، محیط زیست (Stock, 1980)، عوامل شیمیایی و حشرهکشها هستند و احتمالاً با سلامت و کیفیت جنس نر ارتباط مستقیم دارد (Cordoba-Aguilar and Cordero-Rivera, 2005). بررسیهای مولکولی اخیر نیز منشأ آنها را نشان داده است Misof et al., 2000)؛ Weekers et al., 2001؛ (Dumont et al., 2005. بنابراین، مورفولوژی یک معیار برجسته برای طبقهبندی این گروه نیست (Dumont et al., 2005).
برای مدتهای طولانی، در مطالعه سیستماتیک حشرات بالدار از ویژگیهای الگوی رگبندی و مورفومتری بال استفاده شده است، گرچه ممکن است این ویژگی برای پی بردن به روابط تبارشناسی به اندازه کافی قابل اعتماد نباشد (Grimaldi, 2001). ترتیب خاص رگبالها در چندین تحقیق برای مطالعه گونههای سنجاقکشکلان استفاده شده است Wootton, 1991)؛ (Zeng, 1996، اما تحقیقی در این زمینه در مورد گونه Calopteryx splendens وجود ندارد. برخلاف مطالعات اکولوژیک و مورفولوژیک متعدد روی گونه C. splendens، مطالعه کمتری در مورد تبارشناسی (phylogeny) و جغرافیای تباری (phylogeography) این گونه و زیرگونههای آن انجام شده است. وجود تعداد زیادی هیبرید و زیرگونه بر اساس تنها یک ویژگی مورفولوژیک که چندان هم قابل اطمینان نیست، یعنی اندازه لکه تیره بال در این گونه، پژوهشگران را به کاربرد شاخصهای زیستی دیگر از جمله نشانگرهای مولکولی (Misof et al., 2000؛ Weekers et al., 2001؛ Dumont et al.,2005) و روشهای نوین مورفولوژیک هدایت کرده است تا بتوانند تنوع جغرافیایی جمعیتها در این گونه را مقایسه کنند. شایان ذکر است که Sadeghi و همکاران (2009) با همین روش، 10 جمعیت اروپایی
C. splendensرا مقایسه کردند.
اخیراً عرصه نسبتاً جدید مورفومتری هندسی که بر اساس مختصات لندمارکهای همولوگ است، فوایدی را در مقایسه با روش مورفومتری سنتی ارایه کرده است. مورفومتری هندسی یک روش کمّی برای مطالعه مقایسهای ریختی و روش آماری مهمی برای آنالیز تنوع اغلب پدیدهها (ساختارهای زیستی) است که توسط Bookstein (1991) تهیه شده و خلاصه آموزندهای از آن توسط Rohlf و Marcus (1993) ارایه شده است. این روش نتایج ارزشمندی را در بسیاری از عرصههای مورفومتری کلاسیک به دست داده است.
انواع مختلفی از فرآیندهای زیستی، از جمله سازش با عوامل جغرافیایی محلی، تکوین اونتوژنیک یا گوناگون شدن تکاملی بلند مدت، ریخت افراد یا بخشهایی از بدن آنها را تحت تأثیر قرار میدهد. ریخت (shape) به عنوان « تمام اطلاعات هندسی که پس از کنارگذاشتن آثار ناشی از اندازه (size)، گردش (translation) و چرخش (rotation) یک جسم باقی میماند » تعریف میشود (Kendall, 1977). اطلاعات مورفومتریک به طور سنتی شامل اندازهگیری فواصل خطی، زوایا و نسبتهای مختلف بین اجزا میشود. این متغیرهای خطی به طور تیپیک تنها بخشی از اطلاعاتی را ارایه میکنند که ممکن است از موقعیت نسبی نقاط مورد اندازهگیری به دست آید (Rohlf, 2002). این اطلاعات شامل اطلاعاتی در مورد ارتباطات فضایی بین متغیرهای اندازهگیری شده نیست و به نظر میرسد که برای تحلیلها و تفسیرهای تکاملی مناسب نباشد (Bogdanowicz et al., 2005).
ریخت بال به ویژه در خانواده Calopterygidae از این جهت جالب است که ویژگیهای آن در حد بالایی حفظ (conserve) شده است. بنابراین، مطالعه ریخت بال در گونههای این خانواده و از جمله گونه
C. splendens میتواند با اطلاعات ارزشمند خود افق جدیدی را در بررسیهای مورفولوژی و تکاملی به ویژه در مورد گونههای پیچیده فراهم کند.
رویکرد اصلی مطالعه حاضر کمّیسازی و تفسیر ریخت بال جلویی چپ در 10 جمعیت آسیایی ازگونه C. splendens با استفاده از روش مورفومتری هندسی بوده است. در این تحقیق، تنوع و گوناگونی لکه بالی به طور کلی نادیده گرفته شد و تنها الگوی بال چپ جلویی در جنس نر ارزیابی شد تا هر گونه اختلاف مربوط به تفاوت بالهای نر و ماده و همچنین اختلاف ناشی از عدم تقارن بالهای چپ و راست حذف شده باشد.
مواد و روشها.
نمونههای مورد مطالعه: در مطالعه حاضر بال جلویی نرها از 10 جمعیت Calopteryx splendens که در غرب و مرکز آسیا زندگی میکنند بررسی شد. تعداد 242 نمونه نر از 10 زیستگاه از کشورهای مختلف شامل: آذربایجان، ازبکستان، ایران (با دو جمعیت از سه زیستگاه)، تاجیکستان، ترکمنستان (هرکدام با یک جمعیت)، ترکیه، روسیه، قرقیزستان و قزاقستان مطالعه شد (جدول 1). مکانهای جمعآوری از نظر جغرافیایی به خوبی از هم جدا بوده است که برای سهولت بیشتر، جمعیت هر محل را با نام کشور مربوط ذکر میکنیم (شکل 1).
جدول 1- تعداد حشرات نر جمعآوری شده گونه Calopteryx splendens برای هر جمعیت و موقعیت محلهای نمونهبرداری
Population (country) |
No. of males |
Localities |
Azerbaijan (Az) |
20 |
Kura valley, Agsu, Azerbaijan, 8.09.2004 (40:34 N, 48:23 E) |
Northern Iran (I6) |
20 |
Golestan, Minudasht, 27.06.2003 ( 37:15 N, 55:28 E) |
Mazandaran, Alamdeh, 22.07.1996 ( 36:33 N, 51:57 E) |
||
Southern Iran (I9) |
26 |
Kohgiluyeh and Boyer Ahmad, Cheshmeh Belgheis, 14.08.2003, (30:43 N, 50:31 E) |
Kazakhstan (Ka) |
20 |
Irtysh River, China border, 11.07.2004 |
Tajikistan (Ta) |
30 |
Gissar, 40 km W of Dushanbe, 22.07.1991 ( 38:30 N, 68:36 E) |
Turkmenistan (Tm) |
27 |
Kara- Kala, W Kopet- Daq, 08.06.1988 (38:29 N, 56:18 E) |
Uzbekistan (Uz) |
26 |
Tahskent, Boshkizilsay River, 4.07.2005 (41:16 N, 69:10 E) |
Kyrgyzstan (Ki) |
20 |
Sary- Tschelek, natural reservoir, 31.07.1985 (39:45 N, 73:10E) |
Russia (Ru) |
31 |
Ismaylowsky Park, Moscow, 12.07.1989 (57:44 N, 37:37 E) |
Turkey (Tu) |
22 |
Golderesi-kemer, Fethiye, Mugla, 21.07.87 (36:39 N, 29:22 E) |
شکل 1- موقعیت جغرافیایی جمعیتهای نمونهبرداری شده از گونه Calopteryx splendens
.جمعآوری اطلاعات و تفسیر ریخت بال: نمونههای بررسی شده به یکی از دو شکل خشک شده یا نگهداری شده در اتانول 70 درصد (preserve) بود که بخشی توسط نگارندگان جمعآوری شده است و بخشی دیگر از مجموعه پروفسور هنری دومنت (دانشگاه گنت، کشور بلژیک) به امانت دریافت و استفاده شد. نمونههای آسیب دیده حذف و بقیه با دقت و با دستگاه اسکنر کتابی (مدل Agfa SNAPSCAN 1236) اسکن و به صورت تصاویر دیجیتال با کیفیت 400 dpi ذخیره شدند. سپس، با استفاده از نرمافزار tps Dig2 (Rohlf, 2010) 19 نقطه شاخص یا لندمارک هم خاستگاه روی تصویر هر بال نشانهگذاری شد (شکل 2).
شکل 2- موقعیت 19 لندمارک استفاده شده روی بال جلویی سمت چپ در Calopteryx splendens.مشخصات لندمارکها:
(1 Costa (C)-Subcosta (Sc)؛ (2 Nodus (N)؛ (3 محل تلاقی Radius2 (R2) و حاشیه بال؛ (4 رأس انتهایی بال، (5 محل تلاقی Media (M) و حاشیه بال؛ (6 محل تلاقی Cubitus1 (Cu1) و حاشیه بال؛ (7 محل تلاقی Cubitus1 (Cu2) و حاشیه بال؛ (8 رأس شکمی anal triangle؛
(9 رأس پروکسیمال anal triangle؛ (10, 11, 12 رأسهای زوایای دیستال arculus؛ (13 محل منشأ رگبال M؛ (14 محل منشأ رگبال Cubitus1؛ (15 زاویه دیستال anal triangle؛ (16 محل منشأ رگبال Cubitus2؛ (17 محل تلاقی Ru+Rs and R+M؛ (18 محل منشأ رگبال IR3؛ (19 محل منشأ رگبال Radius. نامگذاری رگبالهای اقتباس از Dumont (1991).
خطای لندمارکگذاری (digitizing error) با استفاده از پروتکل مربوط به آن (Adriaens, 2007) تعیین شد. خطای تخمین زده شده در بررسی حاضر حدود 1/7 درصد از تنوع مشاهده شده بود که خطای قابل چشمپوشی به شمار میآید. برای تعیین تغییرات مختصات نقاط و روی هم قرار دادن پیکربندی نقاط نشانهگذاری شده و تبدیل آن به یک شکل میانگین مشترک (common mean shape)، تحلیل پروکراست (procrustes) یا GPA استفاده شد (Adams et al., 2004؛ (Zelditch et al., 2004. اختلافها در ریخت بال بین جمعیتها بر حسب آن چه که توسط شبکههای تغییر شکل thin-plate (thin-plate deformation grids) به تصویر کشیده شده بود، توصیف گردید (Rohlf, 2004). سپس مختصات نقاط لندمارکها برای تخمین اندازه بال که centroid size نامیده میشود Bookstein, 1996)؛ (Slice and Bookstein, 2007 با نرمافزار tps Relw (Rohlf, 2007) به دست آمد. نرمالیتی و همگون بودن واریانس برای اندازه بال جمعیتها به ترتیب با آزمونهای Shapiro-Wilk و Leven’s بررسی شد (Milliken and Johnson, 1984). اختلافهای اندازه بال میان جمعیتها و درون جمعیتها با استفاده از آزمون ANOVA یک طرفه سنجیده شد و آزمون Tukey نیز برای مشخص کردن اختلاف زوج جمعیتها نسبت به هم به کار برده شد (Sokal and Rohlf, 1995).
اطلاعات مختصات لندمارکهای به دست آمده از نرمافزار tps Dig، برای آنالیز تجزیه به مؤلفههای اصلی (PCA) و آنالیز کانونی تشخیصی (CVA) در نرمافزارهای PAST (Hammer and Harper, 2007) و IMP (Sheets, 2000) وارد شد تا ترتیب الگوی تنوع رگبالها و اختلاف ریختهای بال در بین جمعیتها آشکار شود. آنالیز تجزیه به مؤلفههای اصلی (PCA) روی ماتریس تشابه واریانس-کوواریانس مختصات نمونهها انجام شد. همچنین، اختلافهای موجود در ریخت بال در بین جمعیتها با اجرای آزمون تحلیل واریانس چند متغیره یا MANOVA (Wilk’s lambda به عنوان معیار و Hotelling’s test به عنوان Post-hoc) مقایسه شد و از آزمون Goodall’s Goodall, 1991)؛ (Bookstein, 1997 و آزمون Hotelling’s T2 به عنوان آلترناتیو، به ترتیب توسط نرمافزارهای IMP و PAST برای آزمودن اختلاف ریخت بال بین گروهها استفاده شد. برای تجسم بهتر تنوعات ریخت بال، میانگینی (consensus) از پیکربندی بال در هر جمعیت با استفاده از نرمافزارهای tpsSmall و tpsRelw (Rohlf, 2003, 2007) محاسبه شد و به عنوان نمونه مرجع برای تعیین اختلاف هر نمونه بال مورد استفاده قرار گرفت. اختلافهای ریخت بال با ایجاد شبکههای تغییر شکل فیزیکی ناشی از thin-plate spline در طول نخستین محور (Relative Warp, RW) با استفاده از نرمافزار tpsSplin (Rohlf, 2004) به تصویر کشیده شد. یک تحلیل خوشهای (cluster analysis) نیز با استفاده از ماتریکس فواصل پروکراست بین جفت جمعیتها انجام شد تا ارتباطات فنتیک بین این 10 جمعیت بررسی شود. فنوگرامی با بیشترین همبستگی ریختگی حاصل از ماتریکس فواصل پروکراست با روش UPGMA.(Unweighted Pair-Group Method with Arithmetic means) در نرمافزار NTsyspc نسخه 1/2 (Rohlf, 2000) رسم شد.
نتایج.
تنوع اندازه بال: آزمون Shapiro-wilk توزیع نرمالی را در هر یک از جمعیتها نشان داد (p>0.05) و آزمون Levene’s یکدستی (homogeneity) واریانس را بر اساس میانگین اندازه آشکار ساخت. نتایج آزمون Tukey HSD روی اندازه بال (centroid size) اختلاف بالایی را بین جمعیتها نشان داد (جدول 2). نتایج نشان داد که جمعیتها از نظر اندازه بال به سه گروه مشخص تقسیم میشوند (شکل 3): نخستین گروه با بزرگترین اندازه بال جمعیت شمال ایران (I6) است، دومین گروه که اندازهای متوسط را نشان میدهند شامل جمعیتهای: قزاقستان، قرقیزستان، روسیه، تاجیکستان، ترکمنستان، ازبکستان و ترکیه و سومین گروه که کوچکترین اندازه بال را داراست جمعیت غرب ایران (I9) و آذربایجان را شامل میشود. آزمونهای T و F برای اندازه بال در جفت جمعیتها نشان داد که جمعیتهای شمالی و غربی ایران (I6 و I9)، و جمعیتهای آذربایجان و قزاقستان با دیگر جمعیتها اختلاف معنیداری در اندازه بال (centroid size) داشتند (جدول 2). گرچه بیشتر جمعیتها اختلاف معنیداری در اندازه بال نشان دادند اما نتیجهگیری خاصی از نظر جغرافیایی قابل استخراج نبود.
جدول 2- نتایج آزمون Tukey HSD روی اندازه بال (مقادیری که معنیدار نیستند با علامت * مشخص شدهاند). Az، آذربایجان؛ I6، شمال ایران؛ I9، غرب ایران؛ Ka، قزاقستان؛ Ki ، قرقیزستان؛ Ru، روسیه؛ Ta، تاجیکستان؛ Tm، ترکمنستان؛ Tu، ترکیه؛ Uz، ازبکستان
|
Az |
I6 |
I9 |
Ka |
Ta |
Tm |
Uz |
Ki |
Ru |
Tu |
Az |
- |
|
|
|
|
|
|
|
|
|
I6 |
108.91 |
- |
|
|
|
|
|
|
|
|
I9 |
29.41 |
138.31 |
- |
|
|
|
|
|
|
|
Ka |
83.83 |
25.08 |
113.23 |
- |
|
|
|
|
|
|
Ta |
37.01 |
71.90 |
66.41 |
46.82 |
- |
|
|
|
|
|
Tm |
61.03 |
47.87 |
90.45 |
22.79 |
24.3 |
- |
|
|
|
|
Uz |
36.00 |
72.90 |
65.41 |
47.82 |
1.00 * |
25.04 |
- |
|
|
|
Ki |
54.11 |
54.79 |
83.52 |
29.72 |
17.11 * |
6.93 * |
18.11 * |
- |
|
|
Ru |
45.75 |
63.16 |
75.15 |
38.08 |
8.74 * |
15.29 * |
9.74 * |
8.36 * |
- |
|
Tu |
25.06 |
83.85 |
54.47 |
58.77 |
11.95 * |
35.98 |
10.94 * |
29.05 |
20.69 |
- |
شکل 3- میانگین اندازه بال (centroid size) در 10 جمعیت آسیایی Calopteryx splendens. Az، آذربایجان؛ I6، شمال ایران؛ I9، غرب ایران؛ Ka، قزاقستان؛ Ki ، قرقیزستان؛ Ru، روسیه؛ Ta، تاجیکستان؛ Tm، ترکمنستان؛ Tu، ترکیه؛ Uz، ازبکستان
تنوع ریخت بال: نتایج حاصل از آنالیز مؤلفههای اصلی (PCA) برای مختصات نقاط که با روش پروکراست استاندارد شدهاند نشان داد که دو مؤلفه اصلی PC1 و PC2 بخش اعظمی از تنوع موجود در ریخت بال را در بر میگیرد (54/49= PC1 و
65/12= PC2) این مقدار بیش از 62 درصد تنوع را برای ماتریس واریانس- کوواریانس آشکار ساخت. تصاویر اغلب نمونهها بر اساس ارتباط PC1 در مقابل PC2 نشان داد که جمعیتهای شمال ایران (I6) و ترکمنستان (Tm) تمایل به گروهبندی در سمت منفی PC1 و PC2 و جمعیتهای تاجیکستان، ازبکستان، قزاقستان و ترکیه تمایل به گروهبندی در سمت مثبت PC1 و PC2 داشتند (شکل 4). نتایج آنالیز واریانس چند متغیره MANOVA Wilk’s lambda) به عنوان معیار و Hotelling’s test به عنوان Post-hoc) نیز به همراه نتایج آزمون Goodall’s در جدول 3 خلاصه شده است. نتایج تحلیل واریانس چند متغیره MANOVA اختلاف معنیداری
(000/0 =P، 00056/0 =Wilk’s lambda) بین جمعیتهای مطالعه شده نشان داد و نتایج مقایسه زوج جمعیتها با آزمون Hotelling مشخص کرد که جمعیتهای قزاقستان و آذربایجان در سطح معنیداری 05/0 و سایر جمعیتها در سطح معنیداری بسیار پایینتر اختلاف ریخت بال را نشان دادند. این نتایج همچنین با آزمون Goodall’s F تأیید شد (جدول 3). تصویر حاصل از نتایج آنالیز متغیرهای کانونی یا CVA بر اساس مؤلفههای اول و دوم نیز به خوبی نشان دهنده تفاوت معنیدار بین ریخت بال در جمعیتهای مورد مطالعه بود (شکل 5).
شکل 4- نتیجه تصویری تجزیه به مؤلفههای اصلی 1 و 2 (PC1 در مقابل PC2)
جدول 3- نتایج آزمون Hotelling’s pair-wise (نیمه بالا) و آزمون Goodall’s F (نیمه پایین) و مقادیر pبرای 10 جمعیت آسیایی Calopteryx splendens. ** <0.01 و * <0.05. Az، آذربایجان؛ I6، شمال ایران؛ I9، غرب ایران؛ Ka، قزاقستان؛ Ki ، قرقیزستان؛ Ru، روسیه؛ Ta، تاجیکستان؛ Tm، ترکمنستان؛ Tu، ترکیه؛ Uz، ازبکستان.
|
Az |
I6 |
I9 |
Ka |
Ta |
Tm |
Uz |
Ki |
Ru |
Tu |
Az |
- |
0.0044** |
0.0003** |
0.05* |
9.5E-9** |
2.2E-5** |
2.8E-7** |
0.005** |
1.4E-10** |
0.0001** |
I6 |
12.76** |
- |
1.99E-5** |
0.001** |
4.1E-8** |
0.0005** |
2.7E-8** |
0.003** |
2.8E-7** |
0.0015** |
I9 |
5.66* |
10.00** |
- |
0.0011** |
2.3E-8** |
9.6E-8** |
2.0E-5** |
8.1E-6** |
9.4E-5** |
0.0008** |
Ka |
14.65** |
14.00** |
9.00* |
- |
2.6E-7** |
1.4E-6** |
0.0007** |
0.01** |
0.001** |
0.0071** |
Ta |
35.3** |
29.95** |
25.41** |
8.69** |
- |
1.8E-12** |
4.5E-6** |
1.8E-6** |
3.1E-12** |
1.17E-7** |
Tm |
8.12** |
2.41 |
6.97** |
14.50** |
37.95** |
- |
2.1E-9** |
2.6E-8** |
4.1E-8** |
1.50E-5** |
Uz |
7.50** |
5.35** |
2.82** |
5.88* |
17.90** |
5.71** |
- |
1.2E-9** |
1.4E-5** |
7.09E-5** |
Ki |
16.64** |
18.80** |
10.90** |
3.24** |
10.14** |
19.32** |
5.26** |
- |
4.99E-5** |
2.87E-5** |
Ru |
16.64** |
8.76** |
4.00** |
8.47** |
29.80** |
7.05** |
2.58* |
8.17** |
- |
2.74E-5** |
Tu |
17.88* |
6.98** |
5.06* |
5.79** |
18.07** |
7.96** |
1.80** |
5.53** |
1.82** |
- |
شکل 5- نتیجه تصویری آنالیز کانونی تشخیصی 1 و 2 (CV1 در مقابل CV2)
تصویر دستهبندی جمعیتها بر اساس آرایش فضایی مؤلفههای اصلی (Relative Warp Ordination plot) نشان داد که یک الگوی جغرافیایی وجود دارد که در تمام دستهبندیهای مؤلفهها میتواند در نظر گرفته شود. برای بال مورد آزمایش، موقعیت نسبی میانگین پیکربندی بال در جمعیتهای ازبکستان، ترکیه، غرب ایران (I9) و روسیه با هم دستهبندی شد و جمعیتهای قزاقستان و قرقیزستان گروه دومی را در بین تاجیکستان و گروه اول تشکیل دادند (شکل 6) از سوی دیگر، جمعیتهای ترکمنستان و شمال ایران (I6) سومین گروه مجزا را تشکیل داد و جمعیت آذربایجان که از دیگر جمعیتها دور افتاده بود، در طرف مقابل گروه اخیر و نسبت به جمعیت ترکیه بسیار نزدیکتر به جمعیت روسیه قرار گرفت.
شبکههای تغییر شکل (thin-plate deformation grids) ریخت بال چپ جلویی برای هر جمعیت، که در مقایسه با برآیند اطلاعات مختصاتی همه جمعیتها به دست آمد، اختلاف ریخت بال مشخصی را بین گروههای جمعیتی جدا شده در شکل 7 به تصویر کشیده است. شبکههای تغییر ریخت بال برای جمعیتهای ترکمنستان و شمال ایران (I6) باریکترین ریخت بال را در بین جمعیتهای مطالعه شده مشخص کرد و برای جمعیت آذربایجان پهنترین ریخت بال را نشان داد. از سوی دیگر، جمعیت تاجیکستان نیز به خاطر وجود تفاوت در قسمت عقبی ریخت بال به وضوح از سایر جمعیتها قابل تشخیص بود (شکل 7).
شکل 6- طرح دستهبندی جمعیتها بر اساس relative warp و اطلاعات میانگین ریخت بال در هر جمعیت
شکل 7- فنوگرام UPGMA و طرح تغییر شبکهای ریخت بال (thin-plate deformation grid) در جمعیتهای کشورهای مختلف آسیای مرکزی. Az، آذربایجان؛ I6، شمال ایران؛ I9، غرب ایران؛ Ka، قزاقستان؛ Ki ، قرقیزستان؛ Ru، روسیه؛ Ta، تاجیکستان؛ Tm، ترکمنستان؛ Tu، ترکیه؛ Uz، ازبکستان
بحث و نتیجهگیری.
سازگاری به عوامل زیستی و غیرزیستی در محیط زیست گونهها اغلب به عنوان تنوع جغرافیایی شناخته میشود McPeek, 1990)؛ Ricklefs and Miles, 1994). ریخت بال و حالت طبیعی آن بسیاری از جنبههای پرواز حشرات از جمله توانایی مانور و ایجاد جریانهای گردابی به وسیله زنش بالها را تحت تأثیر قرار میدهد. چون ریخت طبیعی بالها به تطابق آنها با شرایط محیطشان وابسته است بنابراین بررسی ریخت بال ممکن است به یافتن ویژگیهای کلیدی بالها در ارتباط با انجام پرواز منجر شود و از این رو بایستی تحت تأثیر " انتخاب " قرار گرفته باشد (Sadeghi et al., 2009). در برخی از سنجاقکها و آسیابکها (Outomuro et al., 2013 a,b) و همچنین مگس سرکه (Klingenberg and Gidaszewski, 2010)، ریخت بال یک ساختار فیلوژنتیک را نشان داده است. تحلیلهای ما در مورد تنوع ریخت بال نشان داد که منابع متنوع و احتمالاً مستقل محیطی ریخت بال را در Calopteryx splendensتحت تأثیر قرار داده است. نتایج تحقیق حاضر نشان داد که محیطهای زیست مختلف که به نحوی به صورت مجزا یا ایزوله جمعیتهایی از C. splendens را در خود جای دادهاند عمدتاً اختلاف معنیدار و مشخصی را در ریخت بال جلوییشان نشان میدهند. به بیان دیگر، به استثنای اختلاف بین دو جمعیت شمال ایران (I6) و ترکمنستان (Tm) که معنیدار نبود و انتظار آن نیز وجود داشت (چون این دو جمعیت از نظر جغرافیایی و منشأ به یکدیگر بسیار نزدیک هستند) سایر جمعیتها به طور معنیدار ریخت بال متفاوت نسبت به یکدیگر داشتند. این اختلاف در فنوگرام UPGMA به دست آمده نیز منعکس بود (شکل 7) به طوری که به خوبی سه دسته مشخص و جدا را در فنوگرام نشان داد. نخستین دسته شامل جمعیت تاجیکستان، شرقیترین جمعیت مورد مطالعه است، دومین انشعاب شامل جمعیتهای شمال ایران (I6) و ترکمنستان (Tm) میشود، دو جمعیتی که از نظر جغرافیایی نیز به هم وابستهاند. سومین یا بزرگترین خوشه، شامل دیگر جمعیتهاست. در این انشعاب بزرگ، جمعیت آذربایجان وضعیت کاملاً متمایزی را نسبت به دیگر جمعیتها نشان میدهد. در حالی که جمعیتهای روسیه و ترکیه از نزدیک به هم وابستهاند، جمعیت قزاقستان و قرقیزستان به طور جالبی یک گروه وابسته به هم را تشکیل داده که از جمعیت تاجیکستان جدا شده است. افزون بر این، جمعیت ازبکستان به جمعیتهای ترکیه و روسیه نزدیکتر است تا قزاقستان و قرقیزستان، که در واقع جمعیتهای همسایه آن از نظر جغرافیایی هستند.
اندازه بال (centroid size) به عنوان تنها معیار اندازه حشره، در مطالعه حاضر هیچ گونه شیب پیوسته (cline) جغرافیایی را نشان نداد. اختلاف ریخت بال در بین 10 جمعیت مطالعه شده C. splendens معنیدار بود و بخش اعظم تنوع مشاهده شده (بیش از 60 درصد) در دو مؤلفه اول در آنالیز تجزیه به مؤلفههای اصلی (PCA) خلاصه میشود و این امر یک توضیح روشن هندسی در این جمعیتها را ارایه میدهد. بنابراین، این دو مؤلفه بیشترین مسؤولیت را در بروز تنوع ریخت بال در میان جمعیتهای مطالعه شده دارند. آنالیز واریانس چند متغیره MANOVA اختلاف معنیدار ریخت بال را در بین جمعیتهای مطالعه شده آشکار ساخت، به استثنای دو جمعیت ترکمنستان و شمال ایران (I6) که مقدار F حاصل از آزمون Goodall نیز وجود منشأ یکسانی برای این دو جمعیت را تأیید کرد.
در مطالعه حاضر، فنوگرام UPGMA سه انشعاب جداگانه را به وضوح مشخص کرد. نخستین انشعاب شامل جمعیتهای آذربایجان، روسیه، غرب ایران، ترکیه، ازبکستان، قزاقستان و قرقیزستان بود. دومین شاخه، جمعیتهای شمال ایران و ترکمنستان را در بر گرفت و سومین شاخه تنها شامل جمعیت تاجیکستان بود. جمعیتهای ترکیه و روسیه بیشترین ارتباط را با یکدیگر و با جمعیت ازبکستان نشان دادند. این سه جمعیت در ارتباط نزدیک با جمعیت غرب ایران است. این نتایج همچنین با نتایج پیشین حاصل از تحلیل مولکولی با روش AFLP مطابقت دارد (Sadeghi et al., 2010). جمعیتهای قرقیزستان و قزاقستان نیز در ارتباط نزدیک با هم بوده، به عنوان گروه خواهری گروه قبلی قرار گرفتند. این نکته قابل توجه است که جمعیت آذربایجان در این فنوگرام با فاصله زیاد از دیگر جمعیتها قرار گرفته است. به نظر میرسد که جمعیت آذربایجان نسبتاً از سایر جمعیتهای شاخه اول جدا شده است، این نتیجه پیش از این در نتایج حاصل از AFLP مشاهده شده بود (Sadeghi et al., 2010).
همان طور که Johansson و همکاران (2009) اشاره کردهاند بخش اعظم اختلاف ریخت بال در بین جمعیتهای بررسی شده، نشان دهنده این است که ریخت بال به وسیله ویژگیهای محیط زیست و رفتار تولید مثلی تحت تأثیر قرار میگیرد. البته سنجاقکها در بسیاری از جنبههای حیاتی که میتواند ریخت بال را متأثر سازد از جمله رفتار غذایابی (foraging) و رفتار جفتگیری با یکدیگر متفاوت هستند. اما در بین جمعیتهای یک گونه واحد، برخی از این رفتارهای حیاتی یکسان است و بیشتر اختلافهای قابل مشاهده در ریخت بال در واقع به اختلافهای اکولوژیک ناشی از وضعیت جغرافیایی خاصشان نسبت داده میشود.
نقش برخی رودخانهها از جمله ولگا، اورال، آمودریا، سیر دریا و حتی رودخانه ایرتیش در آسیای مرکزی در پراکنش حشرات آبزی مهم بوده است، چون همه این رودخانهها به حوضه آبی دریای آرال، دریای خزر و دریای سیاه وارد میشوند. بنابراین، پس از جدایی این دریاها از یکدیگر تصور این که جمعیتهای آسیای مرکزی نیز به چهار جمعیت جداگانه، تقسیم شده باشند دور از ذهن نیست. گستردهترین این جمعیتها در بین رودخانه ولگا و دریای خزر و مرز غربی قرقیزستان و تاجیکستان و حتی تا رودخانه ایرتیش یعنی بخش اعظم فلات ایرانی-تورانی مستقر شده است. این منطقه شامل بخش اصلی فلات ایران نیز میشود، کوههای کپه داغ و صحرای قره قوم در جنوب این ناحیه، آن را از دومین ناحیه ایزوله شده که شامل جنوب ترکمنستان و حاشیه جنوبی دریای خزر (در شمال ایران) میشود جدا میکند. علاوه بر کپه داغ، ناحیه ایزوله شده دوم از جنوب به کوههای البرز و بخشی در شمال و شرق محدود به دریای خزر محدود میشود که نسبت به سایر نقاط مطالعه شده دست کم از نظر رطوبت و دما، کاملاً متفاوت است. سومین ناحیه ایزوله شده احتمالاً از مرز غربی تاجیکستان به سمت شرق کشیده شده است. این ناحیه از غرب به دریای سیاه و از شرق به دریای خزر محدود میشود و بین کوههای قفقاز قرار گرفته است. بر اساس اطلاعات پیشین (Sadeghi, et al., 2009, 2010) به نظر میرسد که جمعیت C. splendens در این منطقه، ارتباطی نزدیک و اجدادی با جمعیتهای شمال ترکیه و اروپای مرکزی دارد و این ارتباط با جمعیتهای جنوب ترکیه و ایران بسیار کمتر است.
سپاسگزاری
نگارندگان از زحمات و تلاش دوستان و همکاران سنجاقکشناس که در کشورهای مختلف آسیای مرکزی اقدام به جمعآوری و ارسال نمونهها نمودند و امکان ذکر نام همه آنها در این مقاله ممکن نیست تشکر و قدردانی مینمایند.