שינויים בתרגום חלבונים ובתהליכים קטבוליים של חומצה קרבוקסילית במחלת כליות סוכרתית

Jul 05, 2024

תַקצִיר: מחלת כליות סוכרתית (DKD)נשאר הגורם המוביל למחלת כליות סופניתלמרות עשרות שנות לימוד. שינויים בגלומרולוס וצינוריות כליותשניהם תורמים לפתוגנזה של DKDלמרות שרוב מאמצי החקירה התמקדו בגלומרולוס. ביקשנו לבחון את חתימת הביטוי הדיפרנציאלי של DKD אנושי בגלומרולוס ובצינורית הפרוקסימלית ולאשש את הממצאים שלנו במודל העכבר db/db של סוכרת. ניתוח רשת תעתיק של נתוני RNAseq מגלומרולוס אנושי מיקרו-דיסק בלייזר (LMD) וצינורית פרוקסימלית של DKD וכריתת נפרקטומית התייחסותדגימות חשפו מסלולים מועשרים כולל קולטנים דמויי רודופסין, איתות חוש ריח וריבוזום (תרגום חלבון) בצינורית הפרוקסימלית של דגימות ביופסיה של DKD אנושית. מסלול התרגום הועשר גם בגלומרולוס. תרגום מוגבר בכליה סוכרתיתs אומת באמצעות פרופיל פוליריבוזמלי במודל העכבר db/db של סוכרת. באמצעות רצף RNA גרעיני יחיד (RNA-Seq) של כליות מעכברי db/db, נתנו עדיפות למסלולים נוספים שזוהו ב-DKD אנושי. המסלול החופף העליון שזוהה בצברי האבובים הפרוקסימליים של snRNAseq העכברים ובאבובית ה-LMD הפרוקסימלית האנושיתתא היה קטבוליזם של חומצה קרבוקסילית. באמצעות כרומטוגרפיה נוזלית-ספקטרומטריית מסה אולטרה-ביצועים, נמצא כי מסלול הקטבוליזם של חומצות השומן אינו מווסת גם במודל העכבר db/db. המטבוליט Acetyl-CoA היה מווסת מטה בעכברי db/db, תוך התאמה לביטוי הדיפרנציאלי האנושי של הגנים ACOX1 ו-ACACB. לסיכום, הממצאים שלנו מראים ששינויים צינוריים פרוקסימליים בתרגום חלבון וקטבוליזם של חומצה קרבוקסיליתהם תכונות מפתח הן בבני אדם והןDKD murine. 



מילות מפתח:חתימת ביטוי גנים; רצף RNA גרעיני יחיד;מטבולומיקה; פרופיל ריבוזומלי

1

עשבי תיבול חדשים עבורכליה סוכרתית


1. הקדמה

סוכרת (DM) משפיעה על למעלה מ-400 מיליון אנשים ברחבי העולם [1] ולמעלה מ-30 מיליון בארה"ב [2]. מחלת כליות סוכרתית (DKD) היא אחד הסיבוכים החמורים ביותר של DM, המשפיעה על 19 מיליון אנשים בארצות הברית [3], והיא הגורם העיקרי למחלת כליות סופנית(ESKD) [4]. DKD נתפסת באופן קלאסי כמחלה גלומרולרית; עם זאת, התמקדות אך ורק בפתולוגיה גלומרולרית כאתר הפציעה העיקרי מתעלמת משינויים ביולוגיים מרכזיים בתא הבין-בין-צינורי התורמים לפתופיזיולוגיה של DKD [5]. ואכן, הפרוגנוזה וההתקדמות של DKD לעיתים קרובות במקביל לאטרופיה צינורית והיקף המחלה הבין-סטילית [6,7]. רלוונטית נוספת לנקודה זו היא עדות לכך שהכוונה לאבובית הפרוקסימלית עם מעכבי SGLT2 מאט את התקדמות DKD [8].

המנגנונים הסלולריים המובילים להתפתחות ולהתקדמות של DKD הם אמנם מורכבים. הסביבה הסוכרתית מתאפיינת בהפרעות מטבוליות, דלקת, מתח חמצוני ושינויים בהומאוסטזיס חלבוני. כל ההפרעות הללו תורמים ללחץ בתאים ולהתפתחות DKD. ישנן עדויות הולכות וגדלות לכך שתרגום שונה, השלב הבסיסי בביטוי גנים שיוצר חלבונים פונקציונליים, חשוב הן בתגובה והן בתרומה ללחץ בתאים [9,10]. שינויים בתרגום תורמים לחוסר איזון בהומאוסטזיס של חלבון ומשפיעים על בריאות האדם ומחלות [11-14].

