Method Article

AMEBaS: Automatische Mittellinienextraktion und Hintergrundsubtraktion von ratiometrischen Fluoreszenz-Zeitraffern polarisierter Einzelzellen

DOI:

10.3791/64857

June 23rd, 2023

In This Article

Summary

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Derzeitige Methoden zur Analyse der intrazellulären Dynamik polarisierter Einzelzellen sind oft manuell und nicht standardisiert. Dieses Manuskript stellt eine neuartige Bildanalyse-Pipeline zur Automatisierung der Mittellinienextraktion einzelner polarisierter Zellen und zur Quantifizierung des raumzeitlichen Verhaltens aus Zeitraffern in einer benutzerfreundlichen Online-Oberfläche vor.

Abstract

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Zellpolarität ist ein makroskopisches Phänomen, das durch eine Ansammlung von räumlich konzentrierten Molekülen und Strukturen entsteht, die in der Entstehung spezialisierter Domänen auf subzellulärer Ebene gipfeln. Es ist mit der Entwicklung asymmetrischer morphologischer Strukturen verbunden, die wichtigen biologischen Funktionen wie Zellteilung, Wachstum und Migration zugrunde liegen. Darüber hinaus wurde die Störung der Zellpolarität mit gewebebezogenen Erkrankungen wie Krebs und Magendysplasie in Verbindung gebracht.

Derzeitige Methoden zur Bewertung der raumzeitlichen Dynamik von fluoreszierenden Reportern in einzelnen polarisierten Zellen beinhalten oft manuelle Schritte, um eine Mittellinie entlang der Hauptachse der Zellen zu verfolgen, was zeitaufwändig und anfällig für starke Verzerrungen ist. Obwohl die ratiometrische Analyse die ungleichmäßige Verteilung von Reportermolekülen mit Hilfe von zwei Fluoreszenzkanälen korrigieren kann, sind die Techniken der Hintergrundsubtraktion häufig willkürlich und nicht statistisch abgesichert.

Dieses Manuskript stellt eine neuartige Rechenpipeline vor, um das raumzeitliche Verhalten einzelner Zellen unter Verwendung eines Modells der Zellpolarität zu automatisieren und zu quantifizieren: Wachstum der Pollenröhren/Wurzelhaare und zytosolische Ionendynamik. Es wurde ein dreistufiger Algorithmus entwickelt, um ratiometrische Bilder zu verarbeiten und eine quantitative Darstellung der intrazellulären Dynamik und des Wachstums zu extrahieren. Der erste Schritt segmentiert die Zelle vom Hintergrund und erzeugt eine binäre Maske durch eine Schwellwerttechnik im Pixelintensitätsraum. Im zweiten Schritt wird ein Weg durch die Mittellinie der Zelle durch eine Skelettierungsoperation nachgezeichnet. Der dritte Schritt schließlich liefert die verarbeiteten Daten als ratiometrischen Zeitraffer und ergibt einen ratiometrischen Kymographen (d.h. ein räumliches 1D-Profil über die Zeit). Daten aus ratiometrischen Bildern, die mit genetisch kodierten fluoreszierenden Reportern aus wachsenden Pollenschläuchen aufgenommen wurden, wurden verwendet, um die Methode zu vergleichen. Diese Pipeline ermöglicht eine schnellere, weniger verzerrte und genauere Darstellung der raumzeitlichen Dynamik entlang der Mittellinie polarisierter Zellen und erweitert so das quantitative Instrumentarium, das zur Untersuchung der Zellpolarität zur Verfügung steht. Der AMEBaS-Python-Quellcode ist verfügbar unter: https://github.com/badain/amebas.git

Introduction

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Die Zellpolarität ist ein grundlegender biologischer Prozess, bei dem die konzertierte Wirkung einer Ansammlung von räumlich konzentrierten Molekülen und Strukturen in der Etablierung spezialisierter morphologischer subzellulärer Domänen gipfelt1. Zellteilung, Wachstum und Migration hängen von solchen Polaritätsstellen ab, während ihr Verlust mit Krebs bei epithelialen Gewebeerkrankungen in Verbindung gebracht wurde2.

