DeepOmicsAE: ייצוג מודולי איתות במחלת אלצהיימר עם ניתוח למידה עמוקה של פרוטאומיקה, מטבולומיקה ונתונים קליניים

Summary

DeepOmicsAE היא זרימת עבודה המתמקדת ביישום שיטת למידה עמוקה (כלומר, מקודד אוטומטי) כדי להפחית את המימדיות של נתוני multi-omics, מתן בסיס למודלים חיזוי ומודולי איתות המייצגים שכבות מרובות של נתוני omics.

Abstract

מערכי נתונים גדולים של omics הופכים זמינים יותר ויותר למחקר על בריאות האדם. מאמר זה מציג את DeepOmicsAE, זרימת עבודה המותאמת לניתוח מערכי נתונים רב-אומיים, כולל פרוטאומיקה, מטבולומיקה ונתונים קליניים. זרימת עבודה זו משתמשת בסוג של רשת עצבית הנקראת autoencoder, כדי לחלץ קבוצה תמציתית של תכונות מנתוני הקלט הרב-מימדיים הגבוהים. יתר על כן, זרימת העבודה מספקת שיטה לייעול הפרמטרים העיקריים הדרושים ליישום המקודד האוטומטי. כדי להציג את זרימת העבודה הזו, נותחו נתונים קליניים מעוקבה של 142 אנשים שהיו בריאים או אובחנו עם מחלת אלצהיימר, יחד עם הפרוטאום והמטבוליזם של דגימות המוח שלהם לאחר המוות. התכונות המופקות מהשכבה הסמויה של המקודד האוטומטי שומרות על המידע הביולוגי המפריד בין חולים בריאים לחולים. בנוסף, התכונות הבודדות שחולצו מייצגות מודולי איתות מולקולרי נפרדים, שכל אחד מהם מקיים אינטראקציה ייחודית עם התכונות הקליניות של הפרטים, ומספק אמצעי לשילוב הפרוטאומיקה, המטבולומיקה והנתונים הקליניים.

Introduction

חלק הולך וגדל מהאוכלוסייה מזדקן והנטל של מחלות הקשורות לגיל, כגון ניוון עצבי, צפוי לעלות בחדות בעשורים הקרובים¹. מחלת אלצהיימר היא הסוג הנפוץ ביותר של מחלות נוירודגנרטיביות². ההתקדמות במציאת טיפול הייתה איטית בהתחשב בהבנה הלקויה שלנו של המנגנונים המולקולריים הבסיסיים המניעים את הופעת המחלה והתקדמותה. רוב המידע על מחלת האלצהיימר מתקבל לאחר המוות מבדיקת רקמת המוח, מה שהפך את ההבחנה בין סיבות ותוצאות למשימה קשה³. פרויקט מחקר המסדרים הדתיים / זיכרון והזדקנות (ROSMAP) הוא מאמץ שאפתני להשיג הבנה רחבה יותר של ניוון עצבי, הכולל מחקר של אלפי אנשים שהתחייבו לעבור בדיקות רפואיות ופסיכולוגיות מדי שנה ולתרום את מוחם למחקר לאחר מותם⁴. המחקר מתמקד במעבר מתפקוד תקין של המוח למחלת אלצהיימר². במסגרת הפרויקט, דגימות מוח לאחר המוות נותחו עם שפע של גישות אומיקס, כולל גנומיקה, אפיגנומיקה, טרנסקריפטומיקס, פרוטאומיקה⁵ ומטבולומיקה.

