Prevalence of Honeybee Colony Collapse Disorder and its Relation to Nosema Spp. and Climate in Apiaries of Iran

Document Type : Epidemiology & Public Health

Authors

1 1Graduated From the Faculty of Veterinary Medicine, University of Tehran, Tehran, Iran, Department of Animal Sciences, Kurdistan Agricultural and Natural Resources Research and Education Center, AREEO, Sanandaj, Iran

2 2Department of Food Hygiene and Quality Control, Faculty of Veterinary Medicine, University of Tehran, Tehran, Iran

3 3Department of Honey Bee, Silk Worm and Wildlife, Razi Vaccine and Serum Research Institute, Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran

4 4Department of Parasitology, Faculty of Veterinary Medicine, University of Tehran, Tehran, Iran

5 5DVM, Expert of Directorate of Health and Control of Bird, Honey Bee and Silk Worm, Iran Veterinary Organization, Tehran, Iran

Abstract

BACKGROUND: Colony Collapse Disorder is a mysterious phenomenon in which worker bees abruptly disappear from a beehive.
Objectives: The aim of this study was to estimate the prevalence of CCD and its relation to climate and Nosema spp infections.
Methods: This Cross sectional study was done from April to September 2016. With respect to different climatic zones of the country, a total of 183 ‎apiaries were selected. In each apiary, 5 percent of the colonies were ‎randomly sampled. Adult bee samples were then examined for the presence of Nosema infections using PCR. Data were analyzed by Chi-square using SPSS version 21.0.
Results: The results showed the prevalence of colony collapse disorder in the studied apiaries at 26.8%. The CCD prevalence was ‎20.5% in ‎humid, ‎16.1% in semi humid‎, ‎22.7% in very humid, 38.2% in arid, ‎43.8% in semi-‎arid and ‎16% in Mediterranean conditions. Comparing CCD phenomenon in different climatic regions, there were significant differences (P<0.05). The prevalence of Nosema ceranae infection was 85 (46.4%), however, infection with Nosema apis was not observed in the samples either in pure form or as ‎associated infection. There was no statistical significant difference between symptomatic and asymptomatic apiaries with colony collapse disorder in terms of presence or absence of N.ceranea (P>0.05).
 Conclusions: The results suggest that climate could influence the prevalence of Colony Collapse Disorder. It may be due to different foraging resources in under studied area. According to findings of this study it seems that N.ceranea alone cannot be the cause for CCD. Further studies are needed to clarify the interactions between climate and other possible causes. 

Keywords


مقدمه

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

کاهش کلنی‌های زنبورعسل در مقیاس وسیع برای صنعت زنبورداری جدید نیست، از سال 1869 میلادی گزارش‌های بین‌المللی تلفات غیر معمول موجود است؛ در برخی از آن‌ها توصیف کاهش کلنی مشابه مواردی است که در چندسال اخیر گزارش شده است. کاهش قابل توجه زنبورعسل در گذشته با نام بیماری می(May Disease) ، ناپدید شدن (Disappearing Disease) و بیماری سقوط

 (Fall dwindle Disease) گزارش شده است، برای مثال در بیماری می که در سال‌های 1891 و 1896 میلادی در ایالت کلرادو امریکا اتفاق افتاد؛ دسته‌هایی از زنبوران به‌طور کامل ناپدید شدند یا در مدت زمان کوتاهی جمعیت آن‌ها به‌طور قابل ملاحظه‌ای کاهش یافت (2).

در سال‌های اخیر بسیاری از کلنی‌های مفقود شده با عنوان «ناپدید شده» گزارش شدند، این کلنی‌ها دارای تعداد کمی زنبور زنده یا بدون زنبور بودند (5). این پدیده در آمریکا اختلال فروپاشی کلنی (Colony Collapse Disease) و در اروپا سندرم کاهش کلنی زنبورعسل (Colony Depopulation Syndrome) نامیده شد (13).