Apikal wachsende Zellen sind ein dramatisches Beispiel für Polarität, bei der sich die Polaritätsstelle an der Spitze typischerweise auf extrazelluläre Signale umorientiert3. Dazu gehören die Entwicklung von Neuriten, Pilzhyphen, Wurzelhaaren und Pollenschläuchen, bei denen mehrere zelluläre Prozesse ausgeprägte Unterschiede von der Zellspitze zum Schaft hin zeigen. Insbesondere in Pollenschläuchen sind die Aktinpolymerisation, der Vesikeltransport und die Ionenkonzentrationen deutlich polarisiert und zeigen spitzenfokussierte Gradienten4. Pollenschläuche sind die männlichen Gametophyten von Blütenpflanzen und sind dafür verantwortlich, die Samenzellen in die Eizelle zu bringen, indem sie ausschließlich an der Spitze der Zelle mit einer der schnellsten Wachstumsraten wachsen, die für eine einzelne Zelle bekannt sind. Die spitzenfokussierten Gradienten von Ionen wie Calcium 5 (Ca2+) und Protonen 6 (H+) spielen eine wichtige Rolle bei der Aufrechterhaltung des Pollenschlauchwachstums, das für die Erfüllung seiner biologischen Hauptfunktion, die in einer doppelten Befruchtung gipfelt, unerlässlich ist 5,6. Daher sind quantitative Methoden zur Analyse der raumzeitlichen Dynamik entlang der Mittellinie apikal wachsender Zellen unerlässlich, um die zellulären und molekularen Mechanismen zu untersuchen, die dem polarisierten Wachstum zugrunde liegen 7,8,9. Forscher verwenden häufig Kymographen, d. h. eine Matrix, die die Pixelintensitäten der Mittellinie der Zelle (z. B. Spalten) im Laufe der Zeit (z. B. Zeilen) darstellt und es ermöglicht, das Zellwachstum und die Zellwanderung in der Diagonalen zu visualisieren (Abbildung 1). Trotz ihrer Nützlichkeit werden Kymographen häufig durch manuelles Nachzeichnen der Mittellinie extrahiert, was anfällig für Verzerrungen und menschliche Fehler ist und gleichzeitig ziemlich mühsam ist. Dies erfordert ein automatisiertes Verfahren der Mittellinienextraktion, das das erste Merkmal der hier vorgestellten Pipeline mit dem Namen AMEBaS ist: Eineutomatische M-Idline-E-Xtraktion und eine Ba ckground-S-Traktion von ratiometrischen Fluoreszenzzeitraffern polarisierter Einzelzellen.

