$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Het repertoire van open source pakketten ontworpen om te analyseren scRNA-Seq datasets dramatisch gestegen40 met de meerderheid van deze pakketten gebruik R gebaseerde talen3. Hier, representatieve resultaten met behulp van twee van deze pakketten worden gepresenteerd: beoordeling van ongecontroleerde groepering van enkele cellen op basis van genexpressie en bestellen van afzonderlijke cellen langs een traject om te lossen cel heterogeniteit en deconstrueren biologische processen.
Figuur 4 illustreert het gebruik van Seurat voor voorbehandeling kwaliteitscontroles en downstream bioinformatics analyse. Ten eerste, filtratie en verwijdering van afwijkende cellen uit de analyse is noodzakelijk om de controle van de kwaliteit. Dit werd gedaan met behulp van viool (figuur 4a) en scatter percelen (figuur 4b) te visualiseren van het percentage van de mitochondriale genen, aantal genen (nGene) en aantal UMI (nUMI) om cel paren en uitschieters te identificeren. Een cel met een duidelijke uitbijter aantal genen, UMI of percentage van mitochondriale genen is verwijderd met behulp van Seurat de FilterCells functie. Aangezien Seurat gebruikt belangrijkste component (PC) analyse scores aan clusters cellen, bepalen van statistisch significante PC's op te nemen is een cruciale stap. Elleboog percelen (Figuur 4 c) werden gebruikt voor de selectie van de PC, in welke PC's buiten het plateau van de "standaarddeviatie van PC' as werden uitgesloten. De resolutie van clustering werd ook gemanipuleerd waaruit blijkt dat het aantal clusters kan worden gewijzigd, variërend van 0.4 (lage resolutie leidt tot minder cel clusters, Figuur 4 d) tot en met 4 (hoge resolutie leidt tot hogere cel clusters, figuur 4e ). Met een lage resolutie is het waarschijnlijk dat elk cluster een bepaald celtype vertegenwoordigt, overwegende dat bij hoge resolutie dit ook subtypen of tijdelijke Staten van een celpopulatie betekenen kan. In dit geval werden met een lage resolutie cluster instellingen gebruikt voor het verder analyseren expressie heatmaps (met behulp van Seurat DoHeatmap functie) om de hoogst uitgedrukte genen in een bepaald cluster (figuur 4f) te identificeren. In dit geval werden de hoogst uitgedrukte genen geïdentificeerd door het beoordelen van differentiële expressie in een bepaald cluster ten opzichte van alle andere clusters gecombineerd, waaruit blijkt dat elk cluster een unieke werd vertegenwoordigd door gedefinieerde genen. Bovendien kunnen individuele kandidaat-genen worden gevisualiseerd op percelen van de tSNE met behulp van Seurat de FeaturePlot functie (Figuur 4 g). Dit toegestaan voor ontcijferen of er clusters die macrofagen vertegenwoordigd waren. Met behulp van FeaturePlot, vonden we dat beide cluster 2 en 4 waren Cd68 - een marker van de pan-macrofaag uitdrukken.
De Monocle-pakket werd gebruikt voor ondersteunend cel clusters geïdentificeerd in Seurat, en voor het bouwen van cel trajecten, of door pseudotemporal te bestellen, te recapituleren biologische processen (Figuur 5). Bestellen van de pseudotemporal kan worden gebruikt voor monsters waar eencellige expressieprofielen verwachting te volgen van een biologische tijdsverloop. Cellen kunnen besteld worden langs een continuüm van pseudotemporal om op te lossen tussentijdse staten, bifurcatie punten van twee alternatieve cel lot, en gene handtekeningen ten grondslag liggen aan de verwerving van elk lot te identificeren. Ten eerste, vergelijkbaar met de Seurat filtratie, kwalitatief slechte cellen werden verwijderd, zodat de verdeling van mRNA over alle cellen log normaal en viel tussen de bovenste en onderste grenzen was zoals aangegeven in figuur 5a. Vervolgens met de Monocle newCellTypeHierarchy functie, afzonderlijke cellen werden ingedeeld en geteld met behulp van bekende afstamming markers (Figuur 5b, 5 c). Bijvoorbeeld, werden cellen, uiting geven aan de PDGF receptor Alfa- of Fibroblast specifieke eiwitten 1 toegewezen aan de cel Type #1 maken een criterium voor het definiëren van fibroblasten. Vervolgens werd deze populatie (cel Type #1) beoordeeld om te ontcijferen fibroblast trajecten. Om dit te doen, werd de Monocle differentiële GeneTest functie gebruikt, die ten opzichte van de cellen die de extreme lidstaten binnen de bevolking vertegenwoordigt en differentiële genen gevonden voor het bestellen van de resterende cellen onder de bevolking (Figuur 5 d). Spruitstuk leermethoden (een type van vermindering van de niet-lineaire dimensionaliteit) toe te passen op alle cellen, was een coördinaat langs een pseudotemporal pad toegewezen. Dit traject werd vervolgens gevisualiseerd door cel staat (figuur 5e) en pseudotime (figuur 5f).

Figuur 1: stroomschema. Stappen van hele dieren voorbereiding naar het analyseren van één cel RNA-Seq datasets tot het indienen van de definitieve datasets in een openbaar archief. Gel kralen in emulsie (edelstenen) verwijzen naar kralen met barcoded oligonucleotides, die weer duizenden afzonderlijke cellen inkapselen. Klik hier voor een grotere versie van dit cijfer.