آن‌چه امروزه فروپاشی کلنی نامیده می‌شود برای نخستین بار در سال 2006 میلادی توسط یک زنبوردار به‌نام Dave Hackenberg ساکن فلوریدا گزارش شد (2). پس از آن تا سال 2007 میلادی، بیش از 22 ایالت اختلال فروپاشی کلنی را گزارش کردند و برخی زنبورداران تا 95 درصد کندوهای خویش را از دست دادند. روند افزایش فروپاشی کلنی تا سال 2008 میلادی ادامه داشت؛ به‌طوری که 35 ایالت، این اختلال را گزارش کردند (27).

فروپاشی کلنی‌ها با مجموعه‌ی گسترده‌ای از علائم شامل ناپدید شدن ناگهانی زنبوران بالغ به‌ویژه زنبوران کارگر از کندوهایی که ذخیره‌ی غذایی کافی و مراحل مختلف نوزاد دارند، نبود زنبوران تلف‌شده در اطراف و داخل کندوها، تأخیر غارت کندوهای تحت تأثیر قرار گرفته توسط آفت‌های ثانویه یا زنبوران دیگر شناسایی می‌شود (5).

کاهش ناگهانی کلنی‌ها محدود به آمریکا نبوده و از سایر کشورها هم گزارش شده است. نخستین کشور اروپایی درگیر CCD جمهوری ایرلند بود؛ به‌طوری که تا 50 درصد تلفات را در سال 2007 میلادی گزارش کرد و به‌تدریج کشورهای اسپانیا،  پرتغال، ایتالیا، دانمارک و آلمان وجود بیماری را گزارش کردند. در آسیا تایوان نخستین کشوری بود که این بیماری را گزارش کرد، اما امروزه ردپای این پدیده­ در اغلب کشورهای آسیایی مشاهده می‌شود (16). در ایران نیز از سال 1388 گزارش‌های پراکنده‌ای از کاهش ناگهانی جمعیت کلنی‌های زنبورعسل از مناطق مختلف کشور وجود دارد.

از آنجایی‌که توزیع جغرافیایی، فراوانی و وسعت کاهش ناگهانی زنبوران عسل رو به افزایش است، فرضیه‌های مختلفی برای توجیه این پدیده مطرح شده است (19). اگرچه علت دقیق فروپاشی کلنی‌های زنبورعسل در بسیاری از موارد ناشناخته مانده است، اما عوامل بسیاری به‌عنوان علل آن مطرح می‌باشند، به‌طوری که در سال‌های اخیر مواردی از کاهش جمعیت کلنی به آلودگی با نوزما سرانه نسبت داده شده است (4، 8، 9، 14، 18و19). هم‌چنین در برخی مطالعه‌ها اثرات اقلیم بر فروپاشی کلنی‌ها مطرح شده است (6). بر این اساس، این پژوهش با هدف برآورد شیوع این پدیده در زنبورستان‌ها و بررسی ارتباط آن با آلودگی به گونه‌های نوزما و اقلیم با همکاری دفتر بهداشت و مدیریت بیماری‌های طیور، زنبورعسل و کرم ابریشم سازمان دامپزشکی کشور در تعدادی از زنبورستان‌های کشور انجام شد.

 

مواد و روش کار

