Summary

Detectie van Copy Number Wijzigingen met behulp van Single Cell Sequencing

Published: February 17, 2017
doi:

Summary

Enkele cel sequencing is een steeds populairder en toegankelijk instrument voor de aanpak van genomische veranderingen op hoge resolutie. Wij bieden een protocol dat enkele cel sequencing gebruikt om het aantal kopieën veranderingen te identificeren in enkele cellen.

Abstract

Detectie van genomische veranderingen op één enkele cel resolutie is belangrijk voor het karakteriseren van genetische heterogeniteit en evolutie in normale weefsels, kanker, en microbiële populaties. Traditionele methoden voor de beoordeling van de genetische heterogeniteit zijn beperkt door lage resolutie, lage gevoeligheid en / of lage specificiteit. Enkele cel sequencing is uitgegroeid tot een krachtig hulpmiddel voor het detecteren van genetische heterogeniteit met hoge resolutie, hoge gevoeligheid en, indien passend geanalyseerd hoge specificiteit. Hier bieden we een protocol voor de isolatie, hele genoom amplificatie, sequencing en analyse van enkele cellen. Onze aanpak maakt het mogelijk om betrouwbare identificatie van megabase-scale copy number varianten in enkele cellen. Echter aspecten van dit protocol kan ook worden toegepast op andere typen genetische veranderingen in individuele cellen te onderzoeken.

Introduction

Veranderingen in DNA copy number kan variëren in grootte van een aantal basenparen (copy aantal varianten) (CNV) en om hele chromosomen (aneuploïdie). Kopie aantal veranderingen beïnvloeden grote gebieden van het genoom kan belangrijke fenotypische gevolgen hebben door het veranderen van de expressie tot duizenden genen 1, 2. CNV aanwezig in alle cellen van een populatie kan worden gedetecteerd door sequentiebepaling of massa-microarray gebaseerde werkwijzen 3, 4 zijn. Echter kunnen ook populaties genetisch heterogeen zijn, waarbij bestaande CNV in een subset van de populatie of zelfs individuele cellen. Genetische heterogeniteit is gebruikelijk bij kanker, tumor evolutie drijven, en ook in normale weefsels, met onbekende gevolg 5, 6, 7, 8, 9 <sup>, 10.

Traditioneel werd genetische heterogeniteit beoordeeld worden door cytologisch benaderingen of bulk sequencing. Cytologisch benaderingen, zoals fluorescentie in situ hybridisatie (FISH), chromosoom spreads, en spectrale karyotypering (SKY), hebben het voordeel van het identificeren van veranderingen aanwezig zijn in individuele cellen, maar hebben een hoge foutenpercentages als gevolg van artefacten van hybridisatie en verspreiden van 11. Deze benaderingen zijn ook beperkt in hun resolutie alleen het onthullen van duplicaties over meerdere megabasen. Sequencing of microarrays van bulk DNA, hoewel hoger in nauwkeurigheid en resolutie, is minder gevoelig. Om genetische heterogeniteit te detecteren door populatie gebaseerde benaderingen moeten de varianten aanwezig in een aanzienlijke fractie van cellen in de populatie. De opkomst van methoden om de genomische DNA van enkelvoudige cellen amplificeren heeft het mogelijk gemaakt om het genoom van afzonderlijke cellen sequentie. Enkele cel Sequéncing heeft de voordelen van hoge resolutie, een hoge gevoeligheid, en, indien van toepassing voor kwaliteitscontrole methoden worden toegepast, hoge nauwkeurigheid 12.

Hier beschrijven we een werkwijze voor het detecteren megabase schaal kopieaantal veranderingen in individuele cellen. We isoleren enkele cellen door microaspiration, genomisch DNA te amplificeren met behulp linker-adapter PCR bibliotheken bereiden voor de volgende generatie sequencing en detecteren kopieaantal varianten zowel verborgen Markov model en circulaire binaire segmentatie.

Protocol

