Questo rapporto descrive il protocollo SLIC-CAGE completo per ottenere librerie pronte per il sequenziamento da nanogrammi di materiale RNA totale iniziale (Figura 1). Per ottenere il mix di portatori di RNA sintetici, in primo luogo, i modelli di portante PCR devono essere preparati e purificati in gel per eliminare i prodotti laterali PCR (Figura 2A). Ogni modello PCR (dieci in totale) viene prodotto utilizzando un inoltro comune, ma un'inversione diversa (Tabella 2), che porta a lunghezze diverse del modello PCR per consentire la variabilità delle dimensioni dei portatori di RNA sintetici. Una volta purificati, i modelli PCR vengono utilizzati per la trascrizione in vitro delle molecole portatrici. Un singolo prodotto portarna è previsto se i modelli sono purificati in gel (vedi analisi rappresentativa del gel nella Figura 2B). La preparazione del supporto può essere scalata a seconda della necessità, e quando preparato, mescolato e congelato a -80 gradi centigradi per un uso futuro.
Utilizzando la quantità minima raccomandata di RNA totale campione (10 ng) combinato con 16-18 cicli di amplificazione PCR, è possibile ottenere librerie SLIC-CAGE ad alta complessità. Il numero di cicli PCR necessari per amplificare la libreria finale dipende molto dalla quantità totale di RNA di input utilizzato (il numero previsto di cicli è indicato nella tabella4).
Dopo il primo ciclo di degradazione, nei risultati qPCR (passaggio 17), la differenza prevista tra i valori Ct ottenuti utilizzando adaptor_f1 o carrier_f1 primer è 1-2, con i valori Ct ottenuti con adaptor_f1 inferiore rispetto a carrier_f1.
La distribuzione delle lunghezze dei frammenti nella libreria finale è compresa tra 200-2.000 bp con la dimensione media del frammento di 700-900 bp (basata sull'analisi della regione utilizzando il software Bioanalyzer, Figura 4B,D). Frammenti più brevi, come illustrato nella figura 4A,C, devono essere rimossi da ulteriori cicli di esclusione di dimensioni (passaggi 20-21). Questi brevi frammenti sono artefatti di amplificazione PCR e non la libreria di destinazione. Si noti che i frammenti più brevi si raggruppano meglio sulle celle di sequenziamento e possono causare problemi di sequenziamento.
La quantità prevista di materiale di libreria ottenuto per campione è compresa tra 5-50 ng. Importi significativamente inferiori sono indicativi di perdita di campione durante il protocollo. Se la quantità bassa ottenuta è sufficiente per la sequenza (è necessaria 2-3 ng delle librerie in pool), le librerie potrebbero essere di minore complessità (vedere di seguito).
A seconda della macchina di sequenziamento, potrebbe essere necessario ottimizzare la quantità della libreria caricata nella cella di flusso. Utilizzando un Illumina HiSeq 2500, il caricamento di 8-12 pM librerie SLIC-CAGE dà in media 150-200 milioni di letture, con >80% delle letture che superano il punteggio di qualità Q30 come soglia.
Le letture ottenute vengono quindi mappate al genoma di riferimento [per letture da 50 bp, Bowtie212 può essere utilizzato con parametri predefiniti che consentono zero mancate corrispondenze per sequenza di semi (22 bp)]. L'efficienza di mappatura prevista dipende dalla quantità totale di ingresso dell'RNA e sono presentate nella tabella 5. Le letture mappate in modo univoco possono quindi essere caricate nell'ambiente di elaborazione grafico e statistico R13 ed elaborate utilizzando CAGEr (pacchetto Bioconductor14). La vignetta pacchetto è facile da seguire e spiega il flusso di lavoro e l'elaborazione dei dati mappati in dettaglio. Un facile controllo visivo della complessità della libreria è la distribuzione della larghezza del promotore, in quanto le librerie a bassa complessità avranno promotori artificialmente stretti (Figura 5A, libreria SLIC-CAGE derivata da 1 ng di RNA totale, per i dettagli si vedano pubblicazione10). Tuttavia, anche le librerie SLIC-CAGE a bassa complessità consentono l'identificazione di CTSS reali, con maggiore precisione rispetto ai metodi alternativi per il mapping TSS a basso/medio input (Figura 5B,C).

Figura 1: Passaggi del protocollo SLIC-CAGE. L'RNA campione viene miscelato con la miscela portante dell'RNA per ottenere 5 g di materiale totale di RNA. cDNA è sintetizzato attraverso la trascrizione inversa e il cappuccio viene ossidato utilizzando periodato di sodio. L'ossidazione consente l'attaccamento della biotina al cappuccio utilizzando l'idrato di biotina. La biotina si attacca all'estremità di 3, in quanto viene anche ossidata usando il periodato di sodio. Per eliminare la biotina dagli ibridi mRNA:cDNA con cDNA in completamente sintetizzato e dalle estremità di 3o mRNA, i campioni vengono trattati con RNase I. cDNA che ha raggiunto l'estremità 5' dell'mRNA viene quindi selezionato per purificazione di affinità su perline magnetiche streptavidin ( cap-trappola). Dopo il rilascio del cDNA, i linker da 5 a 3 vengono legati. Le molecole della biblioteca che provengono dal vettore sono degradate usando le endonucleasi i-SceI e I-CeuI e i frammenti vengono rimossi utilizzando perline magnetiche SPRI. La libreria viene quindi amplificata PCR. Fare clic qui per visualizzare una versione più grande di questa figura.

Figura 2: Analisi congiunta gel dei modelli PCR portante e trascrizioni carrier in vitro. (A) Modelli di portaggio PCR prima della purificazione del gel: il primo pozzo contiene il marcatore da 1 kbp, seguito dai modelli PCR portanti 1, 1-10. (B) Trascrizioni carrier in vitro: il primo pozzo contiene il marcatore da 1 kbp, seguito dalle trascrizioni del vettore 1-10. Le trascrizioni dei portante sono state denaturate riscaldando per 5 minuti a 95 gradi centigradi prima del caricamento. Fare clic qui per visualizzare una versione più grande di questa figura.

Figura 3: Traccia rappresentativa della qualità del DNA (chip di DNA ad alta sensibilità) di SLIC-CAGE prima del primo ciclo di degradazione del vettore. Fare clic qui per visualizzare una versione più grande di questa figura.

Figura 4: Tracce rappresentative di qualità del DNA (chip di DNA ad alta sensibilità) delle librerie SLIC-CAGE dopo l'amplificazione PCR. (A) libreria SLIC-CAGE che richiede un'ulteriore selezione delle dimensioni per la rimozione di brevi frammenti. (B) Libreria SLIC-CAGE dopo la selezione delle dimensioni utilizzando 0,6x perline SPRI al rapporto campione. (C) Libreria SLIC-CAGE di quantità di output inferiore che richiede la selezione delle dimensioni per la rimozione del frammento corto. (D) Libreria SLIC-CAGE di quantità di uscita inferiore dopo la selezione delle dimensioni utilizzando 0.6:1 Spri perline al rapporto campione. Fare clic qui per visualizzare una versione più grande di questa figura.

Figura 5: Convalida delle librerie SLIC-CAGE. (A) Distribuzione delle larghezze interquantili del tag nelle librerie SLIC-CAGE preparate a partire da 1, 5 o 10 ng di RNA totale S. cerevisiae, e nella libreria nAnT-iCAGE preparata da 5 g di S. cerevisiaeRNA totale. Una quantità elevata di cluster di tag stretti nella libreria SLIC-CAGE 1 ng ne indica la bassa complessità. (B) Curve ROC per l'identificazione CTSS nelle librerie S. cerevisiae SLIC-CAGE. Tutti i CTS S. cerevisiae nAnT-iCAGE sono stati utilizzati come un vero set. (C) Curve ROC per l'identificazione CTSS nelle librerie s. cerevisiae nanoCAGE. Tutti i CTS S. cerevisiae nAnT-iCAGE sono stati utilizzati come un vero set. Il confronto delle curve ROC mostra che SLIC-CAGE supera fortemente nanoCAGE nell'identificazione CTSS. Sono stati utilizzati i dati di ArrayExpress E-MTAB-6519. Fare clic qui per visualizzare una versione più grande di questa figura.
Tabella 1: Sequenza del gene sintetico portatore. I siti I-SceI sono in grassetto e corsivo in viola, e i siti di riconoscimenti I-CeuI sono verdi. Fare clic qui per scaricare questo file.
| portatore | primer inverso 5'-3' | Lunghezza del prodotto PCR / bp |
| 1 : il nome del | PCR_N6_r1: NNNNNNCTACGTGTCGCAGAGAATT | 1034 |
| 2 Il nome del sistema | PCR_N6_r2: NNNNNNTATCCAGATCGTTGAGCTGC | 966 |
| 3 (COM del nome | PCR_N6_r3: NNNNNNCACTGCGGCTCTCTCT | 889 889 |
| 4 DEL psu' | PCR_N6_r4: NNNNNNGCCGTCGATAACTTGTTGTTGT | 821 (in questo stato del |
| 5 Del numero 3( | PCR_N6_r5: NNNNNNAGTTGACCGCAGAAGTCTCTC | 744 Del 144 |
| 6 È possibile: | PCR_N6_r6: NNNNNNNNGTGAAGAATTCTTCTCTCCA | 676 del sistema |
| 7 (in questo stato | PCR_N6_r7: NNNNNNCTCGCGGCTCTCTCATAAAC | 599 |
| 8 (IN vio | PCR_N6_r8: NNNNNNTATACGCGATGTCGTCGTAC | 531 |
| 9 (in vie | PCR_N6_r9: NNNNNNACCGGCGCCTCTCTGCAGG | 454 |
| 10 del sistema | PCR_N6_r10: NNNNNNCAGGACGTTTTTGCCCAGCA | 386 mila |
| Il primer di inoltro è lo stesso per tutti i modelli di vettori. Sottolineato è la sequenza promotore T7. PCR_GN5_f1: TAATACGACTCACTATAGNNNCAGCGTTCGCT | |
Tabella 2: Primer per l'amplificazione del modello carrier. Forward primer è lo stesso per tutti i modelli di vettore. Sottolineato è la sequenza promotore T7. PCR_GN5_f1: TAATACGACTCACTATAGNNNNNCGCGTTCGCT. Utilizzando diverse primer inverse, vengono prodotti modelli PCR e quindi RNA portante di lunghezza diversa.
| portatore | lunghezza | non ancora | con tappo/g |
| 1 : il nome del | 1034 | 3.96 (in inglese) | 0,45 (in linguaggio da 1,05) |
| 2 Il nome del sistema | 966 | 8.36 (in inglese) | 0,95 (in questo da fwlinkin base al nome) |
| 3 (COM del nome | 889 889 | 4.4 (in questo stato del documento) | 0,5 0,5 |
| 4 DEL psu' | 821 (in questo stato del | 6.6 | 0,75 (in questo 0,05) |
| 5 Del numero 3( | 744 Del 144 | 4.4 (in questo stato del documento) | 0,5 0,5 |
| 6 È possibile: | 676 del sistema | Ore 3.08 | 0,35 (in questo da fwlinka che) |
| 7 (in questo stato | 599 | 4.4 (in questo stato del documento) | 0,5 0,5 |
| 8 (IN vio | 531 | 3.96 (in inglese) | 0,45 (in linguaggio da 1,05) |
| 9 (in vie | 454 | 2.64 | 0.3 0.3 |
| 10 del sistema | 386 mila | 2.2 2 Il suo | 0,25 |
Tabella 3: Mix di portatori di RNA. In totale 49 g del mix portante 0,3-1 kbp: non con tappo : 44 g, con tappo di 5 g.
| Ingresso TOTALE rna /ng | Cicli PCR |
| 1 ng | 18 mi lato |
| 2 ng | 17 mi lato |
| 5 ng | 16 |
| 10 ng | 15-16 anni |
| 25 ng | 14-15 anni |
| 50 ng | 13-15 anni |
| 100 ng | 12-14 anni |
Tabella 4: Numero previsto di cicli PCR in dipendenza dall'input totale di RNA campione. Il numero approssimativo di cicli si basa su esperimenti eseguiti utilizzando Saccharomyces cerevisiae, Drosophila melanogaster e Mus musculus total RNA.
| Ingresso/ng totale di RNA | % mappata complessivamente | % mappata in modo univoco | % vettore |
| 1 ng | 30 milio | 20-30 anni | 30 milio |
| 2 ng | 60 del sistema | 20-50 anni | 10 del sistema |
| 5 ng | 60-70 anni | 40-60 | 5-10 anni |
| 10 ng | 60-70 anni | 40-60 | 5-10 anni |
| 25 ng | 65-80 | 40-70 | Da 0 a 5 |
| 50 ng | 65-80 | 40-70 | 0-3 |
| 100 ng | 70-85 anni | 40-70 | 0-2 |
Tabella 5: Efficienza di mappatura prevista e dipendenza della quantità totale di ingresso dell'RNA. I numeri approssimativi sono presentati e basati su esperimenti eseguiti utilizzando Saccharomyces cerevisiae e Mus musculus total RNA.