توالی‌یابی کل اگزوم (WES) در ژن نگار آیندگان

١. مقدمه و زمینه بالینی

توالی‌یابی کل اگزوم (Whole Exome Sequencing - WES) یک فناوری ژنومیک تحول‌آفرین است که به‌طور انتخابی نواحی کدکننده پروتئین (اگزون‌ها) در ژنوم انسان را تعیین توالی می‌کند. اگرچه اگزوم تنها توالی حدود ۱.۵ تا ۲ درصد از کل ژنوم انسان را تشکیل می‌دهد، اما حامل بالغ بر ۸۵ درصد از جهش‌های شناخته‌شده و بیماری‌زای ژنتیکی می‌باشد که این امر آن را به یک ابزار تشخیصی بسیار کارآمد تبدیل می‌سازد. WES بالینی، پارادایم‌های تشخیصی برای اختلالات نادر مندلی، ناهنجاری‌های عصبی-تکاملی و فنوتیپ‌های پیچیده کودکان را به طور بنیادی تغییر داده است. با ارزیابی ده‌ها هزار ژن به صورت همزمان، این رویکرد نامتناهی محدودیت‌های آزمایش‌های سنتی و زمان‌بر تک‌ژنی را دور می‌زند.

ارزش بالینی WES در توانایی آن برای بررسی همزمان تغییرات کدکننده در سراسر ژنوم نهفته است؛ این امر اساساً به سرگردانی‌های طولانی‌مدت تشخیصی برای بسیاری از بیماران مبتلا به بیماری‌های نادر پایان می‌دهد. تولید داده‌های با توان عملیاتی بسیار بالا با خطوط لوله بیوانفورماتیکی (Bioinformatics Pipelines) پیشرفته ترکیب می‌شود تا فیلتر، حاشیه‌نویسی (Annotation) و اولویت‌بندی متغیرها (Variants) بر اساس فنوتیپ بالینی دقیق بیمار انجام شود. این رویکرد به ویژه در بیماری‌هایی با ناهمگنی بالینی و ژنتیکی بالا، مانند ناتوانی ذهنی (Intellectual Disability - ID)، تأخیر تکاملی (Developmental Delay - DD) و بیماری‌های متابولیک ارثی بسیار قدرتمند عمل می‌کند.

در آزمایشگاه ژن نگار آیندگان، چارچوب تشخیصی WES بالینی برای دستیابی به حداکثر حساسیت و ویژگی جهت تشخیص واریانت‌های پاتوژن تثبیت‌شده بهینه‌سازی شده است. این آزمایشگاه با بهره‌گیری از فناوری‌های پیشرفته در زمینه کپچر (Capture) و پلتفرم‌های توالی‌یابی، خوانش با دقت بالایی از توالی‌های اگزونی و مرزهای اینترونی مجاور آن‌ها را تضمین می‌کند. یکپارچگی عمیق بالینی با تطبیق ویژگی‌های فنوتیپی بیمار، که با استفاده از اصطلاحات استاندارد HPO (Human Phenotype Ontology) تعریف می‌شوند، با الگوریتم‌های پیشرفته اولویت‌بندی متغیرها حاصل می‌گردد و بازده تشخیصی را به طور قابل‌توجهی افزایش می‌دهد.

یکی از ویژگی‌های حیاتی WES تشخیصی، انعطاف‌پذیری در آنالیز داده‌ها به صورت انفرادی (Singleton) یا به صورت پویا در کنار اعضای خانواده، به ویژه از طریق آنالیز تریو (Trio-WES) (بیمار و والدین بیولوژیک) است. آنالیز تریو با قدرت بالایی واریانت‌های جدید (de novo) را جداسازی کرده و با کاهش چشمگیر تعداد متغیرهای با اهمیت نامشخص (Variants of Uncertain Significance - VUS)، رسیدن به تشخیص قطعی را تسریع می‌بخشد. استفاده از دستورالعمل‌های معتبر و مبتنی بر شواهد جهت تفسیر نتایج، اطمینان حاصل می‌کند که گزارش نهایی، بینشی دقیق، قابل استناد و مفید در رابطه با معماری ژنتیکی بیمار ارائه خواهد داد.

٢. مروری بر تکنولوژی (پشتوانه علمی)

توالی‌یابی نسل جدید (Next-Generation Sequencing - NGS) چشم‌انداز فناوری تشخیص‌های مولکولی را از بررسی‌های سریال و نقطه‌ای به توالی‌یابی موازی و با حجم بالا (Massively Parallel Sequencing) تغییر داده است. WES بالینی با جداسازی و غنی‌سازی (Enrichment) DNA هدف در مناطق اگزونی از طریق استراتژی‌های گیراندازی (Capture)، از این فناوری بهره می‌برد. بر خلاف توالی‌یابی کل ژنوم (Whole Genome Sequencing - WGS) که تمام ژنوم را خوانش می‌کند، WES از پروب‌های اختصاصی (Baits) از جنس RNA یا DNA برای هیبریداسیون هدفمند این نواحی استفاده می‌کند و هزینه‌های عملیاتی را به شدت کاهش می‌دهد.