1. Het isoleren van enkele cellen Bereid de microaspirator Verwijder de doorzichtige plastic uiteinde van de aspirator buisstelsel en steek het smalle uiteinde in een uiteinde van een 1 cm lange PVC buis met 3/16 "binnendiameter. Plaats de uitlaat van een 0,2 urn spuitfilter in het andere uiteinde van de 1 cm lange PVC buis met 3/6 "binnendiameter. Steek de inlaat van de 0,2 urn spuitfilter in een einde van een 6 inch-lange PVC buis met 5/16 "binnendiameter. Optioneel: Snijd de 6 inch-lange PVC buis met 5/16 "binnendiameter in de helft en plaats een in-line water val tussen de twee helften. Breek een 5 ml plastic serologische pipet op 1 ml graduatie en plaats het gebroken uiteinde in het open uiteinde van de PVC buis met 5/16 "binnendiameter. Plaats de uitlaat van het 5 mL plastic serologische pipet in het flexibele uiteinde van een aspirator eenheid (waarbij de doorzichtige plastic eindewas verwijderd). Voor de opslag, hebben betrekking op de rode spreekbuis van de aspirator buis montage met een steriele plastic buis. Trekken uit een glazen capillair met een inwendige diameter van 10-30 urn doordat zij het midden van de buis bij een Bunsen brander vlam en toepassen spanning op beide uiteinden van de buis. Breek de buis getrokken in het midden om twee potentiële aspirator naalden te creëren. Herhaal deze stap met een aantal buizen, aangezien slechts enkele buizen zal eindigen met een inwendige diameter die geschikt is voor het plukken van enkele cellen is. Voor de opslag, plaats twee stroken van boetseerklei in een 15 cm petrischaal en zet de aspirator naalden over de strips. Pick-cellen Bereid een enkele cel schorsing afhankelijk van het type cel en experiment. Bijvoorbeeld om hechtende cellen te bereiden zoals humane fibroblast cellijnen oogst cellen door behandeling met trypsine en overdracht cellen om een ​​conische buis met geschikte media. Voor, tijdens offter de voorbereiding van de enkele cel suspensie (afhankelijk van hoe lang het duurt om de enkele celsuspensie te bereiden), het instellen van de kap voor de hele genoom amplificatie. Spray beneden het oppervlak van de kap, pipetpunt dozen, en de pipet tips met 10% bleekwater en veeg met een papieren handdoek. Herhaal deze stap met 70% ethanol. Voeg 8 pi van water uit het hele genoom amplificatie (WGA) kit afzonderlijke putjes van een 96-well PCR plaat, één putje voor elke cel waarvan de sequentie. Bedek de 96-well PCR plaat met een deksel van een 96-well weefselkweekplaat en plaats op ijs. Voeg 1000 cellen aan 10 ml media of fosfaatgebufferde zoutoplossing (PBS) in een 15 cm petrischaal en plaats de schotel op ijs te voorkomen dat de cellen hechten aan de schaal. Breng de cellen in een petrischaal en de 96-well PCR plaat met deksel op ijs om een ​​licht microscoop met 10x objectief. Verhoog de opening van de aspirator naald door voorzichtig de uitgesponnen eind van Th tikkene aspirator naald op een hard oppervlak, zodanig dat de punt afbreekt. Steek het brede uiteinde van de aspirator naald in de heldere einde van de microaspirator. Plaats de petrischaal met cellen op de microscoop podium. Plaats de rode mondstuk van de aspirator in de mond. Met één hand de naald aspirator en anderzijds naar de petrischaal met cellen bewegen. Identificeer enkele cellen waarvan de sequentie. Met behulp van mond zuigkracht, trek een enkele cel in de aspirator naald samen met ~ 1-2 pi van media of PBS. Breng de cel in de 8 pi water in een enkel putje van de 96-well PCR plaat. Voorkomen dat er luchtbellen bij de overdracht van de cel. Herhaal stap 1.2.7 tot het gewenste aantal cellen zijn geïsoleerd. Houd de PCR-plaat op ijs en bedekt met deksel in tussen plukken cellen. Markeren de putten die cellen hebben ontvangen. Terwijl het oppakken enkele cellen, ontdooien de 10X enkele cel lysis en fragmentatie buffer uit de hele genome amplificatie kit. Na het kiezen van de gewenste aantal cellen, gaat onmiddellijk op hele genoom amplificatie. 2. Whole Genome Amplification Om besmetting tijdens de gehele genoom amplificatie te voorkomen, voeg alle reagentia in een weefselkweek kap, gebruiken pipet tips met filters, en verander de pipet tip in tussen de putten. Bereid een werkende lysis en fragmentatie bufferoplossing uit het hele genoom amplificatie (WGA) kit. Voor elke set van maximaal 32 cellen combineren 32 ui 10x cel lysis en fragmentatie buffer en 2 pl proteinase K oplossing in een microcentrifugebuis. Vortex de buis om de oplossingen te mengen. Voeg 1 pl werken lysis en fragmentatie bufferoplossing aan elk putje en pipet op en neer te mengen. Bedek de plaat en sluit alle putjes met een plastic folie. Kort centrifuge de plaat in een mini plaat spinner. Thermische cyclus als follows: 50 ° C, 1 h en 99 ° C, 4 min. Koel de plaat op ijs en kort centrifuge de plaat in een mini-plaat spinner. Bereid een werkende bibliotheek voorbereiding bufferoplossing uit het hele genoom amplificatie kit. Voor elke cel, combineer 2 pl 1x cel bibliotheek bereiding buffer en 1 pi van bibliotheek stabilisatie oplossing in een microcentrifugebuis. Bereid de werkoplossing meerdere cellen in dezelfde microcentrifugebuis. Verwijder de plastic film en voeg 3 pl werkbibliotheek preparaat bufferoplossing aan elk putje. Pipetteer de inhoud van de put op en neer om te mengen. Vervang de kunststoffolie. OPMERKING: De kunststoffilm kan worden hergebruikt gedurende het gehele genoom amplificatieproces totdat de putjes gaan boringen in de film. Wanneer dit gebeurt, overschakelen naar een nieuwe plastic folie. Kort centrifuge de plaat in een mini spinner plaat en incubeer bij 95 ° C gedurende 2 minuten. Koel de plaat op ijs en in het kortcentrifuge de plaat in een mini plaat spinner. Verwijder de plastic folie, voeg 1 pl bibliotheek enzympreparaat aan elk putje en pipet de inhoud van de put op en neer te mengen. Houd de bibliotheek voorbereiding enzym op ijs of in een koude blok in deze stap. Vervang de kunststoffolie. Kort centrifuge de plaat in mini plaat spinner en thermische cyclus als volgt: 16 ° C, 20 min, 24 ° C, 20 min, 37 ° C, 20 min, 75 ° C, 5 min, en houden op 4 ° C. Koel de plaat op ijs en kort centrifuge de plaat in een mini-plaat spinner. Bereid een werkende versterking mix van het hele genoom amplificatie kit. Voor elke cel, combineer 48,5 pi water, 7.5 pi 10x versterking master mix, en 5 pi WGA DNA polymerase in een microcentrifugebuis. De werkoplossing mix meerdere cellen in dezelfde microcentrifugebuis. Houd de werkende mix op ijs. Verwijder de plastic folie, voeg 61 ul werken versterking mixaan elk putje, en de inhoud van de pipet goed op en neer om te mengen. Houd de werkende versterking mix op ijs of in een koude blok in deze stap. Vervang de kunststoffolie. Kort centrifuge de plaat in een mini plaat spinner en thermische cyclus als volgt: 95 ° C, 3 min, 25 cycli van 94 ° C, 30 sec en 65 ° C, 5 min, en houden op 4 ° C. Overdracht 60 pi van elk monster een aparte microcentrifugebuis en opslag bij – 20 ° C voor gebruik in stap 3. Voeg DNA laadkleurstof het resterende volume van elk monster (dat wil zeggen 3 pi van 6x DNA laadkleurstof tot 15 pl monster) en in werking 5 pi van elke reactie op een 1% agarosegel. OPMERKING: Cellen die met succes versterkt zal verschijnen als een uitstrijkje van 250 tot 1000 bp. Monsters die geen uitstrijkje hebben laten zien of alleen laten een zwakke uitstrijkje is onwaarschijnlijk dat bruikbare sequencing data opleveren. 3. Sequencing Zuiver monsters met paramagnetische kralen. Thaw DNA-monsters op ijs. Vortex paramagnetische kralen tot de oplossing homogeen en incubeer de korrels bij kamertemperatuur gedurende ten minste 30 min. Transfer 20 pi van elk DNA monster in een schone microcentrifugebuis. De rest van het monster kan worden opgeslagen bij -20 ° C. Voeg 30 pi (1,5x) van paramagnetische kralen aan elk DNA-monster en vortex te mengen. Incubeer de monsters bij kamertemperatuur gedurende 10 min. Laat de voorraadoplossing van kralen bij kamertemperatuur voor gebruik in stap 3,4. Plaats de buizen op een magnetische strip buis gedurende 2 minuten of totdat de korrels vormen een pellet en het supernatant helder. Met behulp van een pipet P200, verwijder zoveel van de supernatant mogelijk zonder de kralen. Voeg 180 ui 80% ethanol om het DNA-korrel mengsel. meermaals opzichte van de magneet "spoelen" de kralen door de ethanoloplossing draaien de buizen. Met behulp van een pipet P200, verwijder zoveel mogelijk van de ethanol wassen als mogeble zonder verstoring van de kralen. Herhaal 3.1.6 en 3.1.7. De lucht drogen de kralen gedurende ongeveer 10 minuten of tot er geen ethanol meer zichtbaar is. Naar de volgende stap bij de parels verschijnen gebarsten of de wand van de buis vallen. Verwijder de monsters van de magneetstrip. Voeg 40 ul van 10 mM Tris pH 8,0 met het DNA van de bolletjes elueren en vortex de monsters om de kralen resuspenderen. Incubeer gedurende 2 minuten bij kamertemperatuur. Kort Centrifugeer de monsters om de vloeistof op de bodem van de reageerbuis te verzamelen en plaats de buizen op een magnetische strip buis gedurende ten minste twee minuten. Breng de eluent in schone microcentrifugebuizen zonder verstoring van de kralen. sample normalisatie Kwantificeren van de concentratie van elk monster onder toepassing van een spectrofotometer. De concentraties moet liggen van 10 tot 30 ng / ul. Verdun elk monster tot 0,2 ng / ul gebruikt 10 mM Tris pH 8,0. De volgende stappen zullen een minimum ik vereisennput van 5 pi van deze verdunning. voorbereiding bibliotheek Bereid sequencing bibliotheken door de instructies van de bibliotheek voorbereiding kit 13. eindschoonmaak Reinig de monsters volgens stap 3,1. Gelezen lengte van 50 bp, 75 bp of 150 bp, gebruiken een verhouding van kraag naar volume 1,5x, 1x of 0.6x monster, respectievelijk. Elueer met 15 ui 10 mM Tris pH 8,0. Run de monsters op een fragment analysator om de grootteverdeling van de bibliotheek 14 controleren. De grootteverdeling wordt gelijkmatig verdeeld 150-900 bp. Kwantificeren monsters met behulp van qPCR Met behulp van een bibliotheek kwantificering kit, opzetten van een qPCR reactie volgens de kit instructies 15. Laat een kolom leeg om positieve normen (DNA normen uit de kit) en negatieve normen (water of andere blanco-oplossing) toe te voegen. Thermische cyclus als follows: 95 ° C, 5 min (toenamesnelheid van 4,8 ° C / sec) en 35 cycli van 95 ° C, 30 sec (toenamesnelheid van 4,8 ° C / sec) en 60 ° C, 45 sec (toenamesnelheid van 2,5 ° C / sec). pooling Bepalen hoeveel rijstroken op de doorstroomcel nodig zijn en verdeel de monsters in groepen, waarbij één groep per rijstrook. Voor iedere groep, kiest het monster met de laagste concentratie op basis van de qPCR gegevens uit stap 3.5.2 en normaliseren alle monsters in die groep die concentratie waarbij 10 mM Tris pH 8,0. de monsters in elke groep samen te bundelen. sequencing Load zwembaden op next-generation sequencing machine volgens standaardprocedures 16. 4. Data Analysis OPMERKING: Een Unix-gebaseerde omgeving is vereist om de programma's en scripts in dit gedeelte werking. Installeer de in het protocol genoemde software na hun in stallatie gidsen. Alle scripts zijn te vinden op https://sourceforge.net/projects/singlecellseqcnv/. Knip de leest 40-nt met behulp van fastx_trimmer van FastX-Toolkit versie 0.0.13 17. fastx_trimmer -Q33 -i example.fastq l 40 -o example.trim.fastq Lijn de leest om de juiste verwijzing genoom (MM9 voor muis, hg19 voor menselijke) met behulp van BWA versie 0.6.1 met standaardopties 18. BWA ALN mm9.fa example.trim.fastq> example.sai BWA Samse mm9.fa example.sai example.trim.fastq> example.sam Verwijder Alignmenten Chrm en willekeurige chromosomen, dan sorteren en de index van de resulterende BAM-bestand met SAMTools versie 0.1.19 19. grep -v -w Chrm example.sam | grep -v willekeurige> example.filtered.sam samtools bekijken -uSh example.filtered.sam | samtools sort – example.filtered samtools index example.filtered.bam Ren HMMcopy om CNVs detecteren= "Xref"> 20 Gebruik gcCounter in HMMcopy een GC percentage referentiebestand voor het genoom te genereren. Gebruik de optie '-w 500000 "om het venster formaat op te geven. Gebruik dezelfde versie van het FASTA referentiebestand gebruikt in stap 4.2, maar zorg ervoor dat Chrm en willekeurige chromosomen worden verwijderd uit de fasta bestand. gcCounter -w 500000 mm9.fa> mm9_gc.wig Gebruik generateMap.pl in HMMcopy een wiebelen bestand voor mappability genereren. Gebruik de optie "-w 40" om aan te geven te lezen lengte. Gebruik dezelfde FASTA referentiebestand gebruikt in stap 4.4.1. generateMap.pl -b mm9.fa generateMap.pl -w 40-i mm9.fa mm9.fa -o mm9.bigwig Gebruik mapCounter in HMMcopy een mappability referentiebestand voor het genoom te genereren. Gebruik de optie '-w 500000 "om het venster formaat op te geven. mapCounter -w 500000 mm9.bigwig> mm9_map.wig Gebruik readCounter in HMMcopy een wiebelen bestand voor elk BAM bestand te genereren. readCounter -w 500000 example.filtered.bam> input.wig </li> Wijzig de paden naar de referentie-bestanden in de R scripts (run_hmmcopy.mm9.r of run_hmmcopy.hg19.r), gebruikt u de bovenstaande bestanden gegenereerd in 4.4.1 (bijv mm9_gc.wig) en 4.4.3 (bijv mm9_map. pruik) voor de variabelen GFile en mfile, die verwijzen naar de GC percentage referentiebestand en mappability referentiebestand, respectievelijk. Dan start de meegeleverde R script "run_hmmcopy.mm9.r" of "run_hmmcopy.hg19.r". R CMD BATCH run_hmmcopy.mm9.r of R CMD BATCH run_hmmcopy.hg19.r Voor de batchverwerking van een reeks van BAM-bestanden, zet de BAM-bestanden in één map en voer de meegeleverde script "HMMpipe.pl" in het pakket aan segmenten te bellen en te berekenen variabiliteit scores (VS). Volg het formaat: perl HMMpipe.pl Folder_to_BAMfiles MM9 Ren DNAcopy om CNVs 21 detecteren. Gebruik SAMtools versie 0.1.19 te extraheren unieke wijze in kaart gebracht leest van elk BAM bestand. samtools bekijken -h-F 0x0004 example.filtered.bam | egrep -i "^ @ | XT: A: U" | samtools bekijken -shū -> example.mapped.bam Gebruik MarkDuplicates in Picard versie 1.94 markeren PCR duplicaten in het BAM-dossier 22. java-jar MarkDuplicates.jar INPUT = example.mapped.bam uitgang = example.nondup.bam METRICS_FILE = example.dup REMOVE_DUPLICATES = true Gebruik coveragebed in bedtools versie 2.17.0 te tellen in kaart gebracht leest in elk van de vooraf gedefinieerde dynamische 500kb mappable ramen in de verstrekte BED bestand "mm9.500k.dynamic.win.bed" of "hg19.500k.dynamic.win.bed" 23 . coverageBed -abam example.nondup.bam -b mm9.500k.dynamic.win.bed -counts | sort -k1,1 -k2,2n> example.nondup.counts Gebruik de meegeleverde Perl-script "normalizeGC.pl" om de read tellingen op basis van de GC inhoud referentiebestand "mm9.500k.dynamic.win.fa.gc.txt" of "hg19.500k.dynamic.win.fa normaliseren. gc.txt ". perl normalizeGC.pl mm9.500k.dynamic.win.fa.gc.txt example.nondup.counts> example.bam.norm.counts Gebruik de meegeleverde R script "dnacopy.r" om de segmenten te bellen. R CMD BATCH dnacopy.r Identificeer CNVs Uitsluiten cellen waarvoor de VS berekende 4.4.6 groter is dan 0,26. Filter de segmenten genoemd HMMcopy in 4.4.6 tot alleen die waarvan de gemiddelde log 2 verhouding groter is dan 0,4 (vermeende gain) of kleiner dan -0,35 (vermeende verlies). Filter de segmenten genoemd DNAcopy 4.5.5 daarvan alleen die waarvoor het segment gemiddelde groter dan 1,32 of kleiner dan 0,6. Overlappen de segmenten van 4.6.2 en 4.6.3, roepen CNVs alleen in gebieden waar beide HMMcopy en DNAcopy identificeren van een vermeende winst of vermeende verlies.