חקירה של דגימות ביופסיה של כליות אנושיות משתמשת באופן שגרתי במיקרוסקופ אור, מיקרוסקופ אלקטרוני שידור ופעולות חיסוניות כדי לאפשר לפתולוג לתת אבחנה ופרוגנוזה. עם זאת, לכלי אבחון אלו יש מגבלות בחשיפת המנגנונים המולקולריים העומדים בבסיס מחלת כליות כולל DKD. כתוצאה מכך, מודלים של בעלי חיים של DM שימשו באופן היסטורי לחקר המנגנונים הפתופיזיולוגיים של DKD; עם זאת, אף אחד מהמודלים הללו לא מדמה כראוי DKD אנושי [15]. למרות שהמאמצים האחרונים ליישם טכניקות חקירה מולקולריות מתקדמות יותר על דגימות ביופסיה של כליה אנושית מעודדים [16], אנו מאמינים ששימוש במודלים של בעלי חיים במקביל לרקמת כליה אנושית נותרה גישה ניסיונית חשובה לספק תובנות ייחודיות לגבי היסודות המולקולריים של DKD.

ביקשנו לבחון את חתימת הביטוי הדיפרנציאלי של DKD אנושי ב- glomerulus ו-proximal tubule. שיערנו שהחתימה התעתיק הצינורי הפרוקסימלי, בנוסף לזו של הגלומרולוס, קובעת גם את תוצאות ה-DKD. התייחסנו להשערה זו על ידי חקירת דגימות ביופסיה של כליה אנושית עם מיקרודיסקציה בלייזר של אזורים צינוריים גלומרולריים ופרוקסימליים כדי לזהות מסלולים הקשורים ל-DKD ולהתקדמותו. לאחר מכן, מסלולי מפתח אושרו במודל עכבר סוכרתי באמצעות גישה רב-מודאלית של רצף RNA גרעיני בודד (RNA-Seq), פרופיל פוליריבוזום ומטבולומיקה. בנוסף להפרעות בתרגום חלבונים, הממצאים שלנו תומכים בתפקיד הפוטנציאלי של חילוף החומרים של חומצה קרבוקסילית כגורם מפתח בהתקדמות מחלת כליות סוכרתית.

HERB CISTANCHE FOR DIABETIC KIDNEY

2. חומרים ושיטות

2.1. סובייקטים אנושיים

