ملاحظة: يتم عرض جميع البرامج المستخدمة هنا في جدول المواد. كانت البيانات المستخدمة في هذه الدراسة لأغراض العرض التوضيحي من مشروع الشبكة العصبية البشرية (HCP: http://www. humanconnectome.org)15. تمت الموافقة على جميع الإجراءات التجريبية من قبل مجلس المراجعة المؤسسية (IRB) في جامعة واشنطن. تم الحصول على بيانات التصوير في مجموعة بيانات HCP باستخدام ماسح ضوئي 3T Siemens Skyra معدل مع ملف رأس 32 قناة. تم تفصيل معلمات الحصول على الصور الأخرى في ورقة سابقة16. تم تنزيل الحد الأدنى من البيانات المعالجة مسبقا للعرض التوضيحي ، والذي انتهى من خطوات المعالجة المسبقة التالية: تصحيح تشوه التدرج ، وتصحيح الحركة ، والمعالجة المسبقة للخريطة الميدانية ، وتصحيح التشوه المكاني ، والتطبيع المكاني لمساحة معهد مونتريال للأعصاب (MNI) ، وتطبيع الشدة ، وإزالة مجال التحيز. يمكن أيضا استخدام بيانات التصوير بالرنين المغناطيسي الوظيفي في حالة الراحة من مشاريع الباحثين.
1. المعالجة المسبقة للبيانات
- تحقق من جودة البيانات واستبعد المشاركين الذين لديهم بيانات إعادة اختبار مفقودة وحركة رأس مفرطة (3 مم في الترجمة و 3 درجات في الدوران).
ملاحظة: تمت إزالة خمسة مشاركين ، وتم تضمين 38 شابا (22-35 عاما) في التحليلات الرئيسية.
- افتح مربع أدوات تحليل الشبكة النظرية للرسم البياني (GRETNA)17 في MATLAB لتنفيذ المزيد من خطوات المعالجة المسبقة. انقر فوق دفعة FC Matrix Construction. حدد مسار مجموعة البيانات الوظيفية لتحميل مستندات NIFTI وتنفيذ الخطوات التالية ، كما هو موضح في خيار خط الأنابيب في الشكل 1:
- قم بإزالة أول 10 صور بالنقر المزدوج فوق رقم النقطة الزمنية لإزالتها في إزالة الصور الأولى وإدخال 10.
- سلس مكانيا (العرض الكامل بنصف الحد الأقصى [FWHM] = [4 4 4] بالنقر المزدوج فوق FWHM (مم) في السلس المكاني وإدخال [4 4 4]).
- انحدار المتغيرات المشتركة. اختر إشارات المادة البيضاء وإشارات CSF وحركة الرأس على أنها TRUE. حدد القناع المناسب وفقا لحجم voxel الفعلي ، على سبيل المثال ، قناع ب 2 مم هنا ، واختر معلمات Friston-24 لحركة الرأس.
- تصفية مؤقتا. أدخل قيمة TR وفقا لوقت تكرار فحص التصوير بالرنين المغناطيسي (على سبيل المثال ، 720 مللي ثانية هنا) وقم بإزالة الضوضاء عالية التردد ومنخفضة التردد بالنقر المزدوج فوق النطاق (هرتز) وإدخال [0.01 0.1].
ملاحظة: النتائج مع وبدون انحدار إشارات الدماغ بأكملها معروضة أدناه. عند استخدام البيانات غير المعالجة مسبقا ، يوصى أيضا باستخدام خطوط أنابيب راسخة مثل fMRI-prep18 ومساعد معالجة البيانات للتصوير بالرنين المغناطيسي الوظيفي في حالة الراحة (DPARSF) 19 .

الشكل 1: المعالجة المسبقة للتصوير بالرنين المغناطيسي الوظيفي Rs-fMRI وتقدير اتصال الشبكة الوظيفي. إعدادات المعالجة المسبقة (إزالة أول 10 صور ، والتنعيم المكاني باستخدام FWHM من 4 مم ، وإلغاء الاتجاه الخطي زمنيا ، وتراجع إشارات المادة البيضاء ، وإشارات السائل النخاعي (CSF) ، وحركة الرأس مع 24 معلمة ، وتصفية النطاق من 0.01-0.1 هرتز) والارتباط الثابت مع تحويل فيشر Z. الاختصارات: Rs-fMRI = التصوير بالرنين المغناطيسي الوظيفي في حالة الراحة ؛ FWHM = العرض الكامل عند نصف الحد الأقصى ؛ CSF = السائل النخاعي. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
2. بناء الشبكة وتحليلها
ملاحظة: يتم تلخيص سير العمل العام لبناء وتحليلات شبكة الملاحة في ثلاث خطوات رئيسية (الشكل 2).

الشكل 2: سير العمل العام لبناء وتحليل شبكة الملاحة. (أ) اختر التنقل كمصطلح للبحث عنه في قاعدة بيانات Neurosynth. (ب) يمكن إنشاء قائمة بإحداثيات التنشيط. (ج) قم بإجراء تحليل تلوي باستخدام وظائف من Neurosynth للحصول على العديد من خرائط الدماغ. (د، ه) من خلال دمج الخريطة التحليلية التلوية وأطلس تقسيم الدماغ بالكامل (AICHA) ، يمكن إنشاء العقد (ROI). (و) بناء شبكة ملاحة باستخدام عقد الملاحة الناتجة وتوصيلها الوظيفي (تقدير التوصيلية وتحليل الشبكة). الاختصارات: ROI = منطقة الاهتمام ؛ AICHA = أطلس الاتصال الجوهري للمناطق المتجانسة. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
- تعريف عقدة الشبكة
- قم بتنزيل أحدث قاعدة بيانات Neurosynth (neurosynth.org)20 عن طريق كتابة الأمر في Python:
>استيراد neurosynth ك ns
>ns.dataset.download (المسار = './' ، فك الضغط = صحيح)
ملاحظة: يحتوي أرشيف مجموعة البيانات ('current_data.tar.gz') على ملفين: 'database.txt' و'features.txt'. تحتوي هذه على جميع إحداثيات التنشيط من مقالات التصوير العصبي وعلامات التحليل التلوي التي تحدث بتردد عال في تلك المقالة ، على التوالي.
- قم بإنشاء مثيل Dataset جديد من database.txt وإضافة ميزات إلى هذه البيانات عن طريق كتابة الأمر:
> من استيراد مجموعة بيانات neurosynth.base.dataset
> dataset = Dataset ('data/database.txt')
> dataset.add_features ("بيانات/features.txt")
- قم بإجراء تحليل تلوي باستخدام مصطلح الاهتمام (أي "التنقل") عن طريق كتابة الأمر:
> معرفات = dataset.get_ids_by_features ("التنقل" ، الحد = 0.01)
> ma = ميتا. التحليل التلوي (مجموعة البيانات ، المعرفات)
> ma.save_results('.'، 'الملاحة')
ملاحظة: ينتج عن التحليل التلوي العديد من خرائط الدماغ بتنسيق NIFTI. تم تطبيق عتبة معدل الاكتشاف الخاطئ (FDR) البالغة 0.01 للتحكم في المعدل الإيجابي الخاطئ. هناك حاجة إلى المعرفة المودعة في هذه الخطوة لضمان تضمين المناطق المبلغ عنها بشكل شائع في خريطة التحليل التلوي. يمكن تطبيق خطوات مماثلة لإجراء التحليلات التلوية للوظائف المعرفية الأخرى مثل اللغة والذاكرة.
- حدد مجموعات الاهتمام من خلال دمج الخريطة التحليلية التلوية وأطلس تقسيم الدماغ بالكامل عن طريق كتابة الأمر من FSL:
>fslmaths navigation_0.01.nii.gz -bin navi_bin.nii.gz
>fslmaths navi_bin.nii.gz -mul AICHA/AAL.nii.gz navi_label_aicha/aal.nii.gz
>FSLMATHS navi_label_aicha/aal.nii.gz -thr n -uthr n التسمية _n.nii.gz
>تسمية المجموعة -i _n.nii.gz -t 0.2 -o cluster_n.nii.gz
>fslmaths cluster_n.nii.gz -thr m -uthr m cluster_n_m.nii.gz
>fslmaths cluster_n_m.nii.gz -بن مول س node_x.nii.gz
>fslmaths node_1.nii.gz -إضافة ... -إضافة node_x.nii.gz navi_AICHA/AAL_mask.nii.gz
ملاحظة: تم استخدام أطلسين هنا: AAL و AICHA. AAL هو الأطلس الذي تم استخدامه في الدراسة الأصلية لتعريف العقدة12. تم إنشاء هذا الأطلس بناء على الملامح التشريحية. أطلس الاتصال الجوهري للمناطق المتجانسة (AICHA)21 له دقة أدق بكثير وتجانس وظيفي أعلى. حددنا مناطق الاهتمام باستخدام كل من الأطلس.
- اكتب البرامج النصية في Python للتحقق من حجم كل منطقة في الخريطة:
>بالنسبة لي في np.arange(n)+1:
>____region_list.الملحق (ط)
>____size1_list.append(np.sum(img_dat==i))
>____size2_list.append(np.sum(aicha_img_dat==i))
>____pct_list.append(np.sum(img_dat==i)/np.sum(aicha_img_dat==i))
ملاحظة: يشير العدد الصحيح n في البرنامج النصي إلى العدد الإجمالي للمناطق داخل تقسيم AICHA و AAL (384 و 128 على التوالي). لتجنب آثار العناقيد الزائفة ، يقترح إزالة المجموعات ذات الأحجام الصغيرة نسبيا (على سبيل المثال ، 100 فوكسل). يتم إنشاء أطلس AICHA المستخدم هنا باستخدام بيانات الاتصال الوظيفي ، حيث تظهر كل منطقة تجانس النشاط الزمني الوظيفي داخل نفسها.
- تقدير اتصال الشبكة
ملاحظة: يتم استخدام صندوق أدوات GRETNA لتقدير الاتصال وتحليل الشبكة.- انقر فوق دفعة FC Matrix Construction. قم بتحميل بيانات rs-fMRI المعالجة مسبقا عن طريق تحديد مسار مجموعة البيانات الوظيفية. انقر فوق خيار الارتباط الثابت . قم بتحميل العقدة التي تم الحصول عليها في الخطوة السابقة كأطلس لحساب الارتباط الثابت لإشارات rs-fMRI لكل زوج من المناطق ونقلها إلى درجات Fisher's z لتحسين الحالة الطبيعية.
ملاحظة: العملية التفصيلية موضحة في الشكل 1. سيتم الحصول على مصفوفات شبكة الملاحة N × N (N تمثل عدد العقد) لكل مشارك في .txt الشكل.
- احصل على شبكة موجبة وموزونة بالخطوات التالية ، كما هو موضح في الشكل 3.
- انقر فوق دفعة تحليل الشبكة. أضف مصفوفات الشبكة إلى نافذة مصفوفة اتصال الدماغ واختر دليل إخراج واحد للتحضير.
- بالنسبة لخيار خط الأنابيب لتكوين الشبكة ، حدد موجبا في علامة المصفوفة ، والتي ستضبط الاتصالات السالبة في مصفوفة اتصال الوظيفة على 0 وتزيل الاتصالات الغامضة22. اختر نوع الشبكة على أنه مرجح للحصول على الشبكة المرجحة غير الموجهة.
ملاحظة: إلى جانب الشبكات المرجحة ، يمكن للمرء أيضا تقسيم الشبكات لإنشاء شبكات ثنائية للتحليلات اللاحقة (مع مناهج مختلفة) ، ولكن غالبا ما يعتبر أن الشبكات المرجحة تظهر موثوقية أعلى23,24.
- تحليل الشبكة
- أضف العالم الصغير ، والكفاءة العالمية ، ومعامل التجميع ، وأقصر طول مسار ، ومركزية الدرجة ، والكفاءة المحلية إلى خط أنابيب التحليل المتري لشبكة GRETNA ، كما هو موضح في الشكل 3.
ملاحظة: العالم الصغير والكفاءة العالمية هما مقياسان للشبكة العالمية. على وجه التحديد ، يمكن للشبكة ذات العالم الصغير أن تزيد من كفاءة نقل المعلومات بتكلفة أسلاك منخفضة نسبيا. تعكس الكفاءة العالمية المعلومات الموازية لكفاءة النقل في شبكة النقل. بالنسبة لمقاييس الشبكة العقدية ، تقيس مركزية الدرجة عدد الروابط المتصلة بالعقدة. أقصر طول مسار ، كما هو اسمه ، هو أساس لقياس التكامل. يشير معامل التجميع إلى درجة ترابط جيران العقد مع بعضهم البعض. الكفاءة المحلية هي كفاءة التواصل مع العقدة وجيرانها (تظهر الصيغة التفصيلية والاستخدام في هذه الأوراق) 17,25. يمكن أيضا استخدام صندوق أدوات اتصال الدماغ (BCT)25 وصناديق الأدوات الأخرى لحساب مقاييس الشبكة.
- حدد تباين الشبكة في طريقة العتبة لاستبعاد التأثيرات المربكة للاتصالات الزائفة ، وأدخل مجموعة من تسلسلات العتبة (على سبيل المثال ، 0.05 ، 0.1 ، 0.15 ، 0.2 ، 0.25 ، 0.3 ، 0.35 ، 0.4 ، 0.45 ، 0.5 يستخدم هنا) لتحديد العتبة المناسبة وفقا للنتائج الإحصائية.
ملاحظة: تعرف نسبة الحواف إلى الحد الأقصى لعدد الحواف في شبكة مع نموذج عدد العقد باسم عتبة التناثر. تضمن عتبة التباين أن الأفراد المختلفين لديهم نفس عدد الحواف. اخترنا استكشاف عتبات مختلفة للتحقق من الصحة ، والتي يمكن أن توفر بيانات مفيدة لاختيار عتبة مثالية في الدراسات المستقبلية.
- قم بتعيين رقم الشبكة العشوائي على أنه 1,000 لإنشاء شبكات عشوائية باستخدام خوارزمية أسلاك ماركوف26. انقر فوق تشغيل لتشغيل خط الأنابيب في GRETNA بعد إعداد كافة الخطوات.
ملاحظة: على غرار شبكات الدماغ الحقيقية ، تحافظ الشبكات العشوائية على نفس العدد من العقد والحواف وتوزيع الدرجات. لتحديد ما إذا كانت مبنية طوبولوجيا بشكل غير عشوائي إلى حد كبير ، ستتم مقارنتها بشبكات الدماغ. بعد تشغيل خط الأنابيب ، سيتم الحصول على مجموعة من الدرجات لمقاييس الشبكة لكل من العتبات لإجراء مزيد من التحليلات الإحصائية.
- حدد العدد الأمثل للوحدات في الشبكة في أربع خطوات.
- احسب متوسط شبكة الملاحة. انقر فوق حمام مقارنة المقاييس واختر اتصال. قم بتحميل مصفوفات الشبكة التي تم الحصول عليها أعلاه واختر العملية المتوسطة (الوظيفية ). حدد اتجاه الإخراج للحفاظ على مصفوفة الشبكة المتوسطة ؛ انظر الشكل 4 لمزيد من التفاصيل.
- قسم متوسط الشبكة التي تم الحصول عليها من الخطوة أعلاه إلى وحدات 2 و 3 و 4 و 5 باستخدام العنقود الطيفي للوظيفة في MATLAB.
- احسب نسبة العقد المقسمة إلى نفس الوحدة في REST 1 و REST 2 بعد محاذاة أقسام الوحدة باستخدام البرنامج النصي procrustes_alignment.m. استخدم النسبة كمؤشر لتكرار قسم الوحدة النمطية.
- حدد عدد الوحدات ذات أعلى قابلية للتكرار.
- التحليلات الإحصائية
ملاحظة: التحليلات التالية هي أساسا للتحقق من الصحة ولن تكون ضرورية عند تطبيق هذا البروتوكول على دراسات التباين الفردية.- افحص تشابه مقاييس الشبكة هذه بين شبكتين لهما أنواع مختلفة من الاستراتيجيات لتعريف العقدة (أي أن الشبكة الجديدة التي تم إنشاؤها في الدراسة الحالية ، والتي يطلق عليها اسم NaviNet_AICHA والسابقة من Kong et al. ، تسمى NaviNet_AAL)12. احسب ارتباط بيرسون باستخدام الدالة corrcoef في MATLAB وكرر التحليلات لكل عتبة تفاوت.
ملاحظة: بعد استخراج مقاييس الشبكة ، يمكن للمرء إجراء أي تحليلات إحصائية يهتم بها.
- تحقق من موثوقية الاختبار وإعادة الاختبار لمقاييس الشبكة هذه باستخدام الدالة ICC في MATLAB27,28 ، والتي تنفذ حساب معامل الارتباط داخل الفئة.
ملاحظة: تم الإبلاغ عن قيم p الأصلية غير المصححة في قسم النتائج التمثيلية. 0.2 < يتم تفسير ICC < 0.4 على أنه مؤشر على موثوقية عادلة للاختبار وإعادة الاختبار ويتم تفسير ICC > 0.4 على أنه موثوقية متوسطة إلى جيدة للاختباروإعادة الاختبار 29,30. تم تعيين درجات ICC السلبية إلى الصفر ، بالنظر إلى حقيقة أن وجود ICCs السلبية لا معنى له ويصعب تفسيره31.

الشكل 3: تحليل مقاييس الشبكة. يحدد هذا التحليل الشبكات الإيجابية المرجحة ب 10 عتبات. احسب مقياسين للشبكة العالمية للكلمة الصغيرة والكفاءة ، وأربعة مقاييس للشبكة العقدية لمعامل التجميع ، وأقصر طول مسار ، والكفاءة ، ومركزية الدرجة. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 4: حساب متوسط شبكات الملاحة. تساعد العملية المتوسطة (الوظيفية) في حساب متوسط الشبكات لجميع المشاركين. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.