Dit protocol is ontwikkeld en getest met behulp van monsters die zijn verzameld als onderdeel van studies goedgekeurd door de UCSF Institutional Review Board.

Figuur 1: Cerebrospinale vloeistof (CSF) metagenomische next-generation sequencing (mNGS) monsterverwerking en analytische workflow voor pathogeendetectie. Dit diagram toont de belangrijkste stappen in de CSF mNGS wet-lab en bio-informatische routes. Klik hier om een grotere versie van deze figuur te bekijken.
1. Overwegingen bij het behandelen van CSF-monsters
OPMERKING: mNGS-resultaten zijn optimaal als de voorbereiding van de CSF-bibliotheek plaatsvindt tijdens het verzamelen, of als CSF-monsters worden bewaard via snap-freezing in vloeibare stikstof, gevolgd door langdurige opslag bij -80 °C15. Als alternatief kan het gebruik van conserverende oplossingen zoals DNA/RNA-schild bij het verzamelen van monsters langere tijden bij kamertemperatuur mogelijk maken, namelijk16,17. Dit kan te verkiezen zijn in situaties waarin er zorgen zijn over de betrouwbaarheid van de koelketen. Daarnaast kan het vermijden van onnodige vries-dooicycli ook helpen om de integriteit van nucleïnezuren te behouden, en het kan raadzaam zijn om monsters vooraf te aliquoteren op het moment van verzameling als er meerdere analytische workflows worden verwacht voor de CSF-monsters18,19.
2. Extractie van CSF-nucleïnezuur
VOORZICHTIG: Voer uit in een geschikte bioveiligheidskast met handschoenen en een laboratoriumjas om het risico op overdracht van infectieuze organismen te verminderen. Reagentia zijn schadelijk als ze worden doorgeslikt, ingeademd of in contact komen met huid. Verwijder gebruikte materialen als gevaarlijk chemisch afval in overeenstemming met lokale regelgeving.
- De Quick-DNA/RNA Kit werkt over het algemeen goed voor monsters met een lage hoeveelheid nucleïnezuren.
- Ontdooi de CSF-monsters op ijs, doe vervolgens 100 μL–1 mL van elk monster over in een 1,5 mL buis en centrifugeer op 16.000 x g bij 4 °C gedurende 10 minuten.
- Pipette voorzichtig het supernatant af, waarbij de pellet (die vaak onzichtbaar is voor het blote oog) en 100 μL supernatant achterblijft.
- Voeg 100 μL DNA/RNA Shield toe aan de resterende 100 μL supernatant en pellet, en meng grondig met de pipet.
- Voeg vervolgens 600 μL DNA/RNA lysisbuffer en pipetmengsel toe totdat de oplossing homogeen is.
- Volg de rest van het DNA/RNA-extractieprotocol zoals beschreven door de fabrikant in aanvullende bijlagen A, met de volgende belangrijke aanpassingen specifiek voor CSF:
- Voer alle centrifugatiestappen uit gedurende 1 minuut op maximale snelheid (~21.000 x g) behalve de laatste DNA/RNA-wasstap, die 3 minuten op maximale snelheid moet duren.
- Eluteer het monster door 22 μL (in plaats van 25 μL) nucleasevrij water direct aan de kolommatrix toe te voegen, en incubeer vervolgens 3 minuten bij kamertemperatuur. Het hogere volume resulteert doorgaans in een grotere totale hoeveelheid nucleïnezuur dat uit de kolom wordt geëluerd.
- Na eluatie pipetteer je de doorstroom-door terug op de kolommatrix en centrifugeer je opnieuw een minuut op maximale snelheid. Breng het monster over naar een nieuwe 1,5 mL opslagbuis en ga direct over tot bibliotheekvoorbereiding of opslag bij -80 °C.
3. CSF RNA-bibliotheekvoorbereiding
VOORZICHTIG: Voer uit in een geschikte bioveiligheidskast, met handschoenen en een laboratoriumjas. Reagentia zijn schadelijk als ze worden doorgeslikt, ingeademd of in contact komen met huid. Verwijder gebruikte materialen als gevaarlijk chemisch afval in overeenstemming met lokale regelgeving.
OPMERKING: Belangrijke overwegingen voorafgaand aan het begin: rRNA-depletie heeft de voorkeur boven polyadenylatieverrijking voor isolatie van mRNA, aangezien de rRNA-depletie geen extra kraalreinigingsstappen vereist (wat de opbrengst van bruikbaar mRNA kan verminderen), en sommige neuroinvasieve RNA-virustranscripten geen polyadenylatie hebben. ERCC (External RNA controls consortium) RNA Spike-In kan worden toegevoegd als kwaliteitscontrolemaatregel, hoewel dit niet strikt noodzakelijk is. Als het RNA na extractie bij -80 °C is opgeslagen, ontdooi dan op ijs voordat het begint.
- RNA-fragmentatie en priming: Pipet-mix het volgende in een polymerasekettingreactie (PCR) buis: 3,5 μL RNA-monster, 0,5 μL 1:2500 verdund ERCC spike-in, 4 μL eerste-streng reactiebuffer (5x), 1 μL willekeurige primers en 1 μL 1:100 verdund rRNA (Table of Materials). In een thermocycler zet je het verwarmde deksel op 105 °C en incubeer je de monsters op 75 °C gedurende 2 minuten, 70 °C voor 2 minuten, 65 °C voor 2 minuten, 60 °C voor 2 minuten, 55 °C voor 2 minuten, 37 °C voor 5 minuten en 25 °C voor 5 minuten.
- Eerstestreng cDNA-synthese: Pipette mengt het gefragmenteerde en geprimeerde RNA (10 μL) met 8 μL nucleasevrij water en 2 μL eerstestrengs synthese-enzymmix. In een thermocycler zet je het verwarmde deksel op 105 °C en becubeer je de monsters op 25 °C gedurende 10 minuten, 42 °C gedurende 15 minuten en 70 °C gedurende 15 minuten.
- Tweede-streng cDNA-synthese: Pipettenmix het eerste-streng gesynthetiseerde DNA (20 μL) met 8 μL tweede-streng synthese reactiebuffer, 4 μL tweede-streng synthese enzymmix en 48 μL nuclease-vrij water. In een thermocycler incubeer je de monsters 1 uur bij 16 °C, met het verwarmde deksel uitgeschakeld.
- Voer een magnetische kraalzuivering uit (aanvullende bijlagen B) met een verhouding van 1,8x SPRI magnetische kralen (144 μL). Eluteer in 53 μL nucleasevrij water en breng 50 μL van de uiteindelijke supernatant (met gezuiverde dubbelstrengs cDNA) over in een schone nucleasevrije PCR-buis. Op dit punt kunnen de monsters veilig worden ingevroren bij -20 °C 's nachts indien nodig.
- cDNA-bibliotheek eindvoorbereiding: Als het monster 's nachts bij -20 °C is bewaard, ontdooi dan op ijs voordat het opnieuw wordt gestart. Pipette mengt het gezuiverde ds-cDNA (50 μL) met 7 μL eindprep-reactiebuffer en 3 μL eindprep-enzymmix.
- In een thermocycler incubeer je de monsters bij 20 °C gedurende 30 minuten en 65 °C gedurende 30 minuten, met het verwarmde deksel gesloten.
- Adapterligatie (voer deze stap uit op ijs): Maak een verdunning van de adapter van 1:100 in nucleasevrij water. Pipetteer het eindpreparatreactiemengsel (60 μL) in een buis met 30 μL ligatie-meestermix, 1 μL ligatieversterker en 2,5 μL 1:100 verdunde adapter. Zorg ervoor dat de adapter apart wordt toegevoegd (bijvoorbeeld niet vooraf gemengd met de ligatie-mastermix of ligatie-enhancer) om de vorming van adapterdimeren te voorkomen.
- In een thermocycler incubeer je de monsters bij 20 °C gedurende 15 minuten, met het verwarmde deksel uitgeschakeld.
- Ga direct over tot magnetische kralenzuivering (Aanvullende bijlage 2) met een verhouding van 0,9x SPRI magnetische kralen (87 μL).
- Eluteer in 17 μL nucleasevrij water en breng 15 μL van het uiteindelijke supernatant over naar een schone nucleasevrije PCR-buis
- Barcoding PCR: Pipett het gezuiverde, adapter-geligeerde cDNA (15 μL) in een buis met 3 μL USER-enzym, 25 μL Q5 Master Mix en 10 μL unieke barcodeprimers. In een thermocycler zet je het verwarmde deksel op 105 °C, breng je de monsters uit op 37 °C gedurende 15 minuten, daarna 98 °C gedurende 30 s, en voer je vervolgens 19 cycli uit van 98 °C voor 10 s en 65 °C voor 75 s. Rond de PCR af door 5 minuten te incuberen op 65 °C.
- Voer een laatste magnetische kraalzuivering uit met een verhouding van 0,8x magnetische kralen (43 μL). Eluteer in 23 μL nucleasevrij water en breng na de laatste parelscheidingsstap 20 μL over in een schoone, nucleasevrije PCR-buis. Op dit punt kunnen de monsters veilig worden ingevroren bij -20 °C 's nachts indien nodig.
4. CSF DNA-bibliotheekvoorbereiding
VOORZICHTIG: Werk in een geschikte kast met een bioveiligheidsniveau met handschoenen en een laboratoriumjas. Reagentia zijn schadelijk als ze worden doorgeslikt, ingeademd of in contact komen met huid. Verwijder gebruikte materialen als gevaarlijk chemisch afval in overeenstemming met lokale regelgeving.
- Als het geëxtraheerde DNA 's nachts bij -20 °C is opgeslagen, ontdooi dan op ijs voordat je opnieuw begint.
- DNA-fragmentatie en eindvoorbereiding: Meng de buffer voor de reactie van de eerste streng grondig door te vortexen en de oplossing te mengen met een pipet om eventuele neerslag opnieuw te laten ophangen. Pipetmix bevat 3,5 μL van het geëxtraheerde DNA, 22,5 μL nucleasevrij water, 7 μL eerste-streng reactiebuffer en 2 μL eerste-streng enzymmix.
- In een thermocycler zet je het verwarmde deksel op 105 °C en incubeer je de monsters op 37 °C gedurende 5 minuten, daarna op 65 °C gedurende 30 minuten.
- Het adapterligatie en het barcode PCR-proces is identiek voor de DNA- en RNA-bibliotheken. Herhaal secties 7–10 van de "CSF RNA Library Preparation" met het gefragmenteerde en eindvoorbereide DNA om de DNA-bibliotheekvoorbereiding te voltooien.
5. Kwaliteitscontrole en pooling van bibliotheken
- Kwantificeer de concentratie van de monsterbibliotheken via de DNA-kwantificatiekit en de bijbehorende fluorometer. Watercontroles moeten aanzienlijk lagere (of niet kwantificeerbare) concentraties hebben vergeleken met de rest van de monsters.
- Beoordeel de grootte van de monsterbibliotheek door de monsters te laten draaien op een geautomatiseerde capillaire elektroforesemachine.
- Als er als alternatief geen geautomatiseerde capillaire elektroforesemachine beschikbaar is, versterk dan 1-2 μL van het monster verder door een PCR-reactie met Illumina universele primers uit te voeren, en vervolgens een elektroforese van 10 μL van het eindproduct op een 2% agarosegel. Verwijder gebruikte materialen als gevaarlijk chemisch afval in overeenstemming met lokale regelgeving.
OPMERKING: De optimale bibliotheeklengte is ongeveer 400–500 basenparen, dus toekomstige fragmentatie-incubaties moeten dienovereenkomstig worden aangepast als de bibliotheekgroottes groter of korter zijn dan gewenst. Daarnaast is het cruciaal om te beoordelen op de aanwezigheid van adapterdimeren, die ongeveer 150 basenparen lang zijn en optreden wanneer twee adaptermoleculen per ongeluk met elkaar samengaan zonder insertsequentie. Zelfs als ze in kleine hoeveelheden aanwezig zijn, clusteren en sequencen adapterdimers efficiënt op de flowcel, waardoor het aandeel bruikbare reads tijdens een sequencingrun afneemt. Als er adapterdimere aanwezig zijn, voer dan een extra magnetische kraalzuivering uit met een verhouding van 0,8x SPRI-magnetische kralen. (Opmerkelijk is dat als >10% van het monster adapterdimeer is, het mogelijk 2-3 rondes magnetische parelzuivering kan kosten om de meeste adapterdimeren te verwijderen).
- Om een vergelijkbare diepte van sequencing over monsters binnen dezelfde sequencingrun te garanderen, is het noodzakelijk om equimolaire hoeveelheden van elk monster samen te voegen. Als de bibliotheekgroottes vergelijkbaar zijn over alle monsters, voeg dan simpelweg equivalente massa's van elke RNA-bibliotheek samen en herhaal dit afzonderlijk voor elke DNA-bibliotheek.
6. Volgordeoverwegingen
- De gepoolde monsters zijn nu klaar voor sequencing op Illumina-sequencingapparaten. Als je intern sequenceert (in plaats van via een gespecialiseerde sequencing core-faciliteit), denatureer en verdun dan zorgvuldig de gepoolde bibliotheken volgens de instructies van het sequencing-apparaat voordat je de pool op de reagentcartridge laadt (de typische concentratie is 4 nM voor de meeste Illumina-sequencingplatforms).
- Bereken de geschatte sequencingdiepte door de verwachte dataopbrengst voor de Illumina-flowcel te delen door het aantal gemultiplexeerde monsters in de run.
OPMERKING: Het bepalen van de gewenste sequencingdiepte hangt af van verschillende factoren, waaronder kosten, het gebruikte sequencing-apparaat en de doelstellingen van het experiment. Twee voorbeelden van verschillende sequencingdieptes worden hieronder gepresenteerd in de sectie Representatieve Resultaten.
7. Data-analyse voor het opsporen van ziekteverwekkers
- Maak gebruik van het open-source webgebaseerde platform Chan Zuckerberg ID (CZID, https://www.czid.org, Illumina mNGS Pipeline v8.3) om de metagenomische analyse voor het opsporen van ziekteverwekkers uit te voeren. CZID bevat ook opties om samples te uploaden via de commandoregelinterface en directe overdracht van een Illumina BaseSpace-account naar CZID. Een gedetailleerd overzicht van de analyse van CZID mNGS-gegevens is beschikbaar op https://chanzuckerberg.zendesk.com/hc/en-us/articles/13770737266196-Guide-to-mNGS-Data-Analysis
- Log in bij CZID, klik op "upload", selecteer "metagenomics" voor analysetype en selecteer de input fastq-bestanden uit de sequencing-run. Voer de benodigde informatie in voor elk monster, inclusief monsternaam, monstertype (CSF) en nucleotide (RNA of DNA).
OPMERKING: De geüploade bestanden ondergaan geautomatiseerde verwerking via de CZID-pijplijn, die verschillende computationele tools bevat zoals STAR, Bowtie2, Trimmomatic, PriceSeq en GSNAP, onder andere, om lage kwaliteit leesopdrachten en menselijke leessequenties te verwijderen. De pijplijn lijnt vervolgens de gefilterde sequenties uit en assembleert ze met behulp van de nucleotide (NT) en eiwit (NR) databases van het National Center for Biotechnology Information (NCBI) via respectievelijk de softwareprogramma's Minimap2 en DIAMOND. De uiteindelijke CZID-output omvat NT/NR-taxontellingen en contig-tellingen (aaneengesloten, overlappende segmenten die tijdens de assemblage worden gegenereerd) voor elk monster.
- Na verwerking klik je op de projectmap met de voorbeelden en ga je naar het samenvattingsdashboard.
- Scroll naar beneden om het aantal reads per monster te bekijken (en vergelijk het met het verwachte aantal reads op basis van de gebruikte Illumina-flowcel en het aantal monsters dat in de sequencing run is gemultiplexed), het percentage reads dat door kwaliteitscontrolefiltering komt, en de dubbele compressieverhouding om te beoordelen op biased PCR-overamplificatie.
- De aanwezigheid van besmetting door nucleïnezuurfragmenten uit de omgeving - verkregen tijdens het CSF-monsterverzamelingsproces, hetzij tijdens extractie/bibliotheekvoorbereiding - is een veelvoorkomend verschijnsel. Om te bepalen of bacteriële, schimmel- en parasitaire reads echte positieven zijn, kunnen mogelijke methoden worden overwogen die gebruikmaken van de genormaliseerde readcounts (reads per million, of rPM):
- Maak een achtergrondmodel (bijvoorbeeld met watermonsters die zijn geëxtraheerd, voorbereid en gesequenced) direct binnen CZID om taxonen met rPM's die vergelijkbaar zijn met of lager zijn in hoeveelheid ten opzichte van het achtergrondmodel (en dus waarschijnlijk verontreiniging vertegenwoordigen) eruit te filteren. Controleer binnen de projectmap elk van de watercontrolemonsters en klik vervolgens op "achtergrondmodel" om het aan te maken. Selecteer bij het analyseren van elk monster het achtergrondmodel en klik op "threshold filters" om een "NTZ Score" > 1 (of hoger) toe te passen om waarschijnlijke verontreinigingen eruit te filteren.
- Gebruik vooraf bepaalde rPM-drempels, zoals rpM-monster > 10* (gemiddelde rPM-controles), of log10 (rPM-monster) minstens 1 log groter dan log 10 (gemiddelde rPMhele cohort). Download eerst de taxon-specifieke rPM voor elk monster door elk monster op de projectpagina te controleren, op downloaden te klikken en "Gecombineerde Steekproef Taxon Resultaten" te selecteren. Het resulterende .csv-bestand kan naar wens worden geanalyseerd met statistische software zoals R of Stata.
OPMERKING: Gezien de lagere relatieve hoeveelheid virale pathogenen, moet elk virus met bekend neuroinvasief potentieel en ten minste één read die aansluit bij het virale genoom als positief worden beschouwd (en idealiter bevestigd via herhaalde mNGS of pathogen-specifieke klinische tests). Daarnaast vereist de analytische pijplijn dat elke lezing in het fastq-bestand computationeel is uitgelijnd met het meest waarschijnlijke referentiegenoom in de NCBI-database. Daardoor kunnen sommige reads relatief niet-specifiek zijn voor een bepaald pathogeen en theoretisch kunnen ze zijn uitgelijnd op een breed scala aan organismen in plaats van het specifieke taxon dat in CZID is toegewezen. Gezien dit risico moet handmatige bevestiging van positieve calls worden uitgevoerd door de uitlijning te bevestigen via de webgebaseerde NCBI BLAST-tool (de specifieke reads voor een taxon kunnen direct van CZID naar NCBI BLAST worden verzonden door op het blast-icoon naast het betreffende taxon te klikken).