מחקר זה אושר על ידי מועצת הביקורת המוסדית של בית הספר לרפואה של אוניברסיטת אינדיאנה (IRB מס' 190657223). עשרים ושבע דגימות כליה התקבלו מקונסורציום פרויקט רפואת הדיוק של הכליות ומקוהורט הביופסיה של אינדיאנה [17]. ביופסיות כליה היו ביופסיות אינדיקציה לפרוטאוריה או לירידה בתפקוד הכליות בנבדקים עם סוכרת. 18 דגימות נרכשו ממבוגרים עם מחלת כליות סוכרתית וללא נגע גלומרולרי שני. תשע דגימות ייחוס ללא עדות היסטולוגית ל-DKD נרכשו מפרנכימה לא נאופלסטית של דגימות כריתת כליה או תורמים שנפטרו. משתנים קליניים והיסטופתולוגיים חולצו מתיק הבריאות האלקטרוני. קצב הסינון הגלומרולרי המוערך (eGFR) הוגדר כערך הקרוב ביותר למועד הביופסיה, אלא אם כן ירידה של 20% או יותר בהפניה לביופסיה של eGFR. בתרחיש זה, eGFR בתוך שנה אחת לפני הביופסיה ולפני הירידה של 20% נקבע כקו הבסיס לחישובי התקדמות.


2.2. מיקרודיסקציה בלייזר של רקמת כליה אנושית וריצוף RNA

חתכים קריוסקים בעובי 12 מיקרומטר נחתכו מגושי רקמה קפואים שנשמרו במדיום אופטימלי של טמפרטורת חיתוך, מוצמדים לשקופיות ממברנת פוליפנילן גופרתי (PPS), ועובדו באמצעות פרוטוקול Rapid Stain כפי שתואר קודם [18,19]. שטח מינימלי של 500,000 µm2 נותח עבור תאים צינוריים גלומרולריים ופרוקסימליים באמצעות לייזר UV דופק על Leica LMD6500 ואובייקטיבי 20×. רקמה מנותחת נאספה בשפופרת סטרילית נטולת RNAse המכילה RNA Extraction Buffer ו-RNA בודד לפי הוראות היצרן (Arcturus PicoPure RNA Isolation Kit, ThermoFisher, Waltham, MA, ארה"ב). איכות ה-RNA נקבעה באמצעות Agilent 2100 Bioanalyzer.

ריצוף בוצע במרכז האוניברסיטאי של אינדיאנה לליבה לגנומיקה רפואית. RNA Ribosomal היה מדולדל באמצעות פרוטוקול RiboGone-Mammalian Kit (Cat #634847, Takara Bio USA, Mountain View, CA, ארה"ב). פרוטוקול SMARTer Universal Low Input RNA Kit v2 (Cat #634938, Takara Bio USA, Mountain View, CA, ארה"ב) שימש עבור סינתזת cDNA ובניית ספרייה. ריצוף בוצע עם תצורת קצה זוגי של 2 × 75 bp ב-Illumina HiSeq 4000 באמצעות ערכת HiSeq 3000/4000 PE SBS, והנתונים הרצפים מופו לגנום hg38 באמצעות STAR. קריאות רצף ממופות באופן ייחודי הוקצו לגנים של גנום הייחוס hg38 באמצעות Rsubread featureCounts [20,21].

HERB CISTANCHE FOR DIABETIC KIDNEY



2.3. אישורי מחקר בבעלי חיים

כל הפרוטוקולים של בעלי החיים אושרו על ידי ועדת הטיפול בבעלי חיים המוסדית של אוניברסיטת אינדיאנה ותואמים למדריך הבריאות הלאומי לטיפול ושימוש בחיות מעבדה. עבור רצף RNA בתפזורת, פרופיל פוליריבוזום ומטבולומיקה, זן עכברים זכרים C57BL/6J (#000664, בקרת רקע/עכברי C57) ו-B6.BKS–Leprdb/J (#000697, עכברים db) הושגו ממעבדת ג'קסון ( בר הארבור, ME, ארה"ב). עכברים היו בני 9-12 שבועות ושקלו ~30 גרם (עכברי C57) או 40-50 גרם (עכברים db). עכברים הוקרבו וכליות נקצרו (n=4 או 5 לקבוצה).


2.4. בידוד של רקמת כליות עכבר ורצף RNA בתפזורת

הכליות הוקפאו במהירות ו-RNA חולץ באמצעות ערכת QIAGEN RNeasy Plus Midi עם עמודה להסרת DNA גנומי. איכות ה-RNA נקבעה באמצעות Agilent 2100 Bioanalyzer. ערכת הכנה לספריית mRNA של Illumina TruSeq Stranded שימשה לבניית ספרייה. ריצוף בוצע עם תצורת קצה מזווג של 2 × 75 bp על Illumina HiSeq 4000 באמצעות ערכת HiSeq 3000/4000 PE SBS, והנתונים הרצפים מופו לגנום mm10 באמצעות STAR. קריאות רצף ממופות באופן ייחודי הוקצו לגנים של הגנום הייחוס mm10 באמצעות ספירת תכונות Rsubread.



2.5. בידוד של רקמת כליה של עכבר ורצף RNA חד-גרעיני

עכברים בני 5, 8, 11, 16 ו-20 שבועות הוערכו (n=5 סה"כ, BKS.Cg-Dock7m+/+Leprdb/J #000642). הטרוזיגוטים מאותו זן שימשו כבקרות רקע (n=3). אלבומין שתן נמדד על ידי ערכת אלבומין ELISA של עכברים (Bethyl Laboratories, E99-134, Montgomery, TX, USA) וקריאטינין בשתן נמדד על ידי ערכת בדיקת קריאטינין של QuantiChrom (Bioassay System, DICT-500, Hayward קליפורניה, ארה"ב). על הקרח, חלק מכל כליה נוספה ל-1 מ"ל של חיץ תמוגה של Nuclei EZ בתוספת 2% מעכב פרוטאז (Thermofisher), ו-1% Superase (Thermofisher). לאחר הטחינה, נוספו 1 מ"ל של חיץ ורקמות הומוגגו באמצעות מטחנת הרקמות KONTES Dounce (Kimble Chase, Rockwood, TN, ארה"ב). הומוגנט היה מעורב בעדינות עם 2 מ"ל של חיץ והודגר על קרח למשך 5 דקות. ההומוגנט סונן דרך מסננת בגודל 40 מיקרומטר וטופל במאגר תמוגה בתוספת 1% Superbase, ואז סונן שוב דרך30 µמסננת מ', ולאחר מכן הונחה מחדש ב-PBS בתוספת של 2% BSA ו-1%Superase. תרחיף זה סונן לאחר מכן דרך 5µממסנן והריכוזהותאם ל-1-3 מיליון גרעינים/מ"ל לפני ההגשה לרצף.
רצף RNA גרעיני יחיד 30 בוצע במרכז אינדיאנה לרפואת גנומיקה ליבה באמצעות מערכת Chromium תא יחיד גרסה 3 (10× Genomics, סן פרנסיסקו, קליפורניה, ארה"ב) NovaSeq6000 סיקוונסר (אילומינה, סן דייגו, קליפורניה, ארה"ב). Cell Ranger 4.0 שימש ליצירת קבצי FASTQ ספציפיים לדגימה וקריאות היו מיושרות לגנום ההתייחסות mm10 באמצעות STAR. Seurat v 3.0.1 שימש לשילוב דגימות באמצעות מדדי בקרת האיכות הבאים: ספירת גנים כללה בין 200-3000 ואחוז גן מיטוכונדריאלי נמוך מ-50%. בסך הכל, 58,405 תאים נשמרו לניתוח במורד הזרם. עיבוד מקדים סטנדרטי, בחירת תכונה, הפחתת ממדים (20 רכיבים עיקריים), זיהוי עוגנים בין דגימות ואינטגרציה היומבוצע על פי שיטות עיגון Seurat v3.אשכולות צוינו בהתבסס עלסמני ביטוי גנים נפוצים [22] ושימוש בסייר תאי כליה [23]. 

HERB CISTANCHE FOR DIABETIC KIDNEY


2.6. בידוד של רקמות כליות עכבר ופרופיל פוליריבוזום

בקרת רקע ועכברי db הוקרבו בגיל 9 ו-12 שבועות. זלוף לב דרך החדר השמאלי בוצע עם 6 מ"ל של cycloheximide (100 מיקרוגרם/מ"ל ב-PBS). הכליות שנקטפו הונחו מיד במאגר תמוגה המורכב מ-1% Triton X-100, 0.1% deoxycholate, 20 מ"מ Tris-HCl, 100 מ"מ NaCl, 10 מ"מ MgCl2, טבלית קוקטייל מעכבי פרוטאז ללא EDTA (Roche, Penzberg) , גרמניה), 50 מיקרוגרם/מ"ל ציקלוהקסימיד ו-RNAsin (דילול 1:500). רקמות הומוגגו ב-4 מעלות צלזיוס באמצעות הומוגנית רקמות של Precellys (Precellys, Montigny-le-Bretonneux, צרפת). הומוגניות של רקמות הודגרו על קרח למשך 10 דקות, ולאחר מכן עברו צנטריפוגה ב-9600×g למשך 10 דקות. הסופרנטנט נוסף לחלק העליון של שיפוע סוכרוז שנוצר על ידי BioComp Gradient Master (10% סוכרוזעל גבי 50% סוכרוז ב-20 מ"מ Tris-HCl, 100 מ"מ NaCl, 5 מ"מ MgCl2, ו-50 מ"ג/מ"לcycloheximide) וצנטריפוגה ב-284,000× g למשך שעתיים בשעה 4◦ג. השיפועים נקצרומלמעלה בקציר Biocomp (Biocomp Instruments, Fredericton, NB, קנדה),ותכולת ה-RNA של חלקים ריבוזומליים שנפלטו נוטרה ברציפות באמצעות UVספיגה ב-254 ננומטר.



2.7. מטבולומי

כדי להעריך חתימות מטבולומיות, הוכנו דגימות של בקרת רקע ועכברי db (n=5 לכל קבוצה) באמצעות מערכת MicroLab STAR® האוטומטית (Hamilton Company, Reno, NV, ארה"ב). ניתוח מטבולומי בוצע ב- Metabolon Inc. בקצרה, רקמות כליה קפואות עברו עיבוד לפי שיטת המיצוי הסטנדרטית של Metabolon לפי פרוטוקול החברה. התמציות נותחו ב- Waters ACQUITY Ultra-performance נוזל כרומטוגרפיה (UPLC) עם עמודת C18 (Waters UPLC BEH C18-2.1 × 100 מ"מ, 1.7 מיקרומטר) ו-Thermo Scientific Q-Exactive ברזולוציה גבוהה/מדויק ספקטרומטר מסה המממשק עם מקור יינון אלקטרוספריי מחומם (HESI-II) ומנתח מסה Orbitrap המופעל ברזולוציית מסה של 35,000.

נתונים גולמיים חולצו, זוהו שיא ועובדו QC באמצעות החומרה והתוכנה של Metabolon (Durham, NC, ארה"ב). הפסגות נכמתו באמצעות שטח מתחת לעקומה. ניתוחים סטטיסטיים סטנדרטיים עם מבחן t דו-מדגם של Welch בוצעו ב-ArrayStudio על נתונים שעברו טרנספורמציה ביומן.




2.8. ביטוי גנים דיפרנציאלי, ניתוח נתיבים וסטטיסטיקה

נתוני ביטוי של LMD אנושי ושל רקמת כליות עכבר מרצף RNA נורמלו כמותית וביטוי דיפרנציאלי נקבע באמצעות בדיקה מדויקת ב-edgeR עם ערכי p < 0.05 נחשב מובהק סטטיסטית לאחר Benjamini-Hochberg false תיקון בדיקות מרובות של שיעור גילוי (FDR). עבור דגימות כליות של עכברים שעברו רצף snRNA, ביטוי דיפרנציאלי נקבע באמצעות מבחן Wilcoxon Rank Sum עם ערכי p מותאמים של Bonferroni < 0.05 שנחשבו מובהקים סטטיסטית. עם זאת, ערכי p לא מתוקנים מסופקים בכתב היד. נתוני ביטוי גנים נחקרו עבור מסלולים מועשרים ב- Gene Ontology, Kegg ו-Reactome באמצעות ניתוח רשת טרנסקריפטוגרמה כפי שתואר בעבר [19,24-26]. בקצרה, כל הגנים המקודדים לחלבון, ללא קשר לערך p, מסודרים לפי רשתות אינטראקציה בין חלבון לחלבון. אלגוריתם מונטה קרלו שימש לאיסוף גנים לפי פונקציות ביולוגיות משותפות. המובהקות נקבעת באמצעות מבחן t דו זנבי של Welch לאחר 500 תמורות אקראיות כדי להעריך את שיעור הגילוי השגוי עם ערכי p < 0.05 הנחשבים מובהקים סטטיסטית. מכיוון שרק תת-קבוצה של גנים באה לידי ביטוי בכל תא נתון עבור נתוני snRNA-seq, שיטת ה-transcriptogram אינה אפשרית עבור סוג נתונים זה. ניתוח נתיב עבור snRNA-seq כלל העשרה של גנים המובעים בצורה דיפרנציאלית באונטולוגיה של גנים, קג וראקטום על פי הבדיקה המדויקת של פישר, מבלי לקחת בחשבון את האינטראקציות של רשת חלבון וחלבון. כל צבירי התאים הצינוריים הפרוקסימליים בכל נקודות הזמן מוזגו כדי לזהות גנים מועשרים בין עכברי בקרה ו-db. גם אשכול הפודוציטים הבודד אוחד על פני נקודות זמן. שיטה זו יושמה על נתוני RNA-seq אנושיים גם בעת הערכת חפיפה עם מערך הנתונים של snRNA-seq. ערכי p מתוקנים ב-FDR עבור מסלולים שימשו לזיהוי חפיפה בין מקטעי משנה תואמים. PAM (מחלקים מסביבmedoids) שימש לאיסוף דגימות סוכרת אנושיות בהתאם לגנים המעורבים בכל מסלול חופף. מאפיינים קליניים דיפרנציאליים נקבעו באמצעות בדיקת Wilcoxon Rank Sum בין קבוצות מקובצות.














אולי גם תרצה