1. اختيار مجموعات البيانات والدليل
- لتحديد مجموعات البيانات للتحليل، انقر فوق الزر بحث عن مجموعات البيانات ( الشكل2A، المربع الأحمر 1) للانتقال إلى المجلد الذي يتم إيداع البيانات فيه (على سبيل المثال، مجلد RawData). ملفات البيانات تلقائيا بملء ملفات البيانات كقائمة. يمكن أن يكون هناك أكثر من مجموعة بيانات واحدة في المجلد.
- انقر على زر اختيار الدليل وحدد الدليل (على سبيل المثال، اختبار) حيث سيتم إيداع الملفات المحللة (الشكل 2A، المربع الأحمر 2) . سيتم إنشاء مجموعة من المجلدات وملفات التحليل النهائية بالإضافة إلى الصور المؤقتة المتوسطة في هذا الدليل عند الضغط على زر التحليل. سيتم إنتاج أخطاء إذا لم يتم اختيار دليل.
2. تعيين حجم بكسل ومعدل الإطارات
- ملء معدل الإطار المناسب وحجم بكسل من الصور في المناسبة "Framerate" أو "حجم بكسل" مربع (الشكل 2A، مربع أخضر). إذا لم يتم توفير أية قيم (يتم تعيينها إلى الافتراضي "0")، ثم سيقوم البرنامج بالبحث في بيانات التعريف المقترنة بالملف لحجم معدل الإطارات والبكسل. إذا تعذر العثور على هذه القيم، فسيتخلف البرنامج عن كل بكسل للقياس وفي كل إطار للنقاط الزمنية.
ملاحظة: في المثال المتوفر، معدل الإطارات هو 100، وحجم بكسل هو 0.08.
3. اختيار أو جعل أقنعة
- استخدم الزر Maker قناع لإنشاء أقنعة الخلية تلقائيا للبيانات في قائمة مجموعة بيانات الملف (الشكل 2A، المربع الأزرق). عند استخدام الزر Maker قناع سيتم إنشاء مجلد جديد في الدليل المختار يسمى MaskFiles. يتحول "مؤشر التشغيل" إلى اللون الأصفر أثناء التشغيل ويعود إلى اللون الأخضر عند اكتماله.
ملاحظة: سيتم إنشاء قناع لكل ملف في قائمة ملفات البيانات من الإطارات 10 الأولى من ملف الصورة (أو من كافة الإطارات إذا كان أقل من 10 في ملف الصورة) وإيداعها في المجلد MaskFiles باستخدام نظام التسمية المناسبة (الموضحة أدناه).
- تأكد من أن ملفات القناع تلقائيا بملء قائمة "ملفات قناع". يمكن للمستخدم المتابعة مباشرة إلى التحليل.
ملاحظة: تحقق من ملفات القناع بشكل مرئي دائما وتأكد من التقاطها لمنطقة الاهتمام بأكملها. يتم عرض الإطار الأول من ملف البيانات وملف القناع على واجهة المستخدم عند تحديدها. (الشكل2ب). قد ينتج "صانع القناع" أخطاء في حالة انخفاض الإشارة إلى الضوضاء، لذا فإن التحقق من صحة أن ملفات القناع مناسبة أمر بالغ الأهمية لمراقبة الجودة.
- كبديل لاستخدام Mask Maker، إذا كانت الإشارة إلى ضجيج الصور غير كافية، قم بإنشاء أقنعة يدويا في ImageJ.
- أولا، افتح ملف الصورة الخام في ImageJ (الشكل 3A).
- انقر فوق الزر تحديد المضلع وانقر فوق لرسم قناع حول الخلية. بمجرد الانتهاء، انقر نقرا مزدوجا فوق النقطة الأخيرة لإكمال المضلع.
- بمجرد الانتهاء، انتقل إلى تحرير | | التحديد إنشاء قناع (الشكل 3B). سيتم إنشاء قناع مقلوب جديد استنادا إلى الرسم المضلع. حفظ أقنعة في مجلد MaskFiles المعينة في الدليل المختار. يجب أن يطابق نظام تسمية ملف القناع ملفات البيانات الفردية المطابقة متبوعا ب "_mask_file". على سبيل المثال، إذا كان ملف بيانات يسمى "VAMP2_488_WT_1.tif"، يجب أن يكون اسم ملف القناع المطابق "VAMP2_488_WT_1_mask_file.tif".
- استخدم الزر بحث عن Maskfiles للانتقال إلى المجلد المختار من ملفات القناع المخصص المودعة. سوف الأقنعة ملء "ملفات قناع" كقائمة.
ملاحظة: من المهم أن يكون قناع لكل ملف بيانات قبل تشغيل التحليل.
4. تحليل واستخراج ميزة
- بعد اختيار الدليل وملء قائمة ملفات البيانات وقائمة ملفات القناع، انقر فوق الزر تحليل (الشكل 4A). إذا كان تصنيف الأحداث exocytic مطلوبا، انتقل إلى الخطوة 5.
ملاحظة: سيؤدي النقر فوق الزر تحليل سلسلة من المهام التلقائية لتحليل البيانات. فإنه سيتم إنشاء مجلدات فردية في الدليل المختار لإيداع البيانات التي تم تحليلها. أثناء التشغيل، سيتغير "مؤشر التشغيل" من الأخضر إلى الأصفر (الشكل 4A، المربع الأحمر). بعد الانتهاء من التحليل، سيتغير هذا مرة أخرى إلى اللون الأخضر.
- العثور على مجلد DataFiles (الشكل 4B) مع مجموعة كاملة من ملفات التحليل (وكذلك ملفات استخراج ميزة، لاستخدامها في التصنيف لاحقا) اسمه وفقا لكل ملف البيانات (الشكل 4C).
ملاحظة: يتم تضمين وصف ملفات التحليل هذه أدناه في المقطع نتائج الممثل.
5. تصنيف الأحداث الإكستية
- لإجراء تصنيف متزامن للأحداث الإكسوسيتيكية مع الكشف التلقائي، تحقق من خانة الاختيار تصنيف قبل النقر على زر التحليل. لكل حدث exocytic تعيين درجة احتمال بين 0-1 لكل فئة. يعتبر الحدث exocytic مصنفا كواحد من الفئات الأربع إذا كانت درجة الاحتمال لتلك الفئة > 0.5.
ملاحظة: بمجرد اكتمال التحليل، سيظهر مجلد ملف بيانات جديد ضمن الدليل المختار. سيحتوي المجلد على ملفات تحليل مطابقة لكل ملف صورة.
6. تحليل سباتيوتمبورالي للنصرة باستخدام قيم ريبلي K
- إنشاء ملف قناع "neurite" و "soma" منفصل. أولا، قطع سوما من النيوريت. لا توجد طريقة غير متحيزة لتقسيم سوما من neurites ، لذلك يجب أن يكون المستخدم أعمى عن تجربة التكييف واستخدام أفضل حكم . ويقترح بيضاوي مع عدم وجود ملحقات neurite واضحة.
- فتح ImageJ / فيجي.
- اسحب ملف القناع وأسقطه أو استخدم ملف | افتح ثم حدد ملف القناع.
- استخدام أداة منتقي الألوان وانقر على أي من بكسل أسود في خلفية القناع لتعيين اللون إلى الأسود.
- استخدم أداة تحديد المضلع أو حر لرسم مخطط تفصيلي حول سوما، وفصله عن neurites. وهذا يتطلب اتخاذ قرارات يدوية.
- انقر على تحرير | | التحديد جعل قناع. سيتم فتح صورة جديدة مع سوما دائرة مجزأة من بقية الصورة. سيؤدي هذا إلى إنشاء ملف قناع سوما.
- دون تحريك المنطقة المرسومة، انقر فوق رأس ملف القناع الأصلي.
- انقر على تحرير | ملء لملء سوما دائرة بحيث تبقى فقط neurites القناع.
- بمجرد الحصول على ملفات neurite وقناع منفصلة حفظ كل من ملفات القناع.
- فتح "neurite_2D_network" في ماتلاب.
- في MATLAB، انتقل إلى الدليل مع كافة بيانات التحليل.
- تغيير المسار "اسم القناع" إلى اسم قناع neurite، أي"MaskFiles/VAMP2pHluorin_488_wt_4_mask_file_neurite.tif"
- تغيير "csv_file_name" إلى حيث يوجد ملف fluorescent_traces.csv، أي"MaskFiles/VAMP2pHluorin_488_wt_4_mask_file_neurite.csv".
- انقر فوق تشغيل لهياكل عظمية ملف قناع neurite. يؤدي هذا إلى إنشاء إصدار skeletonized من ملف قناع neurite و إيداعه كملف CSV ضمن المجلد maskfiles.
- بعد ذلك، إنشاء ملف CSV ل soma. فتح "CSV_mask_creator.m" في ماتلاب.
- وضع في المسار ل "اسم قناع" إلى اسم قناع سوما، أي، "MaskFiles / VAMP2phLuorin_488_wt_4_mask_file_soma.tif
- تغيير "writematrix" إلى اسم ملف csv ليتم إنشاؤها، أي" MaskFiles / VAMP2pHluorin_488_wt_4_mask_file_soma.csv "
- انقر فوق تشغيل. يؤدي هذا إلى إنشاء ملف VAMP2pHluorin_488_wt_4_mask_file_soma.csv جديد.
- كرر 6.14 لكل ملف قناع.
7. إعداد RStudio
- فتح Rstudio وفتح ripleys_k_analysis. R ملف.
- تثبيت الحزمة "spatstat" في RStudio عن طريق الانتقال إلى أدوات | تثبيت الحزم والكتابة في "spatstat" متبوعا بالنقر فوق تثبيت.
ملاحظة: هذا يحتاج فقط إلى تنفيذ مرة واحدة لكل تثبيت Rstudio.
- تشغيل مكتبة spatstat في بداية كل جلسة عمل.
- انتبه إلى متغيرين رئيسيين في هذه الحالة: "neuron_mask" و "neuron_datapoints". قناع الخلايا العصبية يشير إلى مجموعة من ملفات قناع لتشغيل، أي، سوما قناع الملفات.
ملاحظة: تشغيل كافة ملفات قناع سوما معا، بشكل منفصل عن ملفات قناع Neurite، والعكس بالعكس.
- قراءة في .csv من ملفات قناع لكل من الخلايا العصبية ليتم تحليلها.
- تشغيل ملفات متعددة في وقت واحد عن طريق نسخ neuron_mask_n+1 إضافية. سيسمح هذا بتجميع تحليل ريبلي معا ، باستخدام السيناريو الموصوف في القسم 8.
- الآن تحقق من المتغير الثاني ، "neuron_datapoints".
- اقرأ في. X_fluorescent_traces.csv" الملف الذي تم إنشاؤه بواسطة برنامج التحليل (بواسطة ميزات كافة extracted_R ملف) مع المواضع x,y,t الأحداث exocytic، وكذلك ملف neurite الخاصة بمواضع x,y لشبكة الاتصال 2D. هذا يذهب في موقف "neuron_datapoints".
- في RStudio حدد رمز | تشغيل | المنطقة تشغيل كافة. وهذا يولد العديد من المؤامرات، بما في ذلك قيم ريبلي K المجمعة وكذلك قطع الكثافة.
- حفظ قطع الأراضي عن طريق الذهاب إلى تصدير | حفظ الصورة كما | حدد تنسيق الصورة المناسب والدليل، ثم أدخل اسم الملف المناسب، ثم انقر فوق حفظ.
ملاحظة: يتم استدعاء دالة رسم خرائط الحرارة في البرنامج النصي باستخدام "plot(density(soma_data,0.4)"." يمثل الرقم "0.4" هنا مدى سلاسة وظيفة الكثافة. يمكن تغييره ليتناسب مع بيانات المستخدم بطريقة ذات معنى، ولكن إذا كان سيتم إجراء المقارنات بين خرائط الحرارة المختلفة، يجب أن يكون الرقم هو نفسه بينهما.
-
تصدير أو حفظ الصورة من Rstudio. إذا تطلبت خريطة الحرارة المزيد من التحرير، فاختر نوع ملف مناسب (SVG أو EPS).
8. تحليل ريبلي
ملاحظة: Ripleys_k_analysis. R الملف أيضا تلقائيا ولدت ريبلي ك قيمة المؤامرات. تشغيل البرنامج النصي بأكمله سيتم تلقائيا تشغيل الوظائف المذكورة أدناه، ولكن يتم تضمينه بالتفصيل إذا كان أحد يرغب في تشغيل كل جزء من البرنامج النصي بشكل فردي أو إجراء تغييرات على التحليل.
- أولا، تشغيل دالة المغلف لكل خلية. هذه الدالة محاكاة عشوائية المكانية كاملة (CSR) لاختبار القيمة K ريبلي من نمط نقطة الحدث exocytic ضد.
Data_envelope_1 = المغلف (soma_data_1، Kest، nsim = 19، savefuns = TRUE)
Data_envelope_2 = المغلف (soma_data_2، Kest، nsim = 19، savefuns = TRUE)
- بعد ذلك، تجمع هذه المغلفات معا وإنشاء تقدير واحد من المسؤولية الاجتماعية للشركات للمجموعة:
Pool_csr = تجمع (Data_envelope_1، Data_envelop_2,...)
بعد ذلك، قم بتشغيل الدالة K ريبلي لكافة نقاط البيانات.
Data_ripleys_k_1 = Kest(soma_data_1، نسبة = TRUE)
Data_ripleys_k_2 = Kest(soma_data_2، نسبة = TRUE)
- بمجرد الانتهاء، تجمع قيم ريبلي K معا وbootstrap فترات الثقة الخاصة بهم مع الأمر التالي:
تجمع البيانات = تجمع (Data_ripleys_k_1، Data_ripleys_k_2,...)
- Bootstrap مع الأمر التالي:
Final_Ripleys_K = varblock (متعة = Kest ، Data_pool)
- رسم البيانات.
- إذا كان هناك حاجة إلى اختلاف إحصائي ثابت، يتم تضمين اختبار التباديل Studentised في حزمة spatstat لاختبار الفرق بين مجموعات من أنماط النقاط:
Test_difference = studpermu.test (all_points_to_test ، exocytic_events ~ المجموعة ، nperm = NP).