Figuur 2: creëren van levensvatbare interne celsuspensie van zenuwweefsel. (a) cartoon overzicht van kwaliteitscontroles. (b) cellen en puin met cellen nog steeds opgenomen in puin (rode pijlen). (c) cellen vrijgelaten uit puin (rode pijlen). (d) cel isolatie door FACS. P0: puin breuk; P1: cel-achtige breuk; P3: uitsluiting van duplets; P4: levensvatbaarheid kleurstof (Sytox Orange) negatieve breuk. (e) geen levensvatbaarheid kleurstof controle. (f) afbeelding P0 fractie vertegenwoordigen geïsoleerd puin. (g) afbeelding P4 fractie vertegenwoordigen geïsoleerd levensvatbare cellen (rode pijlen). (b) (c) (f) en (g) had nucleaire kleurstof toegevoegd 20 minuten voor de beeldvorming. Schaal Bars: 80 µm. Klik hier voor een grotere versie van dit cijfer.

Figuur 3: ondiepe sequencing voorspelt het aantal herstelde cellen in 10 X verwerkte monsters. (a) een voorbeeld (monster 1.6) voor MiSeq-gegenereerd csv aanbieding cel barcodes en de bijbehorende UMI wordt geteld zoals bepaald door vertrouwen toegewezen leest. (b) Barcode rangschikking perceel voor monster 1.6 toont een aanzienlijke daling in UMI tellen als een functie van de cel barcodes. De onderbroken en ononderbroken lijnen geven de cutoff tussen cellen en achtergrond zoals bepaald door de visuele inspectie. (c) cel barcodes waargenomen met behulp van de Cell Ranger pijpleiding post-HiSeq onthult ondiepe sequencing benaderd nauwkeurig het aantal cellen voor monster 1.6. (d) een voorbeeld van een set-up van de stroom-cel op basis van ondiepe sequencing afgeleide ramingen van de cel. Monster 1.6, aangezien ondiepe sequencing voorspeld 3480 cellen, 1,17 rijstroken waren toegewezen om > 100.000 leest per cel sequencing dekking in HiSeq. Opmerking: Alle rijstroken moeten toevoegen aan 100%. Klik hier voor een grotere versie van dit cijfer.

Figuur 4: kwaliteitscontrole en bio-informatica van eencellige RNA-Seq dataset met behulp van Seurat R package. (a) percelen van kwaliteitscontrole metrics, waaronder een aantal genen, aantal unieke moleculaire identificatoren (UMIs) en het percentage van afschriften toewijzen aan het mitochondriaal genoom. (b) monster gene percelen detecting cellen met afwijkende niveaus van mitochondriale afschriften en UMIs. (c) monster elleboog plot gebruikt voor ad hoc bepaling van statistisch significante PC's. De onderbroken en stip onderbroken lijnen geven de cutoff wanneer een duidelijke "elleboog" zich manifesteert in de grafiek. PC afmetingen vóór deze elleboog zijn opgenomen in de downstream-analyse. (d, e) Grafiek gebaseerde cel clusters gevisualiseerd met twee verschillende resoluties in een laag-dimensionale ruimte met behulp van een perceel van tSNE. (f) de bovenste markering genen (geel) voor elk cluster gevisualiseerd op een expressie heatmap met behulp van Seurat de DoHeatmap functie. (g) visualiseren marker expressie van, bijvoorbeeld, Cd68 gen vertegenwoordigen macrofagen (paars) met behulp van Seurat de FeaturePlot functie. Dit suggereert dat cluster 2 en 4 (in deelvenster d) van deze dataset vertegenwoordigt macrofagen. Klik hier voor een grotere versie van dit cijfer.

Figuur 5: cel categorisatie en bestellen langs peudotemporal traject met Monocle toolkit. (a) de verdeling van mRNA (afgeleid van de graven van de UMI) over alle cellen in een monster van de inspectie. Alleen cellen met mRNA tussen 0 - ~ 20.000 werden gebruikt voor downstream-analyse. (b, c) Toewijzen en tellen celtypes gebaseerd op bekende afstamming cel markeringen. Bijvoorbeeld, cellen, uiting geven aan de PDGF receptor Alfa- of Fibroblast specifieke eiwitten 1 werden toegewezen aan de cel Type #1 pan-fibroblasten vertegenwoordigen Monocle de newCellTypeHierarchy functie. Aantal verschillende soorten cellen kan worden gevisualiseerd, een cirkeldiagram (b) en als een tabel (c). (d) het gebruikmaken van cel Type #1 (fibroblasten) als voorbeeld, de genen voor het bestellen van cellen kunnen worden gevisualiseerd met behulp van een scatterplot waaruit gene dispersie vs. gemiddelde expressie gebruikt. De rode curve toont de cutoff voor genen gebruikt voor het bestellen van berekend door het gemiddelde-variantie-model met behulp van de Monocle estimateDispersions functie. Genen die voldoen aan deze cutoff werden gebruikt voor downstream pseudotime bestellen. (e, f) Visualisatie van cel trajecten in een verminderde twee-dimensionale ruimte gekleurd door de cel "staat" (e), en de Monocle-toegewezen "Pseudotime" (f). Klik hier voor een grotere versie van dit cijfer.