חקר יעדים לאי ספיקת לב באמצעות רנדומיזציה מנדלית, קו-לוקליזציה, רנדומיזציה מנדלית מבוססת נתוני סיכום, וחיזוי תרופות

Huiling Zhu; Suxin Luo

doi:10.3791/71544

Method Article

חקר יעדים לאי ספיקת לב באמצעות רנדומיזציה מנדלית, קו-לוקליזציה, רנדומיזציה מנדלית מבוססת נתוני סיכום, וחיזוי תרופות

DOI:

10.3791/71544

⸱

July 10th, 2026

Huiling Zhu¹ , Suxin Luo¹

¹Department of Cardiovascular Medicine, Cardiovascular Research Center, The First Affiliated Hospital of Chongqing Medical University

Summary

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

Abstract

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

אי ספיקת לב (HF) היא מחלת לב ושט רווחת הפוגעת באופן משמעותי באיכות החיים בשלביה המתקדמים. למרות מגוון אסטרטגיות טיפוליות עכשוויות, הנטל המחלתי של אי ספיקת לב נותר משמעותי. לכן, חקירת מטרות טיפוליות חדשות נחוצה בדחיפות. מחקר זה השתמש ברנדומיזציה מנדלית (MR) כדי להעריך את היחסים הסיבתיים בין גנים ניתנים לטיפול ואי ספיקת לב. ניתוח קו-לוקליזציה ורנדומיזציה מנדלית מבוססת נתוני סיכום (SMR) נערכו כדי לאמת עוד יותר את הקשר ביניהם. לבסוף, זוהו 8 מטרות טיפוליות פוטנציאליות לאי ספיקת לב, כולל 4 גורמי סיכון (CYP11A1, GALT, KCNH2, ו-METRN) ו-4 גורמי הגנה (APOM, CHD4, IL11RA, ו-LPAR5). ניתוח העשרת GO, העשרת KEGG ורשת אינטראקציות חלבון-חלבון (PPI) הצביעו על כך שהגנים הללו מעורבים בעיקר בוויסות מטבולי ובתהליכים אלקטרופיזיולוגיים לבביים. יתר על כן, תרופות פוטנציאליות המכוונות למטרות אלו זוהו באמצעות חיזוי תרופות ודוקינג מולקולרי, כולל מיטונאן, אדהל, בנזופוראנים, 1,3-די-או-טולילגואנידין, 96-69-5 וברומונול לקטון. יתר על כן, תוצאות ה-PCR הראו גם שמיטונאן עשוי להיות קשור ל-CYP11A1 ו-KCNH2. תרופות המתוכננות על סמך מטרות טיפוליות פוטנציאליות אלו עשויות להציע שיעורי הצלחה גבוהים יותר ולשפר את התוצאות הקליניות.

Introduction

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

אי ספיקת לב (HF) הוא תסמונת לב וכלי דם נפוצה המוגדרת על ידי מילוי חדרי לקוי או התכווצות שריר הלב, וכתוצאה מכך תפוקת לב לא מספקת כדי לעמוד בדרישות המטבוליות¹^,²^,³. מבחינה קלינית, אי ספיקת לב מאופיינת בגודש ריאתי ומערכתי, היפופרפוזיה של רקמות ותסמינים כמו קוצר נשימה, עייפות ובצקות היקפיות¹^,⁴^,⁵. למרות ההתקדמות הטיפולית, אי ספיקת לב נותרה גורם מוביל לתחלואה ותמותה ברחבי העולם, ומשפיעה על למעלה מ-64 מיליון אנשים על פי מחקר העומס הגלובלי של המחלות (GBD)⁶^,⁷^,⁸. נטל המחלה הגבוה מדגיש את הצורך הדחוף באסטרטגיות טיפוליות חדשות.