טכנולוגיות Omics המציעות קריאות פונקציונליות של מצבים תאיים (כלומר, פרוטאומיקה ומטבולומיקה)^6,7 הן המפתח לפענוח מחלה ^8,9,10,11,12, בשל הקשר הישיר בין חלבון ושפע מטבוליטים ופעילויות תאיות. חלבונים הם המבצעים העיקריים של תהליכים תאיים, בעוד מטבוליטים הם המצע והתוצרים של תגובות ביוכימיות. ניתוח נתונים מולטי-אומיקה מציע את האפשרות להבין את היחסים המורכבים בין נתונים פרוטאומיים ומטאבולומיים במקום להעריך אותם בנפרד. מולטי-אומיקס היא דיסציפלינה החוקרת שכבות מרובות של נתונים ביולוגיים בממדים גבוהים, כולל נתונים מולקולריים (רצף גנום ומוטציות, תעתוק, פרוטום, מטבוליזם), נתוני הדמיה קלינית ותכונות קליניות. בפרט, ניתוח נתונים מולטי-אומיקס נועד לשלב שכבות כאלה של נתונים ביולוגיים, להבין את הוויסות ההדדי שלהם ואת דינמיקת האינטראקציה, ולספק הבנה הוליסטית של התפרצות המחלה והתקדמותה. עם זאת, שיטות לשילוב נתונים מולטי-אומיים נותרו בשלבים המוקדמים של הפיתוח¹³.

Autoencoders, סוג של רשת עצבית¹⁴ ללא פיקוח, הם כלי רב עוצמה לשילוב נתונים multi-omics. שלא כמו רשתות עצביות מפוקחות, מקודדים אוטומטיים אינם ממפים דגימות לערכי מטרה ספציפיים (כגון בריאים או חולים), והם אינם משמשים לחיזוי תוצאות. אחד היישומים העיקריים שלהם טמון בהפחתת ממדיות. עם זאת, מקודדים אוטומטיים מציעים מספר יתרונות על פני שיטות פשוטות יותר להפחתת ממדיות, כגון ניתוח רכיבים עיקריים (PCA), הטבעה של שכן סטוכסטי מבוזר t (tSNE), או קירוב סעפת אחידה והקרנה (UMAP). שלא כמו PCA, מקודדים אוטומטיים יכולים ללכוד קשרים לא ליניאריים בתוך הנתונים. שלא כמו tSNE ו- UMAP, הם יכולים לזהות יחסים היררכיים ורב-מודאליים בתוך הנתונים מכיוון שהם מסתמכים על שכבות מרובות של יחידות חישוביות שכל אחת מהן מכילה פונקציות הפעלה לא ליניאריות. לכן, הם מייצגים מודלים אטרקטיביים כדי ללכוד את המורכבות של נתונים multi-omics. לבסוף, בעוד היישום העיקרי של PCA, tSNE ו- UMAP הוא זה של קיבוץ הנתונים, מקודדים אוטומטיים דוחסים את נתוני הקלט לתכונות מחולצות המתאימות היטב למשימות חיזוי במורד הזרם^15,16.

בקצרה, רשתות עצביות מורכבות מכמה שכבות, שכל אחת מהן מכילה יחידות חישוביות מרובות או "נוירונים". השכבה הראשונה והשכבה האחרונה נקראות שכבות הקלט והפלט, בהתאמה. מקודדים אוטומטיים הם רשתות עצביות בעלות מבנה שעון חול, המורכב משכבת קלט, ואחריה שכבה אחת עד שלוש שכבות נסתרות ושכבה "סמויה" קטנה המכילה בדרך כלל בין שניים לשישה נוירונים. חציו הראשון של מבנה זה מכונה מקודד והוא משולב עם מפענח המשקף את המקודד. המפענח מסתיים בשכבת פלט המכילה את אותו מספר תאי עצב כמו שכבת הקלט. מקודדים אוטומטיים לוקחים את הקלט דרך צוואר הבקבוק ומשחזרים אותו בשכבת הפלט, במטרה לייצר פלט המשקף את המידע המקורי קרוב ככל האפשר. זה מושג על ידי מזעור מתמטי של פרמטר המכונה "אובדן שחזור". הקלט מורכב מסט של תכונות, אשר ביישום המוצג כאן יהיו שפע חלבונים ומטבוליטים, ומאפיינים קליניים (כלומר, מין, השכלה וגיל בעת המוות). השכבה הסמויה מכילה ייצוג דחוס ועשיר במידע של הקלט, שניתן להשתמש בו ליישומים הבאים כגון מודלים לחיזוי^17,18.

