$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Quantifizierung des Kernvolumens in Tumorsphäroiden in 3D
Die automatisierte 3D-Mikroskopie (z-Stacks) ermöglicht es Wissenschaftlern, komplexe mehrzellige Systeme wie Organoide sichtbar zu machen. Um morphologische und Fluoreszenzsignaleigenschaften auf Einzelzellebene zu quantifizieren, ist häufig eine Zellsegmentierung in 3D erforderlich. Das folgende Beispiel zeigt, wie Cell-ACDC die Zellkerne in Organoiden, die Tausende von Zellen aus dem Kernfärbekanal enthalten, segmentieren und ihr Volumen quantifizieren kann (Daten aus Ref.9). Die Segmentierung wurde mit einem benutzerdefinierten Cellpose-Modell durchgeführt, das trainiert und in Ref.9 veröffentlicht wurde. Das trainierte Cellpose-Modell kann direkt in Cell-ACDC verwendet werden, indem der Pfad der Gewichtungsdatei in der GUI angegeben wird, wenn die Modellparameter für Cellpose v2 ausgewählt werden. Die Segmentierung wurde mit dem zweiten Modul von Cell-ACDC durchgeführt, um alle Bilder im Stapelverfahren zu verarbeiten (Abbildung 1A, Modul "Segmentieren und verfolgen"). Als nächstes wurde das Ergebnis in der GUI des dritten Moduls visualisiert (Abbildung 1A, Modul "Visualisieren und korrigieren"). Dieses Modul wurde optimiert, um Tausende von Einzelobjekten mit mehreren Anmerkungsoptionen (z. B. Konturen, überlagerte Segmentierungsmasken oder Text-IDs) zu visualisieren. Nachdem die Messungen aus den 3D-Masken berechnet wurden, wurden alle verfügbaren Messungen direkt in der GUI dokumentiert. Sie sind zugänglich, indem Sie in der oberen Menüleiste (Abbildung 6, "Menüleiste") das Menü Messungen und dann Messungen einstellen auswählen. Im Popup-Dialogfeld können Benutzer über die Info-Schaltflächen messspezifische Informationen abrufen und auswählen, welche Messungen gespeichert werden sollen. Für die repräsentativen Ergebnisse in Abbildung 10 wurde "cell_vol_fl_3D" verwendet, d. h. das Volumen jedes Objekts, das durch Multiplikation der Gesamtvoxel im Objekt mit der Pixelgröße zum Quadrat und der Voxeltiefe berechnet wird. Diese Eigenschaften werden entweder automatisch aus der Mikroskopie-Rohdatei extrahiert oder vom Benutzer bereitgestellt. Die Berechnung der Messungen über diese GUI erfordert das Laden der Rohbilder. Um den Prozess zu rationalisieren und die Stapelverarbeitung zu ermöglichen, können die Messungen auch über das Menü "Dienstprogramme " (Abbildung 1B, "Dienstprogramme") berechnet werden, indem Sie zum Untermenü "Messungen " und dann zu "Messungen berechnen" für ein oder mehrere Experimente gehen. Schließlich wurde die nukleare Volumenverteilung als repräsentatives Ergebnis aufgetragen (Abbildung 10). Diese Analyse zeigt einen signifikanten Anteil kleiner Zellkerne, wahrscheinlich aufgrund von Segmentierungsartefakten. Sie können leicht entfernt werden, indem kleine Objekte aus der Segmentierungsmaske herausgefiltert werden. Gleichzeitig könnten die sehr großen Kerne auf verschmolzene Kerne während der Segmentierung zurückzuführen sein. Es wird immer empfohlen, die Volumenverteilung von Objekten (z. B. einzelnen Zellen) zu zeichnen, um Artefakte zu identifizieren und zusätzliche biologische Informationen über die Zellgröße zu extrahieren.
Quantifizierung von Zeitraffer-Mikroskopie-Daten
Mit der Zeitraffermikroskopie kann die zelluläre Dynamik direkt auf Einzelzellebene beobachtet werden. Neben der Zellsegmentierung erfordert die Extraktion der zeitlichen Dynamik zusätzliche Analysen, einschließlich Zellverfolgung und Annotation von Zellstammbäumen. Aufgrund der gegenseitigen Abhängigkeit dieser Analysephasen können Fehler, die früh in der Pipeline auftreten, auf spätere Schritte übertragen werden. Infolgedessen ist eine kontinuierliche Visualisierung und Korrektur von Segmentierungs-, Tracking- und Annotationsfehlern erforderlich. Im Folgenden wird gezeigt, dass Cell-ACDC für diese Aufgaben geeignet ist. Es wurden zwei Datensätze von zwei verschiedenen Modellorganismen ausgewählt: 1) knospende Hefe (Stamm DCY001-1 aus Ref.35, wo die beiden Histon-H2B-Proteine Htb1 und Htb2 mit mCitrin markiert sind) und 2) embryonale Stammzellen der Maus (mESCs, Daten aus Ref.22).
Diese beiden Datensätze zeigen zwei Annotationsmodi, die in Cell-ACDC verfügbar sind: asymmetrische und symmetrische (d. h. symmetrische Zytokinese oder "normale") Zellteilung. Da sich die Art der Teilung unterscheidet, benötigen die beiden Organismen einen unterschiedlichen Tracking- und Annotationsrahmen.
Bei der asymmetrischen Teilung bildet die Mutterzelle eine Knospe, die wächst und sich schließlich trennt, um zur Tochterzelle zu werden. Nach der Teilung behält die Mutterzelle ihre ursprüngliche Zell-ID bei, und ihre Generationsnummer erhöht sich um eins, während der Tochterzelle eine neue Zell-ID zugewiesen wird und ihre Generationsnummer auf eins gesetzt wird. Die Knospenphase entspricht den S/G2/M-Phasen des Zellzyklus und wird als solche in Cell-ACDC annotiert. Diese Annotationsoptionen helfen dabei, typische biologische Fragen im Zusammenhang mit dem Zellzyklus in knospenden Hefen zu beantworten.
Bei der "symmetrischen" Teilung (z.B. Säugetierzellen) teilt sich die Mutterzelle in zwei Tochterzellen. Die Mutterzelle, also ihre ID, verschwindet bei der Teilung, und die beiden Tochterzellen erhalten neue IDs. Zusätzlich wird die Generationszahl der Tochterzellen um eins relativ zur Mutterzelle erhöht. Cell-ACDC verfolgt auch die übergeordnete ID, die Root-ID (die ursprüngliche Vorfahrzelle am Anfang einer Abstammungslinie) und die Schwester-ID.
Für beide Annotationsmodi wurde ein innovatives Framework zur Korrektur von Annotationsfehlern entwickelt, bei dem die Korrektur automatisch auf alle vergangenen und zukünftigen relevanten Zeitpunkte übertragen wird. Die Visualisierung, Annotation und Korrektur erfolgte in der GUI des dritten Moduls (Abbildung 1A, Modul "Visualisieren und korrigieren" und Abbildung 6). Um die Zellen in Datensatz 1 zu segmentieren und zu verfolgen, wurde das Modell YeaZ_v226 auf den Phasenkontrastkanal angewendet, während für den Kernkanal (Histon) das StarDist25-Modell verwendet wurde. Anschließend wies Cell-ACDC mit dem Dienstprogramm Tracking and lineage > Tracking and/or count sub-cellular objects (Abbildung 1B, Menüleiste "Dienstprogramme ") jedem Zellkern die Zell-ID der entsprechenden Zelle zu, um die Konsistenz zwischen den aus Zell- und Zellkernmasken generierten Tabellen sicherzustellen.
Für Datensatz 2 wurde das Segmentierungsmodell DeepSea22 verwendet. Alle drei Modelle sind bereits in Cell-ACDC verfügbar, was den Vorteil der Integration mehrerer Segmentierungsmodelle in die Software zeigt.
Nach der Korrektur von Segmentierungs- und Tracking-Fehlern wurden Zellstammbäume annotiert, numerische Merkmale berechnet und eine nachgelagerte Analyse durchgeführt. Für Datensatz 1 wurden die Spalte TaYFP_amount_autoBkgr und die Anzahl der segmentierten Kerne (Abbildung 11A) gegen die Zeit aufgetragen. "TaYFP" ist der Name des nuklearen Kanals. "amount_autoBkgr" ist ein Proxy für die gesamte zelluläre Proteinmenge, die aus Epifluoreszenzbildern extrahiert wird36. Sie wird berechnet als die Differenz zwischen der mittleren Fluoreszenzintensität in jeder Zellmaske und dem Hintergrundmedian, multipliziert mit der Zellfläche (in Pixel). Hier wird der Hintergrundmedian aus allen Pixeln berechnet, die nicht als Zellen segmentiert sind. Wie erwartet, beginnt die H2B-Menge mit dem Austreiben der Knospen zu steigen (Abbildung 11A-ii) und erreicht vor der Kernteilung einen konstanten Wert (Abbildung 11A-iii). Dies ist eine wichtige Qualitätskontrolle bei der Homöostase von Histonproteinen, da davon ausgegangen wird, dass die Menge der Histonproteine vom Zellzyklus abhängt. Darüber hinaus bestätigt die Darstellung der Anzahl der Kerne im Laufe der Zeit, dass die Mengen an Histonproteinen das Maximum ungefähr um die Kernteilung herum erreichen.
Für Datensatz 2 wurde die Zellfläche im Zeitverlauf für eine ausgewählte Zelle, die sich in der Zellteilung befindet, dargestellt. Erwartungsgemäß nimmt die Zellfläche bis zum Erreichen eines Maximalwertes zu (Abbildung 11B-i). Dann nimmt sie bis zur Zellteilung (Abbildung 11B-ii) ab, wenn sich die Zelle zusammenzieht. Schließlich beginnt der Zyklus für die beiden Tochterzellen neu. Dies ist eine weitere empfohlene Analyse, da die Überprüfung der Zellgrößenänderungen im Laufe des Zellzyklus unerlässlich ist, um zu bestätigen, dass die Zellen wie erwartet wachsen und sich teilen (oder nicht im Falle bestimmter Mutanten).

