Jul 31, 2023

Public workspaceImage processing to investigate NEMO recruitment and involvement in mitophagy and inflammatory signaling

  • 1Department of Physiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104;
  • 2Aligning Science Across Parkinson’s (ASAP) Collaborative Research Network, Chevy Chase, MD, 20815
Icon indicating open access to content
QR code linking to this content
Protocol CitationOLIVIA HARDING, Erika L.F. Holzbaur 2023. Image processing to investigate NEMO recruitment and involvement in mitophagy and inflammatory signaling. protocols.io https://dx.doi.org/10.17504/protocols.io.n2bvj61xxlk5/v1
License: This is an open access protocol distributed under the terms of the Creative Commons Attribution License,  which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Protocol status: Working
We use this protocol and it's working
Created: June 23, 2022
Last Modified: May 31, 2024
Protocol Integer ID: 65163
Keywords: Confocal imaging, Deconvolution, Image processing, Image analysis, Fluorescent intensity analysis, Subcellular imaging, Organelle morphology, Segmentation, ASAPCRN
Funders Acknowledgement:
Aligning Science Across Parkinson’s
Grant ID: Mechanisms of mitochondrial damage control by PINK1 and Parkin (ASAP-000350)
Abstract
Beautiful images are not sufficient to robustly characterize and interrogate a cellular mechanism. In order to show recruitment of NEMO and its relationship to OPTN and other mitophagy factors, we processed hundreds of still and timelapse images and extracted quantifiable data in order to perform statistical analysis comparing different conditions. In many cases, we used software to deconvolve confocal fluorescent images, allowing us to algorithmically surpass the resolution limit. This was especially useful for live cell images that had been collected with low power settings to preserve cell health. We also employed machine learning software to generate binary segmentations and carry out particle analysis on putatively overlapping structures. Finally, in some cases we simply identified fluorescent structures by hand and measured their fluorescent intensities. We approached image analysis in a multitude of creative, effective ways, and importantly we maintained consistency of analysis within experiments in order to present the results with integrity and reproducibility.
Attachments
Guidelines
  • This protocol addresses multiple image analysis pipelines used in the corresponding manuscript to investigate NEMO and its role in mitophagy.
  • We developed several methods for quantitative analysis with our varying experiments.
Materials

Equipment/software

  • ImageJ/FIJI
  • Ilastik
  • Deconvolution software such as Huygens
  • Excel
Before start

  • For each experiment, collect images for at least 10 cells from at least three biological replicates.
  • Maintain consistent imaging parameters by saving and/or recording laser powers and exposure times.
Blinded whole-cell recruitment assessment
Blinded whole-cell recruitment assessment

Note
This method was used to assess NEMO recruitment in cells before and after AntA/OligA treatment, in the presence or absence of Parkin or p62.

Start with Z stack images of live cells collected with a 60X objective.
Max project the Z stack fields of view by ~2 um to capture the center to upper half of the cell.
Crop each cell and save only the NEMO channel with a file name to identify the conditions.
Copy cropped NEMO channel images to a new folder.
Use a random number generator such as Random.org/strings/ to generate a number of strings equivalent to the total number of cells.
Copy the strings to a document and save.
Rename the copied images in the new folder as the random strings.
Arrange files in the folder by file name in order to randomly mix images.

Open newly named files and judge whether there is NEMO recruitment or not, recording YES or NO for each image by its random string file name.
After categorizing every image, unblind the results and determine the number of cells marked YES for each condition compared to the total for that condition.

Mitochondrial recruitment assay
Mitochondrial recruitment assay

Note
This method was used to assess percentages of mitochondria that had recruited NEMO, OPTN, GABARAPs, and/or LC3B in various conditions.