ستون فقرات فنی توالی‌یابی اگزوم مدرن در درجه اول بر پایه شیمی توالی‌یابی از طریق سنتز (Sequencing by Synthesis - SBS) استوار است که به‌طور عمده در پلتفرم‌های با توان عملیاتی بالای سری Illumina NovaSeq اجرا می‌شود. در تکنیک SBS، نوکلئوتیدهای فلورسنت محدودشده به‌طور سیستماتیک در رشته‌های DNA مکمل جایگزین می‌شوند. سیستم‌های نوری با وضوح بالا در سراسر فلوسل (Flow Cell) سیگنال‌ها را ضبط کرده و میلیون‌ها خوانش کوتاه (Short Reads) دوطرفه (Paired-end)، عموماً با طول ۱۵۰ جفت‌باز، تولید می‌کنند.

روش‌های کپچر نقش بسیار مهمی در یکنواختی و کامل بودن توالی‌یابی اگزوم ایفا می‌کنند. طراحی‌های روز دنیا برای غنی‌سازی هدف (Target Enrichment) توسط ارائه‌دهندگانی مانند Agilent یا Twist به دقت مهندسی شده‌اند تا علاوه بر توالی‌های کدکننده اصلی (CCDS)، نواحی اتصال مهم (Splice-sites) و مناطق تنظیم‌کننده با اهمیت بالینی که در پایگاه‌هایی نظیر ClinVar ذکر شده‌اند را نیز شناسایی کنند. یکنواختی پوشش اطمینان می‌دهد که نواحی با درصد بالای GC که مستعد افت پوشش هستند، به درستی تکثیر و خوانده شوند.

دستیابی به عمق پوشش مطلوب (Depth of Coverage) -به‌طور معمول بیش از 100x در محیط‌های بالینی- برای فراخوانی دقیق واریانت‌ها (Variant Calling) ضروری است، خصوصاً در تشخیص موزاییسم سوماتیک با سطح پایین یا شناسایی حالت‌های هتروزیگوت. در ژن نگار آیندگان، پارامترهای با عمق خوانش مناسب با هم‌ترازی دقیق بیوانفورماتیکی با ژنوم مرجع (GRCh37/GRCh38) همراه می‌شوند. شاخص کیفی Phred بالا در کنار پوشش عمیق، امکان نقشه‌برداری دقیق متغیرهای تک نوکلئوتیدی (SNV) را فراهم می‌آورد و اعتبار علمی گزارشات نهایی را تضمین می‌سازد.

٣. اندیکاسیون‌های بالینی و معیارهای انتخاب بیمار

ادغام توالی‌یابی کل اگزوم در رویه‌های بالینی نیازمند انتخاب دقیق بیماران بر اساس دستورالعمل‌های تدوین‌شده توسط کالج ژنتیک و ژنومیک پزشکی آمریکا (ACMG) می‌باشد. WES به عنوان یک رویکرد تشخیصی خط اول، یا آزمایش خط دوم سریع برای بیمارانی که دارای ناهنجاری‌های مادرزادی متعدد (Multiple Congenital Anomalies - MCA)، تأخیر تکاملی عمومی (GDD) یا ناتوانی ذهنی سندرمیک هستند، به شدت توصیه می‌گردد. در اختلالات عصبی-تکاملی با پیشرفت سریع، درخواست فوری این پروفایل ژنومیک می‌تواند مسیر درمانی را کاملاً تغییر دهد.

خطاهای متابولیک ارثی (Inborn Errors of Metabolism - IEM)، آنسفالوپاتی‌های صرعی کودکان و سندرم‌های نقص ایمنی پیچیده نشان‌دهنده کانون‌های اصلی هدف با هم‌پوشانی فنوتیپی شدید هستند. هنگامی که غربالگری‌های بیوشیمیایی استاندارد یا پنل‌های محدود نتیجه‌ای ندارند، WES گسترده‌ترین چشم‌انداز را برای یافتن پاتولوژی تک‌ژنی (Monogenic) ارائه می‌دهد. علاوه بر این، ارزیابی بیمارانی با فنوتیپ‌های متقاطع مشکوک به تشخیص‌های مولکولی دوگانه بسیار کارآمد خواهد بود.

صلاحیت کاندیداها به نسبت هزینه-فایده و جلوگیری از سرگردانی در راهروهای پاراکلینیکی (Diagnostic Odyssey) بستگی دارد. برای افرادی که سال‌ها درگیر بررسی‌های پراکنده و بی‌نتیجه بوده‌اند، اگزوم یک پایگاه نهایی در نظر گرفته می‌شود. با این حال، ثبت دقیق فنوتیپ و تکمیل پارامترهای معاینه بالینی (استفاده از HPO terms) شامل داده‌های نورولوژیک، دیس‌مورفیک، رادیولوژیک و بیوشیمیایی برای عملکرد صحیح تیم تفسیر کننده به شدت الزامی است.

