The pattern of intrapopulational and interpopulational changes of Betula pendula in Iran, based on leaf morphological traits

Authors

1 Department of Forestry, Faculty of Natural Resources, University of Guilan, Rasht, Iran

2 Department of Agronomy and Plant Breeding, Faculty of Agriculture, University of Guilan, Rasht, Iran

Abstract

Despite superior medicinal properties of Birch, there aren’t any studies reported on its morphological diversity. For the first time, leaf samples were collected from four different sites, each site 10 trees and each tree 30 fertile and sterile shoots and leaves (totally 2400). In analyzing the mean values, the Siahmarz kuh showed the maximum and the Marmishoo minimum size mean. Angle of leaf base and number of veins had minimum and fluctuating asymmetry and distance border had a maximum coefficient of variation, minimum and maximum morphological diversity respectively. Average size of sterile leaf traits was larger than average size of fertile leaf traits. It seemed that Siahmarz Kuh is had the most interapopulation diversity and could apply a transition source to sink. The slightest change in the veins of the population in birch was more strongly controlled by genes. Multivariate analyses of variance indicated high statistical significance of all factors: leaf type, population, type population, tree nested population and type × tree nested population. The cluster analysis identified three clusters, the first cluster contained only Marmishoo trees. Because the majority of the population, geographic distribution of genotypes was not consistent with the grouping that could be due to the influence of environmental factors on morphological traits, so it is recommended to use molecular markers. Canonical discriminant analysis showed that petiole length ratio and fluctuating asymmetry could be used in future research on morphological diversity among populations of B. pendula.

Keywords

Main Subjects


مقدمه

گونه غان (Betula pendula) دارای خواص ممتاز دارویی (ضد سرطان قوی و پر کاربرد در درمان ایدز و هپاتیت) است (Zyryanova, 2010). این گونه می‌تواند در حفاظت از خاک و مقابله با بهمن در ارتفاعات استفاده شود و قابلیت کاشت در مناطق سردسیر کشور را نیز دارد. علیرغم نادر بودن رویشگاه‌ها، بومی بودن، در معرض تهدید بودن، ارزش بالای چوب در صنایع مختلف و حفظ سیمای طبیعی کشور توسط این گونه ارزشمند جنگلی، اطلاعاتی در مورد تنوع مورفولوژی درون جمعیتی و بین جمعیتی این گونه وجود ندارد. این امر ممکن است به علت پراکنش محدود آن در ایران و گسترش این گونه در مناطق صعب‌العبور، با آب و هوای نامساعد باشد. به طوری که به علت عدم شناخت محل دقیق رویشگاه، نبودن جاده و دوری راه نمونه‌گیری در این پژوهش با دشواری همراه بود. از اهداف زیربنایی سازمان جنگل‌ها، مراتع و آبخیزداری، دستیابی به بهترین روش توسعه جنگل در ارتفاعات بالای جنگل‌های شمال متصل به مراتع کوهستانی است که یقیناً موفقیت این برنامه‌های دراز مدت مستلزم شناخت کافی از ساختار تنوع ژرم‌پلاسم موجود است. به طور کلی، از نظر مدیریتی حفظ توده جنگلی به صورت پایدار، نیازمند دریافت اطلاعات مستمر در زمینه سطح تنوع ژنتیکی درون و بین جمعیت‌ها است (Thormann et al., 1994). تعیین جمعیت برتر در توسعه جنگل با عنایت به انتقال source-sink (انتقال از جمعیت‌های اهدا کننده به جمعیت‌های در معرض خطر که نوعی برنامه حفاظتی است به طوری که حتی یک فرد می‌تواند با وارد کردن ژن‌های جدید به جمعیت دیگر باعث نجات ژنتیکی و افزایش نرخ مهاجرت شود) و ارزیابی استعداد آنها برای استفاده در کاربری‌های تخصصی مانند قرق، حفاظت و بازسازی از دیگر ضروریات جدید در تحقیقات جنگلداری است چرا که مدیریت مؤثر منابع ژنتیکی مستلزم داشتن اطلاعات جامع در خصوص اندازه و ساختار جمعیت، پراکنش جغرافیایی و محیط رویشگاه و تنوع بین و درون آنها است.