Representative Results

De geassembleerde aspirator moet vergelijkbaar zijn met die in figuur 1A weergegeven. De naald moet zodanig worden opgesteld dat het voldoende breed om een ​​enkele cel tegemoet te komen, maar niet zo breed dat een groot volume wordt opgesteld met de single cell. Opzuigen enkele cellen eenvoudigst wanneer er één tot vijf cellen in een 10X gebied (Figuur 1B). Als hele genoom amplificatie succesvol is, wordt het monster verschijnen als een uitstrijk op een agarosegel (Figuur 2, lanen 1, 2, 4, 5, 6 en 7). Een zwak of afwezig uitstrijkje geeft een foute amplificatiereactie en het monster moet niet worden gesequenced (Figuur 2, lanen 3 en 8). Na bibliotheek voorbereiding moet het fragment grootteverdeling van de monsters beoordeeld door capillaire elektroforese op een fragment analyse.Met succes bereide bibliotheken een tamelijk gelijkmatige verdeling van fragment maten 150-900 bp (Figuur 3, A, B). Mislukt bibliotheek preparaat resulteert in een scheve fragment grootteverdeling en dergelijke bibliotheken niet worden gesequenced (figuur 3C). Het verwerken van sequencing data door middel van de verborgen Markov model (HMMcopy) en circulaire binaire segmentatie (DNAcopy) zal het genoom van elke cel ontleden in segmenten van de geschatte aantal kopieën. Deze segmenten kunnen dan worden gefilterd die voor een aantal kopieën overeenstemming met winst of verlies in een enkele cel (tabel 1) te identificeren. Deze gefilterd segmenten van HMMcopy en DNAcopy moet dan worden overlapt te identificeren met veel vertrouwen CNVs. Figuur 1. Enkele cel isolement. ( A) een samengestelde microaspirator. (B) A 10X veld tonen gescheiden cellen (pijlen) en microaspirator naald (rechtsonder). Opzuigen enkele cellen eenvoudigst wanneer er één tot vijf cellen in een 10x veld. Klik hier om een grotere versie van deze figuur te bekijken. Figuur 2. Hele genoom versterking. Agarose gelelektroforese van 5 pl gehele genoom amplificatieproducten. Monsters die met succes amplificeren zal verschijnen als heldere vlekken van 100 bp tot 1 kb (lanen 1, 2, 4, 5, 6 en 7) en kan worden gesequenced. Monsters die niet met succes kan amplificeren zal zwakke vlekken of geen uitstrijkje (lanen 3 en 8) produceren en mag niet worden gesequenced.large.jpg "target =" _ blank "> Klik hier om een ​​grotere versie van deze figuur te bekijken. Figuur 3. bereiding Library. Representatieve resultaten van een fragment analysator. De grafieken tonen fragmentgrootte (in bp) op X-as en de relatieve fluorescentie-eenheden (RFU) op de Y-as. Rechts van elke grafiek is een gesimuleerde gellaan. (A) de resultaten van een ideaal voorbeeld, met een gelijkmatige verdeling tussen de 150 en 900 bp en zonder scherpe pieken of neiging naar een kant. Dit monster is aanvaardbaar voor sequencing. (B) Resultaten uit een goed monster, met een grootteverdeling scheef naar lagere fragment maten. Hoewel dit niet optimaal, kan het monster nog worden gesequenced. (C) de resultaten van een mislukte monster, met overwegend kleine fragment maten. Dit wordt waarschijnlijk veroorzaakt door langdurige incubatie tijdens de tagmentation stapbibliotheek preparaat. Het monster moet niet worden gesequenced. Klik hier om een grotere versie van deze figuur te bekijken. HMMcopy Monster Segment chr Begin Einde Staat Mediaan D15-4998 23 chr8 144.500.001 146.500.000 6 0.5008794 D15-4998 29 chr10 67.000.001 134.500.000 6 0.4031945 D15-4998 52 chr19 1 20000000 6 0.4616884 D15-4998 57 Chry 1 59500000 2 -1,506532 DNAcopy Monster Chrom Start bin end bin Begin Einde Seg gemiddelde Seg med D15-4998 chr10 88 197 62.612.945 129.971.511 1,4688 0,157 D15-4998 chr19 0 31 0 28.416.392 1,4141 0,1674 D15-4998 chrX 77 126 </td> 51.659.160 95.343.369 1,3548 0,1874 D15-4998 Chry 0 14 0 23.805.358 -2,7004 0,3591 Tabel 1. Gegevensanalyse. Gefilterd segmenten gegenereerd door HMMcopy en DNAcopy van een enkele cel. Overlappende Resulteert onthult een winst op chromosoom 10 67-130 Mb en een winst op chromosoom 19 van 0 tot 20 Mb. We hebben gevonden dat de voordelen op het proximale gedeelte van chromosoom 19 zijn de invloed van enkele cel 12 sequencing.

