The authors are grateful to FAPESP grants 2015/22308-2, 2019/23343-7, 2019/26129-6, 2020/06744-5, 2021/05363-0, CNPq, NIH R01 grant GM131043 and the NSF grants MCB1714993, MCB1930165 for financial support. Root hair data were produced with the infrastructure and under the supervision of Prof. Andrea Bassi and Prof. Alex Costa.

1. Interactive notebook protocol
The Jupyter notebook can be used directly on the web using Google Colab at https://colab.research.google.com/github/badain/amebas/blob/main/AMEBAS_Colab.ipynb, where the instructions below were based. Alternatively, the Jupyter notebook is available at https://github.com/badain/amebas, where it can be downloaded and configured to run locally in Jupyter (Anaconda can provide an easy and cross-platform installation process). A complete test data ca...

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

The novel method presented here is a potent tool to streamline and automate the analysis of fluorescence microscopy image stacks of polarized cells. Current methods described in the literature, such as ImageJ Kymograph plugins, require manual tracing of the midline of the polarized cell of interest, a task that is not only time consuming but also prone to human errors. Since the definition of the midline in this pipeline is supported by a numerical method18,19 th...

Cell polarity is a fundamental biological process in which the concerted action of a collection of spatially concentrated molecules and structures culminates in the establishment of specialized morphological subcellular domains1. Cell division, growth and migration rely on such polarity sites, while its loss has been associated with cancer in epithelial tissue-related disorders2.
Apically growing cells are a dramatic example of polarity, where the polarity site at the tip typically reorients to extracellular cues3. These include developing neurites, fungal hyphae, root hairs, and pollen tubes, where multiple cellular processes show pronounced differences from the tip of the cell toward the shank. In pollen tubes, in particular, actin polymerization, vesicle trafficking, and ionic concentrations are markedly polarized, showing tip-focused gradients4. Pollen tubes are the male gametophytes of flowering plants and are responsible for delivering the sperm cells to the ovule by growing exclusively at the apex of the cell at one of the fastest growth rates known for a single cell. The tip-focused gradients of ions such as calcium5 (Ca2+) and protons6 (H+) play a major role in sustaining pollen tube growth, which is essential to accomplish its main biological function that culminates in a double fertilization5,6. Thus, quantitative methods to analyze the spatiotemporal dynamics along the midline of apically growing cells are essential to investigate the cellular and molecular mechanisms underlying polarized growth7,8,9. Researchers often use kymographs, i.e., a matrix that represents the pixel intensities of the cell&#39;s midline (e.g., columns) through time (e.g., rows), which allows visualizing cell growth and migration in the diagonal (Figure 1). Despite their usefulness, kymographs are frequently extracted by manually tracing the midline, being prone to biases and human errors while also being rather laborious. This calls for an automated method of midline extraction that is the first feature of the pipeline introduced herein named AMEBaS: Automatic Midline Extraction and Background Subtraction of ratiometric fluorescence time lapses of polarized single cells.
In terms of experimental procedures, quantitative imaging of ions/molecules/species of interest in single cells can be achieved with genetically encoded fluorescent probes10. Among the ever-expanding choices, ratiometric probes are one of the most accurate since they emit different fluorescence wavelengths when bound/unbound to the molecules of interest11. This allows for correction of the spatial heterogeneity in the intracellular concentration of the probe by using the ratio of two channels with their channel specific background subtracted. However, estimating the background threshold for each channel and time point can be a complex task since it often varies in space due to effects like shading, where the corners of the image have luminosity variation relative to the center, and in time due to fading of the fluorophore (photobleaching)12. Although there are multiple possible methods, this manuscript proposes determining the background intensity automatically using the segmentation threshold obtained with the Isodata algorithm13, which is then smoothed across frames through polynomial regression as a standard. Spatial components stemming from fluorescence heterogeneity unrelated to the target cell removed in12, however, were ignored by this method. Automatic thresholding can be performed by several methods, but the Isodata algorithm produced the best results empirically. Thus, automatic background value subtraction and ratiometric calculation are the second main feature of AMEBaS (Figure 1), which, taken together, receives as input a stack of dual-channel fluorescence microscopy images, estimates the cell&#39;s midline and the channel-specific background, and outputs kymographs of both channels and their ratio (main output #1) after background subtraction, smoothing, and outlier removal, together with a stack of ratiometric images (main output #2).
AMEBaS was tested with fluorescence time lapses of growing Arabidopsis pollen tubes obtained under a microscope, either with Ca2+ (CaMeleon)8 or pH (pHluorin)6 ratiometric sensors expressed under the pollen-specific LAT52 promoter. Images from each channel were taken every 4 s coupled to an inverted microscope, a front-illuminated camera (2560 pixels &#215; 2160 pixels, pixel size 6.45 &#956;m), a fluorescence illuminator, and a water immersion objective lens 63x, 1.2NA. Filters settings used for CaMeleon were: excitation 426-450 nm (CFP) and 505-515 nm (YFP), emission 458-487 nm (CFP) and 520-550 nm (YFP), while for pHluorin, excitation 318-390 nm (DAPI) and 428-475 nm (FITC), emission 435-448 nm (DAPI) and 523-536 nm (FITC). A complete data set was added for testing at Zenodo (DOI: 10.5281/zenodo.7975350)14.
In addition, the pipeline was tested with root hair data, where imaging was performed with a light sheet microscope (SPIM) as previously described15,16 with Arabidopsis root hairs expressing the genetically encoded Ca2+ reporter NES-YC3.6 under the control of the UBQ10 promoter17. The home-made LabView software that controlled the camera acquisition, sample translation and shutter of the light sheet microscope permitted the observation of the two cpVenus and CFP channels, but also the visualization of their ratio in real time. Every ratio image of the time-lapse represented a maximum intensity projection (MIP) between the cpVenus and CFP fluorescent channels images obtained from 15 slices of the sample spaced 3 &#956;m apart. The time-lapse cpVenus/CFP ratio of MIPs was saved and directly used for the AMEBaS analysis.
Although this pipeline can work with multiple types of growing and migrating cells, it was specifically designed to analyze growing cells that grow exclusively at the tip, such as pollen tubes, root hairs, and fungal hyphae, where there is a correspondence of the non-growing cytoplasmic regions between frames. When such a correspondence is not present, the user should choose the complete_skeletonization option in step&#160;1.3.1.1 (see the Discussion section for more details).
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/64857/64857fig01.jpg" /> 
Figure 1: An overview of the pipeline workflow. The AMEBaS pipeline analyses and processes microscopic time lapses in three main steps: Single-Cell Segmentation, Midline Tracing, and Kymograph Generation. <a href="https://www.jove.com/files/ftp_upload/64857/64857fig01large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Github</td>
			<td>Github</td>
			<td></td>
			<td>https://github.com/badain/amebas</td>
		</tr>
		<tr>
			<td>Google Colab</td>
			<td>Google</td>
			<td></td>
			<td>https://colab.research.google.com/github/badain/amebas/blob/main/AMEBAS_Colab.ipynb</td>
		</tr>
	</tbody>
</table>

amebas: automatic midline extraction and background subtraction of ratiometric fluorescence time-lapses of polarized single cells

Cell polarity is a macroscopic phenomenon established by a collection of spatially concentrated molecules and structures that culminate in the emergence of specialized domains at the subcellular level. It is associated with developing asymmetric morphological structures that underlie key biological functions such as cell division, growth, and migration. In addition, the disruption of cell polarity has been linked to tissue-related disorders such as cancer and gastric dysplasia.
Current methods to evaluate the spatiotemporal dynamics of fluorescent reporters in individual polarized cells often involve manual steps to trace a midline along the cells&#39; major axis, which is time consuming and prone to strong biases. Furthermore, although ratiometric analysis can correct the uneven distribution of reporter molecules using two fluorescence channels, background subtraction techniques are frequently arbitrary and lack statistical support.
This manuscript introduces a novel computational pipeline to automate and quantify the spatiotemporal behavior of single cells using a model of cell polarity: pollen tube/root hair growth and cytosolic ion dynamics. A three-step algorithm was developed to process ratiometric images and extract a quantitative representation of intracellular dynamics and growth. The first step segments the cell from the background, producing a binary mask through a thresholding technique in the pixel intensity space. The second step traces a path through the midline of the cell through a skeletonization operation. Finally, the third step provides the processed data as a ratiometric timelapse and yields a ratiometric kymograph (i.e., a 1D spatial profile through time). Data from ratiometric images acquired with genetically encoded fluorescent reporters from growing pollen tubes were used to benchmark the method. This pipeline allows for faster, less biased, and more accurate representation of the spatiotemporal dynamics along the midline of polarized cells, thus advancing the quantitative toolkit available to investigate cell polarity. The AMEBaS Python source code is available at:&#160;https://github.com/badain/amebas.git

The AMEBaS pipeline automates the extraction of midline dynamics of polarized single cells from fluorescence microscopy image stacks, making it less time consuming and less prone to human errors. The method quantifies these time lapses by generating kymographs and ratiometric image stacks (Figure 1) in growing single cells. It can be adjusted to work on migrating single cells, but further experiments are necessary. AMEBaS is implemented in Python as an interactive Jupyter Notebook (described...

Watch this Scientific Journal Video about AMEBaS: Automatic Midline Extraction and Background Subtraction of Ratiometric Fluorescence Time-Lapses of Polarized Single Cells at JoVE.com

AMEBaS: Automatic Midline Extraction and Background Subtraction of Ratiometric Fluorescence Time-Lapses of Polarized Single Cells

Current methods for analyzing the intracellular dynamics of polarized single cells are often manual and lack standardization. This manuscript introduces a novel image analysis pipeline for automating midline extraction of single polarized cells and quantifying spatiotemporal behavior from time lapses in a user-friendly online interface.

Cell polarity is a macroscopic phenomenon established by a collection of spatially concentrated molecules and structures that culminate in the emergence ...

amebas-automatic-midline-extraction-background-subtraction

Research

Education

JoVE Core

JoVE Journal

Electrical Engineering

Developmental Biology

Unit 1

Sampling

SCIENTISTS IN ACTION

343 조회수.  University of São Paulo. Current methods for analyzing the intracellular dynamics of polarized single cells are often manual and lack standardization. This manuscript introduces a novel image analysis pipeline for automating midline extraction of single polarized cells and quantifying spatiotemporal behavior from time lapses in a user-friendly online interface.