This method describes an immunofluorescence protocol and quantification pipeline for evaluating protein distribution with varied nuclear organization patterns in human T lymphocytes. This protocol provides step-by-step guidance, starting from sample preparation and continuing through the execution of semi-automated analysis in Fiji, concluding with data handling by a Google Colab notebook.
Various nuclear processes, such as transcriptional control, occur within discrete structures known as foci that are discernable through the immunofluorescence technique. Investigating the dynamics of these foci under diverse cellular conditions via microscopy yields valuable insights into the molecular mechanisms governing cellular identity and functions. However, performing immunofluorescence assays across different cell types and assessing alterations in the assembly, diffusion, and distribution of these foci present numerous challenges. These challenges encompass complexities in sample preparation, determination of parameters for analyzing imaging data, and management of substantial data volumes. Moreover, existing imaging workflows are often tailored for proficient users, thereby limiting accessibility to a broader audience.
In this study, we introduce an optimized immunofluorescence protocol tailored for investigating nuclear proteins in different human primary T cell types that can be customized for any protein of interest and cell type. Furthermore, we present a method for unbiasedly quantifying protein staining, whether they form distinct foci or exhibit a diffuse nuclear distribution.
Our proposed method offers a comprehensive guide, from cellular staining to analysis, leveraging a semi-automated pipeline developed in Jython and executable in Fiji. Furthermore, we provide a user-friendly Python script to streamline data management, publicly accessible on a Google Colab notebook. Our approach has demonstrated efficacy in yielding highly informative immunofluorescence analyses for proteins with diverse patterns of nuclear organization across different contexts.
The organization of the eukaryotic genome is governed by multiple layers of epigenetic modifications1, coordinating several nuclear functions that can occur within specialized compartments called nuclear bodies or condensates2. Within these structures, processes such as transcription initiation3, RNA processing4,5,6, DNA repair7,8, ribosome biogenesis9,10,11, and heterochromatin regulation12,13 take place. The regulation of nuclear bodies adjusts over both spatial and temporal dimensions to accommodate cellular requirements, guided by principles of phase separation14,15. Consequently, these bodies function as transient factories where functional components assemble and disassemble, undergoing changes in size and spatial distribution. Hence, understanding the characteristics of nuclear proteins by microscopy, including their propensity to form bodies and their spatial arrangement in different cellular conditions, offers valuable insights into their functional roles. Fluorescence microscopy is a widely used method for studying nuclear proteins, allowing their detection through fluorescent antibodies or directly expressing targets with a fluorescent protein reporter16,17.
In this context, nuclear bodies appear as bright foci or puncta, with a notable degree of sphericity, making them easily distinguishable from the surrounding environment16,18. Super-resolution techniques like STORM and PALM, by providing improved resolution (up to 10 nm)19, enable more precise characterization of the structure and composition of specific condensates20. However, their accessibility is limited by equipment expenses and the specialized skills needed for data analysis. Therefore, confocal microscopy remains popular due to its favorable balance between resolution and wider usage. Such popularity is facilitated by the inherent removal of out-of-focus light, which diminishes the requirement for extensive post-processing procedures for accurate segmentation, its widespread availability in research institutes, its effective acquisition time, and sample preparation that is typically efficient. However, accurately measuring protein distribution, assembly, or diffusion using immunofluorescence assays across diverse cellular conditions poses challenges, as many existing methods lack guidance on selecting suitable parameters for proteins with varying distribution patterns21. Moreover, handling the resulting large data volume can be daunting for users with limited experience in data analysis, potentially compromising the biological significance of the results.
To address these challenges, we introduce a detailed step-by-step protocol for immunofluorescence preparation and data analysis, aiming to provide an unbiased method for quantifying protein staining with various organization patterns (Figure 1). This semi-automated pipeline is designed for users with limited expertise in computational and imaging analysis. It combines the functionalities of two established Fiji plugins: FindFoci22 and 3D suite23. By integrating the precise foci identification capability of FindFoci with the object identification and segmentation features in 3D space offered by 3D suite, our approach generates two CSV files per channel for each field of acquisition. These files contain complementary information that facilitates the calculation of metrics suitable for various types of signal distribution, such as the count of foci per cell, the distance of foci from the nuclear centroid, and the inhomogeneity coefficient (IC), which we have introduced for diffuse protein staining. In addition, we acknowledge that data extrapolation can be time-consuming for users with limited data handling skills. To streamline this process, we provide a Python script that automatically compiles all collected measurements into a single file for each experiment. Users can execute this script without the need to install any programming language software. We provide an executable code on Google Colab, a cloud-based platform that allows the writing of Python scripts directly in the browser. This ensures that our method is intuitive and readily accessible for immediate use.
We demonstrate the effectiveness of our protocol in analyzing and quantifying alterations in signal distribution of two nuclear proteins: Bromodomain-containing protein 4 (BRD4) and Suppressor of zeste-12 (SUZ12). BRD4 is a well-documented coactivator protein within the Mediator complex known to form condensates associated with polymerase II-dependent transcriptional initiation24,25. SUZ12 is a protein component of the Polycomb Repressive Complex 2 (PRC2) responsible for regulating the deposition of H3K27me3 histone modification26,27. These proteins exhibit different patterns within two distinct cell types: freshly isolated human CD4+ naïve T cells, which are quiescent and exhibit slow rates of transcriptional activity, and in vitro differentiated TH1 CD4+ cells, which are specialized, proliferating effector cells showing increased transcription28.
The use of human samples for research purposes was approved by the Ethics Committees of the Fondazione Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS) Cà Granda Ospedale Maggiore Policlinico (Milan), and informed consent was obtained from all subjects (authorization numbers: 708_2020). The protocol is organized into three primary sections: immunofluorescence execution, image acquisition, and image analysis. On average, it necessitates 4 working days to be completed (Figure 1).
1. Immunofluorescence preparation
NOTE: This immunofluorescence protocol can be easily customized for various cell types and protein targets by adjusting the fixation and permeabilization conditions. Immunofluorescence preparation typically takes less than 3 days to be completed, with the duration of primary antibody incubation varying based on antibody quality and the target protein (Figure 1).
2. Image acquisition
NOTE: The duration of image acquisition depends on the instrument and selected settings.
3. Image analysis
The outlined protocol in this method facilitates the visualization and quantification of alterations in nuclear protein staining within human primary T cells, and it can be customized for diverse cell types and protein targets. As case studies, we conducted and analyzed the staining of BRD4 and SUZ12 in naïve and TH1 CD4+ cells.
BRD4 displays a well-dotted staining pattern in both quiescent naïve and differentiated TH1 CD4+ cells, enabling the use of designed parameters for foci identification for both cell types (Figure 3A). Our pipeline allows for the quantification of various parameters, including the number, volume, fluorescence signal (expressed as mean fluorescence intensity (MFI)), coverage, and spatial arrangement of BRD4 foci within the nucleus (Figure 3B-F). Additionally, we have introduced an automatic calculation to determine the foci's percentage of the nuclear volume, facilitating comparison across cells with different nuclear sizes (Figure 3E).
Finally, we have provided the option to map the foci's distances from the nuclear center, aiding in understanding their nuclear positioning (Figure 3F). The results obtained from the quantification highlight notable alterations in BRD4 foci between TH1 and quiescent CD4+ T cells, encompassing an increase in foci number, size, brightness, and volume, as well as differences in their distribution and localization. These findings are consistent with the heightened transcriptional activity of activated/proliferating T cells compared to their quiescent counterpart28.
In contrast, SUZ12 does not participate in condensate formation. In our investigation, we observed a significant disparity in the staining pattern of SUZ12 between the two cellular states: a punctate pattern in naïve T cells transitions to a diffuse pattern in TH1 CD4+ cells (Figure 4A). This significant shift in protein behavior hampers the comparability of foci features using identical parameters since dot identification in TH1 CD4+ cells is largely unsuccessful (Figure 4A). While it is generally recommended to maintain consistency in both acquisition and analysis parameters when comparing immunofluorescence staining across distinct cellular conditions, adjustments become necessary in cases of substantial changes in the target of interest.
In this specific scenario, as illustrated in Figure 4B, despite experimenting with various parameter combinations, including adjustment to background parameters, Gaussian blur, and search parameters (as suggested in protocol step 3.4.5), FindFoci encounters difficulties in distinguishing true signal peaks from insignificant ones due to the inherent distribution nature of the protein. Moreover, we emphasize that conventional measurements like nuclear MFI can be misleading (Figure 4C), failing to consider alterations in signal distribution.
Hence, to measure such changes, we suggest the use of the coefficient of variation, a metric we have introduced in this protocol to characterize protein diffusion, referred to as the inhomogeneity coefficient (IC). The IC accurately captures the diffusion of SUZ12 foci in TH1 CD4+ cells without forcing their identification with incorrect segmentation approaches (see Figure 4D). Finally, by using DAPI as a control, which remains uniformly distributed across both cell conditions, we further validate the efficacy of this parameter (Figure 4E).
Figure 1: Schematic representation of immunofluorescence preparation and data analysis. The protocol for staining and analyzing nuclear protein immunofluorescence consists of three primary phases: sample preparation, image acquisition, and analysis, spanning approximately 4 working days. The script utilizes two plugins, 3D Suite to conduct nuclei segmentation and FindFoci to process the protein channel. Spatial and quantitative information regarding the foci and the nucleus are collected in separate csv files. A script executable on Google Colab links foci measurements with corresponding nuclei using the nucleus bounding box and computes derived measurements based on collected data. Please click here to view a larger version of this figure.
Figure 2: Representation of nuclear protein quantification pipeline settings. (A) FindFoci popup window illustrating parameters for foci identification (protocol section 3.4.3). (B) Popup window representing the macro recording the selected parameters of FindFoci. The copied string is pasted within the start-up window (protocol steps 3.4.6-3.5.3.4). (C) Panel showing the nuclei segmentation quality control step for manually removing or modifying nuclear regions of interest using the 3D manager (protocol steps 3.5.5.1 -3.5.5.3). Please click here to view a larger version of this figure.
Figure 3: Comparisons between distinct cellular conditions based on BRD4 foci analyses. (A) Representative confocal fluorescence microscopy images of immunofluorescence staining of BRD4 (grey) in human primary naïve CD4+ T cells and TH1 CD4+ cells. Nuclei are counterstained with DAPI (blue). Original magnification 63x; scale bar = 5 µm. Bottom, BRD4 foci, identified by the pipeline, are marked with white arrowheads and countered in yellow (bottom). (B) Box plot representing BRD4 foci number/nucleus in naïve and TH1 CD4+ cells (n = 2 individuals). (C) Box plot representation of volume (µm3) of BRD4 foci in naïve and TH1 CD4+ cells (n = 2 individuals). (D) Box plot representation of foci MFI of BRD4 foci in naïve and TH1 CD4+ cells (n = 2 individuals). (E) Box plot representing the percentage of the nuclear volume occupied by BRD4 foci in relation to the total nucleus volume (n = 2 individuals). (F) Representation of distance frequencies of BRD4 foci from nuclear centroid to nuclear periphery (n = 2 individuals). Please click here to view a larger version of this figure.
Figure 4: Inhomogeneity coefficient, a metric to quantify the transition from a punctate to a dispersed pattern in SUZ12 immunofluorescence signals. (A) Representative confocal fluorescence microscopy images of immunofluorescence staining for SUZ12 (magenta) in human primary naïve CD4+ T cells and TH1 CD4+ cells. Nuclei are counterstained with DAPI (blue). Original magnification 63x; scale bar = 5 µm. Bottom, SUZ12 foci identified by the pipeline are marked with a white arrowhead and countered in orange (B). Examples of wrong foci identification in TH1 CD4+ cells shown by white arrows and orange mask with two different parameter settings. (C) Box plot representation of nuclear MFI of SUZ12 in naïve and TH1 CD4+ cells (n = 2 individuals). (D) Box plot representation of SUZ12 IC in naïve and TH1 CD4+ cells (n = 2 individuals). (E) Box plot representation of DAPI IC in naïve and TH1 CD4+ cells (n = 2 individuals). Abbreviations: DAPI = 4',6-diamidino-2-phenylindole; MFI = mean fluorescence intensity; IC = inhomogeneity coefficient. Please click here to view a larger version of this figure.
Solution | Composition | Comments/Description | |
TH1 medium | RPMI with GlutaMAX-I, 10% (v/v) Fetal Bovine Serum (FBS), 1% (v/v) non-essential amino acids, 1 mM sodium pyruvate, 50 IU/mL penicillin, 50 μg/mL streptomycin, 20 IU/mL recombinant IL-2, 10 ng/mL recombinant IL-12, 2 mg/mL neutralizing anti-IL-4. | Step 1.1.4. | |
Coating solution | 0.1% poly-L-lysine in ddH2O | Step 1.2.1. | |
PBS-T | 1x PBS/0.1% TWEEN 20 pH 7.0 | Step 1.3. | |
PFA solution for fixation | 3% paraformaldehyde (PFA) diluted in 0.1% PBS-TWEEN | Step 1.2.4. | |
TPBS | 0.05% Triton X-100 diluted in 1x PBS | Step 1.2., 1.3. | |
Permeabilization solution | 0.5% TPBS diluted in 1x PBS | Step 1.2., 1.3. | |
Storage solution | 20% glycerol/1x PBS | Step 1.2.7., 1.3.1. | |
Antibody dilution buffer | 0.1% PBS-TWEEN/2% Goat Serum/1% BSA | Step 1.3.6., 1.3.8. |
Table 1: Composition of media and buffers used in this protocol.
Supplemental File 1: "convert_to_TIFF.py". This file contains the script for converting any image file (e.g., ND2, LIF, etc.) into TIFF files required for accurate quantification processes. Please click here to download this File.
Supplemental File 2: "nuclear_prot_q.py". This file contains the script that allows to measure and quantify images containing at least two channels (nucleous staining, nuclear protein staining). Please click here to download this File.
Supplemental File 3: "final_nuclear_protein metrics.ipynb". This file contains a Jupyter notebook that extracts and compiles a summary of all relevant parameters into a single Excel file, using the output of "nuclear_prot_q.py" as input. Please click here to download this File.
In this study, we present a method for performing immunofluorescence experiments on nuclear proteins in human T lymphocytes. This method offers flexibility for use with various cell types through minor modifications in fixation and permeabilization steps, as described previously30,31.
Our imaging workflow builds upon established techniques outlined in the literature, specifically FindFoci and 3D Suite22,23. Unlike previous publications, we provide an executable Jython code integrated into Fiji for a semi-automated pipeline, allowing meticulous parameter tracking and visual quality checks. Notably, we introduced the optional ability to regulate the quality of nuclei segmentation in protocol section 3.5.3.5 and implemented a corrective step in protocol section 3.5.5. With the script executable on Google Colab, users can retrieve comprehensive measurements and choose the most suitable ones for their specific biological scenarios.
The provided analysis pipeline is written in Jython (Java + Python) and requires the installation of Fiji for execution, along with two distinct plugins: 3D Suite and FindFoci (Figure 1). The pipeline begins with image conversion into TIFF format and preprocessing of the channel containing the nuclei with a Gaussian blur filter to enhance nuclei segmentation. We direct 3D Suite to capture both i) spatial information, such as nuclear volume, centroid and bounding box coordinates (representing the smallest box enclosing the object and used to associate foci to the right nucleus), and ii) quantitative information, including MFI and standard deviation of the signal. Subsequently, the pipeline employs FindFoci to identify local maxima in the foci signal, eliminating the dependence on global thresholding methods (which often rely on user discretion and may lack precision), and records both spatial and quantitative information relative to foci. The output folder will thus contain two CSV files per channel, encompassing spatial (“M_.csv”) and quantitative (“Q_.csv”) metrics.
The final file, containing all collected and derived measurements, is generated by executing a Python script on Google Colab. This file includes metrics such as the number of foci per nucleus, foci MFI, inhomogeneity coefficient (IC), nuclear MFI, foci distance from the nuclear centroid, and percentage of occupied nuclear volume. We estimate a processing time of approximately 1 minute per image including visual quality control steps. This processing time can be further expedited through the optimization of immunofluorescence settings and analysis parameters.
The pipeline computes parameters for nuclear proteins exhibiting both well-dotted and diffused patterns, allowing users to select the most appropriate parameter. We recommend utilizing the inhomogeneity coefficient to describe protein foci dissolution.
A critical step to ensure proper foci identification involves the setting of parameters to discriminate genuine signal peaks from background noise. It is crucial to conduct thorough testing on random images before running the pipeline. While our study primarily focuses on a single channel, the pipeline supports the introduction of various dye combinations targeting different proteins in immunofluorescence experiments. It does not support simultaneous analysis of multiple channels and must be run separately for each channel. Furthermore, we have opted to present specific measurements in this protocol; however, additional metrics, such as integrated intensity for assessing absolute changes in fluorescence signal or physical parameters like circularity and compactness, can also be incorporated into the pipeline in section 3.3.2. The inclusion of metrics not currently integrated into this pipeline necessitates proficiency in Python programming to integrate them in the final metrics file; otherwise, they will be automatically included in the spatial (“M_.csv”) and quantitative (“Q_.csv”) information files.
The timing of this pipeline depends on factors such as image size and quantity, number of channels per analysis, and microscope settings. Additionally, machine performance influences execution time, with visual steps potentially slowing down the process.
It is worth noting that while confocal microscopy offers valuable insights, it remains a diffraction-limited technology that may limit the precision of observed nuclear structures or protein foci. To enhance the robustness of the analysis, we suggest increasing the number of nuclei and biological samples examined. As the number of nuclei analyzed rises, any notable changes in nuclear morphology or protein distribution become more evident, irrespective of resolution limitations.
In terms of potential future applications, this method exhibits versatility and adaptability to various cell types and staining protocols. For example, it can be readily employed in radiobiology experiments, where accurate quantification of ionizing radiation-induced foci (IRIFs) serves as a direct indicator of DNA damage32. While our primary emphasis lies on nuclear analysis, the protocol can be easily customized by incorporating a cytoplasmic dye and substituting it for the nuclear channel within the script. Furthermore, the analysis can seamlessly integrate into other experimental methodologies such as DNA and RNA FISH30,31, leveraging default parameters for nuclear topology analysis.
The authors have nothing to disclose.
We acknowledge the scientific and technical assistance of the INGM Imaging Facility, in particular, C. Cordiglieri and A. Fasciani, and the INGM FACS sorting facility in particular M.C Crosti (Istituto Nazionale di Genetica Molecolare 'Romeo ed Enrica Invernizzi' (INGM), Milan, Italy). We acknowledge M. Giannaccari for his technical informatic support. This work was funded by the following grants: Fondazione Cariplo (Bando Giovani, grant nr 2018-0321) and Fondazione AIRC (grant nr MFAG 29165) to F.M. Ricerca Finalizzata, (grant nr GR-2018-12365280), Fondazione AIRC (grant nr 2022 27066), Fondazione Cariplo (grant nr 2019-3416), Fondazione Regionale per la Ricerca Biomedica (FRRB CP2_12/2018,) Piano Nazionale di Ripresa e Resilienza (PNRR) (grant nr G43C22002620007) and Progetti di Rilevante Interesse Nazionale (PRIN) (grant nr 2022PKF9S) to B.B.
1.5 mL Safe-Lock Tubes | Eppendord | #0030121503 | Protocol section 1 |
10 mL Serological pipettes | VWR | #612-3700 | Protocol section 1 |
20 µL barrier pipette tip | Thermo Scientific | #2149P-HR | Protocol section 1 |
50 mL Polypropylene Conical Tube | Falcon | #352070 | Protocol section 1 |
200 µL barrier pipette tip | Thermo Scientific | #2069-HR | Protocol section 1 |
antifade solution – ProlongGlass – mountingmedia | Invitrogen | #P36984 | Step 1.3.12 |
BSA (Bovine Serum Albumin) | Sigma | #A7030 | Step 1.3.6., 1.3.8. |
CD4+ T Cell Isolation Kit | Miltenyi Biotec | #130-096-533 | Step 1.1.2. |
DAPI (4,6-diamidino-2-phenylindole) | Invitrogen | Cat#D1306 | Step 1.3.10. |
Dry ice | Step 1.3.1. | ||
Dynabeads Human T-activator anti-CD3/anti-CD28 bead | Life Technologies | #1131D | magnetic beads step 1.1.4. |
EtOH | Carlo Erba | #4146320 | Step 1.2.1.1. |
FACSAria SORP | BD Bioscences | Step 1.1.3. Equipped with BD FACSDiva Software version 8.0.3 | |
FBS (Fetal Bovine Serum) | Life Technologies | #10270106 | Step 1.1.4 |
FICOLL PAQUE PLUS | Euroclone | GEH17144003F32 | Step 1.1.1. |
FIJI Version 2.14.0 | – | – | Protocol section 3 |
Glass coverslip (10 mm, thickness 1.5 H) | Electron Microscopy Sciences | #72298-13 | Step 1.2.1. |
Glycerol | Sigma | #G5516 | Step 1.2.7-1.3.1. |
Goat anti-Rabbit AF568 secondary antibody | Invitrogen | A11036 | Step 1.3.8. |
HCl | Sigma | #320331 | Step 1.3.4. |
human neutralizing anti-IL-4 | Miltenyi Biotec | Cat#130-095-753 | Step 1.1.4. |
human recombinant IL-12 | Miltenyi Biotec | Cat#130-096-704 | Step 1.1.4. |
human recombinant IL-2 | Miltenyi Biotec | Cat#130-097-744 | Step 1.1.4. |
Leica TCS SP5 Confocal microscope | Leica Microsystems | – | Protocol section 2, Equipped with HCX PL APO 63x, 1.40 NA oil immersion objective, with an additional 3x zoom. Pinhole size : 0.8 AU. Line average 2×. Frame size 1024×1024 pixel. |
MEM Non-Essential Amino Acids Solution | Life Technologies | #11140035 | Step 1.1.4. |
Microscope Slides | VWR | #631-1552 | Step 1.3.12. |
Mouse monoclonal anti-Human CD4 APC-Cy7 (RPA-T4 clone) | BD Bioscience | #557871 | Step 1.1.3. |
Mouse monoclonal anti-Human CD45RA PECy5 (5H9 clone) | BD Bioscience | #552888 | Step 1.1.3. |
Mouse monoclonal anti-Human CD45RO APC (UCHL1 clone) | Miltenyi Biotec | #130-113-546 | Step 1.1.3. |
Multiwell 24 well | Falcon | #353047 | Protocol section 1 |
Normal Goat Serum | Invitrogen | PCN5000 | Step 1.3.6., 1.3.8. |
PBS | Life Technologies | #14190094 | Protocol section 1 |
Penicillin/Streptomycin solution | Life Technologies | #15070063 | Step 1.1.4. |
PFA | Sigma | #P6148 | Step 1.2.4. |
poly-L-lysine | Sigma | #P8920 | 1.2.1. |
Primary antibody – BRD4 | Abcam | #ab128874 | Step 1.3.6. |
Primary antibody – SUZ12 | Cel Signalling | mAb #3737 | Step 1.3.6. |
RPMI 1640 W/GLUTAMAX-I | Life Technologies | #61870010 | Step 1.1.4. |
Sodium Pyruvate | Life Technologies | #11360039 | Step 1.1.4. |
Triton X-100 | Sigma | #T8787 | Step 1.2., 1.3. |
TWEEN 20 | Sigma | #P9416 | Step 1.3. |
Tweezers | – | – | Protocol section 1 |
.