Ce protocole a été développé et testé à partir d’échantillons collectés dans le cadre d’études approuvées par le comité d’examen institutionnel de l’UCSF.

Figure 1 : Traitement des échantillons de séquençage de nouvelle génération métagénomique du liquide céphalorachidien (LCR) et flux de travail analytique pour la détection des agents pathogènes. Ce diagramme illustre les principales étapes des voies de laboratoire humide et bioinformatique du MFC mNGS. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
1. Considérations sur la manipulation des échantillons de LCR
REMARQUE : les résultats du mNGS sont optimaux si la préparation de la bibliothèque de LCR a lieu au moment de la collecte, ou si les échantillons de LCR sont conservés par congélation instantanée dans de l’azote liquide, suivie d’un stockage à long terme à -80 °C 15. Alternativement, l’utilisation de solutions de conservation telles que la protection ADN/ARN au moment de la collecte de l’échantillon peut permettre un stockage à température ambiante pendant des périodes pluslongues 16,17. Cela peut être préférable dans des situations où il existe des inquiétudes quant à la fiabilité de la chaîne du froid. De plus, éviter des cycles de congélation-dégel inutiles peut également aider à préserver l’intégrité des acides nucléiques, et il peut être conseillé de pré-pré-aliquoter les échantillons au moment de la collecte si plusieurs flux de travail analytiques sont anticipés pour les échantillonsde LCR 18,19.
2. Extraction de l’acide nucléique LCR
ATTENTION : Agir dans une armoire de niveau de biosécurité appropriée, avec des gants et une blouse de laboratoire afin de réduire le risque de transmission d’organismes infectieux. Les réactifs sont nocifs s’ils sont avalés, inhalés ou entrent en contact avec la peau. Éliminer les matériaux usagés comme déchets chimiques dangereux conformément à la réglementation locale.
- Le kit Quick-DNA/ARN fonctionne généralement bien pour des échantillons à faible abondance d’acides nucléiques.
- Décongelez les échantillons de LCR sur glace, puis transférez 100 μL–1 mL de chaque échantillon dans un tube de 1,5 mL et centrifugez à 16 000 x g à 4 °C pendant 10 minutes.
- Faites soigneusement une pipette sur le surnageant, en laissant derrière elle la pellete (souvent invisible à l’œil nu) et 100 μL de surnageant.
- Ajoutez 100 μL de bouclier ADN/ARN aux 100 μL restants de surnageant et de pellets, et faites un mélange de pipettes soigneusement.
- Ensuite, ajoutez 600 μL de mélange tampon de lyse ADN/ARN et pipette jusqu’à ce que la solution soit homogène.
- Suivez le reste du protocole d’extraction ADN/ARN tel que décrit par le fabricant dans les annexes complémentaires A, avec les modifications importantes suivantes spécifiques au LCR :
- Effectuez toutes les étapes de centrifugation pendant 1 minute à vitesse maximale (~21 000 x g) sauf l’étape finale de lavage ADN/ARN, qui devrait durer 3 minutes à vitesse maximale.
- Éluez l’échantillon en ajoutant 22 μL (au lieu de 25 μL) d’eau sans nucléase directement sur la matrice colonne, puis incubez pendant 3 minutes à température ambiante. Le volume plus élevé entraîne généralement une plus grande quantité totale d’acide nucléique élue de la colonne.
- Après élution, pipettez le flux à travers la matrice de colonne et centrifugez à nouveau pendant une minute à vitesse maximale. Transférez l’échantillon dans un nouveau tube de stockage de 1,5 mL, puis passez immédiatement à la préparation de la bibliothèque ou stockez à -80 °C.
3. Préparation de la bibliothèque ARN du LCR
ATTENTION : Agir dans un cabinet de biosécurité approprié, en portant des gants et une blouse de laboratoire. Les réactifs sont nocifs s’ils sont avalés, inhalés ou entrent en contact avec la peau. Éliminer les matériaux usagés comme déchets chimiques dangereux conformément à la réglementation locale.
REMARQUE : Considérations importantes avant de commencer : l’épuisement de l’ARNr est préféré à l’enrichissement par polyadénylation pour l’isolement de l’ARNm, étant donné que la déplétion de l’ARNr ne nécessite pas d’étapes supplémentaires de nettoyage des billes (ce qui peut réduire le rendement d’ARNm utilisable), et que certains transcrits de virus à ARN neuroinvasifs ne présentent pas de polyadénylation. Le SPIKE-In de l’ARN (ERCC (External RNA controls consortium) peut être ajouté comme mesure de contrôle qualité, bien que ce ne soit pas strictement nécessaire. Si l’ARN a été stocké à -80 °C après l’extraction, dégelez sur la glace avant le début.
- Fragmentation et amorçage de l’ARN : Pipettez mélanger les éléments suivants dans un tube de réaction en chaîne polymérase (PCR) : échantillon d’ARN de 3,5 μL, spike-in dilué 1:2500 de 0,5 μL, tampon de réaction de première brin de 4 μL (5x), amorces aléatoires de 1 μL, et 1 μL d’ARNr dilué 1:100 (Table of Materials). Dans un thermocycleur, réglez le couvercle chauffé à 105 °C et faites couver les échantillons à 75 °C pendant 2 minutes, 70 °C pendant 2 minutes, 65 °C pendant 2 minutes, 60 °C pendant 2 minutes, 55 °C pendant 2 minutes, 37 °C pendant 5 minutes et 25 °C pendant 5 minutes.
- Synthèse d’ADNc de première brin : Pipette mélange l’ARN fragmenté et amorcé (10 μL) avec 8 μL d’eau sans nucléase et 2 μL de mélange enzymatique de synthèse de première brin. Dans un thermocycleur, réglez le couvercle chauffé à 105 °C et faites couver, les échantillons à 25 °C pendant 10 minutes, 42 °C pendant 15 minutes et 70 °C pendant 15 minutes.
- Synthèse d’ADNc de second brin : Pipette mélangez l’ADN synthétisé de première brin (20 μL) avec un tampon de réaction de synthèse de deuxième brin de 8 μL, un mélange enzymatique de synthèse de deuxième brin de 4 μL, et 48 μL d’eau sans nucléase. Dans un thermocycleur, incubez les échantillons pendant 1 heure à 16 °C, avec le couvercle chauffé éteint.
- Effectuer une purification des billes magnétiques (Annexes supplémentaires B) en utilisant un ratio de 1,8x des billes magnétiques SPRI (144 μL). Éluvez dans de l’eau sans nucléase de 53 μL et transférez 50 μL du surnageant final (contenant de l’ADNc double brin purifié) dans un tube PCR propre sans nucléase. À ce stade, les échantillons peuvent être congelés en toute sécurité à -20 °C toute la nuit si nécessaire.
- Fin de la préparation de la bibliothèque d’ADN : Si l’échantillon a été conservé à -20 °C toute la nuit, décongelez sur glace avant de recommencer. Pipette, mélangez le ds-cDNA purifié (50 μL) avec 7 μL de tampon de réaction de préparation finale et 3 μL de mélange enzymatique de préparation terminale.
- Dans un thermocycleur, incubez les échantillons à 20 °C pendant 30 minutes, et à 65 °C pendant 30 minutes, avec le couvercle chauffé éteint.
- Ligation de l’adaptateur (effectuer cette étape sur glace) : Créer une dilution 1:100 de l’adaptateur dans de l’eau sans nucléase. Pipeter le mélange de la réaction de préparation finale (60 μL) dans un tube contenant un mélange maître de ligation de 30 μL, un amplificateur de ligature de 1 μL et un adaptateur dilué 1:100 de 2,5 μL. Assurez-vous que l’adaptateur est ajouté séparément (par exemple, pas prémélangé avec le mélange maître de ligation ou l’amplificateur de ligatura) pour éviter la formation de dimers adaptateurs.
- Dans un thermocycleur, incubez les échantillons à 20 °C pendant 15 minutes, avec le couvercle chauffé éteint.
- Procéder immédiatement à la purification des billes magnétiques (Annexe supplémentaire 2) en utilisant un ratio 0,9x des billes magnétiques SPRI (87 μL).
- Éluez dans de l’eau sans nucléase de 17 μL, et transférez 15 μL du surnageant final dans un tube PCR propre sans nucléase
- PCR par code-barres : Pipeter l’ADNc purifié, ligué par un adaptateur (15 μL) dans un tube contenant l’enzyme USER de 3 μL, le Master Mix Q5 de 25 μL et des amorces uniques codées à barres de 10 μL. Dans un thermocycleur, réglez le couvercle chauffé à 105 °C, incubez les échantillons à 37 °C pendant 15 minutes, puis à 98 °C pendant 30 s, puis effectuez 19 cycles à 98 °C pendant 10 s et 65 °C pendant 75 s. Terminez la PCR en incubant à 65 °C pendant 5 minutes.
- Effectuer une purification finale des billes magnétiques en utilisant un ratio 0,8x des billes magnétiques (43 μL). Éluez dans de l’eau sans nucléase de 23 μL, puis après la dernière étape de séparation des billes, transférez 20 μL dans un tube PCR propre sans nucléase. À ce stade, les échantillons peuvent être congelés en toute sécurité à -20 °C toute la nuit si nécessaire.
4. Préparation de la bibliothèque d’ADN du LCR
ATTENTION : Intervenez dans un cabinet de niveau de biosécurité approprié, avec des gants et une blouse de laboratoire. Les réactifs sont nocifs s’ils sont avalés, inhalés ou entrent en contact avec la peau. Éliminer les matériaux usagés comme déchets chimiques dangereux conformément à la réglementation locale.
- Si l’ADN extrait a été stocké à -20 °C toute la nuit, dégelez sur de la glace avant de recommencer.
- Fragmentation de l’ADN et préparation finale : Mélangez soigneusement le tampon de réaction de première brin en faisant du vortex et en pipette, mélangeant la solution pour resuspendre tout précipité. Mélangez 3,5 μL de l’ADN extrait, 22,5 μL d’eau sans nucléase, 7 μL de tampon de réaction première brin, et 2 μL de mélange enzymatique de première brin.
- Dans un thermocycleur, réglez le couvercle chauffé à 105 °C et faites éclore les échantillons à 37 °C pendant 5 minutes, puis à 65 °C pendant 30 minutes.
- Le processus de ligation de l’adaptateur et de code-barres PCR est identique pour les bibliothèques d’ADN et d’ARN. Répétez les sections 7 à 10 de la « Préparation de la bibliothèque d’ARN du LCR » avec l’ADN fragmenté et préparé à l’extrémité pour compléter la préparation de la bibliothèque d’ADN.
5. Contrôle qualité des bibliothèques et mise en commun
- Quantifiez la concentration des bibliothèques d’échantillons via le kit de quantification ADN et le fluoromètre correspondant. Les contrôles de l’eau devraient avoir des concentrations nettement inférieures (ou non quantifiables) par rapport aux autres échantillons.
- Évaluez la taille de la bibliothèque d’échantillons en les utilisant sur une machine automatisée d’électrophorèse capillaire.
- Sinon, si une machine automatisée d’électrophorèse capillaire n’est pas disponible, amplifiez davantage 1 à 2 μL de l’échantillon en effectuant une réaction PCR avec des amorces universelles Illumina, puis effectuez une électrophorèse de 10 μL du produit final sur un gel d’agarose à 2 %. Éliminer les matériaux usagés comme déchets chimiques dangereux conformément à la réglementation locale.
REMARQUE : La longueur optimale de la bibliothèque est d’environ 400 à 500 paires de bases, donc les futures incubations de fragmentation doivent être ajustées en conséquence si les tailles de la bibliothèque sont plus grandes ou plus courtes que souhaitées. De plus, il est crucial d’évaluer la présence de dimers adaptateurs, qui font environ 150 paires de bases et se produisent lorsque deux molécules adaptatrices se liguent accidentellement sans séquence d’insertion. Même s’ils sont présents en faible quantité, les dimères adaptateurs ont tendance à se regrouper et à séquencer efficacement sur la cellule de flux, réduisant la proportion de lectures utilisables lors d’un séquençage. Ainsi, si des dimères adaptateurs sont présents, effectuer une purification supplémentaire des billes magnétiques avec un ratio 0,8x des billes magnétiques SPRI. (À noter, si >10 % de l’échantillon est un dimère adaptateur, il peut falloir 2 à 3 cycles de purification magnétique par billes magnétiques pour retirer la plupart des dimères adaptateurs).
- Pour garantir une profondeur de séquençage similaire entre les échantillons au sein d’une même séquençage, il est nécessaire de regrouper les quantités équimolaires de chaque échantillon. Si les tailles des bibliothèques sont comparables sur tous les échantillons, il suffit de regrouper les masses équivalentes de chaque bibliothèque d’ARN et de répéter séparément pour chaque bibliothèque d’ADN.
6. Considérations de séquençage
- Les échantillons regroupés sont désormais prêts pour le séquençage sur les dispositifs de séquençage Illumina. Si le séquençage est effectué en interne (plutôt que via une installation spécialisée de séquençage), dénaturez et diluez soigneusement les bibliothèques mises en pool selon les instructions du dispositif de séquençage avant de charger le pool sur la cartouche de réactif (la concentration typique est de 4 nM pour la plupart des plateformes de séquençage Illumina).
- Calculez la profondeur de séquençage estimée en divisant le rendement attendu pour la cellule de flux Illumina par le nombre d’échantillons multiplexés dans la course.
REMARQUE : La détermination de la profondeur de séquençage souhaitée dépend de plusieurs facteurs, notamment le coût, le dispositif de séquençage utilisé et les objectifs de l’expérience. Deux exemples de profondeurs de séquençage différentes sont présentés ci-dessous dans la section Résultats représentatifs.
7. Analyse des données pour la détection des pathogènes
- Utilisez la plateforme web open source Chan Zuckerberg ID (CZID, https://www.czid.org, Illumina mNGS Pipeline v8.3) pour effectuer l’analyse métagénomique de détection des pathogènes. CZID inclut également des options pour télécharger des échantillons via une interface en ligne de commande et transférer directement d’un compte Illumina BaseSpace vers CZID. Un aperçu détaillé de l’analyse des données CZID mNGS est disponible sur https://chanzuckerberg.zendesk.com/hc/en-us/articles/13770737266196-Guide-to-mNGS-Data-Analysis
- Connectez-vous à CZID, cliquez sur « télécharger », sélectionnez « métagénomique » pour type d’analyse, puis sélectionnez les fichiers fastq d’entrée depuis la séquençage. Saisissez les informations nécessaires pour chaque échantillon, y compris le nom de l’échantillon, le type d’échantillon (LCR) et le nucléotide (ARN ou ADN).
REMARQUE : Les fichiers téléchargés sont traités automatisés via le pipeline CZID, qui intègre plusieurs outils de calcul tels que STAR, Bowtie2, Trimmomatic, PriceSeq et GSNAP, entre autres, afin de supprimer les lectures de faible qualité et les séquences de lecture humaines. Le pipeline aligne ensuite et assemble les séquences filtrées en utilisant les bases de données nucléotides (NT) et protéines (NR) du National Center for Biotechnology Information (NCBI) via les logiciels Minimap2 et DIAMOND, respectivement. La sortie finale du CZID inclut les décomptes de taxons NT/NR et les comptes contigus (segments contigus et superposés générés lors de l’assemblage) pour chaque échantillon.
- Une fois traité, cliquez sur le dossier projet contenant les échantillons et naviguez vers le tableau de bord résumé.
- Faites défiler vers le bas pour examiner le nombre de lectures par échantillon (et comparez au nombre attendu de lectures en fonction de la cellule de flux Illumina utilisée et du nombre d’échantillons multiplexés lors de la séquence), le pourcentage de lectures passant le filtrage de contrôle qualité, et le taux de compression en double afin d’évaluer une suramplification PCR biaisée.
- La présence de contamination par des fragments d’acides nucléiques environnementaux – soit acquis lors du processus de collecte d’échantillons de LCR, soit lors de l’extraction/préparation de la bibliothèque – est un phénomènecourant 20. Pour déterminer si les lectures bactériennes, fongiques et parasites sont de véritables positives, des méthodes potentielles utilisant les décomptes de lectures normalisées (lectures par million, ou rPM) peuvent être envisagées :
- Créer un modèle de fond (par exemple, en utilisant des échantillons d’eau extraits, préparés et séquencés) directement dans CZID pour filtrer tout taxon avec des rPM qui sont similaires ou inférieurs en abondance par rapport au modèle de fond (et donc susceptibles de représenter une contamination). Dans le dossier du projet, vérifiez chacun des échantillons de contrôle de l’eau puis cliquez sur « modèle d’arrière-plan » pour créer. Lors de l’analyse de chaque échantillon, sélectionnez le modèle de fond et cliquez sur « filtres de seuil » pour appliquer un score « NTZ » > 1 (ou supérieur) afin de filtrer les contaminants probables.
- Utilisez des seuils rPM prédéfinis, tels quel’échantillon rpM > 10* (contrôles rpm moyens), ou log10 (échantillon rPM) d’au moins 1 logarthum supérieur au log 10 (rPMmoyen de la cohorte entière). Tout d’abord, téléchargez le rPM spécifique au taxon pour chaque échantillon en vérifiant chaque échantillon sur la page du projet, en cliquant sur Télécharger et en sélectionnant « Résultats combinés des taxons d’échantillon ». Le fichier .csv résultant peut être analysé à l’aide de logiciels statistiques tels que R ou Stata selon le souhait.
REMARQUE : Compte tenu de la faible abondance relative des agents pathogènes viraux, tout virus ayant un potentiel neuroinvasif connu et au moins une lecture correspondant au génome viral doit être considéré comme un résultat positif (et idéalement, confirmé par des tests répétés de mNGS ou des tests cliniques spécifiques à chaque pathogène). De plus, le pipeline analytique exige que chaque lecture dans le fichier fastq soit alignée de manière computationnelle sur le génome de référence le plus probable dans la base de données NCBI. En conséquence, certaines lectures peuvent être relativement non spécifiques à un pathogène particulier et pourraient théoriquement s’aligner sur un large éventail d’organismes plutôt qu’avec le taxon spécifique assigné dans la CZID. Compte tenu de ce risque, la confirmation manuelle des appels positifs doit être effectuée en confirmant l’alignement via l’outil NCBI BLAST basé sur le web (les lectures spécifiques d’un taxon peuvent être envoyées directement du CZID au NCBI BLAST en cliquant sur l’icône blast à côté du taxon concerné).