- روش نمونه‌برداری: این مطالعه  به‌صورت مقطعی از فروردین تا مهر ماه 1395 به مدت 6 ماه انجام شد. حجم نمونه‌ی مورد نیاز برای انجام این مطالعه با احتساب شیوع CCD برابر 44 درصد (بر اساس مطالعه­ی Mahmodi و همکاران در سال 2011)، سطح معنی‌داری 5 درصد و دقت برابر20 درصد شیوع، 122 زنبورستان بدست آمد. بر اساس نمونه­گیری خوشه­ای هر استان یک خوشه و زنبورستانها بعنوان عناصر خوشه در نظر گرفته شدند. ضریب اثر طرح خوشه­ای( 1/5) در حجم نمونه فوق (122) ضرب و تعداد کل نمونه­ها 183 زنبورستان برآورد شد. تعداد نمونه­ها به نسبت زنبورستانهای هر استان به آن اختصاص یافت. از زنبورستان­های استان­های آذربایجان­غربی، خراسان­شمالی، اردبیل، چهارمحال­و­بختیاری، گلستان، خراسان­جنوبی، سمنان، فارس، خراسان­رضوی، قزوین، قم، کردستان، لرستان، کرمان­جنوب، گیلان، مازندران و یزد نمونه­گیری شد (جدول 1). با استفاده از داده­های سیستم اطلاعات جغرافیایی (داده­های سازمان دامپزشکی) پراکنش زنبورستانهای مورد مطالعه تعیین و اقلیم زنبورستانهای نمونه­گیری شده بر اساس نقشه طبقه­بندی اقلیمی دومارتن (20) مشخص شد. به منظور نمونه­گیری، ضمن بازدید کلنی‌های زنبورستان مورد مطالعه توسط کارشناس بیماری‌های زنبورعسل اداره‌های دامپزشکی استان‌ها، اطلاعات مربوط به زنبورستان و مشاهده یا عدم مشاهده­ی علائم CCD در فرم مربوط به هر زنبورستان ثبت شد. سپس بر اساس دستورالعمل سازمان دامپزشکی، به‌صورت تصادفی از 5 درصد کلنی‌های هر زنبورستان نمونه­ی زنبور بالغ اخذ (29) و اطلاعات مربوط به نمونه روی ظروف نمونه‌گیری ثبت شد. نمونه‌ها تحت شرایط سرد برای تشخیص آلودگی به گونه‌های نوزما (نوزما آپیس و سرانه) به آزمایشگاه تشخیص بیماری‌های زنبورعسل مؤسسه­ تحقیقات واکسن و سرم‌سازی رازی منتقل شدند.

-آزمایش مولکولی تشخیص اسپور نوزما

 آماده‌سازی نمونه‌ها: ابتدا تعداد 20 عدد زنبورعسل از نمونه­ی هر زنبورستان شمارش و جدا شد. سپس قسمت شکمی زنبوران جدا و در mm 10 آب مقطر کاملاً له شد، سوسپانسیون حاصل فیلتر و در 800 دور به‌مدت 6 دقیقه سانتریفیوژ شد.

    استخراج DNA : با استفاده از کیت استخراج DNA ساخت شرکت تکاپوزیست مطابق دستورالعمل کارخانه‌ی سازنده DNA استخراج شد.

PCR:غنی‌سازی با استفاده از کیت PCR شرکت سیناکلون بر اساس دستورالعمل OIE (22) در ترموسایکلر ساخت کارخانه­ی اپندورف آلمان انجام شد. سپس محصولات PCR با الکتروفورز بر روی ژل آگاروز 1 درصد جدا، رنگ‌آمیزی و با اشعه‌ی ماوراءبنفش مشاهده شد. در این روش از پرایمر اختصاصی 16S rRNA  برای شناسایی نوزما آپیس و نوزما سرانه استفاده شد (جدول 2). کنترل‌های مثبت برای نوزما آپیس و نوزما سرانه در مؤسسه­ تحقیقات واکسن و سرم‌سازی رازی تهیه شدند.

-جمع‌آوری و آنالیز داده‌ها:اطلاعات نمونه‌های اخذ شده و نتایج آزمایش مولکولی نمونه‌ها ثبت و داده‌ها با استفاده از آزمون‌ غیر پارامتریک مربع کای تحلیل شدند. تمامی آنالیزهای آماری با نرم‌افزار SPSS نسخه‌ی 21 انجام شد. در تمام آنالیزها (p≤0.05) به عنوان سطح معنی‌داری در نظر گرفته شد.

