ENHANCED ridotto Rappresentanza bisolfito sequenziamento per la valutazione della metilazione del DNA in Risoluzione coppia di basi

Sommario

Il sequenziamento del bisolfito a rappresentazione ridotta potenziata è un metodo per la preparazione di librerie di sequenziamento per l'analisi della metilazione del DNA basata sulla digestione dell'enzima di restrizione combinata con la conversione del bisolfito di citosina. Questo protocollo richiede 50 ng di materiale di partenza e produce dati di risoluzione delle coppie di basi in regioni genomiche ricche di GC.

Abstract

DNA metilazione mapping è fortemente studiato in tessuti normali e patologici. Una varietà di metodi sono stati stabiliti per interrogare i modelli citosina metilazione nelle cellule. Rappresentazione ridotta di sequenziamento dell'intero genoma bisolfito è stato sviluppato per individuare coppie di basi modelli risoluzione citosina metilazione quantitativi a ricchi di GC loci genomici. Questo si ottiene combinando l'uso di un enzima di restrizione seguita da conversione bisolfito. Maggiore ridotto Rappresentanza bisolfito Sequencing (ERRBS) aumenta la loci genomici biologicamente rilevanti coperti ed è stato usato per profilare citosina metilazione del DNA da umano, mouse e altri organismi. ERRBS inizia con enzimi di restrizione del DNA per generare frammenti a basso peso molecolare per la redazione biblioteca. Questi frammenti sono sottoposti alla costruzione della libreria standard per la prossima generazione di sequenziamento. Conversione bisolfito di cytosines non metilato prima della amplificati finalesul gradino permette per la risoluzione di base quantitativa dei livelli di metilazione citosina in coperta loci genomici. Il protocollo può essere completata nel giro di quattro giorni. Nonostante bassa complessità nei primi tre basi in sequenza, le biblioteche ERRBS producono dati di alta qualità quando si utilizza una corsia di controllo sequenza designato. Mappatura e bioinformatica analisi viene quindi eseguita e dati produce che può essere facilmente integrato con una varietà di piattaforme genoma. ERRBS può utilizzare piccole quantità di materiale in entrata rendendo così possibile elaborare campioni clinici umani e applicabile in una vasta gamma di applicazioni di ricerca. Il video prodotto dimostra passaggi critici del protocollo ERRBS.

Introduzione

La metilazione del DNA a citosina (5-metilcitosina) è un marchio epigenetico critica in cellule di mammifero per una varietà di processi biologici, compreso ma non limitato a imprinting, inattivazione del cromosoma X, lo sviluppo, e la regolazione dell'espressione genica ^1-8. Lo studio della metilazione del DNA in maligne e altri disturbi ha determinato malattie modelli specifici e contribuito alla comprensione della patogenesi della malattia e potenziali scoperte biomarker ^9-17. Ci sono molti protocolli che interrogano dell'epigenoma per lo stato di metilazione del DNA. Questi possono essere suddivisi in base affinità, enzima di restrizione-based, e saggi basati conversione-bisolfito che utilizzano le piattaforme di microarray o sequenziamento valle. Inoltre, ci sono alcuni protocolli che colmano queste categorie generali compresi, ma non limitatamente a, combinato bisolfito Restriction Analysis ¹⁸ e ridotto Rappresentanza bisolfito Sequencing (RRBs ^19).

RRBs stato originariamente descritto da Meissner et al. ^19,20. Il protocollo ha introdotto un passo per arricchire regioni genomiche ricche di GC seguiti da sequenziamento bisolfito, che ha portato i dati risoluzione di base-pair quantitativo, che è economicamente vantaggioso ^21,22. Le regioni ricche di GC sono mirati dal MspI (C ^ CGG) enzima di restrizione, e citosina metilazione si risolve con la conversione bisolfito di cytosines (deamination di cytosines modificati a uracile), seguita dalla reazione a catena della polimerasi (PCR). RRBs coperto la maggioranza dei promotori genici e isole CpG in una frazione del sequenziamento richiesto per un intero genoma; tuttavia RRBs avevano copertura limitata delle coste CPG e altre regioni intergeniche di rilevanza biologica. Diversi gruppi hanno pubblicato aggiornati protocolli RRBs rispetto alla relazione originale che migliorano la metodologia e la copertura risultante di queste regioni genomiche ^23-25. Maggiore rappresentanza ridotta Bisulfite Sequscere (ERRBS) include modifiche preparazione di libreria e di un allineamento di dati alternativo approccio ²⁶ quando rispetto al RRBs. ERRBS provocato un maggior numero di CPGs rappresentata nei dati generati e aumentato la copertura di tutte le regioni genomiche interrogato ^26. Questo metodo è stato utilizzato per risolvere schemi di metilazione del DNA in paziente umano e altri esemplari di animali ^26-30.

Il protocollo ERRBS descritto particolari offerte su tutti i passi necessari per il completamento e dati sono stati generati utilizzando rappresentante DNA umano (i campioni sono stati ottenuti da precedentemente riportati, i campioni dei pazienti de-identificati ^31, e CD34 + midollo osseo da un campione normale donatore umano). Il protocollo comprende una procedura automatica di selezione del formato, che riduce il tempo di elaborazione per campione e consente una maggiore precisione nella selezione della libreria dimensioni. Il protocollo combina una serie di tecniche di biologia molecolare stabiliti. Alto peso molecolare DNA viene digerito won una metilazione-insensitive enzima di restrizione (MspI) seguita da fine-riparazione, A-tailing, e legatura di adattatori denaturato. Dimensione selezione dei frammenti ricchi di GC è seguita da conversione bisolfito e PCR di amplificazione prima sequenziamento. Conversione bisolfito è stato descritto in precedenza ³² e dettagliata revisione di analisi e applicazioni dati è oltre la portata di questo articolo, tuttavia consigli e riferimenti sono inclusi per l'uso dei lettori. Il protocollo può essere eseguita in quattro giorni ed è suscettibile di piccolo ingresso (50 ng o meno) importi rilevanti. Il protocollo, come descritto produce dati con elevata copertura per sito CpG sufficiente non solo per differenziali sito metilazione e regione determinazioni, ma anche per il rilevamento del polimorfismo epigenetico come descritto da Landan, et al. ^33.

