$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
L'esempio di immagini generate con il reverse-transfezione siRNA protocollo di screening sono stati preparati e analizzati con il software Profiler cellulare. La risultante dati grezzi numerica è tale che ogni cellula è rappresentato individualmente, riconducibile a sua immagine e ben di origine e misurato per diversi parametri di intensità di fluorescenza (Figura 6A). Per ogni cella individuata l'intensità di fluorescenza nucleare media per la RB1 anticorpi P-S780 e l'intensità del DNA integrato per le maschere del DNA nucleare dye-definiti sono determinati. Valori di intensità GFP media per nucleo e nel citoplasma regioni di ogni cella sono anche registrati consentendo il calcolo del nucleare rispetto a fluorescenza citoplasmatica del reporter GFP-CDK2. A valle di questi algoritmi di fluorescenza misure di intensità si fa uso di questi dati delle celle singole per definire porte per due saggi, colorazione anticorpo nucleare e GFP-CDK2 giornalista. Annotazione successiva delle celle the base del risultato del test e l'uso di queste etichette per consentire sottopopolazioni specifiche per essere ulteriormente caratterizzato da una terza misura (contenuto di DNA nucleare) è descritto.
Istogramma dei dati di intensità di fluorescenza grezzi raccolti per ogni test sono un modo efficace di valutare come sottopopolazioni di cellule si comportano in condizioni diverse. Gli istogrammi in figura 6B mostrano le distribuzioni di popolazione di dati delle singole celle di pozzi tripla per ogni condizione di atterramento RNAi. A sinistra sono i dati di intensità anticorpo nucleare e sulla destra sono i dati corrispondenti per il reporter GFP-CDK2. I dati anticorpi RB1 P-S780 rivela che esistono sostanzialmente le cellule in due popolazioni riguardo a questo modificazione post-traduzionale e che popolazioni di cellule con perdita di RB1 fosforilata su S780 possono essere distinti come un picco sinistro intensità nucleare che si arricchisce quando CDK6 viene abbattuto da siRNA. Questo stesso picco sinistrasi vede quando RB1 è di per sé l'obiettivo RNAi, che riflette la rimozione definitiva della proteina e quindi P-S780 RB1 colorazione. Al contrario, le stesse condizioni sperimentali per le stesse cellule, quando osservato attraverso il saggio giornalista GFP-CDK2, mostrano una dinamica differente nei dati delle singole celle. Una distribuzione continua si osserva, con un unico picco, ma siRNA che disturba il ciclo cellulare (siCDK6) e provoca accumulo nei risultati di fase G1 in un'estensione della spalla destra di tale distribuzione (cioè indicando una maggiore presenza di cellule che mostrano un aumento del rapporto GFP nucleare / citoplasma, rappresentato sullo asse X).
Anche mostrato gli istogrammi della figura 6B sono i valori di gate (barre verticali) che vengono scelti sulla base delle distribuzioni di entrambe le serie di dati del saggio. La regola utilizzata per i dati anticorpi RB1 P-S780 è quello di definire la posizione del cancello come mezza altezza: posizione di massima larghezza sulla spalla sinistra del principale(A destra) di picco quando si considerano i dati della cella di controllo negativi (non-targeting siRNA). I dati evidenziati in rosso sono cellule con ridotta e assente RB1 P-S780, che si identifica con questa porta. Un cancello simile posizionato sulla spalla opposta della distribuzione valore del rapporto è utilizzato per il reporter GFP-CDK2. I risultanti cellule sottopopolazione alta rapporto, che mancano o funzionalità ridotta attività CDK2, vengono visualizzati in verde. Per illustrare l'analisi multiplex di due saggi figura 6C mostra l'implementazione di entrambi i valori di gate utilizzando lo script perl 2_gate_classifier.pl per convertire i dati grezzi (Figura 6A) nel file annotato sotto. Questo nuovo file contiene i dati originali accanto una nuova colonna di etichette di classe per ogni cella ei due valori di gate utilizzati per distinguerli (in questo caso porte di 0,004 per i dati anticorpi e 1.5 per il reporter GFP-CDK2 sono stati utilizzati, rispettivamente) .
Dopo aver classificato le singole celle From ciascuna condizione atterramento sulla base dei due test è ora possibile utilizzare queste etichette di classe per aiutare l'annotazione dei grafici dei dati del test. Figura 7 mostra disperdono grafici dei dati delle celle individuali per il P-S780 RB1 e GFP saggi CDK2 Dall'esempio di dati per tutte e tre le condizioni RNAi. Numeri annotazione quadranti sui dispersione mostrano le percentuali relative di ciascuna sottopopolazione gated all'intera per quel contesto knockdown e si generano in R utilizzando le etichette di classe sopra descritti. Queste trame rivelano che, rispetto alle cellule trasfettate con non-targeting siRNA (Figura 7B), cellule trasfettate con siCDK6 rivelare una rete di distribuzione di dati sia spostata verso il basso sulla Y (che indica assenza di RB1 fosforilazione a serina 780) ea destra sulla asse X (che indica bassa attività CDK2, Figura 7C). Entrambi questi cambiamenti sono attesi per atterramento di questo obiettivo. In contrasto con questo, i dati di SirCellule trasfettate B1 (Figura 7A) mostra una perdita della colorazione anticorpale in linea con perdita dell'epitopo, ma poco effetto nella distribuzione dei dati per il reporter CDK2 rispetto ai controlli trasfettate con non-targeting siRNA, suggerendo nessun grande effetto sulla GFP giornalista -CDK2 nasce da RB1 atterramento.
Per esplorare ulteriormente l'uso dei dati di singole cellule, la classificazione sottopopolazione e test multiplexing Figura 8 mostra la rappresentazione della dispersione dei dati siCDK6 dei profili di istogramma Figura 7C accanto appaiati per intensità DNA integrato. Le coppie di istogrammi riguardano opposte metà di tutta la popolazione, diviso sulla base di una intensità anticorpo (destro della dispersione) o GFP-CDK2 valori del rapporto Reporter (sopra la dispersione). La quantificazione dell'intensità DNA nucleare per queste popolazioni mostra due picchi caratteristici di 2N e 4N contenuto di DNA come picchi sinistro e destro, rispettivamente. La intentions delle porte illustrate nelle figure 6, 7 e 8 sono tali che le cellule identificate partire per P-S780 RB1 (contrassegnato: P-S780-) o con un valore elevato rapporto dal reporter GFP-CDK2 (etichettati: G1) sarà essere in fase G1 del ciclo cellulare. In effetti, gli istogrammi profilo DNA per sottopopolazioni identificate con uno di questi saggi prevalentemente contengono cellule con contenuto di DNA 2N. Profili DNA della popolazione opposta gated (etichettati: P-S780 + o Non-G1) contiene cellule con le distribuzioni che vanno da 2N a 4N, in linea con tali cellule adottando una serie di ciclo cellulare posizioni fase post-G1.
Sebbene il fuoco qui è la generazione e l'analisi dei dati delle singole celle da immagini macchiate modo fluorescente, è anche utile essere in grado di prendere questi dati e riassumere ciascun saggio su base ben by-ben-per monitorare la variabilità tra le repliche e le prestazioni di tutti i pozzetti per un esperimento di un intero piatto di dati. a) i dati RB1 P-S780 e b) i dati giornalista GFP-CDK2. I valori tracciati in A e B sono prodotti da due ulteriori script Perl forniti con questo manoscritto; 'Antibody_fluorescence_summary.pl' e 'G1assay_summary.pl', rispettivamente. Questi script utilizzano i dati grezzi creati da Cell Profiler (Nuclei.csv) e dati del rapporto per così come i) cellule totali misurate per pozzetto, ii) il numero di cellule all'interno del cancello, iii) le cellule per cento all'interno del cancello e iv) l'aritmetica media dei misurati, dati grezzi per bene. Questo è incluso come opzione adatta per guardare attraverso grandi insiemi di dati di analisi, prima di concentrarsi sui dati di trattamento individuali utilizzando la valutazione multiplex dei dati delle singole celle come illustrato in 8. I grafici visualizzati qui plot 'iii) le cellule per cento entro il cancello' per entrambi i test, che si adattano alle distribuzioni di dati non normali visti per i dati P-S780 RB1 e GFP-CDK2 nei istogrammi in figura 6B. Questi script anche calcolare 'iv) la media aritmetica dei misurati, dati grezzi per quel bene', che sarebbe adatto analisi dei dati per le risposte popolazione omogenea e la distribuzione dei dati normale, prima e dopo la perturbazione sperimentale.

