$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
מיקובקטריום בוביס זן BCG (bacilli de Calmette et Guérin) 1173P2 שעבר גדילה אקספוננציאלית היה נתון לסדרת זמן (0, 4, 10 ו-20 יום) של רעב תזונתי, ואחריו החייאה של 6 ימים בתווך עשיר בחומרים מזינים כפי שהוצג קודם לכן בהו ואחרים.7. רנ"א קטן בודד מתרבית חיידקים, עם שלושה העתקים ביולוגיים, בכל אחת מחמש נקודות הזמן שנקבעו. ספריות Illumina נבנו באמצעות תהליך העבודה של הכנת ספריית AQRNA-seq שתואר לעיל (איור 1), ולאחר מכן ריצוף על רצף במרכז BioMicro של המכון הטכנולוגי של מסצ'וסטס. לאחר מכן נתוני הרצף עובדו באמצעות צינור ניתוח הנתונים AQRNA-seq (איור 2) שהותאם לכימות שפע tRNA.
לאחר הגברת ה-PCR של ספריית cDNA עם פריימרים לריצוף, נצפתה נוכחות של מוצרי PCR בגודל של 175 זוגות בסיסים (bp) בכל הדגימות (איור 3A), מה שמרמז על היווצרות דימרים ראשוניים. כדי למתן את הנשיאה של דימרים פריימרים, מוצרי PCR שגודלם עולה על 195 bp הוצאו מהג'ל וטוהרו (איור 3B).
קריאות רצף מסוננות וגזוזות מופו לספריית רצפי ייחוס מותאמת אישית הכוללת את 45 איזוקבלי tRNA, התקן הפנימי ורצפי בקרה (כלומר, 23S rRNA, 16S rRNA, 5S rRNA, rnpB ו-ssr). איזוקבלטורים של tRNA היוו 10.5% עד 40.2% מכלל הקריאות שמופו של מדגם נתון, והראו שפע גבוה בהרבה מרצפי הבקרה (איור 4). חשוב לציין כי ניתן לייחס את פרופורציות הקריאה הנמוכות יחסית של איזוקבלי tRNA לשפע היחסי הגבוה יותר של התקנים הפנימיים. לכן, פרופורציות הקריאה של איזוקבלי tRNA ביחס לתקנים פנימיים (איור 4, בלוקי צבע ורודים לעומת ירוקים) ניתנות לשליטה על-ידי המפעיל, באמצעות כוונון עדין של כמות התקן הפנימי המוקפץ לתוך התגובה.
נתוני שפע ה-tRNA הגולמיים נורמלו באמצעות שיטת חציון היחסים שיושמה בגרסת חבילת DESeq2 (להלן v) 1.36.016 בסביבת תכנות סטטיסטית R (להלן R) v 4.2.117. לאחר הנורמליזציה, מושג נוף כמותי של איזוקבלטורים של tRNA ב-Mycobacterium bovis BCG במהלך תקופה של רעב והחייאה של חומרים מזינים (איור 5).
כדי לחשוף אשכולות נפרדים של דגימות עם פנוטיפים שונים בהתבסס על דפוסים בשפע האיזוקבילטורים של tRNA, ניתוח רכיבים עיקריים (PCA) בוצע על נתוני שפע tRNA מנורמלים באמצעות חבילת הנתונים הסטטיסטיים v 4.2.117 in R (איור 6). הניתוח הבחין בין דגימות של יום הרעב 0 ויום ההחייאה 6 לבין דגימות של ימי הרעב 4, 10 ו-20, מה שמצביע על הבדל ניכר בנוף ה-tRNA של Mycobacterium bovis BCG שגדל במדיום חסר חומרים מזינים ובמדיום עשיר בחומרים מזינים.
כדי לאפיין את הדינמיקה של השפע של כל איזוקבלור tRNA על פני חמש נקודות הזמן המיועדות, בוצע ניתוח ביטוי דיפרנציאלי על נתוני שפע tRNA מנורמל באמצעות חבילת DESeq2 v 1.36.0 ב- R (איור 7). הניתוח גילה כי 17 מתוך 20 משפחות האיזוקבילטורים הכילו איזוקבלים שהתבטאו באופן דיפרנציאלי (כלומר, מווסתים באופן משמעותי למעלה או למטה) לפחות באחת מנקודות הזמן, מה שמרמז על תפקיד פוטנציאלי של ויסות מאגר tRNA במצב המתמשך של Mycobacterium bovis BCG במהלך שחפת.

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