Protocollo

Tutte le procedure eseguite sono approvati dalla Scuola di Medicina Istituzionale Animal Care and Use Committee Indiana University e seguono Istituto Nazionale di linee guida di salute.

1. Tecnica chirurgica

1. Mantenere una tecnica asettica durante questa procedura utilizzando guanti sterili, strumenti, e un campo chirurgico sterile, secondo le linee guida NIH ^25. Sterilizzare gli strumenti prima di iniziare l'intervento chirurgico da loro autoclave (vedi Tabella di reagenti specifici / Attrezzature per l'elenco completo). Utilizzare uno sterilizzatore perle di vetro per sterilizzare strumenti durante l'operazione.

2. Anestesia e preparazione

1. Anestetizzare il mouse in una scatola anestesia con una miscela di 0,9 L / min ossigeno e 2,5% isoflurano utilizzando un sistema vaporizzatore isoflurano veterinaria. Assicurarsi che il mouse non risponde ai cambiamenti di posizione del corpo prima di rimuoverlo dalla scatola.
2. Applicare una pomata oftalmica al mouGli occhi di SE per proteggerli dalla disidratazione.
3. Passare il flusso di gas dalla casella a cono. Posizionare il mouse esattamente sul lato sinistro su un tappetino riscaldato coperto da un tappetino chirurgica e carta assorbente con panca il naso e la bocca all'interno del cono. Monitorare costantemente il ritmo del respiro del mouse e la velocità e regolare i livelli isoflurano come necessario (tra 2,5-3% isoflurano) per mantenere un adeguato livello di anestesia, e utilizzare la punta pizzico reflex per confermare sedazione totale.

3. Approccio chirurgico

1. Allineare e mettere a fuoco lo stereoscopio con il campo operatorio. Regolare il cono e il nastro verso il basso in modo che sia posizionato lungo il bordo del campo visivo.
2. Con il mouse sdraiato sul suo lato sinistro, il nastro il bordo dell'orecchio destro al cono, esponendo l'area dietro l'orecchio in cui verrà effettuata l'incisione. Assicurarsi che la vena auricolare posteriore percorre orizzontalmente l'orecchio. Si noti che il corretto posizionamento di tegli animali e taping dell'orecchio sono fondamentali al fine di trovare rapidamente il nervo facciale.
3. Bagnare il pelo sopra e dietro l'orecchio con il 70% di etanolo e radere il sito chirurgico con una lama di rasoio o bisturi. Pre-bagnare il pelo rende la rasatura più facile in questa posizione anatomica.
4. Pulire la pelle con una soluzione di iodio, come Betadine lavaggio chirurgico (7,5% povidone-iodio), seguito da 70% di etanolo. Ripetere questa pulizia due volte di più per disinfettare accuratamente la zona.
5. Per determinare dove fare l'incisione, tracciare la vena auricolare posteriore dall'orecchio caudalmente alla zona posteriore alla protuberanza orecchio. Utilizzando forbici primavera, fare un 4 millimetri un'incisione 2 - 3 mm posteriormente alla protuberanza.
6. Sezionare attraverso il grasso sottocutaneo e la fascia con smussa. Evitare il taglio diretto con le forbici, perché i vasi sanguigni o tessuto muscolare possono essere facilmente danneggiati.
7. In caso di emorragia, applicare pressione al sito chirurgico con un tampone di cotone sterileper almeno 30 secondi. Se si verifica una significativa perdita di fluido, iniettare il mouse per via intraperitoneale con fino a 0,5 ml di soluzione salina sterile allo 0,9% utilizzando un ago G 25 o 27.
8. Utilizzare diversi punti di riferimento chiave, il nervo accessorio spinale, canale uditivo, e anteriore del muscolo digastrico (descritto di seguito), per individuare il nervo facciale. Sezionare intorno a questi punti di riferimento fino a quando sono visualizzati i rami del nervo facciale. Il nervo apparirà come una significativa struttura di solido bianco quando viene rivelato e uno strato di fascia aderisce agli strutture sottostanti.
   1. Trova il nervo accessorio spinale, che viaggia dalla porzione caudale del cranio per innervare il muscolo trapezio, una volta che il grasso sottocutaneo e la fascia sono stati sezionati. Il nervo facciale è profondo al nervo accessorio spinale.
   2. Trovare il canale uditivo cartilagineo che sembra bianco perlato e può essere visto rostrale al nervo facciale.
   3. Trova il ventre muscolare del muscolo digastrico anteriore, che si trova in cima e caudal al nervo facciale.
