تقدم هذه الدراسة مسار تحليل معلوماتي حيوي قابل للتكرار، يدمج بين تحليل النسخ (transcriptomics)، والتعلم الآلي، وتحليل الارتشاح المناعي لتحديد المؤشرات الحيوية التشخيصية المحتملة والشبكات التنظيمية في مرض الانسداد الرئوي المزمن (COPD).
يتطلب اشتراك JoVE لعرض هذا المحتوى. تسجيل الدخول أو ابدأ نسخة تجريبية مجانية
مقالة بحثية
* These authors contributed equally
تقدم هذه الدراسة مسار تحليل معلوماتي حيوي قابل للتكرار، يدمج بين تحليل النسخ (transcriptomics)، والتعلم الآلي، وتحليل الارتشاح المناعي لتحديد المؤشرات الحيوية التشخيصية المحتملة والشبكات التنظيمية في مرض الانسداد الرئوي المزمن (COPD).
يُعد مرض الانسداد الرئوي المزمن (COPD) اضطراباً تنفسياً مترقياً يتميز بتقيد مستمر في تدفق الهواء والتهاب مزمن، ومع ذلك، لا يزال دور تعديل الحمض الريبي النووي N4-acetylcytidine (ac4C) في نشوء هذا المرض غير مستكشف إلى حد كبير. هدفت هذه الدراسة إلى إجراء مسح منهجي للجينات المرتبطة بـ ac4C (ac4C-RGs) من قاعدة بيانات منشورة، وبحث شبكاتها التنظيمية في مرض الانسداد الرئوي المزمن، وبذلك تحديد المؤشرات الحيوية المحتملة لمزيد من الدراسات الميكانيكية دون افتراض وجود علاقة تنظيمية مباشرة بين أي جين محدد وتعديل ac4C. تم تحديد الجينات ذات التعبير التفاضلي (DEGs) من الملفات النسخية، وطُبق تحليل شبكة التعبير الجيني المشترك الموزون (WGCNA) للكشف عن وحدات التعبير المشترك الرئيسية. وأُجري تحليل متقاطع بين الجينات ذات التعبير التفاضلي، والوحدات الهامة، والجينات المرتبطة بـ ac4C. تم فحص الجينات الرئيسية باستخدام خوارزميات LASSO regression وXGBoost والغابات العشوائية (random forest)، تلا ذلك بناء نموذج تشخيصي يعتمد على الانحدار اللوجستي. تم تقييم أداء النموذج عن طريق تحليل منحنى خصائص التشغيل للمستقبل (ROC)، والمساحة تحت المنحنى (AUC) مع فواصل ثقة 95%، وتقييم منحنى المعايرة، وتحليل منحنى القرار (DCA). تم تحديد ما مجموعه 160 جيناً متداخلاً، واختارت خوارزميات التعلم الآلي الثلاثة باستمرار ستة جينات مركزية هي (PTRF، PRKCDBP، UPP1، TOR3A، FAM168B، وB4GALT2). أظهر النموذج التشخيصي أداءً تمييزياً جيداً، حيث بلغت قيم AUC 0.766 و0.759 و0.723 في مجموعات التدريب، والاختبار الداخلي، والتحقق الخارجي، على التوالي. وأشار تحليل الشبكة التنظيمية إلى محاور ceRNA وتفاعلات عوامل النسخ المحتملة، بينما كشف توصيف الارتشاح المناعي عن ارتباطات معنوية بين الجينات الرئيسية ومجموعات فرعية متعددة من الخلايا المناعية. وأشار تحليل التفاعل بين الدواء والجين والالتحام الجزيئي إلى أن fluorouracil وcapecitabine و5-benzylacyclouridine قد تُظهر ألفة ارتباط متوقعة مواتية مع UPP1. وفي الختام، تم تحديد PTRF وPRKCDBP وUPP1 وTOR3A وFAM168B وB4GALT2 كمؤشرات حيوية محتملة مرتبطة بـ ac4C في مرض الانسداد الرئوي المزمن، والتي قد تكون مشتركة في التنظيم المناعي والأيضي، مما يوفر أساساً للتحقيقات الوظيفية والاستكشافات العلاجية المستقبلية.
يُعد مرض الانسداد الرئوي المزمن (COPD) اضطراباً تنفسياً مزمناً وغير متجانس، يتميز بحدوث تقييد تدريجي لتدفق الهواء نتيجة تشوهات في بنية الأسناخ الرئوية والمجاري التنفسية1,2. وتشكل عدم التوازن بين البروتياز ومضادات البروتياز، والإجهاد التأكسدي، والالتهاب المزمن، والشيخوخة الخلوية الآليات المرضية الفسيولوجية الأساسية لمرض الانسداد الرئوي المزمن، مما يؤدي إلى تدمير بنيوي وضعف وظيفي في أنسجة الرئة3,4. علاوة على ذلك، يتأثر مرض الانسداد الرئوي المزمن بعوامل خطر متعددة، بما في ذلك التدخين لفترات طويلة، والتلوث البيئي، والتعرض المهني، والعدوى التنفسية، والاستعداد الوراثي5,6. وقد فرض مرض الانسداد الرئوي المزمن عبئاً كبيراً على الاقتصاد العالمي؛ حيث من المتوقع أن يمثل 0.111% من الناتج المحلي الإجمالي العالمي سنوياً من عام 2020 إلى 20507. ورغم أن الاستراتيجيات العلاجية الحالية، مثل موسعات الشعب الهوائية، والكورتيكوستيرويدات المستنشقة، وإعادة التأهيل الرئوي، والعلاج بالأكسجين طويل الأمد، يمكنها تخفيف الأعراض، إلا أن وقف تطور المرض لا يزال يشكل تحدياً. ونظراً للتباين السريري الملحوظ، فإن نتائج المرضى تختلف بشكل كبير8. وبناءً على ذلك، هناك حاجة ملحة إلى مؤشرات حيوية تشخيصية ومؤشرات إنذارية جديدة لتعزيز إدارة مرض الانسداد الرئوي المزمن وتحسين فرص بقاء المرضى على قيد الحياة.
يشير تعديل RNA إلى التغيير الكيميائي لجزيئات RNA، والذي يمكن أن يغير بنية RNA ووظيفته لتنظيم التعبير الجيني9,10. وتشمل تعديلات RNA الشائعة N6-methyladenosine (m6A)، وpseudouridine (Ψ)، و5-methylcytosine (m5C)، وac4C11. ويلعب تعديل ac4C دورًا محوريًا في الحفاظ على استقرار mRNA وتعزيز ترجمة mRNA12,13. ويُعد NAT10 الإنزيم الحقيقي النواة الوحيد المعروف الذي يحفز تعديل ac4C، ونشاطه ضروري لتكوين هذا التعديل14. وقد أظهرت الدراسات وجود ارتباط قوي بين الإجهاد التأكسدي، والشيخوخة الخلوية، والالتهاب، وتعديل ac4C. فعلى سبيل المثال، تظهر الخلايا اللمفاوية T من نوع CD4+ في أنسجة القولون لدى الأفراد المصابين بمرض الأمعاء الالتهابي (IBD) مستويات مرتفعة بشكل ملحوظ من NAT1015. ويعزز NAT10 عملية أسيتلة ac4C للكيموكينات CCL2 وCXCL1، مما يؤدي إلى تعزيز تسلل البلاعم والعدلات وتفاقم الضرر الالتهابي16. وقد تشمل الاستجابة الخلوية للإجهاد التأكسدي تعديل ac4C، كما يتضح من الزيادة الملحوظة في مستويات ac4C تحت ظروف الإجهاد التأكسدي. بالإضافة إلى ذلك، يعزز NAT10 التليف الرئوي الناجم عن PM2.5 عن طريق تثبيت mRNA الخاص بـ TGFB1 عبر تعديل ac4C، مما يؤدي إلى تحفيز الانتقال الظهاري الميزنشيمي17. ومع ذلك، لا يزال دور تعديل ac4C في مرض الانسداد الرئوي المزمن (COPD) غير مستكشف إلى حد كبير، مما يسلط الضوء على الحاجة إلى مزيد من الأبحاث في هذا المجال.
تم استخدام قاعدة بيانات GEO في هذه الدراسة لتحديد الجينات ذات التعبير المتمايز (DEGs) بين مرضى COPD والمجموعات الضابطة. ثم تم دمج قائمة منشورة من الجينات المرتبطة بـ ac4C (ac4C-RGs)، والتي جُمعت من بيانات متعددة الأوميكات18، مع الـ DEGs لتحديد المرشحين المشتركين. ونظرًا لأن تعديل ac4C معروف بتأثيره على الالتهاب والإجهاد التأكسدي - وهما عمليتان محوريتان في COPD - فقد افترضنا أن الجينات المرتبطة بشبكة تنظيم ac4C قد تكون مختلة التنظيم في COPD. ومع ذلك، لم تُعتبر الجينات المشتركة ركائز مباشرة لـ NAT10 أو جينات يتم تنظيمها مباشرة بواسطة ac4C؛ بل اعتبرت مرشحين مرتبطين بالشبكة المتعلقة بـ ac4C. تم فحص الجينات الرئيسية باستخدام خوارزميات تعلم آلي متعددة، تلا ذلك بناء نموذج تشخيصي والتحقق من صحته. لاحقًا، تم تحديد المسارات ذات الصلة المشاركة في إمراضية COPD، وتنبؤ الأدوية المستهدفة المحتملة. وأخيرًا، أُجريت فحوصات RT-qPCR لتأكيد مستويات تعبير الجينات الرئيسية، مما يوفر رؤى أولية حول الفسيولوجيا المرضية لـ COPD والسبل المحتملة للاستكشاف العلاجي.
الوصول مقيد. يرجى تسجيل الدخول أو بدء فترة تجريبية لعرض هذا المحتوى.
بيان مجلس مراجعة المؤسسة
أُجريت هذه الدراسة وفقًا لإعلان هلسنكي. وقد تمت الموافقة على البروتوكول من قبل لجنة الأخلاقيات في مستشفى شينزين لوهو للطب الصيني التقليدي (رقم الموافقة 2024-LHQZYYYXLL-KY-039)، وتم الحصول على موافقة كتابية مستنيرة من جميع المشاركين قبل الالتحاق. تتوفر تفاصيل أدوات البحث والمواد المستخدمة في هذا البروتوكول في جدول المواد.
مصدر البيانات ومعالجتها
تم الحصول على مجموعات بيانات التعبير الجيني المرتبطة بمرض الانسداد الرئوي المزمن (COPD) من قاعدة بيانات Gene Expression Omnibus (GEO). استُخدمت مجموعة البيانات GSE54837 كمجموعة بيانات النسخ (transcriptome)، بينما عملت مجموعة البيانات GSE112811 كمجموعة تحقق (الجدول 1). جُمعت الجينات المرتبطة بـ ac4C-RGs من الأدبيات العلمية18. تم تحديد الجينات ذات التعبير التفاضلي (DEGs) بين مجموعتي COPD والمجموعة الضابطة باستخدام حزمة limma في لغة R. واعتُبرت الجينات ذات التعبير التفاضلي ذات دلالة إحصائية إذا كانت |log2FC| > 0 وكانت قيمة p < 0.05. كما تم إنشاء مخططات بركانية (Volcano plots) لتصور التوزيع العام للتغيرات في التعبير الجيني.
بناء تحليل شبكة الارتباط الجيني المرجح (WGCNA)
أُجري تحليل شبكة共-تعبير الجينات المرجحة (WGCNA) على مجموعة البيانات GSE54837 باستخدام لغة R لتحديد الوحدات المرتبطة بمرض الانسداد الرئوي المزمن (COPD). وقبل بناء الشبكة، تم تحديد العينات الشاذة وإزالتها من خلال تحليل التجميع الهرمي باستخدام دالة hclust مع طريقة الربط المتوسط ومقياس المسافة الإقليدية. وتم تحديد القوة المثلى للعتبة اللينة (β تم اختيار (= 10) لتحقيق مؤشر ملاءمة طوبولوجيا خالية من المقياس R2 ≥ 0.85، وذلك للموازنة بين الطوبولوجيا الخالية من المقاييس ومتوسط الاتصال. تم بناء مصفوفة تجاور وتحويلها إلى مصفوفة تداخل طوبولوجي (TOM). وقد تم تحديد وحدات الجينات باستخدام خوارزمية القطع الشجري الديناميكي (deepSplit = 2, minClusterSize = 50). أما الوحدات ذات ارتباطات الجين الممثل (eigengene) > دمجت القيم 0.75 لاحقاً باستخدام دالة mergeCloseModules. ثم تم إيجاد الارتباط بين الجينات الذاتية للوحدات (Module eigengenes) والسمات السريرية (حالة مرض الانسداد الرئوي المزمن COPD، والعمر، والجنس، وحالة التدخين) باستخدام معاملات ارتباط بيرسون لتحديد الوحدات المرتبطة بمرض الانسداد الرئوي المزمن من أجل التحليل اللاحق.
الفحص، وتحليل الإثراء، وتحليل شبكة PPI للجينات المتداخلة
تم إنشاء مخطط فين باستخدام حزمة R ggvenn لتحديد الجينات المشتركة بين الجينات المعبر عنها تفاضلياً (DEGs)، وجينات وحدة MEsalmon، والجينات المرتبطة بـ ac4C-RGs. وأُجري تحليل الإثراء الوظيفي للجينات المشتركة باستخدام قواعد بيانات الوجود الجيني (GO) وموسوعة كيوتو للجينات والأجهزة الجينومية (KEGG) بواسطة حزمة R clusterProfiler. كما تم الحصول على معلومات التفاعل بين البروتينات (PPI) من قاعدة بيانات STRING (https://string-db.org/) لتحليل التفاعلات على مستوى البروتين بين الجينات المشتركة. واستُخدم برنامج Cytoscape لتصوير شبكة PPI الناتجة.
تحديد الجينات الرئيسية من خلال التعلم الآلي
تم تطبيق ثلاث تقنيات للتعلم الآلي: انحدار مشغل الاختيار والتقليص المطلق الأدنى (LASSO)، والتعزيز المتدرج الشديد (XGBoost)، والغابة العشوائية (RF). تم تنفيذ انحدار LASSO باستخدام حزمة glmnet مع التحقق المتقاطع بمقدار 10 طيات لتحديد المعلمة الجزائية المثلى λ. تم ضبط المعلمة type.measure على "deviance"، والمعلمة family على "binomial". تم اختيار λ المثلى باستخدام معيار λmin، الذي يقلل من الانحراف المتقاطع، مما أدى إلى تحديد 17 جيناً. تم تنفيذ XGBoost باستخدام حزمة xgboost مع المعلمات الفائقة التالية: nrounds = 100، وmax_depth = 6، وeta = 0.3، وsubsample = 0.8، وcolsample_bytree = 0.8، وeval_metric = "logloss". تم ترتيب أهمية الميزات وفقاً لمقياس الاكتساب (gain)، واختيار أفضل 30 جيناً. أما الغابة العشوائية فقد نُفذت باستخدام حزمة randomForest مع ntree = 200، وتم ترتيب أهمية الميزات وفقاً لمتوسط الانخفاض في معامل جيني (Gini)، واختيار أفضل 30 جيناً. ثم تم تقاطع الجينات المختارة بواسطة طرق التعلم الآلي الثلاث لتحديد الجينات الرئيسية للتحليلات اللاحقة.
بناء وتقييم نموذج الانحدار اللوجستي للتنبؤ بالمخاطر
قُسمت مجموعة بيانات GSE54837 عشوائياً إلى مجموعة تدريب (70%) ومجموعة اختبار (30%). تم بناء نموذج انحدار لوجستي على مجموعة التدريب باستخدام دالة glm من حزمة MASS، مع استخدام مستويات تعبير الجينات الرئيسية كميزات إدخال. وقُيم أداء النموذج باستخدام منحنيات ROC التي تم إنشاؤها بواسطة حزمة pROC. كما حُسبت فواصل الثقة 95% لقيمة AUC عبر 2,000 تكرار بطريقة bootstrap. وصُيغت معايرة النموذج باستخدام منحنيات المعايرة التي تم إنشاؤها بواسطة 1,000 عينة bootstrap متكررة (حزمة rms). أُجري تحليل DCA باستخدام حزمة dca لتقييم صافي الفائدة السريرية عبر مجموعة من احتمالات العتبة. وأُنشئ مخطط نوโมغرام (nomogram) باستخدام دالة nomogram من حزمة rms لتسهيل تقدير المخاطر الفردية.
كانت معادلة الانحدار كالتالي:
logit(P) = 0.5823 + 0.6010 × UPP1 - 0.6563 × PTRF + 0.3853 × B4GALT2 - 0.3972 × FAM168B + 0.1848 × PRKCDBP - 0.4787 × TOR3A. (1)
هنا، يمثل P الاحتمالية المتوقعة للإصابة بـ COPD، ويمثل كل معامل مساهمة قيمة التعبير الجيني المقابلة في لوغاريتم احتمالات الإصابة بـ COPD.
تحليل التعبير الجيني، وشبكة GeneMANIA، والشبكة التنظيمية الجزيئية
تمت مقارنة مستويات التعبير الجيني بين مجموعتي COPD والمجموعة الضابطة في مجموعة بيانات GSE54837 باستخدام اختبار ويلكوكسون لمجموع الرتب (Wilcoxon rank-sum test). تم إنشاء المخططات الصندوقية (Box plots) باستخدام حزمة ggplot2 لتصور توزيع مستويات التعبير، مع تراكب الوسيط، والمدى الربيعي (IQR)، ونقاط البيانات الفردية. واستُخدم GeneMANIA لبناء الشبكات الجينية والتنبؤ بالتفاعلات الوظيفية؛ حيث أُجري البحث باستخدام المعلمات الافتراضية: النوع = Homo sapiens، والحد الأقصى للجينات ذات الصلة = 20. تم تنزيل الشبكة الناتجة وتصورها، مع الإشارة إلى أنواع التفاعلات من خلال ألوان الحواف. كما تم بناء شبكة RNA تنافسية داخلية (ceRNA) لبحث آليات التنظيم ما بعد النسخ. تُم التنبؤ بـ miRNAs التي تستهدف الجينات الستة الرئيسية باستخدام قاعدتي بيانات مستقلتين: DIANA-microT (بدرجة ≥ 0.8) و miRanda (بدرجة ≥ 140، وطاقة ≤ −20 kcal/mol). واستُخدم تقاطع miRNAs المحددة بواسطة كلتا القاعدتين لبناء أزواج miRNA-mRNA. لاحقاً، تم التنبؤ بـ lncRNAs التي تستهدف هذه الـ miRNAs باستخدام قاعدة بيانات StarBase. ثم تم بناء وتصور شبكة تنظيمية من نوع lncRNA-miRNA-mRNA باستخدام برنامج Cytoscape. كما تم التنبؤ بعلاقات التنظيم النسخي باستخدام ChIP-X Enrichment Analysis Version 3 (ChEA3). وبالنسبة لكل جين رئيسي له عوامل نسخ (TFs) متوقعة، تم اختيار أعلى 10 عوامل نسخ ذات أعلى درجات إثراء. وفي الختام، تم بناء شبكة تنظيمية من نوع TF-target في برنامج Cytoscape.
تحليل إثراء مجموعة الجينات وتقييم تسلل الخلايا المناعية
أُجري تحليل إثراء مجموعة الجينات (GSEA) باستخدام حزمة clusterProfiler لاستقصاء الوظائف البيولوجية لكل جين رئيسي. وبالنسبة لكل جين رئيسي، قُسمت العينات إلى مجموعات ذات تعبير عالٍ ومنخفض بناءً على القيمة الوسيطة. أُجري تحليل التعبير التفاضلي بين المجموعتين باستخدام limma، ورُتبت قائمة الجينات الناتجة وفقاً للـ log₂ fold-change الموقّع. نُفذ تحليل GSEA باستخدام دالة gseGO لمصطلحات العمليات البيولوجية في GO ودالة gseKEGG لمسارات KEGG، وفق المعاملات التالية: minGSSize = 10، وmaxGSSize = 500، وpvalueCutoff = 0.05، وnPerm = 1,000. كما قُدِّرت الوفرة النسبية لـ 28 نوعاً من الخلايا المناعية باستخدام تحليل إثراء مجموعة الجينات لعينة واحدة (ssGSEA) المطبق في حزمة GSVA. تم الحصول على مصفوفة توقيع مجموعة جينات منقحة تضم الجينات العلامة لـ 28 نوعاً من الخلايا المناعية من الدراسات السابقة19. وبالنسبة لكل عينة، طُبقت دالة gsva مع method = "ssgsea"، وssgsea.norm = TRUE، وkcdf = "Gaussian". وحُسبت معاملات ارتباط Spearman بين درجات إثراء ssGSEA ومستويات التعبير للجينات الرئيسية الستة باستخدام دالة cor.test. كما عُدِّلت قيم p لاختبارات متعددة باستخدام طريقة Benjamini-Hochberg. وعُرضت مصفوفة الارتباط بصرياً في شكل خريطة حرارية باستخدام حزمة pheatmap.
التنبؤ بالأدوية، والرسو الجزيئي، وتحليل الارتباط بالأمراض
تم تحديد المركبات العلاجية المحتملة التي تستهدف الجينات الرئيسية باستخدام قاعدة بيانات DrugBank. كما تم إنشاء شبكة تفاعلات "العقار المستهدف للجين الرئيسي" في برنامج Cytoscape لتصور تفاعلات العقاقير والجينات المتوقعة. وأجريت عملية الالتحام الجزيئي باستخدام منصة CB-Dock2 لتقييم ألفة الارتباط. تم استرجاع البنية البروتينية ثلاثية الأبعاد لبروتين UPP1 البشري من بنك بيانات البروتينات (PDB ID: 7B8T). كما تم الحصول على البنى الجزيئية للعقاقير (بصيغة SMILES) من PubChem. وقد نُفذ الالتحام باستخدام محرك AutoDock Vina، وصُنفت المخرجات وفقاً لطاقة الارتباط الحرة (ΔG, بـ kcal/mol). وتم تصور معقدات الالتحام باستخدام برنامج PyMOL. كما تم استقصاء الارتباطات بين الجينات الرئيسية والأمراض البشرية المرتبطة بالتعرضات البيئية باستخدام قاعدة بيانات السموم الجينية المقارنة (CTD). حيث تم البحث عن كل جين بشكل فردي، واستُخرجت الأمراض العشرة الأكثر ارتباطاً بها وتم تصويرها باستخدام المخططات الرادارية.
بروتوكول RT-qPCR
تم جمع عينات من الدم الوريدي الطرفي من ثمانية مرضى مصابين بمرض الانسداد الرئوي المزمن (COPD) وثمانية من الأصحاء كمجموعة ضابطة في مستشفى شينزين لوهو للطب الصيني التقليدي. تم تشخيص مرض COPD وفقًا لمعايير المبادرة العالمية لمرض الانسداد الرئوي المزمن (GOLD)، والتي تُعرف بأن تكون نسبة FEV1/FVC بعد استخدام موسعات الشعب الهوائية < 0.70. شملت المجموعة الضابطة متطوعين أصحاء متطابقين في العمر والجنس، وليس لديهم تاريخ من الأمراض التنفسية، مع نتائج طبيعية في اختبارات وظائف الرئة (FEV1% المتوقعة ≥ 80% و FEV1/FVC ≥ 0.70). تظهر المعلومات الأساسية للمرضى في الجدول 2. تم استخلاص إجمالي الـ RNA من عينات دم مرضى COPD باستخدام طقم استخلاص الـ RNA من الدم. ولتخليق الـ cDNA، تم إجراء نسخ عكسي لـ 500 ng من إجمالي الـ RNA باستخدام طقم تخليق الـ cDNA مع إزالة الحمض النووي الجينومي وفقًا للبروتوكول المرفق. ثم خُفف الـ cDNA الناتج ليصل تركيزه إلى 150 ng/μL.
أُجري تفاعل RT-qPCR باستخدام مزيج رئيسي (master mix) لـ qPCR قائم على SYBR Green في نظام PCR في الوقت الحقيقي. احتوى كل تفاعل بحجم 10 μL على 5 μL من المزيج الرئيسي 2x SYBR Green، و0.5 μL من كل من البادئات الأمامية والعكسية (10 μM)، و1 μL من cDNA المخفف (15 ng/μL)، و3 μL من ماء خالٍ من النيوكلييز. كانت ظروف الدورات هي مسخ أولي عند 95 °C لمدة 5 min، تلتها 40 دورة عند 95 °C لمدة 10 s و60 °C لمدة 30 s، مع تحليل نهائي لمنحنى الانصهار من 60 °C إلى 95 °C للتحقق من خصوصية التضخيم. أُجريت جميع التفاعلات في ثلاث مكررات تقنية. استُخدم β-actin كجين مرجعي داخلي. تم التحقق من كفاءة البادئات لكل جين مستهدف باستخدام سلسلة تخفيفات للمنحنى القياسي، وتراوحت الكفاءة من 90% إلى 110%. تمت معايرة مستويات التعبير الجيني بالنسبة لـ β-actin، وحُسب التعبير النسبي باستخدام طريقة 2-ΔΔCt. وأُجريت المقارنات الإحصائية بين مجموعات COPD والمجموعات الضابطة باستخدام اختبار Mann-Whitney U.
التحليل الإحصائي
تم إنشاء التصورات الشبكية باستخدام برنامج Cytoscape، وأُجريت التحليلات الإحصائية باستخدام برنامج R. وما لم يُنص على خلاف ذلك، استُخدم اختبار Mann-Whitney U للبيانات ذات التوزيع غير الطبيعي، واختبار Student's t-test للبيانات ذات التوزيع الطبيعي لمقارنة مجموعتين. واعتُبرت القيمة p < 0.05 ذات دلالة إحصائية.
الوصول مقيد. يرجى تسجيل الدخول أو بدء فترة تجريبية لعرض هذا المحتوى.
تحديد الجينات المتداخلة، وتحليل الإثراء، وبناء شبكة التفاعلات البروتينية-البروتينية (PPI)
تَمَّ تحديد 3,371 جينًا من الجينات ذات التعبير التفاضلي (DEGs) من مجموعة بيانات GSE54837، شملت 1,675 جينًا ذا تعبير مرتفع و1,696 جينًا ذا تعبير منخفض. وتظهر أهم 10 جينات ذات أعلى مستويات زيادة ونقصان في التعبير في الشكل 1أتم إجراء التجميع الهرمي لبيانات GSE54837 (الشكل التكميلي 1-أ)، وطُبقت قوة عتبة لينة قدرها 10 لضمان طبولوجيا شبكة خالية من ...
الوصول مقيد. يرجى تسجيل الدخول أو بدء فترة تجريبية لعرض هذا المحتوى.
يعمل تعديل RNA المحفوظ ac4C، الذي يحدث بشكل أساسي في RNA المرسال (mRNA) وRNA الناقل (tRNA)، على تعزيز استقرار mRNA وكفاءة الترجمة20. ويعد NAT10 هو إنزيم RNA acetyltransferase الوحيد المعروف الذي يتوسط تعديل ac4C21. وقد أظهرت الدراسات أن NAT10 يزداد تعبيره في الخلايا الظهارية الرئوية لدى مرضى COPD. كما يؤدي تثبيط NAT10 إلى تعطيل الوظيفة الميتوكوندرية والاستجابات النسخية22. وبناءً على تحليل multiomics متكامل، حددت هذه الدراسة ستة جينات رئيسية مرتبطة ار...
الوصول مقيد. يرجى تسجيل الدخول أو بدء فترة تجريبية لعرض هذا المحتوى.
يقر المؤلفون بعدم وجود تضارب في المصالح. تم الحصول على موافقة مستنيرة من جميع الأشخاص المشاركين في الدراسة.
نود أن نعرب عن تقديرنا لمستشفى شينزين لوهو للطب الصيني التقليدي لتوفير المرافق السريرية والدعم الإداري الضروريين لهذه الدراسة. وأخيراً، نتوجه بالشكر إلى جميع المرضى والمتطوعين الأصحاء الذين شاركوا في هذا البحث؛ حيث كانت مساهمتهم لا غنى عنها في هذا العمل. تم دعم هذا العمل من قبل مشروع سانمينغ الطبي في شينزين (رقم SZZYSM202401018)، وصناديق تخصصات منطقة لوهو ذات الأولوية (رقم LX202402021)، وصناديق تخصصات منطقة لوهو ذات الأولوية (رقم LX202302064).
الوصول مقيد. يرجى تسجيل الدخول أو بدء فترة تجريبية لعرض هذا المحتوى.
| الاسم | الشركة | رقم فهرسي | التعليقات |
|---|---|---|---|
| βبادئات الأكتين (actin primers) | تسينكيه | نظرًا لعدم توفر نص مصدر، يرجى تقديم المحتوى الذي ترغب في ترجمته لضمان الدقة العلمية والالتزام بالمعايير الأكاديمية المحددة. | إلى الأمام: 5’-CATGTACGTTGCTATCCAGGC-3’ عكس: 5’-CTCCTTAATGTCACGCACGAT-3’ |
| بادئات B4GALT2 | تسينغكيه | نظرًا لعدم توفر نص مصدر للترجمة، يرجى تزويدي بالنص الإنجليزي الذي ترغب في ترجمته إلى اللغة العربية وفقًا للمعايير العلمية والأكاديمية المحددة. | إلى الأمام: 5’-GGGCAGACTGCTGATCGAG-3’ عكس: 5’-CCGGTGTCTAAAGGGGATGAT-3’ |
| CB-Dock2 | لاب شير (LabShare) | عبر الإنترنت | الرسو الجزيئي |
| clusterProfiler | Bioconductor | v4.14.6 | تحليل الإثراء |
| سيتوسكيب (Cytoscape) | اتحاد سايتوسكيب (Cytoscape Consortium) | v3.8.3 | التصور الشبكي |
| DrugBank | جامعة ألبرتا | عبر الإنترنت | التنبؤ بالأدوية |
| بادئات FAM168B | تسينجكه | نظرًا لعدم وجود نص مصدر، يرجى تزويدي بالنص الذي ترغب في ترجمته وسأقوم بتحويله إلى اللغة العربية وفقًا للمعايير العلمية والأكاديمية المحددة. | إلى الأمام: 5’-TCTGGGGTTCCCTATGCAAAT-3’ عكس: 5’-GTAGGATTCGCTCCAGGATACA-3’ |
| glmnet | شبكة CRAN | v4.1 | انحدار لاسو (LASSO regression) |
| تحليل تباين مجموعة الجينات (GSVA) | Bioconductor | v1.52.3 | تحليل ssGSEA |
| مزيج فائق لتخليق الشريط الأول من cDNA بواسطة Hifair III | يايسن | 11141ES | تخليق الحمض النووي الريبوزي المتمم (cDNA) |
| مزيج Hieff RTPCR SYBR Green الرئيسي | ياسين | 11201ES | تضخيم تفاعل البوليميراز المتسلسل الكمي (qPCR) |
| limma | Bioconductor | v3.54.0 | التعبير التفاضلي |
| نظام LightCycler 480 II | Roche | LightCycler 480 II | تفاعل البوليميراز المتسلسل في الوقت الحقيقي (Real-time PCR) |
| بادئات PTRF | تسينغكه | بما أنه لم يتم توفير نص مصدري للترجمة، يرجى تزويدي بالنص المراد ترجمته من الإنجليزية إلى العربية، وسأقوم بتنفيذ ذلك وفقاً للمعايير المهنية والعلمية المحددة في التعليمات. | تقديم: 5’-GGGCCGTAGACCAGATCCA-3’ عكس: 5’-CTTGCTCACCGTATTGCTCGT-3’ |
| بادئات PRKCDBP | تسينغكي | يرجى تزويدي بالنص المصدر الذي ترغب في ترجمته. | الأمام: 5’-CACGTTCTGCTCTTCAAGGAG-3’ عكس: 5’-TGTACCTTCTGCAATCCGGTG-3’ |
| برنامج R | مؤسسة R | v4.4.2 | الحوسبة الإحصائية |
| الغابة العشوائية | CRAN | v4.7 | الغابة العشوائية |
| طقم استخلاص الحمض النووي الريبوزي من الدم MolPure Blood RNA Kit | ييزين | 19241ES50 | استخلاص الحمض الريبي النووي (RNA) |
| قاعدة بيانات STRING | EMBL | عبر الإنترنت | شبكة التفاعلات بين البروتينات |
| بادئات TOR3A | تسينغكي | نظرًا لعدم توفر نص مصدر، يرجى تزويدي بالمحتوى الذي ترغب في ترجمته من الإنجليزية إلى العربية وفقًا للمعايير والتعليمات المذكورة. | تقديم: 5’-CCCTTGCTCTGTCGTTCCAC-3’ عكس: 5’-CCCGTCCCGATACAGGTTC-3’ |
| بادئات UPP1 | تسينغكي | نظرًا لعدم توفر نص مصدر للترجمة، يرجى تزويدي بالنص الذي ترغب في ترجمته من الإنجليزية إلى العربية وفقًا للمعايير العلمية والأكاديمية المحددة. | أمامي: 5’-CTGTCAGTCATGGTATGGGCA-3’ عكس: 5’-GAGCACCGGGCATAGTACA-3’ |
| تحليل شبكة الارتباط الجيني المرجح (WGCNA) | شبكة أرشيف R الشاملة (CRAN) | v1.72 | شبكة التعبير المشترك |
| XGBoost | شبكة أرشيف R الشاملة (CRAN) | الإصدار 1.7 | خوارزمية XGBoost |
الوصول مقيد. يرجى تسجيل الدخول أو بدء فترة تجريبية لعرض هذا المحتوى.