Figura 1:. Panoramica delle fasi del flusso di lavoro per analizzare quantitativamente i dati di immagine al microscopio fluorescente Il flusso di lavoro è qui rappresentato come quattro fasi (A) In primo luogo, è necessario preparare sperimentalmente cellule per l'imaging a fluorescenza.. L'esempio qui descritto è quella di unschermo in cui aderenti cellule tumorali umane siRNA trattati vengono coltivate per 48 ore, fissate e colorate su una piastra di coltura tissutale a 96 pozzetti. Sono presenti in triplice copia in pozzetti separati all'interno della piastra diverse condizioni RNAi. Le cellule sono colorate con un colorante DNA, un anticorpo specifico per RB1 fosforilata sul CDK4 e 6 selettivo sito di destinazione serina-780 (P-S780 RB1) e sono anche stabilmente esprimono un GFP-CDK2-giornalista, riportando l'uscita del ciclo cellulare G1. Collettivamente queste sonde fluorescenti costituiscono due saggi valutati separatamente all'interno del flusso di lavoro. (B) immagini al microscopio parallelo per ciascuna sonda fluorescente (canale) vengono generati e denominate tale da includere i dettagli con cui il software di analisi dell'immagine può organizzare i dati. (C) La file di immagine vengono caricati nel software Profiler Cell, che identifica algoritmicamente singole cellule e le coppie associati di nuclei e citoplasma prima di cedere misure di intensità per le tre sonde fluorescenti rilevanoed in ciascuno di essi. (D) Infine, uno script Perl viene utilizzato per organizzare i dati quantitativi prime prodotte. Questo passaggio si applica cancelli i dati di intensità di fluorescenza per ogni cella, in modo efficace binning le cellule in sottopopolazioni, che può essere inserito, rintracciato e interrogato. Clicca qui per vedere una versione più grande di questa figura.

