Method Article
מאמר זה מתאר את המדידה החצי-אוטומטית של המשרעת וההשהיות של חמש הפסגות והשקתות הראשונות בצורת הגל של תגובת גזע המוח השמיעתית. שגרה נוספת אוספת ומבארת את הנתונים לגיליון אלקטרוני לצורך ניתוח נסיינים. שגרות מחשב חינמיות אלה מבוצעות באמצעות החבילה הסטטיסטית בקוד פתוח R.
דיווחים רבים ב-15 השנים האחרונות העריכו שינויים בצורת הגל של תגובת גזע המוח השמיעתי (ABR) לאחר עלבונות כגון חשיפה לרעש. השינויים הנפוצים כוללים הפחתות במשרעת השיא 1 ובאיחורים היחסיים של הפסגות המאוחרות יותר, וכן עלייה ברווח המרכזי, המתבטאת בעלייה יחסית באמפליטודות של הפסגות המאוחרות יותר בהשוואה לאמפליטודה של שיא 1. נסיינים רבים מזהים את הפסגות והשקתות באופן חזותי כדי להעריך את הגבהים וההשהיות היחסיים שלהן, תהליך מייגע כאשר צורות הגל נאספות במרווחים של 5 dB בכל טווח השמיעה עבור כל תדר ומצב. מאמר זה מתאר שגרות חופשיות שעשויות להתבצע בפלטפורמת הקוד הפתוח R עם ממשק RStudio כדי להפוך את המדידות של הפסגות והשקתות של צורות גל של תגובת גזע המוח השמיעתית (ABR) לאוטומטיות. הרוטינות מזהות את המשרעת וההשהיות של פסגות ושקתות, מציגות אותן על צורת גל שנוצרה לבדיקה, אוספות ומבארות את התוצאות לגיליון אלקטרוני לניתוח סטטיסטי, ומייצרות צורות גל ממוצעות למספרים. במקרים שבהם התהליך האוטומטי מזהה באופן שגוי את צורת הגל של ABR, קיים כלי נוסף המסייע בתיקון. המטרה היא להפחית את הזמן והמאמץ הדרושים לניתוח צורת הגל של ABR, כך שחוקרים נוספים יכללו ניתוחים אלה בעתיד.
תגובת גזע המוח השמיעתית (ABR) משמשת לעתים קרובות לקביעת ספי שמיעה אצל נבדקים בעלי חיים ותינוקות אנושיים. מכיוון שה-ABR הוא תיעוד אלקטרואנצפלוגרם (EEG) של התגובות הראשונות של מערכת העצבים לגירויים שמיעתיים, הוא נושא מידע נוסף המשקף את הירי המתואם של נוירוני גנגליון ספירליים שבלוליים ועיבוד אותות מוקדם בגזע המוח השמיעתי, כולל עיבוד דו-צדדי1. תגובות אלה עשויות להיות מושפעות מטראומת רעש. לדוגמה, חשיפה לרעש שמספיקה כדי לגרום לשינוי סף זמני בעכברים יכולה גם להפחית לצמיתות את המשרעת של שיא ABR 12. יתר על כן, טראומה כזו עשויה להפחית את ההשהיות הבין-פסיכיות ולהגדיל את המשרעת היחסית של הפסגות המאוחרות3, אולי בגלל אובדן של ויסות מעכב4. בנוסף לממצאים אלה, מוטציות גנטיות ספציפיות הוכחו כמשנות את צורת הגל של ABR בהיעדר טראומה 5,6,7. לפיכך, הניתוח השגרתי של צורות גל ABR יכול לספק תובנה על מערכת השמיעה במודלים ניסיוניים.
יש גם עניין בשימוש בצורות גל ABR ככלי אבחון לחולים. דוחות קודמים העריכו אם שיא ABR 1 מופחת בחולים אנושיים לאחר חשיפה לרעש או בחולי טינטון 8,9. יש לציין כי דווח כי התקפי מיגרנה מגבירים באופן זמני את העיכובים הבין-קוליים במשך מספר שבועות, ולאחר מכן צורת הגל של ABR חוזרת לקדמותה אצל אנשים מושפעים10. דווח כי COVID-19 גורם לשינויים ארוכי טווחב-ABR interpeak latencies 11,12, אם כי מחקר אחר דיווח על תוצאות שונות 13. ליקוי שמיעה הוא לעתים קרובות תחלואה משותפת עם דמנציה בהזדקנות, ואנשים עם ליקוי שמיעה גדול יותר נוטים לחוות דמנציה שמתקדמת מהר יותר14. חוקרים חקרו שינויים בצורת הגל של ABR במחלות נוירודגנרטיביות, כגון מחלת פרקינסון (נסקרה ב-Jafari et al.15) ומחלת אלצהיימר (שנסקרה ב-Swords et al.16), כמו גם בהזדקנות רגילה 17. ככל שיותר חוקרים וקלינאים חוקרים ליקויים חושיים כסמנים ביולוגיים למחלות נפוצות בהזדקנות, טכניקות כמו ABR עשויות להפוך לשגרה בתחום הבריאות.
בחינה של פרקי השיטות בספרות מגלה כי מעבדות כותבות לעתים קרובות סקריפטים מותאמים אישית ב- MatLab לניתוח צורות גל של ABR. לפלטפורמת ABR המיוצרת על ידי מערכות שמיעה חכמות יש פונקציה לניתוח צורות גל, אך היא דורשת מהמפעיל לבחור ידנית את הפסגות והשוקת. כאן, כתבנו שגרות ניתוח חצי אוטומטיות עבור הקוד הפתוח, הסביבה הסטטיסטית הזמינה באופן חופשי R וממשק RStudio. דו"ח זה משווה את הנתונים המתקבלים באמצעות השגרה שלנו לנתונים המתקבלים על ידי כך שהנסיין מזהה באופן ידני את הפסגות והשפלים ומראה כי הנתונים משתי השיטות מתואמים מאוד. חשוב לציין שהרוטינות משלבות פונקציה מסנוורת, שבה המטא-נתונים של הדגימות ממוקמים בקובץ נפרד שאינו משולב עד הסוף. פונקציות אלה ייעלו את ניתוח צורות הגל עבור המעבדה שלנו.
כל ההליכים שבוצעו בבעלי חיים אושרו מראש על ידי הוועדה למחקר בבעלי חיים של אוניברסיטת רוצ'סטר. הנבדקים בניסוי היו 12 עכברים זכרים ונקבות מסוג F1 בגיל חודש. עכברי F1 אלה הם תוצר של הזדווגות סכר CBA/CaJ ו-C57BL/6J. העכברים גודלו ושוכנו במתקן הוויבריום עם מחזור אור/חושך סטנדרטי של 12 שעות, מזון ומים ללא הגבלה ואספקת קינון בשפע. לא יותר מחמישה אחים חד-מיניים שוכנו יחד בכלוב אחד.
1. קבלת נתונים לניתוח
הערה: שלב זה חייב לעמוד בהנחיות המוסדיות ולהיות מאושר מראש על ידי הוועדה המוסדית לרווחת בעלי חיים. התהליך המפורט להפקת נתוני ABR מעכברים תואר במקום אחר18.
2. התקנת החבילות הנדרשות וטעינת הנתונים למחשב העובד
3. קבלת הניתוח הראשוני עם FindPeaks.R
4. אימות הניתוח הראשוני
הערה: ברמות צליל נמוכות, חלקים מצורת הגל עשויים להיות קשים להבחנה מרעש, ו- FindPeaks.R עשוי לזהות באופן שגוי את הפסגות או השפלות בהשוואה לדעת הנסיין. אם יש אי התאמה, אז את הקובץ .csv ניתן לתקן עם נתונים שהתקבלו מן הסקריפט See_trace_click.R.
5. הידור והדמיה של מערך הנתונים
בדקנו את הרוטינות על תגובות בצורת גל של ABR לסדרת קליקים, החל מ-75 dB וירידה במרווחים של 5 dB ל-5 dB. נתונים אלה התקבלו כפי שתואר קודםלכן 19. בדקנו את הכלי גם על נתוני טון פיפ והשגנו תוצאות דומות. נתוני ABR מרוב מערכות ABR ניתנים לייצוא כקבצי .txt (ASCII). טענו קבצי ABR ASCII למחשב ופתחנו אותם ב- RStudio כמתואר בפרוטוקול. לאחר הרצת שגרת FindPeaks.R בצורת אצווה, השגנו צורות גל לדוגמה עם תיוג אוטומטי (איור 1) וקובץ .csv עם התוצאות. התוצאות נבדקו כדי להסיר פסגות לא רלוונטיות. כדי לאמת את התיוג האוטומטי, השתמשנו גם ביכולת התוכנית ABR כדי לתייג באופן ידני את חמש הפסגות והשפלים הראשונים בכל צורת גל שהושגה עם סדרת הלחיצות שתוארה לעיל. הנסיין שביצע משימה זו היה בעל ניסיון של שנתיים בהקלטה וניתוח של נתוני ABR. איור 2 מציג השוואה זו, כאשר הנתונים האוטומטיים של FindPeaks.R מופיעים באדום והנתונים המתקבלים באופן ידני בשחור. כל מעקב מייצג את הנתונים מעכבר יחיד. מוצג גם הממוצע עבור שתי השיטות עם סטיית תקן אחת. התוצאות שהתקבלו על ידי FindPeaks.R מתואמות מאוד עם התוצאות המתקבלות באופן ידני (ראו איור משלים S4).
איור 1: תגובה מייצגת בצורת גל לגירוי קליק של 75 dB עבור עכבר F1 צעיר. השהיה באלפיות השנייה משורטטת על ציר ה-x, והמשרעת במיקרו-וולטים משורטטת על ציר ה-y. הפסגות זוהו באופן אוטומטי עם FindPeaks.R ומסומנות באדום, ואילו השקתות מסומנות בכחול. אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור 2: השוואה של הנתונים שהתקבלו מהפסגות שזוהו באופן ידני לנתונים שסופקו על ידי ניתוח FindPeaks.R. (A,C,E,G,I) אמפליטודות במיקרו-וולטים ו-(B,D,F,H,J) באלפיות שניה משורטטות עבור רמות קול בין 5 dB ל-75 dB (ציר x, כל הגרפים) עבור פסגות I-V בצורות הגל המתקבלות עבור גירויי קליק המוצגים ל-12 עכברים. הערכים המתקבלים באופן ידני (שחור) מושווים לאותם מערכי נתונים שנותחו באמצעות FindPeaks.R (אדום). הממוצעים מתווים כקווים כבדים, כאשר האזור המוצלל מייצג סטיית תקן אחת. לא נצפו הבדלים בין השיטות כאשר הוערכו במבחן סכום הדירוג של קרוסקל-וואליס (A, הפרש = 0.0547977 ± 0.0010028, מקסימום = 0.96, p = 0.9216; B, הפרש = −0.0001734 ± 0.0001214, מקסימום = 0.04, p = 0.8289; C, הפרש = −0.0212209 ± 0.0006806, מקסימום = 0.92, p = 0.9687; D, הפרש = −0.0011047 ± 0.0001556, מקסימום = 0.06, p = 0.771; E, הפרש = −0.0323077 ± 0.0006169, מקסימום = 0.66, p = 0.899; F, הפרש = −0.0072189 ± 0.0001460, מקסימום = 0.04, p = 0.8644; G, הפרש = 0.201754 ± 0.0007407, מקסימום = 0.64, p = 0.9312; H, הפרש = −0.0007018 ± 0.0001717, מקסימום = 0.09, p = 0.8013; I, הפרש = 0.0347561 ± 0.0007343, מקסימום = 1.05, p = 0.8856; J, הפרש = −0.0078049 ± 0.0002762, מקסימום = 0.16, p = 0.886), והתוצאות היו מתואמות מאוד (ערכי צ'י בריבוע: A, 0.009696; ב, 0.046684; C, 0.0015395; ד, 0.084742; E, 0.016102; F, 0.029153; ז, 0.0074604; H, 0.063322; I, 0.020699; י, 0.020544; הבדלים המוצגים כממוצע ± SEM; מקסימום = הבדל מקסימלי מוחלט). אנא לחץ כאן כדי להציג גרסה גדולה יותר של נתון זה.
איור משלים S1: ניתוח עם FindPeaks.R. (A) בחירת ספריית העבודה (ראה פרוטוקול שלב 3.2); (B) טעינת התוכנית (ראה פרוטוקול שלב 3.3). אנא לחץ כאן כדי להוריד קובץ זה.
איור משלים S2: פלט סקריפט ופקודות לניתוח צורות הגל. פקודות עבור (A) קבצים בודדים ו- (B) עיבוד אצווה. (C) פלט קובץ PDF המציג צורות גל עם פסגות ושקתות מסומנות. ראה פרוטוקול שלב 3.4. אנא לחץ כאן כדי להוריד קובץ זה.
איור משלים S3: אימות הניתוח. (A) טעינת נתוני צורת הגל (ראה שלב 4.1 בפרוטוקול). (ב) מיקום לחצן ' הפעל אפליקציה' . קובץ הנתונים לדוגמה מצוין גם הוא. (C) חלון מבריק עם צורת גל. במקרה זה, רמת הצליל היא 75 dB, כפי שהוזן בחלון העליון. לחיצה על שיא רצוי ועל השוקת הבאה מתעדת את הנתונים עבור המשרעת וההשהיה בטבלה (שלב פרוטוקול 4.6). נתוני שיא 3 מוצגים. אנא לחץ כאן כדי להוריד קובץ זה.
איור משלים S4: השוואה של הנתונים הבודדים שהתקבלו מהפסגות שזוהו באופן ידני לנתונים שסופקו על ידי ניתוח FindPeaks.R. (א,ג,ה,ז,א) אמפליטודות במיקרו-וולטים ובאיחורים (B,D,F,H,J) באלפיות שניה משורטטות עבור רמות קול בין 5 dB ל-75 dB (ציר x, כל הגרפים) עבור פסגות I-V בצורות גל המתקבלות עבור גירויי קליקים המוצגים ל-12 עכברים. כל חיה מסומנת בצבע ייחודי, כפי שמוצג במקרא. הנתונים המתקבלים עם FindPeaks.R מסומנים בצבעים אחידים, בעוד שהנתונים המתקבלים באופן ידני מסומנים בגרסאות פחות רוויות של אותם צבעים. בעוד ששני מערכי הנתונים מתווים באיור זה, כאשר הם זהים, רק שורה אחת נראית לעין. אנא לחץ כאן כדי להוריד קובץ זה.
הפרוטוקול המתואר בפרסום זה אמור לסייע בייעול רכישת נתונים המתארים יחסי משרעת מתח ומרווחי השהיה עבור ABRs ללחיצות ופיפס טון. על-ידי שימוש בפקודות בודדות ב- RStudio, נסיינית עשויה לחלץ, לקמפל ולהציג מידע זה במסמך יחיד לצורך ניתוח סטטיסטי. על ידי הפיכת ניתוח זה לשגרה, אנו מקווים שהתחום יגלה דרכים חדשות שבהן ניתן לשנות את ה-ABR בהתפתחות, בהזדקנות או על ידי עלבון במינים שונים. מידע כזה יכול להיות בעל ערך לזיהוי מנגנונים חשובים הדומים לסינפטופתיה מרעש2. לעכברים הצעירים ששימשו לניסוי הזה היו תגובות משתנות מאוד, ככל הנראה משום שגזע המוח השמיעתי עדיין מתבגר בגיל20. עם זאת, שתי שיטות הכימות הראו קורלציות חזקות מאוד (איור 2).
הסקריפט משתמש בקובץ שנקרא "Time.csv" כדי להגדיר מרווחי זמן בתוך הנתונים לזיהוי שיא. בקצרה, משרעת מתח מקסימלית המתרחשת במרווח זמן מוגדר מסומנת "שיא 1", מתח מינימלי המתרחש במרווח הבא מסומן "שוקת 1", וכן הלאה. בחרנו את המרווחים כדי להקיף את ההשהיות של תגובות הקליק ופיפס הטון עבור עכברי CBA/CaJ בני חודש עד 12 חודשים באמצעות תדרים המשתרעים על פני 8 קילוהרץ עד 32 קילוהרץ. השתמשנו בכלי בהצלחה גם כדי למדוד תגובות של פיפ טון בעכברים. גם למינים אחרים, כולל בני אדם, יש תגובות ABR בתוך חלונות דומים, ואנו צופים שהכלי הזה יכול לשמש גם לנתונים ממינים אחרים. אנו ממליצים להשתמש בשיטת ABR המקבילה החדשה לבני אדם21, המייצרת צורות גל מצוינות. מגבלת מרווח הזמן מגבילה את השימוש בכלי זה להערכת תגובות ABR מיידיות. עם זאת, נציין כי המשתמשים יכולים לשנות את נתוני המרווח בקובץ זה כדי להפוך את המדידות של תגובות ABR לדיבור או של פוטנציאלים הקשורים לאירוע (ERPs) המתרחשים באופן אופייני בזמנים שונים בתגובה לצליל.
כמה תכונות של הטיפול הסטטיסטי של נתונים אלה ראויים להדגשה. למיטב ידיעתנו, בתחום אין טיפול סטנדרטי להבחנה בין התקדמות משרעת. מחקרים מוקדמים השתמשוב-ANOVA 22,23. הנתונים מסדרת הקליקים כאן (איור 2) היו לא פרמטריים, מה שהוביל לשימוש במבחן סכום הדירוג של קרוסקל-וואליס. בדומה ל- ANOVA, מבחן סכום הדירוג של קרוסקל-וואליס מעריך הבדלים בערכים המתקבלים ברמה נתונה של גירוי; כלומר, הוא משווה את הקווים המתקבלים על הגרף. עם זאת, טיפולים אחרים אפשריים גם כן. מבחינה ביולוגית, התקדמות המשרעת משקפת גיוס נוסף של נוירונים בעלי סף גבוה יותר ככל שרמת הגירוי עולה. זה מצביע על כך שהשטח שמתחת לעקומה, המייצג את האינטגרלים של הקווים, יכול להיות המדד הרלוונטי יותר. משוואות הערכה כלליות (GEE) עשויות לשמש למודל נתונים בודדים לניתוח אינטגרלי, כמו ב- Patel et al.5. יש לציין כי ניתוח GEE יכול לקחת בחשבון את תכנון המדדים החוזרים ונשנים של ניסויים אלה. ככל שיותר חוקרים דנים בשיטות ניתוח הנתונים, אנו צופים את הופעתו של קונצנזוס על שיטות עבודה מומלצות.
לסיכום, מאמר זה מציג כלים חינמיים וקלים לשימוש למדידה, הידור והדמיה של צורות גל ABR. כלים אלה יכולים לשמש תלמידים מתחילים של RStudio על ידי ביצוע פרוטוקול זה, והם משלבים צעד מסנוור לשיפור הקפדנות והשכפול. אנו צופים כי ניתוח שגרתי של צורת הגל של ABR יאפשר גילוי של עלבונות, וריאנטים גנטיים וטיפולים אחרים שיכולים להשפיע על תפקוד השמיעה.
למחברים אין ניגודי עניינים.
עבודה זו נתמכה על ידי שני מענקים מה- NIDCD ל- PW: R01 DC018660 ופרס משלים מנהלי, R01 DC014261-05S1.
Name | Company | Catalog Number | Comments |
C57BL/6J mice | Jackson Labs | 664 | |
CBA/CaJ mice | Jackson Labs | 654 | |
E-series PC | Dell | n/a (this equipment was discontinued) | This runs the IHS system. |
Mini-anechoic chamber | Industrial Acoustics Company | Special order number 104306 | This enclosure reduces noise levels for auditory testing of animals. |
Optiplex 7040 | Dell | i5-6500 | Rstudio may also be run on a Mac or Linux system. |
Universal Smart Box | Intelligent Hearing Systems | n/a (this equipment was discontinued) | Both TDT and IHS can output hearing data as ASCII files. |
Request permission to reuse the text or figures of this JoVE article
Request PermissionThis article has been published
Video Coming Soon
Copyright © 2025 MyJoVE Corporation. All rights reserved