در دسترس بودن اعضای بیولوژیک خانواده به خصوص والدین تأثیر عمیقی در بازده این روش دارد. بررسی خون پدر، مادر و فرزند بیمار که تحت عنوان Trio-WES شناخته می‌شود، همواره به عنوان استاندارد طلایی (Gold Standard) در تفسیر مسیرهای توارث معرفی شده است. اگرچه آزمایش انفرادی (Singleton WES) نیز پذیرفته بوده و غالباً انجام می‌شود، فیلتر کردن تغییرات خانوادگی بی‌خطر بستگی زیادی به بررسی همزمان ژنوم والدین دارد.

• ناتوانی ذهنی متوسط تا شدید (ID) یا تأخیر تکاملی عمومی (GDD).
• ناهنجاری‌های مادرزادی متعدد (MCA) و ویژيگی‌های واضح دیس‌مورفیک (شکل ظاهری).
• آنسفالوپاتی‌های صرعی با شروع زودرس و اختلالات تشنج مقاوم به درمان.
• شک به الگوهای نامنظم اختلالات متابولیک ارثی (IEM) و میتوکندریایی.
• طیف اوتیسم سندرمیک (Syndromic ASD) همراه با پس‌رفت چشمگیر عصبی-تکاملی.
• پایان دادن به سرگردانی‌های طولانی تشخیصی (Diagnostic Odyssey).

٤. شرایط نمونه‌گیری و ملاحظات پیش‌تحلیلی

مرحله پیش‌تحلیلی (Pre-analytical) عامل تعیین‌کننده‌ای جهت تضمین صحت داده‌های NGS است که با استخراج اسیدهای نوکلئیک با کیفیت آغاز می‌شود. خون محیطی وریدی در لوله‌های استاندارد حاوی EDTA به عنوان نمونه مرجع (Gold Standard) برای استخراج DNA ژنومی باقی می‌ماند. برای اطفال، نوزادان، یا در شرایطی که خون‌گیری ناممکن است، کیت‌های جمع‌آوری بزاق یا سواب‌های دهانی حاوی بافرهای نگه‌دارنده جایگزین‌های مناسبی هستند. حفظ زنجیره سرد و کنترل مناسب شرایط انتقال از افت کیفیت نمونه‌ها جلوگیری می‌کند.

پس از مرحله پذیرش، DNA ژنومی پیش از آماده‌سازی کتابخانه (Library) تحت ارزیابی کمی و کیفی دقیقی قرار می‌گیرد. روش‌های فلورومتریک (نظیر سیستم‌های Qubit) برای محاسبه دقیق غلظت مورد استفاده قرار می‌گیرند، در حالی که دستگاه‌های اسپکتروفتومتر نسبت خلوص را در طول‌موج‌های ۲۶۰/۲۸۰ و ۲۶۰/۲۳۰ ارزیابی می‌نمایند. کیفیت و عدم حضور قطعات به شدت تخریب‌شده (Degraded) یا آلودگی‌های شیمیایی نقش شایانی در موفقیت فرایند هیبریداسیون‌های بعدی و بهره‌وری کپچر خواهد داشت.

علاوه بر الزامات بیوشیمیایی، WES نیازمند جمع‌آوری دقیق داده‌های انفورماتیکی پیش‌تست است. پزشک ارجاع‌دهنده موظف به ارائه مستندات جامعی است که ترجیحاً باید در قالب اصطلاحات استاندارد فنوتیپی (HPO Terms) بیان شده باشند. تجزیه و تحلیل حداقل سه نسل از خانواده در قالب شجره‌نامه (Pedigree)، بررسی وضعیت ازدواج خویشاوندی (Consanguinity) و ثبت دقیق نتایج قبلی بررسی‌شده توسط (CMA, MLPA) برای کالیبره کردن الگوریتم‌های تفسیری ضروری می‌باشد.

روند اخذ رضایت‌نامه (Informed Consent) که یکی دیگر از ارکان اصلی در مرحله پیش‌تحلیلی است باید به صراحت به بررسی امکان گزارش یافته‌های ثانویه (Secondary Findings) طبق دستورالعمل‌های به‌روزرسانی شده بپردازد. ژن نگار آیندگان دریافت رضایت آگاهانه در خصوص یافته‌های تصادفی را الزامی می‌داند. مشاوره‌ ژنتیکی قبل از تست به منظور آماده‌سازی ذهن بیمار برای برخورد با واریانت‌های با اهمیت نامشخص (VUS) اکیدا توصیه می‌گردد.

٥. جریان کار آزمایشگاهی (آماده‌سازی کتابخانه، توالی‌یابی، QC)

