$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
La disponibilità di sequenziamento profondo veloce, affidabile e conveniente ha rivoluzionato molti aspetti nel campo delle scienze della vita, consentendo grande profondità nelle analisi di sequenziamento. Un restante sfida il design innovativo e la creazione del rappresentante librerie di sequenziamento. Qui descriviamo un protocollo per catturare molecole di cDNA virale nascente, in particolare gli intermedi del processo di trascrizione inversa di HIV-1.
La fase più critica in questa strategia è la legatura di un adattatore per le aperte 3'-estremità in maniera quantitativa e imparziale. Efficienze di legature tra due termini di ssDNA, sia inter- e intramolecolari, sono stati studiati e ottimizzati per vari applicazioni11,26,27,28,29. La scelta di utilizzare un adattatore di tornante con T4 DNA ligasi nelle condizioni descritte al punto 3.3 è il risultato dell'ottimizzazione empirica in cui abbiamo valutato diverse ligasi, adattatori e reagenti per la legatura dei oligonucleotides sintetici che rappresentano Sequenze di HIV-1 (tabella 2) (dati non mostrati). In queste reazioni in vitro test, abbiamo confermato che la ligasi del DNA T4 mediato legatura dell'adattatore tornante, come descritto da Kwok et al. 11, ha una distorsione molto bassa e raggiunge vicino a legatura completa di molecole acceptor quando l'adattatore viene utilizzato in eccesso. L'efficienza di legatura era inalterata tramite l'aggiunta di sequenza nucleotidica per eseguire il rendering l'adattatore compatibile per il sistema multiplex primer (Vedi Figura 4). In confronto, abbiamo scoperto che una ligasi del DNA/RNA di termostabile 5' (Vedi "A-ligasi" Tabella materiali per esatte ligasi rispetto qui), che è una ligasi derivati dal RNA che è stata sviluppata in parte per migliorare l'efficienza di legatura con ssDNA come accettore 27, era infatti più efficace a legare due molecole di ssDNA di ligasi di RNA ("ligasi B"), ma aveva una distorsione significativa, con forti differenze di efficienza di legatura anche tra oligonucleotidi con differenze di lunghezza base singola [tabella 2 ; HTP con metà G (a) e (b)]. Inoltre, abbiamo trovato solo una minima distorsione in reazioni con "Ligasi C" combinato con un adattatore che trasportano un randomizzato 5'-termini (una strategia utilizzata per compensare la polarizzazione del nucleotide noto di "Ligasi C"; Vedi ad esempio Ding et al. 30). Tuttavia, la "ligasi C"-mediate intermolecolari legature erano incompleti, rendendo il sistema ligasi del DNA T4 la scelta migliore.
Diversi passaggi di controllo di qualità sopra il corso del protocollo e l'inclusione di controlli positivi e negativi permettono per la rilevazione dei potenziali problemi prima continuazione di dosaggio e forniscono indicazioni per la risoluzione dei sforzi. Le quantificazioni di qPCR ai punti 2.2.2 e 2.3.12 garantiscono che la quantità del materiale in ingresso è sufficiente. Numeri di copia cDNA tipico nella gamma di eluizione (dal punto 2.1) 200 µ l da circa 10.000 a 300.000 per µ l. Il passo di cattura di ibrido può comportare una perdita di HIV-1 complessivo cDNA quantità ma deve risultare in un forte arricchimento del cDNA di HIV-1 specifico su DNA cellulare, che può essere determinato utilizzando appositi primer per quantificare il DNA di genomic prima e dopo arricchimento di qPCR o misurando la concentrazione di DNA totale. Recuperati cDNA di HIV-1 dopo l'ibrido catturare passi dovrebbe essere almeno il 10% dell'input. Basso materiale di partenza può spiegare altrimenti un controllo positivo del oligonucleotide di successo (Vedi punto 3.3.2), ma solo limitato letture raggiunta nei campioni. Bassa legge numeri complessivi potrebbero essere spiegati anche da sovrastima della concentrazione libreria a causa della presenza di specie di DNA irrilevante senza adattatori MiSeq. Questo si tradurrebbe in densità bassa cluster e può essere migliorato determinando la concentrazione delle sequenze di HIV-1 nella libreria di qPCR oltre all'importo totale del DNA di analisi fluorometriche. A causa della natura altamente sensibile del metodo, si dovrebbe prestare particolare attenzione per evitare la contaminazione anche a basso livello, sia da altri campioni (in particolare, provenienti dalle scorte del oligonucleotide controllo alta concentrazione) così come da attrezzature di laboratorio. Lavorare in un UV sterilizzazione workstation PCR è utile a questo proposito. L'elettroforesi automatizzata della biblioteca finale (punto 6.1.2) è un'ulteriore misura di controllo di qualità. La gamma di dimensioni di acido nucleico in genere osservata è tra 150 e 500 nt. primer che può essere rilevato nel controllo facoltativo dopo la PCR, e prima di purificazione (Vedi nota al punto 5.2) dovrebbe essere assente. In un risultato rappresentativo, la curva di intensità di campione ha un picco circa 160-170 nt e un secondo più nitida picco circa 320 a 350 nt. Ciò probabilmente riflette l'abbondanza superiore spesso visto in sia relativamente breve (da 1 a 20 nt inserto) inversione trascritti e forte-fermata Full-Length (lunghezza di 180-182 nt inserto) (Figura 3b).
Mentre il protocollo presentato e selezionato gli iniettori sono specifici per primi costrutti di trascrizione inversa di HIV-1, il metodo è generalmente applicabile a qualsiasi studio con l'obiettivo di determinare aperti 3'-termini di DNA. Le principali modifiche richieste in altri contesti sarà il metodo per la cattura di ibrido e la strategia di disegno di primer. Ad esempio, se la destinazione è essere adattata alla fine trascrizioni di HIV-1, un numero maggiore di oligonucleotidi biotinilati diversi acquisizione ricottura su tutta la lunghezza del cDNA sarebbe consigliabile e probabilmente diminuirà la perdita nel passaggio cattura ibrida. Come accennato nell'introduzione, è importante considerare le limitazioni quando si progetta l'intervallo nel quale devono essere rilevato per evitare diverse fonti di bias 3'-stazione termini. Primo, ci può essere una distorsione nelle reazioni di PCR se i modelli con l'adattatore sono di lunghezza notevolmente variabile. In secondo luogo, la piattaforma di sequenziamento usata qui (per esempio, MiSeq) ha una gamma di lunghezza comodo inserto per il clustering ottimale, e significativamente più brevi e più a lunghi i prodotti non possono essere sequenziati con la stessa efficienza. In parte, questo può essere affrontato dal punto di vista, come è stato fatto calcolando il fattore di correzione per bias lunghezza lineare (Vedi Figura 4, grafico in basso). Tuttavia, se la regione di cui 3'-termini mappatura sono volute è lunga (> 1000 nt), è più consigliabile suddividere le reazioni con le trascrizioni dei e utilizzare più a Monte degli iniettori per valutare 3'-termini a sezioni.
Il programma di analisi è stato scritto in-House per lo scopo specifico di analizzare sia l'ultimo nucleotide della sequenza HIV-1 adiacente alla sequenza fissa adattatore così come la variazione di base di tutte le basi per identificare eventuali mutazioni. I singoli passaggi comprendono quanto segue: in primo luogo, le sequenze di adattatore vengono tagliate usando il toolkit fastx-0.0.13; quindi, vengono rimossi eventuali sequenze che vengono duplicate (cioè sequenze identiche tra cui il codice a barre). Tutte le rimanenti letture uniche quindi sono allineate alla sequenza di HIV-1 con Papillon (http://bowtie-bio.sourceforge.net/index.shtml) con il disallineamento massimo fissato a tre basi. La sequenza di modello è composto dal primo 635 nt del cDNA di HIV-1 (ceppo NL4.3), che comprende la sequenza - sss e il primo prodotto di trasferimento di filo fino alla pista di polipurinico (U5-R-U3-PPT; Vedi Figura 1). In tal modo, il software fornito e modelli sono direttamente adatti solo se il metodo viene utilizzato per la stessa applicazione (rilevazione delle trascrizioni precoce inversione del HIV-1NL4.3). Le regolazioni dovranno prevedere altre sequenze di destinazione. Le posizioni della 3'-termini per ogni lettura sono state determinate dalla posizione nell'allineamento. Chiamate di base per ogni posizione sono registrate e tassi di mutazione sono calcolati dalla copertura totale di ogni base, che varia, come letture sono di diverse lunghezze e inserti lunghi potrebbero non essere interamente coperti mediante il sequenziamento di 125-base in Read2.
Per concludere, riteniamo che il metodo descritto per essere uno strumento prezioso per molti tipi di studi. Ovvie applicazioni includono le indagini dei meccanismi che l'inibizione di trascrizione d'inversione attraverso farmaci antiretrovirali o fattori di restrizione cellulare. Tuttavia, solo relativamente piccoli aggiustamenti dovrebbero essere necessari adattare il sistema a 3'-termini mapping all'interno di altri singoli filamenti intermedi DNA virale, che sono presenti, ad esempio, nella replica di parvovirus. Inoltre, il principio del metodo, in particolare il suo passo di legatura ottimizzato, possa fornire una parte fondamentale della progettazione di preparazione di librerie per la caratterizzazione di eventuali estensioni di 3'-DNA, tra cui allungamenti catalizzati da cellulare DNA double-stranded polimerasi.