مطالعه بیان ژن افتراقی زنبور عسل ملکه، نر و کارگر با استفاده از داده‌های rnaseq

این پژوهش با هدف مطالعه پروفایل بیان ژن و تعیین ژن‌های شاخص در تمایز و تکامل ملکه، نر و کارگر با مقایسه تفریقی آن‌ها در سن 5 لاروی یا همان سن تمایز انجام شد. لذا ترانسکریپتوم (توالی کل mrna) 15 نمونه از زنبور عسل نژاد ایتالیایی (a. m. ligustica) شامل 5 زنبور نر، 5 کارگر و 5 ملکه از طریق همردیفی و مکان یابی خوانش های rnaseq بر روی ژنوم مرجع زنبور عسل نسخه amel_4.5 مرتب شد. سپس آنالیز بیان افتراقی ژن صورت گرفت و در نهایت 15962 ژن و 31297 ایزوفرم بر روی ترانسکریپتوم این نمونه‌ها مشخص شد و نهایتاً 465 ژن با بیان متفاوت بین تیمارهای نر و ملکه، 495 ژن بین تیمارهای کارگر و ملکه و 764 ژن بین تیمارهای نر و کارگر بود (p<0.000005)، (p<0.0001) و بیشترین تفاوت بیان ژنی بین تیمارهای کارگر و نر، ژن‌های gb45614 و gb42053 با لگاریتم 2 fold change به‌ترتیب با مقادیر 10 و 11.5 به‌دست آمدند و در مقایسه تیمارهای نر و ملکه ژن‌های gb45614 و gb48020 به‌ترتیب با مقادیر 11.7 و 11.8 و در مقایسه تیمارهای ملکه و کارگر ژن‌های gb43508 و gb42053 به‌ترتیب با مقادیر 6.6 و 9 با بیشترین تفاوت بیان ژنی مشاهده شدند. آنالیز ماهیت ژن (go) و مسیرهایی که در آن ها درگیر هستند نشان داد، که بررسی دقیقی و جامعی بر روی تعداد زیادی از این ژن‌ها انجام نشده است، ولی تعداد زیادی از این ژن‌های شاخص در مقایسه تیمارهای: ملکه و نر مرتبط با اجزای ساختاری و جدایی‌ناپذیر غشا، پردازش متابولیکی کیتین، بازدارنده‌های تریپسین، فعالیت گیرنده گاسترین، باندشونده با هورمون جوانی، پاسخ‌های دفاعی نسبت به باکتری‌ها، پروتئین‌های کوتیکول شفیره هستند درصورتی‌که در مقایسه ژن‌های شاخص متمایز در تیمارهای ملکه و کارگر مرتبط با مسیرهای متابولیکی، متابولیسم داروها و سایر آنزیم ها، متابولیسم لیپیدها، باندشونده با نوکلئیک اسید، انتقال‌دهنده کلسترول داخل سلولی، پردازش‌های متابولیکی کیتین، بیوسنتز یوبیکینون و سایر ترپنوئید کینون‌ها و ژن‌های متمایز در مقایسه کارگر و نر مرتبط ا اجزای ساختاری و جدایی‌ناپذیری از غشای متابولیسم و انتقال آمینو اسیدها، باند شونده با یون کلسیم، باندشونده با فلز یون، انتقال‌دهنده کلسترول کیتین درون سلول، متابولیسم بودند.

Study of differential gene expression in queen, drone and worker honey bee using RNAseq data

The aim of this study was to investigate the gene expression profile and find the important genes in the differentiation and evolution of the queen, worker and drone honey bee at Stage 5 Larvae. Therefore, transcripts (total mRNA sequence) of 15 samples of Italian honey bee (A. m. ligustica) including 5 drone, 5 worker and 5 queen bees were aligned to the reference genome of the honey bee. In gene expression analysis of RNASeq data, 15962 genes and 31297 isoforms were identified. In our differential gene expression analyses, 465 genes between drone and queen bees, 495 genes between worker and queen bees and 764 genes between drone and worker bees, were expressed differently (P <0.000005). The largest difference in expression of genes was observed between drone and workers were for GB45614 and GB42053, with log2 fold change that was 10 and 11.5, respectively. In drone and queen bees comparison, GB45614 and GB48020 genes with log2 fold change 7.11 and 8.11, and in queen and worker bees comparison, GB43508 and GB42053 with log2 fold change 6.6 and 9, had the largest difference in gene expression, respectively. The analysis of the gene ontology (GO) and the pathways involved showed that the function of many of these genes has yet to be found. However, a large number of the genes expressed defiantly in drone and queen bees were related to integral component of membrane, calcium ion binding, carboxypeptidase, cholecystokinin receptor,chitin metabolic process, chymotrypsin inhibitor, haemolymph juvenile hormone binding and pupal cuticle protein, while differentially expressed genes in queen and worker bees comparison were related to metabolic pathways, enzymes metabolism, Pyrimidine metabolism, lipid metabolism, protein kinase, ATP binding and nucleic acid, intracellular cholesterol transport, chitin metabolic process, nitrogen compound metabolic process, hydrolase activity and chitin binding. The genes expressed at different levels in worker bees and drones were related to structural elements, the metabolism membrane and transfer of amino acids, calciumionbound, ionbinder, intracellular cholesterol transport and chitin metabolism.