We thank Pedro Fernandez-Funez for the use of the microscope and camera used in this protocol. We would also like to thank Ava Schapman for providing feedback on the clarity of the protocol. Financial support was provided by The Wallin Neuroscience Discovery Fund to Nam Chul Kim.

1. Preparing for quantification
<ol>
	<li>Download and install ilastik11, ImageJ, and the ImageJ plugin Flynotyper1. See Table of Materials for website links to installation. Download the appropriate packages according to the computer operating system (Mac, Windows, Linux, etc). Follow the installation instructions exactly as stated. 
	NOTE: Following the directions exactly as stated in the provided links is imperative. The computer used needs a 64-bit operating system; otherwise, there are no other required specifications. See the Table of Materials for the computer specifications of the system used for this protocol.</li>
	<li>Use a standard image to train the machine learning model. Take an image of a healthy fly eye using a light microscope with a ring light (such as w1118 x GMR-GAL4), using the center of the eye as the focal point. Analyze males and females separately; therefore, acquire a standard image for males and for females. 
	NOTE: Settings will vary widely depending on the microscope and camera. See the second paragraph in the discussion for more general information on fly preparation and image acquisition.</li>
	<li>Take images of the experimental fly eyes using the same settings used for the standard image acquisition. 
	NOTE:Good orientation is paramount for accurate quantification. See Figure 4 as an example for proper orientation.</li>
</ol>
2. Using ilastik to detect ommatidia from fly eye images
<ol>
	<li>Before starting analysis, make sure all experimental groups, regardless of experimental condition, are in the same file folder to make the following processes easier. Make sure files are named appropriately to identify experimental groups and distinguish between males and females. 
	NOTE:&#160;We recommend printing out the protocol and using the computer to zoom in on figures, as this protocol is most easily followed visually.</li>
	<li>Open the software. 
	NOTE:&#160;It is advised not to have other applications open while running the software; otherwise, the computer might run very slowly, particularly in older or less robust systems.</li>
	<li>Click on Create New Project | Other Workflows, choose Cell Density Counting, and select a location to save the project file (Figure 1A).</li>
	<li>Click on 1. Input Data | Add New | Add Separate Image(s). Choose the standard image (Figure 1B).</li>
	<li>Click on 2. Feature Selection and select features to use. Select the features shown in Figure 1C; choose more sigma values to increase sensitivity (i.e., 10, 15, 20, etc). 
	NOTE: If too many boxes are checked, the program might run for a long time.</li>
	<li>Click on 3. Counting and set the Foreground sigma value to 5.00. Mark 50 or more individual ommatidium using the Foreground tool on the standard image. 
	NOTE:&#160;Fifty markings are typically adequate for analysis; however, the more ommatidia marked, the more accurate the results will be (Figure 1D).</li>
	<li>Click on 3. Counting and use the Background tool to mark the outside of the ommatidia (Figure 1E). Draw green lines in a lattice shape around the ommatidia. 
	NOTE: Similar to more markings on the ommatidia, the more green lines drawn, the more accurate the results will be.</li>
	<li>Click on 4. Density Export | Choose Export Image Settings to set image settings (Figure 1F). Click on Output File Info | select Format and tif for further analysis in Flynotyper.</li>
	<li>Click on 5. Batch Processing | Select Raw Data Files and select all experimental fly photos to be analyzed (Figure 1G). The list of imported images to be analyzed is seen under Select Raw Data Files (Figure 1H). Click Process all Files at the bottom of the 5. Batch Processing section. 
	NOTE: As a reminder, male and female analyses should be done separately. Depending on the computer used, this step can take 10 min or longer, especially if there are many imported images.</li>
	<li>After the software has finished, check the same file folder containing the experimental fly eye images; there will be black images named &#34;YourSampleName_Probabilities&#34; (Figure 1I).</li>
</ol>
3. Using ImageJ to prepare photos for Flynotyper
<ol>
	<li>Open the ilastik-generated .tif file (the black boxes) in ImageJ. 
	NOTE:&#160;Each ilastik generated .tif file needs to be handled separately.</li>
	<li>Use the Rectangle Select tool to make a rectangle around just the eye. After the rectangle has been drawn, crop the area by either choosing Ctrl + X or Ctrl + C (Figure 2A). Wait for a black box in the shape of the rectangle outline to appear. 
	NOTE: The computer operating system will determine whether to use &#39;Ctrl + X&#39; or &#39;Ctrl + C&#39;. Usually, &#39;Ctrl + X&#39; is for Windows and &#39;Ctrl + C&#39; is used for Mac.</li>
	<li>Open the now cut fly eye using ImageJ; click on File | New | Internal Clipboard (Figure 2B).</li>
	<li>Save the cut fly eye as a JPEG by clicking on File | Save As | Jpeg. The cut images are now in the same file folder as the original experimental images. 
	NOTE: To make the next step easier, it is recommended to create a new file folder called &#39;cut&#39; and put all of the cut images into that folder.</li>