از نظر مدیریتی نیز حفظ توده جنگلی به صورت پایدار، نیازمند اطلاعاتی در زمینه سطح تنوع بین و درون جمعیت‌ها است Thormann et al., 1994)؛ Salehi Shanjani et al., 2002). بررسی‌ها نشان داده است که مورفولوژی برگ می‌تواند دید مناسبی از روابط تکاملی و فیلوژنی ارائه دهد که این امر گامی مؤثر در حفاظت از اکوسیستم‌های طبیعی است (Thorssn et al., 2007). بررسی‌ها در موضوع مورفولوژی برگ فراوان و وسیع است و پژوهشگران از روش‌ها و شاخص‌های آماری متفاوتی استفاده کرده‌اند Wilsey et al., 1998)؛ Kovacic and Simic, 2001؛ Franiel and Wieski, 2005).

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

 

مواد و روش‌ها

مواد گیاهی و نمونه‌گیری

در اواخر تابستان 1390 از چهار رویشگاه اصلی غان در استان‌های آذربایجان غربی، البرز، گلستان و مازندران از هر رویشگاه 10 درخت با فاصله تقریباً 100 متر از یکدیگر به طور تصادفی انتخاب و از هر درخت 30 برگ از شاخه‌ زایا و 30 برگ از شاخه نازا از تمامی جهات تاج و در مجموع 2400 برگ سالم جمع‌آوری شد (جدول 1).

درخت غان دارای شاخه‌های زایا و نازا است: در شاخه زایا، برگ‌ها به صورت متقابل با یکدیگر و شاتون ماده در انتهای آن قرار دارد و در شاخه نازا دو برگ به صورت مخالف یکدیگر قرار گرفته‌اند، بنابراین در این تحقیق برگ‌ها از شاخه‌های زایا و نازا جداگانه جمع‌آوری شد. پس از طی دوره خشک شدن برگ، تعدادی صفت مطابق روش Kovacic و Simic (2001) مربوط به کشور کرواسی انتخاب و به صورت دقیق اندازه‌گیری شد. برای تشخیص این که کدام صفت مورفولوژیک برگ غان در شرایط رویشگاهی ایران بهتر از همه قابلیت ایجاد تمایز و تشخیص دارد، تلاش شد چند صفت دیگر به آن اضافه گردد. صفات انتخاب شده عبارت بودند از: طول دمبرگ (ptl به میلی‌متر)، طول پهنک (bl به میلی‌متر) و عرض پهنک (bw به میلی‌متر) که از ویژگی‌های مربوط به اندازه برگ هستند. فاصله از قاعده برگ تا عریض‌ترین نقطه پهنک (dlb به میلی‌متر)، زاویه قاعده (alb به درجه)، اختلاف دو لبه برگ (db به میلی‌متر) و عدم تقارن پهنک (fa) از ویژگی‌های مربوط به شکل برگ هستند. تعداد رگبرگ در قسمت راست پهنک (nvr)، تعداد دندانه بین رگبرگ دوم و سوم طرف راست پهنک (ntr)، عدم تقارن پهنک (fa) و طول دمبرگ نسبی (ptlr).

fa = |wl - wr| / (wr + wl) .

wl= اندازه نیمه چپ برگ

wr= اندازه نیمه راست

ptlr = ptl/bl (شکل 1).

 

 

جدول 1- ویژگی‌های مناطق بررسی شده

جمعیت

استان

ارتفاع از سطح دریا (متر)

طول جغرافیایی

عرض جغرافیایی

سیاه‌مرزکوه

گلستان

2344

′01 °55

′38 °36

سنگنده

مازندران

2579

′13 °53

′59 °35

شهرستانک

البرز

2404

′23 °51

′44 °35

مارمیشو

آذربایجان غربی ( مرز ایران و ترکیه)

1741

′35 °44

′34 °37

 

شکل 1- صفات اندازه‌گیری شده بر روی برگ غان: طول دمبرگ (ptl)، طول پهنک (bl)، عرض پهنک (bw)، فاصله از قاعده برگ تا عریض ترین نقطه پهنک (dlb)، زاویه قاعده (alb)، اختلاف دو لبه برگ (db)، تعداد رگبرگ در قسمت راست پهنک (nvr) و تعداد دندانه بین رگبرگ دوم و سوم طرف راست پهنک (ntr).

 

 

تجزیه آماری