איור 2: סכמה של צינור ניתוח הנתונים AQRNA-seq. שלבי המפתח המתוארים בקו הצינור מפורטים במרכז הסכימה ומחוברים לאיורים הגרפיים המתאימים שלהם באמצעות קווים מקווקוים. התיאור המפורט של כל שלב זמין ב- GitHub (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git). אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.

איור 3: אלקטרופורזה של ג'ל אגרוז של מקטעי cDNA לאחר הגברה של PCR עם פריימרים לריצוף. (A) תמונה של הג'ל לפני מיצוי וטיהור הג'ל. נתיבים 7 ו-14 מצד שמאל מכילים 5 μL של סולם DNA של 50 bp, ואילו הנתיבים האחרים מכילים 20 μL של כל אחת מ-15 הדגימות. לוקליזציה של גודל מוצרי PCR מצביעה על הריכוז הגבוה ביותר שלהם בטווח שבין 175 bp (דימרים פריימר) ל 300 bp (שני פריימרים + 120 bp 5S rRNA). (B) תמונה של הג'ל לאחר מיצוי וטיהור הג'ל. עבור כל דגימה, בלוק הג'ל בין 200 bp ל 400 bp נכרת כדי למזער את הזיהום של דימרים פריימר בספריית הריצוף. אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.

איור 4: מספר קריאות הרצף שמופו בהצלחה לספריית רצפי הייחוס. ציר ה-x מציג את שמות הדגימות (לדוגמה, D18-69XX) מקובצות לפי נקודת זמן (לדוגמה, יום הרעב 0). עבור כל דגימה, ספירת הקריאה המשויכת לקטגוריות נושא יעד שונות מיוצגת באמצעות גושי צבע המוערמים זה על גבי זה. מספרים הממוקמים במרכז גושי הצבע מייצגים את יחסי הקריאות המתאימים לנושאי היעד המתאימים במדגם נתון. אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.

איור 5: נוף כמותי של איזוקבלי tRNA של Mycobacterium bovis BCG בנקודות זמן שונות לאורך מסלול זמן הרעב וההחייאה. נתוני שפע tRNA גולמיים נורמלו באמצעות שיטת חציון היחסים. כאן, כל שורה מתארת שפע tRNA מנורמל (ציר y) כממוצע ± שגיאת תקן עבור 3 העתקים ביולוגיים בכל נקודת זמן. על ציר ה-x, איזוקבלטורים מאותה משפחה קובצו יחד ותויגו עם חומצת האמינו המתאימה. אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.

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

איור 7: ביטוי דיפרנציאלי של איזוקבלי tRNA על פני נקודות זמן שונות. שפע tRNA מנורמל סוכם כאמצעי (קשרי קו) ± שגיאת תקן (פסי שגיאה) על פני 3 עותקים ביולוגיים. בשל מגבלת מקום, התנאים קוצרו כדלקמן: S0-S20 = ימי רעב 0-20; R6 = יום החייאה 6. ניתוח ביטוי דיפרנציאלי בוצע עבור כל איזוקבלור tRNA, תוך השוואה בין נקודות זמן שונות באופן זוגי באמצעות מבחן יחס הסבירות ומבחן Wald. אותיות קומפקטיות שימשו לייצוג מובהקות סטטיסטית, כאשר השפע של איזוקנטור tRNA נתון בנקודות זמן החולקות לפחות אות משותפת אחת לא היה שונה באופן משמעותי זו מזו. לדוגמה, השפע של tRNA-Lys-CTT-1-1 (בלוח ליזין) היה מווסת באופן משמעותי מ-S0 ל-S4 ומ-S4 ל-S10, אך לא מ-S10 ל-S20. לאחר מכן הוא הוסדר משמעותית מ-S20 ל-R6. אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.
טבלה 1: אוליגונוקלאוטידים המעורבים בתהליך הכנת ספריית AQRNA-seq. התקן הפנימי הוא RNA, בעוד שכל שאר האוליגונוקלאוטידים הם DNA. פריימרים PCR ופריימרים לריצוף מותאם אישית המפורטים הם ספציפיים לפלטפורמות הריצוף. ניתן לעצב פריימרים נוספים ל-PCR עם רצפי אינדקס חדשניים. אנא לחץ כאן כדי להוריד טבלה זו.