Discussion

Traditioneel, het identificeren van CNVs en aneuploïdie op het niveau van enkele cellen nodig cytologisch methoden zoals FISH en SKY. Nu heeft enkele cel sequencing ontpopt als een alternatieve aanpak voor dergelijke vragen. Enkele cel sequencing heeft voordelen ten opzichte van FISH en SKY want het is zowel genoom-brede en hoge resolutie. Bovendien, wanneer de juiste kwaliteitscontrole methoden worden toegepast, één cel sequencing kan een meer betrouwbare beoordeling van CNVs en aneuploidy bieden want het is niet gevoelig voor de hybridisatie en het verspreiden van artefacten die inherent zijn aan vis en SKY. Echter, veel van de recente toepassingen van enkele cel sequencing niet zijn onderbouwd door een grondige evaluatie van de gevoeligheid en specificiteit van de methoden en analyses. Inderdaad, enkele analytische methoden dat door andere studies zijn geassocieerd met een hoge frequentie (> 50%) van valspositieve CNV roept 12. De aanpak die we beschrijven is uitvoerig getest met behulp van cellen van kneigen CNV last om ware en valse ontdekking tarieven te bepalen en optimaliseren van de gevoeligheid en specificiteit van CNV detectie 12. Met de kwaliteitscontrole en analytische methoden beschreven in dit protocol, ongeveer 20% van 5 Mb winst, 75% van 5 Mb verliezen, en alle CNV dan 10 Mb kunnen worden gedetecteerd. Hoewel de bepaling van de valse ontdekking tarief van eencellige sequencing moeilijk hebben geschat dat minder dan 25% bedragen. Dit protocol kan worden toegepast op cellen van verschillende bronnen, kan de scripts worden aangepast aan de resolutie van CNV detectie passen en het protocol kan worden aangepast aan andere soorten genomische veranderingen te identificeren.