برای دو نوع تیپ برگ‌ شاخه‌های زایا و نازا در جمعیت‌های مختلف، میانگین X، انحراف معیارSD و ضریب تغییرات CV برای هر صفت محاسبه شد. از نرم‌افزار SAS برای تحلیل واریانس یک متغیره و چند متغیره استفاده شد. شاخص‌های بررسی شده برای آنالیز واریانس عبارت بودند از: تیپ برگ، جمعیت، اثر متقابل تیپ و جمعیت، درخت آشیانه شده در جمعیت و اثر متقابل تیپ و درخت آشیانه شده در جمعیت. در تجزیه خوشه‌ای با روش Ward (1963) برای شناسایی گونه‌های مورد مطالعه و گروه‌بندی آنها از نظر تمام صفات اندازه‌گیری شده از 40 درخت مورد مطالعه استفاده شد. تجزیه تشخیص کانونیک، ترکیباتِ خطیِ صفاتی که بیشترین تفاوت را بین گروه‌ها ایجاد می‌کند را معرفی می‌کند (Dillon and Goldstein, 1984). از این تجزیه، در شناسایی صفاتی که در ایجاد تمایز مؤثر هستند، به کار گرفته شد. برای تعیین همبستگی بین صفات مورفولوژیک برگ غان با عامل ارتفاع از سطح دریا از ضریب همبستگی اسپیرمن استفاده شد (Sokal and Rohlf, 1989).

 

نتایج