פרוטוקול זה מציג זרימת עבודה, DeepOmicsAE, הכוללת: 1) עיבוד מקדים של פרוטאומיקה, מטבולומיקה ונתונים קליניים (כלומר, נורמליזציה, קנה מידה, הסרה חריגה) כדי להשיג נתונים בקנה מידה עקבי לניתוח למידת מכונה; 2) בחירת תכונות קלט מקודד אוטומטי מתאימות, שכן עומס יתר על תכונות עלול לטשטש דפוסי מחלה רלוונטיים; 3) אופטימיזציה ואימון של המקודד האוטומטי, כולל קביעת המספר האופטימלי של חלבונים ומטבוליטים לבחירה, ושל נוירונים לשכבה הסמויה; 4) חילוץ תכונות מהשכבה הסמויה; ו-5) ניצול התכונות המחולצות לפרשנות ביולוגית על ידי זיהוי מודולי איתות מולקולרי והקשר שלהם עם תכונות קליניות.

פרוטוקול זה נועד להיות פשוט וישים על ידי ביולוגים עם ניסיון חישובי מוגבל שיש להם הבנה בסיסית של תכנות עם Python. הפרוטוקול מתמקד בניתוח נתונים מולטי-אומיקס, כולל פרוטאומיקה, מטבולומיקה ותכונות קליניות, אך ניתן להרחיב את השימוש בו לסוגים אחרים של נתוני ביטוי מולקולרי, כולל תעתוק. יישום חדשני חשוב אחד שהוצג על ידי פרוטוקול זה הוא מיפוי ציוני החשיבות של תכונות מקוריות על נוירונים בודדים בשכבה הסמויה. כתוצאה מכך, כל נוירון בשכבה הסמויה מייצג מודול איתות, המפרט את האינטראקציות בין שינויים מולקולריים ספציפיים לבין המאפיינים הקליניים של החולים. הפרשנות הביולוגית של מודולי האיתות המולקולרי מתקבלת באמצעות MetaboAnalyst, כלי זמין לציבור המשלב נתוני גנים/חלבונים ומטבוליטים כדי להפיק מסלולי איתות מטבוליים ותאים מועשרים¹⁷.

Protocol

הערה: הנתונים ששימשו כאן היו נתוני ROSMAP שהורדו מפורטל הידע של AD. אין צורך בהסכמה מדעת כדי להוריד את הנתונים ולעשות בהם שימוש חוזר. הפרוטוקול המוצג כאן משתמש בלמידה עמוקה כדי לנתח נתונים מולטי-אומיים ולזהות מודולי איתות המבחינים בין מטופלים ספציפיים או קבוצות מדגם המבוססות, למשל, על האבחנה שלהם. הפרוטוקול גם מספק קבוצה קטנה של תכונות שחולצו שמסכמות את הנתונים המקוריים בקנה מידה גדול, וניתן להשתמש בהן לניתוח נוסף, כגון אימון מודל חיזוי באמצעות אלגוריתמים של למידת מכונה (איור 1). עיין בקובץ משלים 1 ובטבלת החומרים לקבלת מידע אודות גישה לקוד והגדרת הסביבה החישובית לפני ביצוע הפרוטוקול. יש לבצע את השיטות לפי הסדר המפורט להלן.

איור 1: סכמה של זרימת העבודה של DeepOmicsAE. ייצוג סכמטי של זרימת העבודה לניתוח נתוני multi-omics באמצעות זרימת העבודה. בתיאור המקודד האוטומטי, מלבנים מייצגים שכבות של הרשת העצבית ומעגלים מייצגים נוירונים בתוך שכבות. אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.

