دوره 34، شماره 5 - ( مرداد 1405 )                   جلد 34 شماره 5 صفحات 10234-10223 | برگشت به فهرست نسخه ها

Ethics code: IR.IAU.MSHD.REC.1402.077
Clinical trials code: ندارم


XML English Abstract Print


Download citation:
BibTeX | RIS | EndNote | Medlars | ProCite | Reference Manager | RefWorks
Send citation to:

bababashi M, Baharara2 J, صالحی م, salehi M. Determining the Molecular Mechanisms Regulated by RAD51 and MRE11 Genes in Oocyte Maturation Using Bioinformatics Studies. JSSU 2026; 34 (5) :10223-10234
URL: http://jssu.ssu.ac.ir/article-1-6644-fa.html
باباباشی مرجان، بهار ارا جواد، Shahroshabadi1 Khadije، صالحی محمد. تعیین مکانیسم های مولکولی تنظیم شده توسط ژن های RAD51 و MRE11 در بلوغ تخمک با استفاده از مطالعات بیوانفورماتیکی. مجله علمي پژوهشي دانشگاه علوم پزشكي شهید صدوقی يزد. 1405; 34 (5) :10223-10234

URL: http://jssu.ssu.ac.ir/article-1-6644-fa.html


متن کامل [PDF 766 kb]   (28 دریافت)     |   چکیده (HTML)  (19 مشاهده)
متن کامل:   (9 مشاهده)
Original Article

Determining the Molecular Mechanisms Regulated by RAD51 and MRE11 Genes in Oocyte Maturation Using Bioinformatics Studies
Marjan Bababashi1 , Javad Baharara1       , Khadije Shahroshabadi1 
, Mohammad Salehi2  

1 Department of Biology, Mashahad Branch, Islamic Azad University, Mashahad, Iran
2 Cellular and Molecular Biology Research Center, Shahid Beheshti University of Medical Sciences, Tehran, Iran

 Corresponding author:
Cellular and Molecular Biology Research Center, Shahid Beheshti University of Medical Sciences, Tehran, Iran

Email: msalehi78@gmail.com
Tel: +989003372897
Postal Code: 8915887857

Received: February 16, 2026
Revised: May 05, 2026
Accepted: May 13, 2026

Citation:
Bababashi M, Baharara J, Shahroshabadi KH, Salehi M. Determining the Molecular Mechanisms Regulated by RAD51 and MRE11 Genes in Oocyte Maturation Using Bioinformatics Studies. J Shahid Sadoughi Uni Med Sci 2026; 34(5): 10223-234.
    Introduction: Female infertility presents a major challenge within the field of reproductive health, since it is influenced by multiple genetic and molecular factors. This study was designed to investigate the role of RAD51 and MRE11 genes in oocyte maturation and their involvement in molecular networks.
Methods: The protein–protein interaction network was retrieved from the STRING database and subsequently analyzed with the help of Cytoscape. The MCODE algorithm was then applied to identify functional clusters, while pathway enrichment tools were used to determine relevant biological pathways. Regulatory miRNAs that may target these genes were examined through the miRTarBase database.
Results: Network analysis revealed that BRCA1 and BARD1 function as central hub genes, whereas ATM and ATRX were identified to be present within the key functional clusters identified by MCODE. Enriched pathways indicated that these genes are involved in homologous recombination, the DNA damage response, cellular senescence, and the regulation of the cell cycle. Furthermore, the miRTarBase analysis showed that several miRNAs, including miR-155 and miR-34a, are capable of targeting RAD51 and potentially regulating its expression.
Conclusion: These findings collectively highlight the critical role of DNA repair–related gene networks in maintaining oocyte quality and suggest that RAD51, MRE11, and their associated pathways could serve as potential therapeutic targets for female infertility.

Keywords: DNA repair; Oocyte maturation; RAD51 and MRE11; Bioinformatics analysis.



مقاله پژوهشی
تعیین مکانیسم های مولکولی تنظیم شده توسط ژن‌های RAD51 و MRE11 در بلوغ تخمک با استفاده از مطالعات بیوانفورماتیکی
مرجان باباباشی 1   ، جواد بهارآرا 1   ، خدیجه شاهرخ آبادی1  ، محمد صالحی2  

1 گروه زیست‌شناسی، واحد مشهد، دانشگاه آزاد اسلامی، مشهد، ایران
2 مرکز تحقیقات زیست شناسی سلولی و مولکولی ، دانشگاه علوم پزشکی شهید بهشتی، تهران، ایران

 نویسنده مسئول:
محمد صالحی 
مرکز تحقیقات قلب و عروق، دانشگاه علوم پزشکی شهید صدوقی، یزد، ایران

پست الکترونیک: msalehi78@gmail.com
تلفن: 09003372897
صندوق پستی: 8915887857

تاریخ دریافت: ............
تاریخ بازنگری: ............
تاریخ پذیرش: ............

ارجاع: 
مرجان باباباشی، جواد بهارآرا، خدیجه شاهرخ آبادی، محمد صالحی. تعیین مکانیسم های مولکولی تنظیم شده توسط ژن‌های RAD51 و MRE11 در بلوغ تخمک با استفاده از مطالعات بیوانفورماتیکی. مجله علمی پژوهشی دانشگاه علوم پزشکی شهید صدوقی یزد 1405؛ 34 (5): 234-10223.
    مقدمه: ناباروری زنان یکی از چالش‌های مهم سلامت باروری است و عوامل ژنتیکی و مولکولی متعددی در بروز آن نقش دارند. این مطالعه با هدف بررسی نقش ژن‌های RAD51 و MRE11 در بلوغ تخمک انجام گرفت.