Figura 2: I dati sperimentali a essere ottenuta tramite analisi di immagine I, siRNA-trattati, cellule fluorescente fissate dai dati di esempio sono stati impostati ripreso e misure di intensità corrispondenti prese per cella.. Dati di immagine rappresentativi sono riportati per ciascun parametro registrata durante l'analisi dell'immagine intensità (A) DNA nucleare:. L'intensità della colorazione di colorante DNA nucleare viene utilizzato per produrre una misura del DNA per nucleo (B) Intensità nucleare di fosfo-RB1:.-Immuno colorazione specifica per P-S780 RB1 con un (nero) anticorpo primario e fluorescente contrassegnati anticorpo secondario (rosso) consentire una misurazione dell'intensità RB1 fosforilazione a . S780 per nucleo (C) GFP-CDK2 giornalista: Le cellule utilizzate stabilmente esprimere una proteina reporter GFP-tagged che trasloca tra il nucleo e il citoplasma in un modello di serie con il ciclo cellulare. Doppia misurazione dell'intensità GFP nucleare e citoplasmatica accoppiati per ogni cella permette di calcolare un rapporto per cellula che può essere utilizzata per distinguere fase G1 dal resto del ciclo cellulare. Tre obiettivi siRNA saranno utilizzati per illustrare l'analisi; un non-targeting siRNA controllo negativo; CDK6 siRNA come controllo positivo in perturbando RB1 fosforilazione e il progresso del ciclo cellulare; RB1 siRNA per stabilire specificità anticorpale.k "> Clicca qui per vedere una versione più grande di questa figura.