1. עיבוד מקדים של נתונים

הערה: מטרת סעיף זה היא לעבד מראש את הנתונים, כולל טיפול בנתונים חסרים; נרמול וקנה מידה של נתונים פרוטאומיים, מטבוליים וקליניים; והסרת חריגים. הפרוטוקול מיועד למערך נתונים הכולל נתוני פרוטאומיקה המבוטאים כיומן₂(יחס); נתוני מטבולומיקה המתבטאים בשינוי קיפול; ומאפיינים קליניים הכוללים מאפיינים רציפים וקטגוריים. יש לקבץ את החולים או הדגימות על סמך אבחנה או פרמטרים דומים אחרים. דגימות או מטופלים צריכים להיות לאורך השורות והתכונות לאורך העמודות.

1. כדי להתחיל מופע חדש של Jupyter Notebook בדפדפן, פתח חלון מסוף חדש, הקלד את הטקסט הבא והקש Enter.
   מחברת יופיטר
2. בדף הבית של Jupyter בדפדפן, לחץ על המחברת M01 - expression data pre-processing.ipynb כדי לפתוח אותה בכרטיסייה חדשה (קובץ משלים 2, שלב 1.1).
3. בתא השני של המחברת, הקלד את שם קובץ ערכת הנתונים במקום your_dataset_name.csv.
4. בתא האחרון של המחברת, הקלד את השם הרצוי של קובץ נתוני הפלט במקום M01_output_data.csv.
5. בתא החמישי של המחברת, ציין את מיקום העמודות עבור כל סוג נתונים באופן הבא: נתוני פרוטאומיקה (cols_prot), נתוני מטבולומיקה (cols_met), נתונים קליניים רציפים (לדוגמה, גיל) (cols_clin_con), נתונים קליניים בינאריים (לדוגמה, מין) (cols_clin_bin). הזן את אינדקס העמודות הראשון עבור כל סוג נתונים במקום col_start ואת אינדקס העמודות האחרונות במקום col_end; לדוגמה: cols_prot = פרוסה(0, 8817). ודאו שהערכים שצוינו בעצמי הפרוסה תואמים לאינדקסים של העמודה הראשונה והאחרונה המתאימים לכל סוג נתונים. השתמש בפקודה בתא הרביעי של אותה מחברת (df.iloc[:, :]) כדי לקבוע את מיקום ההתחלה והסיום עבור כל סוג נתונים (קובץ משלים 2, שלב 1.2).
6. בחר תא | הפעל הכל משורת התפריטים ב- Jupyter כדי ליצור את קובץ נתוני הפלט בתיקייה שצוינה (קובץ משלים 2, שלב 1.3).
   הערה: נתונים אלה ישמשו כקלט עבור הפרוטוקולים המתוארים בסעיפים 2, 3 או 4.

2. אופטימיזציה מותאמת אישית של זרימת העבודה (אופציונלי)

הערה: סעיף 2 הוא אופציונלי מכיוון שהוא דורש צורך רב במחשב. משתמשים צריכים לדלג ישירות לסעיף 4 אם הם מחליטים לא לבצע את סעיף 2. פרוטוקול זה ינחה את המשתמש דרך אופטימיזציה של זרימת העבודה באופן אוטומטי. באופן ספציפי, השיטה מזהה את הפרמטרים המספקים את הביצועים הטובים ביותר של המקודד האוטומטי במונחים של יצירת תכונות שחולצו שמפרידות היטב בין קבוצות המדגם. הפרמטרים הממוטבים הנוצרים כפלט כוללים את מספר התכונות לשימוש לבחירת תכונות (k_prot ו- k_met) ואת מספר תאי העצב בשכבה הסמויה של המקודד האוטומטי (סמוי). לאחר מכן ניתן להשתמש בפרמטרים אלה בפרוטוקול המתואר בסעיף 3 כדי ליצור את המודל.

