تم تنفيذ الخطوات من 2 إلى 9 من التحليل باستخدام برنامج R (انظر الشكل التكميلي 1 والملف التكميلي 1) ، بينما تم تنفيذ الخطوات المتبقية على المنصات عبر الإنترنت. يتم توفير تفاصيل قواعد البيانات المستخدمة في هذا البروتوكول (جنبا إلى جنب مع روابط الويب) في جدول المواد.
1. الحصول على البيانات
- الوصول إلى قاعدة بيانات التعبير الجيني الشامل (GEO) المتاحة للجمهور في المركز الوطني لمعلومات التكنولوجيا الحيوية.
- ابحث في بيانات GEO المسماة مرض الزهايمر في مربع البحث.
- حدد الكائنات الحية العليا مثل الإنسان العاقل على الجانب الأيمن.
ملاحظة: كانت نتائج البحث بيانات عن مرض الزهايمر في الإنسان العاقل.
- بعد تصفية المعلومات التي تم البحث عنها ، قم بتنزيل ملفات البيانات GSE167490 و GSE183068 ، والتي تشمل features.tsv و barcode.tsv و matrix.mtx لكل عينة نوى مفردة. وتألفت مجموعات البيانات من 34 عينة من مرض الزهايمر من القشرة الأمامية، مع توزيع متساو ل 17 عينة من الذكور و 17 عينة من الإناث (الجدول التكميلي 1).
2. دمج العينة
- قم بتكوين مسارات البيانات وأسماء العينات وفقا لذلك على الكمبيوتر. قم باستيراد 34 عينة تم تنزيلها وقم بتعيين أسماء خاصة بالجنس للعينات باستخدام أسماء الوظائف.
- قم بإنشاء كائنات Seurat لجميع العينات بطريقة معالجة الدفعات باستخدام قائمة الوظائف و Read10X ، مع تحديد المعلمات على أنها min.cells = 3 و min.features = 200.
- استخدم الدالة RenameCells لإضافة معرفات نموذجية كبادئات إلى الرموز الشريطية للخلية للحفاظ على الرموز الشريطية للخلية أثناء عملية الدمج. وقد ضمن ذلك احتفاظ كل خلية بهويتها الفريدة ويمكن إرجاعها إلى مصدر العينة الأصلي بعد الدمج.
3. مراقبة الجودة (QC)
- استخدم الدالة PercentFeatureSet لحساب نسب جينات الميتوكوندريا ونسب جين كرات الدم الحمراء ونسب جين الريبوسوم لكل خلية.
- قم بتخزين هذه النسب المحسوبة في بيانات التعريف باستخدام عامل التشغيل [[ ]] لإرفاق هذه المعلومات مباشرة ببيانات تعريف كل خلية.
- استخدم دالة المجموعة الفرعية لإجراء ترشيح الخلية ، وتحديد المعلمات على النحو nFeature_RNA > 200 ، nFeature_RNA < 10000 ، nCount_RNA < 60000 ، percent.mt < 10 ، percent.rb < 5 ، والنسبة المئوية. HB < 75.
- استبعاد GSM5106107 من التحليل.
4. فحص تأثير الدفعة
- إجراء معالجة البيانات.
- تطبيع البيانات باستخدام الدالة NormalizeData .
- حدد أهم 2000 معلم متغير في مجموعة البيانات باستخدام الدالة FindVariableFeatures .
- إجراء تحليل المكونات الرئيسية (PCA)17 على البيانات باستخدام RunPCA، مع الاحتفاظ ب 50 مكونا رئيسيا.
- قم بإنشاء مخطط كوع باستخدام وظيفة ElbowPlot لتحديد العدد الأمثل للأبعاد للتحليل اللاحق. النظر في أول 50 بعدا.
- قم بتوسيع نطاق البيانات باستخدام ScaleData للتأكد من أن جميع الميزات على نطاق قابل للمقارنة.
- حدد أقرب الجيران باستخدام FindNeighbors استنادا إلى 30 بعدا.
- قم بتطبيق خوارزمية UMAP باستخدام RunUMAP لتقليل أبعاد البيانات إلى 30 بعدا.
- تصور البيانات المعالجة باستخدام الدالة DimPlot مع تعيين معلمة الاختزال على umap وتعيين معلمة group.by على orig.ident.
ملاحظة: قد تؤدي هذه الخطوة إلى إنشاء مخطط تصور البيانات في مساحة UMAP المخفضة، مجمعة حسب هويات الخلايا الأصلية. عند فحص مؤامرات UMAP ، أصبح من الواضح أن هناك وجود تأثير دفعة. يشير التجميع المميز أو فصل الخلايا بناء على دفعتها أو أصلها التجريبي إلى أن الدفعات التجريبية قد أثرت على ملامح التعبير الجيني.
5. تكامل البيانات
- تطبيع البيانات وتوحيدها باستخدام وظيفة SCTransform .
- قم بتطبيق خوارزمية الانسجام18 لدمج البيانات ال 33 المتبقية أحادية النواة. استخدم مقايسة SCT للتكامل واضبط الحد الأقصى لعدد تكرارات الانسجام على 20.
- استخدم الدالة FindClusters مع تعيين معلمة دقة إلى 0.07 لتحديد مجموعات مميزة داخل البيانات.
- استخدم الدالة RunUMAP مع عدد محدد من الأبعاد (dims = 30) لتقليل أبعاد البيانات بشكل أكبر وتصور المجموعات في مساحة ذات أبعاد أقل.
6. التعليق التوضيحي لنوع الخلية
- جمع جينات الواسمات (الجدول التكميلي 2) للخلايا من خلال مراجعة شاملة للأدبيات الموجودة.
- بعد تحديد عدم تجانس الكتلة الخلوية ، صنف نوع كل خلية عنقودية حسب جينات العلامة المعبر عنها على وجه التحديد.
- تقديم أنواع خلوية مختلفة مع تصور UMAP باستخدام حزمة ggplot2 ، حيث تم تسليط الضوء على oligodendrocyte برمز اللون #DB7093 ، والخلايا العصبية المثيرة مع #FF69B4 ، والخلايا النجمية مع #1874CD ، والخلايا الدبقية الصغيرة مع #63B8FF ، وخلية سلائف oligodendrocyte مع #DB7093 ، والخلايا العصبية المثبطة مع #FFC0CB ، والخلية البطانية مع #FF69B4.
- احسب نسب كل نوع خلية طبقيا حسب الجنس.
7. استخراج بيانات الخلايا الدبقية
- استخراج بيانات الخلايا النجمية من البيانات المجمعة المتكاملة باستخدام وظيفة المجموعة الفرعية.
- استخراج بيانات الخلايا الدبقية الصغيرة من البيانات المجمعة المتكاملة باستخدام وظيفة المجموعة الفرعية.
- استخراج بيانات oligodendrocyte من البيانات المجمعة المتكاملة باستخدام وظيفة مجموعة فرعية.
8. التقاط الجينات المعبر عنها تفاضليا الخاصة بالجنس الدبقي (DEGs)
- حدد DEGs الخاصة بالجنس للخلايا النجمية باستخدام وظيفة FindMarkers (ident.1 = ذكر ، ident.2 = أنثى ، group.by = group.sum ، مقايسة = RNA) مع قيم العتبة: قيمة p < 0.05 و |avg_log2FC| > 30. قم بتسمية DEGs المنظمة لأعلى على أنها أعلى ، و DEGs المنظمة لأسفل على أنها لأسفل ، والباقي على أنها مستقرة.
- تصور DEGs باستخدام وظيفة ggplot ، حيث يمثل المحور x الفرق في النسبة المئوية بين شرطين (pct.1 - pct.2) ، والمحور y الذي يصور avg_log2FC. تم تسليط الضوء على الجينات المنظمة باستخدام اللون PaleVioletRed ، والجينات المنظمة مع Pink ، والجينات المستقرة مع DodgerBlue3.
- حدد DEGs الخاصة بالجنس للخلايا الدبقية الصغيرة باستخدام دالة FindMarkers (ident.1 = ذكر ، ident.2 = أنثى ، group.by = group.sum ، مقايسة = RNA) مع قيم العتبة: قيمة p < 0.05 و |avg_log2FC| > 1. قم بتسمية DEGs المنظمة لأعلى على أنها أعلى ، و DEGs المنظمة لأسفل على أنها لأسفل ، والباقي على أنها مستقرة.
- تصور DEGs باستخدام وظيفة ggplot ، حيث يمثل المحور x الفرق في النسبة المئوية بين شرطين (pct.1 - pct.2) ، والمحور y الذي يصور avg_log2FC. تم تسليط الضوء على الجينات المنظمة باستخدام اللون OrangeRed ، والجينات المنظمة مع LightSalmon ، والجينات المستقرة مع SteelBlue1.
- حدد DEGs الخاصة بالجنس من oligodendrocyte باستخدام دالة FindMarkers (ident.1 = ذكر ، ident.2 = أنثى ، group.by = group.sum ، مقايسة = RNA) مع قيم العتبة: قيمة p < 0.05 و |avg_log2FC| > 10. قم بتسمية DEGs المنظمة لأعلى على أنها أعلى ، و DEGs المنظمة لأسفل على أنها لأسفل ، والباقي على أنها مستقرة.
- تصور DEGs باستخدام دالة ggplot ، حيث يمثل المحور x الفرق في النسبة المئوية بين شرطين (pct.1 - pct.2) ، والمحور y الذي يصور avg_log2FC. تم تسليط الضوء على الجينات المنظمة باستخدام اللون DeepPink ، والجينات المنظمة مع HotPink ، والجينات المستقرة مع DeepSkyBlue3.
9. تحليلات الإثراء الوظيفي ل DEGs الخاصة بالجنس
- إجراء تحليل إثراء أنطولوجيا الجينات (GO) على DEGs الخاصة بالجنس لكل نوع من الخلايا الدبقية باستخدام وظيفة enrichGO . قم بتعيين المعلمات التالية: OrgDb = org. Hs.eg.db ، keyType = SYMBOL ، ont = ALLL ، pAdjustMethod = BH ، pvalueCutoff = 0.01 ، و qvalueCutoff = 0.05.
- تحويل رموز الجينات إلى معرفات الجينات المقابلة باستخدام دالة bitr. إجراء تحليل إثراء موسوعة كيوتو للجينات والجينوم (KEGG) على DEGs الخاصة بالجنس لكل نوع من أنواع الخلايا الدبقية باستخدام وظيفة enrichKEGG . اضبط الإعدادات على النحو التالي: الكائن الحي = لديه ، نوع المفتاح = kegg ، pAdjustMethod = BH ، pvalueCutoff = 0.01 ، و qvalueCutoff = 0.05.
10. إحصائيات التردد ل DEGs الدبقية في مسارات go و kegg ، ومخططات venn لكل DEGs الخاصة بالجنس الدبقي ، وبناء شبكة PPI
- احسب تواتر DEGs الخاصة بالجنس الدبقي في مسارات GO و KEGG باستخدام الرسم البياني للتردد.
- الوصول إلى قاعدة بيانات STRING لإنشاء شبكات PPI.
- اختر البروتينات المتعددة. ابحث عن قائمة الأسماء في مربع البحث. تعيين "الكائنات الحية" كما الإنسان العاقل.
- راجع قائمة البروتينات التي تم الحصول عليها من البحث. انقر فوق "متابعة " للمتابعة.
- قم بتصدير شبكات PPI عن طريق تحديد خيار التنزيل ، ويفضل أن يكون ذلك بتنسيق PNG بدقة أعلى.
- تصور توزيع التعبير المشترك للجينات الخاصة بالجنس باستخدام مخططات فن.
- حدد الجين (الجينات) المشتركة كجين (جينات) رئيسي في الدراسة بناء على تحليل مخطط فين.
11. بناء شبكة تنظيمية متعددة العوامل
- قم بالوصول إلى محلل الشبكة.
- انقر فوق إدخال قائمة الجينات وحدد الكائن الحي على أنه H. العاقل (الإنسان). قم بتعيين نوع المعرف كرمز جيني رسمي. أدخل اسم الجين في حقل البحث ، ثم انقر فوق تحميل ومتابعة.
- حدد تفاعلات Gene-miRNA واختر miRTarBase v8.0. قم بتأكيد التحديد بالنقر فوق موافق.
- انتقل إلى تفاعلات TF-gene وحدد قاعدة بيانات ENCODE . انقر فوق "موافق " لتأكيد التحديد.
- بعد ذلك ، انتقل إلى شبكة التنظيم المشترك TF-miRNA وانقر فوق "موافق " للمتابعة.
- أخيرا ، اختر متابعة لإنشاء شبكة تنظيمية متعددة العوامل تتضمن تفاعلات الجينات miRNA وتفاعلات TF-gene.
12. تحليل الجينات والطب الصيني التقليدي المستهدف
- الوصول إلى قاعدة بيانات Coremine Medical عبر الإنترنت.
- أدخل اسم الجين المحدد وحدد الجين المقابل مع الجين / البروتين اللاحق ، الإنسان في مربع البحث ضمن قسم الاستكشاف .
- انتقل إلى قسم الأدوية وحدد الطب الصيني التقليدي المرتبط بالأدوية التي تم البحث عنها.
ملاحظة: تم تمييز الأدوية ذات الدلالة الإحصائية باللون الأزرق.
- حدد أعلى خمسة مقاييس للطب الصيني التقليدي بناء على قيمتها "الأهمية" كمقاييس علاجية للطب التقليدي.
13. ملخص بحثي لمكونات الطب الصيني التقليدي في استهداف الجينات الرئيسية
- الوصول إلى منصة الأبحاث التكاملية القائمة على علم الأدوية للطب الصيني التقليدي (TCMIP) وقاعدة بيانات علم الأدوية وأنظمة الطب الصيني التقليدي ومنصة التحليل (TCMSP). أدخل أسماء الأعشاب في شريط البحث لاسترداد المكونات المقابلة لها.
- استرجع المكونات الموجودة في قاعدة بيانات PubMed مع حد زمني حتى 10 أبريل 2023. تضمنت مصطلحات البحث المستخدمة اسم الجزيء في TCMSP والمكونات الكيميائية في TCMIP كمصطلحات بحث ، واقتصرت على المقالات المنشورة باللغة الإنجليزية.
- تلخيص وتحليل الأعشاب والمكونات المقابلة لها التي تعمل على AD.
14. تأكيد استهداف وظيفة علاج الطب الصيني التقليدي في خصوصية جنس مرض الزهايمر
- استيراد الأعشاب إلى TCMIP وانتقل إلى صفحة الوصف المقابلة.
- استخدم وظيفة تصدير البيانات وحدد تنسيق CSV لتنزيل شروط التخصيب الخاصة ب GO - العملية البيولوجية ، GO - المكون الخلوي ، GO - الوظيفة الجزيئية ، ومسار المفاعل.
- تصور شروط الإثراء التي تم تنزيلها لكل عشب باستخدام المخططات الشريطية.