Figura 3: Organizzazione di file di immagine prima dell'analisi immagini Le immagini prelevate dalla piastra di coltura tissutale sono denominati sistematicamente per consentire al software di analisi di immagine per correlare i dati di immagine sul contesto sperimentale originale.. Questa informazione è collocato all'interno del nome del file per ogni immagine. (A) In ciascun pozzetto della piastra esperimento può corrispondere a diversi target RNAi o trattamenti, la parte forme di indirizzo e del nome del file. (B) Il numero di telaio è parte del nome del file . come ogni pozzetto è ripreso a raccogliere più, non si sovrappongono frames (C) sonde fluorescenti da ciascun frame vengono esposte separatamente; di conseguenza, i nomi dei file devono anche riflettere il canale ogni immagine si riferisce a. Cliccate qui per vedere una versione più grande di questa figura.

Figura 4:. L'uso del cellulare Profiler per misurare il DNA nucleare e anticorpi colorazione Con le impostazioni nel file fornito gasdotto (3_channels_pipeline.cppipe), il profiler cellulare misure di software di analisi delle immagini di fluorescenza valori di intensità per il DNA nucleare e anticorpi vincolante relativo alle singole celle . (A) I nuclei sono identificati il 'blu' immagine del canale di DNA tinto. (B) < / Strong> Le posizioni dei DNA macchiato nuclei sono tenuti temporaneamente in un 'maschera Nuclei'. La maschera Nuclei viene quindi sovrapposto (C) le immagini blu e rosso di canale (DNA e anticorpi fluorescenza dati, rispettivamente) ed i valori di fluorescenza di segmenti di immagine che si sovrappongono con la maschera sono registrati contro ogni cella identificata. Individuazione di successo di nuclei vicini, separati può essere valutato visivamente l'aspetto della maschera Nuclei. Per l'illustrazione, mostrata cerchiato in questa immagine maschera, sono esempi in cui le impostazioni scelte per l'algoritmo hanno mis-identificati nuclei vicini come un unico nucleo. Regolazione delle impostazioni algoritmo per ridurre al minimo questi eventi viene introdotti nella sezione di discussione. Cliccate qui per vedere una versione più grande di questa figura.
/ftp_upload/51882/51882fig5highres.jpg "/>
Figura 5: Uso di Profiler cellulare per misurare le intensità GFP nucleari e citoplasmatici Il CDK2 giornalista GFP-tag trasloca tra il nucleo e il citoplasma in relazione alla posizione del ciclo cellulare delle cellule.. Allo stesso tempo, che Profiler cellulare calcola i DNA e anticorpi intensità nucleari per cella (Figura 4), calcola anche il nucleare rapporto citoplasma di intensità GFP per ogni cella. (A) I dati colorante DNA per ogni immagine viene utilizzato per generare una maschera Nuclei. Profiler (B) Cell utilizza la maschera Nuclei congiuntamente dell'immagine GFP dal reporter GFP-CDK2 per inizializzare la posizione di ogni cella con e poi si espande il perimetro di ogni cella per stimare l'intera impronta di ciascuna cella. Questo diventa un nuovo, 'maschera cellulare'. (C) La maschera nuclei viene sottratto dalla maschera cella di cedere una serie ciambella-come citoplasma contorni, che diventanola 'maschera citoplasma'. (D) La maschera e il citoplasma maschera nuclei sono utilizzati da Profiler cellulare per misurare coppie di valori GFP nucleari e citoplasmatici. Questi valori associati vengono quindi utilizzati per Profiler cellulare per calcolare indici, che informano per la posizione di ogni cella nel ciclo cellulare. Cliccate qui per vedere una versione più grande di questa figura.