ההתערבויות התרופתיות והמבוססות על מכשירים הנוכחיות לאי ספיקת לב מוגבלות על ידי יעילות ותגובות לוואי תת-אופטימליות, ותורמות לשיעורי אשפוז חוזרים גבוהים ולאיכות חיים לקויה⁹^,¹⁰^,¹¹^,¹². לכן, זיהוי מטרות מולקולריות חדשות עם קישורים סיבתיים לפתוגנזה של אי ספיקת לב הוא קריטי לשיפור התוצאות הקליניות. שילוב גנומיקה בגילוי תרופות מייצג אסטרטגיה יעילה מאוד לשיפור היעילות, מכיוון שטיפולים שאושרו גנטית מדגימים שיעור הצלחה גבוה משמעותית בניסויים קליניים¹³^,¹⁴^,¹⁵. יתרה מכך, חלבונים המקודדים על ידי גנים פעילים אלו משמשים כמטרות טיפוליות מרכזיות הן למולקולות קטנות והן לנוגדנים מונוקלונליים¹⁶^,¹⁷.

התקדמות חדשה בגנטיקה האנושית הקלה על עדיפות מטרות תרופות באמצעות ניתוח אינטגרטיבי של מחקרי קישור גנומי (GWAS) ונתוני לוקוס תכונות קוונטית של ביטוי (eQTL). רנדומיזציה מנדלית (MR) מנצלת גרסאות גנטיות כמשתנים אינסטרומנטליים כדי להעריך השפעות סיבתיות, וממזערת בלבול וסיבתיות הפוכה הטבועים במחקרים תצפיתיים¹³^,¹⁴. בשילוב עם קולוקליזציה ונתוני MR מבוססי סיכום (SMR), גישה זו מאפשרת עדיפות חזקה של גנים שאותם אפשר להשפיע על סיכון למחלה, ומציעה אסטרטגיה מעוגנת גנטית לגילוי מטרות¹⁸^,¹⁹^,²⁰. בהשוואה למחקרי קישור קונבנציונליים, שילוב של MR, SMR וקולוקליזציה מתגבר על בלבול וסיבתיות הפוכה כדי להקים מסקנות סיבתיות חזקות. יתרה מכך, מסגרת זו מסננת אותות שווא מתחום השוויון על ידי הבטחה שביטוי הגנים וסיכון לאי ספיקת לב חולקים את אותו משתנה סיבתי בסיסי, ומייעלת את עדיפות המטרות. ראוי לציין שבעוד שרקמות לב וכלי דם מציגות ספציפיות גבוהה, מחקר זה משתמש בעיקר בנתוני eQTL מרחבי גנום מדם בקנה מידה גדול כדי למקסם את הכוח הסטטיסטי. למרות שגישה זו עשויה שלא ללכוד באופן מלא את המנגנונים הספציפיים לרקמות הטבועים באי ספיקת לב, גרסאות שמקורן בדם משמשים כסמנים ביולוגיים נגישים ומספקים תובנות חזקות לגבי הגורמים הגנטיים המערכתיים של המחלה.

במחקר זה, נעשה שימוש במחקרי GWAS ובנתוני eQTL בקנה מידה גדול כדי להעריך את ההשפעות הסיבתיות של גנים שניתן לטפל בהם על סיכון לאי ספיקת לב. אסוציאציות משמעותיות אושרו עוד יותר באמצעות ניתוחי קולוקליזציה ו-SMR. כדי להבהיר את הפונקציות הביולוגיות של המטרות המועדפות, בוצעו ניתוחי העשרה של Gene Ontology (GO) ו-Kyoto Encyclopedia of Genes and Genomes (KEGG), יחד עם ניתוח רשת אינטראקציות חלבון-חלבון (PPI). הפוטנציאל הטיפולי הוערך לאחר מכן באמצעות חיזוי ת

Protocol

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

1. עיצוב המחקר