9. Quando i principali rami del nervo facciale sono visualizzati, li risalire dorsale di trovare la loro origine dal forame stylomastoid. Utilizzando punta fine Dumont pinze # 5/45 di tenere il sito chirurgico aperto, far avanzare le punte primavera forbice seguito il percorso del nervo, quindi spostare la pinza dorsale per mantenere la nuova area Advanced Open.
10. Visualizza il tronco del nervo facciale con zigomatico, buccale, e rami mandibolari marginali a questo punto.
    NOTA: Il ramo temporale si troverà più vicino al forame. I rami marginali nervo mandibolare nelle sue parti superiore e inferiore più vicino alla mascella, pertanto quei rami nervose non saranno visibili a questo livello.
    1. Se si esegue una resezione del nervo, di stabilizzare il nervo delicatamente con le pinze punta fine e tagliare il nervo con le forbici a molla. Evitare di applicare troppa trazione al nervo con le pinze per evitare avulsing nervo dal tronco encefalico. Spingerei monconi distanti, o tagliare e rimuovere una porzione del nervo distale per garantire che nessun riconnessione può verificarsi.
    2. Se si esegue una lesione da schiacciamento, utilizzare Dumont # 5/45 pinze per comprimere il nervo per 30 secondi con una pressione costante a recidere tutti gli assoni, quindi ripetere questa cotta in un secondo angolo perpendicolare al primo sito cotta. Evitare l'applicazione di quantità variabili di pressione durante la 30 sec cotta, altrimenti il ​​danno non sarà coerente tra gli animali.

4. Chiusura e recupero

1. Riposizionare il grasso e muscoli sopra le strutture sottostanti.
2. Approssima i bordi dell'incisione e chiudere la ferita utilizzando un clip ferita 7,5 millimetri. Suture o colla sono accettabili per la chiusura della ferita. Analgesici postoperatoria possono essere forniti in questo momento.
3. Rimuovere il nastro dal orecchio del mouse. Spegnere il flusso isoflurano e lasciare il mouse per respirare ossigeno puro per 30 secondi a 1 minuto. Plasso il mouse in una gabbia vuota senza biancheria da letto per recuperare da anestesia.
4. Quando il mouse è recuperato, esaminare il suo comportamento per i segni di conferma della paralisi facciale. I baffi saranno paralizzati e angolati verso la guancia, il naso viene deviato, e l'occhio non lampeggia in risposta ad un soffio d'aria.
5. Animali Casa congiuntamente dopo l'intervento chirurgico se sono di sesso femminile. Evitare ospitano topi maschi congiuntamente perché sono più aggressivi e tendono a rimuovere forzatamente clip ferita del loro cagemate, che porta a infezioni. Fornire analgesici postchirurgiche in questo momento, se necessario.
6. Monitorare i topi una volta al giorno per diversi giorni dopo l'operazione per assicurarsi che nessuna infezione o altre complicanze si verifica dopo l'intervento. Rimuovere clip ferita 7 - 10 giorni dopo l'intervento chirurgico, se non sono caduti fuori da sè.
7. Applicare lubrificante pomata occhio l'occhio interessato al giorno per prevenire le complicanze corneali, sia fino a quando il riflesso corneale occhio è recoperti o fino eutanasia.

Risultati

La figura 1 fornisce una panoramica di ERRBS, evidenziando passaggi chiave, che sono spiegate in tutto il protocollo descritto. Biblioteche ERRBS sono stati preparati utilizzando 50 ng DNA di ingresso.

Valutare la qualità delle biblioteche preparati. Produzione Libreria produce abitualmente formati frazione di 150-250 bp e 250-400 bp (Figura 3A-C). Sono attesi Lievi differenze nella distribuzione delle librerie di dimensioni tra i campioni. Si noti che in entrambe le frazioni biblioteca inferiori e superiori non sono formati DNA molto intense, indicativi di arricchimento di una particolare sequenza. MspI risultati digestione l'arricchimento di una famiglia di DNA ripetitivo sequenze presenti nel genoma umano a 190 bp, 250 bp e 310 bp nelle biblioteche ERRBS. Questi tre ripetizioni rappresentano una firma caratteristica di una libreria ERRBS ²⁰ (vedi figure 3A-C e 3G). Biblioteche rappresentativi sono stati sequenziati su una sequenza di prossima generazioner utilizzando single-end legge. Quando si carica alla concentrazione biblioteca consigliata su Illumina HiSeq 2500 sequencer, sono attesi densità di cluster di 500,000-700,000 per mm ^2. A questa densità clustering, 81,6% ± 3,14% (n = 81) dei cluster filtro passa (Figura 4A). A causa della fine complessità basso degli inserti di libreria (sito di riconoscimento MspI: C ^ CGG), i valori di intensità e punteggi di qualità registrati durante il sequenziamento sono altamente variabili nei primi tre basi (Figura 4B-C), tuttavia, se una corsia di controllo indipendente è incluso (vedi la discussione), l'85% delle basi avranno punteggi di qualità di 30 o superiore (valori Q30; Figura 4D).

L'allineamento dei dati e determinazione citosina metilazione come descritto nei dati risoluzione rendimenti protocollo base-pair (Tabella 7). Per il genoma umano, un 51-ciclo unico lettura run sequenziamento di una libreria ERRBS in una corsia di un HiSeq 2500 in modalità ad alto rendimento normalely genera 153.194.882 ± 12.918.302 totale si legge che dopo il filtraggio di qualità e l'adattatore taglio rendimenti 152.231.183 ± 13.189.678 legge per l'ingresso nella pipeline di analisi. Rendimento medio mappatura per una libreria ERRBS è tipicamente 62.95% ± 5,92% con la rappresentazione di 3.183.594 ± 713.547 CPGs con una copertura minima per CpG di 10x e una copertura media per CpG di 84,94 ± 16,29 (n = 100).

