Pour déterminer la performance de notre protocole, nous avons utilisé le 10x Genomics Chromium Single Cell 3' Gene Expression Kit v3.1 pour la préparation de la bibliothèque et analysé les données de snRNA-seq avec le paquet Seurat12,13.
La figure 6 montre les résultats d’une bibliothèque représentative de snRNA-seq. Pour évaluer la qualité de nos noyaux, nous avons tracé le nombre de gènes par rapport au nombre de transcrits (définis par des identifiants moléculaires uniques (UMI)) colorés par la fraction de lectures mitochondriales (Figure 6A). Les noyaux de bonne qualité présentent généralement un nombre plus élevé de lectures, corrélant le nombre d’UMI et de gènes, et de faibles fractions de lecture mitochondriales.
Pour l’analyse ultérieure, les noyaux avec moins de 500 ou plus de 4000 gènes comptés, ou plus de 5% d’ARN mitochondrial ont été exclus (n = 828). Seuls les gènes exprimés dans un minimum de trois noyaux ont été inclus. Nous avons détecté environ 20 000 gènes au total dans les 6 000 noyaux restants avec 1 600 gènes médians et 2 800 UMI médians par noyau (Figure 6B).
Le regroupement était basé sur des gènes très variables. Nous avons identifié un total de 18 grappes. Les identités cellulaires ont été annotées en fonction de gènes marqueurs connus (non représentés). Les sous-grappes d’un type cellulaire ont été résumées en une seule grappe, ce qui donne un total de 11 types de cellules distinctes : podocytes, tubule proximal (PT), membre mince (tL), membre ascendant épais (TAL), tubule convoluté distal (DCT), tubule de connexion (CNT), principal du canal collecteur et cellules intercalées (CD-PC, A-IC, B-IC), épithélium médullaire profond du bassin (DMEP) et endothélium. Les profils d’expression génique des marqueurs enrichis en grappes ont été visualisés dans un diagramme à points (figure 6C) et les groupes de types cellulaires dans un diagramme d’incorporation de voisins stochastiques distribués en t (t-SNE) (figure 6D).
Pour évaluer les distributions des types cellulaires dans notre échantillon, le pourcentage de chaque type de cellule a été calculé (figure 6E) et utilisé pour déterminer le rapport du PT au TAL. Le PT est principalement situé dans le cortex rénal et souvent surreprésenté dans les ensembles de données sur les cellules uniques du rein, car les cellules du PT sont faciles à dissocier et très abondantes dans les échantillons de reins entiers. Le TAL, quant à lui, s’étend sur toute la moelle extérieure14. Ainsi, le rapport des fractions PT et TAL représente une bonne mesure de l’enrichissement des types de cellules médullaires dans un ensemble de données unicellulaires rénales. En général, le rapport PT/TAL dans les ensembles de données sur les reins entiers unicellulaires variait de 8 (données non publiées provenant de tissus rénaux entiers traités à la protéase froide) à 45 pour les tissus dissociés enzymatiquement10,14,15. Dans l’ensemble de données snRNA-seq présenté ici, nous avons pu atteindre un rapport PT/TAL de 2. Ce résultat illustre que l’élimination de l’excès de cortex pendant la dissection tissulaire combinée à un séquençage de l’ARNn entraîne une amélioration frappante de la représentation du type de cellules rénales.

