ژن‌های مرتبط با هورمون‌های جنسی در استئوآرتریت: شناسایی بیوانفورماتیک و ارتباط با نفوذ ایمنی

مقدمه

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

هدف مطالعه

مطالعه منتشرشده در PLOS One با استفاده از روش‌های بیوانفورماتیک و یادگیری ماشین تلاش کرد تا ژن‌های نامزد مرتبط با هورمون‌های جنسی در OA را شناسایی کند و ویژگی‌های عملکردی و تعاملات ایمنی آن‌ها را بررسی نماید. این مطالعه بر اساس تحلیل داده‌های بیان ژنی عمومی و پایگاه‌های اطلاعاتی ژنی انجام شده است و هدف آن ارائه شواهد انجمنیاتی (associative) برای نقش احتمالی این ژن‌ها در OA بود.

خلاصه روش‌شناسی

داده‌ها و انتخاب ژن‌ها

برای تحلیل، داده‌های بیان ژنی OA از پایگاه GEO استخراج شد و فهرستی از ژن‌های مرتبط با هورمون‌های جنسی از پایگاه GeneCards گردآوری شد. این ترکیب داده‌ها امکان پیدا کردن ژن‌هایی را فراهم کرد که هم در OA و هم در ارتباط با هورمون‌های جنسی مطرح‌اند.

شناسایی ژن‌های متفاوت بیان

با استفاده از بسته آماری limma فهرست ژن‌های متفاوت‌بیان (DEGs) بدست آمد. سپس از تقاطع میان این DEGs و ژن‌های مرتبط با هورمون جنسی، ژن‌هایی تحت عنوان SADEGs (sex hormone-associated DEGs) استخراج شد.

تحلیل شبکه‌ای و هم‌بیان

برای تعیین گروه‌های ژنی مرتبط با OA، از تحلیل شبکه هم‌بیان وزن‌دار یا WGCNA استفاده شد. ماژول یا گروه‌هایی که بیشترین همبستگی را با شاخص‌های بالینی یا وضعیت OA داشتند، انتخاب و ژن‌های مرتبط به عنوان OA-SADEGs تعیین شدند.

غنی‌سازی عملکردی

برای تفسیر عملکردی ژن‌ها از تحلیل‌های GO و KEGG استفاده شد تا مسیرهای زیستی و فرآیندهای بیولوژیکی غالب برای این ژن‌ها مشخص شود. نتایج تمرکز بر مسیرهایی مانند بازسازی ماتریکس خارج‌سلولی، سیگنال‌دهی گیرنده و ارائه آنتی‌ژن را نشان داد.

شناسایی ژن‌های هاب با یادگیری ماشین

سه مدل یادگیری ماشین برای یافتن ژن‌های هاب که بیشترین توان تفکیک را دارند به کار گرفته شد. ژن‌هایی که در هر سه مدل به‌عنوان کاندیدای مهم ظاهر شدند به‌عنوان هاب نهایی گزارش شدند.

تحلیل نفوذ ایمنی

با استفاده از ابزار xCell میزان نفوذ انواع سلول‌های ایمنی در نمونه‌ها برآورد شد و ارتباط بین امتیازهای نفوذ ایمنی و بیان ژن‌های هاب بررسی شد.

اعتبارسنجی خارجی و آزمایشی

برای اعتبارسنجی نتایج، بیان ژن‌های هاب در دیتاست‌های مستقل بررسی شد و در نمونه‌های خون محیطی بیماران OA با روش qRT-PCR سنجیده شد تا تغییر بیان این ژن‌ها در نمونه‌های بالینی ارزیابی گردد.

نتایج کلیدی

از تحلیل‌ها ۳۲ ژن مرتبط با هورمون‌های جنسی در زمینه OA شناسایی شد که در مسیرهایی مانند بازسازی ماتریکس خارج‌سلولی، سیگنال‌دهی گیرنده و ارائه آنتی‌ژن غنی شدند. سه ژن به‌عنوان کاندیدای هاب معرفی شدند: LOXL1، HLA-DRA و CYBB. بیان هر سه ژن در نمونه‌های OA در مقایسه با نمونه‌های کنترل افزایش نشان داد.

تحلیل xCell نشان داد که ۱۳ نوع سلول ایمنی در بافت‌های مرتبط با OA به‌صورت متفاوت بین گروه‌ها توزیع شده‌اند و ارتباط آماری میان برخی از این انواع سلول‌ها و بیان ژن‌های هاب وجود داشت. اعتبارسنجی در دیتاست‌های مستقل و نتایج qRT-PCR از افزایش بیان این سه ژن در نمونه‌های بالینی حمایت کرد.

تفسیر بیولوژیک

معنی یافتن LOXL1، HLA-DRA و CYBB

• LOXL1: عضو خانواده لیزیل اکسیدازها که در کراس‌لینکینگ و بازسازی ماتریکس خارج‌سلولی نقش دارد. افزایش بیان آن می‌تواند نشان‌دهنده تغییرات در ساختار و خواص مکانیکی بافت‌های مفصلی باشد.
• HLA-DRA: مولکولی مرتبط با ارائه آنتی‌ژن نوع II که نشان‌دهنده فعال‌سازی مسیرهای ایمنی و پیام‌رسانی التهابی در میکروسکوپ مفصل است.
• CYBB: زیرواحدی از آنزیم NADPH اکسیداز (NOX2) که در تولید گونه‌های واکنشی اکسیژن و پاسخ‌های التهابی نقش دارد؛ افزایش آن می‌تواند به استرس اکسیداتیو و آسیب بافتی منجر شود.