Il protocollo ERRBS è suscettibile di multiplexing (vedi file Supplemento 1: Protocollo di adattamento per il sequenziamento multiplex). I dati di sequenziamento rappresentante corre è riassunto nella figura 5 Dati da piste di sequenziamento multiplex (51-ciclo di esecuzione singola lettura sequenziamento;. N = 128 per due biblioteche per corsia; n = 11 per tre biblioteche per corsia; n = 11 per quattro biblioteche per corsia) sono stati confrontati con una corsia sequenza completa di una libreria ERRBS (51-ciclo piste sequenziamento singoli-letto; n = 100) e downsampling una sola corsia per SIMULAZIONEe 50%, 33% e 25% di letture per corsia (2, 3, e 4 campione multiplexing per corsia rispettivamente; n = 3). Poiché il numero di letture per campione diminuisce con il fattore di multiplexing, il numero di CPGs coperte ad una copertura minima di 10x e la copertura per CpG diminuisce pure (Figura 5 e Tabella 8). Significa tassi di conversione dei siti non-CpG attesi sono 99,85% ± 0,04% (n = 400). I tassi di conversione più bassi del 99% possono indicare non ottimale di conversione bisolfito che può causare alti tassi di livelli di metilazione falsi.

I dati di una libreria ERRBS preparato da un DNA rappresentante umano genomico è stato analizzato in R 2.15.2 ⁴⁵ utilizzando il pacchetto methylKit ²⁶ (vedi file di codice supplementare 1 per i dettagli di comando). I dati possono essere visualizzati nei browser genoma di uso comune (Figura 6A). I dati citosina metilazione è equamente derivato da entrambi i filamenti (Figura 6B) e va tuttospettro di potenziali livelli di metilazione citosina (Figura 6c). Analisi delle repliche tecniche da un alto rappresentante concordanza rendimenti campione di DNA umano tra i risultati dei dati (Figura 6D) e copre CPGs in un ampio spettro di loci genomici (Figura 6 E e F e, come descritto in precedenza ^26). Mentre repliche tecniche produrranno alta R ² valori (superiore al 97%), repliche biologica produrrà R ² valori compresi 0,92-0,96 ^26, e confrontando i differenti tipi di cellule umane produrrà R ² valori inferiori a 0.86 (dati non mostrati).

Figura 1: Diagramma di flusso delle fasi di protocollo ERRBS. Grafico rappresenta passi, che può essere completato in un giorno di lavoro tradizionale. * Indica un potenziale punto di pausa (immediatamente seguire ing legatura ripulire e prima selezione del formato, protocollo punto 5) in cui i campioni possono essere congelati a -20 ° C prima di procedere con la durata di protocollo.

Figura 2: protocollo selezione Size. (A) Schermata di impostazioni utilizzate nel protocollo ERRBS Pippin Prep (vedi sezione protocollo 5.1.2 - 5.1.6): (1) Selezionare il tipo di cassetta. (2) Selezionare tipo da impiegare. (3) Selezionare la modalità di raccolta per ogni corsia. (4) Inserire i campi di raccolta bp. (5) Salvare il protocollo (B) Fasi di estrazione gel manuale utilizzato nella sezione del protocollo 5.2:. (1) scale gel visualizzati. (2) Contrassegnato formati per la selezione dimensioni utilizzando una lama di rasoio. (3) L'immagine dei campioni asportati (frazione inferiore: 150-250 bp e superiori frazione: 250-400 bp)."> Clicca qui per vedere una versione più grande di questa figura.

Figura 3: i risultati dei controlli di qualità per biblioteche rappresentative ERRBS preparati da campioni di DNA umano con una macchina Bioanalyzer. (A) immagine Gel-come mostra una scala standard (1), frazione libreria inferiore (135-240 bp frazione da Pipino Prep); 2) e la frazione biblioteca superiore (240-410 bp frazione da Pipino Prep); . 3) (B) Bioanalyzer elettroferogramma della frazione biblioteca previsto minor (C) Bioanalyzer elettroferogramma della frazione biblioteca più alta prevista D -.. F) i dati rappresentativi da una povera libreria qualità prep. Immagine Gel-like (D) della scala standard (1), frazione inferiore library (2) e la frazione biblioteca superiore (3). La banda a 150 bp contrassegnati con arrow indica una quantità eccessiva di adattatore. Elettroferogramma del minore (E) e le frazioni più alte di libreria (F) con i picchi adattatore in eccesso a 150 bp (segnato con frecce). (G) Bioanalyzer elettroferogramma di una libreria ERRBS pool per il sequenziamento. Traccia rossa rappresenta una libreria pool di alta qualità con una pari rappresentanza di frazioni superiori e inferiori. Blu traccia rappresenta una libreria pool non adeguato per il sequenziamento a causa della mancanza di pari rappresentanza delle frazioni più alte e più basse. Clicca qui per vedere una versione più grande di questa figura.

Figura 4: Sequencing classifica per un rappresentante ERRBS 51-ciclo di esecuzione sequenziamento single-letto su un HiSeq 2500 sequencer in alto rendimentomodalità. (A) densità Cluster (K / mm ² = 1.000 cluster per millimetro quadrato; blu). E densità di cluster passaggio filtrante (verde) in due corsie con librerie ERRBS (B) intensità tipiche visto nei primi 30 cicli in una corsia con una biblioteca ERRBS. Nota la firma CGG da MspI digestione nelle intensità dei primi tre cicli. (C) Percentuale di basi con un punteggio di qualità di 30 o superiore (%> Q30) per ogni ciclo in una corsia ERRBS. (D) Distribuzione punteggio di qualità per tutti i cicli di una corsia ERRBS. Blu = meno di Q30, Verde = maggiore o uguale a Q30. In questa corsia, 84,7% di basi avevano punteggi di qualità di 30 o superiore.