فاز آزمایشگاهی WES با قطعه‌قطعه کردن (Fragmentation) مکانیکی یا آنزیمی DNA تصفیه‌شده آغاز می‌شود. استفاده از دستگاه‌های اولتراسونیک و فراصوتی منجر به تکه‌تکه شدن تصادفی و بدون سوگیری طول‌های DNA در بازه ۲۰۰ تا ۳۰۰ جفت‌باز می‌گردد. متعاقب آن، آنزیم‌ها اقدام به افزودن یک باز آدنین در انتهای قطعات جهت تشکیل دم (A-tailing) می‌کنند. این گام حیاتی اجازه می‌دهد تا آداپتورهای اختصاصی (Adapters) که حاوی شاخص‌های مولکولی یکتا (UMIs) و بارکدهای دوگانه هستند، به‌خوبی جهت تفکیک نمونه‌ها به قطعات متصل شوند.

کتابخانه‌های اتصال‌یافته با آداپتور، تحت پروتکل‌های کپچر هدفمند (Targeted Capture) قرار می‌گیرند تا مناطق ژنتیکی مورد علاقه غنی‌سازی شوند. قطعات مذکور در طول یک شبانه‌روز با مجموعه پیچیده‌ای از پروب‌های هیبریدی بیوتین‌دار RNA یا DNA که برای کل اگزوم انسانی طراحی شده‌اند، مواجه می‌شوند. دانه‌های مغناطیسی پوشیده‌شده با استرپتاویدین، مجموعه‌های هدف را جدا کرده در حالی که فرآیندهای شستشوی شدید، به شستن قطعات بدون هدف و غیر ضروری غیر اگزونی می‌پردازند. سپس با یک مرحله PCR، محصول غنی‌سازی‌شده و نهایی ساخته می‌شود.

پیش از بارگذاری در فلوسل، کتابخانه‌های غنی‌شده، کنترل کیفی نهایی چندگانه‌ای (QC) را طی می‌کنند. دستگاه‌های آنالیزگر قطعات متصل‌نشده را جدا کرده و اندازه صحیح قطعات (Fragment Sizing) را رصد می‌نمایند، سپس با استفاده از تکنولوژی qPCR، مولاریته قطعی نمونه‌های آماده‌شده ارزیابی می‌گردد. با ادغام ده‌ها نمونه درون استخرهای ترکیبی هم‌مولار (Equimolar Pooling)، توزیعی متناسب برای پردازش در ماشین‌های قدرتمندی نظیر NovaSeq شکل می‌گیرد.

توالی‌یابی روی دستگاه از طریق ایجاد کلاسترهای میلیونی با روش تقطیر پل (Bridge Amplification) و سنتز مستقیم انجام می‌پذیرد. فاکتورهای اولیه کیفیت دستگاه به بررسی درصد کلاسترهای تایید شده (%PF)، شاخص کیفی پایه Phred (بالاتر از %۸۵ با Q30) و نسبت‌های جداسازی می‌پردازد. تخصیص بالای پوشش منطقه خوانش نسبت به هدف و مقادیر یکنواختی، نشانگر اجرای بی‌نقص شیمیایی در این مرحله بوده و سلامت پایه فرایند بالینی ژن نگار آیندگان را بازگو می‌کند.

٦. جریان کار بیوانفورماتیکی (هم‌ترازی، فراخوانی واریانت، حاشیه‌نویسی)

ترجمه داده‌های خام اپتیکال NGS به فرمت‌های ژنتیکی قابل تفسیر توسط مجموعه‌ای غول‌پیکر و پیوسته و چندمرحله‌ای از خطوط لوله بیوانفورماتیک (Bioinformatics Pipelines) صورت می‌پذیرد. در ابتدا، داده‌های پردازش تصویر به فایل‌های توالی در قالب FASTQ بازنویسی و تبدیل می‌گردند. سپس هم‌ترازی قطعات کوتاه (Reads) بر روی ژنوم مرجع (مانند GRCh37/GRCh38) به وسیله ابزار دقیق Burrows-Wheeler Aligner (مانند BWA-MEM) انجام می‌شود. پس از تهیه نقشه هم‌ترازی اصلی در فایل‌های BAM، خوانش‌های تکراری حاصل از خطای PCR توسط ابزاری مانند Picard علامت‌گذاری و حذف می‌گردند.

برای اطمینان از قرارگیری کدهای بااهمیت، کالیبراسیون مجدد امتیاز کیفیت باز (Base Quality Score Recalibration - BQSR) وارد عمل می‌شود تا خطاهای سیستماتیکی که در طول خوانش سنتزی بوجود آمده‌اند را از لحاظ ریاضیاتی اصلاح کند. در ادامه، نرم‌افزارهای بسیار دقیق جستجوی واریانت که با ساختار سلول‌های زایا کالیبره شده‌اند (نظیر برنامه‌های پکیج قدرتمند GATK HaplotypeCaller)، وارد صحنه می‌شوند تا انواع تفاوت‌های موجود از جمله تغییر تک‌نوکلئوتیدی (SNV) و جهش‌های کوچکتر درج/حذف (Indel) را شناسایی کرده و فایل جامع (VCF) را ایجاد نمایند.

