DOI: 10.5061/dryad.76hdr7t64
title: Data from: Phagocytosed photoreceptor outer segment particles within the retinal pigment epithelium show diurnal rhythmicity and variation between cone subtypes in larval zebrafish

description: During phagocytosis of retinal rod and cone outer segments (OSs) by the
 retinal pigment epithelium (RPE), RPE bites off and internalizes the
 oldest parts of the OS tips to prevent the accumulation of harmful
 compounds in the photoreceptors. Increases in the number of phagocytosed
 OS particles in the RPE appear rhythmically once or twice a day, depending
 on the animal species. However, the variation of this rhythmicity between
 the distinct photoreceptor types is not well understood. To compare the OS
 phagosome numbers and their daily rhythms between the different cone
 subtypes, we used zebrafish larvae as an animal model. We collected
 zebrafish larvae at seven different time points during a 24 h circadian
 cycle and prepared histological sections of their eyes. We immunolabelled
 the different cone opsins together with the RPE-specific proteins in the
 histological sections. The sections were then imaged with a laser scanning
 confocal microscope, after which the internalized OS particles were
 quantified from the images. To make the quantification of the OS particles
 more efficient, we developed a semi-automated analysis tool for
 ImageJ-Fiji. The analysis tool allows the user to select the RPE area from
 which it then segments and counts the detected OS particles. After
 quantification, the results revealed the presence of OS particles of all
 cone subtypes in the RPE throughout the day in larval zebrafish. However,
 the phagosome numbers from UV and blue cone OSs increased significantly at
 two time points, whereas the number of green cone OS particles was more
 constant. The numbers of red cone OS particles were continuously lower
 compared to other cone types and showed only one increase in the latest
 time point during the dark phase of the day. We also investigated whether
 the rhythmicity is affected by external light by keeping the larvae in
 constant darkness for at least 24 h before sample preparation. We found
 that the complete darkness condition decreased the phagosome numbers of
 all cone subtypes, and additionally, the daytime peaks in the UV and blue
 cones disappeared. This indicates that the rhythmicity is strongly
 affected by the external light in the larval zebrafish. Our findings
 provide new understanding on the cone OS phagocytosis and its rhythmicity.