Figura 5: Sequencing risultati di output. Trame Box di dati sperimentali provenienti da campione multiplex e singola per corsia sequenziamento ru ns (visualizzati come scatole verdi) e di dati ricavati da downsampling simulato dal sequenziamento piste di tre ERRBS librerie (visualizzato come scatole blu, campione di cinque volte per ogni ciclo di sequenziamento) da 51 a ciclo sequenziamento singola lettura corre Il fattore multiplexing corrisponde. il numero di biblioteche ERRBS sequenziato per corsia. 1 = corsia intero o 100% di legge e rappresenta i dati di una singola libreria ERRBS per corsia; 2 = 50% di corsia e rappresenta i dati da due ERRBS biblioteche per corsia; 3 = 33% di una corsia e rappresenta i dati da tre ERRBS biblioteche per corsia; e, 4 = 25% di una corsia e rappresenta i dati provenienti da quattro ERRBS librerie per corsia. (A) I conti di lettura, o il numero di sequenze analizzate, per ogni fattore di multiplexing. (B) Il numero di CpG di coperta dai dati di sequenziamento per multiplexing factor. (C) La copertura media per CpG per fattore multiplexing._blank "> Clicca qui per vedere una versione più grande di questa figura.

Figura 6: i dati rappresentativi da una libreria ERRBS preparata dal DNA genomico umano (A) University of California, Santa Cruz (UCSC) Browser genoma ⁴³ immagine di dati rappresentativi di una corsia di ERRBS sequenziamento.. La barra di scala asse y rappresenta 0-100% metilazione ad ogni citosina coperto con un minimo di 10x. La pista su ordinazione superiore rappresenta il filo in avanti e la pista personalizzata inferiore rappresenta il filo opposto. Viene mostrato chr12:.. 6,489,523-6,802,422 (hg19) comprensivo di geni RefSeq e isole CpG all'interno di questa regione genomica istogrammi (B) Distribuzione di copertura CpG lungo in avanti e indietro i fili in un CD34 + umane del midollo osseo campione rappresentativo (C) istogramma di distribuzionedei livelli di metilazione CpG lungo entrambi i filamenti in un CD34 umana + midollo osseo campione rappresentativo. (D) Correlazione appezzamento di livelli di metilazione CpG da una replica tecnico rappresentativo di un campione di DNA umano. grafico (E) Pie che illustra le proporzioni di CPGs coperti in ERRBS che annotata in isole CpG (verde chiaro), CpG rive (grigio) e in altre regioni (bianco) in un campione rappresentativo preparata dal DNA genomico umano. grafico (F) Pie che illustra le proporzioni di CPGs coperti in ERRBS che annotati a promotori di geni (rosso ), esoni (verde), introni (blu) e regioni intergeniche (viola). Clicca qui per vedere una versione più grande di questa figura.

Reagente	Volume	Commento
10x T4 DNA Ligase Buffer Reaction	10 microlitri
Trifosfato deossinucleotidi (dNTP) Soluzione Mix	4 pl	mix di 10 mm di ogni nucleotide
T4 DNA Polymerase	5 ml	3.000 unità / ml
DNA polimerasi I Large (Klenow) Fragment	1 ml	5.000 unità / ml
T4 Polynucleotide Kinase	5 ml	10.000 unità / ml
Acqua priva di DNase	45 microlitri

Tabella 1:. End reagenti reazione riparazione nomi dei reagenti e le quantità utilizzate nella reazione di riparazione fine (protocollo punto 2.1).

Reagente	Volume	Commento
Tampone di reazione 10x	5 ml	per esempio, NEBuffer 2
1 mM 2'-deossiadenosina 5'-trifosfato (dATP)	10 microlitri
Klenow Fragment (3 '→ 5' eso)	3 ml	5.000 unità / ml

Tabella 2:. A-tailing reagenti reazione nomi dei reagenti e le quantità utilizzate nella reazione A-tailing (protocollo punto 3.1).

Reagente	Volume	Commento
15 micron adattatori ricotto in acqua priva di DNase	3 ml	PE adattatore 1.0 e PE adattatore 2.0; vedi tabella 4 per le sequenze e riferimento
10x T4 DNA Ligase Buffer Reaction	581; l
T4 DNA Ligase	1 ml	2.000.000 unità / ml
Acqua priva di DNase	31 microlitri

Tabella 3: Adattatore legatura reagenti reazione nomi dei reagenti e le quantità utilizzate nella reazione adattatore legatura (protocollo punto 4.2)..

Tabella 4:. Oligos utilizzati nel protocollo ERRBS Elenco dei oligonucleotidi utilizzati nel protocollo ERRBS nella reazione legatura (protocollo punto 4) e fasi di amplificazione PCR (protocollo fase 7).

Reagente	Volume	Commento
10x FastStart High Fidelity Reaction Buffer with 18 mM cloruro di magnesio	20 l
10 mM dNTP Mix Solution	5 ml
25 micron PCR Primer PE 1.0	4 pl	Vedere Tabella 4
25 micron PCR Primer PE 2.0	4 pl	Vedere Tabella 4
FastStart High Fidelity Enzyme	2 pl	5 unità / ml FastStart Taq DNA polimerasi
Acqua priva di DNase	125 ml

Tabella 5: reagenti reazione PCR nomi reagenti e le quantità usate nella reazione di amplificazione PCR (protocollo punto 7.1)..

Passo Protocol	Reattivo / detail protocollo	Importo DNA Input
		5-10 ng	25 ng	50 ng
1	MspI enzima	1 ml	2 pl	2 pl
	Volume di reazione MspI digest	50	100	100