پس از تولید فایل پایه VCF، حاشیه‌نویسی محاسباتی و ادغام اطلاعات (Annotation Framework) تغییرات ژنومی خشک و سخت را به حقایق بیولوژیکی مرتبط تفسیر می‌نماید. نرم‌افزارهای این حوزه با اتصال به پایگاه‌های اطلاعاتی نظیر gnomAD و 1000 Genomes Project، فراوانی جمعیتی هر آلل مدنظر (MAF) را استخراج می‌کنند. همچنین دسته‌بندی‌های قبلی ثبت‌شده در سیستم‌های ClinVar و HGMD به محتویات اضافه می‌شوند. ابزارهای درون-رایانه‌ای که به اصطلاح in silico نامیده می‌گردند (مانند CADD ،REVEL و SpliceAI)، جهت محاسبه تخمین‌های صدمه و بیماری‌زایی، داده‌های خود را الحاق می‌سازند.

به طور تکمیلی در فاز بیوانفورماتیک، تجزیه و تحلیل عمق خوانش (Read Depth Analysis) امکان یافتن تغییرات تعداد کپی (CNVs) در سرتاسر اگزوم را توسط تکنیک‌های هم‌ترازی مانند ExomeDepth یا ابزارهای GATK-gCNV در دسترس قرار می‌دهد. محاسبات انطباقی، افت یا خیزش دقیق خوانش هر قسمت را با بانک اطلاعات مرجع بررسی می‌نمایند. نتیجه نهایی از این اسمبلینگ پیچیده الگوریتمی، یک شبکه تفکیک شده، توصیف شده و منسجم از هزاران یافته جهت شروع گام نهایی ارزیابی در بخش ژنتیک پزشکی و تشخیصی است.

٧. تفسیر و گزارش‌دهی (طبقه‌بندی ACMG/AMP)

حیاتی‌ترین و تخصصی‌ترین مرحله WES بالینی، تفسیر واریانت‌ها است؛ فرایندی طاقت‌فرسا که سنتز داده‌های محاسباتی با تخصص دقیق انسانی را می‌طلبد. آزمایشگاه ژن نگار آیندگان این مرحله حیاتی را به طور سیستماتیک و منطبق بر دستورالعمل‌های توافقی استاندارد منتشر شده توسط کالج علوم ژنتیک و ژنومیک پزشکی آمریکا (ACMG) و انجمن آسیب‌شناسی مولکولی (AMP) اجرا می‌نماید. این فریم‌ورک از روش چندوجهی بررسی‌های عملکردی، مدل‌های آسیب‌پذیری زیستی و داده‌های پراکندگی خانوادگی بهره‌برداری می‌کند.

بر اساس معماری تعیین شده توسط ACMG/AMP سال ۲۰۱۵، واریانت‌های استخراج شده به طور قطعی در پنج کلاس اصلی سطح‌بندی می‌شوند: پاتوژن (P)، محتمل بر پاتوژن (LP)، واریانت‌های با اهمیت نامشخص (VUS)، محتمل بر بی‌خطر (LB) و بی‌خطر (B). این پروسه عمیقاً از ترم‌های اختصاری HPO که ویژگی‌های بالینی بیمار را شرح داده‌اند برای فیلترینگ استفاده می‌کند؛ همچنین ابزارهایی مانند Exomiser وظیفه مقایسه با فنوتیپ‌های پایگاه داده را دارند. تأکید اولیه همیشه بر روی جهش‌های پاتوژن تطابق یافته با سیمای بیماری مراجعه‌کننده می‌باشد.

تشخیص همه‌گیر و پرتکرارِ متغیرهای با اهمیت نامشخص (VUS) همواره بزرگترین مانع و چالش در مسیر ارائه بیماری‌های مدل مندلی محسوب می‌شود. یک VUS نماینده اختلالی در عملکرد است اما مدارک متقن برای رسیدن به حد نصاب دسته‌بندی قطعی به عنوان پاتوژن را تا این لحظه ندارد. رفع این ابهامات نیازمند تحلیل عمیق پراکندگی جهش در خانواده (Segregation)، مطالعات بافت‌های RNA یا برنامه بازتحلیل سیستماتیک (Reanalysis) نتایج پس از بازه‌های طولانی خواهد بود که بانک اطلاعاتی رشد کرده باشد.

در گزارش بالینی نهایی، کشفیات ژنومیک در قالب کاربردی درمی‌آیند؛ تا روشن‌ترین توجیهات جهت پزشک ارائه شود. گزارش‌های مثبت علاوه بر نام‌گذاری رسمی منطبق بر HGVS شامل نوع زیگوسیتی، میزان خوانش دقیق و پیشینه مقالات می‌باشند. همچنین در صورتی که رضایت کتبی وجود داشته باشد، گزارش یافته‌های ثانویه (مطابق با آپدیت‌های لیست ACMG v3.2) حول محور حساسیت‌های توارثی قلبی، عروقی و انکولوژی فارغ از دلیل اصلی ارجاعِ بیمار ارائه خواهند شد.