1. בדף הבית של Jupyter בדפדפן, לחץ על המחברת M02 - DeepOmicsAE model optimization.ipynb כדי לפתוח אותה בכרטיסייה חדשה (קובץ משלים 2, שלב 2.1).
2. בתא השני של המחברת, הקלד את שם קובץ הקלט במקום M01_output_data.csv. הקלט לפונקציה זו הוא נתוני הפלט מסעיף 1.
3. בתא החמישי של המחברת, ציין את מיקום העמודות עבור כל סוג נתונים באופן הבא: נתוני פרוטאומיקה (cols_X_prot), נתוני מטבולומיקה (cols_X_met), נתונים קליניים (cols_clin; כולל כל הנתונים הקליניים), כל נתוני הביטוי המולקולרי, כולל נתוני פרוטאומיקה ומטבולומיקה (cols_X_expr). הזן את אינדקס העמודות הראשון עבור כל סוג נתונים במקום col_start ואת אינדקס העמודות האחרונות במקום col_end; לדוגמה, cols_prot = פרוסה(0, 8817). ודא שהערכים שצוינו באובייקטי הפרוסה תואמים לאינדקס העמודות הראשונה והאחרונה המתאים לכל סוג נתונים, והשתמש בפקודות בתא השלישי והרביעי של המחברת כדי לחקור את הנתונים ולקבוע את מיקום ההתחלה והסיום של כל סוג נתונים. ציין את שם העמודה הכוללת את משתנה היעד במקום y_column_name כ- y_label (קובץ משלים 2, שלב 2.2).
   הערה: הערכים של האינדקסים שצוינו ב- cols_X_prot, cols_X_met, cols_clin ו- cols_X_expr יהיו שונים מאלה המשמשים בסעיף 1 עקב העיצוב מחדש של מסגרת הנתונים המתרחש במהלך עיבוד מקדים של נתונים.
4. בתא השישי של המחברת, ציין כמה סבבי מיטוב יש לבצע על-ידי הקצאת ערך ל - n_comb. זמני העיבוד הם כ 4-5 דקות במשך 10 סיבובים; 20 דקות ל-50 סיבובים, ו-40 דקות ל-100 סיבובים (קובץ משלים 2, שלב 2.3).
5. בחר תא | הפעל הכל משורת התפריטים ביופיטר.
   הערה: משתני הפלט kprot, kmet ו - latent יאוחסנו וניתן יהיה לגשת אליהם מהמחברות האחרות, אשר ישמשו להמשך זרימת העבודה האנליטית. AE_optimization_plot.pdf העלילה ייווצר ויישמר בתיקייה המקומית (איור 2).

3. יישום זרימת עבודה עם פרמטרים מותאמים אישית

הערה: בצע פרוטוקול זה רק לאחר מיטוב השיטה (סעיף 2). אם משתמשים בוחרים שלא לבצע מיטוב שיטה, דלג ישירות לסעיף 4. פרוטוקול זה ינחה את המשתמש ביצירת מודל באמצעות הפרמטרים המותאמים אישית הנגזרים מסעיף 2. המקודד האוטומטי 1) ייצור קבוצה של תכונות שחולצו המשחזרות את הנתונים המקוריים ו-2) יזהה את התכונות החשובות המניעות כל נוירון בשכבה הסמויה, תוך ייצוג יעיל של מודולי איתות ייחודיים. מודולי האיתות יפורשו באמצעות הפרוטוקול המפורט בסעיף 5.