Er zijn verscheidene middelen dissociëren verse weefsels in enkele cellen en vele publicaties beschrijven procedures geoptimaliseerd voor specifieke weefsels, zoals huid en hersenen 24 25. Wij verkiezen enkele cellen isoleren microaspiration omdat hiermee fo r visuele beoordeling van elke cel te sequeneren. Het is echter ook mogelijk om enkele cellen te isoleren door fluorescentie-geactiveerde celsortering (FACS) 26 en 27 microfluïdische inrichtingen. Als één cel isolatie en gehele genoom amplificatie handmatig wordt uitgevoerd, is het redelijk te isoleren en amplificeren tot veertig cellen in een enkele zitting. Om een ​​hoge kwaliteit cel sequentiedata te verkrijgen, is het cruciaal dat de amplificatie van eencellige genomen uniform en volledig. We zien dat de kwaliteit van afzonderlijke cellen en die de efficiëntie van lysis en amplificatie een aanzienlijke invloed op de kwaliteit van sequentiebepalingsgegevens. Als zodanig moet de cellen worden geoogst uit hun natieve omgeving juist voorafgaande aan isolatie en gehele genoom amplificatie dient onmiddellijk na cellen geïsoleerd. Bovendien moet de lysis en fragmentatie stap precies gevolgd zoals beschreven in stappen 2.3-2.6.