نتایج

   در این مطالعه از تعداد 183 زنبورستان نمونه‌برداری شده، 49 زنبورستان (26/8 درصد، فاصله‌ی اطمینان 95 درصد: 33/3-20/3) دارای علائم CCD بودند. میزان شیوع فروپاشی کلنی‌های زنبور عسل در اقلیم‌های متفاوت کشور در جدول شماره 3 نشان داده شده است. در مقایسه‌ی میزان شیوع CCD با استفاده از آزمون مربع کای، شیوع آلودگی به CCD در مناطق اقلیمی کشور تفاوت آماری معنی‌داری داشت (0.05>p).

بر اساس نتایج آزمایش مولکولی آلودگی با نوزما سرانه در همه‌ی اقلیم‌های مورد مطالعه وجود داشت (جدول شماره 4) در حالی­که آلودگی با نوزما آپیس در هیچ‌یک از نمونه‌ها مشاهده نشد. شیوع نوزما سرانه در مناطق اقلیمی کشور از نظر آماری تفاوت معنی‌داری داشت (0.05>p).

همان‌طور که در جدول شماره 5 مشاهده می‌شود، بین آلودگی با نوزما سرانه و شیوع CCD ارتباط آماری معنی‌داری دیده نشد (p>0.05) به‌طوری که آلودگی با این انگل شانس ابتلا به پدیده‌ی CCD را تغییر نداده است.

بحث

بر اساس یافته‌های این مطالعه شیوع CCD در زنبورستان‌های مورد مطالعه 26/8 درصد بود. در مطالعه‌ی Mahmodi و همکاران در سال 2011، شیوع فروپاشی کلنی در زنبورستان‌های ایلام 43/3 درصد گزارش شده است (18). در این مطالعه با توجه به شیوع بیشتر این پدیده در مناطق نیمه خشک (43/8 درصد) و خشک (38/2 درصد) نسبت به سایر مناطق، احتمال تأثیر گرما بر وقوع بیشتر پدیده­ی CCD مطرح است. بر اساس مطالعه‌‍‌های انجام شده اقلیم محلی که کندو در آن قرار گرفته بر سلامت کندو تأثیر دارد. تغییرات اقلیمی در سطوح مختلف بر زنبورعسل اثر می‌گذارند، این تغییرات ممکن است بر رفتار و فیزیولوژی زنبور تأثیر مستقیم داشته باشد، کیفیت محیط گیاهان را تغییر دهد، توان تولید کلنی را کاهش یا افزایش دهد، توزیع زنبوران را تعیین کند و باعث افزایش رقابت در بین گونه‌ها و نژادهای زنبور و نیز در بین پارازیت‌ها و پاتوژن‌های آن‌ها شود. بنابراین تغییرات محیطی و استرس‌های مربوط به تغییرات اقلیمی ممکن است در فروپاشی کلنی‌ها تأثیر داشته باشند (6).

بنابر نظر Conte و همکاران در سال 2008 اقلیم بر توزیع محصولات کشاورزی، رشد گل‌ها و تولید شهد و گرده گیاهان اثر می‌گذارد و این‌ها بر فعالیت چرای زنبورها تأثیر می‌گذارند. تغییرات اقلیم می‌تواند روی گسترش و حدت آفت‌ها و پاتوژ‌ن‌ها اثر بگذارد. هم‌چنین حساسیت زنبور عسل به پاتوژن‌ها و پارازیت‌ها تحت تاثیر هوا قرار می‌گیرد؛ مثلاً در زمستان‌های گرم مانند سال‌های 2012-2011 کاهش کلنی‌ها در آمریکا به‌طور قابل ملاحظه‌ای کمتر بوده است (6).

 یافته‌های این مطالعه نشان دهنده‌ی گسترش نوزما سرانه در 85 نمونه از 183 نمونه مطالعه شده (46/4 درصد) و عدم حضور نوزما آپیس است. نوزما سرانه نخستین بار در ایران در سال 2011 میلادی توسط Nabian و همکاران گزارش شد (21). به دنبال آن آلودگی با این انگل  توسط مطالعه‌های  دیگری تأیید شد. به‌طوری که بر اساس یافته‌های Razmaraii و همکاران در سال 2011، نوزما سرانه تنها گونه جدا شده در آذربایجان شرقی بود (26). هم‌چنین بر اساس گزارش Modirrousta و همکاران در سال 2014، نمونه‌های جمع آوری شده‌ی 5 استان در فواصل زمانی 2013-2004 از نظر آلودگی به نوزما سرانه مثبت بودند (20). هم‌چنین در مطالعه Aroee و همکاران در سال 2017، آلودگی نمونه‌های زنبورعسل استان‌های اصفهان، چهارمحال و بختیاری و فارس با نوزما سرانه گزارش شده است (1).

