- Introduction on Fluorescence Lifetime Imaging Microscopy (FLIM)
- Specifications
- General pipeline
- Principles
- Phasor transform and autofluorescence filtering
- Note on the multi-component exponential fit
- Prepulse and afterpulse in the measured IRF
- Normalization of event counts
- How to cite this work?
- Acknowledgments
Jung et al., 2011 on the fluorescence lifetime of fluorescent proteins
Malacrida et al., 2021 on the phasor plot, especially literature [45, 46, 54, 59, 60, 61] cited in this review.
This toolkit analyzes raw FLIM data (in .iss-tdflim files acquired by VistaVision) obtained by an Alba v5 scanning confocal microscope equipped with a 20-MHz pulsed excitation light source. The microscope is hosted at the Single Molecule Analysis in Real-Time (SMART) Center at the University of Michigan, Ann Arbor. To operate the microscope, please read the manual. Photons are detected by an SPCM-AQRH-15 avalanche photodiode (APD). The time-correlated single photon counting module (which registers detected photon events to their corresponding excitation pulses) is an SPC-830.
(Note: steps 1-3 are only needed to be done once per experiment.)
- Open MATLAB. In the Command Window, call the
calculateIRF
function and pick the .fcs file containing the instrument response function (IRF) measurement data in the pop-up UI. Note that this requires another repository of mine, readHeader, to function. The normalized IRF is then automatically saved into a .mat file with the same filename as the original .fcs input file. Load it to continue later steps.
>> calculateIRF('Early'); % For the green channel, specify the early pulse.
- Visually examine the IRF curve. Optional: if you see a major prepulse before the main IRF spike, remove it manually by assigning those
IRFProb
values belonging to the prepulse to 0. The reason is explained in a section below. However, if you do this, make sure to renormalize theIRFProb
.
>> plot(25/4096:50/4096:25-25/4096, IRFProb, '.-'); % A ClockFrequency of 20 MHz and an ADCResolution of 4096 were used
...
>> % IRFProb = IRFProb / sum(IRFProb); % Normalization again
- Calculate
IRFTransform
. Do not use other values for$\omega$ because the empirical standards (hard-coded in the third step of the pipeline) to identify autofluorescent pixels depend on$\omega$ . For more details, see another section below. SaveIRFProb
,Omega
, andIRFTransform
into a single .mat file.
>> Omega = calculateBestOmega(2, 3); % ~ 0.4082, which optimally resolves fluorescence lifetimes in the 2-3 ns range.
>> IRFTransform = calculateIRFTransform(IRFProb, 25/4096:50/4096:25-25/4096, Omega); % A ClockFrequency of 20 MHz and an ADCResolution of 4096 were used
- Load the .mat file containing
IRFProb
,Omega
, andIRFTransform
in the third step of the pipeline. This pipeline lays out all the steps in a typical data analysis: intensity thresholding (for localized fluorophores), phasor plot-based pixel filtering, region exclusion (manual correction), and fitting (using the MATLAB nonlinear optimization functionfmincon
). UseRun Section
to perform your analysis in a guided, step-by-step manner. To minimize$\chi^2$ (the correct way for the Poisson process but numerically problematic due to the low and noisy event counts at the two tails) during fitting, setFittingOption = Fitting2
(for a 2-component exponential decay fit) orFittingOption = Fitting1
(for a mono-exponential decay fit) in the fifth step. To minimize the Manhattan distance between the fitted curve and raw data (a practical way employed in our data analysis, which is further explained here), setFittingOption = Fitting2S
(by default; for a 2-component exponential decay fit) orFittingOption = Fitting1S
(for a mono-exponential decay fit) in the fifth step.
All fitting parameters are automatically saved into a .mat file and two associated plots (including an overlay of the raw FLIM data in black and the fitted curve in red, as well as a residual plot tiled together with a plot of the auto-correlation function of residuals) are also automatically saved as individual .fig files.
Exemplary output figure 1: An overlay of the raw FLIM data in black and the fitted curve in red. The "microtime" refers to the interval between the arrival of a detected photon and its corresponding excitation pulse.
Exemplary output figure 2: A residual plot tiled together with a plot of the auto-correlation function (ACF) of the residuals in the exemplary output figure 1 above.
To demonstrate how fluorescence lifetime measurements can quantify the FRET efficiency, consider a large number of donor fluorophore molecules with a lifetime of
The total donor fluorescence intensity is
wherein the pre-exponential factor
Here,
The effective lifetime of the donor fluorophore (which can be measured through FLIM) becomes
and the total donor fluorescence intensity becomes
wherein
Below is an exemplary FLIM experiment (performed by Dr. Ajit Joglekar) from the reference study where a high FRET efficiency was observed. Here in the experimental group (red), the nuclear pore complex protein NUP50 is tandemly tagged by mNeonGreen and mScarlet-I. In the control group (black), NUP50 is only tagged by mNeonGreen. We can see that the lifetime of mNeonGreen in NUP50-mNeonGreen-mScarlet-I is greatly shortened compared to the lifetime of mNeonGreen in NUP50-mNeonGreen. This results from the highly efficient FRET between mNeonGreen (the donor fluorophore) and mScarlet-I (the acceptor fluorophore) that are closely linked. And by fitting the FLIM data with (multi-component) exponential decays, we can quantify the FRET efficiency using the equation above.
The scatter plots show the mean FLIM data (PDF: the empirical distribution function of the photon arrival microtime relative to the corresponding excitation pulse) of mNeonGreen (red: NUP50-mNeonGreen-mScarlet-I; black: NUP50-mNeonGreen) and each error bar represents the standard deviation of each microtime bin (the number of cells
The phasor transform is a normalized Fourier transform that converts time-resolved emission data into a single point in the complex plane. For a donor fluorophore with an exponential decay
On the complex plane, the corresponding phasor
For an ensemble of fluorophores with different exponential decay lifetimes, it is easy to derive that the phasor is a linear combination of the phasors of the composing species
Since the fractional intensity
The actual FLIM signal
If we ignore the background baseline term and the noise term, the phasor transform of the raw time-resolved emission data
wherein the denominator
and the numerator
Therefore,
The calculateIRFTransform
function calculates
A critical application of the phasor plot in the pipeline is to identify pixels with mostly autofluorescence, without performing tedious fitting pixel by pixel. Autofluorescent substances in the green channel in mammalian cells (see Aubin, 1979 and FPbase Fluorescence Spectra Viewer) typically feature shorter lifetimes (see Horilova, Cunderlikova, and Chorvatova, 2014, Kalinina et al., 2021, and Morrow et al., 2024) than the green fluorescent protein chosen for the FLIM experiment does. For example, in the reference study, we used mNeonGreen (with a reported lifetime of ~ 3 ns) and mScarlet-I as the donor and the acceptor fluorophores (based on McCullock, MacLean, and Kammermeier, 2020), respectively. Reflected on the phasor plot, the phasor transforms of time-resolved emission data registered to pixels with mostly autofluorescence are well separated from those registered to pixels with mostly fluorescent proteins.
... Without careful consideration of the nature of the problem, deconvolution as an information-improving device can easily become an exercise in self-delusion. — Knight and Selinger (1971)
The multi-component exponential fit is intrinsically flexible. Regarding this, Grinvald and Steinberg (1974) raised two educational examples which are reproduced here:
Therefore, one should always be vigilant when translating the parameters derived from a multi-component exponential fit into actual physical parameters. Regardless, the average lifetime of a good multi-component exponential fit weighted by the corresponding FittingOption = Fitting1S
or FittingOption = Fitting2S
, and
[As a side note, Leica TauSense uses the sum of a continous spectrum of exponential decays for FLIM fitting, which is not implemented here. Related info can be found in James and Ware, 1986 and Mérola, Rigler, Holmgren, and Brochon, 1989.]
Our current protocol uses a mirror on the sample plane to measure the IRF. The emission filter is removed and internal reflection at lenses is observed as a prepulse in the measured IRF. This prepulse is an artifact due to the removal of the emission filter and should be manually removed in postprocessing. Additionally, the APD detector has an afterpulse feature (see Ziarkash et al., 2018). This is intrinsic to the detector and an integral part of the IRF that should NOT be removed in postprocessing.
For a scanning microscope, the amplification factor is determined by the setup and the objective but not by the scale of the FOV set for scanning. Therefore, the power of the excitation light on the sample plane and the corresponding area on the sample plane of the APD detector are not affected by the FOV set for scanning. Because the pixel dwell time is fixed, the event count per pixel is directly comparable, regardless of the scale of the FOV set for scanning.
This toolkit is licensed under . If you have used any of the codes in your research, please kindly consider citing the following reference study:
Chen, C., Piano, V., Alex, A., Han, S.J.Y., Huis in ’t Veld, P.J., Roy, B., Fergle, D., Musacchio, A., and Joglekar, A.P. "The structural flexibility of MAD1 facilitates the assembly of the Mitotic Checkpoint Complex." Nature Communications 14 (2023): 1529. (https://doi.org/10.1038/s41467-023-37235-z)
BibTeX:
@article{chen2023,
title={The structural flexibility of MAD1 facilitates the assembly of the Mitotic Checkpoint Complex},
author={Chen, Chu and Piano, Valentina and Alex, Amal and Han, Simon JY and Huis in ’t Veld, Pim J and Roy, Babhrubahan and Fergle, Daniel and Musacchio, Andrea and Joglekar, Ajit P},
journal={Nature Communications},
volume={14},
pages={1529},
year={2023}
}
I would like to thank Dr. J. Damon Hoff (the SMART Center at the Univerisity of Michigan, Ann Arbor) for his suggestions on the manual and his ground-laying contributions to scripts for parsing raw data files.