$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
تمت الموافقة على هذا البروتوكول من قبل لجنة الأخلاقيات في كلية الطب بجامعة إيجي. تم جمع جميع عينات الدم الطرفية بعد الحصول على موافقة خطية مستنيرة من المرضى أو أوصيائهم القانونيين، وفقا لإعلان هلسنكي. تم تصميم هذه الدراسة ذات المركز الواحد بأثر رجعي لتقييم تكرار المتغيرات المرتبطة بالاعتلال التداخلي بين المرضى الذين يخضعون لتحليل جيني لمرض يشتبه في الالتهاب الذاتي ولتحقيق علاقتها بالتهاب الأوعية الدموية. بين عامي 2022 و2025، تم تضمين عينات دم محيطية أرسلت إلى مختبر الطب الجزيئي بجامعة إيجي مع تشخيص أولي لمرض التهاب ذاتي. خضعت ما مجموعه 1,204 عينات لتحليل جيني. تم استخدام CES/WES كنهج تسلسل من الجيل القادم (NGS) يستهدف الإكسونات والتسلسلات الإنترونيكية المجاورة. تم تصنيف المتغيرات باستخدام معايير ACMG وتحليل المعلوماتية الحيوية. تم استبعاد المتحورات الحميدة والمحتملة الحميدة. من بين هذه الدفعة، تم تضمين 132 حالة أطفال تم تحديد متغيرات مرتبطة بالاعتلال التداخلي واستوفين معايير الإدراج في المجموعة النهائية للدراسة.
ضمن هذه المجموعة من طب الأطفال، تم استرجاع البيانات السريرية من قاعدة بيانات عيادات روماتيزم. بشكل عام، تم تحديد 92 متغيرا جينيا فريدا بين 132 مريضا، بما في ذلك 13 متحورا ممرضا أو محتملا تم الإبلاغ عنها سابقا وتم توثيقها في قواعد بيانات عامة للمتغيرات (مثل ClinVar) و79 متحورا جديدا غائبا عن قواعد البيانات العامة حتى فبراير 2026. كانت معظم المتحورات متغاير الشكل ومرتبطة بأنماط شبيهة ب CAPS أو بعروض التهاب الأوعية الدموية الذاتي المعقدة. استخدم WES نهج NGS مستهدفا المناطق المشفرة للبروتين (الإكسونات) والتسلسلات الإنترونيكية المجاورة للكشف عن متغيرات موقع التوصيل والتنظيم. على الرغم من أن WES يغطي حوالي 1٪–2٪ من الجينوم البشري، إلا أنه يلتقط ما يقرب من 85٪ من الطفرات المعروفة المسببة للأمراض، مما يجعله أداة فعالة للغاية للكشف عن المتحورات النادرة.
شمل سير العمل المخبري استخراج الحمض النووي الجيني من الدم المحيطي، وتجزئة الحمض النووي، وتحضير المكتبة المرتبطة بالمكيلات، وتخصيب الإكسوم باستخدام مجسات التقاط تعتمد على التهجين عبر سحب البيوتين-ستريبتافيدين. تم إجراء تسلسل مزدوج عالي الإنتاجية على منصة تسلسل كرات الحمض النووي النانوية، محققا عمق تغطية متوسط يتراوح بين 100–200x لضمان اكتشاف المتحورات بشكل موثوق. شمل التحليل المعلوماتي الحيوية محاذاة الجينوم المرجعي، واستدعاء المتغيرات، والتعليق باستخدام قواعد بيانات متغيرة مختاراة، مع إجراء تصنيف المتغيرات وفقا لإرشادات الكلية الأمريكية لعلم الوراثة الطبية والجينوميات (ACMG). تم تصنيف المتغيرات كمتغيرة تم الإبلاغ عنها سابقا أو جديدة. تم التحقق من صحة المتغيرات المختارة بواسطة تسلسل سانجر.
بالإضافة إلى خطوط الأنابيب المعلوماتية الحيوية القياسية، تم تطبيق مبادئ معالجة الإشارات القائمة على الفيزياء والنمذجة الحاسوبية المبسطة لتقييم عمق التسلسل، وتوحيد التغطية، وخصائص الإشارة إلى الضوضاء، مما يوفر طبقة استكشافية وتكميلية لتحليل عمق التسلسل الخام والتغطية الإقليمية. علاوة على ذلك، تم استخدام مفاهيم النمذجة الفيزيائية الحيوية لتقييم التأثير التراكمي للمتغيرات الجينية المرتبطة بالتهاب الأوعية الدموية.
يعمل هذا الإطار متعدد التخصصات كنموذج مختبري أولي لدراسة خصائص تصفية المتغيرات في سياق بحثي استكشافي، لا سيما لتوجيه استراتيجيات علاجية مضادة للإنترفيرون. تشمل القيود انخفاض الحساسية للمتغيرات الفسيفسائية منخفضة المستوى والحاجة إلى التحقق الوظيفي من النتائج الجديدة. تشمل التطبيقات المحتملة التشخيص المبكر لالتهاب الأوعية الدموية لدى الأطفال، واتخاذ قرارات علاجية شخصية، وتوسيع قواعد بيانات المتحورات ذات الصلة بالمناعة والروماتيزم والطب الوعائي.
تجزئة الحمض النووي
تم نقل عشرين ميكرولتر من الحمض النووي المخفف إلى أنابيب PCR معقمة جديدة. إلى كل أنبوب، أضيف 2 ميكرولتر من مخزن التجزئة/الأدينيل و3 ميكرولتر من خليط إنزيم التجزئة/الأدينيل؛ كان حجم التفاعل الكلي 25 ميكرولتر. تم خلط الخليط برفق، وتم طرد مركزه لفترة وجيزة، وتم تشغيل برنامج التجزئة على جهاز الدورة الحرارية. فور الانتهاء، تم تطرد المركزي للأنابيب لفترة وجيزة ووضعها على الثلج. كانت النتيجة المتوقعة حجم متوسط لجزء DNA يتراوح بين 200–300 قطعة حراسة.
ربط المحول
إلى كل منتج تجزئة، أضيف 2.5 ميكرولتر من محول التسلسل، يليه 10 ميكرولتر من خليط الربط الرئيسي (بدون دوامة). تم تفريغ الخليط بلطف لضمان التجانس، وتم تفريغه لفترة وجيزة، وحضنه عند 20 درجة مئوية لمدة 15 دقيقة مع فتح غطاء جهاز الدورة الحرارية. بعد الحضانة، تم تجميع الأنابيب لفترة وجيزة وتركها على الثلج.
التنقية القائمة على الخرزات
تمت إضافة ثلاثين ميكرولتر من الخرزات المغناطيسية إلى كل عينة وخلطها جيدا عن طريق التوصيل حتى يتم الحصول على تعليق متجانس. تم حضن الخليط في درجة حرارة الغرفة لمدة 5 دقائق، ثم وضعت الأنابيب على حامل مغناطيسي لمدة 3 دقائق. بعد أن تنظف السائل الفائق، تم إزالته بعناية. تم غسل الحبيبة مرتين ب 100 ميكرولتر من 80٪ إيثانول، وتمت إزالة الإيثانول المتبقي بعد الغسل النهائي. ثم تجف الخرزات بالهواء على الحامل المغناطيسي لمدة تصل إلى 5 دقائق. تمت إضافة تسعة ميكرولترات من الماء الخالي من النوكليز، وأعيد تعليق الحبيبة عن طريق الأنابيب وحضنتها لمدة دقيقتين في درجة حرارة الغرفة. بعد فصل مغناطيسي لمدة 3 دقائق، تم نقل 7.5 ميكرولتر من المادة الفائقة إلى أنبوب PCR جديد معنون.
تضخيم PCR
تم تخفيف بادئات النسخ الصناعي (IVT) (20 ميكرولتر من مخزون البرايمر + 80 ميكرولتر من الماء الخالي من النوكلياز). لتفاعل PCR (الحجم الإجمالي 17.5 ميكرولتر)، تمت إضافة 2.5 ميكرولتر من برايمر IVT 1، و2.5 ميكرولتر من برايمر IVT 2، و12.5 ميكرولتر من خليط التضخيم الرئيسي للمكتبة. تم وضع الخليط بلطف على القناة، وتم تنفيذ برنامج PCR-3.
التنقية بعد PCR
تمت إضافة خمسة وعشرين ميكرولتر من الخرزات المغناطيسية إلى كل منتج PCR. بعد 5 دقائق حضانة بدرجة حرارة الغرفة وفصل مغناطيسي، تم غسل الخرزات مرتين ب 100 ميكرولتر من الإيثانول بنسبة 80٪. بعد التجفيف بالهواء لمدة تصل إلى 5 دقائق، أضيف 11 ميكرولتر من الماء الخالي من النوكليز. تم نقل عشرة ميكرولتر من المادة الفائقة إلى أنبوب جديد. تم قياس تركيز المكتبة؛ الهدف: >25 نانوغرام/ميكرولتر.
التجميع والتهجين (عينات لكل مجموعة)
تم تجميع المكتبات مع ثماني عينات مرضى في كل مجموعة (93.75 نانوغرام لكل عينة، بإجمالي 750 نانوغرام). تم تعديل الحجم إلى 12.5 ميكرولتر باستخدام مياه خالية من النوكلياز إذا لزم الأمر. تم حضان خليط التهجين مسبقا عند 65 درجة مئوية لمدة 15 دقيقة. تمت إضافة محلول الحجب، وأوليجونكليوتيدات الحجب الشامل، ومسبار التقاط الإكسوم، والماء الخالي من النوكلياز، ومحسن التهجين بشكل متسلسل. تم إطلاق برنامج التهجين لمدة 16 ساعة (18:00–10:00).
التقاط الأهداف المهجنة على خرز الستربتافيدين
تم غسل خرز الستريبتافيدين ثلاث مرات باستخدام حاجز الربط. بعد التهجين لمدة 16 ساعة، يضاف الخليط إلى الخرزات ويحتضن عند حرارة 25 درجة مئوية لمدة 30 دقيقة مع خلط لطيف كل 5 دقائق. تم الغسل باستخدام مخزن غسل الالتقاط 1 في درجة حرارة الغرفة، تلاها ثلاث غسلات باستخدام مخزن غسل الالتقاط المسبق (48 درجة مئوية) 2، بما في ذلك الحضانة عند 48 درجة مئوية. بعد الغسل النهائي، أعيد تعليق الحبيبة في 23 ميكرولتر من الماء الخالي من النوكلياز.
تحضير الحمض النووي أحادي السلسلة (ssDNA)
تمت إضافة أربعة وعشرين ميكرولتر من مخزن TE، وتم تغيير الخليط عند 95 درجة مئوية لمدة 3 دقائق، ثم وضع فورا على الثلج. تمت إضافة خليط رئيسي يحتوي على حاجز ربط الجبيرة وليغاز DNA سريع، وتم تشغيل برنامج SS-2 (37 درجة مئوية، 30 دقيقة) لتدوير الحمض النووي أحادي السلسلة. تمت إضافة حاجز الهضم وخليط إنزيمات الهضم، وتم تشغيل برنامج SS-3 (37 درجة مئوية، 30 دقيقة). ثم أضيف 3.75 ميكرولتر من حاجز توقف الهضم. تمت إضافة خمسة وثمانين ميكرولترا من الخرزات المغناطيسية، تلتها تنقية خرزات قياسية. تم نقل خمسة عشر ميكرولتر من المادة الفائقة إلى أنبوب جديد؛ كان التركيز المتوقع 0.8–2 نانوغرام/ميكرولتر.
تحضير كرة نانوية للحمض النووي (DNB)
تم استخدام منتجات ssDNA لتكوين DNB. تمت إضافة خلطات مخازل تكوين DNB، ومخزن TE منخفض EDTA، وخلطات إنزيم DNB 1 و2. تم تشغيل برامج DNB-1 و DNB-2 بشكل متسلسل. بعد الانتهاء، أضيف 20 ميكرولتر من مخزن توقف تفاعل DNB وخلطه بلطف (5–8x) باستخدام رؤوس ماصة عريضة السبطانة؛ كان التركيز المتوقع من 8 إلى 40 نانوغرام. تم تحميل DNBs المحضرة على منصة تسلسل كرات الحمض النووي النانوية للتسلسل عالي الإنتاجية.
معالجة بيانات التسلسل وتحليل الإشارات
تم التحكم في جودة قراءات التسلسل الخام باستخدام FastQC (v0.11.9) وfastp (v0.23.1). تم تحسين عتبة نسبة الإشارة إلى الضوضاء (SNR) البالغة 20 ديسيبل باستخدام تحليل منحنى خصائص تشغيل المستقبل (ROC) مقابل مجموعة بيانات مرجعية للمتغيرات المعروفة الالتهابية الذاتية، مع موازنة معدل الاكتشاف الكاذب (FDR) البالغ <1٪ وحساسية مستهدفة تبلغ >95٪ للمتغيرات منخفضة التردد؛ تم التخلص من القراءات التي تسقط تحت هذا العتبة ال20 ديسيبل. تم إجراء محاذاة مع الجينوم المرجعي لاتحاد الجينوم Human Build 38 (GRCh38) المرجعي واستدعاء المتغيرات الأولي باستخدام BWA-MEM (v0.7.17) ومجموعة أدوات تحليل الجينوم (GATK، v4.2.6).
بالنسبة لخطوات معالجة الإشارات القائمة على الفيزياء، تم تحويل تسلسلات النيوكليوتيدات إلى إشارات رقمية ثنائية (0 تمثل البيورينات؛ 1 تمثل البيريميدينات). لحساب DFT، تمت معالجة الإشارة الثنائية بنافذة منزلقة N = 512 زوج قواعد وتداخل بنسبة 50٪ (حجم خطوة 256 بيتاب) للحفاظ على دقة جينومية محلية. تم تعريف DFT على أنها:

لتصفية التشويشات في التسلسل عالي التردد دون الإفراط في تنعيم متغيرات النيوكليوتيد الأحادي الحقيقية (SNVs)، والتي تظهر كانتقالات محلية حادة وعالية التردد، تم معايرة مرشح رقمي منخفض التمرير برمجيا. تم تحديد تردد القطع المطبع الأمثل (fc) بشكل تكراري عن طريق مسح نطاق 0.05–0.25 دورة/قاعدة. اختارت خوارزمية التحسين fc = 0.15 دورات/قاعدة، وتعرف بأنها نقطة التحول التي يحتفظ فيها طيف قوة الإشارة ب ≥85٪ من التباين الكلي لمتغيرات التحكم المعروفة ذات الإيجابية الحقيقية مع إزالة الضوضاء التقنية الخلفية. تم تطبيق التصفية المعتمدة على DFT باستخدام سكريبتات مخصصة مكتوبة بلغة بايثون (الإصدار 3.9)، مع استخدام مكتبات NumPy (الإصدار 1.23.0) وSciPy (الإصدار 1.9.1)، لتقليل الضوضاء عالية التردد مع الحفاظ على الميزات الطيفية المرتبطة بالطفرات. تم معايرة معلمات التصفية بشكل تكراري برمجيا لمنع التلميع المفرط للإشارات النادرة المتغيرة. يتم تلخيص سير العمل المتكامل لمعالجة الإشارات والنمذجة الفيزيائية الحيوية في الشكل 2.
النمذجة الفيزيائية الحيوية والعشوائية لإشارات الإنترفيرون
للتحقيق في النتائج الوظيفية للمتغيرات المحددة، تم نمذجة ديناميكيات مسار JAK–STAT باستخدام المعادلات التفاضلية العادية (ODEs) المعرفة كما يلي:

تم تنفيذ محاكاة المعادلات التفاضلية الحتمية باستخدام COPASI (محاكي المسارات المعقد، الإصدار 4.36) وتم تأكيدها بواسطة سكريبتات بايثون مخصصة باستخدام وحدة scipy.integrate.solve_ivp. تم تحديد ظروف أولية واقعية بيولوجيا، وأجريت تحليلات حساسية على ثوابت معدل التضخيم باستخدام SALib (مكتبة تحليل الحساسية بلغة بايثون، v1.4.5). تم دمج حركيات الإشارة IFN-α التي بوساطة TYK2 لمحاكاة تضخيم STAT2 المدفوع بزيادة الوظيفة.
تم إدخال التأثيرات العشوائية باستخدام صياغة لانجفين:
dX = f(X) dt + g(X)dW
بالإضافة إلى ذلك، أجريت محاكاة عشوائية لجيليسبي باستخدام مكتبة بايثون GillesPy2 (v1.7.0) لنمذجة الانفجار النسخي وتنشيط الجين المحفز بالإنترفيرون (ISG) المحفز بالإنترفيرون الناتج عن IFN β غير متجانس. تم تنفيذ دوال التأخير الموزعة التي تمثلها صيغ نواة غاما في بايثون باستخدام التكامل العددي لمحاكاة آليات التغذية الراجعة النسخية المؤجلة. تم تنفيذ جميع خطوط الأنابيب الحسابية، بما في ذلك تصفية الإشارات والنمذجة الرياضية، على بيئة حوسبة عالية الأداء (HPC) تعتمد على لينكس.