Figure 1 : Vue d’ensemble schématique du flux de travail. Le protocole comprend quatre étapes principales qui comprennent la dissection tissulaire suivie de l’isolement des noyaux, du tri des noyaux et d’une évaluation finale de la pureté et de la concentration. Barre d’échelle = 100 μm. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 2 : Dissection rénale et préparation tissulaire. (A) Image représentative d’un rein entier disséqué. Les lignes pointillées indiquent les coupes nécessaires pour obtenir une tranche médiane de 1-2 mm avec une représentation de tous les types de cellules rénales. (B) Image représentative de la tranche du milieu obtenue. Les lignes pointillées indiquent les coupes pour la coupe du cortex sur le côté. (C) Image représentative d’un morceau de rein central avec cortex coupé. Le cortex (C), la moelle externe (OM) et la moelle interne (IM) sont clairement visibles. Barre d’échelle = 500 μm. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 3 : Homogénéisation tissulaire et purification des noyaux. (A) Image représentative montrant du tissu rénal suffisamment haché. Barre d’échelle = 500 μm. (B) Homogéner après la première étape d’homogénéisation (25 coups avec le pilon A, tube de meulage de 2 mL). (C) Homogéner après deuxième étape d’homogénéisation (15 coups avec le pilon B, tube de meulage de 2 mL). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 4 : Stratégie de contrôle pour le tri des noyaux. (A) Une première porte P1 a été définie sur la base de la diffusion vers l’avant (FSC) par rapport à la diffusion latérale (SSC) pour exclure les débris et les agrégats. (B) Une porte ultérieure basée sur la zone DAPI (DAPI-A) par rapport à la hauteur DAPI (DAPI-H) excluait les noyaux et les multiplets vides ou endommagés. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 5 : Suspension de noyaux avant et après tri des noyaux. (A, C) Noyaux colorés DAPI (bleu). (B, D) Superposition de DAPI et de canal en fond clair (BF). Avant le tri (panneau supérieur), la suspension de noyaux contient des débris cellulaires et des agrégats (marqués avec des pointes de flèches blanches). Après tri (panneau inférieur), la suspension des noyaux semble beaucoup plus propre. Des exemples de noyaux colorés DAPI sont étiquetés avec des pointes de flèches noires. Les noyaux de bonne qualité apparaissent ronds et lisses avec une membrane intacte et sont bien séparés, tandis que les noyaux de mauvaise qualité apparaissent ridés et présentent une perte de la membrane nucléaire. Barre d’échelle = 250 μm. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 6 : Contrôle de la qualité et analyse d’un ensemble de données représentatif de l’ARNn-seq. (A) Nombre de gènes (nGene) tracé en fonction du nombre d’identificateurs moléculaires uniques (nUMI) colorés par la fraction de lectures mitochondriales (percent.mt). Les noyaux de faible qualité correspondent au quadrant inférieur gauche du graphique (n = 828) et ont été exclus de l’analyse ultérieure. (B) Distribution et médiane du gène nGene et du nUMI détectés par noyau dans l’ensemble de données snRNA-seq, représentant 6 177 noyaux (> 500 gènes). Les bibliothèques ont été séquencées à une profondeur médiane de ~ 8 200 lectures mappées par noyau. (C) Diagramme à points montrant les profils d’expression génique des marqueurs enrichis en grappes (axe des x) pour les types de cellules individuels (axe des y). La taille du point correspond à la proportion de cellules exprimant le gène indiqué. La couleur correspond à l’expression moyenne. (D) Tracé d’incorporation stochastique de voisins distribués en T (t-SNE) de types cellulaires identifiés. (E) Distribution des types cellulaires dans l’ensemble de données snRNA-seq. PT, tubule proximal; tL, membre mince; TAL, membre ascendant épais, DCT, tubule convolué distal; CNT, tubule de raccordement; CD-PC, cellules principales de collecte des conduits; A-IC, cellules intercalées de type A; B-IC, cellules intercalées de type B; DMEP, épithélium médullaire profond du bassin. Veuillez cliquer ici pour voir une version agrandie de cette figure.
| Réactif | Concentration finale | Volume (mL) |
| a) 4 % BSA / PBS | | |
| Solution saline tamponnée au phosphate (PBS) contenant 10 % d’albumine bovine | 4% | 2 |
| PBS (solution saline tamponnée au phosphate) 1X sans calcium ni magnésium | - | 3 |
| b) 0,04 % BSA / PBS | | |
| 4 % BSA / PBS | 0.04 % | 0.5 |
| PBS (solution saline tamponnée au phosphate) 1X sans calcium ni magnésium | - | 49.5 |
| (C) Tampon de lyse des noyaux 1 (NLB1) | | |
| Tampon de lyse nucléaire EZ | - | 4 |
| Inhibiteur de la RNase RiboLock (40 U/μL) | 1 U/μL | 0.1 |
| Complexe ribonucléoside-vanadyle (200 mM) | 10 mM | 0.2 |
| d) Tampon de lyse des noyaux 2 (NLB2) | | |
| Tampon de lyse nucléaire EZ | - | 4 |
| Inhibiteur de la RNase RiboLock (40 U/μL) | 1 U/μL | 0.1 |
| e) Tampon de suspension de noyaux (NSB) | | |
| 0,04 % BSA / PBS | - | 2 |
| Inhibiteur de la RNase RiboLock (40 U/μL) | 1 U/μL | 0.05 |
| f) Tampon de gradient de saccharose (10 % saccharose) | | |
| Poids 1 g de saccharose | | |
| Dissoudre dans 6 ml de tampon de lyse EZ nucléaire | | |
| Remplir jusqu’à 10 mL avec le tampon de lyse EZ nucléaire | | |
| Filtrer à travers un filtre à seringue de 0,2 μm dans un tube neuf | | |
Tableau 1 : Recettes de solution : (A) Préparation de 4 % de BSA/1x PBS. Filtrer à l’aide d’un filtre à seringue à membrane SFCA de 0,2 μm et conserver sur la glace jusqu’à utilisation. (B) Préparation de 0,04% BSA/1 x PBS. Filtrer à l’aide d’un filtre à seringue à membrane SFCA de 0,2 μm et conserver sur la glace jusqu’à utilisation. c) Préparation du tampon de lyse des noyaux 1 (NLB1). Les volumes indiqués sont fournis par échantillon. Conserver sur la glace jusqu’à utilisation. d) Préparation du tampon de lyse des noyaux 2 (NLB2). Les volumes indiqués sont fournis par échantillon. Ajouter l’inhibiteur de la RNase RiboLock à NLB2 directement avant utilisation, comme indiqué dans le protocole. Conserver sur la glace jusqu’à utilisation. e) Préparation du tampon de suspension de noyaux (NSB). Les volumes indiqués sont fournis par échantillon. Ajouter l’inhibiteur de la RNase RiboLock au NSB directement avant utilisation, comme indiqué dans le protocole. Conserver sur la glace jusqu’à utilisation. f) Préparation du tampon de gradient de saccharose. Filtrer à l’aide d’un filtre à seringue à membrane SFCA de 0,2 μm et conserver sur la glace jusqu’à utilisation.