Start with Z stack images of cells collected with 60X objective.
Crop and save each cell visible in the field of view.
Deconvolve images with Huygens or similar deconvolution software.
Max project each cell ~2 um.
Split channels and save each channel.
Import 5-7 channels displaying mitochondria into the Ilastik software for segmentation.
Train algorithm to recognize “mitochondria” or “not mitochondria”.
Note
If there are cells with poor expression or labeling of the mitochondria, these may be discarded.

Generate binary images for all mitochondria channels.
Import segmented Ilastik results to FIJI.
Use the original image file to draw an ROI outlining the cell.
Clear Outside the cell ROI on the segmented image.
Run the threshold function with threshold set to (255,255).
Save binary image.
Repeat previous 3 steps for NEMO puncta and OPTN rings/puncta.
GABARAPs and LC3B antibodies do not produce high enough signal-to-noise to segment with Ilastik. For these experiments, only segment NEMO and mitochondria. In FIJI, draw ROI’s around GABARAPs- or LC3B-positive mitochondria by hand. Save ROIs.
To determine the % of mitochondria that recruited NEMO,
Open the binary mitochondria image in FIJI.
Use the Analyze Particles function to generate a mask of binary image, filtering particles to exclude those fewer than 5 pixels.
Summarize and add to manager.
Record the total number of mitochondria in Excel.
Open binary image of NEMO channel.
Use the Analyze Particles function to generate a mask of binary image, filtering particles to exclude those fewer than 5 pixels, and add particles to manager.
Project NEMO particles to Mito binary channel and multimeasure.
Copy results to Excel.
The count of particles with an average intensity of >127.5 is the number of mitochondria positive for NEMO.
Divide this number by the total number of mitochondria to calculate the proportion.
To determine the % of mitochondria that recruited NEMO and/or OPTN.

Repeat the first four steps above to record total number of mitochondria.
Load the binary OPTN channel to FIJI and Fill Holes.
Project mitochondrial particle ROIs onto OPTN binary image and multi-measure.
Record results in Excel.
The count of particles with an average intensity of >127.5 is the number of mitochondria positive for OPTN.
Delete mitochondria particles from the manager and add NEMO particles.
Project NEMO particles onto mitochondria and perform the same calculation.
Use the image calculator to add the mitochondria channel to the OPTN channel with Filled Holes.
Use Math > Subtract to subtract 255 and float the resulting image.
Convert to an 8-bit image.
Project NEMO particles to this new image and multi measure.

Record results in Excel.
The count of particles with an average intensity greater or equal to 64 is the count of mitochondria that are positive for both NEMO and OPTN.
To determine the % of mitochondria that recruited NEMO and/or GABARAPs or LC3B.
Use the Image Calculator to add NEMO particles to Mitochondria particles.
Use Math > Subtract to subtract 255 and float the resulting image.
Convert to an 8-bit image.
Project hand-drawn ROIs to resulting image and multi measure.
Record results in Excel.
The count of particles with an average intensity greater or equal to 64 is the count of mitochondria that are positive for both NEMO and GABARAPs or LC3B.
Use the R code library(eulerr) to generate Euler diagrams of the resulting data.
Recruitment assay
Recruitment assay
Note
This method was used to quantify the extent of NEMO, OPTN, and p62 recruitment over time.

Start with confocal timelapse images in which the timepoint with added AntA/OligA is known.
Choose several events per cell in which the mitochondria stays in the field of view for a majority of the timecourse and the protein of interest is recruited.

For each time point, use the mitochondria channel to draw a generous ROI around the mitochondria.
This will capture rings or puncta that are recruited.
Measure the intensity of NEMO and/or OPTN or p62 for every ROI and record in Excel.

Calculate the 5-frame moving average across the timecourse for each channel.
Calculate the background intensity by averaging the moving average intensity of the first ten frames of the timecourse.
Subtract the background from every frame.
Determine the maximum intensity for the event.
Normalize the timecourse intensities to the max intensity so that the intensity of each frame is a percentage of the max intensity.
Calculate the half-max timepoint based on when the normalized intensity surpasses 50%.