الشكل 2. إطار متكامل لمعالجة الإشارات والنمذجة الفيزيائية الحيوية لإشارات JAK-STAT. 1: مرحلة معالجة الإشارة: رسم خريطة ثنائية لتسلسلات النيوكليوتيدات (البيورينات = 0، البيريميدينات = 1) تليها ترشيح قائم على DFT. لاحظ عتبة SNR الحرجة عند 20 ديسيبل لاستدعاء المتغيرات بدقة. 2: النمذجة الفيزيائية الحيوية: محاكاة معدلات التضخيم القائمة على المعادلات التبسيطية حيث تؤدي الاضطرابات الصغيرة في الظروف الأولية إلى تحولات عتبة في الإشارات المستمرة. 3: النمط الظاهري البيولوجي: نمذجة متغيرات كسب الوظيفة في TYK2 في مسار IFN-α، مما يؤدي إلى تفعيل STAT2 المتزايد وأنماط الاعتلال التداخلي الناتجة. الاختصارات: DFT = تحويل فورييه المتقطع؛ IFN-α = الإنترفيرون ألفا؛ JAK-STAT = محول ومحفز وفعال النسخ؛ ODE = المعادلة التفاضلية العادية؛ نسبة SNR = نسبة الإشارة إلى الضوضاء؛ STAT2 = محول الإشارة ومنشط النسخ 2؛ TYK2 = كينازين التيروسين 2. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
يلخص الشكل 2 دمج معالجة الإشارات المعتمدة على DFT مع مكونات النمذجة التعددية والمتغيرة والاحتمالية المستخدمة في سير العمل الحسابي الاستكشافي. يدمج هذا البروتوكول الشامل علم الوراثة الجزيئي المتقدم، وتقنيات التسلسل عالي الإنتاجية، وخوارزميات معالجة الإشارات، والنمذجة الفيزيائية الحيوية لدعم اكتشاف وتفسير المتغيرات المرتبطة بالاعتلال التداخلي في التهاب الأوعية الدموية الذاتي الالتهابي لدى الأطفال.
الإطار الإحصائي وتوزيع احتمالية رسم الإشارة
لتقنين رسم الخرائط الرقمية لتسلسلات الجينوم قبل ترشيح الإشارات، تم إنشاء إطار تحويل ثنائي يعتمد على الكيمياء الحيوية للنوكليوتيدات. لكل نافذة جينومية هيكلية بطول N، يتم تعيين البيورينات {A, G} إلى قيمة رقمية 0، ويتم تعيين البيريميدينات {C, T} إلى 1. وفقا للفرضية الصفرية (H0) لتوزيع جينوم خلفي غير متحيز ومنتظم، يتبع هذا التحويل الثنائي إطار تجربة برنولي. تعرف دالة كتلة الاحتمالية (PMF) للإشارة الخريطة X كما يلي
P(X = x) = px(1 - p)1-x ل x ∈ {0,1}
حيث يمثل p = 0.5 احتمال مواجهة بقايا بيريميدين عبر مسار خلفية غير محدد. عند مقياس هذا التحويل عبر قواعد متسلسلة لحساب كثافة الطيف الطاقي (PSD) عبر DFT، يتصرف توزيع الضوضاء الخلفية التراكمي كمسار عشوائي، يتقارب إلى توزيع ضوضاء بيضاء غاوسي وفقا لمبرهنة الحد المركزي. وبالتالي، فإن طيف القدرة المعدل لهذا التوزيع الصفري يتبع توزيعا كاي تربيع (χ2) بدرجتين من الحرية. للحفاظ على عتبة دلالة إحصائية صارمة (α = 0.05)، تم حساب عتبة شدة القدرة الحرجة لتعريف ارتفاع إشارة متغير ممرض حقيقي تحليليا باستخدام تكامل كثافة الاحتمالات التالي:
العتبة = - In(a) ×a 2
حيث يمثل σ2 التباين التشغيلي لأرضية الضوضاء الجينومية الخلفية المحلية. أي ذروة طيفية تتجاوز هذا الحد (p < 0.05، ما يعادل SNR > 20 ديسيبل) كانت تعطى أولوية للبوابات الترشيحية في الأسفل في السيليكو، مما يضمن أن الاستدعاء المتنوع يعتمد على البيانات وأقل تأثرا بضوضاء التسلسل العشوائي.
الكفاءة التشغيلية وتحليل التكاليف (الشكل 3)
من منظور التنفيذ الانتقالي والسريري، تم مقارنة كفاءة العمليات، ووقت الإنجاز السريري (TAT)، والجدوى الاقتصادية لهذا البروتوكول المتكامل مقابل المسارات التشخيصية التقليدية، مثل تسلسل سانجر المتسلسل أو لوحات الجينات المستهدفة المقيدة. بينما تمتد الرحلات التشخيصية التقليدية لالتهاب الأوعية الدموية للأطفال أو اعتلالات التداخلية من النوع الأول المشتبه بها غالبا من 8 إلى 12 أسبوعا بسبب الاختبارات التكرارية لجين واحد، فإن سير العمل المبسط—الذي يشمل WES عالي الإنتاجية، والتهجين المحسن لمدة 16 ساعة، وترشيح إشارات DFT المتوازي على مجموعة الحوسبة عالية الأداء (HPC)—يحقق إجمالي TAT سريري يبلغ 10–14 يوما من استلام العينة الأولي حتى التقرير الجزيئي النهائي. علاوة على ذلك، وبفضل تعدد العينات الفعال (تجميع ثماني عينات مرضى لكل كتلة هجينة)، تبلغ تكلفة الكاشف الأساسي والتسلسل حوالي 250–300 دولار لكل مريض، مقارنة باللوحات الشاملة التقليدية التي غالبا ما تتجاوز 1200–1800 دولار. يشير هذا الدمج لكل من الجدول الزمني التشخيصي وأطر التكلفة إلى أن سير عمل المعلومات الحيوية المقترح المستوحى من الفيزياء قد يكون قابلا للتنفيذ والتوسع لمختبرات الوراثة السريرية الروتينية.