در این مطالعه‎ بیشتر بودن شیوع نوزما سرانه در مناطق نیمه مرطوب (71 درصد)، خیلی مرطوب (68/2درصد) و مرطوب (53/8 درصد) نسبت به سایر مناطق می‌تواند بر تأثیر رطوبت در گسترش آلودگی تأکید داشته باشد. نتایج مطالعه‌ی انجام شده در ترکیه نیز نشان می‌دهد نوزما سرانه در مناطقی که رطوبت بالایی دارند، گسترش بیشتری دارد و میزان رطوبت به‌دلیل تأثیر آن بر دما، یک عامل حیاتی برای گونه‌های نوزما توصیف شده است (23).

در این مطالعه نسبت شانس CCD در زنبورستان‌های آلوده به نوزما سرانه (1/59-0/42) 0/82 با 0/56p=، نشان می‌دهد آلودگی به نوزما سرانه تأثیر آماری معنی‌داری در شانس ابتلا به CCD ندارد. بنابراین آلودگی با نوزما به‌عنوان یک عامل مؤثر بر CCD مطرح نیست. مشابه یافته‌ی اخیر در مطالعه‌ی انجام شده در اردن توسط Haddad  در سال 2011، ارتباطی بین فروپاشی کلنی و نوزما یافت نشده است (10). هم‌چنین مطابق مطالعات متعددی نوزما به تنهایی علتCCD  در ایالات متحده امریکا، آلمان و کانادا نبوده است (15). در مطالعه‌ی Cox-Foster و همکاران در سال 2007، CCD با نوزما سرانه ارتباط علیتی نداشته است. در این مطالعه ویروس فلجی حاد اسرائیلی خطر فروپاشی کلنی را افزایش داده (Odds=65)  و تأثیر این ویروس در نمونه‌های آلوده با نوزما آپیس افزایش داشته در حالی­که نوزما سرانه تأثیر ویروس فلجی حاد اسرائیلی را بر CCD تغییر نداده است (7).

در برخی مطالعه‌ها تأثیر آلودگی به نوزما بر فروپاشی کلنی‌های زنبورعسل گزارش شده است، به‌طور مثال در مطالعه‌ی Mahmodi و همکاران در سال 2011، از 13 زنبورستان دارای علائم CCD، 9 زنبوررستان (69/23 درصد) آلوده به اسپور نوزما بوده، به‌طوری که ارتباط آماری معنی‌داری بین CCD و نوزما وجود داشته است (18).

بر اساس 3 منبع منتشر شده در مورد کلنی‌های زنبور در اسپانیا، نوزما سرانه در مدت 18 ماه آلودگی کندو، موجب CCD می‌شود (13، 19). اگرچه به‌نظر می رسد نوزما سرانه نقش کلیدی در فروپاشی کلنی در اسپانیا داشته باشد (12) اما جای تعجب است که با وجود نوزما سرانه در قاره اروپا از سال 1998 میلادی، تلفات بیشتری از اروپا گزارش نشده است. بر اساس نظر پژوهشگران ممکن است آپیس ملیفرای اسپانیایی نسبت به سایر نژادهای زنبور به نوزما سرانه حساس‌تر باشد یا سویه‌ی نوزما سرانه‌ی اسپانیا حدت بیشتری داشته باشد (25).