4	Adattatori in reazione legatura	1 ml	2 pl	3 ml
	Volume di reazione legatura	20 l	25 microlitri	50 microlitri
5	Protocollo selezione Size	Solo gel Manuale	Pipino Prep o gel manuali	Pipino Prep o gel manuali
7	Concentrazione dei primer PCR	25 micron	25 micron	10 mM per 14 cicli; 25 mM per 18 cicli
	Numero di cicli di PCR	18	18	14-18

Tabella 6:. Protocol modifiche passo per ingresso quantità dei materiali che vanno 5-50 ng Diversi passi in tutto il protocollo richiedono modifiche delle quantità di reagenti utilizzati per generare le librerie di alta qualità provenienti da diverse quantità di materie prime. Modifiche alla quantità di reagenti chiave sono inclusi qui. Regolare tampone e volumi d'acqua nelle reazioni di conseguenza.

Chr	Base	Filo	Copertura	freqC	freqT
chr1	10564	R	366	85.52	14.48
chr1	10571	F	423	91.25	8.75
chr1	10542	F	432	91.2	8.8
chr1	10563	F	429	94.64	5.36
chr1	10572	R	366	96.99	3.01
chr1	10590	R	370	88.11	11.89
chr1	10526	R	350	92	8
chr1	10543	R	368	92.93	7.07
chr1	10525	F	433	91.92	8.08
chr1	10497	F	435	88.74	11.26

Tabella 7: dati ERRBS rappresentativi. Dopo l'allineamento dei dati e citosina determinazione metilazione, si ottiene dati coppia di basi Per ciascuna CpG coperto, il protocollo di allineamento come descritto determina la genomica coordinate (colonne: CHR = cromosoma, Base e Strand)., Il tasso di copertura del locus specifico (Coverage ), e il tasso di citosina rilevamento rispetto timidina come percentuale (freqC e freqT rispettivamente).

Numero di biblioteche ERRBS per corsia	Numero medio di allineato unicamente legge	Il numero medio di CPGs coperto	Copertura per CpG media
1	152.231.184 ± 13.189.678	3.183.594 ± 713.547	85 ± 16
2	77.680.837 ± 7.657.058	2674823± 153.494	49 ± 9
3	49.938.156 ± 2.436.865	2.552.186 ± - 76.624	39 ± 2
4	34.457.208 ± 4.441.686	1.814.461 ± 144.339	28 ± 4

Tabella 8:. Parametri rappresentativi di sequenziamento librerie singoli e multiplex ERRBS riportati sono dati per corsia da 51 a ciclo piste sequenziamento singoli-letto: media e deviazioni standard di allineate in modo univoco si legge, il numero di CPGs coperto e la copertura per sito CpG ottenuto da sequenziamento singolo biblioteche ERRBS per corsia (n = 100), due ERRBS biblioteche per corsia (n = 128), tre ERRBS biblioteche per corsia (n = 11), e quattro ERRBS biblioteche per corsia (n = 11).

Discussione

I dati di risoluzione rendimenti base-pair protocollo presentato di citosina metilazione a biologicamente rilevanti regioni genomiche. Il protocollo come scritto è ottimizzato per 50 ng di materiale di partenza, tuttavia, può essere adattato per gestire una gamma di materiali di ingresso (5 ng o più) ^26. Ciò richiederà aggiustamenti di alcuni dei passaggi del protocollo, come visto nella tabella 6. Le librerie ERRBS sono suscettibili di sequenziamento fine associato e ulteriormente la copertura genomica può anche essere realizzato mediante sequenziamento legge più di 51 cicli. Sequenziamento Multiplexed offrirà un protocollo minor costo per campione, tuttavia, questo comporta copertura ridotta per sito CpG rappresentato nei dati (Figura 5 e la Tabella 8), e non produrrà sufficiente profondità di copertura effettuare analisi che richiedono elevata copertura per sito CpG (ad esempio come descritto da Landan et al. ^33). Infine, questo protocollo (o qualsiasi protoco basata bisolfito-l) non può distinguere tra metil-citosina e idrossimetil-citosina ^46,47. Tuttavia, i dati generati può essere integrato con altri protocollo risultati ^48,49 per delineare le diverse modifiche, e altre modifiche citosina recentemente riportato ^50, dovrebbero essere di interesse.

Appariranno librerie di alta qualità come mostrato nella figura 3A-C, e una volta aggregate per sequenziamento produce una traccia, come mostrato nella figura 3G (traccia rossa) rappresenta contributi molari uguali da due frazioni di libreria. Insufficienza preparazione Libreria può derivare da qualsiasi punto durante la procedura. Se il DNA degradato viene elaborato il risultato sarà in librerie che non sono arricchiti in frammenti MspI e quindi nella copertura bassa CpG utilizzando i parametri di sequenziamento descritti in questo protocollo. Se un enzima è non funzionale o inavvertitamente esclusi da una delle reazioni, il protocollo non produrrà la libreria previsto. Se la rea ​​legaturaction è inefficiente, adattatori sono ad una maggiore concentrazione del previsto, e / o la concentrazione di primer utilizzato è un reagente limitante per le fasi finali di amplificazione, può verificarsi insufficienza libreria. Adattatori in eccesso (visto come un picchi a ~ 150 bp a risultati Bioanalyzer, Figura 3D-F) nella libreria sarà anche interferire con il sequenziamento causa del raggruppamento indiscriminata sia della biblioteca e gli adattatori in eccesso. Mentre una tale libreria può sequenziare apparentemente normalmente, una parte significativa della legge sarà sequenze semplicemente adattatore. Se gli adattatori in eccesso si osservano in una libreria, è meglio ripetere la preparazione biblioteca se il materiale è disponibile usando il materiale in ingresso ottimale all'adattatore rapporti quantitativi. Infine, per garantire efficienza di amplificazione PCR delle librerie, le frazioni biblioteca inferiori e superiori sono mantenuti come campioni separati per tutta la conversione bisolfito e fasi di arricchimento PCR. In caso contrario si produce efficienza differenziale dell'amplificazione durante la PCR reazione di frazioni superiori e inferiori (come si vede nella Figura 3G traccia blu) e il potenziale per la rappresentazione disuguale dei rispettivi loci genomici coperto in ogni frazione libreria durante sequenziamento. L'utente può scegliere di includere una PCR quantitativa passo immediatamente dopo la conversione bisolfito per ulteriore titolazione di cicli PCR ottimali necessari per amplificare le librerie generati.