"> De algoritmes kan worden aangepast om de resolutie van CNV detectiewijziging, met tegengestelde effecten op de gevoeligheid en specificiteit 12. Het is ook mogelijk om het bij te stellen om globale chromosoom aneuploidie te detecteren in de omgeving van tetraploidy 11. We zien echter dat onze aanpak is beperkt tot het opsporen van CNVs groter dan 5 Mb, als ruis geïntroduceerd tijdens hele genoom amplificatie bemoeilijkt de opsporing van kleinere varianten 12. Toekomstige verbeteringen in de hele genoom amplificatie benaderingen moeten uiteindelijk verbeteren van de resolutie van CNV detectie met behulp van enkele cel sequencing.

Enkele cel sequencing maakt onderzoek niet alleen kopieaantal veranderingen maar één nucleotide verschillen 28, 29 en 30 structuurvariatie. Ons protocol voor cel isolatie kan worden toegepast op deze andere que beantwoordenstions. De keuze van hele genoom amplificatie is afhankelijk van de specifieke toepassing. De werkwijze in dit protocol, dat is gebaseerd op polymerasekettingreactie beschreven, is het meest geschikt voor het detecteren kopienummer veranderingen, omdat het wordt geassocieerd met lagere amplificatie vertekening 32. Voor onderzoek naar andere soorten genomische veranderingen, zoals single nucleotide polymorphisms, zijn andere werkwijzen gehele genoom amplificatie aangenomen geschikter 31, 32 te zijn.