در مطالعه‌ی دیگری در اسپانیا در 2 زنبورستان واقع در دو منطقه که از نظر آب و هوایی و شرایط پرورش زنبور از جمله نوع پوشش گیاهی موجود برای زنبوران و روش‌های زنبورداری کاملاً متفاوت بودند، خسارت‌های زیاد کلنی با ناپدید شدن زنبوران بالغ در زمستان سال 2006 میلادی گزارش شده است. بر اساس یافته‌های این مطالعه، تنها پاتوژن موجود در همه‌ی نمونه‌های زنبور عسل نوزما سرانه بوده است (23). هم‌چنین مرگ کلنی‌ها درلهستان، اتریش و فرانسه و کاهش چشمگیر کلنی در اروپای شمالی به وجود نوزما سرانه مربوط بوده است (28).

امروزه نوپدید بودن نوزما سرانه برای زنبور آپیس ملیفرا و اثرات آن بر کلنی‌ها مشخص است. این‌که این پاتوژن به تنهایی یا همراه سایر عوامل از جمله از بین رفتن سکونت‌گاه‌های زنبور و کمبود منابع گیاهی موجب فروپاشی کلنی‌ها می‌شود؛ هنوز نامشخص است و نیاز به بررسی‌های بیشتر دارد.

نتیجه‌گیری

بنا به اطلاعات نویسندگان، این پژوهش نخستین بررسی گسترده‌ی کلنی‌های زنبورعسل در کشور است که شواهدی از پدیده‌ای به‌نام CCD را فراهم نموده است. تأثیر اقلیم بر این پدیده در این پژوهش نیز می‌تواند ناشی از اثرات متقابل عوامل دیگر باشد. به هرحال تغییرات اقلیمی در سال‌های اخیر و تأثیر آن بر منابع چرای در دسترس زنبوران، بر رفتار، عملکرد و بیماری‌های زنبورعسل مؤثر بوده است. در این مطالعه اگرچه نوزما سرانه در تعدادی از موارد فروپاشی کلنی دیده شد، اما آلودگی با این انگل تأثیر معنی‌داری در شانس ابتلا به CCD نداشت. بنابراین ممکن است نوزما سرانه همراه با سایر انگل‌ها و عوامل بیماری‌زا در فروپاشی کلنی‌ها نقش داشته باشد. هم‌چنین عدم مشاهده نوزما آپیس و پراکنش نوزما سرانه در زنبورستان‌ها می‌تواند نشانه‌ی جایگزینی نوزما سرانه با نوزما آپیس باشد. بر اساس این یافته ها و با توجه به این‌که بر اساس مطالعه‌های اخیر پژوهشگران تأثیر متقابل عوامل محیطی، مدیریتی، پاتوژن‌ها و پارازیت‌ها در پدیده‌ی فروپاشی کلنی محتمل است، بررسی نقش این عوامل در فروپاشی کلنی­ها در زنبورستان‌های کشور مفید خواهد بود.

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

بدین‌وسیله ازهمکاری سازمان دامپزشکی کشور برای انجام هماهنگی‌های لازم برای نمونه‌گیری در استان‌های مختلف قدردانی می‌شود. از تمامی دامپزشکان که در کار فیلدی این مطالعه همکاری نمودند و همه‌ی زنبوردارانی که از زنبورستان‌های آن‌ها نمونه‌گیری شد، تشکر و قدردانی می‌شود. هم‌چنین از زحمات کارکنان آزمایشگاه تشخیص بیماری‌های زنبورعسل مؤسسه تحقیقات واکسن و سرم‌سازی رازی به‌دلیل انجام آزمایش‌های مربوط قدردانی می‌شود.

 

 

