$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
L'exemple d'images généré en utilisant l'inverse-transfection siRNA protocole de dépistage ont été préparées et analysées avec le logiciel Cell Profiler. Les données brutes numérique résultant est telle que chaque cellule est représentée individuellement, traçables à son image et bien d'origine et mesurée pour plusieurs paramètres d'intensité de fluorescence (figure 6A). Pour chaque cellule identifiée l'intensité moyenne de fluorescence nucléaire pour l'anticorps RB1 P-S780 et l'intensité de l'ADN intégré pour les masques de protection nucléaire de colorant défini ADN sont déterminés. Les valeurs moyennes d'intensité de la GFP pour les régions de noyau et le cytoplasme de chaque cellule sont également enregistrés permettant le calcul de l'énergie nucléaire par rapport fluorescence cytoplasmique du rapporteur GFP-CDK2. En aval de ces fluorescence algorithmique mesures d'intensité on utilise ces données de cellules individuelles pour définir des grilles pour les deux essais, la coloration d'anticorps nucléaire et CDK2-GFP rapporteur. Annotation ultérieure des cellules sur ee base des résultats dosage et l'utilisation de ces étiquettes pour permettre à des sous-populations spécifiques devant être en outre caractérisé par une troisième mesure (teneur en ADN nucléaire) est décrit.
parcelles histogramme des données d'intensité de fluorescence brutes recueillies pour chaque essai sont un moyen efficace d'évaluer comment les sous-populations de cellules se comportent dans des conditions différentes. Les histogrammes de la figure 6B montrent les distributions de population de données cellulaires individuels de puits en triple pour chaque condition knockdown ARNi. Pour la gauche sont les données pour l'intensité d'anticorps nucléaire et sur le droit sont les données correspondantes pour le rapporteur GFP-CDK2. Les données d'anticorps RB1 P-S780 révèle que les cellules existent largement dans deux populations à l'égard de cette modification post-traductionnelle et que les populations de cellules avec une perte de RB1 phosphorylée sur S780 peuvent être distingués en tant que pic de gauche de l'intensité nucléaire qui est enrichie quand CDK6 est renversé par siRNA. Cette même pic gaucheest vu quand RB1 se est la cible ARNi, reflétant la suppression pure et simple de la protéine et donc P-S780 RB1 coloration. En revanche, les mêmes conditions expérimentales pour les mêmes cellules, lorsque observée par l'intermédiaire du dosage de rapporteur GFP-CDK2, présentent une dynamique différente dans les données de cellules individuelles. Une distribution continue est observée, avec un seul pic, mais siRNA qui perturbe le cycle cellulaire (siCDK6) et provoque l'accumulation dans les résultats de la phase G1 dans une extension de l'épaule droite de cette distribution (ce est à dire indiquant la présence accrue des cellules présentant une augmentation du ratio de la GFP nucléaire / cytoplasme, tracé sur l'axe X).
On voit également sur les histogrammes de la figure 6B sont des valeurs de grille (barres verticales) qui sont choisis sur la base des distributions de deux ensembles de données d'essai. La règle utilisée pour les données d'anticorps RB1 P-S780 est de définir la position de grille comme la demi-hauteur: la position de largeur maximale sur l'épaule gauche de la principale(À droite) de pointe lors de l'examen des données négatives de cellules de contrôle (non-ciblage siRNA). Données souligné rouge sont des cellules avec réduite et absent RB1 P-S780, qui sont identifiés avec cette porte. Une porte similaire placé sur l'épaule opposée de la distribution de valeur de rapport est utilisé pour le rapporteur GFP-CDK2. Les cellules de la sous-population à ratio élevé résultant, qui manquent ou la fonction réduction de l'activité de CDK2, sont indiquées en vert. Pour illustrer l'analyse multiplexée des deux dosages Figure 6C montre la mise en œuvre des deux valeurs de grille en utilisant le script 2_gate_classifier.pl perl pour convertir les données brutes (figure 6A) dans le fichier annoté ci-dessous. Ce nouveau fichier comprend des données d'origine le long d'une nouvelle colonne d'étiquettes de classe pour chaque cellule et les deux valeurs de grille utilisés pour les distinguer (dans ce cas des portes de 0,004 pour les données d'anticorps et 1,5 pour le rapporteur GFP-CDK2 ont été utilisés, respectivement) .
Ayant classé les cellules individuelles From chaque condition knockdown sur la base des deux essais, il est maintenant possible d'utiliser ces étiquettes de classe pour aider l'annotation de parcelles de données de titrage. La figure 7 montre dispersent parcelles de données cellulaires individuelles pour le P-S780 RB1 et GFP dosages de CDK2 de l'ensemble pour les trois conditions ARNi données d'exemple. Numéros annoter les quadrants sur les diagrammes de dispersion montrent les pourcentages relatifs de chaque sous-population fermée à l'ensemble de ce contexte d'effet de choc et sont générés dans la R en utilisant les étiquettes de classe décrites ci-dessus. Ces courbes montrent que, par rapport aux cellules transfectées avec des non-ciblage siRNA (figure 7B), les cellules transfectées avec siCDK6 révèlent une distribution de données net déplacé à la fois vers le bas sur l'axe des Y (indiquant l'absence d'RB1 phosphorylation sur la serine 780) et à droite, sur l'axe X (ce qui indique une faible activité de CDK2, Figure 7C). Ces deux changements sont attendus pour effet de choc de cet objectif. Contrairement à cela, les données de SirB1 cellules transfectées (figure 7A) montrer une perte de la coloration d'anticorps conforme à la perte de l'épitope, mais peu d'effet dans la distribution de données pour le journaliste de CDK2 par rapport aux témoins transfectées avec non-ciblage siRNA, suggérant l'absence de grand effet sur la GFP journaliste -CDK2 découle de RB1 knockdown.
Pour explorer davantage l'utilisation des données de la cellule individuelle, la classification des sous-population et le multiplexage dosage figure 8 montre le diagramme de dispersion pour les données siCDK6 de profils d'histogramme Figure 7C aux côtés appariés pour l'intensité de l'ADN intégré. Les paires d'histogrammes concernent moitiés opposées de l'ensemble de la population, répartis sur la base de l'intensité soit d'anticorps (de droite du diagramme de dispersion) ou GFP-CDK2 valeurs de rapport rapporteurs (au-dessus du nuage de points). La quantification de l'intensité de l'ADN nucléaire de ces populations montre deux pics caractéristiques de 2N et 4N teneur en ADN comme des pics de gauche et de droite, respectivement. Le intentions des grilles représentées sur les figures 6, 7 et 8 sont telles que les cellules identifiées comme faible pour P-S780 RB1 (marqué P-S780-) ou avec une forte valeur de rapport à partir du rapporteur GFP-CDK2 (marqués: G1) seront être en phase G1 du cycle cellulaire. En effet, les histogrammes de profils d'ADN pour les sous-populations identifiées avec l'une de ces essais contiennent principalement des cellules avec l'ADN contenu 2N. profils ADN de la population opposée fermée (étiquetés: P-S780 + ou non-G1) contient des cellules avec des distributions allant de 2n à 4N, en accord avec ces cellules adoptant une gamme de cycle cellulaire positionne phase post-G1.
Bien que l'accent ici a été la génération et l'analyse des données de cellules individuelles à partir d'images colorées par fluorescence, il est également utile d'être en mesure de prendre ces données et résumer chaque essai sur une base bien par puits pour surveiller la variabilité entre les répétitions et la performance de tous les puits pour un essai donné dans l'ensemble une plaque de données. A) les données de RB1 et B P-S780) les données rapporteurs GFP-CDK2. Les valeurs tracées dans A et B sont produits par deux scripts Perl supplémentaires fournies avec ce manuscrit; 'Antibody_fluorescence_summary.pl »et« G1assay_summary.pl', respectivement. Ces scripts utilisent les données brutes créées par Profiler cellulaire (Nuclei.csv) et les données du rapport ainsi que par i) cellules totales mesurées par puits, ii) le nombre de cellules dans la porte, iii) pourcentage de cellules dans la porte et iv) la moyenne arithmétique moyenne des données brutes, mesurées pour ce puits. Ce est inclus comme une option appropriée pour la recherche sur de vastes ensembles de données d'essai, avant de se concentrer sur les données individuelles de traitement utilisant l'évaluation des données multiplexé de cellules individuelles comme illustré dans 8. Les graphiques affichés ici complot »iii) pourcentage de cellules au sein de la porte» pour les deux dosages, qui conviennent les distributions de données non-normales observées pour les données de P-S780 RB1 et GFP-CDK2 dans les histogrammes de la figure 6B. Ces scripts calculent aussi 'iv) la moyenne arithmétique des données brutes, mesurées pour ce bien », qui conviendrait à l'analyse des données pour les réactions des populations homogènes et de la distribution normale des données avant et après la perturbation expérimentale.

Figure 1:. Vue d'ensemble des étapes du flux de travail pour l'analyse quantitative des données d'image de microscope marquées par fluorescence Le workflow est représenté ici comme quatre étapes (A) d'abord, il est nécessaire de préparer expérimentalement des cellules pour l'imagerie de fluorescence.. L'exemple décrit ici est celui d'uneécran dans lequel les cellules tumorales humaines adhérentes siARN traités sont cultivées pendant 48 heures, fixées et colorées sur une plaque de culture de tissu à 96 puits. Conditions ARNi différents sont présents en triple exemplaire dans des puits séparés dans la plaque. Les cellules sont colorées avec un colorant d'ADN, un anticorps spécifique pour RB1 phosphorylée sur la CDK4 et 6 sélectif site cible sérine-780 (P-S780 RB1) et ils ont aussi exprimer de manière stable un GFP-CDK2-reporter, les rapports de sortie du cycle cellulaire de G1. Collectivement, ces sondes fluorescentes constituent deux dosages évalués séparément dans le flux de travail. (B) images de microscope parallèle pour chaque sonde fluorescente (canal) sont générés et nommés de manière à inclure des détails par lequel le logiciel d'analyse d'image peut organiser les données. (C) Le fichiers d'images sont chargées dans le logiciel Cell Profiler, qui identifie algorithmique cellules individuelles et les paires associées de noyaux et le cytoplasme avant de céder mesures d'intensité pour les trois sondes fluorescentes détectented dans chaque. (D) Enfin, un script Perl est utilisé pour organiser les données quantitatives brutes produites. Cette étape se applique portes aux données d'intensité de fluorescence pour chaque cellule, binning efficacement les cellules en sous-populations, qui peuvent être tracées, suivis et contre-interrogé. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 2: Les données expérimentales à obtenir par analyse d'image, Les siRNA-traitées, les cellules marquées par fluorescence fixes de l'ensemble ont été imagées et des mesures d'intensité correspondantes prises par cellule données d'exemple.. Des données d'image représentatives sont indiquées pour chaque paramètre enregistré lors de l'analyse de l'image d'intensité (A) ADN nucléaire. L'intensité de la coloration du colorant de l'ADN nucléaire est utilisée pour donner une mesure de l'ADN par noyau (B) Intensité nucléaire de phospho-RB1:. Immuno-coloration spécifique pour P-S780 RB1 l'aide d'un (noir) anticorps primaire et fluorescence marqués anticorps secondaire (rouge) permettent une mesure d'intensité de RB1 phosphorylation au . S780 par noyau (C) GFP-CDK2 journaliste: Les cellules utilisées exprimer de façon stable une protéine rapporteur GFP-tagged que translocation entre le noyau et le cytoplasme dans un modèle d'ensemble avec le cycle cellulaire. Double mesure de l'intensité de la GFP nucléaire et cytoplasmique couplé pour chaque cellule permet le calcul d'un rapport par cellule qui peut être utilisée pour distinguer la phase G1 du reste du cycle cellulaire. Trois cibles d'ARNsi seront utilisés pour illustrer l'analyse; un non-ciblage siRNA contrôle négatif; CDK6 siRNA comme témoin positif dans RB1 perturber la phosphorylation et la progression du cycle cellulaire; RB1 siRNA pour établir la spécificité des anticorps.k "> Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 3: Organisation des fichiers image avant l'analyse d'image Les images prises à partir de la plaque de culture de tissus sont nommés systématiquement afin de permettre le logiciel d'analyse d'image de relier les données d'image replacer dans le contexte expérimental original.. Cette information est placée dans le nom du fichier pour chaque image. (A) Comme chaque puits sur la plaque d'expérimentation peut correspondre à différentes cibles ou traitements ARNi, les formes d'adresse et une partie du nom de fichier. (B) Le numéro de l'image fait partie du nom de fichier . que chaque puits est imagée pour recueillir multiples, non-chevauchement des cadres (C) sondes fluorescentes de chaque trame sont imagés séparément; par conséquent, les noms de fichiers doivent également refléter le canal chaque image se rapporte. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 4:. L'utilisation de cellules Profiler pour mesurer l'ADN nucléaire et la coloration des anticorps Avec les paramètres dans le fichier fourni de pipeline (3_channels_pipeline.cppipe), les valeurs d'intensité des mesures de logiciels d'analyse d'images de fluorescence cellulaire de Profiler pour l'ADN nucléaire et anticorps contraignant relatif à des cellules individuelles . (A) Les noyaux sont identifiés dans l'image de canal «bleu» de l'ADN coloré. (B) < / Strong> Les positions des noyaux d'ADN colorées sont détenues temporairement dans un «masque noyaux». Le masque est ensuite noyaux superposé à (C) les images bleues et rouges canal (ADN et de fluorescence de l'anticorps de données, respectivement) et les valeurs de fluorescence à partir de segments d'image qui se chevauchent avec le masque sont enregistrés contre chaque cellule identifiée. Identification réussie des noyaux voisins séparés, peut être évaluée visuellement à l'apparence du masque noyaux. Pour illustration, montré encerclé dans cette image de masque, sont des exemples où les paramètres choisis pour l'algorithme ont mis-identifiés noyaux voisins comme un seul noyau. Réglage des paramètres de l'algorithme pour minimiser ces événements est introduite dans la section Discussion. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.
/ftp_upload/51882/51882fig5highres.jpg "/>
Figure 5: Utilisation de la cellule pour mesurer les intensités de profils GFP nucléaires et cytoplasmiques La CDK2 rapporteur GFP-tagged translocation entre le noyau et le cytoplasme en fonction de la position du cycle cellulaire des cellules.. En même temps que cellule Profiler calcule les ADN et d'anticorps intensités nucléaires par cellule (figure 4), on calcule également le nucléaire à taux de cytoplasme des intensités GFP pour chaque cellule. (A) Les données ADN de colorant pour chaque image est utilisé pour générer un masque noyaux. (B) Cell Profiler utilise le masque de noyaux en liaison avec l'image de la GFP à partir du rapporteur GFP-CDK2 pour ensemencer la position de chaque cellule puis se dilate au périmètre de chaque cellule pour estimer l'ensemble de l'empreinte de chaque cellule. Cela devient une nouvelle, «masque cellulaire '. (C) Le masque de noyaux est soustrait du masque cellulaire pour obtenir une série de forme de beignet de cytoplasme contours, qui deviennentle «masque cytoplasme». (D) Le masque et le cytoplasme masque noyaux sont utilisés par Cell Profiler pour mesurer paires de valeurs de GFP nucléaires et cytoplasmiques. Ces valeurs appariées sont ensuite utilisés par Cell Profiler pour calculer les ratios, qui informent à la position de chaque cellule dans le cycle cellulaire. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 6:. Extraction des données - Traitement des données brutes de cellules individuelles en imposant des portes sur l'évolution des valeurs de dosage biologiques à partir des données de cellules individuelles pour la coloration d'anticorps et des dosages de rapporteur GFP CDK2 sont extraites en utilisant les données bloquées. Histogrammes des données brutes permettent d'identifier des valeurs de grille appropriés. Elles sont ensuite imposées avec un script Perl. (A) Leproduit de l'analyse des fichiers d'image avec les paramètres fournis pour cellulaire Profiler final sont séparées par des virgules valeur (.csv). Ces fichiers C ontain données individuelles de cellules relatives à chacun des différents segments de sous-cellulaire. Le fichier 'Nuclei.csv' contient toutes les mesures sélectionnées relatives à l'utilisation du masque noyaux. Ces mesures comprennent l'intensité des anticorps nucléaire, l'intensité de l'ADN nucléaire, et le rapport de la GFP (noyau / cytoplasme). (B) histogrammes d'intensité d'anticorps nucléaire (à gauche) et les rapports de rapporteurs GFP-CDK2 (à droite) tracée à partir des données de cellules individuelles pour chaque état knockdown siARN . Les barres sur les histogrammes affichés montrent les positions de départ désirées pour ces essais. Des données de couleur sur les histogrammes indiquent les sous-populations gated. (C) Les portes pour les deux essais illustrés dans B sont appliquées aux données brutes en utilisant le script Perl «2_gate_classifier.pl '. Le scriptcrée une copie modifiée du fichier de sortie d'origine Cell Profiler (Nuclei.csv) pour aider tracer par la suite. Les deux valeurs de grille sont enregistrés dans le nouveau fichier (surligné en couleur ici) et une nouvelle colonne «Label» est ajouté. Le bac étiquettes chaque cellule dans l'une des quatre sous-groupes possibles en fonction des deux valeurs de dosage bloquées pour chaque cellule. Ces étiquettes sont utilisées dans les parcelles suivantes qui disposent calculs des contributions de chaque sous-population ainsi que le recoupement des paramètres supplémentaires générés dans Cell Profiler. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 7:. Les diagrammes de dispersion pour chaque condition siRNA représentant les données brutes de cellules individuelles et les positions de grille Scatter parcelles de individouble données cellulaires à partir de toutes les images pour les conditions siARN indiqué: (A) siRB1; (B) de commande de ciblage Sinon négative; (C) siCDK6. Complotèrent contre les axes Y sont des valeurs de fluorescence nucléaire de l'anti-P-S780 RB1 coloration. Comploté contre les X-axes sont les valeurs de rapport correspondantes calculées à partir du rapporteur GFP-CDK2. Les barres rouges et vertes indiquent les positions des portes pour la porte de la RB1 P-S780 et les portes rapporteurs GFP-CDK2, respectivement. Les deux grilles des cellules se divisent en quatre sous-populations et les nombres résultants sur les quadrants représentent le nombre de cellules de pour cent de chacun de ceux-ci. Annotations autour des axes pour A indiquent les quatre étiquettes-éléments possibles appliquées à chaque cellule par le 2_gate_classifier.pl script Perl. Ces étiquettes sont représentés par rapport à leur grille respective de dosage et sont utilisés dans le R-écriture (analysis.r) pour générer les tracés des figures 6, 7 et Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 8: cellulaires sous-populations définies par les deux dosages de transit G1 montrent profils 2N et 4N ADN en accord avec dosage résultats Le nuage de points de données pour les cellules siCDK6 est répété à partir de la figure 7C.. Entourant le nuage de points sont des histogrammes pour intégrée intensité de l'ADN nucléaire relatif à des sous-ensembles de la population. Ceux au-dessus du nuage de points se rapportent à l'essai rapporteur GFP-CDK2. Ceux à la droite de la nuage concernent nucléaire phospho-RB1 anticorps mesures seuls. Les lignes de grille de couleur sont étendues à montrer leur relation avec les histogrammes. Étiquettes porte par laquelle les données de la cellule ont été sélectionnés pour ces parcelles supplémentaires sont également présentés. Les cellules avec perte de RB1 phosphorylée sur des sérine 780 (P-S780-) ou ceux avec une GFP-CDK2 journaliste nucléaires de haute ratio cytoplasmique (indiquant une faible activité de CDK2) montrent des profils d'ADN principalement 2N-like, alors que leurs homologues opposées pour chaque dosage respective montrer une distribution de 2N et 4N, caractéristique d'un, post-G1 population de phase mixte de cellules. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 9: Résumé des parcelles de valeurs de dosage pour chaque gated siRNA état parcelles de données Résumé de gated (a) P-S780 données RB1 et données (B) GFP-CDK2 de puits en triple pour chaque condition knockdown. SiRNA. Les valeurs ont été calculées à partir de la sortie brute de Profiler cellulaire (Nuclei.csv) en utilisant leScripts Perl, 'antibody_fluorescence_summary.pl' (A) ou «G1assay_summary.pl '(B). Les valeurs tracées sont des moyens de le pourcentage de cellules au sein de la grille appliquée à chaque dosage. Les barres indiquent les erreurs-types calculées à partir des puits en triple. , Les valeurs homoscédastiques test t non apparié P pour chaque condition knockdown par rapport à la non-ciblage siRNA sont indiquées ci-dessus les données tracées où P <0,001 (**) et P <0,05 (*). Se il vous plaît cliquer ici pour voir une version plus grande de cette figure.
Figure S1. Mise en place logiciel Cell Profiler pour l'analyse d'image. (A) Screenshot of Cell Profiler avant les paramètres d'analyse d'images sont saisies. (B) Screenshot of Cell Profiler après les détails de l'algorithme contenues dans '3_channels_pipeline.cppipe' ont été chargés. La haute onglet allumé dans le coin supérieur gauche indique que cet écran affiche les paramètres pour l'étape "loadImages» de l'analyse. En cliquant sur les autres parties de la liste ci-dessous cela révéler les détails pour les étapes ultérieures de l'analyse. (C) Screenshot of Cell Profiler avec des détails pour dossier d'entrée et de sortie de dossier est entré. (D) Screenshot of Cell Profiler après la «Analyser le bouton de images a été cliqué pour commencer l'analyse. Superposées sont trois nouvelles fenêtres illustrant les masques algorithmique-produits générés par le logiciel à partir des images en cours d'analyse. Ces fenêtres sont accessibles en cliquant sur les icônes «oculaires» à la position ouverte à côté des mesures pertinentes dans l'analyse, dans le coin supérieur gauche de la fenêtre principale Cell Profiler. Ces vues permettent à l'utilisateur de vérifier si les paramètres générant les masques confettis de couleur en accord avec les données en niveaux de gris, d'accompagnement d'origine.
ent ">
Figure S2. L'utilisation de Perl et RStudio aux données de cellule individuelle de grille et de tracer les sous-populations de cellules qui en résultent. (A) Le panneau de droite montre le dossier choisi de recevoir les dossiers de la sortie (icônes vertes) de l'analyse cellulaire de Profiler. Les scripts Perl fournis avec le manuscrit (icônes bleues) sont copiés dans ce dossier. surbrillance est le script Perl
«2_gate_classifier.pl», qui a été double-cliqué avec la souris pour produire la boîte de dialogue dans le panneau de gauche. indiquées sont les invites et réponses correspondant tapé nécessaires à la porte les données de cellules individuelles à partir du fichier
'Nuclei.csv'. (B) Capture d'écran de RStudio immédiatement après le chargement du script
'analysis.R'. surbrillance sont les commandes de télécharger les données de
A gated dans le logiciel avant de traçage (détails de notes dans les lignes 5 et 6 devra être ajustée en fonction de l'endroit où les données fermée se trouve sur le calculr utilisé pour l'analyse).
(C) Capture d'écran de RStudio fois que les données ont été téléchargées.
(D) Capture d'écran de RStudio démontrant souligné le bloc de code requis pour produire le tracé affiché dans la fenêtre en bas à droite. Codes pour chaque parcelle sont séparés par des lignes vides et regroupés par type de tracé.
| cible siRNA | Plaque de puits adresses |
| Non-ciblage (NT) | E5, F5, G5 |
| Rétinoblastome (RB) | E7, F7, G7 |
| Kinase cycline-dépendante 6 (CDK6) | B2, C2, D2 |
Tableau 1: Well adresses et conditions siARN correspondants utilisés dans l'ensemble de données d'exemple.