٨. ژن‌ها و نواحی که بیشترین فراوانی آنالیز را دارند

بر خلاف پنل‌های تشخیصی محدود، WES تمام نقشه اگزوم را اسکن نموده و بالغ بر ۲۰ هزار ژن متمایز را به طور موازی بررسی می‌نماید. با این وجود، بعضی از حوزه‌های بالینی دارای نرخ موفقیت و مراجعه بسیار بالاتری هستند. شبکه اختلالات تکاملی بخش عظیمی از مراجعات را تشکیل می‌دهد که دربرگیرنده ژن‌های متصل به ناتوانی ذهنی زودرس، مالفورماسیون‌های غشای کورتکس مغز و فنوتیپ‌های مرتبط با اختلالات طیف اوتیسم می‌باشد.

بررسی و کندوکاو در ژن‌های شناخته شده مقصر در بیماری‌های متابولیک ارثی از دیگر بخش‌های موفق تکنیک WES بشمار می‌آیند. خطاهای متابولیکی (IEM) اکثراً در دوران شیرخوارگی علایمی بسیار پراکنده بروز می‌دهند و دستیابی به تشخیص سنتی در آنها به غایت پیچیده است. تست سیستمیک توانایی گنجاندن بررسی‌های مرتبط با مشکلات لیزوزومال، ژن‌های متصل به میتوکندری و سایر ناهنجاری‌ها را همزمان دارد که منجر به استفاده موثر از فرمول‌های رژیم و جایگزین‌های آنزیمی خواهد شد.

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

همواره باید در نظر داشت که به رغم قدرت پوشش بالای WES، محدودیت‌های ذاتی خاصی باعث می‌شود نتوان آنالیز بی نقصی روی تمامی ژن‌ها پیاده‌سازی کرد؛ برخی ژن‌ها نیاز به کاوش‌های پیشرفته دارند تا اثر پنهان‌سازی توسط پسودوژن‌ها (Pseudogenes) یا اشکالات خوانش مناطق سنگین از پیوندهای GC مرتفع شود. جدول زیر شامل لیستی با بیشترین تقاضا برای آنالیز در محیط‌های آزمایشگاهی ژنتیکی است.

٩. نقاط قوت WES

برجسته‌ترین نقطه مثبت و امتیاز ممتاز WES مقیاس آنالیزهای فوق‌گسترده در کنار عمق خوانش متمرکز است. در حالی که پنل‌های ایزوله تنها نقاطی مشخص از ژنوم را پوشش قرار می‌دهند، پایش حدود ۲۰,۰۰۰ ژن، چارچوبی طلایی خصوصاً برای عارضه‌های چندسیستمی ناشناخته ارائه می‌سازد. بدین ترتیب بازده تشخیصی می‌تواند در آنالیزهای کامل خانوادگی (Trio-WES) در مجموعه‌های تأیید شده، تا مرز ۴۰٪ نیز پیشتازی کرده و موقعیت این آزمایش را از لحاظ اقتصادی و دقت علمی بی‌رقیب سازد.

توان بالای سیستم در امکان یافتن لوسای (Loci) و مناطق جهش‌یافته بیماری‌زای جدید نشان می‌دهد که اگزوم نه تنها راه‌حلی برای بخش درمان بالینی بوده بلکه موتور قدرتمندی برای محققین و پزشکی ژنومیک (Translational Genomic) است. وقتی با تست‌های رایج نتایج منفی حاصل می‌گردند، WES با بررسی همزمان پدر و مادر و فرزند دارای ظرفیت بالایی در یافتن متغیرهای دی‌نوو (De novo) بوده و توان تشخیص بالایی برای رخدادهای لحظه‌ای ثبت‌نشده در نسل‌های قبلی ایجاد می‌نماید.

از سایر ارزش‌های راهبردی WES، امکان بازتفسیر طولی بدون نیاز مجدد به فرآیندهای آزاردهنده خون‌گیری است. آرشیو فایل‌ها (نظیر FASTQ/BAM) بخش کاملی از رخداد‌ها را در خود محبوس کرده است. بر پایه مطالعات، بازتحلیل سیستماتیک (Reanalysis) و فیلترگذاری مجدد به طور سالیانه و انطباق داده‌ها با پیشرفت‌های جهانی پایگاه داده (مثل آپدیت‌های ClinVar)، قابلیت افزایش ۱۰ تا ۱۵ درصدی در بازده تشخیص نهایی طی بازه‌ای یک تا دو ساله خواهد داشت.