1. בדף הבית של Jupyter בדפדפן, לחץ על המחברת M03a - יישום DeepOmicsAE עם פרמטרים מותאמים אישית.ipynb כדי לפתוח אותו בכרטיסייה חדשה (קובץ משלים 2, שלב 3.1).
2. בתא השני של המחברת, הקלד את שם קובץ הקלט במקום M01_output_data.csv. הקלט לפונקציה זו הוא נתוני הפלט מסעיף 1.
3. בתא החמישי של המחברת, ציין את מיקום העמודות עבור כל סוג נתונים באופן הבא: נתוני פרוטאומיקה (cols_prot), נתוני מטבולומיקה (cols_met), נתונים קליניים (cols_clin; כולל את כל הנתונים הקליניים). הזן את אינדקס העמודות הראשון עבור כל סוג נתונים במקום col_start ואת אינדקס העמודות האחרונות במקום col_end; לדוגמה: cols_prot = פרוסה(0, 8817). ודא שהערכים שצוינו באובייקטי הפרוסה תואמים לאינדקסים של העמודה הראשונה והאחרונה המתאימים לכל סוג נתונים, והשתמש בפקודות בתא השלישי והרביעי של המחברת כדי לחקור את הנתונים ולקבוע את מיקום ההתחלה והסיום של כל סוג נתונים. ציין את שם העמודה המכילה את משתנה היעד (לדוגמה, 0 או 1, המתאים לבריא או חולה) במקום y_column_name כ- y_label.
   הערה: ערך האינדקסים המצוין ב- cols_X_prot, cols_X_met, cols_clin ו- cols_X_expr יהיה שונה מאלה המשמשים בסעיף 1 עקב העיצוב מחדש של מסגרת הנתונים המתרחש במהלך עיבוד מקדים של נתונים.
4. בחר תא | הריצו הכל משורת התפריטים בצדק כדי ליצור ולשמור את החלקות PCA_initial_data.pdf, PCA_extracted_features.pdf ו-distribution_important_feature_scores.pdf בתיקייה המקומית (איור 3 ואיור משלים S1). בנוסף, רשימות של תכונות חשובות עבור כל מודול איתות מזוהה יאוחסנו בקבצי טקסט בתיקייה המקומית, בשם module_n.txt, כאשר n יוחלף במספר המודול.

4. יישום זרימת עבודה עם פרמטרים מוגדרים מראש

1. עיין בסעיף 3 לקבלת הוראות מפורטות כיצד להפעיל שיטה זו (קובץ משלים 2, שלב 4.1). ההבדל היחיד בין שני הפרוטוקולים האלה הוא שהפרמטרים kprot, kmet ו-latent (בתא השביעי של המחברת) נגזרים מתמטית בהתבסס על תוצאות האופטימיזציה שבוצעה כפי שמוצג באיור 2.
   הערה: אם סעיף 4 מספק הפרדה גרועה בין קבוצות המדגם, דבר המצביע על ביצועי מודל לא אופטימליים, מומלץ לבצע אופטימיזציה של המודל (סעיף 2) באמצעות לפחות 15 חזרות, ואם אפשר, עד 50.

5. פרשנות ביולוגית באמצעות MetaboAnalyst

1. פתח את הדפדפן ונווט לקישור למטה כדי לגשת לפונקציונליות ניתוח המסלול המשותף באתר MetaboAnalyst : https://www.metaboanalyst.ca/MetaboAnalyst/upload/JointUploadView.xhtml.
2. גש לתיקייה שבה נשמרו קבצי הפלט משיטה 3 או שיטה 4 ופתח את קבצי הטקסט module_n.txt עבור כל מודול איתות n שנוצר על-ידי שיטה 3 או על-ידי שיטה 4.
3. אתר את החלבונים בקובצי הטקסט והעתק אותם.
4. הדבק את רשימת החלבונים בחלון גנים/חלבונים עם שינויי קיפול אופציונליים בדף האינטרנט של MetaboAnalyst.
5. חזור על השלב לעיל עבור מטבוליטים והדבק אותם ברשימת התרכובות בחלון עם שינויי קיפול אופציונליים באותו דף אינטרנט.
6. בחר את האורגניזם המתאים ואת סוג המזהה ולאחר מכן לחץ על שלח בתחתית הדף (קובץ משלים 2, שלב 5.1).
   הערה: ודא שהמזהים מזוהים על-ידי MetaboAnalyst. מזהים מוכרים כוללים מזהה Entrez, סמלי גנים רשמיים ומזהה Uniprot לחלבונים; שם תרכובת, מזהה HMDB ומזהה KEGG עבור מטבוליטים. אם המזהים אינם מסוגים אלה, יש צורך בהמרה מתאימה לפני הניתוח.