Disclosures

The authors have nothing to disclose.

Acknowledgements

Wij danken Stuart Levine voor reacties op dit manuscript. Dit werk werd ondersteund door de National Institutes of Health Grant GM056800 en Kathy en Curt Marble Cancer Research Fund om Angelika Amon en gedeeltelijk door de Koch Instituut Ondersteuning Grant P30-CA14051. Angelika Amon is ook een onderzoeker van het Howard Hughes Medical Institute en de Glenn Foundation for Biomedical Research. KAK wordt ondersteund door de NIGMS Training Grant T32GM007753.

Materials

Aspirator tube assemblies for calibrated microcapillary pipettes Sigma A5177
PVC tubing (ID 3/16", OD 5/16", ~1 foot) VWR 89068-500
PVC tubing (ID 5/16", OD 7/16", ~6 inches) VWR 89068-508
Acrodisc Syringe Filter with HT Tuffryn Membrane (diameter 25 mm, pore size 0.2 um) VWR 28144-040
Serological pipette (5 ml) BioExpress P-2837-5 Any plastic 5 mL pipette should suffice
In-Line Water Trap for Oxygen Use Amazon.com 700220210813
Capillary Melting Point Tubes (ID 0.8-1.1 mm, 100 mm length) VWR 34502-99
Modeling clay VWR 470156-850
Petri dish (150 mm diameter) VWR 25384-326 Surface should be non-adherent to facilitate aspiration of single cells
Hard-Shell Full-Height 96-Well Semi-Skirted PCR Plates Bio-Rad HSS9601 Any 96-well PCR plate should suffice
96-well tissue culture plate VWR 62406-081 Any 96-well tissue culture plate should suffice
GenomePlex Single Cell Whole Genome Amplification Kit Sigma WGA4
Microseal 'A' Film Bio-Rad MSA5001
Mini plate spinner Thomas Scientific 1225Z37
Thermal cycler Bio-Rad 1861096  Any thermal cycler should suffice so long as it can accommodate a 96-well PCR plate
Agencourt Ampure Beads Beckman Coulter A63880
Dynamag-2 Magnet Thermo Fisher Scientific 12321D Any similar magnetic tube strip should suffice
Nextera XT DNA Library Preparation Kit Illumina FC-131-1096
Nextera XT Index Kit Illumina FC-121-1012
Complete kit (optimized for Roche LightCycler 480) Kapa Biosystems KK4845  This kit is optimized for the Roche LightCycler 480 real-time PCR machine. If using a different machine, substitute a kit optimized for your machine.