ERRBS protocollo di preparazione biblioteca ha diversi passaggi chiave in cui sono raccomandati reagenti specifici. Nella fase finale di riparazione, l'uso di una miscela dNTP quattro nucleotidi consente end-riparazione di qualsiasi prodotto che non contengono la sporgenza CG, come quelli derivanti da attività enzimatica MspI stella e frammenti di DNA tranciate presenti nel campione di DNA originale. Ciò si traduce in una migliore rappresentazione CpG nei risultati. Al passo legatura è fondamentale utilizzare una ligasi alta concentrazione (2.000.000 unità / ml) e adattatori metilati per garantire che il riatti legaturaon è efficiente e che la conversione bisolfito non influenza le sequenze essenziali per l'allineamento dei dati accurati adattatore. Nella fase PCR, utilizzando una polimerasi grado di amplificare ricchi GC frammenti di DNA bisolfito trattati è necessaria per alta specificità. Infine, per garantire l'eliminazione di adattatori eccesso e primer, SPRI purificazione tallone (per esempio: Agencourt AMPure XP) è consigliato anziché test basati colonna per legatura ed isolamenti prodotto PCR.

Al fine di generare dati di alta qualità, è importante assicurare la conversione bisolfito efficiente. Il controllo ha presentato offre all'utente la possibilità di determinare l'efficienza di conversione prima di sequenziamento. In alternativa, un DNA non umano come DNA lambda può essere utilizzato come controllo interno (picco-in). A causa delle differenze nelle specie, questo tipo di controllo può essere direttamente incluso nella sequenza a valle (es usati per Yu et al. ^34). Tuttavia, se l'i-picco ins utilizzato, non può essere utilizzato per determinare l'efficienza di conversione prima sequenziamento libreria a meno amplificato modo separato ed indipendente sequenziato prima sequenziamento libreria. I tassi di conversione determinati sono basate sullo stato di metilazione in siti non-CpG. Questo può non essere adatto per l'uso nel contesto di alta metilazione citosina in non-CpG contesto (ad esempio cellule staminali embrionali) e campioni paralleli o altri mezzi per valutare l'efficienza di conversione può essere utilizzato per questo scopo.

Ci sono alcune avvertenze da affrontare che sono unici per il sequenziamento del biblioteche ERRBS. I primi tre basi delle frazioni di libreria sequenziate sono quasi uniformemente non casuale dovuta al sito di riconoscimento MspI taglio (C ^ CGG; vedi Figura 4B, C). Il risultato è la possibilità di una significativa perdita di dati a causa di bassa qualità derivanti da letture di localizzazione grappolo scarsa nonostante apparente ad alta densità di cluster durante la sequenza. Per superare questo ostacolo,includere una libreria di complessità in una corsia indipendente (controllo Phix o altro tipo di biblioteca) come una corsia di controllo dedicato. Biblioteche elevata complessità hanno estremità che contengono una rappresentanza equilibrata di A, C, T e G nelle prime quattro basi in sequenza. Corsie di controllo adatti comprendono biblioteche, come RNA-Seq, ChIP-seq, sequenziamento dell'intero genoma, o un controllo offerto dal costruttore della macchina sequenza (ad esempio Phix controllo v3). Quando designato come una corsia di controllo del rispettivo corsa sequenziamento, può servire come base per la generazione della matrice che viene utilizzata durante i primi quattro basi di sequenziamento per rilevare posizioni cluster. La qualità superiore si legge catturato aumenterà la copertura media per sito CpG 5,2 (n = 4). In alternativa, questa difficoltà tecnica può anche essere superato utilizzando un approccio di sequenziamento scuro come descritto in precedenza ^23. Altri criteri di sequenziamento seguire le procedure operative standard per i protocolli del produttore. Infine, la copertura per CpG chosen per l'analisi dei dati sarà guidato dall'utente e in parte dalle domande biologiche di interesse. 10x soglia di copertura offre un approccio di analisi ad alta copertura, ma questa soglia può essere abbassata che dovrebbe essere di interesse.

Una discussione completa di analisi dei dati ERRBS è oltre la portata di questo articolo, tuttavia, differenziale cytosines e regioni denaturato possono essere determinate utilizzando strumenti open source ^31,51-53. Considerazioni ulteriori analisi e approcci sono stati ben descritto ^54,55, e il lettore è incoraggiato a cercare letteratura per strumenti più adeguati per l'analisi pianificata.