در مجموع، این سه ژن ترکیبی از عملکردهای ساختاری، ایمنی و اکسیداتیو را نشان می‌دهند که ممکن است به توضیح برخی از جنبه‌های مولکولی تأثیر هورمون‌های جنسی بر مسیرهای بیماری‌زایی OA کمک کنند.

محدودیت‌های مطالعه

اگرچه این مطالعه با رعایت استانداردهای بیوانفورماتیک انجام شده، چند محدودیت مهم وجود دارد که باید مدنظر قرار گیرد:

• داده‌ها برگرفته از مجموعه‌های عمومی هستند و ممکن است از نظر تکنیکی و جمعیتی همگن نباشند. تفاوت در روش جمع‌آوری، تعریف نمونه و شرایط آزمایشی می‌تواند بر نتایج تأثیر بگذارد.
• یافته‌ها از نوع همبستگی و انجمنیاتی هستند؛ این نتایج علت و معلول را اثبات نمی‌کنند و نیاز به آزمایش‌های مکانیکی و مدل‌های تجربی برای تایید دارد.
• تحلیل نفوذ ایمنی با ابزارهای محاسباتی مانند xCell برآورد نظری از حضور سلول‌ها ارائه می‌دهد و باید با روش‌های تجربی مانند ایمونوهیستوشیمی یا سیتومتری فلوئورسانت تأیید شود.
• اعتبارسنجی qRT-PCR در خون محیطی ممکن است نشان‌دهنده تغییراتی سیستمیک باشد که دقیقاً بازتاب‌دهنده اتفاقات داخل مفصل نیستند؛ اختلاف میان بافت و خون باید با احتیاط تفسیر شود.
• اطلاعات مربوط به نقش مستقیم هورمون‌ها و تأثیرات تفکیکی جنسیتی نیاز به مطالعات طولی و بالینی با داده‌های هورمونی دقیق دارد که در این مطالعه به‌طور مستقیم بررسی نشده است.

کاربردهای بالینی بالقوه و مسیرهای آینده

یافته‌های این مطالعه می‌تواند چند مسیر تحقیقاتی و بالینی را هدایت کند:

• این ژن‌ها ممکن است به عنوان بیومارکرهای تشخیصی یا پیش‌آگهی بالقوه مورد بررسی قرار گیرند؛ با این حال نیاز به اعتبارسنجی در گروه‌های بزرگتر و نمونه‌های بافتی مفصلی است.
• مداخلات دارویی یا مولکولی که مسیرهای مرتبط با LOXL1، HLA-DRA یا CYBB را هدف قرار می‌دهند می‌توانند در مدل‌های پیش‌بالینی ارزیابی شوند تا اثبات مفهومی (proof-of-concept) برای تغییر مسیر بیماری فراهم کنند.
• مطالعات آینده باید به تفکیک نقش هورمون‌های جنسی (مثلاً استروژن، تستوسترون، پروژسترون) در تنظیم بیان این ژن‌ها بپردازند و تفاوت‌های جنسیتی را به‌صراحت بررسی کنند.
• ترکیب داده‌های بیولوژیکی با اطلاعات بالینی مانند شدت درد، یافته‌های تصویربرداری و وضعیت هورمونی ممکن است به توسعه مدل‌های پیش‌بینی‌کننده قابل‌اطمینان‌تر برای مدیریت شخصی‌شده بیماران OA کمک کند.

نکته مهم برای بیماران

اگر شما یا نزدیکانتان مبتلا به استئوآرتریت هستید، مهم است بدانید که این مطالعه بر پایه تحلیل داده‌های زیستی و محاسباتی است و نتایج آن هنوز در مراحل تحقیقاتی قرار دارد. افزایش بیان ژن‌هایی مانند LOXL1، HLA-DRA و CYBB ارتباط احتمالی با مسیرهای ساختاری و ایمنی در OA را نشان می‌دهد، اما این یافته‌ها به هیچ‌وجه به معنی تغییر فوری در مراقبت بالینی یا توصیه درمانی نیستند. هرگونه تغییر درمانی حتماً باید توسط پزشک معالج و بر اساس شواهد بالینی تایید شود.

نکات عملی برای پزشکان و محققان

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

جمع‌بندی

تحلیل بیوانفورماتیک ارائه‌شده در این مطالعه نشان می‌دهد که مجموعه‌ای از ژن‌های مرتبط با هورمون‌های جنسی در استئوآرتریت تغییر بیان دارند و سه ژن LOXL1، HLA-DRA و CYBB به‌عنوان کاندیداهای کلیدی برجسته شده‌اند. این ژن‌ها با مسیرهایی مرتبطند که می‌توانند تلاقی میان فرایندهای ساختاری، ایمنی و استرس اکسیداتیو را در مفصل نشان دهند. هرچند یافته‌ها امیدوارکننده‌اند، اما لازم است این نتایج در مطالعات تجربی و بالینی بزرگتر و با روش‌های مولکولی و سلولی دقیق‌تر تایید شوند تا بتوانند به کاربرد بالینی مانند بیومارکر یا هدف درمانی تبدیل شوند.

منبع

Identification of candidate sex hormone-associated genes and immune infiltration characteristics in osteoarthritis based on bioinformatics analysis and machine learning. PLOS One 2026. DOI: https://doi.org/10.1371/journal.pone.0351556