In Bezug auf experimentelle Verfahren kann die quantitative Abbildung von Ionen/Molekülen/Spezies von Interesse in einzelnen Zellen mit genetisch kodierten Fluoreszenzsondenerreicht werden 10. Unter den ständig wachsenden Möglichkeiten sind ratiometrische Sonden eine der genauesten, da sie unterschiedliche Fluoreszenzwellenlängen emittieren, wenn sie an die interessierenden Moleküle gebunden oder ungebunden sind11. Dies ermöglicht die Korrektur der räumlichen Heterogenität in der intrazellulären Konzentration der Sonde, indem das Verhältnis zweier Kanäle mit ihrem kanalspezifischen Hintergrund subtrahiert wird. Die Schätzung des Hintergrundschwellenwerts für jeden Kanal und jeden Zeitpunkt kann jedoch eine komplexe Aufgabe sein, da er oft räumlich aufgrund von Effekten wie Schattierungen, bei denen die Ecken des Bildes Helligkeitsschwankungen relativ zur Mitte aufweisen, und im Laufe der Zeit aufgrund des Verblassens des Fluorophors (Photobleaching)12 variiert. Obwohl es mehrere mögliche Methoden gibt, schlägt dieses Manuskript vor, die Hintergrundintensität automatisch unter Verwendung der mit dem Isodata-Algorithmus13 erhaltenen Segmentierungsschwelle zu bestimmen, die dann durch polynomiale Regression als Standard über Frames geglättet wird. Räumliche Komponenten, die sich aus der Fluoreszenzheterogenität ergeben, die nicht mit der entfernten Zielzelle inVerbindung steht, wurden bei dieser Methode jedoch ignoriert. Die automatische Schwellenwertbestimmung kann mit verschiedenen Methoden durchgeführt werden, aber der Isodata-Algorithmus lieferte empirisch die besten Ergebnisse. Daher sind die automatische Hintergrundwertsubtraktion und die ratiometrische Berechnung das zweite Hauptmerkmal von AMEBaS (Abbildung 1), das zusammengenommen einen Stapel von Zweikanal-Fluoreszenzmikroskopiebildern als Eingabe empfängt, die Mittellinie der Zelle und den kanalspezifischen Hintergrund schätzt und Kymographen beider Kanäle und ihres Verhältnisses (Hauptausgang #1) nach Hintergrundsubtraktion, Glättung und Ausreißerentfernung ausgibt. zusammen mit einem Stapel ratiometrischer Bilder (Hauptausgabe #2).

AMEBaS wurde mit Fluoreszenzzeitraffern von wachsenden Arabidopsis-Pollenröhrchen getestet, die unter einem Mikroskop aufgenommen wurden, entweder mit ratiometrischen Ca2+ (CaMeleon)8 oder pH (pHluorin)6 Sensoren, die unter dem pollenspezifischen LAT52-Promotor exprimiert wurden. Die Bilder von jedem Kanal wurden alle 4 s aufgenommen, gekoppelt mit einem inversen Mikroskop, einer frontbelichteten Kamera (2560 Pixel × 2160 Pixel, Pixelgröße 6,45 μm), einem Fluoreszenzstrahler und einem Wasserimmersionsobjektiv 63x, 1,2NA. Die für CaMeleon verwendeten Filtereinstellungen waren: Anregung 426-450 nm (CFP) und 505-515 nm (YFP), Emission 458-487 nm (CFP) und 520-550 nm (YFP), während für pHluorin Anregung 318-390 nm (DAPI) und 428-475 nm (FITC), Emission 435-448 nm (DAPI) und 523-536 nm (FITC). Ein vollständiger Datensatz wurde zum Testen bei Zenodo hinzugefügt (DOI: 10.5281/zenodo.7975350)14.

Darüber hinaus wurde die Pipeline mit Wurzelhaardaten getestet, wobei die Bildgebung mit einem Lichtblattmikroskop (SPIM) durchgeführt wurde, wie zuvor beschrieben 15,16, wobei Arabidopsis-Wurzelhaare den genetisch kodierten Ca2+ Reporter NES-YC3.6 unter der Kontrolle des UBQ10-Promotors17 exprimierten. Die selbstgebaute Software LabView, die die Kameraaufnahme, die Probenverschiebung und den Shutter des Lichtblattmikroskops steuerte, ermöglichte die Beobachtung der beiden cpVenus- und CFP-Kanäle, aber auch die Visualisierung ihres Verhältnisses in Echtzeit. Jedes Verhältnisbild des Zeitraffers stellte eine Projektion maximaler Intensität (MIP) zwischen den cpVenus- und CFP-Fluoreszenzkanälen dar, die von 15 Schichten der Probe im Abstand von 3 μm aufgenommen wurden. Das Zeitraffer-cpVenus/CFP-Verhältnis der MIPs wurde gespeichert und direkt für die AMEBaS-Analyse verwendet.

Obwohl diese Pipeline mit mehreren Arten von wachsenden und migrierenden Zellen arbeiten kann, wurde sie speziell entwickelt, um wachsende Zellen zu analysieren, die ausschließlich an der Spitze wachsen, wie z. B. Pollenschläuche, Wurzelhaare und Pilzhyphen, bei denen es eine Entsprechung der nicht wachsenden zytoplasmatischen Regionen zwischen den Frames gibt. Wenn eine solche Entsprechung nicht vorhanden ist, sollte der Benutzer die Option complete_skeletonization in Schritt 1.3.1.1 auswählen (weitere Informationen finden Sie im Abschnitt Diskussion).

figure-introduction-1
Abbildung 1: Ein Überblick über den Pipeline-Workflow. Die AMEBaS-Pipeline analysiert und verarbeitet mikroskopische Zeitraffer in drei Hauptschritten: Einzelzellsegmentierung, Midline-Tracing und Kymograph-Generierung. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Access restricted. Please log in or start a trial to view this content.

Protocol

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

1. Interaktives Notebook-Protokoll

Das Jupyter-Notebook kann direkt im Web verwendet werden, indem Google Colab at https://colab.research.google.com/github/badain/amebas/blob/main/AMEBAS_Colab.ipynb verwendet wird, auf dem die folgenden Anweisungen basieren. Alternativ ist das Jupyter-Notebook unter https://github.com/badain/amebas verfügbar, wo es heruntergeladen und für die lokale Ausführung in Jupyter konfiguriert werden kann (Anaconda kann einen einfachen und plattformübergreifenden Installationsprozess bereitstellen). Vollständige Testdaten sind bei Zenodo (https://doi.org/10.5281/zenodo.7975350) zu finden, die ein- und zweikanalige Daten von Arabidopsis-Pollenröhrchen enthalten, die entweder pH- oder Ca2+ -Reporter exprimieren14. Die Pipeline wurde in Teile unterteilt, in denen jeder Schritt durch Klicken auf den Play-Button ausgeführt werden kann, nachdem die benutzerspezifischen Optionen festgelegt wurden. Die für diese Studie notwendigen Dateien sind im AMEBaS-Haupt-ZIP-Ordner (Supplementary Coding File 1) verfügbar.

  1. Öffnen Sie Jupyter Notebook, und lesen Sie die Zeitrafferdateien.
    1. Navigieren Sie zur Startseite des oben genannten interaktiven Notizbuchs in Google Colab, oder laden Sie das Notizbuch AMEBaS_Local.ipynb von GitHub herunter und öffnen Sie es.
    2. Bereiten Sie die Verzeichniseinrichtung für die Ein- und Ausgabedaten vor:
      1. Wenn Sie die lokale Version verwenden, platzieren Sie den Fluoreszenz-Zeitraffer als TIFF-Datei oder DV-Datei in einem Ordner mit dem Namen data , der sich im Stammordner des Programms befinden muss. Es muss ein Ordner mit dem Namen out erstellt werden, um die generierten Daten zu empfangen. Führen Sie dann den Setup-Codeblock aus.
      2. Wenn Sie das Notizbuch in Google Colab verwenden, führen Sie den Setup-Codeblock aus, um die Daten - und Ausgangsordner automatisch zu generieren.
    3. Führen Sie den Codeblock Dateieingabe aus, um die Zeitrafferdaten zu lesen, indem Sie auf die Wiedergabeschaltfläche klicken. Wenn Sie die Google Colab-Version des Notizbuchs verwenden, klicken Sie auf die Schaltfläche Datei auswählen, um die Zeitrafferdatei direkt in den Datenordner hochzuladen.
      HINWEIS: Die Anzahl der Kanäle wird automatisch anhand der Größe des Bildes erkannt.
    4. Wählen Sie aus, ob zusätzliche Ausgaben für jeden Schritt generiert werden sollen, indem Sie den Parameter "verbose" auf "True " oder " False" setzen.
  2. Erkennen Sie die Hauptzelle und das Segment im Hintergrund (Abbildung 2).
    1. Führen Sie den Codeblock Einzelzellensegmentierung aus, um die gewünschte Zelle automatisch vom Hintergrund zu trennen, indem Sie auf die Wiedergabeschaltfläche klicken.
      HINWEIS: Median- und Gauß-Filter werden als Vorverarbeitungsschritt angewendet, um unerwünschtes Rauschen zu entfernen, bevor der Vordergrund vom Hintergrund durch Isodata-Schwellenwerte segmentiert und der Bereich des größten Bereichs isoliert wird, um unerwünschte Artefakte zu entfernen.
      1. Passen Sie den Sigma-Wert an, der von der Gauß-Funktion in der Variablen "Sigma" verwendet wird, um die Glätte der Segmentierungsmaske fein abzustimmen. Der Standardwert ist 2,0.
      2. Legen Sie die Variablenschätzung auf False fest, um den aus Isodata geschätzten Schwellenwert direkt zu speichern, oder auf True , um ihn mithilfe der lokalen Polynomregression (LOESS) über benachbarte Frames hinweg zu glätten. Passen Sie die Funktion an, indem Sie die n_points Variable ändern. Der Standardwert ist 40.
  3. Zeichnen Sie die Mittellinie entlang der Zellverlängerung nach (Abbildung 3).
    1. Führen Sie den Codeblock Cell Midline Tracing aus, indem Sie auf die Wiedergabeschaltfläche klicken, um die Zelle mithilfe der Lee-Methode18 automatisch zu skelettieren und die Spitze des letzten Skeletts durch lineare Extrapolation zu verlängern.
      1. Wählen Sie aus, ob die Mittellinie nur auf dem letzten Frame oder einmal pro Frame nachgezeichnet werden soll, indem Sie das Argument complete_skeletonization anpassen.
        HINWEIS: Wenn alle Frames skelettiert sind, wird die Extrapolation übersprungen.
      2. Legen Sie den Anteil der Punkte im Skelett fest, der während der Extrapolation interpoliert werden soll, indem Sie die Variable interpolation_fraction anpassen. Der Standardwert ist 0,25.
      3. Wählen Sie die Länge der Mittellinienextrapolation aus, indem Sie die Variable extrapolation_length ändern. Der Standardwert ist -1, wodurch das Skelett bis zur nächsten Kante verlängert wird.
  4. Generieren Sie Kymographen für jeden Kanal ( Abbildung 4).
    1. Führen Sie den ersten Codeblock für die Datenvisualisierung aus, indem Sie auf die Schaltfläche " Wiedergabe" klicken, um automatisch Kymographen für beide Kanäle zu generieren.
      1. Wählen Sie die Größe des Gaußschen Kerns, der für die Glättung verwendet wird, indem Sie die Variable kymograph_kernel anpassen.
        HINWEIS: Sie entspricht der Größe der Nachbarschaft (in Pixeln), über die die Pixelintensitäten gemittelt werden. Der Standardwert ist 3 Pixel x 3 Pixel.
      2. Nicht verlängerte Skelette erzeugen gekappte Kymographen, die eine benutzerdefinierte Farbkarte verwenden müssen, um ihre Intensitäten korrekt anzuzeigen. Wählen Sie den prozentualen Anteil der Intensitäten, die der Hintergrundfarbe Schwarz zugewiesen werden, und passen Sie die Variable shift_fraction an. Der Standardwert ist 0,7.
  5. Berechnen Sie das Verhältnis zwischen den Kanälen (Abbildung 5).
    1. Führen Sie den zweiten Codeblock für die Datenvisualisierung aus, indem Sie auf die Wiedergabeschaltfläche klicken, um automatisch einen ratiometrischen Kymographen und einen ratiometrischen Zeitraffer zu generieren (Abbildung 6).
      HINWEIS: Dieser Schritt ist nur verfügbar, wenn Sie Dual-Channel-Zeitraffer verwenden. Der in Schritt 1.2.1.2 gespeicherte Schwellenwert für die Hintergrundintensität wird von jedem Kanal subtrahiert.
      1. Passen Sie die Variable switch_ratio an, um die Reihenfolge der Kanäle zu ändern, die bei der Berechnung des Verhältnisses als Zähler und Nenner verwendet werden. Der Standardwert ist False.
      2. Wählen Sie aus, ob das Zeitrafferverhältnis mit einem Median-Filterdurchlauf weiter geglättet werden muss, indem Sie die Variable smooth_ratio anpassen. Der Standardwert ist False.
      3. Wählen Sie aus, ob Ausreißer, die durch das Low-Signal des Nennerkanals erzeugt werden, durch Manipulation der Variablen reject_outliers entfernt werden sollen. Der Standardwert ist True und definiert Ausreißer als Werte, die dem 1,5-fachen des Interquartilsabstands über dem dritten Quartil (wo 75 % der Werte liegen) entsprechen.
      4. Wählen Sie aus, ob der Hintergrund in der ratiometrischen Ausgabe exportiert werden soll, indem Sie die Variable background_ratio anpassen. Der Standardwert ist False, wodurch er durch Nullen ersetzt wird.

figure-protocol-1
Abbildung 2: Schritt der Einzelzellsegmentierung. Bildverarbeitungstechniken wie Filterung, Schwellwertbildung und Bereichsbeschriftung werden verwendet, um das interessierende Signal zu isolieren (Schritt 1.2). Diese speziellen Daten enthielten die folgenden Werte für die niedrigste Intensität: 2556, den Median: 3441 und die höchste Intensität: 32125. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

figure-protocol-2
Abbildung 3: Übersicht über die Mittellinienverfolgung - Die Mittellinie der einzelnen Zelle wird durch Berechnung ihres Skeletts (weiß) ermittelt. Die Spitze (Magenta) wird linear von den letzten Punkten am Ende des Skeletts extrapoliert (Schritt 1.3). In dieser Komposition werden sowohl die Mittellinie als auch ihre Spitze über die ursprüngliche Zelle gelegt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

figure-protocol-3
Abbildung 4: Zeitraffer-Kymographen - Vergleich der Kymographen für jeden Kanal, der mit ausgeschaltetem "complete_skeletonization" erzeugt wurde (Schritt 1.4). Die vertikale Achse beschreibt den Verlauf der Zeit, und die horizontale Achse stellt die durchschnittliche Intensität des extrapolierten Mittellinienpfads dar, dem eine einzelne Zelle folgt. Für diese speziellen Daten stellt die Farbkarte die folgenden Werte für Kanal 1 (Niedrigste Intensität): 2886, Median: 3167, Höchste Intensität: 21021 dar. Kanal 2 Niedrigste Intensität: 3030, Median: 3400, Höchste Intensität: 29688. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

figure-protocol-4
Abbildung 5: Glättung des Hintergrundschwellenwerts - Der Schwellenwert für die Hintergrundsegmentierung wird über den Isodata-Algorithmus geschätzt und dann über die lokale polynomiale Regression geglättet (Schritt 1.5). Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

figure-protocol-5
Abbildung 6: Ratiometrische Ergebnisse- (A) Vergleich zwischen dem letzten Frame des ratiometrischen Zeitraffers und dem segmentierten ursprünglichen ersten Kanal. (B) Kymograph, der aus dem ratiometrischen Zeitraffer (Schritt 1.5) generiert wurde. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

2. Protokoll für den Batch-Modus

  1. Laden Sie die Datei pipeline.py aus dem AMEBaS GitHub-Repository herunter und legen Sie sie im selben Verzeichnis wie die Daten ab.
  2. Geben Sie den Speicherort der Datei in der Befehlszeile nach der Programmdatei ein.
  3. Fügen Sie - -v als Positionsargument ein, um die internen Schritte der Pipeline anzuzeigen, falls gewünscht.
  4. Schließen Sie --s ein, um den Sigma-Wert auszuwählen, der im Vorverarbeitungsschritt des Gaußschen Filters in Vorbereitung auf die Zellsegmentierung verwendet wird. Der Standardwert ist 2.
  5. Fügen Sie - -a ein, um die Mittellinie für jedes Bild des Zeitraffers nachzuzeichnen. Standardmäßig verwendet die Pipeline nur den letzten Frame.
  6. Fügen Sie - -f ein, um den Bruchteil [0,1] des Skeletts auszuwählen, der in der Interpolation verwendet werden soll. Der Standardwert ist 0,25.
  7. Fügen Sie -e ein, um die Länge des extrapolierten Skeletts in Pixeln auszuwählen. Der Standardwert ist -1, wodurch das Skelett bis zur nächsten Kante verlängert wird.
  8. Fügen Sie - -sf ein, um den Anteil des Farbbereichs auszuwählen, der in nicht extrapolierten Kymographen in den Hintergrund verschoben wird. Der Standardwert ist 0,7.
  9. Fügen Sie - -k ein, um die Größe des Kernels zu bestimmen, der in der Gaußschen Filterung des kymographen verwendet wird. Der Standardwert ist 3.
  10. Fügen Sie --eb ein, um die globale Intensität des Hintergrundschwellenwerts über die polynomiale LOESS-Regression der framespezifischen Hintergrundschwellenwertintensitäten zu schätzen.
    1. Passen Sie die Anzahl der Punkte an, die bei der LOESS-Glättung der Hintergrundschwellenwerte verwendet werden, indem Sie den Parameter - -n ändern. Der Standardwert ist 40.
  11. Wechseln Sie die Kanäle, die als Zähler und Nenner bei der Berechnung des Verhältnisses verwendet werden, einschließlich - r oder - -switch_ratio, wenn der Zeitraffer zwei Kanäle hat. Standardmäßig ist der zweite Kanal der Zähler und der erste der Nenner.
  12. Wählen Sie aus, ob das Zeitrafferverhältnis mit einem Medianfilter-Durchlauf mit dem Argument - -sm weiter geglättet werden muss. Der Standardwert ist False.
  13. Schließen Sie -o ein, um Pixel mit abnormaler Intensität während der ratiometrischen Zeitraffererzeugung abzulehnen.
  14. Wählen Sie mit dem Argument - -b aus, ob der Hintergrund in der ratiometrischen Ausgabe exportiert werden soll. Der Standardwert ist False, wodurch er durch Nullen ersetzt wird.
  15. Drücken Sie die Eingabetaste, um auszuführen. Die Ausgabe wird im selben Verzeichnis wie die Programmdatei erzeugt.

Access restricted. Please log in or start a trial to view this content.

Results

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Die AMEBaS-Pipeline automatisiert die Extraktion der Mittelliniendynamik polarisierter Einzelzellen aus fluoreszenzmikroskopischen Bildstapeln, wodurch sie weniger zeitaufwändig und weniger anfällig für menschliche Fehler ist. Die Methode quantifiziert diese Zeitspanner durch die Erzeugung von Kymographen und ratiometrischen Bildstapeln (Abbildung 1) in wachsenden Einzelzellen. Es kann so eingestellt werden, dass es an der Migration einzelner Zellen arbeitet, aber weitere Experimente sind no...

Access restricted. Please log in or start a trial to view this content.

Discussion

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Die hier vorgestellte neuartige Methode ist ein leistungsfähiges Werkzeug, um die Analyse von fluoreszenzmikroskopischen Bildstapeln polarisierter Zellen zu rationalisieren und zu automatisieren. Aktuelle Methoden, die in der Literatur beschrieben werden, wie z. B. ImageJ Kymograph-Plugins, erfordern eine manuelle Verfolgung der Mittellinie der polarisierten Zelle von Interesse, eine Aufgabe, die nicht nur zeitaufwändig, sondern auch anfällig für menschliche Fehler ist. Da die Definition der Mittellinie in dieser Pipelin...

Access restricted. Please log in or start a trial to view this content.

Disclosures

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Die Autoren dieses Manuskripts erklären keine konkurrierenden finanziellen Interessen oder andere Interessenkonflikte.

Acknowledgements

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Die Autoren bedanken sich für die FAPESP-Zuschüsse 2015/22308-2, 2019/23343-7, 2019/26129-6, 2020/06744-5, 2021/05363-0, CNPq, NIH R01 Grant GM131043 und die NSF-Zuschüsse MCB1714993, MCB1930165 für die finanzielle Unterstützung. Die Wurzelhaardaten wurden mit der Infrastruktur und unter der Aufsicht von Prof. Andrea Bassi und Prof. Alex Costa erstellt.

Access restricted. Please log in or start a trial to view this content.

Materials

List of materials used in this article
NameCompanyCatalog NumberComments
GithubGithubhttps://github.com/badain/amebas
Google ColabGooglehttps://colab.research.google.com/github/badain/amebas/blob/main/AMEBAS_Colab.ipynb

References

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Drubin, D. G., Nelson, W. J. Origins of cell polarity. Cell. 84 (3), 335-344 (1996).
  2. Wodarz, A., Näthke, I. Cell polarity in development and cancer. Nature Cell Biology. 9 (9), 1016-1024 (2007).
  3. Palanivelu, R., Preuss, D. Pollen tube targeting and axon guidance: parallels in tip growth mechanisms. Trends in Cell Biology. 10 (12), 517-524 (2000).
  4. Portes, M. T., et al. The Pollen Tube Oscillator: Integrating Biophysics and Biochemistry into Cellular Growth and Morphogenesis. Rhythms in Plants: Dynamic Responses in a Dynamic Environment. , Springer. Cham. (2015).
  5. Wudick, M. M., et al. CORNICHON sorting and regulation of GLR channels underlie pollen tube Ca2+ homeostasis. Science. 360 (6388), 533-536 (2018).
  6. Hoffmann, R. D., et al. Plasma membrane H+-ATPases sustain pollen tube growth and fertilization. Nature Communications. 11 (1), 1-15 (2020).
  7. Michard, E., et al. Glutamate receptor-like genes form Ca2+ channels in pollen tubes and are regulated by pistil D-serine. Science. 332 (6028), 434-437 (2011).
  8. Damineli, D. S., Portes, M. T., Feijó, J. A. Oscillatory signatures underlie growth regimes in Arabidopsis pollen tubes: computational methods to estimate tip location, periodicity, and synchronization in growing cells. Journal of Experimental Botany. 68 (12), 3267-3281 (2017).
  9. Li, K., et al. An optimized genetically encoded dual reporter for simultaneous ratio imaging of Ca2+ and H+ reveals new insights into ion signaling in plants. New Phytologist. 230 (6), 2292-2310 (2021).
  10. Sadoine, M., et al. Designs, applications, and limitations of genetically encoded fluorescent sensors to explore plant biology. Plant Physiology. 187 (2), 485-503 (2021).
  11. Grenzi, M., et al. Illuminating the hidden world of calcium ions in plants with a universe of indicators. Plant Physiology. 187 (2), 550-571 (2021).
  12. Munglani, G., Vogler, H., Grossniklaus, U. Fast and flexible processing of large FRET image stacks using the FRET-IBRA toolkit. PLoS Computational Biology. 18 (4), 1009242(2022).
  13. Ridler, T., Calvard, S. Picture thresholding using an iterative selection method. IEEE Transactions on Systems, Man, and Cybernetics. 8 (8), 630-632 (1978).
  14. Portes, M. T., Feijó, J. Growing Arabidopsis pollen tubes expressing genetically encoded reporters for calcium and pH. , (2023).
  15. Candeo, A., Doccula, F. G., Valentini, G., Bassi, A., Costa, A. Light sheet fluorescence microscopy quantifies calcium oscillations in root hairs of Arabidopsis thaliana. Plant & Cell Physiology. 58 (7), 1161-1172 (2017).
  16. Romano Armada, N., et al. In vivo light sheet fluorescence microscopy of calcium oscillations in Arabidopsis thaliana. Methods in Molecular Biology. 1925, 87-101 (2019).
  17. Krebs, M., et al. FRET-based genetically encoded sensors allow high-resolution live cell imaging of Ca2+ dynamics. The Plant Journal. 69 (1), 181-192 (2012).
  18. Lee, T. -C., Kashyap, R. L., Chu, C. -N. Building skeleton models via 3-D medial surface axis thinning algorithms. CVGIP: Graphical Models and Image Processing. 56 (6), 462-478 (1994).
  19. Nunez-Iglesias, J., Blanch, A. J., Looker, O., Dixon, M. W., Tilley, L. A new Python library to analyse skeleton images confirms malaria parasite remodelling of the red blood cell membrane skeleton. PeerJ. 6, 4312(2018).

Access restricted. Please log in or start a trial to view this content.

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

Tags

Cell PolarityRatiometric FluorescenceMidline ExtractionBackground SubtractionSingle Cell AnalysisFluorescence Time LapseCell SegmentationSkeletonizationKymograph GenerationPolarized Cells

Related Articles