description: # Data from: Phagocytosed photoreceptor outer segment particles within the
 retinal pigment epithelium show diurnal rhythmicity and variation between
 cone subtypes in larval zebrafish Dataset DOI:
 [10.5061/dryad.76hdr7t64](10.5061/dryad.76hdr7t64) ## Description of the
 data and file structure This is accompanying data for the research article
 “Phagocytosed photoreceptor outer segment particles within the retinal
 pigment epithelium show diurnal rhythmicity and variation between cone
 subtypes in larval zebrafish” compressed into Dataset_.zip folder. We have
 included raw and processed versions of the images presented in the
 published article and its supplementary files. Additionally, we provide
 the full script of our developed semi-automatized analysis tool for
 imageJ-Fiji macro used for data analysis, as well as numerical data on the
 OS phagosome quantification used to generate the graphs in the main
 article and in the supplementary files. ### Files and variables The
 dataset inside compressed archive named as Dataset_.zip is organized into
 three folders as follows: 1. Imaging data, 2. Quantitative analysis tool,
 3. Quantitative data. ### 1.Imaging data The data in this folder is
 further organized into the subfolders “Raw images” and “Processed images”
 which are organized even further into folders named after the figures of
 the main article and the supplementary material (Figure 1, Figure 2,
 Supplementary Figure S1, etc.). Raw images are in the form of .ND2 files
 and processed images as .tiff files. Imaging data folder also contains a
 .docx file with more detailed information on the naming of the image
 files, used microscope and laser lines, as well as image processing.
 Abbreviations in the file names: UV=UV cones, Blue= Blue cones,
 Green=Green cones, Red=Red cones, ZPR2= antibody that labels RPE tissue,
 RPE65=antibody that labels RPE tissue, ZPR3=antibody that labels green
 cones, UV opsin= antibody that labels UV cones, Blue opsin= antibody that
 labels blue cones, Red opsin= antibody that labels red cones, ZT=zeitgeber
 time ### 2.Quantitative analysis tool This folder contains files that
 provide the full script of the analysis tool for quantitative analysis of
 outer segment (OS) phagosomes in the RPE tissue as .ijm file and .docx
 file. The tool is for ImageJ-Fiji (ImageJ 1.54), and it utilizes standard
 ImageJ-Fiji libraries for general image processing and MorphoLibJ
 (MorphoLibJ_-1.6.2) for all morphological operations. The full script of
 the analysis tool .docx file also contains information on the description
 and key features of the analysis tool. The entire workflow is performed as
 an interactive script that allows the user to configure the used channels,
 adjust parameters and select RPE area from which the OS phagosomes are
 then quantified by the tool. The analysis tool uses confocal z-stack
 images of fluorescently labelled cryosections containing RPE and
 photoreceptor channels as input. The output of the analysis (segmented RPE
 area and detected OS phagosomes) are saved as ROI files and .xlsx files.
 Before further analysis, the stacks are processed into maximum intensity
 projections and Gaussian blur filter with radius 1 is applied. The
 workflow of the analysis tool is divided into two steps: Manual RPE area
 selection by the user and OS phagosome segmentation. OS phagosome
 segmentation is based on the extraction of local bright objects using the
 h-dome transform to detect all local maxima. h-dome transform method
 suppresses low-contrast regions and enhances significant intensity peaks
 bigger than defined h-value. The Maxima Finder is then applied to detect
 the intensity peaks with an additional criterion called prominence that is
 the minimum height difference between a local maximum and its surrounding
 region for the maximum to be considered significant. In addition, the
 folder contains .xlsx file with different numerical parameters used for
 the analysis tool development and testing needed for the evaluation of the
 performance of the analysis tool's peak detection algorithm. A
 randomly chosen set of Z-stacks with each cone subtype (UV, blue and
 green) labelled together with the RPE tissue was used. Empirical testing
 was conducted to determine appropriate values for the parameters
 "h-value" and "prominence" to find the balance between
 sensitivity (detecting enough intensity peaks) and specificity (avoiding
 false positives or irrelevant peaks) and the values are included in
 separate worksheets for each cone subtype (named as Parameter
 testing_UVcones, Parameter testing_Bluecones, Parameter testing_Green
 cones). The test analysis for each image was performed with the analysis
 tool using all combinations of five h-dome values (5, 15, 25, 50, 100) and
 five prominence thresholds (5, 10, 20, 40, 80). The local maxima detected
 by the analysis tool were then matched to the manually annotated local
 maxima (Ground truth) and the detected peaks were defined as True positive
 (TP), False positive (FP) and False negative (FN) (see the definitions in
 the table below). By using the counts of TP, FP and FN, the values for
 Precision, Sensitivity, F1-score and False positive rate (FPR) (see the
 definitions and equations in the table below) were calculated for each
 parameter combination of h-values and prominence thresholds to evaluate
 the algorithm’s performance and its overall detection reliability in the
 current analysis. The "Plotted values" worksheet in the .xlsx
 file contains the average values and standard deviations (SD) of UV, blue
 and green cones for each parameter (F1-score, FPR, Precision and Recall)
 with each prominence and h-value combinations. The parameter combinations
 of all five h-dome values and prominence thresholds between 5-20 produced
 relatively similar values for precision, sensitivity, F1-score and FPR,
 the performance and peak detection reliability with these parameters were
 highly equal. Therefore, in our analysis tool, the h-value was set to 15,
 whereas the prominence threshold was set specifically, always under the
 value of 20, for each analysed image during the OS phagosome
 quantification to reach the most reliable detection of OS phagosomes.
 Missing values in this .xlsx file are marked as NaN. | **Parameter** |
 **Definition** | **Equation** | | :------------------------ |
 :---------------------------------------------------------------------------------------- | :------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ | | Precision | The proportion of detected peaks that were correct | ![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAANsAAABNCAMAAAAYYuEsAAAAAXNSR0IArs4c6QAAAK5QTFRFAAAAAAAAAAA6AABmADo6ADpmADqQAGa2OgAAOgA6OgBmOjoAOjpmOjqQOmZmOmaQOma2OpC2OpDbZgAAZgA6ZgBmZjoAZjo6ZjpmZmY6ZpC2ZpDbZrbbZrb/kDoAkDpmkGYAkGY6kLbbkLb/kNv/tmYAtmY6tpA6tpBmtrZmttvbttv/tv//25A625Bm27Zm27aQ29u22//b2////7Zm/9uQ/9u2/9vb//+2///bg4IXPQAAAAF0Uk5TAEDm2GYAAAAJcEhZcwAAEnQAABJ0Ad5mH3gAAAAZdEVYdFNvZnR3YXJlAE1pY3Jvc29mdCBPZmZpY2V/7TVxAAAED0lEQVR4Xu1a7XbTMAyNt3UEGNAS6GAbAcY6vpa10KYk7/9iyJLSyPEOdbflODX2j6bREs9XsmRdqUkSR9RA1IC7Bgq1GU/grRndnVy6zzDcJ2cHV0mZHi+TlXoNq6w/6cs6Ux+Gu2TXlVUZQEFY5dMr/ZLGqiXaiHs+ymc3eiOClfBbUueH+hIENm2aOkdb4ShTNFiBGzSAUWXgbjzI66qsRbvfAOUGLDSoxTSEUIJGERuwzvUJcPDq535bq109hUYccnuGAI9DI0Ihdwtn4MnNowjG0QiQtJXYnkEYD0/uxt3w5A5mQO543oCZq8NNWAkAYJVB0OcAstLxPyRwAdinPwhBUb6OmnaifFXWiVaWoD8j3GfmXSgf84z231iC+6ygt3dCo3xSUaFRPolNUr7bFM6T9VShUy1OlRrRkTJ/DsfLGJgv/YXuT760ApTgw80U/GZv281pYkH5ivM6P/5+tpxp8n6tJkvKVYHYXwAMXZQhWr+CdKH+ipUMEsATl8T/i7Pf+fF8vFzRH+hpMdrc12lxD3vIpHxwN6b5kBUCVMCGWV2VIRRME2baetVb+GgEmPdxqS0/gikGEWNMylemHOPJYrhCQTSYZ8zUBSuUBPwEYaMptrEtact+vm+U3e4fzt+02f58S0ewTEGc0GDahoqpPQmY+dOFUA2CbZmLaEgS7tSjl7o8LcpoGxvXn9WBNhgL+C26oP+ahNmXvxmUb7MksTbxVTgR7TwS8BMElG4GUdwwc6bWMFTDLV8spQnIxhhEqIDd+B8ajIgkTbHN3R4W/xzfNilfu6RCh4s51uHRjX6dwsIByvrdzQoOh6RECCQAJJMlnAkYYWkK+Fx/9Ey3OpRPOB+cwVwLrK+BFp780Gac6hhSvk9BMNZFGhKAEkDCvSO2ba5GDoXE9enE0QbOj92OhtHDKiAnANvfp8vGzExhomTQtMUU1e55YERFx929y0bMLKmmGMyMzlyV+W9irTj2uHTZCqvI1Zw7lCgZnTlOljxajvI56q1t7bJZ2OQBanXmOMXwB05od3uXzcImK85WZ877CUTehmN7l83Chsuv3tAU3c6cb48zs7LW++/uslnY8KhppN13jA3rYWdKW23vsnWxETNTxAbtzpzn3oHE9q8u293MFlPFOifSYievw8Hm0GXr2o16xTPqhNihYzjYHLpsXWzGjxVsrjgcbA5dti42Y/E2Es/YxEZ06LJ1sBnMzC5t+46TxM/1cOmydbAZzMzuzHlnxs1OdOqymdgMZnZHZ857XsL5ZB9Hq/d8EkJ3Tz8W9M8DuCL2+HbznU0SIuTdjz0W6RB4N4TIcOslj22yOF/UQNRA1EDUQNRA1EDUQNRA1EDUQNRA1EDUQNRA1EDUwH+qgb+OqXrpAxA9LQAAAABJRU5ErkJggg==) | | Ground truth | The manually annotated local intensity maxima |   | | Sensitivity (Recall) | The proportion of ground truth peaks that were successfully detected | ![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAOsAAABNCAMAAABwkWMUAAAAAXNSR0IArs4c6QAAALRQTFRFAAAAAAAAAAA6AABmADo6ADpmADqQAGa2OgAAOgA6OgBmOjoAOjo6OjpmOjqQOmZmOmaQOma2OpC2OpDbZgAAZgA6ZgBmZjoAZjo6ZjpmZjqQZmY6ZpC2ZpDbZra2ZrbbZrb/kDoAkDo6kDpmkGYAkGY6kLbbkLb/kNv/tmYAtmY6tpA6tpBmttv/tv/btv//25A625Bm27Zm27aQ2//b2////7Zm/9uQ/9u2/9vb//+2///buWmIkwAAAAF0Uk5TAEDm2GYAAAAJcEhZcwAAEnQAABJ0Ad5mH3gAAAAZdEVYdFNvZnR3YXJlAE1pY3Jvc29mdCBPZmZpY2V/7TVxAAAFAUlEQVR4Xu1aaXuaQBB2NWnokeay1fRKaJtWbNOgVqDw//9X51pYMNEV1FKe3Q/q6gL7zrXvzNjrueEk4CRwMAmEKh8n8NCAZ6ffDraBAz4o6E96sfds2YvUFTw2+4pvyVDdHnAPB3pUOgRoBDN+PsFnInb8BpXcsRG/eEDDBS3Sp17mD/Ctk1hRdZlPuqQRe6TQkAy6gyMdgrvKYK9NhwX6bgE2DTZEkItRF0MTKc0w2MzHE6d/ft8tbRZoOPTSMM25i3Al9BI0dtfuDmISMsLOOioDNHVpmHMnlUtMQrsrMYnODuC+HzW4uRrkYaqDgNMhHDISkCI8b7oMtoP62yek5A4M4qi7BMSQXeRdLEsJxFq5psNS/KtMmdC0NkIK14qKU3otVsmk0iFnyDI1L6ms2KdFbnvverlutJG1bF6x7U6br4+UpUbNRyVGzvz4FjavaL71re8Apxb4qx6LsVLHcG7NPDjJkhF9hjF/iRW772CyCp2RS3j9iUyzKZ11M69/u7JiRKKUhVtvbtcXRJ461kngVF0ukWCHH377z+YXy4jrHBHwlewHlXr4i5K7ptf3SMSnV5kPP5ZWJO+WzNELT0EWkI8aJtUMPqhPKAilwpmPO8j8o4t84wEqM30LL7FHxEWTcJlCDo2AqNJXWSEoo7aQnCxgn+WUiRUTe4hPMAXqRsQpmZROqPQ7Jx1U36uskJw6WF/YNHW9r8+CAOoZqBdU65+f3jGXIlGBsnHwaaEapGFdrgTH5aNUcugQL6ms4N/SN+3hroSJKjhHZ9wzoMpVXgfI7hTEnbzwoesf5ff4FC2dK31GhQR1HhZq/Xf+KhSCbNCscPBnY8ds0xJ4NIMovWefyY45eBkcA+QohfpmYaXx1aHYJIcjmsSv4TOjlBIsBiXuKIhN0/ezWz2lddkUFV9dQbAv/VaUhgJ1suzpcz/EEDSnZgmjhNfkyySCg6gXE3jAmlw/UGydg0hkirE3GXP4qqwgc67DVhprcfUGs7FnFFmBQUjBVQ4RH4/e+D2uIcIRjyhIofsCfj1FODpNqqygwGV73iTjy91DnB0fsGlJp67NCPtgGjWbrNKMVZrWSVsL6xiLkcELbfZRfw2RKZtBkZ/CQo0mKzdje+kIHwYTTk1m0vKx3IHNLp9eE14lI8tslvMiyyYrh1RzSAk1IBsKzoTbsvdodtcMy6arwbQGdg0XZqZMvTY3WVewmuclRIgbPlJC4dtCeTbtttHvwE7Oixxq7a0M2Vs0WVewmh0B7LuTuHL3aVk/hL2VhkWTdQUrodFMNBo8UF6VDslri7SskeJ2dnG5K1XEkiearCtY6YzU3wIjpaJXTtdKFr6zPde9kalLiyZrFSs3Y4WzZEjT0IiLVa1q/phY1zZZH08dSI2EkSwWdIxGXCSSbcVq02St6pXLHgG7J7grVW5/FTS8rVhtmqxVrKU/w/BJE6hPOjTpv0LVdbAdX2fYsE2TtYq1pDc23Ui9or9B0WiVXg3DtWmyVrCWWg3Cv830ql1xmOsgOKyarBWspWbsXEFRkIw4P7pa9j8Sbbl2TdYy1lIzFicUo4zGQ8t4k/DhHYcB7a6W+cd+nr561z32fw6T52wjKIMRb3PZ5rW6SbF55eFWUF1i92OB/eXWjf+/3tQ6kboNOQk4CTgJOAk4CTgJOAk4CTgJOAk4CTgJOAk4CTgJOAk4CTgJbCWBv4H1mb3jzegBAAAAAElFTkSuQmCC) | | F1-score (F1) | The harmonic mean of precision and sensitivity, providing a single measure of performance | ![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAATsAAABNCAMAAADuDjOdAAAAAXNSR0IArs4c6QAAALpQTFRFAAAAAAAAAAA6AABmADo6ADpmADqQAGa2OgAAOgA6OgBmOjoAOjo6OjpmOjqQOmZmOmaQOma2OpC2OpDbZgAAZgA6ZgBmZjoAZjo6ZjpmZjqQZmY6ZpC2ZpDbZrbbZrb/kDoAkDo6kDpmkGYAkGY6kLbbkLb/kNv/tmYAtmY6tmZmtpA6tpBmtrZmttvbttv/tv//25A625Bm27Zm27aQ29u22//b2////7Zm/9uQ/9u2/9vb//+2///bARSCRQAAAAF0Uk5TAEDm2GYAAAAJcEhZcwAAEnQAABJ0Ad5mH3gAAAAZdEVYdFNvZnR3YXJlAE1pY3Jvc29mdCBPZmZpY2V/7TVxAAAGsklEQVR4Xu1b60LbNhSOoJB0g61AttD1QtaN1WzdqAlr4tR+/9fauelqZThR4kKQfzRG0bno09GRjvJ1MMhPRmBPEFhenu/JSKLDuDu6Thre10/fK3VwNo8pKQ/eY3Oh6DlNM5Tk5Y6E7yfxgXczt1CHt4OmUMcR8MqDj6Sk+U1dDAbLsbrqoLMeH372urUaOihZp0uov23vfzyox8N1bPl9F4RPPRaY3C8XBquCvl2oLnaqUdCr1bC5s1FJ0W9QaNsLe7hqqhGERdLTTNvYNVMdi82UQqkbdkl+JAjbmV6lJN6jCJbJ2i7E4s7OiMxbiSv3sT7LyOz7vq7osUgdlYTUbISK6gktVp3tMOAQNMb3DrssJ4pm6/5SqSPOiTPYcdTZoBrxN/z36R+2gVqws9YggtTZtetKOzYCsUA/b2fgIDvQ3OA7mDq40h7YHhNK7trTpIwH/ulFefYFVidAh8OXNhwJoXhPiJbvYCn//XZeYFK7UefzaoSOwJp/D94AxBKjC/Vu0PyFEtIAPa6x28fy7Zfp8exsvrB5ceHaRYNW2tgIxFr6vXRXv74twfbNRTMFP/10t/xljt+ZFOSMc6M1xXuBwFQc3+JbPbbpjub0FTXD+CG8uC+GI2fFAr2pxwQVLWzKIvXP8I9uEH9J5AVo8FM6TI/Y9aUdG55YS79ZemKP8anHNJ3kkl2cvMp4lwRXI9tkdxTt6gQ7JwKZxc6+ocpqJMmVI44g4FcBlA4yhaKTISGMDdKD/GcNfp6xdj1px4YvFuoXM8aeIFJ9B5bYA9tDwoKWTip2JSQF/cQQ80dp/sKQ+PppdGS9I08Y2XqsJE65QfYZ+mANMiIx7M+PkXZs+GKhfrvy2J78XaIlrwXt0Zf1T7LWUuLO27ptBNnR+KPUpmDtKvXiR6w1nBOOkWp+VzQl0iBS9EFhHqQZJ3JpdCzt2AjFAv3GrrzwR3UKOHotNtRKfQxNwI5Sgn6ql2/06rfYecrNmJ3BO69OFuOlyQ3Sg3Tye5AIrF1xhaRdgFtirn6bPcUB+mh+xdByW0Q5REP1Uqf4BOwKSlXNlLPSlYlC47df0djA4gNz9YNGg/ySFI2bBKccCVp2ULYUtOglAteubDEkLT6gDQFbjkst/dR+hwcBTm74d3Pj5jnTg8Z5PtURk7DPLrjOV2ASXofg4rD5k6ZE7yAzdainSJzi2StxO5jR9FEa+/eSXV++/ryAw8ugomXDDaD7fA5nFtyhefjw7/ID6w3sUl8jbWy4YjH9w8EMIRZ7sLcuL3m7Mh5ID2iCbGk2N38BSGB2+5A7EsSuUHipUCq+PNGRgeA6u7iT/OC0KgcXOIqqg9N/cEYnuEdUb0b6boYb6PgrFzESmlN1xIee0K4rjQdpseGKtfVj+gPEtQPOXiUuUQqlHvC4JWhyXRHD2daz3WbhCfVyc3xyPRsd98PV9ROCy3WVqg15TNVuFiAfspYmIW42SOfMvJmCxyhVXiwn9urEVrN8SQnFPOFK1XPSI/fGSToemzDEF9z2ynM/svfG+giK17yTU74PSXn27/cKOG6/0vXjwP29wj2vfrgOyp4UDJ+BbFDZ+MXUMxh/yhDprGKK3Bx362BJgVaaXcSPO1NBYBUR+1lsHUt715euHRxc8prtPsVUtTf2UNcBO13NPu9PgJjvkgtzMOmAXfeJ2fOe4Y+AOd91n/DwFi/HXWfsWrSLjF1n7PxLSrilKtTQlB+dtTzLjsEl5d0J7Z0na5a0+1fBtoIhlXW3Krr2nnGHA09j3a2Gjn9CeNyMO/QwSO4Rit2uWHcrsOuXcWeLxw1ybMCpi5D8dsu6C13umXGXhJ34/vDvAjti3YXY9cy42wJ23451F2LXI+MOTbew05w6TOfC7OvGuhMuXY+su/aStbdXO2bckekQuwTWHeW2Pll3IXb9Mu7a2FlO3dqsO51t+mLdtTa3/hh38WtYw6lbn3Wnq0+u6HfOumufCyx2fTDu2nGXwLqTH/d7Yt1FjlT9Mu4ie8XGrDvteT+su9hptFfGXWyfFTrt+qw7fTLuh3UXw65Pxl17nyVi94asO53uemHdxYug3hh3bN4/o6Sw7oRih4z23bPu4tjJFrF7xl0MO4dTtzbrTpP8LIO7d9bdHjPucLZ2y7p7uLLe4MbjsYhEWXdbdG4vGXeUW1ew7raIHRC09f/H2abWb69rJetum67t5+8VK1l324Qu68oIZAQyAhmBjEBGICOQEcgIZAQyAhmBjEBGICOQEcgIZAQyAhmBjEBGICOQEcgIPC0E/gNFqVwtVBwfigAAAABJRU5ErkJggg==) | | False Positive Rate (FPR) | The proportion of negative ground truth peaks incorrectly classified as positive | ![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAJUAAABICAMAAADI3EeMAAAAAXNSR0IArs4c6QAAAG9QTFRFAAAAAAAAAAA6AABmADpmADqQAGa2OgAAOgA6OgBmOjpmOmaQOpDbZgAAZgA6ZjoAZpC2ZpDbZrbbZrb/kDoAkGYAkLbbkNv/tmYAttv/tv//25A627Zm27aQ29u22////7Zm/9uQ/9u2//+2///bcHGEfQAAAAF0Uk5TAEDm2GYAAAAJcEhZcwAAEnQAABJ0Ad5mH3gAAAAZdEVYdFNvZnR3YXJlAE1pY3Jvc29mdCBPZmZpY2V/7TVxAAACM0lEQVRoQ+1Za0/DMAxsNmC8yngUKJTHSvv/fyOJ7bbJopmoauYJ2R+mWta009lxrrei0DghBmpDsX4vCkquX6QB9s/m1mLothf2E5Of0jxJw6pXrw5C7bAVmOyMgygZfeVaR0GJPKp2c/49omo3QFIDTZWMnUPQ3UEXbetgxkpsqmA0brIb6mLj4HxuxYe9r2AxYBcxWd28CbIEP92Vlqa+wjnqSm/GJJHhcatxP+FYnUAExw1m7BSCdihCCRJBdDBWQwSJIKjiw6yn1RQkgqh2bhEMsIJEENS/+2lGlgmKNEaWkUirJNYMI8uwFIiRI80KI8uoJIGKkWUk0iSkIyPL8DbtKwGRxsgyKH1tV4/jMLlNOEZGWcLIMhJp95PuPtKos7IMRFqC8vf5W+IZlNhBWQYlmYuekWUwVuGsH2uuGFmGpVpgszOyjEoS24qRZVQS2OyMLBtKXWnOxF/g5q2kGS5VMy0R2AKYLmptzXCp3FmC6aB7LYe1leBSDWYDdqMr7R0LgNpLNLiWt7ZSXKoQVXsFjqTdiPBk9yK8oy164lNcqhCVo8ff0DmsrRSXKkbluyE5rK0UlypG5bcrg7WV5FLFqLw7OIe19YdLdehO9+7gHNZWkksVceUf3BzWVpJLFaHyD24OayvJpYpQ+fxksLbSXKoIlafccijeNJdqH5X9O+dh0AIZrK1El2oPlVVIZnBv1dqaJ9X0W8qAMqAMKAPKgDKgDCgDyoAyoAzMZeAXFe0x0kI1be8AAAAASUVORK5CYII=) | | True positive (TP) | A detected peak that was within three pixels of a ground truth peak |   | | False positive (FP) | A detected peak that could not be matched to any ground truth peak |   | | False Negative (FN) | A ground truth peak that was not matched to any detected peak |   | ### 3.Quantitative data The folder of Quantitative data contains two .xlsx files: for the raw and processed quantitative data shown in the main article and for the data shown in the supplementary files. Quantitative data_OS phagosomes_main article.xlsx file contains numerical data on the quantified OS phagosomes counted by the semi-automatized analysis tool. The data is organized into multiple Excel worksheets, where the cone subtype (UV, Blue, Green) and the light condition (LD= normal light cycle; DD=constant darkness) are indicated. Each worksheet contains the number of outer segment (OS) phagosomes counted by the analysis tool from the RPE tissue of the whole eye cryosection and RPE length, which was manually measured from the RPE-OS interface individually in each cryosection. In addition, the number of OS phagosomes/1 µm RPE, as well as OS phagosomes/10 µm RPE were calculated and shown in the columns. These data are shown at each studied time point as zeitgeber times (ZT1, ZT3, ZT5, ZT10, ZT16, ZT18, ZT23). The h-value was fixed to get value 15, whereas the prominence values were set specifically by the user for each sample in the analysis tool during quantification of the OS phagosomes and are shown in the worksheets. The file also contains worksheets for LD and DD conditions with the averages of OS phagosomes/10 µm of RPE at each time point for UV, blue and green cones shown in the same table one below another in order UV, blue and green cones (named as Averages_UV, Blue, Green_LD and Averages_UV, Blue, Green_DD). Additionally, averaged values were normalized to the largest value of the averages throughout the time points for each cone subtype individually and shown in separate worksheets for LD and DD conditions (named as Normalized_UV,Blue,Green_LD and Normalized_UV,Blue,Green_DD). The order of photoreceptor types is UV, blue and green. Quantitative data_OS_phagosomes_supplementary material.xlsx file contains raw and processed numerical data on the quantified OS phagosomes counted by the analysis tool (rod precursors, LD) and raw and processed numerical data of the manually quantified OS phagosomes (red cones, LD and red cones, DD). Worksheets and columns are named similarly to the ones in Quantitative data_OS phagosomes_main article .xlsx file and the red cone and rod precursor data are shown in the same format. Abbreviations in the excel work sheets: OS phagosome= outer segment phagosome, ZT=Zeitgeber time, LD= normal light cycle, DD=constant darkness, µm= micrometer ## Code/software ### **Semi-automatized analysis tool for quantitative analysis of outer segment (OS) phagosomes in the RPE tissue** We developed the analysis tool for ImageJ-Fiji (ImageJ 1.54). It uses standard ImageJ-Fiji libraries for general image processing and MorphoLibJ (MorphoLibJ_-1.6.2) for all morphological operations. The analysis tool utilizes confocal z-stack images of fluorescently labeled whole eye cryosections containing RPE and photoreceptor channels as input. Before further analysis, the stacks are processed into maximum intensity projections and denoised with Gaussian filter with radius 1. The workflow is performed as an interactive script that allows the user to select the channels from the z-stack image, adjust certain parameters, and select the RPE area from which the OS phagosomes are then counted by the analysis tool. First, the tool enables the user to define which channels to use for analysis, after which the workflow of the analysis tool is divided into two steps: Manual RPE area selection by the user and OS phagosome segmentation. The tool then allows the user to define RPE area after which the selected area is binarized through auto-thresholding using Li's method. A morphological opening with a disc-shaped structuring element of radius 6 pixels is then applied automatically by the tool to remove holes and smooth the selection. The final selection is saved as an ImageJ region of interest (ROI) to be used as the area to search for OS phagosomes in the following steps. OS phagosome segmentation is based on the extraction of local bright objects using the h-dome transform method to find all local maxima. h-dome transform dampens low-contrast regions and boosts significant intensity peaks bigger than defined h-value. The Maxima Finder is applied to find the intensity peaks with an additional criterion called prominence. Prominence is the minimum height difference between a local intensity maximum and its surrounding area for the maximum to be considered significant. In the analysis tool, the prominence value of the analysis tool can be adjusted by the user to get any value between 1 and 20. The tool then quantifies the total number of detected intensity peaks representing the OS phagosomes within the segmented RPE area. The results are saved as ROI images and optionally as an Excel file.