References

  1. Kahlem, P. Transcript Level Alterations Reflect Gene Dosage Effects Across Multiple Tissues in a Mouse Model of Down Syndrome. Genome Res. 14 (7), 1258-1267 (2004).
  2. Torres, E. M. Effects of Aneuploidy on Cellular Physiology and Cell Division in Haploid Yeast. Science. 317 (5840), 916-924 (2007).
  3. Itsara, A. Population Analysis of Large Copy Number Variants and Hotspots of Human Genetic Disease. Am. J. Human Gen. 84 (2), 148-161 (2009).
  4. Sudmant, P. H. Global diversity, population stratification, and selection of human copy-number variation. Science. 349 (6253), aab3761 (2015).
  5. Navin, N. Inferring tumor progression from genomic heterogeneity. Genome Res. 20 (1), 68-80 (2010).
  6. Torres, L., Ribeiro, F. R., Pandis, N., Andersen, J. A., Heim, S., Teixeira, M. R. Intratumor genomic heterogeneity in breast cancer with clonal divergence between primary carcinomas and lymph node metastases. Breast cancer res and treat. 102 (2), 143-155 (2007).
  7. Yates, L. R. Subclonal diversification of primary breast cancer revealed by multiregion sequencing. Nat. Med. 21 (7), 751-759 (2015).
  8. Forsberg, L. A. Age-related somatic structural changes in the nuclear genome of human blood cells. Am. J. Human Gen. 90 (2), 217-228 (2012).
  9. Laurie, C. C. Detectable clonal mosaicism from birth to old age and its relationship to cancer. Nat. Gen. 44 (6), 642-650 (2012).
  10. Jacobs, K. B. Detectable clonal mosaicism and its relationship to aging and cancer. Nat. Gen. 44 (6), 651-658 (2012).
  11. Knouse, K. A., Wu, J., Whittaker, C. A., Amon, A. Single cell sequencing reveals low levels of aneuploidy across mammalian tissues. Proc. Natl. Acad. Sci. 111 (37), 13409-13414 (2014).
  12. Knouse, K. A., Wu, J., Amon, A. Assessment of megabase-scale somatic copy number variation using single-cell sequencing. Genome Res. 26 (3), 376-384 (2016).
  13. . KAPA Library Quantification Kit Available from: https://www.kapabiosystems.com/product-applications/products/next-generation-sequencing-2/library-quantification/ (2016)
  14. . DNAcopy Available from: https://bioconductor.org/packages/release/bioc/html/DNAcopy.html (2016)
  15. Lichti, U., Anders, J., Yuspa, S. H. Isolation and short-term culture of primary keratinocytes, hair follicle populations and dermal cells from newborn mice and keratinocytes from adult mice for in vitro analysis and for grafting to immunodeficient mice. Nat. Protoc. 3 (5), 799-810 (2008).
  16. Brewer, G. J., Torricelli, J. R. Isolation and culture of adult neurons and neurospheres. Nat. Protoc. 2 (6), 1490-1498 (2007).
  17. Navin, N. Tumour evolution inferred by single-cell sequencing. Nature. 472 (7341), 90-94 (2011).
  18. Szulwach, K. E., Chen, P. Single-Cell Genetic Analysis Using Automated Microfluidics to Resolve Somatic Mosaicism. PLoS ONE. 10 (8), e0135007 (2015).
  19. Zong, C., Lu, S., Chapman, A. R., Xie, X. S. Genome-Wide Detection of Single-Nucleotide and Copy-Number Variations of a Single Human Cell. Science. 338 (6114), 1622-1626 (2012).
  20. Wang, Y. Clonal evolution in breast cancer revealed by single nucleus genome sequencing. Nature. 512 (7513), 155-160 (2014).
  21. Zhang, C. Z. Chromothripsis from DNA damage in micronuclei. Nature. 522 (7555), 179-184 (2015).
  22. Macaulay, I. C., Voet, T. Single Cell Genomics: Advances and Future Perspectives. PLoS Genetics. 10 (1), e1004126 (2014).
  23. Gawad, C., Koh, W., Quake, S. R. Single-cell genome sequencing: current state of the science. Nat. Rev. Gen. , 1-14 (2016).

Play Video

Cite This Article
Knouse, K. A., Wu, J., Hendricks, A. Detection of Copy Number Alterations Using Single Cell Sequencing. J. Vis. Exp. (120), e55143, doi:10.3791/55143 (2017).

View Video