7. בדף הבא, בדוק את מיפוי המזהים לפני שתלחץ על המשך כדי לוודא שהמזהים מזוהים.
8. בדף הגדרת פרמטרים , בחר מסלולים מטבוליים (משולבים) או כל המסלולים (משולבים) כדי להציג בהתאמה את תרומת הקלט למסלולים מטבוליים בלבד או לכל מסלולי האיתות (קובץ משלים 2, שלב 5.2). בחלונית בחירת אלגוריתם , בחר Enrichment analysis: Hypergeometric test, Topology measure: degree centrality, and integration method: Combine p values (pathway-level). לחץ על שלח בתחתית הדף.
9. העמוד האחרון הוא תצוגת התוצאות, המציגה את תוצאות ניתוח ההעשרה. מסלולים מועשרים מתווטים על פי השפעתם ומשמעותם, ורשימת המסלולים מסופקת גם בפורמט טבלאי.

תוצאות

כדי להציג את הפרוטוקול, ניתחנו מערך נתונים הכולל פרוטום, מטבוליזם ומידע קליני שמקורו במוחות לאחר המוות של 142 אנשים שהיו בריאים או אובחנו עם מחלת אלצהיימר.

לאחר ביצוע פרוטוקול סעיף 1 כדי לעבד מראש את הנתונים, מערך הנתונים כלל 6,497 חלבונים, 443 מטבוליטים, ושלושה מאפיינים קליניים (?...

Discussion

מבנה מערך הנתונים הוא קריטי להצלחת הפרוטוקול ויש לבדוק אותו בקפידה. יש לעצב את הנתונים כפי שמצוין בסעיף 1 של הפרוטוקול. הקצאה נכונה של מיקומי טורים היא גם קריטית להצלחת השיטה. נתוני פרוטאומיקה ומטבולומיקה מעובדים מראש באופן שונה ובחירת התכונות מתבצעת בנפרד בשל אופיים השונה של הנתונים. לכן...

Materials

Name	Company	Catalog Number	Comments
Computer	Apple	Mac Studio	Apple M1 Ultra with 20-core CPU, 48-core GPU, 32-core Neural Engine; 64 GB unified memory
Conda v23.3.1	Anaconda, Inc.	N/A	package management system and environment manager
conda environment DeepOmicsAE	N/A	DeepOmicsAE_env.yml	contains packages necessary to run the worflow
github repository DeepOmicsAE	Microsoft	https://github.com/elepan84/DeepOmicsAE/	provides scripts, Jupyter notebooks, and the conda environment file
Jupyter notebook v6.5.4	Project Jupyter	N/A	a platform for interactive data science and scientific computing
DT01-metabolomics data	N/A	ROSMAP_Metabolon_HD4_Brain 514_assay_data.csv	This data was used to generate the Results reported in the article. Specifically, DT01-DT04 were merged by matching them based on the individualID. The column final consensus diagnosis (cogdx) was filtered to keep only patients classified as healthy or AD. Climnical features were filtered to keep the following: age at death, sex and education. Finally, age reported as 90+ was set to 91, then the age column was transformed to float64. The data is available at https://adknowledgeportal.synapse.org
DT02-TMT proteomics data	N/A	C2.median_polish_corrected_log2 (abundanceRatioCenteredOn MedianOfBatchMediansPer Protein)-8817x400.csv
DT03-clinical data	N/A	ROSMAP_clinical.csv
DT04-biospecimen metadata	N/A	ROSMAP_biospecimen_metadata .csv
Python 3.11.3	Python Software Foundation	N/A	programming language