روش بررسی: در یک مطالعه بیو انفورماتیکی شبکه تعاملی پروتئین–پروتئین از پایگاه STRING استخراج و تحلیل‌های شبکه‌ای در نرم‌افزار Cytoscape انجام گردید. برای شناسایی کلاسترهای عملکردی از الگوریتم MCODE و برای تحلیل مسیرهای زیستی از ابزارهای غنی‌سازی مسیر استفاده شد. بررسی miRNAهای تنظیم‌کننده نیز با استفاده از پایگاه miRTarBase صورت گرفت.
نتایج: تحلیل شبکه نشان داد ژن‌های BRCA1 و BARD1 به‌عنوان گره‌های مرکزی عمل می‌کنند و براساس الگوریتم MCODE ژن‌های ATM و ATRX در کلاسترهای عملکردی مهم شناسایی شدند. مسیرهای زیستی به‌دست‌آمده نشان‌دهنده مشارکت این ژن‌ها در مسیرهایی همچون Homologous Recombination، پاسخ به آسیب DNA، senescence سلولی و چرخه سلولی بود. همچنین نتایج miRTarBase نشان داد چندین miRNA از جمله miR-155 و miR-34a قادر به هدف‌گیری ژن RAD51 هستند. 
نتیجه‌گیری: یافته‌ها اهمیت تنظیم شبکه‌های ژنی مرتبط با ترمیم DNA را در حفظ کیفیت اووسیت نشان داده و بیان می‌کنند که ژن‌های مورد بررسی می‌توانند اهداف بالقوه‌ای برای مداخلات درمانی در ناباروری زنان باشند.

