פרוטוקול זה מציג שלבים של ניתוח metabarcoding, מיקוד 16S rRNA ו 18S rRNA גנים, לניטור פריחות אצות מזיקות ואת המיקרוביום הקשורים שלהם בדגימות מי ים. זהו כלי רב עוצמה המבוסס על מולקולות אך דורש מספר הליכים, אשר מוסברים חזותית כאן צעד אחר צעד.
ניטור פריחת אצות מזיקות (HABs) יושם ברחבי העולם, וצ’ילה, מדינה המפורסמת בדיג ובחקלאות הימית שלה, השתמשה באופן אינטנסיבי במיקרוסקופיה ובניתוח רעלים במשך עשרות שנים למטרה זו. שיטות ביולוגיות מולקולריות, כגון ריצוף דנ”א בעל תפוקה גבוהה וגישות מבוססות הרכבה חיידקית, רק מתחילות להיות מוצגות בניטור HAB הצ’יליאני, וההליכים טרם הוקרכו. כאן, 16S rRNA 18S rRNA ניתוחי metabarcoding לניטור HABs צ’יליאני מוצגים צעד אחד. על פי השערה עדכנית, אסוציאציה הדדית אצות-בקטריאלית משחקת מערכת יחסים סינרגטית או עוינת קריטית המהווה חניכה, תחזוקה ורגרסיה. לפיכך, ניטור HAB מנקודת מבט אצות-חיידקים עשוי לספק הבנה רחבה יותר של מנגנוני HAB ואת הבסיס לאזהרה מוקדמת. ניתוח Metabarcoding הוא אחד הכלים מבוססי מולקולות המתאימים ביותר למטרה זו כי זה יכול לזהות מידע טקסונומי אצות-חיידקי מסיבי במדגם. ההליכים החזותיים של דגימה לניתוח metabarcoding להלן מספקים הוראות ספציפיות, במטרה להפחית שגיאות ואיסוף נתונים אמינים.
מינים רבים של פיטופלנקטון ימי ידועים כמי שמייצרים רעלנים אנדוגניים, וכאשר מינים אלה מצטברים במספרים מספיקים, הם מזיקים לסביבה הימית. פריחות אצות מזיקות כאלה (HABs) נצפות היום בחופי רוב יבשות העולם1. פיטופלנקטון רעיל מצטבר לראשונה ברקמות דו-שנתיות, מה שמוביל למחלות ומוות ברמות טרופיות גבוהות יותר של אורגניזמים, כולל בני אדם, על העיכול. לאחר מכן, אירועים אלה מביאים השלכות חמורות על הכלכלה המקומית, חברתית-כלכלה, ובריאות הציבור2. הנזק לכלכלה העולמית שנגרם על ידי HABs מוערך במיליוני עד מיליארדי דולרים מדי שנה3. צ’ילה היא אחת מהמדינות הרבות הסובלות מ- HABs תכופים.
צ’ילה היא מדינה ארוכה, המשתרעת על פני 4,300 ק”מ מצפון ומדרום. האדמה המערבית המוארכת פונה לאוקיינוס השקט, מה שמגדיל באופן טבעי את הסיכויים שצ’ילה תחווה HABs. במיוחד בדרום צ’ילה יש הרבה חקלאות ימית סלמון מפורסמת בעולם, ובכל פעם HAB מתרחש באזור, אצות-רעלים לגרום סלמון חקלאי מסיבי לחלות למות4,5,6. בצ’ילה, השנה שבה ה-HAB פגע בכלכלה בצורה הקשה ביותר הייתה 2016, עם הפסד שנתי מוערך של 800 מיליון דולר7,8. מינים אצות רעילים סיבתיים משתנים לפי שנה ומיקום. במקרה של 2016, קומפלקס של אלכסנדריום קטנלה ופסאודוצ’אטלונלה וורוקלוזה גרם ל-HAB נרחב ברוב חופידרום צ’ילה 7,8. HAB האחרון בצ’ילה התרחש עם מינים אצות סיבתיים של Heterosigma Akashiwo במרץ 2021 בקמנג’קה, הממוקם בדרום צ’ילה, שם באזור יש חוות סלמון גדולה.
צ’ילה מנהלת ניטור חופים כבר שנים רבות בשתי שיטות עיקריות; התבוננות במי ים עם מיקרוסקופים כדי לזהות מינים רעילים של אצות באופן קבוע ומדידת רמות הרעלן ברכיכות על ידי בדיקות ביוכימיות9. גילוי מוקדם של אצות רעילות ורמות רעלן ברכיכות אינם מונעים HABs; עם זאת, ניתוחים אלה יכולים להפעיל אמצעי נגד מיידיים ולהפחית את הנזקים לקהילות מקומיות. כדי לחזק עוד יותר את האפקטיביות של אסטרטגיה זו, שיטה אנליטית מבוססת מולקולרית נוספה לאחרונה לתוכנית ניטור HAB הצ’יליאנית הקונבנציונלית שלנו כדי לזהות אצות וקהילות חיידקים קשורות בדגימות מי ים. באופן ספציפי, נבחרה שיטת רצף מקבילי מסיבית באמצעות מטא-ברקודינג המתמקדת בגנים 16S rRNA ו- 18S rRNA. למרות שטכניקה זו דורשת הליכים מסובכים ומכונות יקרות וריאגנטים, זוהי טכנולוגיה מתקדמת שיכולה לזהות אלפי אצות וחיידקים genera / מינים נוכחים בדגימת מי ים בבת אחת.
הסיבות של HAB הם העריך להיות שונים, כגון טמפרטורה ועונה, אבל זה בלתי אפשרי להכליל אותם. הסיבה לכך היא שמינים HAB ותדרים תלויים באזור ובתנאים spatiotemporal, מעורבים תופעות טבע כגון ייחודיות גיאוגרפית, ערבוב חומרים מזינים upwelling, ו נגר אלמנטים מהיבשה עקב לטעות2,10,11,12. בנוסף, גורמים מלאכותיים כגון eutrophication להשפיע על HABs המקומי12,13. בשל הרב-גורמי המורכבים, לא קל לבצע חיזוי מדויק של HAB. בשנים האחרונות, יש השקפה כי אוכלוסיות חיידקים ספציפיות עשויות להיות קשורות להתפתחות HABs כאחד הגורמים, ומחקר לתמיכה בהשערה זו ניכר יותר ויותר14,15,16,17,18. טכניקות ביולוגיה מולקולרית משמשות בדרך כלל לחקר הרכבה חיידקית; עם זאת, שיטה מתוקננת כזו טרם הוקמה ניטור HAB צ’יליאני9. כדי לחקור את הקשר האצות-חיידקי עם HABs, זה הכרחי לבצע במקביל ניתוח metabarcoding עבור תוכניות ניטור החוף הנוכחי צ’ילה. לפיכך, פרוטוקול זה מציג חזותית את תוכנית ניטור HAB הצ’יליאנית שלנו, תוך התמקדות בהליך צעד לניתוח מינים של אצות וחיידקים בדגימות מי ים באמצעות ניתוח metabarcoding.
הפרוטוקול המלא המתאר את תוכנית הניטור של HAB הצ’יליאנית שלנו זמין בירימיזו ואח‘9. הוא כולל את הנהלים של דגימת מי ים, זיהוי מינים אצות מיקרוסקופיים, זיהוי גנים אצות-חיידקיים, ניתוח פיגמנטים, איסוף נתונים מטאורולוגיים, ובוחות תכונות מים פיזיות וכימיות. פרוטוקול הצעד של 16S rRNA וניתוח מטא-בר-רנ”א 18S לזיהוי מינים אצות וחיידקיים זמין כהדפסה מוקדמת19. פרוטוקול זה מדגים במיוחד metabarcoding צעדי ניתוח שכן הוא החלק המורכב ביותר ואת גולת הכותרת של תוכנית ניטור HAB שלנו. פרוטוקול זה כולל גם הקדמה של התוכנית וגילוי מיקרוסקופיה של מינים אצות. בעת ניתוח מינים אצות על ידי metabarcoding, זה חיוני לבצע בו זמנית מיקרוסקופיה כדי לאמת את התוצאות משתי השיטות. פרוטוקול זה אינו כולל כיצד להשתמש בתוכנה להקצאת טקסונומיה, אך המלצת מסד הנתונים נאמרת בקצרה בסוף הסעיף הבא.
פרוטוקול זה ביצע בהצלחה ניתוח מטא-בר-קודינג של גנים 18S rRNA ו-16S rRNA לזיהוי מינים אצות וחיידקיים בדגימות מי ים לניטור HABs צ’יליאני. האצות הדומיננטיות וכמה אצות קלות שזוהו על ידי פרוטוקול זה היו עקביות עם אלה שהושגו על ידי מיקרוסקופיה, המאשר את אמינות הפרוטוקול. גולת הכותרת של הפרוטוקול היא כי הניתוח זיהה אלכסנדריום spp. ו Pseudochattonella spp., שני HAB הבעייתי ביותר גרימת מינים בדרום צ’ילה, ממי הים של מטרי אפילו בשפע נמוך. במיוחד, פסאודוצ’אטלונלה סאפ. הם מינים מאתגרים לזהות תחת מיקרוסקופ כי התאים הם קטנים בקלות נקרע תחת אור או שימוש של מתקנים27,28. מטא-ברקודינג יכול לספק למינים אצות מידע שבקושי קיים בדגימות מי ים שמיקרוסקופיה אינה יכולה לזהות. הפרוטוקול גם זיהה מעל 30 מינים חיידקיים. למרות כיצד חיידקים אלה קשורים לצמיחת האצות אינו ידוע בשלב זה, איסוף נתונים מסיבי של מינים אצות-חיידקים בניתוח סדרת הזמן עשוי לספק תגליות חשובות כאלה. לכן, הוספת ניתוח metabarcoding לתוכניות ניטור HAB צ’יליאני קונבנציונאלי יהיה ללא ספק לחזק את יעילות ניטור HAB הנוכחי. למעשה, פרוטוקול חזותי זה לניתוח metabarcoding יכול להיות מועיל לא רק עבור ניטור אצות-חיידקים של מים צ’יליאניים, אלא גם עבור תוכניות ניטור החוף באזורים אחרים מושפעים HAB ברחבי העולם.
למרות פרוטוקול זה סיפק את היתרונות הנ”ל, יש לדון בחסרונות של השיטה. כפי שניתן לראות בפרוטוקול החזותי, שיטת המטא-ברקודינג גוזלת זמן ומורכבת, ודורשת מכונות יקרות וריאגנטים. יש להכשיר את אנשי המעבדה באופן מיוחד, או שהם יבזבזו חומר, עבודה וזמן. בנוסף, זיהוי אצות על ידי metabarcoding חייב להיות מזווג עם מיקרוסקופיה כדי לוודא כי אותם מינים דומיננטיים מזוהים משתי שיטות אורתוגונל29. מיקרוסקופיה היא כלי לא פולשני, על פי רוב, לזיהוי מינים אצות, כלומר קשה יותר לעשות טעויות כאשר למינים האצות יש צורות ייחודיות לכאורה. כמובן, טעות אנוש יכולה להתרחש אם למין יש צורות דומות מאוד זו לזו. מצד שני, מאז ניתוח metabarcoding דורש צעדים מרובים, זה באופן טבעי מגביר שגיאות. זה יכול להיות מדגם מעורבב, תוספת ריאגנט שגויה, או חסר כמה הליכים. לכן, זה קריטי להשוות את תוצאות metabarcoding עם אלה שהושגו על ידי מיקרוסקופיה. לבסוף, הפרוטוקול ישים רק מבחינה איכותית, ויש לחשב את התוצאות עבור שפע יחסי. כפי שציין Lamb et al. בשנת 2019, המטא-ברקוד הנוכחי מכיוון שהביצועים הכמותיים עדיין מוגבלים, ונדרשים מחקרים נוספים לפני שניתן להשתמש בביטחון במטה-ברקודינג עבור יישומים כמותיים30.
פרוטוקול זה עבר אופטימיזציה עבור 140 – 170 דוגמאות לריצה ממדריך הכנת ספריית הרצף המטגנומי 16S שהונפק על ידי Illumina Inc. הפרוטוקול הממוטב נבדק פעמים רבות ושונה עוד יותר כדי להנפיק גירסה סופית זו. לכן, מומלץ מאוד לעקוב אחר כל צעד בדיוק. החלק הקריטי ביותר בפרוטוקול הוא שנדרש טיפול נוסף כדי למנוע זיהום מדגם. יש לנקות את הפיפטות ואת מכסה המנוע לזרם למינאר עם חשיפה של 70% לאלכוהול ול-UV בכל עת, ויש לעקר חומרים הניתנים לחיטוי אוטומטי. פריימרים ריאגנטים צריך להיות מדולל ישירות מהמלאי בכל ריצה חדשה, או שימוש חוזר של ריאגנטים וחשיפה למספר דילולים יכול להיות גורם זיהום מדגם. הפרוטוקול מציין מתי ניתן לבצע את הפעולה מחוץ למכסה המנוע של זרימת למינאר. אלא אם כן אחרת, יש לטפל בדגימות במכסה זרימה למינארי נקי. כאשר מתמודדים עם רצועות 8-tube מרובות בו זמנית, מומלץ לעבוד על הרצועה הראשונה 8-tube וכובע זה לפני המעבר לרצועת 8 הצינור הבאה. השארת המכסה פתוח במשך זמן רב עלולה לגרום לזיהום מדגם שכן גנים rRNA 16S ו- 18S נמצאים בכל מקום בסביבה. יש לעקוב אחר הזמן האמור, כגון זמן ערבוב וזמן הדגירה, כמתואר בדיוק מכיוון שאורך הזמן הטוב ביותר עבור כל שלב נבחר בהתבסס על ריצות הבדיקה המרובות. עודף או זמן לא מספיק יכול, למשל, להפחית את התפוקה לדוגמה. יש לעצור את התהליך רק בחלק שהפרוטוקול אומר שהוא יכול להיפסק. זה חיוני כדי לקבל שליטה שלילית כפי שהוא מאשר כי התוצאות החיוביות אינן חיוביות שווא מלאכותיות. לבסוף, מומלץ מאוד לתכנן את הימים כי זה דורש לפחות חמישה ימים כדי לעבד את כל השלבים, וחלקים מסוימים לא ניתן לעצור עד סימן העצירה הבאה.
המגבלה של הקצאה טקסונומית עבור נתונים רצפים נרשמת כאן בקצרה. רצפי נוקלאוטידים שהתגלו לאחרונה עבור מינים חיידקיים ואצות מתעדכנים במאגר הנתונים מדי יום. בעוד מינים אצות וחיידקים שנחקרו היטב רשומים באופן אמין, ישנם גם רצפים רבים עבור מינים לא ידועים המתעדכנים במאגר הנתונים. הוא מציין כי הפתרון של רצפי נוקלאוטידים רשומים אינו תמיד מספיק לזיהוי מינים אלא רק לזיהוי סוג. לכן, זיהוי מינים מיקרוסקופיים חייב להתבצע אורתוגונית. ליצירת צינור גישה פתוחה לעבודה זו יש יתרון משמעותי בהתמודדות עם קבוצה גדולה של נתונים רצפים בבת אחת ובחקירה יעילה כאשר מתרחשת שגיאה. בנוסף, הפלט של DADA2 הוא בקובץ טקסט או csv, מה שמקל על ניתוח סטטיסטי נוסף. מצד שני, שרת גדול נדרש לבצע הקצאות טקסונומיות, במיוחד כאשר ערכת נתונים גדולה. כדי להקים את הצינור, מהנדסים נדרשים כדי להגדיר את הצינור, ואנשי המקצוע חייבים להבין פרמטרים, פעולה, לינוקס וביואינפורמטיקה.
מלבד ההיקף ככלי ניטור HAB, ניתן להרחיב את השימוש בפרוטוקול זה למטרות מחקר חקירה. לדוגמה, מתנהלים ויכוחים מתמשכים בין ממשלת צ’ילה לדייגים מקומיים, יחד עם עמותות לשלום סביבתי בנוגע לעלייה במספרם שלHAB 31,32. הממשלה טוענת כי הסיבה לכך היא תופעת טבע כמו התחממות כדור הארץ ואל ניניו, בעוד הצדדים המאוחרים יותר טוענים כי חקלאות ימית סלמון היא הגורם. סלמון אינו מין ילידי צ’ילה ולא היה במים צ’יליאניים לפני חצי מאה. ממשלת צ’ילה באותה תקופה שמה לה למטרה להצמיח את הכלכלה וליצור מקומות עבודה לאזורים עניים על ידי פיתוח עסק סלמון33. עם התערבות בחו”ל לרווח ההון, צ’ילה הצליחה מאוד בפיתוח חקלאות ימית בסגנון עט, ומספר הסלמון גדל להפליא בעשורים האחרונים31,32,33. לאחר מכן, כמות גדולה של מזון לסלמון נזרקה לים, מה שמוביל את המקומיים לחשוד כי חקלאות ימית סלמון היא הסיבה HABs המקומי התכוף31,32. האמת אינה ידועה כעת, אך יש להבין אותה בסופו של דבר, כך שניתן יהיה לתכנן אסטרטגיות להגנה על הסביבה הימית המקומית, הכלכלה ובריאות האדם מפני HABs. הניתוח המולקולרי על מערכת היחסים האצות-חיידקית המקומית עשוי לתרום להבהרה צעד קדימה לנושא כזה. לדוגמה, הטכניקה יכולה לחפש את מיני האצות והחיידקים שלא היו נוכחים באזור הימי היעד לפני כן אך גדלו באופן דרמטי בשנים האחרונות, וזה דומה במקצת למחקר שנעשה על ידי Sakai et al. לגילוי אוכלוסייה לא מתועדת של סלמנדרה יאמאטו (Hynobius vandenburghi) על ידי GIS וניתוח eDNA34. לכן, מעבר לכלי ניטור HAB, פרוטוקול metabarcoding זה יש את השימוש הפוטנציאלי עבור לקוחות פוטנציאליים אחרים.
The authors have nothing to disclose.
מחקר זה נתמך על ידי המענק (JPMJSA1705) למחקר על שותפות מחקר מדע וטכנולוגיה עבור אצות ניטור פיתוח בר קיימא בצ’ילה (SATREPS-MACH). אנו מודים לד”ר סנדרה ריוס (תוני לוס לאגוס) על שאפשרה לנו להשתמש בקליפ. אנו מודים לניל אנדרו הולנד על הגהתו החרוצה של כתב היד. אנו מביעים את הערכתנו הכנה לחברי קבוצת המעבדה ב- CREAN-IFOP בפוארטו מונט, צ’ילה, על כך שייעצו לנו על זיהוי פיטופלנקטון.
1.5 mL single tube | ThermoFisher | Q32856 | sterile |
1.5 mL tubes | ThermoFisher | 2150N | sterile |
8-strip 0.2 mL PCR tubes and flat cap | ThermoFisher | AB0451 & 4323032 | sterile |
96-well 0.2 mL PCR plate | ThermoFisher | N8010560 | |
AccuRuler 100 bp DNA Ladder | MaestroGen | 02001-500 | |
Chelex 100 Chelating Resin | Bio-Rad | 1432832 | |
Chemical Duty Vacuum Pressure Pump | MilliporeSigma | WP6122050 | |
D1000 Reagents | Agilent Technologies | 5067-5583 | |
D1000 sample buffer | Agilent Technologies | 5067-5602 | |
DNA loading dye (Blue/Orange Loading Dye 6X) | Promega | G1881 | |
Ethanol Molecular Biology Grade | E7023 | Merck KGaA | |
Filtration device | ThermoFisher | 09-740-36H, 09-740-23E | |
Fluorescent DNA concentration measurement kit (Qubit dsDNA HS Assay kit) | Thermo Fisher Scientific | Q32851 | |
Fluorometer (Qubit4 Fluorometer) | ThermoFisher | Q33238 | |
Fragment analysis verification matrix (D1000 ScreenTape) | Agilent Technologies | 5067-5582 | |
Fragment Analyzer (Agilent Tapestation 4150) | Agilent | G2992AA | |
GelRed Nucleic Acid Gel Stain | Biotium | 41002 | |
Generic bench vortexer (Vortex Genie 2) | Scientific Industries | SI-0236 | |
Heat block | |||
High performance sequencing-grade DNA polymerase (KAPA HiFi Hotstart ReadyMix), 2X | Roche | KK2602 | |
HT1 Hybridization buffer | Illumina | 20015892 | |
Illumina Nextera XT v2 Index Kit set A, B, C and D | Illumina | FC-131-2001, -2002, -2003, -2004 | |
Laminar hood cabinet | Biobase Co. | Biobase PCR-800 PCR Cabinet | |
Loading tips (Specific for Agilent TapeStation systems) | Agilent Technologies | 5067-5598 | |
Magnetic beads DNA cleanup system (ProNex® Size-Selective Purification System) | Promega | NG2002 | |
Magnetic stand for 96 wells (Magnetic Stand-96) | ThermoFisher | AM10027 | |
Micro-seal film for 96-well plate | ThermoFisher | 4311971 | sterile |
MiSeq Reagent Kit v3 | Illumina | MS-102-3003 | includes pre-filled, ready-to-use reagent cartridges. |
MiSeq system | Illumina | SY-410-1003 | includes pre-installed software |
Molecular grade nuclease-free water | Integrated DNA Technologies | 11-05-01-04 | |
Molecular grade sodium hydroxide (NaOH) 1M | Merck KGaA | 1091371000 | |
Multichannel pipettes | Not specified | Not specified | |
Multi-Purpose Agarose | Cleaver Scientific | CSL-AG500 | |
Ow D2 Wide-Gel Electrophoresis System | ThermoFisher | D2 | |
Ow EC-105 Compact Power Supply | ThermoFisher | 105ECA-115 | |
Pellet Pestle— Cordless Motor | ThermoFisher | K749540-0000 | |
PhiX control Kit v3 | Illumina | FC-110-300 | ready-to-use control library for Illumina sequencing |
Pipette tips | Not specified | Not specified | sterile |
Pipettes | Not specified | Not specified | 2, 10, 20, 100, 1000 μL |
SterivexTM GP 0.22 μm filter unit | MilliporeSigma | SVGP01050 | |
Tabletop centrifuge | Eppendorf 5427R Centrifuge | Eppendorf AG | |
TBE buffer | ThermoFisher | B52 | |
TE buffer (pH 8.0) | ThermoFisher | AM9849 | |
Terr PCR Direct FFPE Kit | Takara Bio USA | 639284 | includes 2X Terra PCR Direct Buffer |
Terra PCR Direct Polymerase Mix, 1.25 U/μL | Takara Bio USA | 639271 | High performance Inhibitor-tolerant DNA polymerase |
ThermalCycler (MiniAMP Plus Thermal Cycler) | ThermoFisher | A37835 | |
TransIlluminator and image capture system | Analytik Jena, AG | GelDoc-ItTS2 Imager | |
Whatman 1.0 m pore-sized membrane | MilliporeSigma | WHA111110 |