תהליך העבודה של מחקר זה מוצג באיור 1. ראשית, נתוני חשיפה הושגו מ-cis-eQTLs של גנים ניתנים לטיפול. שנית, ניתוח MR נערך כדי לחקור את הקשר הסיבתי בין גנים ניתנים לטיפול לבין HF. מטרות טיפוליות מועמדות ל-HF זוהו לאחר מכן על ידי ניתוחי קולוקליזציה ו-SMR. לאחר מכן בוצעו ניתוח העשרה של גנים ובניית רשת PPI. לבסוף, בוצעו חיזוי תרופות ו-docking מולקולרי כדי לחקור את הזיקה של תרופות לגנים.

2. מקור הנתונים

בסך הכל זוהו 6,888 גנים ניתנים לטיפול משתי מקורות. מערך נתונים אחד כלל 4,463 גנים ממחקר שפורסם בעבר ששילב מקורות נתונים מרובים¹⁷. מערך הנתונים השני הכיל 5,012 גנים ממאגר האינטראקציות בין תרופות לגנים (DGIdb)²¹. מידע מפורט מוצג באתר (https://www.dgidb.org/). לאחר מיזוג והסרת כפילויות, נשמרו 6,888 גנים ניתנים לטיפול ייחודיים לניתוח.

נתוני cis-eQTL הושגו ממאגר eQTLGen (https://eqtlgen.org/cis-eqtls.html), שכולל פרופילי ביטוי מ-31,684 דגימות דם ומכסה 16,987 גנים²².

נתוני סיכום HF הושגו מהמהדורה R12 של מאגר FinnGen²³, שכלל 37,653 מקרים ו-462,695 בקרות. האבחון של HF התבסס על סיווג המחלות הבינלאומי (ICD), כולל ICD-10—I11.0, I13.0, I13.2, I50, ICD-9—4029B|428, ו-ICD-8—42700|42710|428|7824. כל המשתתפים היו ממוצא אירופי. מידע מפורט זמין באתר FinnGen והמאמר המקורי (https://www.finngen.fi/en/access_results).

3. ניתוח רנדומיזציה מנדלית

משתנים אינסטרומנטליים (IV) נבחרו על סמך שלושת ההנחות הליבה של MR²⁴. ראשית, IV מקושרים ישירות עם גורמי חשיפה. שנית, IV אינם מקושרים עם גורמי בלבול. שלישית, IV משפיעים על התוצאה באמצעות חשיפה בלבד. בהתחשב בכך ש-eQTLs פרוקסימליים מפעילים שליטה רגולטורית מיידית יותר וחזקה יותר על גנים מטרה, MR השתמשה אך ורק ב-SNPs שנמצאים בתוך אזורים של 100 kb של הגנים כדי למקסם את הקירבה הביולוגית. הועסקו קריטריונים קפדניים לבחירת SNP. כל ה-SNPs שנבחרו היו ממוקמים בטווח של 100 kb במעלה הזרם של אתר התחלת התעתיק ו-100 kb במורד הזרם של אתר סיום התעתיק של כל גן נייטרלי. נכללו רק SNPs עם משמעות גנומית (p < 5×10^-8)²⁵. כדי למזער את אי-שיווי משקל קישור (LD), בוצעה קלומפינג LD באמצעות סף r² של 0.001 ומרחק קלומפינג של 10,000 kb²⁶. נשמרו SNPs עם סטטיסטיקת F >10, שחושבו באמצעות הנוסחה: F = [R²(N-k-1)]/[k(1-R²)]²⁷. מידע על SNP מסופק בטבלה משלמת 1.

שיטות הניתוח העיקריות היו משוקללות הפך ההופכי (IVW) ויחס וולד. ערך p פחות מ-0.05 נחשב למשמעותי סטטיסטית²⁸^,²⁹. שיעור תגלית שווא (FDR) הוחל כדי להתאים את ערך ה-p כדי למזער תוצאות חיוביות שווא³⁰. כאשר היה זמין רק SNP אחד, נעשה שימוש בשיטת יחס וולד; אחרת, IVW הוחל. הניתוחים כללו MR-Egger, מדיה משוקללת, מצב משוקלל ומצב פשוט כאשר היו מספיק SNPs.

נעשה שימוש במבחן Cochran's Q ובניתוח MR-Egger intercept כדי לבחון את ההטרוגניות והפליאוטרופיה. ערך Q-p גדול מ-0.05 ה

Results

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

תורשה מנדלית
תוצאות ה-MR להערכת ההשפעות הסיבתיות בין ביטוי גנים ניתנים לתרופה ו-HF הוצגו בטבלה משלמת 3. ערך P מותאם ל-FDR פחות מ-0.05 נחשב למשמעותי מבחינה סטטיסטית. בסך הכל זוהו 11 גנים עם קשרים סיבתיים משמעותיים עם HF (איור 2). ביניהם, שישה גנים היו קשורים לסיכון מוגבר ל-HF, כולל CDKN1A (P_FDR: 0.03, OR: 1.293), CYP11A1 (P_FDR: 0.034, OR: 1.494), GALT (P_FDR: 0.042, OR: 1.074), KCNH2 (P_FDR: 0.018, OR: 1.163), LST1 (P_FDR: 0.03, OR: 1.112) ו-METRN (P_FDR: 0.018, OR: 1.115). חמשת הגנים הנותרים היו קשורים לסיכון מופחת ל-HF, כולל APOM (P_FDR: 0.03, OR: 0.773), CAMK2G (P_FDR: 0.018, OR: 0.699), CHD4 (P_FDR: 0.018, OR: 0.811), IL11RA (P_FDR: 0.018, OR: 0.947), ו-LPAR5 (P_FDR: 0.021, OR: 0.915).

ניתוח רגישות בוצע במחקר זה, והתוצאות הראו כי לא זוהתה הטרוגניות משמעותית או פליוטרופיה אופקית (טבלה 1), וההערכות של MR לא הונעו על ידי SNP בודד (איור משלים 1).

ניתוח קו-לוקליזציה
ניתוח קו-לוקליזציה בוצע על 11 הגנים שזוהו (איור 3 וטבלה 2). תשעה גנים הדגימו עדות חזקה לשיתוף גרסאות סיבתיות עם HF, כולל APOM (PP.H4: 0.923), CHD4 (PP.H4: 0.942), CYP11A1 (PP.H4: 0.91), GALT (PP.H4: 0.967), IL11RA (PP.H4: 0.967), KCNH2 (PP.H4: 0.815), LPAR5 (PP.H4: 0.938), LST1 (PP.H4: 0.913) ו-METRN (PP.H4: 0.907). שני גנים הושרו מניתוח נוסף מכיוון ש-PP.H4 היה פחות מ-0.8, כלומר, CDKN1A ו-CAMK2G.

ניתוח SMR
כדי לאמת עוד יותר את קשרי הגן-מחלה, בוצע ניתוח SMR עבור תשעת הגנים המשותפים (טבלה 2). ניתוח ה-SMR הראה ששמונה גנים נותרו קשורים משמעותית ל-HF. היו ארבעה גנים הקשורים באופן שלילי ל-HF, כולל APOM (P_SMR: 0.018, P_HEIDI: 0.115, OR: 0.838), CHD4 (P_SMR: 5.45E-05, P_HEIDI: 0.458, OR: 0.732), IL11RA (P_SMR: 1.00E-04, P_HEIDI: 0.388, OR: 0.938), ו-LPAR5 (P_SMR: 3.71E-05, P_HEIDI: 0.168, OR: 0.889). וארבעה היו קשורים באופן חיובי ל-HF, כולל CYP11A1 (P_SMR: 5.68E-04, P_HEIDI: 0.056, OR: 1.81), GALT (P_SMR: 1.02E-04, P_HEIDI: 0.726, OR: 1.085), KCNH2 (P_SMR: 2.76E-04, P_HEIDI: 0.162, OR: 1.236) ו-METRN (P_SMR: 1.55E-05, P_HEIDI: 0.068, OR: 1.103). למרות ש-LST1 (P_SMR: 4.61E-05, P_HEIDI: 0.025, OR: 1.135) הראה קשר עם סיכון מוגבר ל-HF, מבחן HEIDI הצביע על כך שהקשר עשוי להיות נגרם על ידי שיווי משקל קישור. לכן, LST1 הוצא.

ניתוח העשרה
שמונת הגנים שזוהו באמצעות ניתוחי MR, קו-לוקליזציה ו-SMR עברו ניתוח העשרה GO וניתוח העשרה KEGG (איור 4A–E). בקטגוריית BP, הגנים היו מועשרים בפירוק תהליכים מטבוליים (תהליך מטבולי של קורטיזול, תהליך מטבולי של גלקטוז, תהליך ביו-סינתטי של גלוקוקורטיקואידים ותהליך הובלת כולסטרול הפוכה), תהליכים הקשורים לחלקיקי ליפופרוטאין (הרכבת חלקיקי ליפופרוטאין

Discussion

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

מחקר זה העריך באופן שיטתי את ההשפעות הסיבתיות של גנים ניתנים לתרופות על סיכון ל-HF באמצעות MR, וזיהה שמונה מטרות טיפוליות פוטנציאליות. אלו כללו ארבעה גורמי הגנה (APOM, CHD4, IL11RA, ו-LPAR5) וארבעה גורמי סיכון (CYP11A1, GALT, KCNH2, ו-METRN). קשרים אלו נתמכו עוד על ידי ניתוח מיקום משותף ו-SMR. אנוטציה פונקציונלית באמצעות GO, העשרה ב-KEGG, וניתוח רשת PPI הציעו שגנים אלו עשויים להשפיע על HF בעיקר דרך מסלולים הקשורים למטבוליזם ואלקטרופיזיולוגיה. בנוסף, חיזוי תרופות ומחקרי עגינה מולקולרית הדגישו את הפוטנציאל התרגומי של מטרות אלו.

APOM זוהה כגורם הגנה ל-HF. APOM קושר ספינגוזין-1-פוספט (S1P) ומפעיל השפעות אנטי-דלקתיות על ידי הפחתת ביטוי TNF-α, הגבלת הדבקה של מונוציטים-אנדותל, ושמירה על שלמות המחסום האנדותלי, ובכך מפחית נזק מיאוקרדי שנגרם על ידי דלקת⁴²^,⁴³^,⁴⁴. בנוסף, APOM תומך בהישרדות וריבוי קרדיומיוציטים תוך הפחתת אפופטוזיס, ושמירה על מבנה ותפקוד הלב. חסרונו קשור לתוצאות לבביות שליליות, תוך התאמה לממצא זה⁴²^,⁴⁵. עם זאת, המנגנון הבסיסי המקשר בין APOM ל-HF דורש חקירה נוספת.

CHD4, גן הגנה נוסף שזוהה, הוא מרכיב ליבה של קומפלקס הנוקלאוזום ועיצוב דאוקסילאז (NuRD) ומשחק תפקיד חיוני בהתפתחות לבבית וצמיחת מיוקרד⁴⁶^,⁴⁷^,⁴⁸. אובדן CHD4 מוביל לשילוב של חלבוני מיופיברילים לא נכונים לתוך הסרקומרים הלבביים, וכתוצאה מכך חריגות מבניות ופונקציונליות, כולל התרחבות חדרים ופגיעה בהתפתחות חדרית. תוצאה זו מחזקת את התפקיד הקריטי של CHD4 בשמירה על תפקוד לבבי תקין.

גם נמצא ש-IL11RA יכול להפחית את הסיכון ל-HF. IL11RA הוא תת-היחידה הקשורים לליגנד העיקרי של קומפלקס הקולטן IL-11. בעוד ש-IL-11 הוכח כמניע תגובות פיברוטיות בלב, התפקיד הספציפי של IL11RA נותר לא ברור⁴⁹. ההשפעה המגינה שנצפתה במחקר זה מצביעה על כך ש-IL11RA עשוי לממש תפקידים תלויי הקשר או לקיים אינטראקציה עם רכיבי איתות אחרים כדי להגביל פיברוזיס בתנאים מסוימים⁵⁰. הפער הברור הזה מדגיש את הצורך במחקרים מכניסטיים נוספים כדי להבהיר את התפקיד הפונקציונאלי של IL11RA בפתוגנזה של HF.

מחקר זה הראה ש-LPAR5 היה קשור לסיכון מופחת ל-HF. LPAR5 הוא קולטן לחומצה ליזופוספטידית (LPA). למרות ש-LPA נחשב באופן כללי לפרו-דלקתי והוצג כמעורב ברמודלינג לבבי שלילי, לקולטני LPA בודדים עשויים להיות תפקידים מגוונים⁵¹^,⁵². לדוגמה, איתות LPA מווסת את היווצרות סיבי הלחץ הציטושלדיים דרך מסלולי Gα12/13 ו-Gαq/11, וממלא תפקיד חשוב בשמירה על מורפולוגיה, תפקוד חוזי, ותהליך העברת אותות של קרדיומיוציטים⁵³^,⁵⁴. עם זאת, התפקיד המדויק של LPAR5 בפתופיזיולוגיה של HF דורש חקירה נוספת.

חקירה זו הצביעה על כך ש-CYP11A1 היה גורם סיכון ל-HF. CYP11A1 מקודד את האנזים המגביל את הקצב בביוסינתזה של הורמוני סטרואיד⁵⁵^,⁵⁶. ייצור מוגבר של הורמוני סטרואיד, במיוחד אלדוסטרון, תורם לשמירת נתרן ולהעמסה קרדיאלית מוגברת, תכונות מרכזיות בהתקדמות HF⁵⁷^,⁵⁸. בנוסף, CYP11A1 מייצר מינים פעילים של חמצן במהלך פעילותו המטבולית, שעשויים להחמיר את הפגיעה בתאי הל

Disclosures

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

המחברים מצהירים כי אין להם עניינים מתנגשים.

Acknowledgements

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

מוכרים המשתתפים והחוקרים של מחקר FinnGen. קונסורציום eQTLGen, DGIdb וחוקרים אחרים שסיפקו נתונים זמינים לציבור לניתוח זה מוכרים בהכרת תודה. עבודה זו נתמכה על ידי פרויקט חדשנות המחקר של המאסטר של המכללה הקלינית הראשונה של אוניברסיטת רפואית צ'ונגצ'ינג (CYYY-SSCX202516) והאבחון והטיפול הסטנדרטי של אי ספיקת לב (2025cyjstg006).

תרומת הסופר
הוי-לינג ג'ו תרם לתפיסה, רכישת מימון, חקירה, מתודולוגיה, משאבים, תוכנה, פיקוח, אימות, כתיבת הטיוטה המקורית, וביקורת ועריכת המאמר. סוקסין לוו תרם לרכישת מימון, מתודולוגיה, פיקוח וביקורת ועריכת המאמר.

Materials

List of materials used in this article
Name	Company	Catalog Number
2× Universal SYBR Green Fast qPCR Mix	Abclonal	RK21203
AC16	HyCyte	TCH-C119
AC16 cell-specific culture medium	HyCyte	TCH-G119
cDNA reverse kit	Abclonal	RK20400
CFX96™ Real-Time system	Bio-Rad Laboratories
Mitotane	TargetMol	T1199
Palmitic acid	TargetMol	T2908
RNA extraction kit	Abclonal	RK30120

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

חקר יעדים לאי ספיקת לב באמצעות רנדומיזציה מנדלית, קו-לוקליזציה, רנדומיזציה מנדלית מבוססת נתוני סיכום, וחיזוי תרופות

In This Article

Summary

Abstract

Introduction

Protocol

Results

Discussion

Disclosures

Acknowledgements

Materials

Reprints and Permissions

Tags

Related Articles