نوع مقاله : علمی - پژوهشی
نویسندگان
1 گروه زیست شناسی، دانشکده علوم پایه، دانشگاه گیلان، رشت ، ایران
2 گروه زیست شناسی، دانشکده علوم پایه، دانشگاه گیلان،
3 گروه زنان و زایمان، مرکز تحقیقات بهداشت باروری، بیمارستان الزهرا، دانشکده پزشکی، دانشگاه علوم پزشکی گیلان، رشت، ایران.
تازه های تحقیق
-
کلیدواژهها
موضوعات
عنوان مقاله English
نویسندگان English
Introduction: Numerous studies on microRNAs (miRNAs) have indicated that abnormal interactions between miRNAs and mRNAs contribute to the development and worsening of endometriosis. miR-223 influences genes associated with various biological functions such as signaling pathways, gene expression, cellular development, proliferation, angiogenesis, and programmed cell death. Alterations in genes have been demonstrated to be significant in the progression of endometriosis. The p53 gene, recognized as a tumor suppressor, plays a vital role in inhibiting cancer by regulating cell growth and division. Recent investigations suggest a potential connection between the p53 gene and endometriosis. In addition, Forkhead box protein O1 (FOXO1) acts as a cell-specific core transcription factor essential for effective endometrial remodeling throughout the menstrual cycle and may play an important role in the onset of endometriosis.
Aims: This study aimed to explore the expression levels of miR-223-3p and its target genes, TP53 and FOXO1, in cases of ovarian endometriosis through bioinformatic analysis and experimental verification.
Materials and methods: Gene expression data were obtained from the GSE105765 dataset in the Gene Expression Omnibus (GEO) database, which encompasses both normal (eutopic) and ectopic tissue samples from endometriosis patients. Differentially, expressed genes were pinpointed using the "limma" package in R software, with the results visualized as a heat map. Among the various expressed miRNAs, miR-223-3p was chosen for further investigation. Subsequently, the target genes of miR-223-3p were explored utilizing the miRDB database, leading to the identification of FOXO1 and TP53 as target genes. In the experimental component of this study, 40 ectopic tissue samples, 40 eutopic tissue samples from individuals with endometriosis, and 40 normal endometrial samples (control group) were analyzed. Following the extraction of total RNA and synthesis of cDNA, the expression levels of miR-223-3p, TP53, and FOXO1 were assessed using Real-Time PCR.
Results: The average relative expression of miR-223-3p in the ectopic group was measured at 4.13±0.72, while in the eutopic group it was 1.10±0.10, and in the control group it was 1.00±0.09. The relative expression of miR-223-3p was significantly higher in the ectopic group compared to both the eutopic and control groups (p<0.0001), whereas no significant differences were found between the control and eutopic groups (p=0.54). The average expression of FOXO1 in the ectopic group (0.47±0.12) was significantly lower than that in the eutopic group (1.02±0.06) and the control group (1.001±0.11). One-way ANOVA revealed a significant reduction in FOXO1 expression in the ectopic group when compared to the eutopic group and the control group (p<0.0001). However, there was no significant difference in FOXO1 expression between the eutopic and control groups (p=0.53). The relative expression level of the TP53 gene was analyzed in the ectopic group (0.335±0.14), the eutopic group (1.05±0.09), and the control group (1.01±0.12). Statistical evaluations showed that TP53 gene expression was significantly lower in the ectopic group compared to both the eutopic and control groups (p<0.0001). Conversely, no significant difference was observed in TP53 expression between the eutopic and control groups (p=0.28).
Conclusion: The findings indicate that miR-223-3p expression in ectopic samples was significantly elevated compared to both the eutopic and the control samples, while the expression of the FOXO1 and TP53 genes in the ectopic group was notably reduced relative to the eutopic and control groups. These results suggest that the expressions of miR-223-3p, TP53, and FOXO1 genes are changed in the endometrium of individuals with endometriosis, highlighting their potential roles in the disease's pathogenesis and their therapeutic implications.
کلیدواژهها English
آندومتریوز، یک بیماری مزمن، التهابی و وابسته به استروژن است که ناشی از قرارگیری و رشد بافت آندومتر (شامل غدد اپیتلیالی و سلولهای استرومایی) در خارج از حفره رحمی میباشد. بافت آندومتر نابهجا (اکتوپیک) معمولا در حفره لگن یافت میشود، اما ممکن است در تخمدانها، لولههای فالوپ، قسمت بالایی شکم و سایر مناطق بدن نیز دیده شود (1). بهطور کلی، تخمین زده میشود که حدود 10 درصد از زنان در سنین باروری به این بیماری مبتلا هستند (2). آندومتریوز بر اساس آسیب شناسی بافتی و مکانهای بروز خود به سه فنوتیپ مختلف تقسیمبندی میشود که شامل آندومتریوز سطحی صفاقی، آندومتریومای تخمدانی و آندومتریوز عمیق نفوذی میباشد (3). این بیماری منجر به قاعدگی نامنظم، درد مزمن لگنی و ناباروری میشود (4). آندومتریوز به عنوان یک بیماری چند عاملی در نظر گرفته میشود زیرا فرایند ایجاد این بیماری (پاتوژنز) شامل طیفی از ویژگیها نظیر وجود بافت آندومتر نابهجا، اختلال در سیستم ایمنی، عدم تعادل بین تکثیر سلولی و مرگ سلولی برنامه ریزی شده (آپوپتوزیس)، سیگنالدهی غددی غیرطبیعی و عوامل ژنتیکی است (5). بهطور کلی، اعتقاد بر این است که بروز و پیشرفت آندومتریوز یک فرآیند بیولوژیکی پیچیده است که در آن عوامل محیطی و ژنهای زیادی دخیل میباشند. در سالهای اخیر، گروهی از ژنها که محصولات آنها RNAهایی با نقشهای تنظیمی هستند، مورد توجه محققان قرار گرفتهاند. RNAهای غیرکدکننده قطعاتی از RNA هستند که از DNA رونویسی میشوند، اما به پروتئین ترجمه نمیشوند (6). میکروRNAها (miRs=miRNAs)، RNAهای غیرکدکننده کوتاهی هستند که طولی در حدود 20 تا 24 نوکلئوتید دارند که بیان ژن را با کاهش ترجمه یا کاهش پایداری mRNA هدف تنظیم میکنند (7). با توجه به اینکه آندومتریوز یک بیماری چند عاملی میباشد، نتایج مطالعات حاکی از اهمیت نقش RNAهای غیر کدکننده از جمله miRNAها در طی روند ایجاد و پیشرفت ضایعات آندومتریوزی است. تحقیقات نشان دادهاند که میزان بیان برخی از miRNAها در بافت آندومتر یوتوپیک (طبیعی) زنان مبتلا به آندومتریوز، در مقایسه با بافت اکتوپیک، دچار تغییراتی میشود (8). miR-223 یک miRNA بسیار حفاظت شده در مهره داران است و ژن آن در جایگاه q12 کروموزوم X قرار دارد. سنتز miR-223 در مرحله اول منجر به ایجاد رشته miR-223-3p میشود، که miR-223 نیز نامیده میشود (9). این miRNA، ژنهایی که در فرآیندهای بیولوژیکی مختلفی مانند انتقال سیگنال، رونویسی، تکوین سلولی، تکثیر سلولی و آپوپتوزیس دخالت دارند را تنظیم میکند (10). مطالعات اخیر همچنین نشان دادهاند که miR-223 در التهاب، عفونت و پیشرفت سرطان نقش دارد (11). علاوه بر این miR-223 در فرآیند انتقال اپیتلیال به مزانشیمی (Epithelial-mesenchymal transition=EMT) در بیماریهای مختلف نظیر آندومتریوز دخالت دارد (12 و 13).
miRNAها، بیان ژن را با اتصال به mRNA تنظیم میکنند که منجر به تخریب mRNA یا سرکوب ترجمه میشود. این فرآیند معمولا از طریق کمپلکس خاموشسازی القا شده-RNA (RNA-induced silencing complex =RISC) رخ میدهد که miRNA آن را به mRNA هدف خود هدایت میکند. در حالیکه بیشتر تعاملات در ناحیه ترجمه نشده '3 UTR(mRNA) رخ میدهد، اتصال میتواند در 'UTR5، نواحی کدکننده و حتی پروموترهای ژن نیز رخ دهد. بسته به میزان جفت شدن بازها، نتیجه میتواند تخریب mRNA یا کاهش سنتز پروتئین باشد. در برخی موارد، miRNAها همچنین میتوانند ترجمه را فعال یا رونویسی را تنظیم کنند (14). از آنجاییکه عملکرد اصلی miRNAها جلوگیری از بیان ژنهای هدف از طریق مهار ترجمه mRNA است، میتوانند با تاثیر روی بیان ژنهای مختلف در پاتوژنز بیماریهای مختلف از جمله آندومتریوز موثر باشند.
ژنهای زیادی در بیماری آندومتریوز نقش دارند که یکی از آنها ژن TP53 میباشد، که بر روی کروموزوم 17 قرار دارد و حاوی 11 اگزون است و پروتئین سرکوبگر تومور P53 را کد میکند. عملکرد اصلی آن محافظت از یکپارچگی DNA سلول است (15). انواع مختلف استرسهای سلولی منجر به فعالسازی ژن TP53 میشود که این استرسهای سلولی شامل آسیب DNA، فعال شدن انکوژنها، هیپوکسی، استرس اکسیداتیو و تغییرات متابولیک سلولی میباشد (16). فرآیندهای زیستی مهمی که نقش P53 در آنها مشخص شده است، توقف چرخه سلولی، پیری، ترمیم DNA، آپوپتوزیس و پیشگیری از سرطان هستند که در میان آنها، شناخته شده ترین مسیرها توقف چرخه سلولی و آپوپتوز است (17). P53 در القای آپوپتوزیس نقش دارد بهنحوی که اعضای خانواده پرو آپوپتوزیس مانند Bax، Bak، Bid را فعال میکند که منجر به تخریب اعضای خانواده ضد آپوپتوزیس Bcl-2 میشود. (18). جهشهای TP53 میتوانند فرآیندهای سلولی مختلفی از جمله آپوپتوزیس، کنترل چرخه سلولی ترمیم DNA را مختل کرده و با تمایززدایی اپیتلیال مرتبط هستند، فرآیندی که در آن سلولها ویژگیهای تخصصی خود را از دست میدهند که همه آنها با ایجاد و پیشرفت آندومتریوز مرتبط هستند. بهطور خلاصه، جهشهای TP53 و عملکرد تغییر یافته p53 با مختل کردن فرآیندهای سلولی مختلف، تاثیر بر بیان ژن و بهطور بالقوه تاثیرگذاری بر ریزمحیط تومور، در ایجاد و پیشرفت آندومتریوز نقش دارند (19).
یکی دیگر از ژنهای دخیل در این بیماری، FOXO1 است که عضوی از خانواده FOXO میباشد که در استرس اکسیداتیو (20)، آپوپتوزیس سلولی (21)، تنظیم متابولیک و ایمنی (22)، طول عمر و سرطانزایی (23) دخالت دارد FOXO1 نقش مهمی در بازسازی آندومتر و هومئوستاز در طول چرخه قاعدگی و اوایل بارداری دارد. در آندومتر، پروتئینهای FOXO1 بهدلیل توانایی آنها در تنظیم ژنهای پروآپوپتوزیس یا ژنهای دخیل در تمایز، توقف چرخه سلولی، دفاع اکسیداتیو و ترمیم DNA بهعنوان واسطههای کلیدی سرنوشت سلول مطرح شدهاند. FOXO1 در پاییندست مسیر NOTCH عمل میکند و کاهش سیگنالدهی NOTCH در آندومتریوز میتواند منجر به کاهش بیان FOXO1 و اختلال در دسیجوایی شدن گردد. فسفریله شدن FOXO1 در ناحیه سرین 184 باعث مهاجرت سلولهای استرومایی آندومتر شده و در نتیجه منجر به بیماری آندومتریوزیس میشود (24).
بهطور کلی برای ایجاد آندومتریوز فرایندهایی نظیر تکثیر سلولی، مهاجرت و مرگ سلولی و همچنین ایجاد عروق خونی جدید نقش مهمی بازی میکند. لذا مطالعه بیان ژنهایی نظیر TP53 و FOXO1 که در فرایندهای فوق نقش دارند بههمراه بررسی بیان miRNAها که در تنظیم بیان ژنها نقش مهمی بازی میکنند، میتواند به درک مکانیسمهای ایجاد بیماری آندومتریوز کمک کند. با توجه به مطالب فوق و نقش تنظیمی miR-223-3p در بیان ژنها و همچنین در فرایند EMT و از طرفی اهمیت ژنهای TP53 و FOXO1 در فرایندهایی نظیر آپوپتوزیس، تکثیر سلولی و بازسازی آندومتر رحم، لذا در این تحقیق به بررسی بیان miR-223-3p و ژنهای هدف آن یعنی TP53 و FOXO1 در بافت آندومتریوز تخمدانی از طریق تحلیل بیوانفورماتیکی و تایید تجربی پرداخته شده است.
2-مواد و روش
پژوهش حاضر در دو بخش بیوانفورماتیک و تجربی صورت گرفت.
بخش بیوانفورماتیک
در این پژوهش، از اصطلاح " آندومتریوز" بهعنوان یک کلمه کلیدی در پایگاه داده بیان ژنGEO (Gene Expression Omnibus) به آدرس )http://www.ncbi.nlm.nih.gov/geo( استفاده شد و مجموعه داده سری GSE105765 انتخاب شد. سری GSE105765 شامل پروفایل بیانmicroRNA های مختلف در بافتهای آندومتر اکتوپیک و یوتوپیک در بیماران مبتلا به آندومتریوز تخمدانی است. این مجموعه داده بر اساس پلتفرم GPL11154 تهیه شده که مبتنی بر ریز آرایه ژنومی انسانی و بر پایه دادههای توالییابی RNA (RNA-seq) میباشد. نمونهها مربوط به 8 بیمار مبتلا به آندومتریوز تخمدانی بود و از هر بیمار یک نمونه آندومتر اکتوپیک و یک نمونه آندومتریوتوپیک آنالیز شد.
تجزیه و تحلیل ژنهای با بیان افتراقی در دوگروه اکتوپیک و یوتوپیک: فایلهای مربوط به مجموعه داده سری GSE105765 از پایگاه دادهGEO دانلود شد. تجزیه و تحلیل دادههای دانلود شده، بین دو گروه نمونههای اکتوپیک و نمونههای یوتوپیک بیماران مبتلا به آندومتریوز تخمدانی، با استفاده از بستههای آماری limma و voom در محیط RStudio انجام شد. با انجام پردازشهای مربوطه، داده های غیر مرتبط حذف و جدول مربوط به ژنهای با بیان افتراقی بین دو گروه بهدست آمد. در ادامه، نمودارهایی رسم شدند که هر یک به تحلیل و درک بهتر تغییرات بیان ژنها، بین نمونههای اکتوپیک و یوتوپیک کمک کردند. ابتدا نمودار نقشه حرارتی (Heatmap) رسم شد، این نمودار رایجترین روش برای تصویرسازی دادهها بهشمار میرود و میزان بیان هر ژن در هر نمونه به تفکیک گروهها را نشان میدهد. در این نمودار، رنگ قرمز نشاندهنده بیان بالای ژن ها و رنگ آبی نشاندهنده بیان پایین ژنها میباشد. در مرحله بعد نمودار آتشفشانی (Volcano) این ژنها رسم شد. نمودار آتشفشانی نوعی نمودار پراکنده است که اهمیت سطح معنیداری (P-value) را در مقابل بزرگی تغییر (fc=fold change) نشان میدهد. در نهایت بر اساس نتایج این تحلیلها، ژنهایی که تفاوت بیان معنیدار و بیشتری در گروه اکتوپیک نسبت به یوتوپیک داشتند و تاکنون بهصورت گسترده بررسی نشده بودند، شناسایی و از میان آن ها ژن mir-223-3p انتخاب شد.
پیشبینی اهداف miR-223-3p : بهمنظور دست یابی به اهداف پیش بینی شده miR-223-3p در انسان، از پایگاه داده معتبر miRDB به آدرس (https://mirdb.org) استفاده شد. با وارد کردن miRNA مورد نظر و انتخاب جانور مورد مطالعه یعنی انسان، ژنهای مرتبط با miR-223-3p نمایش داده شد و در نهایت، از بین ژنهای هدف miR-223-3p، با بررسیهای صورت گرفته دو ژن FOXO1 و TP53 انتخاب شدند. در مرحله بعدی با بهرهگیری از نرمافزار Cytoscape (نسخه 3.9.1) شبکه ای از ژنهای مرتبط با miR-223-3p، سازماندهی و ترسیم شدند. این شبکه نشان دهنده تعاملات پیشبینی شده بین miR-223-3p و ژنهای هدف مختلف آن میباشد.
بخش آزمایشگاهی
روش جمع آوری نمونهها: جهت تایید یافتههای بیوانفورماتیکی در خصوص تغییر بیان miR-223-3p در نمونههای آندومتریوزی، نمونههای بافتی اکتوپیک و یوتوپیک از ۴۰ بیمار مبتلا به آندومتریوز تخمدانی در بیمارستان آریا شهر رشت جمعآوری شد. گروه کنترل نیز شامل ۴۰ بیمار بدون آندومتریوز بوده که بهدلیل لاپاراسکوپی تشخیصی تحت جراحی قرار گرفتند. میانگین سنی افراد بیمار در این پژوهش 4±5/31 سال و برای گروه کنترل 4±4/30 سال بود. معیارهای انتخاب شامل: تشخیص بیماری از طریق لاپاراسکوپی، عدم مصرف دارو طی سه ماه پیش از لاپاراسکوپی و نمونهبرداری و ارائه رضایتنامه کتبی بود. در مقابل، معیارهای عدم انتخاب شامل عدم ارائه رضایتنامه کتبی، عدم ابتلا به انواع تومورها، دیابت، بیماریهای مزمن کبدی، کلیوی و بیماری های روماتولوژی نظیر آرتریت روماتوئید در نظرگرفته شد. تشخیص بیماری توسط پزشک متخصص زنان و زایمان بیمارستان آریای شهر رشت در طی لاپاروسکوپی بیماران صورت گرفت و کلیه نمونههای دریافت شده از اتاق عمل، بلافاصله توسط ازت مایع فریز شده و سریعا به فریزر 70- درجه سانتیگراد در آزمایشگاه تحقیقاتی گروه زیست شناسی دانشگاه گیلان منتقل شدند. این مطالعه توسط کمیته اخلاق دانشگاه گیلان با شناسه اخلاق IR.GUILAN.REC.1403.006 تایید شد.
استخراج RNA کل از نمونه های بافتی: بهمنظور استخراج RNA کل از نمونههای بافتی، از معرف ترایزول استفاده شد و مطابق با دستورالمعل مربوط به شرکت سازنده ترایزول، استخراج RNA کل صورت گرفت. برای اطمینان از استخراج صحیح RNA از الکتروفورز افقی در حضور ژل 5/1 درصد استفاده شد و سپس سطح کمّی غلظت RNA، بهکمک دستگاه نانودراپ مورد بررسی قرار گرفت. برای اطمینان از وجود DNA، تمام نمونهها تحت تیمار با DNase محصول شرکت (Thermo Scientific, USA) قرار گرفتند.
طراحی آغازگر و سنتز cDNA: با توجه به ناپایداری و امکان تخریب RNA، لازم است تا به مولکول پایدار cDNA تبدیل شود. بدینمنظور ابتدا طراحی آغازگرهای اختصاصی جهت سنتز cDNA توسط نرم افزار Oligo7 انجام شد (جدول 1).
جدول1: ویژگی آغازگرهای مورد استفاده.
|
طول محصول(جفت باز) |
توالی |
اغازگر |
|
63 bp |
GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACCCCAGCTA 5’- ACGCAGAAGGAGCTTACAATCTA-3’ 5’- GTGCAGGGTCCGAGGT-3’ |
hsa-miR-223-3p-STEM
hsa-miR-223-3p-R hsa-miR-223-3p-F |
|
122 bp |
5’- TCAAGAGCGTGCCCTACTTC-3’ 5’- AGTTCCTTCATTCTGCACACGA-3’ |
FOXO1-F FOXO1-R |
|
120 bp |
5’- TCCTCAGCATCTTATCCGAGT-3’ 5’- GTACAGTCAGAGCCAACCTCA-3’ |
TP53-F TP53-R |
|
235 bp |
5’- CATCACCATCTTCCAGGAGCG-3’ 5’- GGAGGCATTGCTGATGATCTTG-3’ |
GAPDH-F GAPDH-R |
|
94 bp |
5’- CTCGCTTCGGCAGCACA-3’ 5’- AACGCTTCACGAATTTGCGT-3’ |
U6-F U6-R |
ویژگیهای پرایمر نظیر طول، درصد GC و دمای ذوب بررسی شدند. برای اطمینان از اختصاصی بودن پرایمرهای طراحی شده، از پایگاه Blast (www.blast.ncbi.nlm.nih.gov/Blast.cgi) استفاده شد. نتایج این تحلیل، تاییدکننده اختصاصی بودن محل اتصال پرایمرهای طراحی شده بود.
سنتز cDNA با استفاده از کیت سنتز cDNA شرکت زیست ویرایش به شماره کاتالوگ ZV06010050 ، ایران و بر اساس دستورالعمل پیشنهادی شرکت سازنده، صورت گرفت. مراحل انکوباسیون شامل 60 دقیقه در دمای 55 درجه سانتیگراد، 10 دقیقه در دمای 95 درجه سانتیگراد و در نهایت نگهداری در دمای 4 درجه سانتیگراد بود. پس از اتمام مراحل، میکروتیوپهای حاوی cDNA برای نگهداری طولانی مدت به فریزر با دمای منفی 20 درجه سانتیگراد منتقل شدند.
بررسی بیان ژن توسط واکنش زنجیرهای پلیمراز در زمان واقعی (Real Time PCR): برای ارزیابی بیان ژنهای miR-223-3p، FOXO1 و TP53 از روش Real Time PCR با استفاده از رنگ SYBR Green (یکتا تجهیز، ایران، شماره کاتالوگ: YT2551) استفاده شد. مخلوط واکنش Real Time PCRبا توجه به پروتکل پیشنهاد شده توسط شرکت سازنده SYBR Green تهیه شد و به استریپهای مخصوص منتقل شد. بهمنظور اطمینان از صحت نتایج، هر نمونه در دو تکرار بهدرون استریپ منتقل شد و در نهایت استریپهای حاوی نمونه به دستگاه Real Time PCR (Roche, Germany) منتقل و واکنش زنجیرهای پلی مراز مطابق شرایط دمایی – زمانی جدول 2 انجام شد. بر اساس نتایج بهدستآمده از Gradient PCR، دمای بهینه اتصال برای تمامی ژنهای مورد مطالعه 60 درجه سانتیگراد تعیین شد. همچنین بهمنظور مقایسه و بررسی درستی بیان ژنهای هدف، از GAPDH بهعنوان ژن مرجع برای ژنهای FOXO1، TP53 و ژن مرجع U6 برای miR-223-3p استفاده شد.
جدول2: برنامه دمایی-زمانی واکنش Real Time PCR
|
تعداد چرخه |
دمای (درجه سانتیگراد) |
زمان (دقیقه) |
مرحله |
|
1 |
95°C |
5 دقیقه |
واسرشتسازی اولیه |
|
40 |
95°C |
30 ثانیه |
واسرشتسازی ثانویه |
|
60°C |
35 ثانیه |
اتصال آغازگرها به رشته الگو |
|
|
72°C |
40 ثانیه |
طویل شدن |
|
|
1 |
72°C |
5 دقیقه |
بسط نهایی |
|
- |
4°C |
نامحدود |
نگهداری |
3-آنالیز آماری
پس از اتمام فرآیند Real-time PCR جهت بررسی بیان نسبی ژنهای GAPDH، FOXO1، TP53، U6 و miR-223-3p، تمام Ct ها در فایل اکسل جمعآوری شده و طبق فرمولهای زیر مقادیر ΔCt،ΔΔCt و Fold change محاسبه شد. سپس آنالیز آماری جهت بررسی میزان تغییرات بیان ژنهای مربوطه، بررسی معنیدار بودن میزان تغییرات بیان و همچنین نمودارها و گرافها توسط نرمافزار GraphPad Prism نسخه 8 رسم شد.
فرمولها:
گروه بیمارΔCt = Ct (ژن مورد نظر) - Ct (ژن کنترل داخلی)
گروه کنترلΔCt = Ct (ژن مورد نظر) - Ct (ژن کنترل داخلی)
ΔΔCt = ΔCt - ΔCt (میانگین گروه کنترل)
بخش بیو انفورماتیک
مقایسه بیان miR-223-3p در دو گروه اکتوپیک و یوتوپیک
در مجموعه داده سریGSE105765، وضعیت بیان miRNA های مختلف در بافت آندومتر یوتوپیک با بافت اکتوپیک در بیماران مبتلا به آندومتریوز تخمدانی مقایسه شدند. تجزیه و تحلیل ژنهای با بیان افتراقی از طریق بستههای مختلف در نرم افزار Rstudiu انجام و نقشه حرارتی این ژنها رسم شد (شکل 1).
شکل1: نقشه حرارتی مربوط به سطوح بیان miRNA های با بیان متفاوت در نمونههای بافت اکتوپیک (EC) در مقایسه با نمونههای بافت یوتوپیک (EU) در بیماران مبتلا به آندومتریوز تخمدانی. رنگ قرمز بیانکننده افزایش بیان و رنگ آبی نشاندهنده کاهش بیان miRNA ها است. از میان miRNA های دارای افزایش بیان، miR-223-3p انتخاب شد.
مطابق شکل 1 برخی از miRNA ها در نمونههای اکتوپیک نسبت به نمونههای یوتوپیک کاهش بیان و برخی از miRNA ها افزایش بیان داشتند. سپس miRNA هایی که در بافت اکتوپیک نسبت به یوتوپیک افزایش بیان داشتند، مورد بررسی قرار گرفتند و در نهایت mir-223-3p که تغییرات بیان آن در دو گروه اکتوپیک و یوتوپیک معنیدار بود، انتخاب شد. سپس نمودار آتشفشانی (volcano) که نشاندهنده معنیدار بودن و تفاوت بیان miRNA های مختلف میباشد، رسم و miRNA منتخب نیز در آن مشخص شد (شکل 2).
شکل2: نمودار آتشفشانی مربوط به .mir-223-3pهر نقطه نشاندهنده یک ژن است. نقاط قرمز نشاندهنده ژنهای با افزایش بیان، نقاط آبی نمایانگر ژنهای با کاهش بیان و نقاط خاکستری ژنهای بدون تغییر بیان را نشان میدهند.
تعیین ژنهای هدف و رسم شبکه miR-223-3p : برای یافتن اهداف پیش بینی شده miR-223-3p از پایگاه داده miRDB استفاده شد. پس از شناسایی ژنهای هدف پیش بینی شده miR-223-3p از این پایگاه داده، با استفاده از نرم افزار Cytoscape (نسخه 3.9.1) برای miR-223-3p شبکهای طراحی شد که ارتباط این miRNA با ژنهای هدف را نشان میدهد. در نهایت از بین ژنهای هدف miR-223-3p، دو ژن FOXO1 و TP53 که در بیماری آندومتریوز نیز نقش داشتند، انتخاب شدند و در شبکه ترسیم شده نیز نمایش داده شدند (شکل 3).
شکل 3: شبکه ژنهای هدف .miR-223-3p شبکه فوق ژنهای هدفmiR-223-3p را نشان میدهد که از پایگاه داده miRDB استخراج شدهاند. miR-223-3p با رنگ آبی مشخص شده است. هر گره نمایانگر یک ژن هدف و رنگ آن نشاندهنده امتیاز پیشبینی شده تعامل است، نقاط پررنگتر بیانگر امتیازات بالاتر در تعامل هستند. از میان ژن های هدف miR-223-3p، دو ژن FOXO1 و TP53 انتخاب شدند.
نتایج بخش آزمایشگاهی
نتایج حاصل از بررسی کمی و کیفی RNA استخراج شده
بررسی نمونهها با دستگاه نانودراپ نشان داد که غلظت RNA در محدوده 600 تا 900 نانوگرم بر میکرولیتر متغیر بود. همچنین نسبت جذب نوری در طول موج های 260 و 280 نانومتر، در محدوده مطلوب 1.8 تا 2 قرار داشت. برای ارزیابی دقیقتر کیفیت RNA های استخراج شده، از الکتروفورز ژل آگارز استفاده شد. مشاهده باندهای مشخص 18s rRNA،28s rRNA و 5.8s rRNA تاییدکننده کیفیت مناسب RNA های استخراج شده بود (شکل 4).
شکل 4: تصویر ژل آگارز 5/1 درصد از RNA کل استخراجشده. ستون اول از سمت چپ نشاندهنده مارکر 100 bp و ستونهای سمت راست نشاندهنده RNA کل استخراجشده هستند. مشاهده سه باند نشان دهنده کیفیت خوب RNA استخراج شده می باشد.M : مارکر100 bp.
نتایج بیان miR-223-3pدر نمونههای اکتوپیک، یوتوپیک و گروه کنترل
جهت بررسی تغییرات بیان ژن miR-223-3p در نمونههای بافت آندومتریوم اکتوپیک و یوتوپیک افراد مبتلا به آندومتریوز تخمدانی و بافت آندومتر طبیعی افراد سالم، Ct های بهدست آمده برای هر نمونه به نرمافزار Excel منتقل و سپس میزان بیان و Fold Change نمونهها محاسبه شد. میانگین بیان نسبی ژن miR-223-3p در گروه اکتوپیک برابر با 72/0±13/4 در گروه یوتوپیک 10/0±10/1 و در گروه کنترل برابر با 09/0±00/1 محاسبه شد. بهعبارت دیگر، بیان نسبی ژن miR-223-3p در گروه اکتوپیک نسبت به گروه یوتوپیک و گروه کنترل بهطور معنیداری افزایش یافته است (0001/0< p(. درحالیکه بین دو گروه کنترل و یوتوپیک تفاوت معنیداری مشاهده نشد (54/0= p) (شکل 5).
شکل5: میزان بیان نسبی miR-223-3p در سه گروه بافت کنترل(آندومتر طبیعی) و بافت اکتوپیک و یوتوپیک بیماران مبتلا به آندومتریوز تخمدانی. Control: گروه کنترل،:Eutopic بافت یوتوپیک، :Ectopic بافت نابهجا(اکتوپیک)، **** نشاندهنده تفاوت معنیدار بین گروهها و ns نشان دهنده عدم تفاوت معنیدار بین گروهها است. میزان بیان miR-223-3p در گروه اکتوپیک، 72/0 ± 138/4 در گروه یوتوپیک، 10/0 ± 100/1 و در گروه کنترل، 09/0±000/1 است. افزایش معنیداری بین گروه اکتوپیک نسبت به گروه یوتوپیک و گروه کنترل دیده شد (0001/0. درحالیکه بین دو گروه کنترل و یوتوپیک تفاوت معنیداری مشاهده نشد (54/0= p).
نتایج بیان FOXO1 در نمونههای اکتوپیک و یوتوپیک و گروه کنترل
در این مطالعه، میزان بیان نسبی ژن FOXO1 در نمونههای بافت اکتوپیک، یوتوپیک بیماران مبتلا به آندومتریوز تخمدانی و گروه کنترل بررسی شد. نتایج نشان داد که میانگین بیان FOXO1 در گروه اکتوپیک (12/0±47/0) بهطور قابل توجهی پایین تر از گروه یوتوپیک (06/0± 02/1) و گروه کنترل (11/0±001/1) بود. تحلیل واریانس یکطرفه (One-way ANOVA) کاهش معنی داری در بیان FOXO1 در گروه اکتوپیک نسبت به گروه یوتوپیک و همچنین نسبت به گروه کنترل نشان داد (0001/0. از سوی دیگر، مقایسه بین گروه یوتوپیک و گروه کنترل تفاوت معنیداری در میزان بیان FOXO1 نشان نداد (53/0= p) (شکل 6).
شکل6: میزان بیان نسبی FOXO1 در سه گروه کنترل (آندومتر طبیعی) و بافت اکتوپیک و یوتوپیک بیماران مبتلا به آندومتریوز تخمدانی. Control: گروه کنترل، :Eutopic بافت یوتوپیک، :Ectopic بافت نابهجا(اکتوپیک)، **** نشان دهنده تفاوت معنیدار بین گروهها و ns نشان دهنده عدم تفاوت معنیدار بین گروهها است. میزان بیان FOXO1 در گروه اکتوپیک، 12/0 ± 475/0 در گروه یوتوپیک، 06/0±026/1 و در گروه کنترل، 11/0±001/1 است. کاهش معنیداری بین گروه اکتوپیک نسبت به گروه یوتوپیک و گروه کنترل دیده شد (0001/0
. درحالیکه بین دو گروه کنترل و یوتوپیک تفاوت معنیداری مشاهده نشد (53/0= p).
نتایج بیان TP53 در نمونههای اکتوپیک و یوتوپیک و گروه کنترل
میزان بیان نسبی ژن TP53 در گروه اکتوپیک (14/0±335/0)، گروه یوتوپیک (09/0±059/1) و گروه کنترل (12/0± 018/1) مورد بررسی قرار گرفت. تحلیل آماری نشان داد که بیان ژن TP53 در گروه اکتوپیک نسبت به گروه یوتوپیک و گروه کنترل بهطور معنیداری کاهش یافته است (0001/0 . از سوی دیگر، تفاوت معنیداری بین بیان TP53 در گروه یوتوپیک و گروه کنترل مشاهده نشد (28/0= p) (شکل7).
شکل7: میزان بیان نسبی TP53 در سه گروه کنترل (آندومتر طبیعی) و بافت اکتوپیک و یوتوپیک بیماران مبتلا به آندومتریوز تخمدانی. Control: گروه کنترل، :Eutopic بافت یوتوپیک، :Ectopic بافت نابهجا(اکتوپیک)، **** نشان دهنده تفاوت معنیدار بین گروهها و ns نشان دهنده عدم تفاوت معنیدار بین گروهها است. میزان بیان TP53 در گروه اکتوپیک، 14/0 ± 335/0 در گروه یوتوپیک، 09/0 ± 059/1 و در گروه کنترل، 12/0 ± 018/1 است. کاهش معنیداری بین گروه اکتوپیک نسبت به گروه یوتوپیک و گروه کنترل دیده شد (0001/0< p(. درحالیکه بین دو گروه کنترل و یوتوپیک تفاوت معنیداری مشاهده نشد (28/0= p (.
آندومتریوز یکی از بیماریهای شایع زنان است که حدود ۱۰ درصد زنان در سنین باروری را در سراسر جهان تحت تاثیر قرار میدهد (۲۵). اگرچه آندومتریوز بهطور کلی یک وضعیت خوشخیم در نظر گرفته میشود، اما خطر تبدیل شدن به بیماری سرطان، از جمله سرطان تخمدان را دارد. شباهتهای فنوتیپی میان سلولهای آندومتریوزی و سلولهای سرطانی، نظیر رشد بیرویه سلولی و کاهش فرآیند آپوپتوزیس، نشاندهنده ویژگیهای مولکولی مشترک آنها است (۲۶). miRNA ها گروهی از RNA های کوچک غیرکدکننده هستند که در تنظیم بیان ژنها نقش کلیدی ایفا میکنند. این مولکولها با اتصال به ناحیه ۳’- UTR در mRNA های هدف خود باعث کاهش ترجمه یا تخریب mRNA میشوند (۲۷) و بدینترتیب در فرآیندهای سلولی مانند تکثیر، تمایز، مرگ برنامهریزی شده سلول و پاسخ به استرسهای سلولی نقش دارند (۲۸). تغییر در بیان miRNA ها با بسیاری از بیماریها از جمله سرطانها (۲۹) و بیماریهای التهابی مانند آندومتریوز مرتبط است (30 _ miR-223-3p) یکی از miRNA های کلیدی است که با هدف قرار دادن ژنهای مختلف، در فرایندهای التهابی، تمایز سلولهای ایمنی، تنظیم تکثیر سلولی و مرگ سلولی نقش دارد (۳۱). مطالعات متعددی در خصوص بیان miR-223-3p در بیماریهای مختلف انجام شده است:
در مطالعهای Wang و همکاران (۱۳)، نشان دادند که بیان miR-223-3p در بافتهای سرطان پستان بهطور چشمگیری افزایش مییابد و یک رابطه فیزیکی و عملکردی بین miR-223-3p و (F-box and WD repeat domain-containing 7) FBXW7 وجود دارد و miR-223-3p با تنظیم منفی بیان FBXW7، بهعنوان یک فاکتور محرک تومور عمل کرده و تهاجم و متاستاز سلولهای سرطان پستان را تقویت میکند. Du و همکاران (۳۲) نشان دادند که miR-223 از طریق فعالسازی مسیر سیگنالدهی Hippo/Yap1 (Yes-associated protein 1) در تکثیر، مهاجرت، تهاجم و EMT سلولهای سرطان پستان نقش دارد. در مطالعهی Fang و همکاران (۳۳) نشان دادند که بیان miR-223 در ردههای سلولی سرطان تخمدان شامل SKOV3، OVCAR3، A2780 و ES2 در مقایسه با رده سلولی تخمدانی غیرسرطانی تغییر نشان میدهد. بیان نابجای miR-223 در این ردههای سلولی، سطح پروتئینهای Ki67 (نشانگر تکثیر سلولی) و PCNA (proliferating cell nuclear antigen) را افزایش میدهد. آنها همچنین نشان دادند که مهار اختصاصی miR-223 در سلولهای سرطان تخمدان، از طریق حذف مهار مستقیم آن بر SRY-box transcription factor 11 (SOX11) منجر به افزایش بیان این پروتئین شده و مهاجرت و تهاجم سلولی را کاهش میدهد.Lu و همکاران (۳۴) نشان دادند که در سرطان اپیتلیال تخمدان بیان circBNC2 (RNA حلقوی BNC2) و LARP4 (پروتئین مرتبط با La-4) کاهش یافته، در حالیکه بیان miR-223-3p افزایش یافته است. circBNC2 بهعنوان یک اسفنج مولکولی برای miR-223-3p عمل میکند و آن را مهار میکند، LARP4 نیز هدف مستقیم miR-223-3p است. بهطورکلی، کاهش بیان circBNC2 و LARP4 همزمان با افزایش بیان miR-223، رشد سرطان اپیتلیال تخمدان را تقویت میکند، بنابراین محور circBNC2/miR-223/LARP4 میتواند بهعنوان یک مکانیسم تنظیمی هدفمند برای درمان سرطان اپیتلیال تخمدان مطرح شود. در طی مطالعاتی که Liu و همکاران (۳۵) انجام دادند، نشان داده شد که miR-223 بهطور مستقیم رونوشت FBXW7 را مورد هدف قرار میدهد، که زیر واحد شناسایی سوبسترا در (SCF (SKP1/CUL1/F-box protein است، بنابراین آپوپتوزیس را سرکوب کرده و از طریق مسیرهای Notch و Akt/mTOR میزان تکثیر سلولی را تغییر میدهد. Cheng و همکاران (۳۶) نشان دادند که miR-223 با دخالت در مکانیسمهای مرگ سلولی در سلولهای سرطان روده بزرگ و افزایش تکثیر سلولی از طریق کاهش بیان BCL2-like 11 (BCL2L11/BIM) از طریق fork head box O3 (FOXO3a)، موجب افزایش تکثیر سلولی میشود. همچنین Ju و همکاران (۳۷) در مطالعهای نشان دادند که در رده سلولی SW620 مربوط به سرطان کولورکتال، miR-223 سیگنالدهی FoxO3a/BIM را بهطور منفی تنظیم میکند و بدینترتیب تکثیر سلولی افزایش مییابد.
FOXO1 و TP53 از مهمترین اهداف miR-223-3p هستند که نقشهای حیاتی در تنظیم فعالیتهای سلولی دارند. Wang و همکاران (38)، نشان دادند که بیان NEK2 در آندومتر اکتوپیک بیماران مبتلا به آندومتریوز افزایش یافته است. در عین حال، NEK2 با FOXO1 تعامل دارد و FOXO1 را در موقعیت Ser184 فسفریله میکند، که باعث کاهش پایداری پروتئین FOXO1 میشود. فسفریلاسیون FOXO1 توسط NEK2 موجب افزایش تکثیر سلولی، مهاجرت، تهاجم و اختلال در واکنش دسیدوایی میشود. همچنین Li و همکاران (39)، نشان دادند که m6A (N6-methyladenosine) تنظیمشده توسط (methyltransferase-like 3) METTL3 باعث تخریب mRNA FOXO1 از طریق YTHDF2 میشود و به اینترتیب بر واکنش دسیدوایی سلولی، کاشت جنین و قدرت باروری زنان مبتلا به آندومتریوز تاثیر میگذارد. Shan و همکاران (40)، در طی مطالعهای نشان دادند که سطوح بیان Rab27b، در ضایعات آندومتریوزی بهطور قابل توجهی بالاتر از آندومتر طبیعی میباشد واز طرفی استروژن سطوح FOXO1 فسفریله شده را در ضایعات آندومتریوزی افزایش میدهد و باعث افزایش بیان Rab27b میشود. Rab27b یکی از اعضای خانواده GTPase های Rab است و ارتباط قوی با رشد، تهاجم و متاستاز انواع مختلف تومورها دارد. در مطالعهی Rezk و همکاران (41) دریافتند که miR-34a ممکن است از طریق SIRT-1 بیان FOXO1 را در بافت آندومتریوزی تنظیم کند و به این ترتیب میتواند با کاهش آپوپتوزیس طبیعی در آندومتر به پاتوژنز آندومتریوز کمک کند. ژن TP53 پروتئین p53 را کد میکند که میتواند بهطور موثر باعث توقف چرخه سلولی، ترمیم DNA، پیری سلولی، آپوپتوزیس و اتوفاژی شود. جهشهای TP53 تقریبا در 50 درصد از تومورها مشاهده میشود (42). Ghafouri- Fard و همکاران (43)، کاهش غیرطبیعی بیان TP53 و lncRNA های مرتبط با آن را در بافتهای سرطان پستان نشان دادند که با نقش حیاتی این ژن سرکوبکننده تومور در پیشگیری از سرطان همخوانی دارد. همچنین Guo و همکاران (44) نشان دادند که TRIM31 بهعنوان یک E3 ubiquitin ligase میتواند با تثبیت و جلوگیری ازتخریب p53 از پیشرفت سرطان پستان جلوگیری کند و شاید در آینده بهعنوان یک روش درمانی جدید برای درمان این نوع سرطان مورد استفاده قرار بگیرد. در مطالعهای که توسط Guo و همکاران (45) انجام شد، مشخص گردید که بیان LEMD1-AS1 (RNA آنتیسنس LEMD1) در بافتهای سرطان تخمدان و ردههای سلولی آن کاهش یافته است و نتایج این تحقیق نشان داد که LEMD1-AS1 میتواند مستقیما با miR-183-5p و پروتئین p53 تعامل داشته باشد. نقش ضدتوموری LEMD1-AS1 در پیشرفت سرطان تخمدان، وابسته به بیان TP53 توسط miR-183-5p است. Duan و همکاران (46)، به بررسی بیان و اهمیت بالینی کیناز دی فسفات نوکلئوتیدآ (Nm23-H1)، p53 و اینتگرین β1 در آندومتریوز، بافتهای آندومتر نرمال و اکتوپیک پرداختند و دریافتند که در مقایسه با بافت آندومتر نرمال، سطح پروتئینهای nm23-H1 و p53 به طور قابل توجهی در بافتهای آندومتر اکتوپیک کاهش یافته، در حالیکه پروتئین اینتگرین β1 افزایش یافته است. بنابراینNm23-H1 و p53 ممکن است پیشرفت آندومتریوز را مهار کنند، در حالیکه اینتگرین β1 اثر ترویجی دارد. Ma و همکاران (47)، نشان دادند که ارتباط معنیداری بین پلی مورفیسم ژن p53 و میزان شیوع آندومتریوز در زنان وجود دارد. بهطوری که کاهش فراوانی ژن GG و افزایش آلل C ممکن است خطر ابتلا به این بیماریها را افزایش دهد. نتایج این مطالعه نشان داد که بیان miR-223-3p در نمونههای اکتوپیک بیماران مبتلا به آندومتریوز تخمدانی نسبت به نمونههای یوتوپیک و کنترل بهطور معنیداری افزایش یافته است، در حالیکه بیان ژنهای FOXO1 و TP53 در نمونههای اکتوپیک بیماران مبتلا به آندومتریوز تخمدانی نسبت به نمونههای یوتوپیک و کنترل بهطور معنیداری کاهش یافته است. از محدودیتهای این مطالعه میتوان به موارد زیر اشاره کرد:
این تحقیق منحصرا روی miR-223-3p و دو ژن هدف آن یعنی FOXO1 وTP53 متمرکز شد. یک miRNA منفرد میتواند چندین ژن را هدف قرار دهد، در حالیکه یک ژن میتواند تحت تاثیر RNA های تنظیمی مختلف مانند microRNA ها و RNAهای حلقوی (circRNAs) قرار گیرد. بنابراین، سایر RNA های طولانی غیر کد کننده (lncRNA) و ژنهای هدف آنها نیز ممکن است در ایجاد آندومتریوز نقش داشته باشند. این مطالعه نشان داد که سطح بیان miR-223-3p در نمونههای بافتی از بیماران مبتلا به آندومتریوز تخمدان بیشتر از گروه کنترل بود. برعکس در بیمارانی که با آندومتریوز تخمدان تشخیص داده شدند، سطوح FOXO1 و TP53 بهطور قابل توجهی کاهش نشان داد. کاهش سطوح FOXO1 و TP53 ممکن است نتیجه افزایش بیان miR-223-3p باشد، زیرا هر دو هدف این miRNA خاص هستند، همانطور که توسط تجزیه و تحلیل بیوانفورماتیک پیشنهاد شده است. لذا بررسی تعداد بیشتری از miRNA ها و ژنهای دخیل در این بیماری میتواند به نتیجهگیری صحیحتر کمک کند.
نتیجهگیری میشود که بیان ژنهای miR-223-3p، FOXO1 و TP53 در بافت آندومتر بیماران مبتلا به آندومتریوز تغییر میکند. همچنین، این نتایج نشان میدهد که این ژنها ممکن است بهعنوان عوامل بالقوه در پاتوژنز (بیماریزایی) آندومتریوز تخمدانی مطرح شوند و دارای پتانسیل درمانی برای این بیماری باشند.
7-تشکر و قدردانی
بدین وسیله از معاونت پژوهشی دانشکده علوم پایه دانشگاه گیلان رشت جهت حمایت و تامین منابع مالی این پژوهش تشکر و قدردانی میشود.
-