</ol>
4. Using Flynotyper to calculate phenotypic scores
<ol>
	<li>Open Flynotyper as an ImageJ plugin by clicking on Plugins | Flynotyper.</li>
	<li>Select Add Genotypes and open the folder containing the cut fly eye images (cut file folder). The folder name will appear under Genotypes (Figure 2C).</li>
	<li>Check the Light microscope and Vertically boxes but adjust accordingly if needed. Check Stability and Distance to the Center boxes under Rank Ommatidia by: (Figure 2C).</li>
	<li>For Number of ranked ommatidia considered, input 200. 
	NOTE: In general,&#160;200 is a good number but adjust accordingly to preferences.</li>
	<li>Click Run and wait for the results of the analysis to appear (Figure 2D). 
	NOTE: This step may take 5 min or longer (or shorter) depending on the processing power of the computer used.</li>
	<li>Copy and paste the Sample File and P-Score into statistical software for further analysis.</li>
</ol>

Unbiased Quantification of Ommatidia in Drosophila Eye Models

<ol>
	<li>Iyer, J., 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=Quantitative+assessment+of+eye+phenotypes+for+functional+genetic+studies+using+Drosophila+melanogaster.">Quantitative assessment of eye phenotypes for functional genetic studies using Drosophila melanogaster.</a> G3. 6 (5), 1427-1437 (2016).</li><li>Thomas, B. J., Wassarman, D. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+fly's+eye+view+of+biology.">A fly's eye view of biology.</a> Trends Genet. 15 (5), 184-190 (1999).</li><li>Roignant, J. -. Y., Treisman, J. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pattern+formation+in+the+Drosophila+eye+disc.">Pattern formation in the Drosophila eye disc.</a> Int J Dev Biol. 53 (5-6), 795-804 (2009).</li><li>Diez-Hermano, S., Ganfornina, M. D., Vegas, E., Sanchez, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Machine+learning+representation+of+loss+of+eye+regularity+in+a+Drosophila+neurodegenerative+model.">Machine learning representation of loss of eye regularity in a Drosophila neurodegenerative model.</a> Front Neurosci. 14 (1), 1-14 (2020).</li><li>Appocher, C., Klima, R., Feiguin, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Functional+screening+in+Drosophila+reveals+the+conserved+role+of+REEP1+in+promoting+stress+resistance+and+preventing+the+formation+of+Tau+aggregates.">Functional screening in Drosophila reveals the conserved role of REEP1 in promoting stress resistance and preventing the formation of Tau aggregates.</a> Hum Mol Genet. 23 (25), 6762-6772 (2014).</li><li>Pandey, U. B., 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=HDAC6+rescues+neurodegeneration+and+provides+an+essential+link+between+autophagy+and+the+UPS.">HDAC6 rescues neurodegeneration and provides an essential link between autophagy and the UPS.</a> Nature. 447 (7146), 860-864 (2007).</li><li>Outa, A. A., 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=Validation+of+a+Drosophila+model+of+wild-type+and+T315I+mutated+BCR-ABL1+in+chronic+myeloid+leukemia:+an+effective+platform+for+treatment+screening.">Validation of a Drosophila model of wild-type and T315I mutated BCR-ABL1 in chronic myeloid leukemia: an effective platform for treatment screening.</a> Haematologica. 105 (2), 387-397 (2019).</li><li>Shirinian, 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=A+transgenic+Drosophila+melanogaster+model+to+study+human+T-lymphotropic+virus+oncoprotein+Tax-1-driven+transformation+in+vivo.">A transgenic Drosophila melanogaster model to study human T-lymphotropic virus oncoprotein Tax-1-driven transformation in vivo.</a> J Virol. 89 (15), 8092-8095 (2015).</li><li>Diez-Hermano, S., Valero, J., Rueda, C., Ganfornina, M. D., Sanchez, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=An+automated+image+analysis+method+to+measure+regularity+in+biological+patterns:+a+case+study+in+a+Drosophila+neurodegenerative+model.">An automated image analysis method to measure regularity in biological patterns: a case study in a Drosophila neurodegenerative model.</a> Mol Neurodegener. 10 (1), 1-10 (2015).</li><li>Yusuff, 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=Drosophila+models+of+pathogenic+copy-number+variant+genes+show+global+and+non-neuronal+defects+during+development.">Drosophila models of pathogenic copy-number variant genes show global and non-neuronal defects during development.</a> PLoS Genet. 16 (6), e1008792-e1008792 (2020).</li><li>Berg, S., 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=ilastik:+interactive+machine+learning+for+(bio)image+analysis.">ilastik: interactive machine learning for (bio)image analysis.</a> Nat Methods. 16, 1226-1232 (2019).</li><li>Chalmers, M. R., Kim, J., Kim, N. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Eip74EF+is+a+dominant+modifier+for+ALS-FTD-linked+VCPR152H+phenotypes+in+the+Drosophila+eye+model.">Eip74EF is a dominant modifier for ALS-FTD-linked VCPR152H phenotypes in the Drosophila eye model.</a> BMC Res Notes. 16 (1), 1-5 (2023).</li><li>Baek, 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=TDP-43+and+PINK1+mediate+CHCHD10S59L+mutation-induced+defects+in+Drosophila+and+in+vitro.">TDP-43 and PINK1 mediate CHCHD10S59L mutation-induced defects in Drosophila and in vitro.</a> Nat Commun. 12 (1), 1-20 (2021).</li></ol>