Figura 6:. L'estrazione dei dati - Elaborazione dati delle singole celle grezzi imponendo porte su valori di analisi delle tendenze biologici dai dati delle celle individuali per la colorazione di anticorpi e saggi giornalista GFP-CDK2 vengono estratti utilizzando dati gated. Istogrammi dei dati grezzi consentono l'identificazione dei valori cancello adatti. Questi sono poi imposti con uno script Perl. (A) Laprodotto finale di analizzare i file di immagini con le impostazioni fornite di Profiler cellule sono separate da virgole valore (.csv) file. Questi file c ontain dati delle celle individuali relativi a ciascuno dei diversi segmenti sub-cellulare. Il file 'Nuclei.csv' contiene tutte le misure selezionate relative all'uso della maschera Nuclei. Queste misure comprendono intensità anticorpo nucleare, intensità DNA nucleare e il rapporto di GFP (nucleo / citoplasma). (B) istogrammi dell'intensità anticorpo nucleare (a sinistra) e rapporti di giornalista GFP-CDK2 (destra) tracciata dai dati delle celle individuali per ogni condizione di atterramento siRNA . Le barre sugli istogrammi visualizzati mostrano le posizioni del cancello desiderati per questi test. Dati colorate su istogrammi indicano le sottopopolazioni gated. (C) Le porte per i due saggi illustrati in B vengono applicati ai dati grezzi usando lo script Perl '2_gate_classifier.pl'. Lo scriptcrea una copia modificata del file di output Profiler cellulare originale (Nuclei.csv) assistere successiva tracciato. I due valori cancello vengono registrati nel nuovo file (evidenziato in colore qui) e viene aggiunta una nuova colonna 'Label'. La etichette bin ciascuna cella in una delle quattro possibili sottogruppi basati sui due valori di analisi gated per ogni cella. Queste etichette sono utilizzati in successivi lotti che presentano calcoli dei contributi di ogni sottopopolazione, nonché il controllo incrociato dei parametri aggiuntivi generati in Profiler cellulare. Clicca qui per vedere una versione più grande di questa figura.

Figura 7:. Diagrammi a dispersione per ogni condizione di siRNA raffigurante i dati grezzi per le singole cellule e le posizioni del cancello Scatter appezzamenti di individati doppi cellulari da tutte le immagini per le condizioni siRNA indicati: (A) siRB1; (B) controllo negativo-Sinon mira; (C) siCDK6. Tramato contro le assi Y sono valori di fluorescenza nucleare da anti-P-S780 RB1 colorazione. Tramato contro le assi X sono i valori del rapporto corrispondenti calcolati dal giornalista GFP-CDK2. Le barre rosse e verdi indicano le posizioni delle porte per il cancello RB1 P-S780 e le porte del reporter GFP-CDK2, rispettivamente. Le due porte dividono le cellule in quattro sottopopolazioni ei numeri sui quadranti risultanti sono il numero percentuale di cellule da ciascuna di esse. Annotazioni intorno agli assi per A indicano le quattro possibili etichette elementi applicati a ciascuna cella dal 2_gate_classifier.pl script Perl. Queste etichette sono mostrati in relazione al rispettivo cancello dosaggio e sono utilizzati in R-script (analysis.r) per generare le trame in figure 6, 7 e Cliccate qui per vedere una versione più grande di questa figura.