Abbildung 1: Cell-ACDC-Module. (A) Überblick über die 4 Hauptmodule, die von der Haupt-Cell-ACDC-Trägerrakete gestartet werden können. Nach dem Segmentieren, Verfolgen und Annotieren von Mikroskopiedaten können numerische Merkmale entweder aus dem dritten Modul ("Visualisieren und korrigieren") oder aus (B) dem Menü "Dienstprogramme " in der oberen Menüleiste berechnet werden. Die "Dienstprogramme" sind die Routinen, die automatisch auf mehreren Datensätzen ohne Benutzereingabe ausgeführt werden können. Neben der Berechnung der Messungen gehören zu den weiteren Dienstprogrammen die Verkettung mehrerer Ausgabetabellen in einer einzigen Tabelle, die Verfolgung subzellulärer Objekte und die Bildvorverarbeitung. Cell-ACDC unterstützt 2D-, 3D- (Z-Stack oder Zeitraffer) und 4D-Daten (Z-Stacks im Zeitverlauf) mit beliebig vielen zusätzlichen Kanälen. Die Ausgabetabelle mit den numerischen Merkmalen kann dann für nachgelagerte Analysen und biologische Entdeckungen verwendet werden (Abbildung 10 und Abbildung 11). Zu diesem Zweck stehen Jupyter-Notebooks auf der GitHub-Seite Cell-ACDC zur Verfügung, die Beispiele für Diagramme enthalten, die aus der Ausgabetabelle abgerufen werden können. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 2: Flussdiagramm der Cell-ACDC-Entscheidung. Ein Flussdiagramm, in dem das Modul beschrieben wird, das je nach Datensatztyp und Analyseanforderungen verwendet werden soll. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 3: Laden Sie Beispieldaten herunter. (A) Offener Willkommensleitfaden. (B) Laden Sie Beispieldaten herunter, die zum Replizieren des Protokolls erforderlich sind. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 4: Laden von Daten für die Datenvorbereitung. (A) Laden Sie Daten in die GUI der Datenvorbereitung. (B) Wählen Sie den zu ladenden Kanal aus. (C) Bearbeiten und bestätigen Sie die Bildmetadaten. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 5: Ausführen des Datenvorbereitungsprozesses. (A) Starten Sie den Vorgang. (B) Positions-ROIs (für das Zuschneiden) und Hintergrund-ROIs. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 6: GUI des dritten Moduls zur Visualisierung und Korrektur der Ergebnisse. Screenshot der GUI des dritten Moduls ("Visualisieren und korrigieren" in Abbildung 1A) mit hervorgehobenen Elementen. Beachten Sie, dass die meisten Schaltflächen in den Symbolleisten über einen Tooltip verfügen (auf den Sie zugreifen können, indem Sie den Mauszeiger über die Schaltfläche bewegen), der erklärt, wie diese spezielle Funktion verwendet wird. Mit dem Moduswähler kann zwischen 5 Modi umgeschaltet werden: "Viewer", "Segmentierung und Tracking", "Zellzyklusanalyse" (für asymmetrisch teilende Zellen), "Normalteilung: Abstammungsbaum" (für symmetrisch teilende Zellen, z. B. Säugetierzellen) und "Benutzerdefinierte Anmerkungen". Beachten Sie, dass in den anderen GUIs häufig eine Symbolleiste oder Menüleiste vorhanden ist (z. B. im Modul "Datenvorverarbeitung", Abbildung 1). Die Bearbeiten-Symbolleiste enthält alle Funktionen, mit denen Segmentierungs- und Tracking-Fehler bearbeitet und korrigiert werden können (z. B. Pinsel, Radiergummi, Editier-ID usw.). Das geladene Bild wird in einer zweiteiligen Ansicht angezeigt, was hilfreich ist, wenn verschiedene Anmerkungsoptionen benötigt werden (z. B. Zellzyklusinformationen auf dem linken Bild und IDs auf dem rechten Bild). Das rechte Bild kann auch ausgeschaltet werden (klicken Sie mit der rechten Maustaste auf das Bild und deaktivieren Sie die Option Gespiegeltes Bild anzeigen). Jedes Bildfenster verfügt über einen LUT-Schieberegler an der Seite, mit dem Sie die Intensitätsstufen schnell anpassen können. Mit einem Rechtsklick auf den LUT-Regler kann der Benutzer verschiedene Farbkarten für die Intensitätsbilder auswählen. Zusätzlich gibt es auf der rechten Seite einen LUT-Selektor für die Farbe der Segmentierungsbeschriftungen, die über die Intensitätsbilder gelegt werden sollen ( Anmerkungsoption namens Segm.-Masken). Auf der linken Seite der Anmerkungsoptionen für das linke Bild gibt es zusätzliche Umschalter, mit denen Sie einige Einstellungen steuern können, wie z. B. automatisches Speichern, Schriftgröße usw. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 7: Visualisieren Sie die Vorverarbeitung in der Haupt-GUI. (A) Öffnen Sie das Dialogfeld für die Vorverarbeitung. (B) Vorverarbeitungsdialog mit Parametern, die im initialisierten Protokoll verwendet werden. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 8: Visualisieren Sie die Segmentierungsausgabe in der Haupt-GUI. (A) Wählen Sie ein Segmentierungsmodell aus. (B) Richten Sie Segmentierungsparameter ein. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 9: Von Cell-ACDC benötigte Ordnerstruktur. Um mit Cell-ACDC arbeiten zu können, müssen die Daten in einer bestimmten Ordnerstruktur angeordnet werden. Obwohl Cell-ACDC ein Modul zur automatischen Generierung dieser Struktur bereitstellt, ist es wichtig zu verstehen, wie diese Struktur aussehen soll. Zunächst müssen sich alle Dateien in einem Ordner mit dem Namen "Bilder" befinden. Als nächstes müssen sie alle mit dem gleichen Namen beginnen, dem sogenannten "basename_". Der Mindestsatz der erforderlichen Dateien besteht aus einer Einkanal-TIFF-Datei (2D, 3D-Z-Stapel oder 3D+Zeit) und einer CSV-Datei, die mit "_metadata.csv" endet. Bei dieser Datei muss es sich um eine Tabelle mit zwei Spalten handeln, wobei die erste Spalte "Beschreibung " und die zweite Spalte "Werte" lauten muss, und sie sollte mindestens Einträge für SizeT und SizeZ für die Anzahl der Frames bzw. Z-Slices enthalten. Wenn eine Bilddatei keine Z-Slices enthält, muss SizeZ auf 1 festgelegt werden. Das Gleiche gilt für Bilder ohne Zeitraffer, bei denen SizeT 1 sein muss. Da es sich um eine CSV-Datei handelt, besteht ein Eintrag aus einer einzelnen Zeile von Description,value, die durch ein Komma getrennt ist, z. B. SizeZ,1. Für mehrere Kanäle muss eine TIFF-Datei pro Kanal generiert werden. Der Ordner Bilder muss dann in einem Ordner mit dem Namen "Position_1" abgelegt werden. Es sind mehrere Positionen zulässig, die mit einer fortlaufenden Nummer benannt werden müssen. Beim Laden von Daten in eines der Cell-ACDC-Module kann der Benutzer entweder einen bestimmten Positionsordner oder den gesamten Experimentordner auswählen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 10: 3D-Quantifizierung von Tumor-Organoiden. Screenshot eines repräsentativen Tumor-Organoids (Daten von9), das in die GUI des dritten Moduls von Cell-ACDC geladen wurde (links), Beispiel-Z-Schichten mit roten Konturen, die die Segmentierungsmasken hervorheben (Mitte), und Histogramm der Zellvolumenverteilung, berechnet aus den 3D-Segmentierungsmasken (rechts). Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 11: Quantifizierung von Zeitraffer-Mikroskopie-Daten. (A) Screenshot von Zeitraffer-Mikroskopiedaten von knospenden Hefezellen, die in die GUI des dritten Moduls von Cell-ACDC geladen wurden (links), und Quantifizierung der Histon-H2B-Proteinmenge über die Zeit in einem repräsentativen Zellzyklus (rechts). Die vergrößerten Bilder zeigen die Beispielzelle und ihre Knospe (weiße Pfeile) zu Beginn des Zellzyklus (i), beim Knospenaufgang (ii) und bei der Kernteilung (iii). Die Daten stammen von Chatzitheodoridou et al.35 (B) Screenshot von Zeitraffermikroskopiedaten von embryonalen Stammzellen der Maus, die in das dritte Modul geladen wurden, GUI von Cell-ACDC (links) und Zellfläche (μm2), aufgetragen als Funktion der Zeit einer repräsentativen Zelle, die sich in der Zellteilung befindet. Die vergrößerten Bilder zeigen die Beispielzelle und ihre Tochterzellen bei maximaler Zellfläche vor der Teilung (i), der Teilung in zwei Tochterzellen (ii) und dem letzten analysierten Frame nach der Teilung (iii). Die Daten stammen von Zargari et al.22,33. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.