The authors have no conflicts of interest to disclose.

The ommatidia of Drosophila comprise a useful system for studying various biological functions and genetic diseases. The regularity of ommatidia is a good measurement to examine the effect of genetic mutations4. Even though several methods for calculating ommatidial regularity exist, such as manual ranking, these methods can be heavily biased. To overcome this biased approach, semi-automatic tools have been developed1,11.
<p c...

The Drosophila melanogaster genome contains ~75% of human disease-related gene orthologs. Additionally, during Drosophila eye development, approximately two-thirds of the genes in the genome are expressed, making the Drosophila eye an outstanding genetic system to investigate various molecular and cellular functions, development, and disease models1,2. Thus, the Drosophila eye system is a useful experimental tool to study various biological processes.
The Drosophila compound eye is a well-structured and comprehensive array of ~800 ommatidia that exhibit a symmetrical and hexagonal pattern3. The regularity of this hexagonal pattern can be used to estimate the effect of introducing mutations and gene expression changes in eye morphology4. Previous studies that require evaluation of eye morphology have heavily relied on manual ranking of the severity of eye phenotypes detected by the naked eye. To rank the eye phenotypes, external eye morphology images are taken by a stereomicroscope5,6. The eye phenotype of each group is assessed by splitting the external eye into four areas and calculating the proportion of degeneration in each area5,6. Then, the values are used to calculate averages that are compared with the values obtained from control flies7. The scoring is based on the extent of fusion, loss of ommatidia, and bristle organization7,8. Fly eye photos taken with a stereomicroscope are acquired by one researcher, and the eye phenotype analysis is performed by another researcher with triple validation sets7,8.
When it comes to ranking weak alterations in eye morphology with the naked eye, there are limitations4. To overcome these limitations, computational approaches such as FLEYE and Flynotyper have been developed1,9. Flynotyper is a novel computational method for quantitatively estimating morphological changes in the Drosophila eye system1. It automatically detects the Drosophila eye and individual ommatidium, calculating Phenotypic Scores (P-Scores) based on the irregularity of theeye1. A higher P-Score indicates the fly eye is more degenerated. This software was successfully used in quantifying the abnormality of Drosophila eyes10. Although Flynotyper ensures an automated process, it still cannot be successfully applied to some eye images taken by various light microscopy methods.
Qualitatively, we prefer a ring light source compared to a single-point light source, as it offers a more accurate representation of each ommatidium. However, when the ring light is used, it generates a ring-shaped shadow at the top of each ommatidium due to the hemispherical shape of the ommatidium. This ring-shaped shadow inhibits accurate ommatidial detection by Flynotyper, leading to incorrect calculation of P-Scores.
To overcome these issues, we implemented ilastik, a machine learning-based tool for various analyses, to classify ommatidia in fly eye images11. We then fed the ilastik-generated images into Flynotyper to calculate P-Scores. This allows us to quantify the Drosophila eye morphological defects unbiasedly1.

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Computer specifications</td>
			<td></td>
			<td></td>
			<td>Ryzen 5, 16 GB RAM, Nvidia RTX 3070 Super, Windows 10</td>
		</tr>
		<tr>
			<td>Flynotyper</td>
			<td>Iyer, J. et al. (2016)</td>
			<td>Download software here: https://flynotyper.sourceforge.net/imageJ.html</td>
			<td>Open source software. Do not use Flynotyper 2.0. At the time of publication, 2.0 was fairly new and this protocol is optimized for the original version of Flynotyper.</td>
		</tr>
		<tr>
			<td>ilastik</td>
			<td>Berg, S. et al. (2019)</td>
			<td>Download software here: https://www.ilastik.org/download.html</td>
			<td>Open source software. Download Version 1.4.0.post1 under Regular Builds corresponding to your computer operating system.</td>
		</tr>
		<tr>
			<td>ImageJ</td>
			<td></td>
			<td>Download software here: https://imagej.net/ij/download.html</td>
			<td>Open source software. Versions 1.53 and 1.54 were used. 1.54 is the updated version and is the default download.</td>
		</tr>
		<tr>
			<td>Leica Application Suite (LAS X)</td>
			<td>Leica Microsystems</td>
			<td>LASX Office 1.4.6 28433</td>
			<td>System and software used for z-stack acquisition.</td>
		</tr>
		<tr>
			<td>Leica Z16 APO microscope with a DMC2900 camera</td>
			<td>Leica Microsystems</td>
			<td>10 447 173, 12 730 466</td>
			<td>Referred to as Z-stack microscope and camera in the text. This product is now archived.</td>
		</tr>
	</tbody>
