הפרוטוקול מספק הוראות לשינוי RNA עם דימתיל סולפט לניסויים בפרופיל מוטציות. הוא כולל בדיקה במבחנה ו – in vivo עם שתי שיטות הכנה חלופיות לספרייה.
תפקידו של מבנה הרנ”א כמעט בכל תהליך ביולוגי נעשה ברור יותר ויותר, במיוחד בעשור האחרון. עם זאת, גישות קלאסיות לפתרון מבנה RNA, כגון קריסטלוגרפיה של RNA או cryo-EM, לא הצליחו לעמוד בקצב של התחום המתפתח במהירות והצורך בפתרונות בעלי תפוקה גבוהה. פרופיל מוטציה עם ריצוף באמצעות דימתיל סולפט (DMS) MaPseq היא גישה מבוססת ריצוף כדי להסיק את מבנה הרנ”א מתגובתיות הבסיס עם DMS. DMS מתילאט את החנקן N1 באדנוסינים ואת N3 בציטוזינים על פני ווטסון-קריק שלהם כאשר הבסיס אינו מזווג. תמלול לאחור של הרנ”א המתוקן עם הטרנסקריפטאז ההפוך של הקבוצה התרמוספטיבית II (TGIRT-III) מוביל לכך שהבסיסים שעברו מתילציה משולבים כמוטציות ב-cDNA. כאשר מריצפים את ה-cDNA המתקבל וממפים אותו בחזרה לתעתיק ייחוס, שיעורי המוטציה היחסיים עבור כל בסיס מעידים על “מצבו” של הבסיס כמזווג או לא מזווג. למרות שלתגובתיות DMS יש יחס אות לרעש גבוה הן במבחנה והן בתאים , שיטה זו רגישה להטיה בהליכי הטיפול. כדי להפחית את ההטיה הזו, מאמר זה מספק פרוטוקול לטיפול ב-RNA עם DMS בתאים ועם רנ”א מתועתק במבחנה .
מאז הגילוי שלרנ”א יש תכונותמבניות של 1,2 וגם תכונותקטליטיות 3, נחשפו בהדרגה חשיבותו של הרנ”א ותפקודו הרגולטורי במגוון רחב של תהליכים ביולוגיים. ואכן, ההשפעה של מבנה הרנ”א על ויסות גנים זכתה לתשומת לב גוברת4. בדומה לחלבונים, לרנ”א יש מבנים ראשוניים, משניים ושלישוניים, המתייחסים לרצף הנוקלאוטידים, למיפוי הדו-ממדי של אינטראקציות זיווג-בסיס, ולקיפול התלת-ממדי של מבנים אלה המזווגים בבסיס, בהתאמה. בעוד שקביעת המבנה השלישוני היא המפתח להבנת המנגנונים המדויקים העומדים מאחורי תהליכים תלויי RNA, המבנה המשני הוא גם אינפורמטיבי מאוד לגבי תפקוד הרנ”א והוא הבסיס לקיפול תלת-ממדי נוסף5.
עם זאת, קביעת מבנה הרנ”א הייתה מאתגרת במהותה בגישות קונבנציונליות. בעוד שעבור חלבונים, קריסטלוגרפיה, תהודה מגנטית גרעינית (NMR) ומיקרוסקופיית אלקטרונים קריוגנית (cryo-EM) אפשרו לקבוע את מגוון המוטיבים המבניים, מה שמאפשר חיזוי מבנה מהרצף בלבד6, גישות אלה אינן ישימות באופן נרחב על רנ”א. ואכן, רנ”א הן מולקולות גמישות עם אבני בניין (נוקלאוטידים) שיש להן הרבה יותר חופש קונפורמציה וסיבובי בהשוואה למקבילותיהן בחומצות אמינו. יתר על כן, האינטראקציות באמצעות זיווג בסיסים הן דינמיות ורב-תכליתיות יותר מאלו של שאריות חומצות אמינו. כתוצאה מכך, גישות קלאסיות הצליחו רק עבור רנ”א קטנים יחסית עם מבנים מוגדרים היטב וקומפקטיים מאוד7.
גישה נוספת לקביעת מבנה הרנ”א היא באמצעות בדיקה כימית בשילוב עם ריצוף מהדור הבא (NGS). אסטרטגיה זו מייצרת מידע על מצב הקשירה של כל בסיס ברצף RNA (כלומר, המבנה המשני שלו). בקצרה, הבסיסים במולקולת RNA שאינם עוסקים בזיווג בסיסים משתנים באופן דיפרנציאלי על ידי תרכובות כימיות קטנות. תמלול לאחור של רנ”א אלה באמצעות תעתיק הפוך מיוחד (RTs) משלב את השינויים בחומצה דאוקסיריבונוקלאית משלימה (cDNA) כמוטציות. מולקולות cDNA אלה מוגברות לאחר מכן על ידי תגובת שרשרת פולימראז (PCR) ומרוצפות. כדי לקבל מידע על “מצבם” כמאוגדים או לא מאוגדים, תדרי המוטציות בכל בסיס ברנ”א של עניין מחושבים ומוכנסים לתוכנת חיזוי מבנה כאילוצים8. בהתבסס על כללי השכן הקרוב ביותר9 וחישובי אנרגיה חופשית מינימליים 10, תוכנה זו מייצרת מודלים מבניים המתאימים ביותר לנתונים הניסיוניים שהתקבלו11,12.
DMS-MaPseq משתמש ב-DMS, אשר מתילציה של חנקן N1 באדנוסינים וחנקן N3 בציטוזינים על פני ווטסון-קריק שלהם באופן ספציפי מאוד13. שימוש בטרנסקריפטאז הפוך (TGIRT-III) מקבוצה תרמוספטיבית II (TGIRT-III) בתעתיק הפוך יוצר פרופילים מוטציוניים עם יחסי אות לרעש חסרי תקדים, ואף מאפשר פירוק של פרופילים חופפים הנוצרים על ידי שני קונפורמציות חלופיות או יותר14,15. יתר על כן, DMS יכול לחדור לקרומי תאים ולרקמות שלמות, מה שהופך את הבדיקה בהקשרים פיזיולוגיים לאפשרית. עם זאת, יצירת נתונים באיכות טובה היא מאתגרת, מכיוון ששינויים בהליך הטיפול יכולים להשפיע על התוצאות. לכן, אנו מספקים פרוטוקול מפורט הן למבחנה והן ל-DMS-MaPseq בתוך התא כדי להפחית את ההטיה ולהנחות מצטרפים חדשים לשיטה דרך הקשיים שהם עלולים להיתקל בהם. במיוחד לאור מגיפת SARS-CoV2 האחרונה, נתונים באיכות גבוהה על נגיפי RNA הם כלי חשוב לחקר ביטוי גנים ומציאת טיפולים אפשריים.
הפרוטוקול כאן מתאר כיצד לחקור RNA במבחנה ובתאים באמצעות ניסויי פרופיל מוטציה של DMS. יתר על כן, הוא נותן הוראות כיצד להכין ספריות לריצוף Illumina כדי ליצור נתונים ספציפיים לגנים ולנתח את קבצי ה- .fastq המתקבלים. בנוסף, ניתן להשתמש בגישות ספרייה כלל-גנומית. עם זאת, RT-PCR ספציפי לגן מפיק את הנתונים האיכותיים והחזקים ביותר. לכן, אם משווים בין דוגמאות, חשוב לוודא שהם מוכנים עם אסטרטגיות רצף זהות, שכן יצירת הספרייה גורמת להטיה מסוימת. תמיד יש למדוד את יכולת השכפול באמצעות משכפלים.
מספר אמצעי זהירות
רנ”א היא מולקולה לא יציבה הרגישה להתפרקות הן באמצעות טמפרטורות גבוהות והן על ידי RNases. לכן, מומלץ להשתמש באמצעים מיוחדים – שימוש בציוד מגן אישי (PPE), חומר נטול RNAse ומעכבי RNAse. והכי חשוב, RNA צריך להישמר על קרח במידת האפשר. זה חל במיוחד על RNA מתילציה, אשר רגיש עוד יותר לטמפרטורות גבוהות.
חשוב לוודא שמבנה הרנ”א אינו רגיש לריכוז ה-DMS ולתנאי החיץ. מאגרים כגון 100 mM Tris, 100 mM MOPS ו- 100 mM HEPES ב- pH 7-7.5 נותנים אות גבוה אך עשויים שלא להספיק כדי לשמור על ה- pH במהלך התגובה21. כאשר DMS עובר הידרוליזה במים, מה שמפחית את ה-pH, חיץ חזק הוא קריטי לשמירה על pH נייטרלי במהלך תגובת השינוי. התוספת של bicine הוכחה כמסייעת לשמור על ה- pH כבסיסי מעט21 אך גורמת לשינוי DMS נמוך ב- Gs ו- Us, אשר יכול להיות אינפורמטיבי אך יש לנתח אותו בנפרד בשל ייצור אות נמוך בהרבה מאשר As ו- Cs ואינו נדון בהמשך פרוטוקול זה.
ב-RT-PCR ספציפי לגנים, הרנ”א המהונדס מתועתק לאחור לתוך הדנ”א ומוגבר בשברים על ידי PCR. בעוד שגודל הרנ”א יכול תיאורטית להיות בלתי מוגבל, מקטעי PCR אלה לא יעלו על אורך של 400-500 זוגות בסיסים (bp) כדי למנוע הטיה במהלך תגובת השעתוק ההפוכה. באופן אידיאלי, השברים צריכים להיות בטווח של הרצת הרצף (כלומר, אם הרצף מתבצע באמצעות תוכנית ריצוף של 150 x 150 מחזור, מקטע יחיד לא יעלה על 300 bp). בעת שימוש בתוכניות ריצוף עם פחות מחזורים, ניתן לפצל את מוצרי ה-PCR באמצעות dsDNase. יתר על כן, מכיוון שרצפים בתוך רצפי הפריימר אינם מכילים מידע מבני כלשהו, השברים חייבים לחפוף כאשר הרנ”א הנחקר כולל מקטע >1. תגובות RT יכולות להכיל פריימרים מרובים של RT עבור מקטעים שונים (עד 10 פריימרים שונים של RT). בהתאם לרצפים, איגום פריימרים RT יכול להפוך את השעתוק ההפוך לפחות יעיל, אך בדרך כלל עובד היטב. כל תגובת PCR צריכה להתבצע בנפרד.
כאשר חוקרים RNA עם DMS, תנאי הניסוי ממלאים תפקיד נוסף, שכן רנ”א רבים אינם יציבים מבחינה תרמודינמית ומשנים את הקונפורמציה שלהם בהתבסס על גורמים סביבתיים כגון טמפרטורה. כדי למנוע אי סדרים, יש לשמור על תנאי הניסוי קבועים ככל האפשר, גם ביחס לזמני התגובה. נראה כי תנאי החיץ ניתנים להחלפה במידה מסוימת 17,20,22,23 כאשר התנאים הבסיסיים נשמרים – יכולת האגירה והנוכחות של יונים חד-ערכיים (Na) ודיוולנטיים (Mg) – כדי להבטיח קיפול תקין של RNA 24.
לגבי הכנת הספרייה של רנ”א שונה, יש לקחת בחשבון מספר היבטים. ראשית, כפי שצוין קודם לכן, רנ”א שעברו שינוי הם פחות יציבים ממקביליהם שלא שונו, כלומר הם עשויים לדרוש אופטימיזציה של זמני הפיצול להתפלגות אופטימלית של גודל השבר. יתר על כן, ערכות הכנה מסוימות של ספריית RNA, כמו גם גישות RNAseq רבות אחרות, משתמשות בפריימרים אקראיים בערכת השעתוק ההפוכה. זה עלול להוביל לכיסוי נמוך יותר של הייחוס, במיוחד ב-3′ של גן, ובסופו של דבר, לעומק כיסוי לא מספיק. אם הכיסוי של אזור מסוים נמוך מדי, ייתכן שיהיה צורך להסיר בסיסים אלה מחיזוי המבנה. מלבד ערכות RT-PCR ו-RNAseq של גנום שלם, ניתן להשתמש בגישות אחרות להכנת ספרייה. פרוטוקולים הכוללים קשירת מתאמי 3′ ו/או 5′ ל-RNA הם יתרון כאשר משתמשים בשברי RNA קטנים או כאשר יש להימנע מאובדן מידע בדיקה באזורי הפריימר.
לבסוף, יש לפרש תמיד בזהירות את ניתוח הניסויים הכימיים. נכון לעכשיו, אין תוכנה המנבאת את מבנה הרנ”א של כל רנ”א מהרצף לבדו בדיוק גבוה. למרות שאילוצים של בדיקה כימית משפרים מאוד את הדיוק, יצירת מודלים טובים עבור רנ”א ארוך (>500 nt) עדיין מאתגרת. מודלים אלה צריכים להיבדק עוד יותר על ידי גישות אחרות ו/או מוטגנזה. תוכנת חיזוי RNA מבצעת אופטימיזציה למספר המרבי של זוגות בסיסים, ובכך מענישה באופן משמעותי קונפורמציות פתוחות, שאולי אינן מייצגות במדויק קיפול RNA5. לפיכך, יש לבחון את מודל המבנה המתקבל על ידי כימות הסכם החיזוי עם נתוני הבדיקה הכימיים הבסיסיים (למשל, על ידי AUROC) ובין שכפולים (למשל, על ידי mFMI), כפי שהודגם על ידי Lan et al.20.
באופן אידיאלי, יש להשתמש במספר ניסויים במערכות שונות כדי לאתגר את מודל המבנה המתקבל כדי לחזק את ההשערה. אלה יכולים לכלול שימוש בגישות חוץ גופיות ובתוך התא, מוטציות מפצות וקווי תאים ומינים שונים. יתר על כן, תגובות גולמיות הן לעתים קרובות אינפורמטיביות באותה מידה או אפילו יותר מאשר תחזיות מבנה, מכיוון שהן מתעדות את תמונת “האמת הקרקעית” של אנסמבל קיפול הרנ”א. ככזה, תגובתיות גולמית מתאימה מאוד ואינפורמטיבית להשוואת שינויי מבנה בין תנאים שונים. חשוב לציין שמבני האנרגיה החופשית הנמוכים ביותר המחושבים באמצעות אילוצי בדיקה כימית עם חיזוי חישובי צריכים לשמש רק כהשערת מוצא לקראת מודל מבנה שלם.
The authors have nothing to disclose.
ללא
1 Kb Plus DNA Ladder | 10787018 | Thermo | |
2-mercaptoethanol | M6250-250ML | Sigma | |
Acid-Phenol:Chloroform, pH 4.5 | AM9720 | Thermo | |
Advantage PCR | 639206 | Takara | |
CloneAmp HiFi PCR Premix | 639298 | Takara | |
DMS | D186309 |
Sigma | |
dNTPs 10 mM each | U151B | Promega | |
E-Gel EX Agarose Gels, 2% | G402022 | Thermo | precast agarose gels |
Ethanol (200 proof) | E7023-4X4L | Sigma | |
Falcon tubes, 15 mL, 50 mL | |||
GlycoBlue | co-precipitant | ||
HCT-8 cells | ATCC #CCL-244 | ||
Invitrogen MgCl2 (1 M) | AM9530G | fisherscientific | |
Isopropanol | 278475 | Sigma | |
Megascript T7 transcription | AM1334 | Thermo | |
NanoDrop spectrophotometer | |||
Novex TBE Gels, 8%, 10 well | EC6215BOX | Thermo | |
OC43 | ATCC #VR-1558 | ||
RiboRuler Low Range RNA Ladder | SM1831 | Thermo | |
RNAse H | M0297L | NEB | |
Sodium Cacodylate, 0.4 M, pH 7.2 | 102090-964 | VWR | |
Sodium hydroxide solution | S8263-150ML | Sigma | |
SuperScript II Reverse Transcriptase for FSB and DTT | 18064014 | Thermo | |
TGIRT-III Enzyme | TGIRT50 | Ingex | |
The Oligo Clean & Concentrator | D4060 | Genesee | |
The RNA Clean & Concentrator kits are RNA clean up kits | R1016 | Genesee | |
TRIzol Reagents | 15596018 | Thermo | RNA isolation reagent |
Water, (For RNA Work) (DEPC-Treated, DNASE, RNASE free/Mol. Biol.) | BP561-1 | fisherscientific | |
xGen Broad-range RNA Library Prep 16rxn | 10009865 | IDT | |
Zymo RNA clean and concentrator columns |