صرفه‌جویی همه‌جانبه اقتصادی نسبت به تکنیک‌های انزواگرا و همچنین جستجوی یافته‌های ثانویه (Secondary Findings) نیز نقاط قوتی کم‌نظیرند. همراه کردن تحلیل ژنتیکِ بیماری با بررسی اتفاقی فنوتیپ‌های سرطانی پنهان یا رخدادهای قلبی (بر اساس تأییدیه‌های مقالات در لیست ACMG SF) باعث قدرت بخشیدن به بیماران شده تا تصمیماتی آگاهانه برای پیشگیری‌های دوره‌ای مرتبط اتخاذ نمایند.

١٠. محدودیت‌ها

علی‌رغم تمام توانمندی‌های گسترده، WES استانداردِ بالینی محدودیت‌های بنیادینی ناشی از چارچوب مکانیکال گیراندازی (Capture) دارد. روش‌های اگزوم به صراحت مناطق پردازشی خود را در اگزون‌ها و لبه‌های اتصالی اطراف نگه می‌دارند. این عمل باعث می‌گردد حجم بالایی از نواحی بی‌کران میان‌ژنی، پروموترهای فاصله دار و اینترون‌های عمیق که همگی نقش‌های اساسی در بیان فرستنده‌های mRNA بازی می‌کنند رها شده و آنالیز نگردند.

فروریزش پتانسیل تحلیلی در ارزیابی واریانت‌های ساختاری بزرگ (Structural Variants - SVs) و نوسانات کپی نامبرهای عظیم از خلأهای تکنیکال دیگر می‌باشد. اگرچه الگوریتم‌های عمق‌سنجی خوانش تا حد قابل اعتمادی محاسبات را پیش می‌برند، ولی دستیابی به اطلاعات دقیق نقاط تقاطع شکست در وارونگی‌های کروموزومی در توالی‌یابی خوانش کوتاه (Short-read Sequence) همواره با نقص‌هایی روبروست؛ لذا تکنولوژی پایداری نظیر CMA غالبا پیش‌نیاز بررسی قبل کار قرار می‌گیرد.

محدودیت ذاتی در حوزه ردیابی گسترش تکرارهای پشت سرهم (Repeat Expansions) و تمایز پسودوژن‌ها نیز کاملاً بارز است. بیماری‌های مهلکی که در آن‌ها الگوهای سه‌تایی تکراری به صورت ناپایدار پدیدار می‌شوند —نظیر بیماری هانتینگتون یا نوع شایع سندرم ایکس شکننده (Fragile X)— با این روند پوشش قطعی نداشته و نیازمند روشی نظیر Repeat-primed PCR می‌باشند. همچنین فرار از کدهایی با تشابه بالای هومولوژی (بمانند SMN1/SMN2) همچنان در پردازش‌های پایه دردسر آفرین است.

در پایان باید به ایجاد تداخل و خطا در شناسایی سوماتیک موزاییسم‌ها و یا جهش‌های هتروپلاسمی میتوکندری نیز اشاره شود. چنانچه عمق خوانش برای تک‌تک بلوک‌ها در محدوده‌های بین 500x الی 1000x کالیبره نباشند، امکان فرق‌گذاری میان خطای دستگاه و وجود جهش‌های موزاییک زیگوتیک ضعیف شده خواهد بود. با این وضع معمولاً بدون اضافه‌سازی پروب‌های غنی شده اختصاصی، امکان پوشش عالی میتوکندریایی در روش پایه میسر نخواهد بود.

• ناتوانی در ارزیابی مناطق ژنومی غیرکدکننده (Non-coding) و عمیق اینترونی.
• رزولوشن ضعیف تصویرسازی برای ترانس‌لوکاسیون‌های بالانس یا وارونگی‌های (Inversions) بزرگ.
• عدم پوشش اختلالات مبتنی بر بسط تکرارهای نوکلئوتیدی (همانند هانتینگتون یا Fragile X).
• چالش‌های جدی نقشه برداری در مناطق سنگین با مقادیر فوق العاده بالای نواحی GC/AT.
• خطای تداخل سیگنالی بین توالی پسودوژن‌ها و توالی ژن استاندارد اصلی (نظیر ژن‌های SMN1/2).
• دقت پایین‌تر در شناسایی موزاییسم (Mosaicism) سوماتیکِ با درصد پایین بدون بررسی در عمق‌های بسیار بالا.

١١. مقایسه با روش‌های جایگزین

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

توالی‌یابی برنامه‌ریزی‌شده و پنل‌های NGS (Targeted Gene Panels) در مناطقی حکم‌رانی می‌کنند که فنوتیپ‌های مشخص منطبق بر لیست ژن‌های کوچکتر وجود دارد. پنل‌ها با ایجاد اعماق خوانش خیره‌کننده بر روی بخشی محدود از ژنوم، WES را در کشف مناطق موزاییسم با وضوح شگفت‌انگیز در سایه قرار می‌دهند و با کاستن از پیچیدگی‌های تحلیل انفورماتیکی مقرون‌به‌صرفه‌تر هستند. البته این امر مشروط است به واضح بودن بیماری؛ هر تلاقی کوچکی میان فنوتیپ‌ها به فروپاشی ارزش تشخیصی پنل‌ها دامن می‌زند.