</table>

quantifying drosophila melanogaster eye phenotypes: a computational approach integrating ilastik and flynotyper

The Drosophila melanogaster compound eye is a well-structured and comprehensive array of around 800 ommatidia, exhibiting a symmetrical and hexagonal pattern. This regularity and ease of observation make the Drosophila eye system a powerful tool to model various human neurodegenerative diseases. However, ways of quantifying abnormal phenotypes, such as manual ranking of eye severity scores, have limitations, especially when ranking weak alterations in eye morphology. To overcome these limitations, computational approaches have been developed such as Flynotyper. The use of a ring light allows for better qualitative images accessing the intactness of individual ommatidia. However, these images cannot be analyzed by Flynotyper directly due to shadows on ommatidia introduced by the ring light. Here, we describe an unbiased way to quantify rough eye phenotypes observed in Drosophila disease models by combining two software, ilastik and Flynotyper. By preprocessing the images with ilastik, successful quantification of the rough eye phenotype can be achieved with Flynotyper.

Author Spotlight: Quantifying Rough Eye Phenotypes in Drosophila Models of Amyotrophic Lateral Sclerosis with Frontotemporal Dementia

Quantifying Drosophila melanogaster Eye Phenotypes: A Computational Approach Integrating ilastik and Flynotyper

In a previous study, we used this protocol to determine genetic modifiers of mutant VCP protein linked to amyotrophic lateral sclerosis with frontotemporal dementia (ALS-FTD)12. Additionally, this method was also used in another paper to assess the toxicity of CHCHD10S59L-mediated ALS-FTD, even when using an older stereomicroscope13. To further validate these results, we used ilastik and Flynotyper in the same ALS-FTD disease mutation as Baek et al.<sup class="xr...

The Drosophila eye system is a useful tool for studying various biological processes, particularly human neurodegenerative diseases. However, manual quantification of rough eye phenotypes can be biased and unreliable. Here we describe a method by which ilastik and Flynotyper are used to quantify eye phenotype in an unbiased way.

The Drosophila melanogaster compound eye is a well-structured and comprehensive array of around 800 ommatidia, exhibiting a symmetrical and hexagonal ...

quantifying-drosophila-melanogaster-eye-phenotypes-computational

Without Section

Research

JoVE Journal

Developmental Biology

Combining ilastik and Flynotyper for Accurate Detection and Quantification of Ommatidia in Drosophila melanogaster Fly Eye Images

472 Views.  University of Minnesota.  The Drosophila eye system is a useful tool for studying various biological processes, particularly human neurodegenerative diseases. However, manual quantification of rough eye phenotypes can be biased and unreliable. Here we describe a method by which ilastik and Flynotyper are used to quantify eye phenotype in an unbiased way. 

Unbiased Quantification of Ommatidia in Drosophila Eye Models

Summary