واژه‌های کلیدی: ترمیم  DNA، بلوغ اووسیت، RAD51  و MRE11، تحلیل بیوانفورماتیکی

 
مقدمه
ناباروری در زنان به ناتوانی در باردار شدن و حفظ بارداری تا پایان دوره طبیعی آن، علی‌رغم داشتن روابط جنسی منظم و محافظت‌نشده، اطلاق می‌شود. این مسئله یکی از چالش‌های مهم جهانی است که حدود 5/17٪ از جمعیت بزرگسال، معادل ۱ نفر از هر ۶ نفر را تحت تأثیر قرار می‌دهد و نیاز فوری به مراقبت‌های باروری در دسترس و مقرون‌به‌صرفه را در سراسر جهان نشان می‌دهد (1). سلول‌های انسانی به‌طور مداوم در معرض عوامل داخلی و محیطی قرار دارند که منجر به ایجاد هزاران آسیب به DNA می‌شوند. در صورت عدم ترمیم این آسیب‌ها، ممکن است منجر به جهش و بی‌ثباتی ژنومی شده و توانایی بقای سلول را کاهش داده یا آن را مستعد بیماری‌هایی مانند سرطان سازند. ترمیم DNA در مراحل ابتدایی رشد جنین از اهمیت بالایی برخوردار است (2). در این مرحله، سلول‌ها سطوح مختلفی از حساسیت را نسبت به عوامل درونی و بیرونی نشان می‌دهند (3). پروتئین‌های ترمیمی در همکاری با تنظیم‌کننده‌های چرخه سلولی، با ایجاد توقف در چرخه، به مکانیسم‌های ترمیم اجازه می‌دهند تا DNA آسیب‌دیده را اصلاح کنند  (4). آسیب به DNA در اووسیت‌ها می‌تواند به ناباروری، از دست رفتن بارداری و ناهنجاری‌های مادرزادی منجر شود و مشکلاتی در تولیدمثل ایجاد کند. شکست‌های دو رشته‌ای در DNA (DSBs) از خطرناک‌ترین انواع آسیب‌ها هستند که می‌توانند پایداری ژنوم را مختل سازند. ترمیم صحیح این نوع آسیب‌ها در اووسیت‌ها برای حفظ کیفیت و کمیت تخمک ضروری است، همان‌طور که در مطالعات انسانی و موشی نشان داده شده است(6, 5) . حداقل پنج مسیر اصلی برای ترمیم DNA وجود دارد: ترمیم با جایگزینی باز (BER)، اتصال غیرهمولوگ انتهاها (NHEJ)، ترمیم عدم تطابق (MMR)، ترمیم برداشت نوکلئوتیدی (NER) و نوترکیبی همولوگ .(HR) این فرآیندها در فازهای مختلف چرخه سلولی عمل می‌کنند تا امکان ترمیم مؤثر آسیب‌های DNA فراهم شود (7). ترمیم شکست‌های دو رشته‌ای DNA از طریق نوترکیبی همولوگ(HR)  در حفظ کیفیت و کمیت تخمک‌ها نقش کلیدی دارد (8). علاوه بر این، کمپلکس پروتئینی MRN متشکل از NBS1، RAD50  و MRE11 نقش مهمی در تشخیص آسیب DNA و فعال‌سازی مسیر ATM دارد و از نخستین پاسخ‌دهندگان به شکست‌های دو رشته‌ای به‌شمار می‌رود (9). ژن RAD51 به‌عنوان یکی از اجزای کلیدی مسیر نوترکیبی همولوگ(HR) ، در بازسازی شکست‌های دو رشته‌ای DNA (DSBs)  نقش دارد. این پروتئین با ایجاد رشته‌های نوکلئوپروتئینی روی DNA تک‌رشته‌ای، فرایند جستجوی توالی همولوگ و تبادل رشته را انجام می‌دهد. این مکانیسم برای حفظ ثبات ژنومی در اووسیت‌ها حیاتی است، به‌ویژه در طول میوز زمانی که تقسیم دقیق کروموزوم‌ها ضروری است (10). MRE11 به‌طور فعال در مسیر ترمیم مشارکت دارد و فعالیت نوکلئازی نیز از خود نشان می‌دهد. نقش MRE11 در ترمیم شکست‌های دو رشته‌ای DNA در مراحل ابتدایی میوز I در موش‌ها به خوبی شناخته شده است. اگرچه تحقیقات محدودی در مورد بیان و عملکرد آن در اووسیت‌های مشتق‌شده از فولیکول‌های اولیه انجام شده، اما مطالعات نشان می‌دهند که MRE11 در شناسایی آسیب‌های DNA در تخمک‌های بالغ حیاتی است. مهار عملکرد MRE11 در این موارد، توانایی ترمیم سلول را مختل می‌سازد (11). با پیشرفت فناوری‌های زیستی، مطالعات بیوانفورماتیکی به‌عنوان ابزاری قدرتمند برای شناسایی الگوهای مولکولی، مسیرهای دخیل، و شبکه‌های تعاملات پروتئینی در فرآیندهای زیستی نظیر بلوغ تخمک، جایگاه ویژه‌ای یافته‌اند. ابزارهایی مانند تحلیل بیان ژن(RNA-seq, microarray)، شناسایی مسیرهای زیستی (Pathway analysis)، تحلیل شبکه‌های تعاملی پروتئین–پروتئین (PPI) و پیش‌بینی تنظیمات اپی‌ژنتیکی، امکان بررسی جامع عملکرد ژن‌هایی نظیر RAD51 و MRE11 را در زمینه تولیدمثل فراهم ساخته‌اند (13, 12).  هدف مطالعه حاضر، استفاده از ابزارهای بیوانفورماتیکی برای شناسایی و تحلیل مسیرهای تنظیمی مرتبط با ژن‌هایRAD51  و MRE11 در فرآیند بلوغ تخمک است. این مطالعه با بهره‌گیری از پایگاه‌های داده عمومی و تحلیل‌های سیستماتیک در سطح ترنسکریپتومی، به درک بهتر نقش این ژن‌ها در تنظیم پاسخ‌های ترمیمی در اووسیت‌های مشتق از IVM کمک می‌کند. 
روش بررسی
جمع‌آوری داده‌ها: مطالعات قبلی نشان داده‌اند که ژن‌های RAD51  و MRE11 در فرآیند بلوغ اووسیت نقش حیاتی دارند و اختلال در عملکرد آن‌ها می‌تواند منجر به اختلال در تقسیمات میوزی و کاهش کیفیت ژنومی تخمک شود (11, 14).  در این مطالعه، جهت بررسی مکانیسم‌های مولکولی و تنظیمی این دو ژن در بلوغ تخمک، از ابزارهای بیوانفورماتیکی بهره‌گیری شد. بررسی شبکه اینترکشن پروتئین پروتئین: جهت شناسایی تعاملات مولکولی ژن‌های RAD51 وMRE11 ، از پایگاه داده STRING (Search Tool for the Retrieval of Interacting Genes/Proteins)  استفاده شد. این پایگاه داده، اطلاعات مربوط به تعاملات عملکردی و فیزیکی پروتئین‌ها را فراهم می‌کند و بر اساس منابع متعددی نظیر داده‌های آزمایشگاهی، هم‌بیانی ژن‌ها و تحلیل‌های محاسباتی ساخته شده است (15). در این تحلیل، نام ژن‌ها در نسخه انسان (Homo sapiens) وارد شده و شبکه تعاملی مربوطه با تنظیم confidence score  برابر با 7/0 (سطح بالا) استخراج شد.
آنالیز شبکه PPI درCytoscape: نرم‌افزار  Network Analyzer  برای شناسایی گره‌هایی با درجه بالا(hub) ، مورد استفاده قرار گرفت. گره‌هایی با  hubبیش از ۱ به‌عنوان گره‌های عملکردی بحرانی درنظر گرفته شدند. هر دو نوع گره یعنی  hubوbottleneck  که دارای مرکزیّت نزدیکی بالا هستند، نقش‌های اساسی در فرآیندهای بنیادی سلولی ایفا می‌کنند (16). شبکه به‌دست‌آمده از STRING وارد نرم‌افزار  Cytoscape (v3.9.1) شد و برای تحلیل ساختار شبکه از افزونه NetworkAnalyzer  استفاده گردید. شاخص‌های توپولوژیکی شبکه شامل:
•    ) Degree centrality تعداد اتصالات یک گره با سایر گره‌ها)
تعیین کلاسترها در شبکه PPI : به‌منظور شناسایی خوشه‌ها (کلاسترها) و نواحی به‌شدت درهم‌تنیده در شبکه تعاملی پروتئین-پروتئین(PPI) ، از نرم‌افزار  Cytoscape  و افزونه  MCODE (Molecular Complex Detection) استفاده شد. این ابزار زیرشبکه‌هایی با تراکم بالای اتصالات را شناسایی کرده و ژن‌های بذر (Seed genes) مرتبط با هر کلاستر را مشخص می‌کند. پارامترهای مورد استفاده در این تحلیل شامل :Node Score Cutoff=0/2، Degree Cutoff=2، K-Core=2 و Max-Depth=100 بودند (16). ژن‌های بذر به‌عنوان گره‌هایی با بالاترین درجه در شبکه درنظر گرفته شدند. این تحلیل به‌منظور تعیین ساختارهای عملکردی مهم در شبکه و شناسایی کلاستر‌های کلیدی مرتبط با ژن‌های هدف انجام گرفت.
پیش‌بینی میکرو RNA‌های تنظیم کننده ژن‌ها: جهت شناسایی miRNAهایی که ممکن است تنظیم‌کننده ژن‌های کلیدی شناسایی‌شده در شبکه باشند، از پایگاه داده miRTarBase  استفاده شد. این پایگاه داده، شامل روابط اثبات‌شده بین ژن‌ها و miRNA ها است که از طریق آزمایش‌های تجربی  مانند reporter assay  و (qPCR) اعتبارسنجی شده‌اند. ژن‌های کاندید وارد پایگاه شده و لیستی از miRNAهای تنظیم‌کننده استخراج گردید (17) . 
آنالیز مسیرهای زیستی (Pathway Enrichment): برای شناسایی مسیرهای عملکردی که تحت تأثیر ژن‌های کلیدی شبکه قرار دارند، از سرور  DAVID (Database for  Annotation, Visualization and Integrated Discovery) استفاده شد. ژن‌های هدف وارد سامانه شده و تحلیل‌های Gene Ontology (GO) وKEGG Pathway  جهت بررسی فرآیندهای زیستی، عملکرد مولکولی و مسیرهای سلولی انجام شد. مسیرهایی با مقدار p-value  کمتر از 05/0 از نظر آماری معنادار تلقی گردیدند (18).
نتایج
نتایج حاصل از پایگاه داده  STRING نشان داد که پروتئین‌هایRAD51  و MRE11 با ۱۸ پروتئین دیگر تعامل دارند و 109 اینترکشن دارند (شکل 1). برای تحلیل دقیق‌تر، شبکه استخراج‌شده به نرم‌افزار Cytoscape منتقل شد و با استفاده از ابزار NetworkAnalyzer، شاخص‌های مرکزیّت گره‌ها بررسی گردید. بر اساس نتایج، ژن‌های B‏ARD1  و BRCA1 دارای بالاترین مقدار  Degree بوده و به‌عنوان گره‌های اصلی (hub) شبکه شناخته شدند (جدول 1). با اجرای افزونه MCODE درCytoscape ، دو کلاستر عملکردی با اتصال داخلی بالا در شبکه PPI شناسایی شدند. ژن‌های ATM  و ATRX در این کلاسترها دارای بالاترین  Score بودند، که نشان‌دهنده نقش محوری آن‌ها در کلاستر‌های ژنی است. ویژگی‌های ساختاری این کلاسترها شامل تعداد گره‌ها، یال‌ها، و ژن‌های درگیر در جدول 2 ارائه شده است. با استفاده از پایگاه داده miRTarBase، میکرو  RNAهای متعددی که قادر به تنظیم ژن‌های کلیدی شناسایی‌شده در شبکه هستند، شناسایی شدند. سپس با تحلیل ترکیبی دو شاخص Degree و Betweenness centrality، زیرشبکه‌ای متشکل از ۱۰٪ برتر گره‌ها از کل شبکه استخراج گردید. شبکه تعامل ژن- miRNA در شکل 3 ارائه شده است.
تحلیل غنی‌سازی مسیرها (Pathway Enrichment) با استفاده از پایگاه  DAVID نشان داد که ژن‌های مورد بررسی در چندین مسیر کلیدی زیستی نقش دارند. مهم‌ترین این مسیرها شامل:
•    نوترکیبی همولوگ (Homologous recombination)، اتصال انتهای غیرهمولوگ (Non-homologous end joining)، پیری سلولی ((Cellular senescence، چرخه سلولی (Cell cycle) و مسیر سیگنال‌دهی p53می باشند. اطلاعات دقیق مربوط به مسیرها، ژن‌های درگیر، شماره دسترسی KEGG و مقادیرp-value  در (جدول 3) آورده شده است.
 


شکل 1. شبکه تعاملی پروتئین–پروتئین مربوط به ژن‌های RAD51 و MRE11 استخراج‌شده از پایگاه STRING 


جدول 1: نتایج آنالیز شبکه (Network analysis) در نرم‌افزار Cytoscape شامل مقادیر Degree برای گره‌های کلیدی
نام ژن‌ها Degree
BARD1 19
BRCA1 19
ATRX 10
PALB2 10
TP53 10
FANCD2 10
BLM 10
RAD52 10
CHEK1 10
BRCA2 10
RAD51 10
ATM 10
BRIP1 10
NBN 10
FEN1 10
EXO1 10
RBBP8 10
RAD50 10
XRCC6 10
MRE11 10


شکل 2: کلاسترهای عملکردی استخراج‌شده از شبکه PPI با استفاده از MCODE در .Cytoscapeنودهایی که به رنگ زرد و به شکل دایره نشان داده شده اند seed هستند.

جدول 2: ویژگی‌های کلاسترهای شناسایی‌شده توسط افزونه MCODE در شبکه پروتئینی: تعداد ژن‌ها، یال‌ها و مقادیر Score




 شکل 3. شبکه تعامل ژن- miRNA مربوط به ژن‌های کلیدی حاصل از آنالیز PPI بر اساس داده‌های پایگاه  miRTarBase

جدول 3 : مسیرهای عملکردی زیستی استخراج‌شده از پایگاه DAVID به همراه ژن‌های شرکت‌کننده و مقادیر آماری

 
بحث
تحلیل بیوانفورماتیکی ژن‌های RAD51 وMRE11  در این مطالعه، چشم‌انداز جدیدی از نقش حیاتی آن‌ها در بلوغ تخمک و حفظ سلامت ژنومی در سلول‌های زایشی ارائه داده است. تحقیقات پیشین نشان داده‌اند که اووسیت‌ها ظرفیت بالایی برای ترمیم شکست‌های دو رشته‌ای DNA (DSBs) دارند و مسیر Homologous recombination (HR)، به‌ویژه توسط RAD51 هدایت می‌شود (19). همچنین، نقل‌شده است که مونوکلئاز MRE11 به‌طور فعال در تشخیص و بازآفرینی اولیهDSB ها مشارکت دارد و برای عملکرد مؤثر کامپلکسMRN  ضروری است (20). در مطالعه حاضر، با به‌کارگیری روش‌های تحلیلی بیوانفورماتیکی، از جمله پایگاه STRING  برای ساختاردهی شبکه‌های تعاملی پروتئینی و ابزارهای آنالیز شبکه مانندCytoscape ، توانستیم شبکه‌ای از تعاملات میان ژن‌های دخیل در ترمیمDNA  و کنترل چرخه سلولی را شناسایی کنیم. سپس با تحلیل MCODE و DAVID، مشخص شد که ژن‌های کلیدی این شبکه، در مسیرهای مرتبط با پایداری ژنومی، بلوغ سلولی، پاسخ به آسیب DNA و کنترل تقسیمات میوزی فعالیت می‌کنند. نتایج تحلیل‌های توپولوژیکی شبکه PPI نشان داد که RAD51 با امتیازات بالای Degree به‌عنوان یکhub  کلیدی در شبکه عمل می‌کند و این امر مقبولیت آن را در حفظ انسجام ساختاری شبکه تأیید می‌نماید. این یافته‌ها نشان می‌دهد که RAD51 رابطی حیاتی بین مسیرهای ترمیم DNA ، عملکرد میتوکندری و تقسیم میوز در سلول‌های زایشی است. به‌علاوه، استنادهای بیوانفورماتیکی و مطالعات تجربی بر نقش گسترده‌تر RAD51  در حفظ عملکرد سلولی در سیستم‌های بیولوژیکی پیچیده تأکید دارند (21, 21, 10). در تحلیل شبکه تعاملی پروتئین–پروتئین (PPI) انجام‌شده با استفاده از پایگاه STRING  و نرم‌افزار Cytoscape، ژن‌های  BARD1 و BRCA1 به‌عنوان گره‌هایی با  degree بالا و  betweenness centrality قابل‌توجه شناسایی شدند. حضور این دو ژن در مرکز شبکه، نشان‌دهنده‌ نقش تنظیمی و ساختاری مهم آن‌ها در تعاملات مولکولی مرتبط با پاسخ به آسیبDNA  و حفظ یکپارچگی ژنومی است. ژن‌هایی باdegree  بالا معمولاً در مرکز شبکه‌های زیستی قرار دارند و حذف عملکرد آن‌ها، موجب اختلال جدی در مسیرهای سلولی می‌شود. ژن‌های BRCA1  و BARD1 پس از تشکیل کمپلکس، یکی از ارکان اصلی مسیر  HR به‌شمار می‌آیند. مطالعه Westermark  و همکاران در سال ۲۰۰۳ نشان داد که کمپلکس BRCA1–BARD1 به‌صورت فیزیکی با پروتئین RAD51  در فرایند ترمیم DSB  مشارکت دارد وعدم وجود این کمپلکس موجب بی‌نظمی کروموزومی و نقص در ترمیم HR می‌شود (23). علاوه بر این، پژوهش Boulton در سال ۲۰۰4 تأیید کرد که دسترسی BRCA1–BARD1 به محل شکست DNA برای تشکیل اتصال‌های هومولوگیک ضروری است و فقدان آن باعث افزایش حساسیت به آسیب‌های ژنتیکی می‌شود (24). این یافته‌ها نشان‌دهنده نقش برجسته BRCA1–BARD1 در پایداری ژنومی اووسیت و جلوگیری از انتقال جهش به نسل بعد است. مطالعات بالینی نشان داده‌اند که زنان حامل موتاسیونBRCA1  دارای سطح پایین‌تر هورمون AMH به‌عنوان نشانگری برای ذخیره تخمدان هستند و این امر نشان‌دهنده کاهش ذخیره اولیه اووسیت است (26, 25). در مدل حیوانی، مطالعه Titus و همکاران نشان داد که کاهش عملکردBRCA1 ، منجر به شکل‌گیری کمتر فولیکول‌های پریمیوردیال و تجمع بیشتر DNA DSB  در اووسیت‌های موش می‌شود، که این امر با کاهش ذخیره تخمک و شروع زودرس یائسگی مرتبط است (27). مطالعات سامانه‌ای و مروری دیگری نیز چنین یافته‌هایی را تأیید می‌کنند: بررسی BRCA و ترمیم DNA نشان می‌دهد که نقص عملکرد BRCA1 موجب پیری زودرس تخمدان، افت کیفیت تولیدمثل و کاهش طول عمر باروری می‌شود  (28). اجرای افزونه MCODE در Cytoscape، ژن  ATM را در یکی از کلاسترهای اصلی شبکه PPI با بالاترین امتیاز مکان‌یابی کرده است. ‌مطالعه‌ای ژنتیکی در اووسیت‌های موش نشان داده که کاهش بیان ATM منجر به اختلال در پیشرفت میوز، افزایش تجمع  DSB، ایجاد آنیوپلوئیدی و بروز آپوپتوز زودرس می‌شود. این یافته‌ها حاکی از آن است که ATM نقش حفاظتی کلیدی در مقابل استرس ژنومی دارد و ضعف عملکرد آن می‌تواند به‌عنوان یک عامل مهم در پیر شدن تخمدان و نازایی تلقی گردد (29). ژن  ATRXنیز در کلاستر دوم با امتیاز بالا قرار گرفت؛ این ژن به‌عنوان یک بازآرایی‌کننده کروماتین، برای تشکیل صحیح هتروکروماتین اطراف سانترومرها ضروری است. مطالعات روی موش نشان داده که فقدان ATRX منجر به اختلال در تراز کروموزومی میتوز دوم و افزایش چشمگیر آنیوپلوئیدی تا حدود ۶۰٪ اووسیت‌ها می‌شود. این ناپایداری کروموزومی با کاهش قدرت باروری همراه است و اهمیت ATRX را در بلوغ صحیح و سلامت اووسیت به‌خوبی نشان می‌دهد (30). یافته‌های بیوانفورماتیکی و مولکولی نشان داد کهATM  وATRX  دو ژن کلیدی در شبکه PPI هستند که نقش عمده‌ای در حفظ کیفیت ژنومی، پیشرفت دقیق میوز و جلوگیری از پیری کروموزومی در اووسیت‌ها دارند. این شناسایی علاوه بر آگاهی علمی، پتانسیل استفاده از این ژن‌ها در طراحی بیومارکرهای تشخیصی یا اهداف درمانی برای بهبود باروری و کنترل ذخیره تخمک را نمایان می‌کند. نتایج تحقیق حاضر تأکید می‌کند که ژن RAD51 نه فقط در پایداری ژنومی، بلکه در بلوغ اووسیت نقش کلیدی ایفا می‌کند. این ژن به‌عنوان نوترکیب‌کننده مرکزی (Recombinase) در مسیر  HR شناخته شده است و مطالعات متعددی نشان داده‌اند که کاهش بیان RAD51 باعث افزایش DSB، آپوپتوز و اختلال در بلوغ میواز در اووسیت‌های موش می‌شود (32, 31). برای مثال، Kim  و همکاران در سال 2016 با استفاده از مهار RNAi در اووسیت‌های موش نشان دادند که خاموش‌سازی RAD51 منجر به توقف میوز در مرحله MI، بی‌نظمی در تراز کروموزومی، آسیب به دوک تقسیم و کاهش عملکرد میتوکندری شد (10) همچنین، Kujjo و همکاران در سال 2010 با میکروتزریق RAD51 ریکامبیننت به اووسیت‌های AKR/J  دریافتند که این افزایش بیان موجب کاهش قابل‌توجهِ شکست‌های DSB و سرکوب مسیر آپوپتوتیک Bax شد و توسعه جنینی را بهبود بخشید (22). مطالعات تجربی نشان داده‌اند که RAD51 نقشی حیاتی در حفظ عملکرد میتوکندری در طول بلوغ اووسیت‌ها ایفا می‌کند. در یک کار تجربی توسط Kim و همکاران در سال ۲۰۱۶، کاهش بیان RAD51 در اووسیت‌های موش موجب کاهش تولید ATP، افت پتانسیل غشایی میتوکندری، تغییر نامناسب در توزیع میتوکندری و افزایش اتوفاژی شد (10).  این یافته‌ها نشان می‌دهند که RAD51  علاوه بر نقش شناخته‌شده‌اش در ترمیم DNA، یک نقش مرتبط با متابولیسم انرژی و بقای سلولی در اووسیت‌ها دارد. تحلیل غنی‌سازی مسیرهای عملکردی با استفاده از DAVID، مسیرهای کلیدی مرتبط با ژن‌های RAD51 و MRE11 را شناسایی کرد. شاخص‌ترین مسیر شناسایی‌شده، مسیر  HR است که شامل ژن‌هایی مانند RAD51،MRE11، RAD50، BRCA1  و BRCA2  می‌شود. این مسیر نقش بسیار مهمی در ترمیم دقیق شکست‌های دو رشته‌ای DNA (DSBs)  دارد و در سلول‌های زایشی، به‌ویژه اووسیت‌ها، مسئول حفظ پایداری ژنومی است. اختلال در این مسیر با افزایش نرخ آنیوپلوئیدی در اووسیت‌های انسانی و موشی و کاهش تعداد فولیکول‌های اولیه مرتبط است (10). مسیر NHEJ  علاوه بر نقش درHR، در مسیرmicrohomology-mediated  end joining (MMEJ)  نیز مشارکت دارد. مطالعات تجربی نشان داده‌اند که فقدان کارکرد MRE11 باعث نقص در پردازش انتهای DNA و کاهش ترمیم MMEJ می‌شود، که بر ناکارآمدی ترمیم DSB تأکید دارد (33). تحلیل مسیرهای زیستی نشان داد که مولکول‌های کلیدی ATM، TP53، CHK1/2  به‌طور معناداری در پاسخ به آسیب DNA فعال شده‌اند. ژن ATM، به‌عنوان آغازگر مسیر پاسخ به شکست‌های دو رشته‌ای DNA، از طریق کمپلکس MRN فعال می‌شود و سپس پروتئین‌هایی نظیر  H2AXγ،  CHK2 و p53  را فسفریله می‌کند. این فرآیند منجر به توقف چرخه سلولی، فعال‌سازی نقاط کنترل و در صورت لزوم، آغاز آپوپتوز می‌گردد. در اووسیت‌ها، نشانه‌هایی از فعال‌شدن CHK1 و CHK2 در مرحله G2/M پس از آسیب DNA  گزارش شده که نشان‌دهنده عملکرد حفاظتی این مسیر در حفظ کیفیت تخمک است (34). 
مسیرهای مرتبط با Cellular senescence  و checkpoints چرخه سلولی نیز در آنالیز غنی‌سازی مشاهده شدند. این مسیرها با حذف سلول‌های دارای آسیب DNA از طریق توقف چرخه در فازهای G2/M یاG1/S، نقشی کلیدی در ارتقای کیفیت اووسیت‌ها دارند. این نتایج با مطالعاتی که نشان می‌دهند DNA damage checkpoint در اووسیت‌ها، مخصوصاً در مرحله میوز I، توقف‌های حساسیت‌زا ایجاد می‌کند هم‌راستا هستند (35, 34). کاهش توان ترمیم DNA در اووسیت‌های مسن‌تر و افزایش DSB ها با علائم پیری سلولی هم گزارش شده است(Changes in DNA repair ). همچنین، ژنMRE11  که معمولاً به‌عنوان یک عامل کلیدی در شناسایی و ترمیم شکست‌های دو رشته‌ای DNA شناخته می‌شود، در مسیر پیری سلولی نیز نقش دارد. فعالیت این ژن به‌همراه  ATM و TP53 موجب فعال‌سازی پاسخ به آسیبDNA ، توقف چرخه سلولی و در صورت لزوم القای آپوپتوز می‌شود که از مکانیسم‌های اصلی پیری سلولی در اووسیت‌ها محسوب می‌گردد (20). از منظر تنظیم‌های اپی‌ژنتیکی، بررسی تعامل ژن–میکرو‌RNA  نشان داد که تعدادی از miRNAها قادر به هدف‌گیری ژن‌های کلیدی شبکه مانندRAD51 ، MRE11 و BRCA1  هستند. برای درک بهتر، ابتدا باید به نقش بنیادی microRNAها در تنظیم پس از رونویسی ژن‌ها اشاره کرد. این مولکول‌های کوچک RNA با اتصال به mRNA هدف، موجب سرکوب ترجمه یا تخریب پیام‌رسان می‌شوند و بدین ترتیب در مسیرهای تمایز، چرخه سلولی و پاسخ به آسیب‌های  DNA نقش دارند (36). حضور این miRNAها به‌عنوان تنظیم‌کننده‌های منفی یا مثبت در مسیرهای ترمیم، اهمیت بالایی در پایداری ژنومی و تنظیم بلوغ اووسیت دارد. تحلیل miRNA ژن نشان داد که چند miRNA کلیدی می‌توانند بیان RAD51 را هدف قرار دهند و از این طریق بر مسیر Homologous Recombination  تأثیر بگذارند:
•     Gasparini و همکاران در سال 2014 نشان دادند که افزایش بیان miR 155 در سلول‌های سرطان سینه منجر به کاهش بیان RAD51  و اختلال در ترمیم HR می‌شود، به‌طوری که miR 155 به‌طور مستقیم روی 3′ UTR RAD51 اثرگذار است و سطح پروتئین آن را کاهش می‌دهد (37).
•     Chen و همکاران در سال 2019 یافتند که miR 34a  از طریق اتصال به mRNA RAD51، بیان آن را کاهش می‌دهد. این مکانیسم باعث اختلال در HR و تأثیر بر مسیر کارکرد p53 شده و منجر به افزایش حساسیت سلولی به آسیب DNA می‌گردد (38).
مطالعات متعددی نشان داده‌اند که miRNAهایی مانند miR 155، miR 34a، miR 506، miR 182، و  miR 103/107  مستقیماً روی RAD51 تأثیر می‌گذارند و باعث کاهش تشکیل فسیل Fociهای RAD51 در پاسخ به آسیب DNA و کاهش کارایی مسیر  Homologous Recombination  می‌شوند (39). همچنین مطالعات مروری در جانوران پستاندار (شامل انسان، گاو و موش) نشان داده‌اند که miRNAها در طول بلوغ اووسیت و پیش‌روی مراحل اولیه جنینی نقش تنظیمی مهمی دارند. این نقش شامل تنظیم دقیق بیان ژن‌های مرتبط با چرخه سلولی، آپوپتوز و ترمیمDNA  است (41, 40). افزون برRAD51، مطالعات نشان داده‌اند که چندین miRNA نظیرmiR-103  وmiR-107  نیز می‌توانند ژن MRE11  را هدف قرار دهند. این تنظیم منفی موجب کاهش تشکیل کمپلکس MRN و تضعیف پاسخ به آسیب‌های DNA می‌شود که به‌ویژه در سلول‌های زایشی می‌تواند منجر به افت کیفیت اووسیت گردد (21). لذا، شایع بودن miRNAهای هدفمند برای RAD51 و MRE11 در این مسیرها، نشان‌دهنده‌ یک لایه تنظیمی حیاتی در کیفیت اووسیت است. یافته‌های این مطالعه، ضمن شناسایی اهداف بالقوه درمانی یا تشخیصی، می‌توانند زمینه را برای طراحی راه‌کارهای بیولوژیکی نوین در بهبود عملکرد باروری، به‌ویژه در بیماران نابارور یا تحت درمان‌های کمک‌باروری فراهم سازند. هم‌چنین، بررسی دقیق‌تر بیان ژن و تنظیم‌های miRNA در مدل‌های تجربی می‌تواند افق‌های جدیدی برای مداخله در مسیرهای مولکولی مؤثر بر سلامت تخمک‌ها ایجاد کند.
نتیجه‌گیری
این مطالعه با بهره‌گیری از رویکرد بیوانفورماتیکی، نقش بنیادی ژن‌هایRAD51  وMRE11  را در فرآیند بلوغ تخمک، ترمیم DNA و حفظ سلامت ژنومی اووسیت‌ها روشن ساخت. تحلیل شبکه PPI و مسیرهای زیستی، حضور محوری ژن‌هایی نظیر BRCA1، BARD1،ATM  و ATRX را در فرآیندهای کلیدی مرتبط با باروری تأیید نمود. همچنین، یافته‌ها نشان دادند که تنظیم اپی‌ژنتیکی از طریق miRNAها، لایه‌ای حیاتی از کنترل بیان ژن‌های مرتبط با ترمیم را فراهم می‌کند. داده‌های حاصل از این پژوهش می‌توانند مبنایی برای توسعه بیومارکرهای تشخیصی یا اهداف درمانی جهت بهبود عملکرد تخمدانی و مدیریت ناباروری فراهم سازند. پیشنهاد می‌شود در مطالعات آینده، اعتبارسنجی تجربی داده‌های استخراج‌شده و بررسی عملکردی  miRNAهای کلیدی در مدل‌های حیوانی و انسانی پیگیری شود.

سپاس‌گزاری
این مقاله برگرفته از رساله دکتری با عنوان «بررسی میزان ترمیم DNA در تخمک‌های موشی حاصل از بلوغ آزمایشگاهی» است.
نویسندگان از همکاری و مشارکت یکدیگر در انجام این مطالعه و نگارش مقاله صمیمانه قدردانی می‌کنند.
حامی مالی
ندارد.
تعارض در منافع
وجود ندارد.
ملاحظات اخلاقی
این مطالعه شامل کار با داده‌های عمومی و تحلیل بیوانفورماتیکی بود و نیازی به دخالت انسانی یا حیوانی نداشت. کلیه روش‌ها مطابق با استانداردهای اخلاقی پژوهش‌های علمی انجام شده است.
مشارکت نویسندگان
تمام نویسندگان در ایده‌پردازی، جمع‌آوری داده‌ها، تجزیه و تحلیل بخش‌های مختلف طرح و همچنین نگارش و ویرایش علمی مقاله مشارکت داشتند و سهم یکسانی در انجام این فعالیت‌ها دارند. 
References:
1.    Thaker N, Dhande R, Parihar P, Dhande R. Role of transvaginal sonography in the diagnosis of female infertility: a comprehensive review. Cureus 2023; 15(12): e50505.  
2.    Khokhlova E, Fesenko Z, Sopova J, Leonova E. Features of DNA repair in the early stages of mammalian embryonic development. Genes. 2020; 11(10): 1138.  
3.    Hamdoun A, Epel D. Embryo stability and vulnerability in an always changing world. Proc Natl Acad Sci USA. 2007; 104(6): 1745-1750.  
4.    Smith ML, Fornace AJ Jr. Mammalian DNA damage-inducible genes associated with growth arrest and apoptosis. Mutat Res Rev Genet Toxicol 1996; 340(2-3): 109-124.  
5.    Oktay K, Moy F, Titus S, Stobezki R, Turan V, Dickler M, et al. Age-related decline in DNA repair function explains diminished ovarian reserve, earlier menopause, and possible oocyte vulnerability to chemotherapy in women with BRCA mutations. J Clin Oncol. 2014; 32(10): 1093-1094.  
6.    Titus S, Li F, Stobezki R, Akula K, Unsal E, Jeong K, et al. Impairment of BRCA1-related DNA double-strand break repair leads to ovarian aging in mice and humans. Sci Transl Med. 2013; 5(172): 172ra21.  
7.    Chatterjee N, Walker GC. Mechanisms of DNA damage, repair, and mutagenesis. Environ Mol Mutagen. 2017; 58(5): 235-263.  
8.    He DJ, Wang L, Zhang ZB, Guo K, Li JZ, He XC, et al. Maternal gene Ooep may participate in homologous recombination-mediated DNA double-strand break repair in mouse oocytes. Zool Res. 2018; 39(6): 387-397.  
9.    Wang X, Liu D, He D, Suo S, Xia X, He X, et al. Transcriptome analyses of rhesus monkey preimplantation embryos reveal a reduced capacity for DNA double-strand break repair in primate oocytes and early embryos. Genome Res. 2017; 27(4): 567-579.  
10.    Kim KH, Park JH, Kim EY, Ko JJ, Park KS, Lee KA. The role of Rad51 in safeguarding mitochondrial activity during the meiotic cell cycle in mammalian oocytes. Sci Rep. 2016; 6(1): 34110.  
11.    Mayer A, Baran V, Sakakibara Y, Brzakova A, Ferencova I, Motlik J, et al. DNA damage response during mouse oocyte maturation. Cell Cycle. 2016; 15(4): 546-558.  
12.    Primig M. The bioinformatics tool box for reproductive biology. Biochim Biophys Acta Mol Basis Dis. 2012; 1822(12): 1880-1895.  
13.    Mulas F, Sacchi L, Zagar L, Garagna S, Zuccotti M, Zupan B, et al. Knowledge-based bioinformatics for the study of mammalian oocytes. Int J Dev Biol. 2012; 56(10-11-12): 859-866.  
14.    Pailas A, Niaka K, Zorzompokou C, Marangos P. The DNA damage response in fully grown mammalian oocytes. Cells. 2022; 11(5): 798.  
15.    Dehghan Z, Mirmotalebisohi SA, Sameni M, Bazgiri M, Zali H. A motif-based network analysis of regulatory patterns in Doxorubicin effects on treating breast cancer, a systems biology study. Avicenna J Med Biotechnol. 2022; 14(2): 137-145.  
16.    Saberi F, Dehghan Z, Noori E, Taheri Z, Sameni M, Zali H. Identification of critical molecular factors and side effects underlying the response to thalicthuberine in prostate cancer: a systems biology approach. Avicenna J Med Biotechnol. 2023; 15(1): 53-61.  
17.    Cui S, Yu S, Huang HY, Lin YCD, Huang Y, Zhang B, et al. miRTarBase. 2025: updates to the collection of experimentally validated microRNA-target interactions. Nucleic Acids Res. 2025; 53(D1): D147-D156.  
18.    Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, et al. DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007; 35(Suppl 2): W169-W175.  
19.    Kim KH, Park JH, Kim EY, Ko JJ, Park KS, Lee KA. The role of Rad51 in safeguarding mitochondrial activity during the meiotic cell cycle in mammalian oocytes. Sci Rep. 2016; 6: 34110.  
20.    Buis J, Wu Y, Deng Y, Leddon J, Westfield G, Eckersdorff M, et al. Mre11 nuclease activity has essential roles in DNA repair and genomic stability distinct from ATM activation. Cell. 2008; 135(1): 85-96.  
21.    Peraza-Vega RI, Valverde M, Rojas E. Interactions between miRNAs and double-strand breaks DNA repair genes, pursuing a fine-tuning of repair. Int J Mol Sci. 2022; 23(6): 3231.  
22.    Kujjo LL, Laine T, Pereira RJ, Kagawa W, Kurumizaka H, Yokoyama S, et al. Enhancing survival of mouse oocytes following chemotherapy or aging by targeting Bax and Rad51. PLoS One. 2010; 5(2): e9204.  
23.    Westermark UK, Reyngold M, Olshen AB, Baer R, Jasin M, Moynahan ME. BARD1 participates with BRCA1 in homology-directed repair of chromosome breaks. Mol Cell Biol. 2003; 23(21): 7926-7936.  
24.    Boulton SJ, Martin JS, Polanowska J, Hill DE, Gartner A, Vidal M. BRCA1/BARD1 orthologs required for DNA repair in Caenorhabditis elegans. Curr Biol. 2004; 14(1): 33-39.  
25.    Wang X, Wang L, Xiang W. Mechanisms of ovarian aging in women: a review. J Ovarian Res. 2023; 16(1): 67.  
26.    Wang ET, Pisarska MD, Bresee C, Chen YDI, Lester J, Afshar Y, et al. BRCA1 germline mutations may be associated with reduced ovarian reserve. Fertil Steril. 2014; 102(6): 1723-1728.  
27.    Lin W, Titus S, Moy F, Ginsburg ES, Oktay K. Ovarian aging in women with BRCA germline mutations. J Clin Endocrinol Metab. 2017; 102(10): 3839-3847.  
28.    Oktay K, Turan V, Titus S, Stobezki R, Liu L. BRCA mutations, DNA repair deficiency, and ovarian aging. Biol Reprod. 2015; 93(3): 67, 1-10.  
29.    Suzuki R, Tan X, Szymanska KJ, Kubikova N, Perez CA, Wells D, et al. The role of declining ataxia-telangiectasia-mutated (ATM) function in oocyte aging. Cell Death Discov. 2024; 10(1): 302.  
30.    Baumann C, Viveiros MM, De La Fuente R. Loss of maternal ATRX results in centromere instability and aneuploidy in the mammalian oocyte and pre-implantation embryo. PLoS Genet. 2010; 6(9): e1001137.  
31.    Kim K, Park J, Kim E, Ko J, Park K, Lee K. Silencing of Rad51 impairs oocyte maturation via mitochondrial dysfunction and mitophagy. Fertil Steril. 2016; 106(3): e188.  
32.    Borde V. The multiple roles of the Mre11 complex for meiotic recombination. Chromosome Res. 2007; 15(5): 551-563.  
33.    Xie A, Kwok A, Scully R. Role of mammalian Mre11 in classical and alternative nonhomologous end joining. Nat Struct Mol Biol. 2009; 16(8): 814-818.  
34.    Carroll J, Marangos P. The DNA damage response in mammalian oocytes. Front Genet. 2013; 4: 117.  
35.    Subramanian GN, Greaney J, Wei Z, Becherel O, Lavin M, Homer HA. Oocytes mount a noncanonical DNA damage response involving APC-Cdh1-mediated proteolysis. J Cell Biol. 2020; 219(4): e201907213.  
36.    Salehi M. The expression of miR-31 and its target gene FOXP3 in recurrent implantation failure patients. Int J Womens Health Reprod. Sci. 2020; 8: 389-395.  
37.    Gasparini P, Lovat F, Fassan M, Casadei L, Cascione L, Jacob NK, et al. Protective role of miR-155 in breast cancer through RAD51 targeting impairs homologous recombination after irradiation. Proc Natl Acad Sci USA. 2014; 111(12): 4536-4541.  
38.    Chen S, Liu R, Wang Q, Qi Z, Hu Y, Zhou P, et al. MiR-34s negatively regulate homologous recombination through targeting RAD51. Arch Biochem Biophys. 2019; 666: 73-82.  
39.    Peraza-Vega RI, Valverde M, Rojas E. Interactions between miRNAs and double-strand breaks DNA repair genes, pursuing a fine-tuning of repair. Int J Mol Sci. 2022; 23(6): 3231.  
40.    Dehghan Z, Mohammadi-Yeganeh S, Rezaee D, Salehi M. MicroRNA-21 is involved in oocyte maturation, blastocyst formation, and pre-implantation embryo development. Dev Biol. 2021; 480: 1-10.  
41.    Dehghan Z, Mohammadi-Yeganeh S, Salehi M. MiRNA-155 regulates cumulus cells function, oocyte maturation, and blastocyst formation. Biol Reprod. 2020; 103(3): 548-559.


 
نوع مطالعه: پژوهشي | موضوع مقاله: زیست شناسی
دریافت: 1404/11/27 | پذیرش: 1405/2/16 | انتشار: 1405/5/15

ارسال نظر درباره این مقاله : نام کاربری یا پست الکترونیک شما:
CAPTCHA

ارسال پیام به نویسنده مسئول


بازنشر اطلاعات
Creative Commons License این مقاله تحت شرایط Creative Commons Attribution-NonCommercial 4.0 International License قابل بازنشر است.

کلیه حقوق این وب سایت متعلق به ماهنامه علمی پ‍ژوهشی دانشگاه علوم پزشکی شهید صدوقی یزد می باشد.

طراحی و برنامه نویسی : یکتاوب افزار شرق

© 2026 CC BY-NC 4.0 | SSU_Journals

Designed & Developed by : Yektaweb