Figura 8: Cellulari sottopopolazioni definite dai due saggi di transito G1 mostrano 2N e 4N DNA profili in linea con il risultato del test La rappresentazione della dispersione dei dati per le celle siCDK6 si ripete dalla figura 7C.. Intorno della dispersione sono istogrammi per intensità DNA nucleare integrata di sottogruppi della popolazione. Quelli sopra la dispersione riguardano il saggio giornalista GFP-CDK2. Quelli a destra della dispersione riguardano nucleare fosfo-RB1 solo misure di anticorpi. Le linee di gate colorati vengono estesi per mostrare il loro rapporto con gli istogrammi. Sono indicate anche le etichette porta per la quale sono stati selezionati i dati della cella per questi complotti supplementari. Le cellule con perdita di RB1 fosforilati su serina 780 (P-S780-) o quelli con un alto GFP-CDK2 giornalista nucleare rapporto citoplasmatica (indicanti un'attività basso CDK2) mostrano profili DNA prevalentemente 2N-like, mentre i loro omologhi opposte per ogni rispettiva analisi mostrano una distribuzione di 2N e 4N, caratteristica di un post-G1 popolazione fase mista di cellule. Cliccate qui per vedere una versione più grande di questa figura.

Figura 9: Riepilogo appezzamenti di valori di analisi gated per ogni siRNA condizione trame dati di sintesi del gated (A) P-S780 dati RB1 e dati (B) GFP-CDK2 da pozzi tripla per ogni condizione di knockdown siRNA.. I valori sono stati calcolati dalla uscita Profiler cellulare grezzo (Nuclei.csv) con ilScript Perl, 'antibody_fluorescence_summary.pl' (A) o 'G1assay_summary.pl' (B). I valori tracciati sono mezzi per la percentuale delle cellule all'interno del cancello applicato a ciascun test. Bar indicano errori standard calcolati da pozzi triplice copia. Spaiati, homoscedastic valori T-Test P per ogni condizione di atterramento rispetto ai non-targeting siRNA sono riportati sopra i dati tracciati dove P <0.001 (**) e P <0,05 (*). Cliccate qui per visualizzare una versione più grande di questo figura.
Figura S1. Installazione di software Profiler cellulare per l'analisi delle immagini. (A) Schermata di Profiler cellulare prima di qualsiasi impostazione di analisi delle immagini sono inseriti. (B) Schermata di Profiler cellulare dopo i particolari algoritmi contenuti in '3_channels_pipeline.cppipe' sono stati caricati. L'alto scheda illuminato in alto a sinistra indica che questa schermata mostra i parametri per la fase 'LoadImages' dell'analisi. Cliccando sulle altre parti della lista di sotto di questo rivelerà i dettagli per le fasi successive di analisi. (C) Schermata di Profiler cellulare con dettagli per cartella di input e la cartella di output è entrato. (D) Schermata di Profiler cellulare dopo il 'Analyze Pulsante immagini "è stato cliccato per iniziare l'analisi. Sovrapposto sono tre nuove finestre che illustrano le maschere algoritmicamente prodotte generati dal software dalle immagini sotto analisi. Queste finestre sono accessibili facendo clic sulle icone 'occhio' alla posizione aperta accanto ai passi importanti nell'analisi, nell'angolo in alto a sinistra della finestra principale Profiler cellulare. Questi punti di vista aiutano l'utente a verificare se le impostazioni che generano le maschere coriandoli color d'accordo con gli originali, in scala di grigi, i dati di accompagnamento.
ent ">
Figura S2. L'uso di Perl e rstudio al cancello dati delle celle individuali e tracciare le sottopopolazioni di cellule risultanti. (A) Il pannello di destra mostra la cartella scelto di ricevere i file di output .csv (icone verdi) dall'analisi Profiler cella. Gli script Perl forniti con il manoscritto (icone blu) vengono copiati in questa cartella. Evidenziato è lo script Perl
'2_gate_classifier.pl', che è stato fatto doppio clic con il mouse per produrre la finestra di dialogo nel pannello di sinistra. Sono indicate le prompt e corrispondenti digitato le risposte necessarie per cancello i dati delle celle singole dal file
'Nuclei.csv'. (B) Schermata di rstudio subito dopo aver caricato lo script
'analysis.R'. Evidenziati sono i comandi per caricare i dati gated da
A in il software prima del tracciato (i dettagli della nota nelle linee 5 e 6 dovrà essere regolato in base alla quale i dati gated si trova sul calcolor utilizzato per l'analisi).
(C) Schermata di rstudio una volta i dati sono stati caricati.
(D) Schermata di rstudio mostra ha evidenziato il blocco di codice necessaria per produrre il grafico mostrato nella finestra in basso a destra. Codici per ciascun lotto sono separati da righe vuote e raggruppati per tipo di trama.
| bersaglio siRNA | Indirizzi pozzetti |
| Non-targeting (NT) | E5, F5, G5 |
| Retinoblastoma (RB) | E7, F7, G7 |
| Ciclina chinasi dipendente 6 (CDK6) | B2, C2, D2 |
Tabella 1: Well indirizzi e corrispondenti condizioni siRNA utilizzato nell'esempio set di dati.