به طور کلی، منطقه سیاه‌مرزکوه در استان گلستان بیشترین و مارمیشو در مرز ایران و ترکیه کمترین خصوصیات مربوط به اندازه برگ bw)، ptl و (bl را دارند. خصوصیات مربوط به شکل برگ alb) و (dlb اغلب رابطه معکوس با یکدیگر دارند. عمدتاً برای همه صفات میانگین صفت در برگ‌ شاخه‌های نازا از زایا بزرگتر بود. صفت db برای برگ‌‌ شاخه‌های نازا و زایا مساوی بود، صفت alb بیشترین انحراف معیار را داشت. تعداد رگبرگ کمترین تغییر را در بین جمعیت‌ها نشان داد (جدول 2). برای برگ‌ شاخه‌های زایا کمترین CV مربوط به nvr جمعیت‌های سیاه‌مرزکوه و شهرستانک و alb در سنگده و مارمیشو مشاهده شد و بیشترین CV مربوط به db در جمعیت‌های سیاه‌مرزکوه و شهرستانک و fa در سنگده و مارمیشو مشاهده شد. برای برگ‌ شاخه‌های نازا کمترین CV برای صفت nvr در جمعیت‌های سیاه‌‌مرزکوه، مارمیشو و شهرستانک و alb در سنگده مشاهده شد و بیشترین CV برای صفت fa در جمعیت‌های سیاه‌مرزکوه، سنگده و شهرستانک و db در مارمیشو بود (جدول 3).

آنالیز واریانس چند متغیره MANOVA اختلاف معنی‌داری را برای همه شاخص‌های تیپ برگ، جمعیت، اثر متقابل تیپ برگ و جمعیت، درخت آشیانه شده در جمعیت و اثر متقابل تیپ و درخت آشیانه شده در جمعیت نشان داد. نتایج ANOVA نشان داد که در تمامی صفات، اثر شاخص‌های درخت آشیانه شده در جمعیت و اثر متقابل تیپ و درخت آشیانه شده در جمعیت معنی‌دار بود. اختلاف جمعیتی بین تیپ‌های برگ اختلاف معنی‌داری را برای ptl، ptlr، dlb و db نشان نداد اما برای جمعیت‌ها در تمامی صفات اختلاف معنی‌داری وجود داشت. رابطه بین تیپ برگ و جمعیت به غیر از صفات ntr، fa و db اختلاف معنی‌داری برای صفات نداشت (جدول 4).


 

 

جدول 2- میانگین X، انحراف معیار SD و ضریب تغییرات CV صفات مورفولوژی برگ غان. طول دمبرگ (ptl)، طول پهنک (bl)، عرض پهنک (bw)، فاصله از قاعده برگ تا عریض ترین نقطه پهنک (dlb)، زاویه قاعده (alb)، اختلاف دو لبه برگ (db)، تعداد رگبرگ در قسمت راست پهنک (nvr) و تعداد دندانه بین رگبرگ دوم و سوم طرف راست پهنک (ntr)، عدم تقارن پهنک (fa) و طول دمبرگ نسبی (ptlr).

صفت

شاخص‌های جمعیت

برگ‌ شاخه‌های نازا

برگ‌ شاخه‌های زایا

سیاه‌‌مرزکوه

سنگنده

شهرستانک

مارمیشو

سیاه‌‌مرزکوه

سنگنده

شهرستانک

مارمیشو

Ptl

X (mm)

8/22

1/19

8/21

0/18

3/20

4/17

0/20

9/16

SD (mm)

9/4

6/3

06/5

3/4

9/3

39/3

1/5

3/4

CV (%)

7/21

4/20

1/23

0/23

1/19

4/19

4/25

5/24

Bl

X (mm)

3/62

4/56

5/59

8/55

9/55

7/53

9/55

4/52

SD (mm)

2/11

0/8

06/9

76/7

16/8

24/9

9/9

05/8

CV (%)

9/17

8/14

2/15

19/13

5/14

8/17

5/17

0/15

Ptlr

X (mm)

9/36

5/33

8/36

9/31

6/36

1/34

5/35

2/33

SD (mm)

5/6

3/6

7/6

2/6

0/6

4/6

51/8

3/6

CV (%)

6/17

9/18

2/18

5/19

3/16

7/18

9/23

1/19

Bw

X (mm)

7/51

8/46

5/47

4/48

6/43

5/43

8/41

2/37

SD (mm)

2/9

8/6

3/7

9/7

1/8

8/8

19/7

09/5

CV (%)

8/17

6/14

48/15

3/16

5/18

3/20

1/17

6/13

Dlb

X (mm)

2/18

2/16

71/17

14/16

43/17

8/16

41/17

0/17

SD (mm)

4/4

0/4

29/4

1/3

51/4

04/4

57/4

5/5

CV (%)

2/42

8/24

2/42

2/19

8/25

9/23

8/24

1/31

Alb

X (mm)

7/168

3/167

159

3/152

1/146

9/146

4/137

2/117

SD (mm)

8/29

2/19

9/27

3/24

48/27

2/16

1/19

0/12

CV (%)

6/17

4/11

6/17

16

78/18

03/11

9/13

2/10

Nvr

X

3/8

2/8

1/8

4/7

9/7

9/7

8/7

8/6

SD

1/1

9/0

0/1

95/0

0/1

06/1

0/1

0/1

CV (%)

1/13

8/11

4/12

8/12

1/13

3/13

0/13

8/14

Ntr

X

7/4

2/4

0/5

1/3

68/3

4/3

64/3

3/2

SD

6/1

37/1

5/1

8/0

18/1

2/1

06/1

49/0

CV (%)

9/33

07/32

08/30

28

2/34

4/31

06/28

4/21

Fa

X (mm)

79/2

0/3

01/3

03/0

1/3

7/2

1/3

04/0

SD (mm)

2/1

4/1

4/1

02/0

38/1

3/1

5/1

02/0

CV (%)

47/45

48/47

4/46

58/45

17/46

96/48

7/41

3/42

Db

X (mm)

5/0

5/0

2/0

5/0

57/0

49/0

1/0

48/0

SD (mm)

2/0

2/0

09/0

2/0

28/0

22/0

06/0

18/0

CV (%)

7/44

9/43

7/44

9/45

9/48

6/46

9/45

0/37


جدول 3- میانگین X، انحراف معیار SD و ضریب تغییرات CV صفات برای 1200 برگ شاخه‌های نازا و 1200 برگ شاخه‌های زایا. طول دمبرگ (ptl)، طول پهنک (bl)، عرض پهنک (bw)، فاصله از قاعده برگ تا عریض ترین نقطه پهنک (dlb)، زاویه قاعده (alb)، اختلاف دو لبه برگ (db)، تعداد رگبرگ در قسمت راست پهنک (nvr) و تعداد دندانه بین رگبرگ دوم و سوم طرف راست پهنک (ntr)، عدم تقارن پهنک (fa) و طول دمبرگ نسبی (ptlr).

صفت

شاخص‌های جمعیت

برگ‌ شاخه‌های نازا

برگ‌ شاخه‌های زایا

Ptl

X (mm)

4/20

9/18

SD (mm)

5/4

2/4

CV (%)

1/22

2/22

Bl

X (mm)

7/58

5/54

SD (mm)

0/9

9/8

CV (%)

3/15

3/16

Ptlr

X (mm)

8/34

9/34

SD (mm)

5/6

8/6

CV (%)

6/18

5/19

Bw

X (mm)

6/48

5/41

SD (mm)

8/7

3/7

CV (%)

1/16

5/17

Dlb

X (mm)

1/17

6/17

SD (mm)

0/4

7/4

CV (%)

1/23

5/26

Alb

X (mm)

9/161

0/137

SD (mm)

4/25

7/18

CV (%)

7/15

5/13

Nvr

X

0/8

6/7

SD

0/1

0/1

CV (%)

6/12

6/13

Ntr

X

3/4

3/3

SD

3/1

0/1

CV (%)

1/31

1/29

Fa

X (mm)

2/2

4/2

SD (mm)

0/1

1/1

CV (%)

2/46

8/44

Db

X (mm)

4/0

40/0

SD (mm)

2/0

2/0

CV (%)

8/44

2/44

 


جدول 4- نتایج MANOVA و ANOVA برای تمام نمونه‌ها در پنج شاخص‌ مختلف. *، ** و *** به ترتیب نشان دهنده سطح معنی‌داری در سطوح 5، 1 و 1 هزارم درصد و ns عدم معنی‌داری است. طول دمبرگ (ptl)، طول پهنک (bl)، عرض پهنک (bw)، فاصله از قاعده برگ تا عریض ترین نقطه پهنک (dlb)، زاویه قاعده (alb)، اختلاف دو لبه برگ (db)، تعداد رگبرگ در قسمت راست پهنک (nvr) و تعداد دندانه بین رگبرگ دوم و سوم طرف راست پهنک (ntr)، عدم تقارن پهنک (fa) و طول دمبرگ نسبی (ptlr).

تجزیه واریانس

 

تیپ برگ

جمعیت

جمعیت * تیپ

درخت (جمعیت)

تیپ * درخت (جمعیت)

MANOVA

درجه آزادی صورت

10

30

30

360

360

درجه آزادی مخرج

27

87

87

22017

22017

سطح معنی‌داری

7/6***

1/62***

2/3**

3/23***

7/7***

ANOVA

درجه آزادی صورت

1

3

3

36

36

درجه آزادی مخرج

36

36

36

2160

2160

ptl

6/2 ns

8/36***

5/1 ns

5/45***

4/5***

bl

9/5*

7/6**

8/0 ns

8/40***

5/14***

ptlr

01/0ns

4/17***

4/1ns

8/22***

1/4***

bw

7/25***

3/5***

9/3*

1/33***

7/12***

dlb

8/0ns

0/4**

9/1ns

5/34***

4/6***

alb

3/26 ***

0/14***

5/1ns

4/57***

1/19***

nvr

5/4*

5/20***

47/0ns

8/36***

1/11***

ntr

1/19***

8/31***

3*

9/31***

6/11***

fa

6/14***

7/429***

6/6***

7/1***

4/2***

db

75/0ns

3/260***

6/7***

3/14***

0/2***

 

 

در بررسی حاضر، آنالیز خوشه‌ای با روش Ward (1963) 40 درخت مطالعه شده را در سه خوشه مختلف قرار داد: خوشه اول شامل 10 گونه، خوشه دوم شامل 9 گونه و خوشه سوم شامل 21 گونه بود. خوشه اول فقط شامل گونه‌های جمعیت مارمیشو، خوشه دوم شامل گونه‌های سیاه‌مرزکوه و شهرستانک و خوشه سوم شامل گونه‌های جمعیت سنگنده و تعدادی هم از سیاه‌مرزکوه و شهرستانک بود (شکل 2).

نتایج حاصل از تجزیه تشخیص کانونیک نشان داد که fa و ptlr هم در برگ شاخه‌های زایا و هم در شاخه‌های نازا در شناسایی صفات مؤثر در ایجاد تمایز بین جمعیتی با موفقیت بیشتری همراه بوده‌اند (جدول 5).

نتایج همبستگی اسپیرمن بین عامل ارتفاع از سطح دریا و متغیرهای مورفولوژیک برگ نشان داد که در مجموع، همبستگی ضعیفی بین ارتفاع از سطح دریا با صفات مهم مورفولوژیک برگ وجود دارد. بیشترین مقدار همبستگی بین ارتفاع از سطح دریا با صفت تعداد دندانه بین رگبرگ دوم و سوم طرف راست پهنک (ntr) مشاهده شد (جدول 6).


 

 

 

شکل 2- آنالیز خوشه‌ای با روش وارد بر مبنای صفات مورفولوژی برگ غان (Si1-Si10 سیاه‌مرزکوه، Sa1-Sa10 سنگنده، Sh1-Sh10 شهرستانک و Ma1-Ma5 مارمیشو)

 

 

جدول 5- ضریب استاندارد کانونیک صفات اندازه‌گیری شده در برگ غان به همراه مقادیر ویژه. طول دمبرگ (ptl)، طول پهنک (bl)، عرض پهنک (bw)، فاصله از قاعده برگ تا عریض ترین نقطه پهنک (dlb)، زاویه قاعده (alb)، اختلاف دو لبه برگ (db)، تعداد رگبرگ در قسمت راست پهنک (nvr) و تعداد دندانه بین رگبرگ دوم و سوم طرف راست پهنک (ntr)، عدم تقارن پهنک (fa) و طول دمبرگ نسبی (ptlr).

متغیر کانونیکی اول

 

برگ‌ شاخه‌های زایا

برگ‌ شاخه‌های نازا

صفت

661/0 -

285/0

ptl

314/0

154/0 -

lbl

721/0

381/0 -

ptlr

142/0

182/0

lbw

101/0

117/0

dlb

196/0

074/0

alb

241/0

254/0

nvr

292/0

339/0

ntr

892/0

782/0

fa

385/0 -

374/0 -

db

5/72

1/73

واریانس

841/0

798/0

همبستگی کانونیک

4/2

75/1

مقدار ویژه

 


جدول 6- همبستگی اسپیرمن بین متغیرهای مورفولوژی برگ و عامل ارتفاع از سطح دریا. طول دمبرگ (ptl)، طول پهنک (bl)، عرض پهنک (bw)، فاصله از قاعده برگ تا عریض ترین نقطه پهنک (dlb)، زاویه قاعده (alb)، اختلاف دو لبه برگ (db)، تعداد رگبرگ در قسمت راست پهنک (nvr) و تعداد دندانه بین رگبرگ دوم و سوم طرف راست پهنک (ntr)، عدم تقارن پهنک (fa) و طول دمبرگ نسبی (ptlr).

صفت

ptl

lbl

ptlr

lbw

dlb

alb

nvr

ntr

fa

db

همبستگی

40/0-

40/0-

2/0

4/0

2/0-

4/0

4/0

8/0

4/0

6/0-

سطح احتمال

6/0

6/0

8/0

6/0

8/0

6/0

6/0

2/0

6/0

4/0

 


بحث

اندازه اغلب صفات اندازه‌گیری شده در این مطالعه از سیاه‌مرزکوه (شرق ایران) به طرف مارمیشو (غرب ایران) در حال کاهش است. از آنجا که CV واحد ویژه‌ای ندارد، برای اندازه‌گیری تنوع صفات، معیار مناسبی است به طوری که CV بیشتر نشان دهنده تنوع فنوتیپی بیشتر است. صفاتی که ضریب تغییرات بالایی دارند، محدوده وسیع‌تری از کمّیت را دارا هستند که دامنه انتخاب بیش‌تری برای آن صفات محسوب می‌شود. در بررسی حاضر، alb و nvr کمترین و fa و db دارای بیشترین CV هستند و به ترتیب کمترین و بیشترین تنوع فنوتیپی را دارند. از آنجا که در هر جمعیت، افراد تشکیل دهنده یکسان بوده، صفاتی که تنوع بیشتری دارند نشان دهنده اثر بیشتر واریانس محیطی در آنها است، یعنی قابلیت توارث‌پذیری آن کمتر است و بیشتر تحت تأثیر محیط قرار می‌گیرند. بیشترین CV صفات مختلف (بیشترین تنوع درون جمعیتی) در سیاه‌مرزکوه وجود داشت و به نظر می‌رسد بهترین جمعیت از لحاظ تنوع درون جمعیتی است و می‌تواند از آن در انتقال source-sink استفاده نمود.

از میان تمام صفات اندازه‌گیری شده در این بررسی، تعداد رگبرگ کمترین تغییر را داشت. برای درخت بارانک (Sorbus aucuparia) نیز این صفت کمترین تغییر را داشته است (Bednorz and Byzia, 2005). چون این صفت برای زنده‌مانی گونه با اهمیت است و در شرایط محیطی گوناگون واریانس کمتری دارد، بنابراین می‌تواند تحت کنترل ژن‌های کیفی و دارای توارث‌پذیری بالایی باشد. صفات اندازه‌گیری شده از روی برگ‌ شاخه‌های زایا عمدتاً کوچکتر از برگ‌ شاخه‌های نازا هستند که این امر ممکن است ناشی از نیاز شاخه زایا برای مصرف انرژی در پر کردن دانه، تولید بذر و شاتون باشد، همان طور که در مشاهدات میدانی مشاهده شد که گونه غان بذرهای فراوانی بر روی شاخه زایا تولید کرده است.

نتایج MANOVA در این پژوهش نشان داد که جمعیت‌های غان تغییرات درون جمعیتی بیشتری از تغییرات بین جمعیتی دارند. این قاعده توسط Kovacic و Simic (2001) نیز ثابت شده بود. معنی‌داری رابطه بین تیپ برگ و درخت آشیانه شده در جمعیت‌ها این موضوع را نشان می‌دهد که رابطه بین برگ‌ شاخه‌های زایا و نازا در بین گونه‌های درون هر جمعیت متفاوت است. گونه‌ها در جنگل طبیعی، تنوع شایان توجهی در خصوصیات مورفولوژیک بروز می‌دهند که نوعی فرآیند سازگاری محلی است (Boratynska, 2002). سازمان جنگل‌ها، مراتع و آبخیزداری می‌تواند گونه غان را که پتانسیل زیادی برای سازگاری محلی در شرایط محیطی مختلف دارد را برای جنگل‌کاری در مناطق دیگر برای تولید چوب و ایجاد پوشش گیاهی به کار گیرد.

دندروگرام نشان داد که تنها گونه‌های منطقه مارمیشو در یک خوشه جداگانه قرار گرفته‌اند، یعنی این جمعیت به دلیل متفاوت بودن با سایر جمعیت‌ها باید بیشتر مورد حفاظت قرار گیرد و سعی شود بذر مناطق دیگر به این منطقه آورده شود تا بر تنوع ژنتیکی آنها افزوده شود. در کلاستر مارمیشو نوعی همخوانی بین الگوی جغرافیایی و ژنوتیپ‌ها مشاهده شد که به علت فاصله جغرافیایی زیاد این منطقه با سایر رویشگاه‌ها است که احتمالاً سبب تغییرات ژنتیکی قابل توجهی در این جمعیت شده است. جمعیت کوهستانی مارمیشو در مرز ایران و ترکیه و خارج از شمال کشور قرار دارد. رویشگاه غان در آن منطقه کوهسری (orophyte) محسوب می‌شود. فارغ از ارتفاع از سطح دریا، شرایط کوهسری در تمامی رویشگاه‌های غان در کشور تا حدودی مشترک است. این شرایط بوم‌شناختی عبارتند از: درجه حرارت محدود کننده، سرعت باد، ارتفاع کل درخت به صورت محدود، قطر درخت به صورت محدود و در برخی موارد رویش قطری در خارج از مرکز و غیره. شاید علت شباهت صفات مورفولوژیکی جمعیت‌های دیگر که با هم در یک خوشه هستند، به خاطر تنظیم ژنتیکی آنها باشد. از آنجا که در اغلب جمعیت‌ها، توزیع جغرافیایی با گروه‌بندی ژنوتیپ‌ها همخوانی ندارد و این امر می‌تواند به دلیل تأثیر عوامل بوم‌شناختی بر صفات مورفولوژیک باشد، بنابراین پیشنهاد می‌شود از نشانگرهای مولکولی هم استفاده گردد. گونه غان در ایران در مناطقی با تنوع جغرافیایی زیاد و به صورت گسسته وجود دارد و عموماً همین تنوع جغرافیایی عامل مؤثر در ایجاد تنوع است (Dunlap and Stettler, 1995). سازگاری محلی می‌تواند در میزان تمایز بین جمعیتی تأثیرگذار باشد، بدین معنی که کاهش جریان ژنی در گونه‌هایی از غان که سازگاری محلی کسب کرده باشند، رخ می‌دهد (Williams and Arnold, 2001). در نهایت، همین نوع سازگاری به تکامل در مولکول‌های زیستی منجر می‌شود (Gillespie, 1991). پژوهش‌های گذشته نشان می‌دهد که عامل توپوگرافی در رویشگاه‌های غان در مناطق مرکزی ژاپن در ایجاد "تمایز ویژه محل" تأثیرگذار بوده است (Tsuda et al., 2010) بنابراین پیشنهاد می‌شود تأثیر توپوگرافی بر تنوع همین صفات بررسی شود.

از داده‌های به دست آمده از تجزیه تشخیص کانونیک می‌توان برای گروه‌بندی درختان به زیر گروه‌های کوچکتر که شباهت زیادی درون آنها وجود دارد، استفاده نمود. بر مبنای ضریب استاندارد کانونیک به دست آمده در این بررسی، پیشنهاد می‌شود که در بررسی‌های مربوط به تنوع ژنتیکی میان جمعیت‌های غان در آینده از دو صفت مهم عدم تقارن پهنک (fa) و طول دمبرگ نسبی (ptlr) نیز استفاده شود.

شایان ذکر است که صفات مورفولوژیک برگ غان نظیر شکل برگ، از عامل محیطی ارتفاع از سطح دریا تأثیر چندانی نمی‌پذیرند (Kovacic and Nikolik, 2005) که با نتایج پژوهش حاضر هماهنگی دارد. گونه‌هایی مانند راش (Fagus orientalis) که در محدوده وسیع‌تری از عامل ارتفاع از سطح دریا رویش دارند (600 تا بیش از 2000 متری ارتفاع از سطح دریا در ایران) در رابطه با عامل ارتفاع از سطح دریا تنوع زیستی بیشتری نشان می‌دهند (Salehi Shanjani et al., 2002).


سپاسگزاری


نگارندگان از اداره کل منابع طبیعی گلستان، شرکت سهامی بهره‌برداری چوب فریم مازندران و دانشگاه ارومیه که نهایت همکاری را برای در اختیار گذاشتن نمونه‌های برگ از رویشگاه‌های صعب‌العبور غان داشته‌اند، سپاسگزاری می‌نمایند.

 

 

منابع
Bednorz, L. and Byzia, A. (2005) Morphological leaf variability of Sorbus aucuparia in the Wielkopolska national park. Roczniki Akademii Rolniczej w Poznaniu 9: 13-22.
Boratynska, K. (2002) Needle variability of Pinus mugo in the West Tatra Mts. Dendrobiology 48: 3-8.
Dillon, W. R. and Goldstein, M. (1984) Multivariate analysis methods and applications. John Wiley and Sons, New York.
Dunlap, J. M. and Stettler, R. F. (1995) Genetic variation and productivity of Populus trichocarpa and its hybrids X, trait correlations in young black cottonwood from four river valleys in Washington. Canadian Journal of Forestry Research 25: 1710-1725.
Franiel, I. and Wieski, K. (2005) Leaf features of silver birch (Betula pendula) variability within and between two populations. Trees 19: 81-88.
Gillespie, J. H. (1991) The causes of molecular evolution. Oxford University press, Oxford.
Kovacic, S. and Nikolic, T. (2005) Relations between Betula pendula leaf morphology and environmental factors in five regions of Croatia. Acta biologica cracoviensia series botanica 47: 1-7.
Kovacic, S. and Simic, D. (2001) Intrapopulational and interpopulational relations of Betula pendula in Croatia, based on leaf morphology. Acta biologica cracoviensia series botanica 43: 87-96.
Salehi Shanjani, P., Paule, L., Gomory, D. (2002) Diversity and differentiation in beech forests of Iran. Iranian Journal of Forest and Poplar Research 11: 35-94 (in Persian).
Sokal, R. R. and Rohlf, F. J. (1989) Biometry. WH Freeman and Company, San Francisco.
Thormann, C. E., Ferreira, M. E., Camargo, L. E. A., Tivanga, J. G. and Osborn, T. C. (1994) Comparison of RFLP and RAPD markers to estimating genetic relationships within and among cruciferous species. Theoretical Appliaed Genetics 88: 973-980.
Thorssn, A. T., Palsson, S., Sigurgeirsson, A. and Anamthawat, J. K. (2007) Morphological variation among Betula nana, B. pubescens and their triploid hybrids in Iceland. Annals of Botany 99: 1183-1193.
Tsuda, Y., Sawada, H., Takafumi, O., Katsuhiro, N., Hiroki, N. and Yuji, I. (2010) Landscape genetic structure of Betula maximowicziana in the Chichibu mountain range, central Japan. Tree Genetics and Genomes 6: 377-387.
Ward, J. H. (1963) Hierarchical grouping to optimize an objective function. Journal of the American Statistical Association 48: 236-244.
Williams, J. H. and Arnold, M. L. (2001) Sources of genetic structure in the woody perennial Betula occidentalis. International Journal of Plant Sciences 162: 1097-1109.
Wilsey, B. J., Haukioja, E., Koricheva, J. and Sulkinoja, M. (1998) Leaf fluctuating asymmetry increase with hybridization elevation tree Birches. Ecology 79: 2092-2099.
Zyryanova, O. (2010) White birch trees as resource species of Russia, their distribution, ecophysiological, features and multiple utilizations. Eurasian Journal of Forest Research 13: 25-40.