در سمت دیگر، توالی‌یابی کل ژنوم (Whole Genome Sequencing - WGS) بدون توسل بر روش هیبریدِ پروب‌ها به طور پیوسته توپولوژی توالی را اسکن می‌کند. WGS از چالش مناطق غنی GC سربلند بیرون آمده، ساختارهای پیچیده‌ی کروموزومی را کشف می‌نماید و اجزای دوردست پروموتور را رمزگشایی می‌سازد. با این حال، محدودیت‌های شدید مالی و مشکلات فنی در استخراج حقایق، غالباً آن را تبدیل به پناهگاهی برای روزهای سخت پس از ناکامی کامل WES قرار داده است.

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

١٢. سناریوهای بالینی / مثال‌های تشخیصی

ترجمه فعالیت‌های آزمایشگاهی در سناریوهای تشخیصیِ ویژه اطفال، توانایی سیستم را کاملاً برجسته می‌سازد. تصور کنید زوجین دارای نسبت فامیلی نزدیک، کودکی با هیپوتونیای شدیدِ عضلانی و آنسفالوپاتی ناشناخته دارند. بررسی‌های سوخت‌وسازی در ابتدا هیچ نقصی پیدا نمی‌کنند. ساختار پردازشی WES منظره تغییرات ژنتیکی را با دقت فیلتر کرده و تطابق می‌بخشد، و در نهایت با استخراج مناطقی با افت هتروزیگوسیتی (AOH) یک متغیر مغلوب بسیار نادر واقع در لوسای مرتبط با دیستروفی عضلانی ناشناخته پیدا می‌نماید.

نمونه مهم دیگر نمایانگر وقوع پاتولوژی حاد به‌صورت منزوی از الگوهای نسل پیشین می‌باشد. یک کودک نوپا تغییراتی سندرمیک در ترکیب با طیف اوتیسم و ناهنجاری اسکلتی مختصری را نمایش می‌دهد، درحالی که نقشه CMA تماما نرمال است. پیاده‌سازی فرمول Trio-WES و تمرکز روی داده‌های مثلث والدین و فرزند، یک نقطه متغیر جهش‌یافته غالب ایزوله‌شده (de novo mutation) و مرتبط با کمپلکس سرکوب‌کننده رونویسی عصبی مشخص می‌سازد؛ این یعنی پایان دادن به سالها تشخیصِ پرهزینه.

علاوه بر تشخیص‌گذاری پایه، ارزش‌های پنهان WES در بازخوانی‌ها نهفته‌اند. طی ارزیابی انفرادی بیماری با تاخیر عصبی مزمن، تغییر ساختاری مشخص در قالب متغیر با اهمیت نامشخص (VUS) به‌دلیل فقر مستندات جهانی، بایگانی می‌گردد. پس از طی گذشت ۲۴ ماه با بازنگری الگوریتمی روی پایگاه داده‌های ارتقایافته‌ی جهانی، وضعیت واریانت مذکور بازنگری گردیده و ارتقاء دسته‌بندی به حالت پاتوژنیکِ قطعی رخ می‌دهد که درهای برنامه‌ریزی اختصاصی توان‌بخشی را بازگشایی می‌نماید.

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

١٣. تضمین کیفیت، اعتباربخشی و زمان پاسخ‌دهی

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

حفاظت از درستی تحلیل‌ها به واسطه استقرار سیستم نظارتِ چندمرحله‌ای بر مسیر ارزیابی تفسیر‌ واریانت‌ها امکان‌پذیر است. تمامی رکوردهای محاسباتی در روند مستقل و بدون مداخله، طبق لیست ارتقا یافته در پایگاه داده‌های کلین‌ژن (ClinGen) مورد بازنگری توسط متخصصین قرار گرفته و آستانه تحمل خطا کاهش می‌یابد. پایش دقیق و دوره‌ای با مدل‌های خارجی امکان برابری کیفیت در بررسی حساسیت‌های تحلیلی را افزایش می‌دهد.

معماری کلان ساختار به پایه‌های نرم‌افزارهای یکپارچه محدود است که سبب گشته مدیریت همگام‌ساز داده‌های gnomAD / ClinVar طی چرخه اجرای تحلیل حفظ گردد. همچنین، تضمین کیفیت بر قابلیت تکرارپذیری فرآیندها در دراز مدت و رهگیری مجدد استوار بوده و فضای حفاظت محرمانه در سرورها مبتنی بر مستندات استتار شده پیش می‌رود.

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

• زمان پاسخ‌دهی (Turnaround Time): در حال تکمیل
• مدارک مورد نیاز جهت پذیرش: در حال تکمیل
• نظارت بر کنترل کیفی فرآیندهای هم‌ترازی و توالی‌یابی به صورت مستمر پیگیری می‌گردد.