Rispetto ad altri metodi pubblicati, ERRBS offre un protocollo di quattro giorni che, quando eseguita come descritto rendimenti elevati tassi di riproducibilità. E 'stato convalidato rispetto al gold standard di MassARRAY EpiTYPER ^26, è conveniente per i dati di copertura alte, ed è adattabile per vari materiale in entrataimporti (favorevoli per l'elaborazione del campione clinico e altri tipi di cellule a bassa frequenza) e approcci di sequenziamento. Offre risoluzione di base-pair a loci biologicamente rilevanti e può essere utilizzato in integrativo analisi con altre tecniche di profilazione genome-wide fattore di trascrizione vincolante, il rimodellamento della cromatina, segni epigenetici e altre modifiche citosina di interesse. ERRBS uso dei dati in tali studi possono contribuire ad un approccio molecolare completa e consentire per alta dimensionale analisi nello studio dei modelli biologici e le malattie umane.

Materiali

Elenco dei materiali utilizzati in questo articolo

Nome	Azienda	Numero di catalogo	Commenti
MspI	New England Biolabs	R0106M	100.000 unità/ml
NEBuffer 2	New England Biolabs	B7002S Tampone	di reazione per l'enzima MspI; fase del protocollo 1.2
Soluzione fenolica	Sigma-Aldrich	P4557	Equilibrato con 10  mM Tris HCl, pH  8.0; vedere le istruzioni di sicurezza e manipolazione su http://www.sigmaaldrich.com/catalog/product/sigma/p4557
Cloroformio	Sigma-Aldrich	C2432	Vedere le istruzioni di sicurezza e manipolazione su http://www.sigmaaldrich.com/catalog/product/sial/c2432
Glicogeno	Sigma-Aldrich	G1767	19-22 mg/ml
NaOAc	Sigma-Aldrich	S7899	3 M, pH 5.2
Etanolo	Sigma-Aldrich	E7023	200 proof, per biologia molecolare
Tampone EB	Qiagen	19086	10 mM Tris-Cl, pH 8,5
tris(idrossimetil)amminometano (Tris)	Sigma-Aldrich	T1503	preparare una soluzione da 1 M, pH 8,5
Acido tris-etilendiamminotetraacetico (TE)	Sigma-Aldrich	T9285	Diluire in 1x soluzione tampone secondo le raccomandazioni
del produttoreTampone di reazione DNA ligasi T4	New England Biolabs	B0202S	Concentrazione 10x
Miscela di soluzioni di deossinucleotide trifosfato (dNTP) New	England Biolabs	N0447L	10 mM per nucleotide
T4 DNA polimerasi	New England Biolabs	M0203L	3.000 unità/ml
DNA polimerasi I, frammento grande (Klenow)	New England Biolabs	M0210L	5.000 unità/ml
T4 Polinucleotide chinasi	New England Biolabs	M0201L	10.000 unità/ml
Kit di purificazione QIAquick PCR	Qiagen	28104	Utilizzato per la purificazione del prodotto DNA nella fase del protocollo 2.3
2'-deossiadenosina 5'-trifosfato (dATP)	Promega	U1201	100 mM
Frammento di Klenow (3'→ 5' exo-)	New England Biolabs	M0212L	5.000 unità/ml
Kit di purificazione MinElute PCR	Qiagen	28004	Utilizzato per la purificazione del prodotto DNA nella fase del protocollo 3.3
T4 DNA Ligasi	New England Biolabs	M0202M	2.000.000 unità/ml
Adattatore di metilazione Oligo Kit	Illumina	ME-100-0010
Agencourt AMPure XP	Beckman Coulter	A63881	Utilizzato nelle sezioni di protocollo che implementano passaggi di purificazione con biglie magnetiche (passaggi 4.3 e 8.2). Equilibrare a temperatura ambiente prima dell'uso.
Cassette Pippin Prep Gel, 2% agarosio, senza coloranti	Sage Science	CDF2010	con standard interni
Certificato Low Range Ultra Agarose	Bio-Rad	161-3106
Tris-Borato-EDTA (TBE) tampone	Sigma-Aldrich	T4415
Soluzione di bromuro di etidio	Sigma-Aldrich	E1510	10 mg/ml
50 bp DNA Ladder	NEB	N3236S
100 bp DNA Ladder	NEB	N3231S
Gel Caricamento Colorante, Arancione (6x)	NEB	B7022S
Bisturi  Lama  No. 11	Fisher Scientific 3120030
Kit di estrazione del gel QIAquick	Qiagen	28704
Kit di metilazione del DNA EZ	Zymo Research	D5001	Utilizzato nel protocollo passo 6.2
Kit di metilazione-fulmine del DNA EZ	Zymo Research	D5030	Alternativa al passaggio 6.2
Standard universale del DNA umano metilato Zymo	Research	D5011	Utilizzato come controllo della conversione del bisolfito
Sistema PCR ad alta fedeltà	FastStart Kit	di
Qubit dsDNA Life	Technologies	Q32854	Un test di quantificazione del DNA basato sulla fluorescenza; utilizzato nelle fasi del protocollo 1.1, 9.1 e 10.1
DynaMag-2 Magnet	Life Technologies	12321D
	Agilent Technologies	5067-4626
2100 Bioanalizzatore	Agilent Technologies
PhiX Control v3	Illumina	FC-110-3001
HiSeq 2500	Illumina
Pippin Prep	Sage Science
Qubit 2.0 Fluorimetro	Life Technologies	Q32872
Kit cluster TruSeq SR v3-cBot-HS	Illumina	GD-401-3001
Kit TruSeq SBS v3-HS	Illumina	FC-401-3002
TruSeq RNA Preparazione del campione	Illumina	RS-122-2001	Adattatori con codice a barre utilizzati per librerie multiplexing; Vedere File supplementare per il protocollo di multiplexing.
Microcentrifuga
Vortex Mixer
Dry Block Heater
Thermal Cycler
Water Bath
Gel electrophoresis system
Alimentazione
Gel doc
			UV o a luce blu