Aroee, F., Azizi, H., Shiran, B., Pirali Kheirabadi, K. (2017). Molecular identification of Nosema species in provinces of Fars, Chaharmahal and Bakhtiari and Isfahan (Southwestern Iran). Asian Pac J Trop Biomed, 7, 10-13. https://doi.org/10.1016/j.apjtb.2016.11.004
Benjamin, A., McCallum, B. (2009). A world without bees. The mysterious decline of the honeybee- and what is means for us. Guardian Books, London, UK.
Biesmeijer, J.C., Roberts, S.P.M., Reemer, M., Ohlemüller, R., Edwards, M., Peeters, T., Schaffers, A.P., Potts, S.G., Kleukers, R., Thomas, C.D., Settele, J., Kunin, W.E. (2006). Parallel declines in pollinators and insect-pollinated plants in Britain and the Netherlands. Science, 313 (5785), 351-354. https://doi.org/ 10.1126/science.1127863
Cepero, A., Martín-Hernández, R., Bartolomé, C., Gómez-Moracho, T., Barrios, L., Bernal, J., Teresa Martín, M., Meana, A., Higes, M. (2015). Passive laboratory surveillance in Spain: pathogens as risk factors for honey bee colony collapse. J Apic Res, 54, 525-531. https://doi.org/10.1080/00218839.2016.1162978
Chensheng, LU., Warchol, KM., Callahan, RA. (2012). In situ replication of honey bee colony collapse disorder. Bull Insectol, 65(1), 99-106.
Conte, Y. Le., Navajas, M. (2008). Climate change: impact on honey bee populations and diseases. Rev Sci Tech, 27 (2), 499-510
Cox-Foster, D.L., Conlan, S., Holmes, E.C., Palacios, G., Evans, J.D., Moran, N.A., Quan, P.L., Briese, T., Hornig, M., Geiser, D.M. Martinson, V., VanEngelsdorp, D., Kalkstein, A.L., Drysdale, A., Hui, J., Zhai, J., Cui, L., Hutchison, S.K., Simons, J.F., Egholm, M., Pettis, J.S., Lipkin, W.I. (2007). A metagenomic survey of microbes in honey bee colony collapse disorder. Science, (New York, N.Y.) 318, 283-287. http://dx.doi.org/10.1126/science.1146498
Cornman, R.S., Chen, Y.P., Schatz, M.C., Street, C., Zhao, Y., Desany, B., Egholm, M., Hutchison, S., Pettis, J.S., Lipkin, W.L., Evans, J.D. (2009). Genomic analyses of the microsporidian Nosema ceranae, an emergent pathogen of honey bees. PLoS Pathogens, 5 (6), 1-14. https://doi.org/10.1371/journal.ppat.1000466
Gajger, I.T., Vugrek, O., Grilec, D., Petrinec, Z. (2010). Prevalence and distribution of Nosema ceranae in Croatian honeybee colonies. Vet Med, 55 (9), 457-462.
Haddad, N. (2011). Honey bee viruses, diseases and hive management in the Middle East and their relation to the Colony Collapse Disorder and bee losses. Uludag Bee J, 11 (1), 17-24.
Higes, M., Martín-Hernández, R., Meana, A. (2006). Nosema ceranae, a new microsporidian parasite in honeybees in Europe. J Invert Pathol, 92, 93-95. https://doi.org/10.1016/j.jip.2006.02.005
Higes, M., Martín-Hernández, R., Botías, C., Garrido Bailón, E., González-Porto, A.V., Barrios, L., Nozal, M.J.D., Bernal, J.L., Jimenez, J.J., Garcia Palencia, P., Meana, A. (2008).  How natural infection by Nosema ceranae causes honeybee colony collapse. Environ. Microbiol, 10(10), 2659-2669. http://dx.doi.org/10.1111/j.1462-2920.2008.01687.x
Higes, M., Martín-Hernández, R., Garrido-Bailón, E., González-Porto, A.V., García-Palencia, P., Meana, A., del Nozal, M.J., Mayo, R., Bernal, J.L. (2009). Honeybee colony collapse due to Nosema ceranae in professional apiaries. Environ Microbiol Rep, 1(2), 110-113. http://dx.doi.org/10.1111/j.1758-2229.2009.00014.x
Higes, M., Meana, A., Bartolomé, C., Botías, C., Martín-Hernández, R. (2013). Nosema ceranae (Microsporidia), a controversial 21st century honey bee pathogen. Environ Microbiol Rep, 5, 17-29. https://doi.org/10.1111/1758-2229.12024
Huang, Z. (2010). Honey Bee Nutrition. Am.  Bee J, 150, 773-776.
Jacobsen, R. (2010). Fruitless fall: The Collapse of the Honey Bee and the Coming Agricultural Crisis. Bloomsbury Publishing. New York, USA. p. 66-80.
Kluser, S., Peduzzi, P. (2007). Global Pollinator Decline: A Literature Review. UNEP/GRID-Europ, Switzerland. p. 1-12.
Mahmodi, M., Basami, Sh., Nabian, S., Bahonar, A.R. (2011). Study on some possible factors in occurrence of colony collaps disorder in Ilam. HBSJ, 5 (5), 4-9.
Martín-Hernández, R., Meana, A., Prieto, L., Salvador, A.M., Garrido-Bailón, E., Higes, M., (2007). Outcome of colonization of Apis mellifera by Nosema ceranae. Appl Environ Microbiol, 73, 6331-6338.
Masodian, A., Kaviani, M. (2007). Climatology of Iran. University of Isfahan press. Isfahan, Iran. p.179.
Modirrousta, H., Moharrami, M., Mansouri, M.A. (2014). Retrospective study of the Nosema ceranae infection of honey bee colonies in Iran (2004-2013). Arch Razi Inst, 69, 197-200. DOI: 10.7508/ARI.2014.02.012
Nabian, S., Ahmadi, K., Shirazi, M.N., Sadeghian, A.G. (2011). First detection of Nosema ceranae, a microsporidian protozoa of European honeybees (Apis mellifera) in Iran. Iran J Parasitol, 6, 89. PMID:  22347302
Office International Des Epizooties (2008). Nosemosis of Honeybees. Manual of Diagnostic Tests and Vaccines for Terrestrial Animals.
Özgör, E., Güzerin, E., Keskin, N. (2015). Determination and comparison of Nosema apis and Nosema ceranae in terms of geographic and climatic factors. Hacettepe J Biol Chem, 43(1): 9-15.
Partap, U., Partap, T., Yoghua, H. (2001). Pollination failure in apple crop and farmers management strategies in Hengduan, China. Acta Hortic, 561: 225-230.
Paxton, R.J. (2010). Does infection by Nosema ceranae cause. Colony Collapse Disorder. in honey bees (Apis mellifera)? J Apic Res, 49(1), 80-84. https://doi.org/10.3896/IBRA.1.49.1.11
Razmaraii, N., Sadegh-Eteghad, S., Babaei, H., Paykari, H., Esmaeilnia, K., Froghy, L. (2013). Molecular identification of Nosema species in East Azerbaijan province, Iran. Arch Razi Inst, 68 (1), 23-27.
Schacker, M. (2008). A spring without Bees: How Colony Collapse Disorder Has Endangered Our Food Supply. Schacker, M. (2008). A spring without Bees: How Colony Collapse Disorder Has Endangered Our Food Supply. Guilford, Conn. Lyons Press, USA. p.19
Tehrani, F., Rajab, A., Tabatabaaei, SM., Forci, M., Ghafouri, S.A., Amirhajlo, S., Khosravi, D., Hashemi, A., Khoshnevisan, Sh., Charkhkar, S., Fallah, M.H., Moghadas, E. (2014). Executive plan of survey and control of Poultry and Bee Diseases. Veterinary Organization Executive plan, p. 684-693.
Topolska, G., Gajda, A., Hartwig, A. (2008). Polish honey bee colony-loss during the winter of 2007/2008, J Apic Sci, 52 (2), 95-104.
VanEngelsdorp, D., Underwood, R., Caron, D., Hayes, J. (2007). An estimate of managed colony losses in the winter of 2006-2007: a report commissioned by the apiary inspectors of America. Am Bee J, 147(7), 599-603