חיידקים היו בתרבית משתן שנאסף מנשים בהסכמה כחלק ממחקרים שאושרו על ידי ועדת הבדיקה המוסדית 19MR0011 (UTD) ו- STU 032016-006 (UTSW).
1. תרבית שתן משופרת שונה
הערה: כל צעדי התרבות חייבים להתבצע בתנאים סטריליים. עקר את כל המכשירים, הפתרונות והמדיה. נקו את אזור העבודה עם 70% אתנול, ואז הקימו מבער בונזן ועבדו בזהירות קרוב ללהבה כדי להפחית את סיכויי הזיהום. לחלופין, ארון בטיחות ביולוגית מסוג II עשוי לשמש לשמירה על סביבה סטרילית. ללבוש ציוד מגן אישי מתאים (PPE) כדי למנוע חשיפה חיידקים פתוגניים פוטנציאליים.
- ציפוי שתן מלא גלצל ובידוד מושבה
- הפשר שתן המצויד בגליצריל בטמפרטורת החדר (RT). לאחר הפשרה, מערבולת המדגם עבור 5 s לערבב. בצינורות מיקרוצנטריפוגה סטריליים, להכין 1:3 ו 1:30 דילול של השתן סטרילי 1x פוספט-חוצץ מלוח (PBS) לנפח סופי של 100 μL.
הערה: שתן מלאי גליצל מוכן על ידי ערבוב 500 μL של שתן לא דולל ו 500 μL של 50% גלילסול סטרילי cryovials ולאחסן ב -80 °C (80 °F).
- צלחות אגר חמות מראש ב 37 °C (5 °F) במשך 15 דקות לפני השימוש. אנא עיינו באיור 1 לסוגי מדיה ולתנאי תרבות המתאימים לסוגי חיידקי השתן הנפוצים. מערבבים היטב את השתן המדולל על ידי צנרת לפני הציפוי, צלחת 100 μL של השתן המדולל על צלחת אגר הרצויה ולהפיץ את המדגם באמצעות חרוזי זכוכית סטריליים. צלחת 100 μL של 1x PBS דילול על צלחת נפרדת כמו בקרת צמיחה.
הערה: אם מנסים לתרבות מינים אורופתיוגנייםנפוצים(למשל, Escherichia coli, Klebsiella spp., Enterococcus faecalis וכו '), מומלץ להשתמש אגר כרומוגניים ( טבלת חומרים ) כפי שהוא מאפשר זיהוי קל של מינים חיידקיים אורופתיגניים (איור 1). חומצה קוליסטין Nalidixic (CNA) או MRS אגר שימושיים לבידוד מינים גראם חיובי בררני (למשל, לקטובצילוס spp.) משתן ידוע להכיל אורופטוגנים גרם שלילי, אשר עשויים להערים על המינים בררניים אגרים שאינם סלקטיביים.
- לדגור את הצלחת הפוכה במצב האטמוספרי הרצוי ב 35 °C (35 °C) לתקופה של 24 שעות עבור אורופטוגנים ו 3-5 ימים עבור חיידקים בררניים(איור 1).
- לאחר תקופת הדגירה, להסיר את הצלחות מן האינקובטור. מכל צלחת, בחר את המושבות המציגות צבע ייחודי, מורפולוגיה או דפוסים המוליטיים.
- פס מחדש את מושבת החיידקים באמצעות לולאה סטרילית על אגר המתאים ודגורת הצלחת הפוכה במשך 2-5 ימים באטמוספירה הרצויה כדי להשיג מושבות מבודדות היטב.
הערה: אם משתמשים BAP לתרבות ראשונית, תיקון מושבות על אגר כרומוגניים עשוי לספק מידע שימושי על ההטרוגניות של אוכלוסיית החיידקים במדגם.
- פולחן במרק נוזלי וגליצל-גרביים חיידקים מבודדים
- לאחר שהמושבות המבודדות התואמות את המורפולוגיה של מושבת ההורים מתקבלות, בחרו מושבה אחת וחסן ל-3 מ"ל של מרק נוזלי באמצעות לולאת חיסון סטרילית. עיין באיור 1 למרק המסוגל לתמוך בצמיחת המיקרוביוטה נפוצה בדרכי השתן. לאטום את צלחות אגר עם parafilm ולאחסן אותם ב 4 °C (4 ימים). לדגור על תרבויות נוזליות בתנאים האטמוספריים הרצויים במשך 1-5 ימים עד שהתרבות היא מעוכבת בעליל.
- לאחר הצמיחה הוא ציין, מערבולת התרבות, ולאחר מכן להוסיף 1 מ"ל של תרבות הלילה כדי 500 μL של גליצרל סטרילי 50% בקריוביאלי 2 מ"ל; לאטום ולערבב בעדינות על ידי היפוך. הכן שתי מניות גליצריל עבור כל מושבה (אחד משמש כגיבוי) ולאחסן ב -80 °C (80 °F).
2. זיהוי מינים חיידקיים על ידי רצף גן סנגר 16S rRNA
הערה: זהות מיקרוביאלית יכולה להיות מאושרת לחילופין באמצעות זמן יינון לייזר בסיוע מטריצה של ספקטרומטר מסה טיסה (MALDI-TOF)20.
- תגובת שרשרת מושבה-פולימראז (PCR)
- הכן 25 μL של תגובת PCR בצינורות PCR על ידי הוספת 12.5 μL של 2x טאק פולימראז מאסטר מיקס, 0.5 μL של 10 μM 8F פריימר, 0.5 μL של 10 μM 1492R פריימר(שולחן החומרים),ו 11.5 μL של מים ללא נוקלאז21.
הערה: אם אתם מבצעים PCR עבור דגימות מרובות, בצעו תערובת ראשית של תערובת טאק פולימראז, פריימרים ומים סטריליים ללא גרעין. ואז aliquot 25 μL לתוך כל צינור PCR.
- כדי לבצע מושבה-PCR, החלק מושבה מבודדת היטב מן הפס מחדש באמצעות קיסם סטרילי או קצה פיפטה. resuspend המושבה בתערובת תגובת PCR מוכן בשלב 2.1.1. מערבבים בעדינות. לאסוף את הנוזל בתחתית הצינור על ידי סיבוב מהיר ב 2000 x g.
הערה: ודאו שהדגימה נקייה מבועות אוויר. כלול דוגמה ללא פקד תבנית (NTC) המכילה את תערובת התגובה של PCR בלבד.
- מניחים את צינורות המדגם בתרמוסקלר ולהפעיל את התוכנית הבאה: 95 °C (55 °F) במשך 3 דקות; 40 מחזורים של: 95 °C (55 °F) ל-30 מעלות צלזיוס, 51 °C (50 °F) ל-30 מעלות צלזיוס ו-72 °C (72°F) למשך דקה 30; 72 °C (72 °F) למשך 10 דקות; החזק ב 10 °C (50 °F).
- מיצוי ג'ל וזיהוי מינים
- עם השלמת ריצת ה- PCR, בדוק את מוצר ה- PCR על ג'ל אגרוז 1% שהוכן במאגר 0.5x Tris-Borate-EDTA (TBE). לפני הליהוק של הג'ל, מוסיפים אתידיום ברומיד (EtBr). לאחר מכן, להטיל את הג'ל באמצעות מסרקים עבור בארות להחזיק לפחות 20 נפח מדגם μL.
זהירות: EtBr הוא סוכן intercalating חשוד להיות מסרטן. תמיד ללבוש כפפות PPE בעת טיפול בו ולהיפטר מחומרים המכילים EtBr על פי הנחיות המוסד.
- כאשר הג'ל מוגדר, מניחים את הג'ל במיכל האלקטרופורזה מלא במאגר TBE 0.5x ומסירים את המסרק. לטעון את סולם 1 kb בבאר הראשונה ו 10-20 μL של תגובת PCR לתוך הבא. הפעל ב 100-140 V עד לפתרון. דמיינו את הג'ל תחת אור UV ואשרו את נוכחותה של רצועה מוגדרת בבירור ב- ~ 1.5 kb הנעדרת בבאר NTC.
אזהרה: קרני UV מזיקות לעור ולעיניים, השתמשו בשומר מתאים בעת הדמיית הג'ל ולובש PPE מתאים.
הערה: PCR המושבה עשוי להיות לא מוצלח עבור חיידקים מסוימים; המשך עם PCR מ gDNA מבודד היא אפשרות חלופית22.
- Excise את רצועות ~ 1.5 kb באמצעות סכין גילוח ולהעביר את ג'ל חתכים לתוך צינורות microcentrifuge נקי. המשך עם פרוטוקול מיצוי ג'ל בהתאם להוראותהיצרן (טבלת חומרים). למדוד את הריכוז של DNA מטוהר על ידי ספקטרופוטומטר microvolume.
הערה: ריכוז >10 ננוגרם /μL רצוי, ו A260/280 בין 1.7-2.0 מקובל.
- הכן שתי תגובות רצף סנגר עבור כל מדגם, אחד באמצעות 8F והשני באמצעות פריימר 1492R במים ללא נוקלאז על פי ההנחיות של כל שירות רצף סנגר שנבחר.
- לאחר קבלת נתוני הרצף, העלה את רצפי הדנ"א לאתר האינטרנט של כלי חיפוש היישור המקומי הבסיסי של NCBI (BLAST) (blast.ncbi.nlm.nih.gov/Blast.cgi), בחר Nucleotide BLAST (blastn), בחר את רצפי ה- RNA rRNA / ITS 16S ריבוזומלי (חיידקים וארכיאה) והפעל את תוכנית Megablast. ניתן לזהות את הבידוד על-ידי הפגיעה באיכות הגבוהה ביותר בהפניה ממסד הנתונים.
הערה: מינים חיידקיים מסוימים מפגינים זהות גבוהה ברצפי rRNA 16S שלהם ועשויים להיות מובדלים בשיטה זו בלבד. דגימה תדרוש הומולוגיה DNA ניתוחים ביוכימיים כדי להבחין בביטחון חברים מאותו סוג23.
3. הפקת דנ"א גנומי (gDNA)
הערה: סעיף זה משתמש ריאגנטים ועמודות ספין המסופקים בערכת החילוץ gDNA המוזכרת בטבלת החומרים להפקת תשואה גבוהה של DNA גנומי איכותי ממיני חיידקים מגוונים. להלן שינויים מומלצים והוראות.
- הכן ריאגנטים של ערכה לפי הוראות היצרן.
- הכן תרבויות של 3-10 מ"ל במרק סטרילי מתאים (איור 1) על ידי חנק חיידקים ממושבות מבודדות היטב לתקשורת ודגרה בטמפרטורה ובלחץ האטמוספרי שצוין באיור 1 עד לצמיחה מספקת.
- לאחר הדגירה, למדוד את הצפיפות האופטית ב 600 ננומטר (OD600)של התרבות באמצעות ספקטרופוטומטר24.
- הכן את המדגם לכימות על ידי דילול תרביות לילה ביחס של 1:10. כלול ריק של מדיה תרבותית סטרילית למדידה גם כן. חשב את הצפיפות האופטית על-ידי חיסור הקריאה הריקה מקריאת המדגם והכפלה בגורם הדילול של עשר.
- באמצעות מדידת OD600 ויחס OD600 ל- CFU / mL שנקבע מראש עבור המינים, לחשב כמה מיליליטרים של תרבות נחוצים כדי להשיג 2 x 109 תאים.
- צנטריפוגה נפח התרבות הנדרש במשך 5 דקות ב 5000 x g לכדור. שאפו את supernatant ו resuspened הכדור ב 200 μL חיץ TE קר (מראש צמרמורת על קרח בתחילת ההליך).
- צנטריפוגה הדגימה במשך 2 דקות ב 5000 x g. הסר את supernatant, ולאחר מכן resuspend הכדור ב 180 μL של חוצץ תמיסת אנזימטית (ELB) ולהוסיף 20 μL של RNase מבושל מראש (10 מ"ג / מ"ל). לקבלת תמוגה יעילה של חיידקים גראם-חיוביים, יש להוסיף 18 מיקרו-אל של מוטנוליזין (25 קילוואט/מ"ל). מערבולת היטב, ולאחר מכן לדגור על הדגימות ב 37 °C (50 °F) על סיבוב במשך 2 שעות.
הערה: מומלץ להשתמש ELB המתואר בפרוטוקול היצרן עבור חיידקים גרם חיובי ו גרם שלילי.
- המשך בהתאם להוראות היצרן.
הערה: חזור על שלבי ההעמקה עבור פעם אחת או פעמיים נוספות כדי להשיג תפוקת gDNA נוספת, אם תרצה.
- להעריך את האיכות של gDNA שחולץ כפי שהורה בסעיף 4 ולאחסן gDNA ב 4 °C (70 °F) אם זה ישמש בתוך שבוע אחד. לחלופין, שמור gDNA ב-20 °C לאחסון לטווח ארוך.
4. הערכת איכות ה-GDNA שחולץ
- כדי להעריך את האיכות על ידי אלקטרופורזה ג'ל, להכין 1% ג'ל אגרוז כמתואר בסערה 2.2. הכן את המדגם בצינור נקי: לערבב 1-2 μL של gDNA שחולץ ו 3 μL של צבע טעינה 2x על parafilm. הפעל את הג'ל פעם טעון, ולאחר מכן לדמיין אותו תחת אור UV.
הערה: חילוץ gDNA מוצלח יתברר על ידי להקה נפרדת בחלק העליון של הג'ל והכפשה מינימלית (איור 2A). מריחה מעידה על גימה. אם לא ברורה רצועת gDNA ו/או מריחה היא משמעותית, חזור על חילוץ gDNA. שקול להפחית את זמני הדגירה ב RNase A ו Proteinase K. אם נצפות שתי רצועות סביב 1.5-3 kb, הדבר מצביע על זיהום RNA(איור 2B). הכן RNase A טריים וחזור על החילוץ.
- כדי להעריך את האיכות על ידי ספקטרופוטומטר מיקרו-וולטום, מדוד את ריכוז ה-gDNA ואת יחס הספיגה A260/280 על ידי ספקטרופוטומטר מיקרו-וולטום. ריכוזים >50 ננוגרם /μL ו- A260/280 בין 1.7-2.0 מקובלים.
הערה: תפוקת gDNA נמוכה עשויה לנבוע מקלט נמוך, קלט גבוה, זיהום של גרקלאזים, תמוגה לא מספקת. יחסי ספיגה מעל הטווח מצביעים על זיהום RNA. מיצוי חוזר אם איכות gDNA ירודה.
- כדי להעריך את האיכות על ידי פלואורומטר, בצע את הוראות היצרן לכמת את ריכוז gDNA באמצעות ערכת בדיקת רגישות גבוהה ומכשיר פלואורומטר(שולחן החומרים). ריכוז >50 ננוגרם / μL רצוי.
5. הדור הבא של הדור הבא בהכנה קצרת קריאה והכנת ספריה
הערה: ניתן לבצע רצף לקריאה קצרה בכלים שונים באורכים ובהכוונה ברורים לקריאה. רצף 150 bp (300 מחזור) מומלץ עבור WGS חיידקי. הן הכנת הספרייה והן הרצף עשויים להיות במיקור חוץ למתקני ליבה או למעבדות מסחריות.
- הכן ספריית רצף בהתאם להוראות היצרן (טבלת חומרים). עקוב אחר ריכוז ספריית הטעינה הסופית המומלץ של היצרן; עם זאת, שינוי מומלץ הוא לטעון את הספרייה המאגדת ב- 1.8 pM ליצירת קריאה אופטימלית במכשירי NextSeq.
- למרות האופציה, השתמש ביואנלייזר(טבלת חומרים)כדי להעריך את התפלגות שברי הספרייה במאגר ולוודא שגודל הקטע הוא 600 bp בממוצע.
6. הכנת ספריית רצף מיניון ננופור
- הכן את ספריית הרצף בהתאם לפרוטוקול היצרן (טבלת חומרים). שימוש בשתי ערכות הרחבת ברקוד מאפשר ריבוי של עד 24 דגימות בתא זרימה יחיד. מומלץ לבצע הכנת ספריה בשני חלקים, 12 דגימות בכל פעם בעת מולטיפלקסינג 24 דגימות. ניתן לאגד את כל 24 הדגימות כמתואר להלן.
הערה: דגימות עשויות להיות מאוחסנות ב 4 °C (5 °F) לילה עם סיום קשירת ברקוד מקורי - זה מספק נקודת עצירה בפרוטוקול, במידת הצורך. בסוף סעיף קשירת הברקוד המקורי של פרוטוקול הכנת הספרייה, מומלץ לאסוף כמויות שווימולריות של כל דגימה עד למסת ה- DNA המרבית (ng) האפשרית.
- כדי לעשות זאת, לכמת את כל הדגימות הבאות קשירת ברקוד באמצעות פלואורומטר (טבלת חומרים) לפי הוראות היצרן. להעריך את נפח המדגם עם ריכוז dsDNA הנמוך ביותר ולאחר מכן לחשב את הסכום הכולל dsDNA נמצא במדגם זה. השתמש במספר זה כדי לקבוע את הכמויות ההזוותיות של כל הדגימות האחרות שיאחדו יחד.
הערה: מכיוון שהחישוב השוומור ימקסם את כמות ה- dsDNA המאוחס ובכך יניב מאגר בנפח גבוה (>65 μL), יש צורך בניקוי כדי לרכז את הבריכה.
- ניקוי וריכוז בבריכת dsDNA
- הוסף נפח של 2.5x של חרוזים פרמגנטיים (שולחן החומרים) למאגר הדנ"א, ולאחר מכן החלק בעדינות את הצינור כדי לערבב את התוכן. מניחים את הצינור ברוטציה במשך 5 דקות ב RT. ספין את המדגם ב 2000 x g וכדור על מגנט.
- מוסיפים 250 μL טרי מוכן 70% אתנול (במים ללא גרעין), דואג לא להפריע לכדור. שאפו את האתנול וחזרו על שטיפת האתנול פעם אחת.
- לאחר השאיפה השנייה, לסובב את המדגם ב 2000 x g ולהניח אותו בחזרה על המגנט. פיפטה את כל אתנול שיורית ולאפשר לדגימה להתייבש במשך כ 30 s.
- הסר את הצינור מן המגנט ו resuspend הכדור ב 60-70 μL של מים ללא גרקלאז. דגירה ב RT במשך 2 דקות. גלם את המדגם על המגנט עד elute הוא ברור, ולאחר מכן להסיר את elute ולהעביר לתוך צינור microcentrifuge 1.5 מ"ל נקי.
- לכמת את הבריכה המרוכזת באמצעות פלואורומטר, ולאחר מכן להכין aliquot להמשיך לשלב קשירת המתאם: להכין 700 ננוגרם של המדגם בנפח הסופי 65 μL. שמור על שאר הבריכה ב 4 °C (7 °F) עבור ריצה שנייה כדי להסתיים לאחר הריצה הראשונה הסתיימה.
- המשך עם קשירת מתאם בהתאם להנחיית היצרן וטען את הדגימה בתא הזרימה. התחל את ריצת הרצף.
הערה: שאיפה אוויר ~ 200 μL של מאגר אחסון מיציאת התחלה של תא הזרימה לפני טעינת המדגם. הדבר קריטי עבור תימוי תאי הזרימה המוצלחים וטעינת הדגימה. השתמש בפיפטה p1000 וטיפים בעת ציור והפקדת פתרונות דרך היציאה הראשונית של תא הזרימה.
- רצף את הספריה בהתאם להוראות היצרן.
- פתח את תוכנת ההפעלה לרצף ולחץ על התחל. הזן שם עבור הניסוי, המינוח המומלץ כולל את תאריך ההפעלה ואת שם המשתמש. לחץ על המשך לבחירת ערכות, בחר בערכת ההכנה המתאימה לספריה ובחבילות הרחבה של ברקוד המשמשות ולאחר מכן לחץ על המשך כדי להפעיל אפשרויות.
- התאם את אורך הריצה ל- 48 שעות אם הוא מתכנן להכין ספריה מספקת להפעלה שנייה (אחרת השאר כברירת מחדל 72 שעות). לחץ על המשך להתקשרות בסיסית.
- בדוק את אפשרות basecalling Config: Basecalling מהיר וודא ש- Barcoding מוגדר כזמין כך שקבצי FASTQ פלט יקצצו את רצפי הברקוד ויהפכו לספריות נפרדות המבוססות על ברקוד. לחץ על המשך לפלט.
- בחר היכן לשמור נתוני ריצוף פלט. צפה לכ-30-50 ג'יגה-בתים של נתונים אם רק שמירת פלט FASTQ ו->500 Gb של נתונים אם גם שומרת את פלט FAST5. בטל את הסימון באפשרות הסינון Qscore: 7 | אורך קריאה: לא מסונן אם הוא מתכנן להמשיך בסינון המתואר בסעיף 7.2, השאר אחרת מסומן ולהתאים את אורך הקריאה ל- 200.
- לחץ על המשך להפעיל את תוכנית ההתקנה וסקור את כל ההגדרות. אם ההגדרות נכונות, לחץ על התחל, אחרת לחץ על אחורה ובצע את כל ההתאמות הדרושות.
- במידת הרצון, ניתן לשטוף את תא הזרימה בהתאם להוראות היצרן ולהיטען מחדש עם המאגר הנותר. חזור על השלבים ב- 6.2 עבור המאגר הנותר לאחר השלמת ההפעלה הראשונה ותא הזרימה נשטף.
הערה: בעת הגדרת הריצה השנייה, התאם את מתח ההטיה ל- -250 mV לפי המלצות היצרן לתאי זרימה ששימשו בעבר בריצות מעל 48 שעות.
7. הערכה והכנת קריאות
הערה: מבנה ספריה מומלץ מתואר באיור 4. צור את הספריות הנמצאות בשולחן העבודה, כלומר, Long_Reads, Short_Reads ו- Trimmed_Reads, לפני שתמשיך עם שלבי החישוב שלהלן.
- קריאות קצרות (איור 3)
הערה: קריאות קצרות נוצרות בתבנית FASTQ. הקבצים מכילים 4000 קריאות מרביות לכל FASTQ. לעתים קרובות הם מכווצים (ארכיון .gz) ומאורגנים בקבצים מרובים. בהתאם לפלטפורמה, ברקודים נחתכים בדרך כלל. תוכניות מסוימות מקבלות קבצים בתבנית מכווצת, אחרות עשויות לדרוש את החילוץ שלהם לפני הייבוא. הקריאות חייבות לעבור שלבי בקרת איכות (QC) כדי להבטיח דיוק נתונים במהלך הרכבת הגנום. אם שולחן העבודה של CLC Genomics אינו זמין, תוכניות חלופיות עשויות לשמש לקיצוץ וקריאה קצרה של QC כגון Trimmomatic25 או חתוך גוון בשפע (https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/) לגיזום ו- FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) להערכת איכות הקריאה. כיסוי קריאה קצר ממוצע, המוערך על ידי הכפלת מספר הקריאות לפי אורך הקריאה הממוצע וחלוקה לפי גודל הגנום, מומלץ להיות >100x.- פתח את תוכנת שולחן העבודהשלGenomics ( טבלת חומרים ) וייבוא את כל קבצי FASTQ קצרי הקצה. קבצים משויבים ייווצרו באופן אוטומטי.
- צור תיקיה חדשה תחת CLC_Data על-ידי לחיצה על חדש בסרגל הכלים העליון ובחירה בתיקיה... כדי לאחסן את הקבצים. תן שם לתיקיה כרצונך, מוסכמה מומלצת משתמשת במזהה לדוגמה. שמור את כל הפלט מהשלבים הבאים לתיקיה זו.
- בסרגל הכלים העליון, לחץ על לחצן ייבוא ובחר אילומינה... נווט אל כל הקבצים שנקראו קצרים ובחרו התואמים לדוגמה. ודא שהאפשרות קריאות משויפות נבחרה ובטל את הסימון באפשרות הסר קריאות שנכשלו. לחץ על הבא, בחר שמורולחץ שוב על הבא. בחר לשמור את הקבצים המיובאים בתיקיה החדשה שנוצרה בשלב הקודם ולחץ על סיום.
- צור רשימת רצף של כל הקבצים המשויים עבור הבידוד; פעולה זו תצרר נתוני קריאה לקובץ יחיד לפשטות הניתוח.
- בסרגל הכלים העליון, לחץ על לחצן חדש ובחר רשימת רצף... ברשימת הספריות מימין, בחר את הקבצים שיש לערוץ והשתמש בחצים כדי להעביר אותם לרשימת הקבצים שנבחרו מימין. לחץ על הבא, בחר שמורולחץ שוב על הבא. בחר לשמור את רשימת הרצף ולחץ על סיום.
- לאחר יצירת רשימת הרצף, שנה את שמה באופן מיידי למזהה לדוגמה.
- הפעל את הכלי QC לקריאות רצף ברשימת הרצף: הליך זה יעריך את פרמטרי האיכות הכוללים של הקריאות שנוצרו על-ידי NGS לקריאה קצרה.
- חפש את הכלי QC לרצף קריאות בתפריט ארגז הכלים (חלון מימין למטה). לחץ פעמיים על הכלי ולאחר מכן בחר את רשימת הרצף שיש לנתח ולחץ על הבא.
- ודא שכל אפשרויות הפלט מסומנות ובחר שמור תחת טיפול בתוצאות. לחץ על הבא וציין כדי לשמור את קבצי הפלט ולאחר מכן לחץ על סיום.
- הפעל את הכלי Trim Reads ברשימת הרצף: החיתוך ייעשה בהתבסס על איכות, אורך ואי בהירות. תהליך זה מניח שהברקודים המשמשים לרצף קוצצו לפני שלב זה.
- חפשו את הכלי חיתוך קריאות בארגז הכלים (חלון מימין למטה). לחץ פעמייםעל חיתוך קריאות ולאחר מכן בחר את רשימת הרצף שיש לנתח ולחץ על הבא.
- חיתוך איכות: הגדר את מגבלת ציון האיכות ל- 0.01 והותיר נוקלאוטידים מעורפלים ב- 2. לחץ על הבא.
הערה: ניתן להתאים פרמטרים לפי שיקול דעת המשתמש; אלה ההגדרות המומלצות.
- בטל את הסימון של חיתוך מתאם קריאה אוטומטי (בצע זאת רק אם מתאמים נחתכו מהקריאות לפני הייבוא ל- CLC). לחץ על הבא וסמן את מחק קריאות מתחת לאורך, השתמש בברירת המחדל 15.
- לחץ על הבא, סמן יצירת דוחולאחר מכן בחר שמור. לחץ על הבא וציין היכן לשמור את קבצי הפלט. לחץ על סיום.
- ייצוא רשימת הרצף החתוך: ההרכבה והניתוח ההיברידיים הבאים יושלמו מחוץ ל- CLC ודורשים ייצוא של קבצים קצוצים לקריאה קצרה.
- מתוך ניווט הספריה מימין למעלה, בחר את הקובץ החתוך שנוצר בשלב 7.1.4 ולאחר מכן לחץ על ייצוא בסרגל הכלים העליון. בחר Fastq עבור סוג קובץ הייצוא ולחץ על הבא. בדוק יצא רשימת רצף מזווג לשני קבצים. לאחר מכן, לחץ על הבא ובחר את הספריה Trimmed_Reads לייצוא הקבצים. לחץ על סיום. ודא שהקבצים החתוך לקריאה קצרה יוצאו בהצלחה כשני קבצים (R1 ו- R2) עם הסיומת .fastq.
הערה: יש לייצא את רשימת הרצף החתוך לשני קבצים, המוגדרים בדרך כלל על-ידי CLC כ- R1 ו- R2. הדבר קריטי מכיוון שהרכבה היברידית במורד הזרם דורשת להגדיר קלט נתונים לקריאה קצרה ככזה.
- שנה את שמות הקבצים המיוצאים, נא להימנע משימוש ברווחים ובתווים מיוחדים בשמות קבצים. לפשטות trimmed_short_file תבנית מומלצת. אר-1.פאסטק.
- קריאות ארוכות (MinION) (איור 3)
הערה: הצינור הבא להכנת קריאות רצף ארוכות (MinION) להרכבה היברידית משתמש בתוכניות NanoFilt ו- Nanostat26 המבוצעות על-ידי שורת הפקודה. התקן את הכלים לפני שתמשיך והכר את היסודות של UNIX כדי לבצע פקודות אלה. מומלץ למסופי ברירת מחדל ולמעטפת Bash. מדריך שיעורים עבור פקודות מסוף נפוצות ושימוש נמצא ב- Software Carpentry27. ההוראות שלהלן מניחות כי הקבצים שנוצרו ייקראו עם המינוח של הברקוד (NB01, NB02 וכו ') ונשמרו בספריה Long_Reads. לחלופין, ניתן לבצע סינון קריאה באמצעות MinKNOW בעת הגדרת הפעלת הרצף. כיסוי קריאה ארוך ממוצע מומלץ להיות >100x. אורך הקריאה הממוצע המומלץ הוא >2000 bp; לכן, מספר הקריאות הארוכות הדרושות נמוך ממספר הקריאות הקצרות.
- צור ספריות חדשות עבור כל ברקוד המשמש בהפעלה (ברקוד01, ברקוד02 וכו') בתוך הספריה Long_Reads(איור 4). העתק את כל קבצי ה- .fastq התואמים לכל ברקוד לתיקיה המתאימה. שלב את כל קבצי .fastq עבור כל ברקוד מכל ריצה.
- פתח את Terminal ונווט אל ספריות הברקוד בספריה Long_Reads באמצעות הפקודה cd: cd Desktop/Long_Reads/ברקוד01
- לכנס את כל קבצי .fastq לכל ברקוד לקובץ .fastq יחיד על ידי ביצוע הפקודה הבאה: cat *.fastq > NB01.fastq
הערה: פקודה זו משלבת את כל הקריאות מכל אחד מקבצי FASTQ ל- FASTQ אחד גדול ובודד בשם NB01.fastq.
- השתמש NanoStat כדי להעריך את איכות הקריאה של המדגם על ידי ביצוע הפקודה הבאה: NanoStat --fastq NB01.fastq
- הקלט את התוצאות על-ידי העתקת הפלט לקובץ טקסט או Word לעיון עתידי.
- השתמש NanoFilt כדי לסנן MinION קורא ביטול קריאות עם Q < 7 ואורך < 200 על ידי ביצוע הפקודה: NanoFilt -q 7 -l 200 bp NB01.fastq | gzip > NB01 _trimmed.fastq.gz
- הפעל NanoStat על הקובץ החתוך שנוצר בשלב 7.2.6 על ידי ביצוע הפקודה: NanoStat -fastq NB01 _trimmed.fastq.gz
- הקלט את התוצאות על-ידי העתקת הפלט לקובץ טקסט או Word והשווה לתוצאות שלב 7.2.4 כדי להבטיח שהסינון הצליח (טבלה 1).
- חזור על שלבים 7.2.2 עד 7.2.8 עבור כל ברקוד המשמש בריצת הרצף.
הערה: הקובץ NB01_trimmed.fastq.gz שנוצר בשלב 7.2.6 ישמש להרכבה היברידית.
8. יצירת הרכבה גנומית היברידית
הערה: צינור ההרכבה הבא משתמש ב- Unicycler19,28,29,30 כדי לשלב קריאות קצרות וארוכות שהוכנו בסעיפים 7.1 ו- 7.2 ( איור3). התקן את Unicycler ואת יחסי התלות שלו ובצע את הפקודות שלהלן. מניחים שקבצים הנקראים זמן קצר שיוצאו בשלב 7.1.5 נקראים trimmed_short_file. אר-1.פאסטק trimmed_short_file. R2.fastq לפשטות.
- ארגן את הקבצים שנקראו זמן קצר ואת הקבצים שנקראו זמן רב בספריה אחת בשם Trimmed_Reads. הספריה חייבת להכיל את הדברים הבאים:
- קובץ .fastq.gz עבור קריאות ארוכות חתוכות (נוצר בשלב 7.2.6).
- שני קבצי .fastq (R1 ו- R2) לקריאות קצרות חתוכות (שנוצרו בשלב 7.1.5).
- נווט אל הספריה Trimmed_Reads המאחסנת את קבצי הקריאה באמצעות הפקודה cd בטרמינל: cd Desktop/Trimmed_Reads
- פעם בספריה הנכונה, תסגור את שני קבצי הקריאה הקצרים כך שהם גם בתבנית .fastq.gz על-ידי ביצוע הפקודה הבאה: gzip trimmed_short_file. R1.fastq
- חזור על שלב 8.2 הן עבור R1 והן עבור R2. ודא שכל קבצי הקריאה נמצאים כעת בתבנית .fastq.gz וודא שכל הקבצים תואמים לאותו בידוד.
- התחל את ההרכבה ההיברידית באמצעות חד-ציקלילר על-ידי הפעלת הפקודה הבאה:
חד-קליט -1 trimmed_short_file. R1.fastq.gz -2 trimmed_short_file. R2.fastq.gz -l NB01 _trimmed.fastq.gz -o unicycler_output_directory
הערה: -o מציין את הספריה שבה תישמר פלט יוניציקלר, Unicycler ייצור ספריה זו לאחר ביצוע הפקודה; אל תיצור את הספריה מראש. זמן הריצה משתנה בהתאם לעוצמה החישובית של המחשב המשמש, כמו גם גודל הגנום ומספר הקריאות. פעולה זו עשויה להימשך בין השעות 4 שעות ל-1 או 2 ימים. פרוטוקול זה בוצע על מכונת CentOS Linux 7 עם זיכרון RAM של 250 ג'יגה בייט, מעבד Intel Xeon (R) עם ליבות מעשיות של 2.5 ג'יגה-הרץ ו-48 ליבות וירטואליות. לחלופין, מחשבים אישיים עם זיכרון RAM של 16 ג'יגה-בתים ומעבדי 6 ליבות ב- 2.6 GHz יכולים לחשב הרכבות אלה בזמן עיבוד ארוך יותר.
- לאחר השלמת ההפעלה, סקור את קובץ .log unicycler כדי להבטיח שאין שגיאות - הקלט את המספר, הגודל והמצב (שלם, לא שלם) של ההסתבפויות שנוצרו.
- אם מזוהים קונטיגים לא שלמים (מסומנים כלא שלמים ביומן יוניציקלר), הפעל מחדש את Unicycler במצב מודגש על-ידי הוספת הדגל הבא לפקודה בשלב 8.4: --mode מודגש.
הערה: מצב מודגש יוריד את סף האיכות המקובל עבור גשרי קריאה ארוכה במהלך ההרכבה; זה עשוי להניב הרכבה מלאה, אבל איכות ההרכבה עשויה להיות מופחתת. מומלץ להשתמש במצב מודגש רק בעת הצורך וכראיה ראשונית להצטרפות קונטיג שיאושרו מאוחר יותר על ידי PCR.
9. הערכת איכות ההרכבה
הערה: הפרוטוקול הבא משתמש בתחבושת31 וב- QUAST32, שתי תוכניות שיש להגדיר לפני השימוש (איור 2 ואיור 4). התחבושת אינה דורשת התקנה לאחר ההורדה ו- QUAST דורש היכרות עם שימוש בסיסי בשורת הפקודה. מומלץ גם להעריך את שלמות הגנום באמצעות אורתולוגיות אוניברסליות בעלות העתק יחיד (BUSCO)33.
- תחבושת: לחץ על קובץ. לאחר מכן, בחר טען גרף ובחר את הקובץ assembly.gfa שנשמר unicycler_output_directory שנוצר על-ידי Unicycler בשלב 8.4. לאחר טעינה, לחץ על לחצן צייר גרף בסרגל הכלים השמאלי ובדוק כיצד הקונטיגים (הנקראים צמתים) מחוברים ומאורגנים כדי להעריך אם ההרכבה הושלמה (איור 5).
הערה: הרכבות שלמות מיוצגות על-ידי קונטיגסים מעגליים בודדים המקושרים בשני קצותיה (איור 5A,B). להרכבות שלא הושלמו יש קונטיגים מרובים המקושרים זה לזה או שהם ליניאריים(איור 5C). קונטיגס ליניארי קטן לא יכול להיות שלם כפי שהם עשויים להצביע על אלמנטים חוץ-כרומוזומאליים ליניאריים. כיסוי, המכונה גם עומק, יצוין בתחבושת ומייצג את השפע היחסי של הרצף לכרומוזום, מנורמל ב Unicycler ל 1x.
- QUAST
- בתוך המסוף, נווט לתיקיה המאחסנת את פלט Unicycler באמצעות הפקודה cd: cd Desktop/Trimmed_Reads/unicycler_output_directory
הערה: רווחים אינם מותרים בנתיב שבו ממוקמת ההרכבה, כלומר, אין ספריות המובילות לפלט יוניציקלר שיכולות לכלול רווחים על שמם. לחלופין, העתק את הקובץ assembly.fasta לשולחן העבודה לקבלת גישה נוחה.
- הפעל את QUAST על-ידי ביצוע הפקודה הבאה: assembly quast.fasta -o quast_output_directory
- סקור את הדוחות שנוצרו על-ידי QUAST בספריה quast_output_directory.
10. ביאור גנום
הערה: צינור הביאורים שלהלן משתמש ב- Prokka34, כלי שורת פקודה שיש להתקין לפני השימוש. לחלופין, השתמש Prokka דרך אוטומטי GUI K-Base(טבלת חומרים)או גנומים ביאורים באמצעות שרת האינטרנט RAST35. אם הפקדת גנומים לתוך NCBI, הם יירו באופן אוטומטי באמצעות צינור ביאור הגנום הפרוקריוטי (PGAP)36.
- נווט בתוך המסוף לתיקיה המאחסנת את פלט יוניציקלר באמצעות הפקודה cd (ראה שלב 9.2.1). לאחר מכן, הפעל את Prokka על ידי ביצוע הפקודה הבאה: prokka -קידומת sample_ID -להקדים את prokka_output_directory ההרכבה.fasta
הערה: --קידומת תציין את כל קבצי הפלט בהתבסס על sample_ID שצוין. --outdir ייצור ספריית פלט בשם שצוין שבו יישמרו כל קבצי הפלט של Prokka; אל תיצור ספריית פלט עבור Prokka מראש.
- סקור את הביאורים על-ידי פתיחת טבלת ה- .tsv ו/או על-ידי העלאת קובץ ה- .gff שנוצר לתוכנת ניתוח רצף כדי להציג באופן חזותי ולנתח את הביאורים (איור 6).
- סוגים מסוימים של ביאורים יכולים להיווצר בהתאם לגורמים גנטיים של עניין. מומלץ להתחיל עם הכלים הידידותיים למשתמש על המרכז לאפידמיולוגיה גנומית (www.genomicepidemiology.org/) שרת אינטרנט לניתוח ראשוני37,38,39,40,41. כלים נוספים לזיהוי מערכות CRISPR-cas ופרואפאז' זמינים (איור 3)42,43.
11. שיטות מומלצות לדמוקרטיזציה של נתונים
- במידת האפשר, הפקד את כל נתוני הקריאה הגולמיים, כמו גם גנומים שהורכבו במאגר ציבורי כגון ארכיון הקריאה של רצף NCBI (SRA) וג'נבנק. הגנומים מובאים באופן אוטומטי באמצעות צינור PGAP במהלך תהליך התצהיר של NCBI.