Translator Disclaimer
13 May 2020 Fluorescence lifetime imaging microscopy: fundamentals and advances in instrumentation, analysis, and applications
Author Affiliations +

Significance: Fluorescence lifetime imaging microscopy (FLIM) is a powerful technique to distinguish the unique molecular environment of fluorophores. FLIM measures the time a fluorophore remains in an excited state before emitting a photon, and detects molecular variations of fluorophores that are not apparent with spectral techniques alone. FLIM is sensitive to multiple biomedical processes including disease progression and drug efficacy.

Aim: We provide an overview of FLIM principles, instrumentation, and analysis while highlighting the latest developments and biological applications.

Approach: This review covers FLIM principles and theory, including advantages over intensity-based fluorescence measurements. Fundamentals of FLIM instrumentation in time- and frequency-domains are summarized, along with recent developments. Image segmentation and analysis strategies that quantify spatial and molecular features of cellular heterogeneity are reviewed. Finally, representative applications are provided including high-resolution FLIM of cell- and organelle-level molecular changes, use of exogenous and endogenous fluorophores, and imaging protein-protein interactions with Förster resonance energy transfer (FRET). Advantages and limitations of FLIM are also discussed.

Conclusions: FLIM is advantageous for probing molecular environments of fluorophores to inform on fluorophore behavior that cannot be elucidated with intensity measurements alone. Development of FLIM technologies, analysis, and applications will further advance biological research and clinical assessments.



Fluorescence microscopy is a core biomedical imaging tool that provides high-resolution images of molecular contrast in living samples. Stokes coined the term “fluorescence” in 1852 for the “remarkable phenomena of light” observed in the materials that emitted light at a different color than the absorbed light.1 At that time, refraction (or refrangibility), internal dispersion, circular dichroism, and other phenomena of light were well studied. Stokes examined flower petals, leaves, turmeric, calcium fluoride, and many other compounds. In the 20th century, fluorescence was redefined as a short-lived emission of photons caused by the incidence of higher energy photons and became a popular tool for studying molecular dynamics and characterizing compounds.

In the mid-20th century, Weber used fluorescence properties of molecules such as depolarization along with absorption and emission spectra to pinpoint molecular dynamics and reveal kinetic parameters for biologically relevant processes such as enzyme binding.2 These experiments advanced fluorescence as a major means for biophysical and biochemical investigation. By the late 20th century, numerous brightly fluorescent small molecules had been categorized and repurposed as markers bound to other molecules. Fluorescence-based targeting provided unique molecular specificity in high-resolution microscopy. For example, mitochondria were identified using a small fluorescent molecule called tetramethyl-rhodamine-ethyl-ester (TMRE), which binds only to the mitochondrial membrane. This discovery has greatly advanced the study of energy distribution in biological systems. A weaker endogenous source of fluorescence was also identified within mitochondria. This inherent ability of many biological systems to fluoresce without the addition of external fluorophores was termed “autofluorescence.”

Along with reduced pyridine nucleotides, oxidized flavins, and other metabolic agents, proteins containing an abundance of amino acids, such as tryptophan, phenylalanine, and tyrosine, are the major endogenous fluorophores in biological systems. Studies in the 1980s identified the presence of a fluorescent protein expressed by jellyfish. This small protein was cloned into a functionally expressible green fluorescent protein (GFP)3 and was genetically expressed in Escherichia coli (E. coli) to create bacteria capable of green fluorescence.4 Thousands of varieties of this protein that fluoresce at different parts of the spectrum have been engineered, and their fluorescence behavior has been modified to cater to unique probing interests. The in vivo imaging capabilities of GFP-tagged proteins within organisms have bolstered fluorescence imaging as a robust and flexible assessment method for biomedical research.

Fluorescence lifetime imaging microscopy (FLIM), which exploits the lifetime property of fluorescence, is a microscopy technique that has gained popularity because of its high sensitivity to the molecular environment and changes in molecular conformation. FLIM has been extensively used in autofluorescent molecular imaging to study cellular metabolism. FLIM of autofluorescent molecules provides unique insights into cellular health in a nondestructive manner and is often used to study live animals and as a contrast mechanism for fluorescence-guided surgery.512 Exogenous fluorescent molecules that are capable of monitoring microenvironmental parameters, such as temperature, viscosity, pH, and ion concentration, are categorized as FLIM-based sensors.1315 Protein–protein interactions can be monitored using Förster resonance energy transfer (FRET) sensors that are specific for cellular signaling, cellular proliferation, cytokinesis, and other molecular interactions.1619 Thus, leveraging both endogenous and exogenous fluorophores, FLIM can monitor numerous processes in cells and tissues, including disease progression and drug efficacy.

In this review, we discuss the principles and theory behind FLIM and its unique advantages over intensity-based fluorescence microscopy methods. Then, we review FLIM instrumentation and FLIM analysis methods including segmentation and population density modeling of cell heterogeneity, and we close with a summary of FLIM applications in vivo and in vitro.


Fluorescence Lifetime

When a molecule in ground state (denoted as S0 in Fig. 1) absorbs light of energy equal or greater than the higher energy levels (S1,S2,,Sn), an electron is excited to a higher energy level for a short period. The electron will undergo vibrational relaxation to the lowest vibrational level of the excited state (denoted as S1) by a nonradiative process called internal conversion. From the S1 electronic state, molecules return to the ground state either by a radiative or nonradiative process. Figure 1 represents the different luminescence phenomena that occur in these levels.

Fig. 1

Schematic of Jablonski diagram.


Fluorescence is a radiative process in which molecules (fluorophores) decay to the ground state by emitting detectable photons (on the timescale of 109  s). The fluorescence emission happens from the lowest excited electronic level (S1). This mandatory emission from the lowest excited electronic level ensures that the emission spectrum remains the same and is independent of the excitation wavelength. The energy of the emitted fluorescence photon is lower (i.e., emission occurs at a longer wavelength than the excitation) due to energy loss in vibrational relaxation and internal conversions. This shift in emission wavelength is referred to as the Stokes shift. Another predominant luminescence process, phosphorescence, occurs when the excited electron energy transitions into a triplet energy level (T1,T2,,Tn) by a process known as intersystem crossing (ISC). Electrons in the triplet state have parallel spins, and these electron transitions are “spin-forbidden,” resulting in a slow transition to ground level by emission of a phosphorescence photon or reversal of the ISC and emission of a delayed fluorescence photon. Phosphorescence occurs on timescales on the order of milliseconds to hundreds of seconds. The Jablonski diagram shown in Fig. 1 concisely illustrates these processes.

The quantum yield of the molecule is defined as the ratio of emitted photons to the absorbed photons. Quantum yields for common fluorescent compounds include 80% for fluorescein,20 60% for eGFP,21 6% for tryptophan,22 and 2% for reduced nicotinamide adenine dinucleotide (NADH).23 This emission efficiency of a molecule depends on (1) its spatial orientation with respect to the incident electromagnetic wave’s electric field orientation (polarization), (2) the electronic energy levels available for absorbing the incoming photon energy (absorption spectrum), (3) the efficiency of rearrangement of vibrational levels (fluorescence lifetime), (4) relaxation back to the ground state electronic energy level (Stokes shift), and (5) the population of vibrational levels within this ground state (emission spectrum). Fluorophores are characterized by their absorption spectrum, fluorescence lifetime, Stokes shift, and emission spectrum.

Conventionally, we define fluorescence lifetime (τ) as the average time that a fluorophore remains in its excited state. In this interval, the intensity I(t) decreases to 1/e or 36.8% of its original value. The decaying intensity at time t is given by a first-order kinetics equation summed across all species, i, in the sample

Eq. (1)

where α is the pre-exponential factor or the amplitude of the exponential function. The mean lifetime (τm) of a multiexponential mixture of species is the sum of each species lifetime (τi) weighted by fractional contribution of each species (αi)

Eq. (2)

In addition, the number of excited molecules at a time t is given as

Eq. (3)

where n(t) is the number of molecules in the excited state at time t.

Fluorescence lifetime can be measured in either the time-domain or frequency-domain, and these methods will be covered in detail in Secs. 2 and 3. Briefly, for time-domain methods, the sample is excited by a short excitation pulse and the decay is calculated either from time-of-arrival of photons that are binned into a histogram or by time-gated detection or pulse sampling techniques. If multiple fluorescent species are present, all species are summed into a single histogram. In frequency-domain methods, each photon is represented by its phase delay with respect to the excitation photon, which is similar to the arrival time histogram. For multiple species, this phase distribution is analyzed in Fourier space to extract the modulation and demodulation parameters that separate multiple species. Both time-domain and frequency-domain offer unique advantages and challenges in different FLIM scenarios including low photon budget imaging, high dynamic range imaging, or high time resolution imaging.


Autofluorescence FLIM Measurements

Biological systems are rich in endogenous fluorophores that are used for autofluorescence molecular imaging in a convenient, label-free manner. Endogenous fluorophores are powerful biomarkers because their emission properties are often influenced by their microenvironment, as well as the morphology, metabolic state, and pathological conditions of the sample. Notable endogenous fluorophores along with their excitation and emission wavelengths and fluorescence lifetimes are listed in Table 1. Imaging endogenous fluorophores is advantageous because it avoids the administration of external fluorescent dyes, thus circumventing complications introduced by these contrast agents including nonspecific binding, toxicity, and interference with the biochemical and physiological functions of the sample. Furthermore, autofluorescence imaging can be easily translated to in vivo monitoring in animal models and in humans for impactful clinical measurements.

Table 1

Spectral characteristics and lifetimes of endogenous fluorophores.

Endogenous fluorophoreExcitation (nm)Emission (nm)Lifetime (ns)Reference
Metabolic coenzymes
NAD(P)H free340 (max)470 (max)0.4 (free), 1 to 5 (bound)2425.26.27.28
FAD, flavin450 (max)535 (max)2.3 to 2.9 (free), <0.1  ns (bound)26, 29, 30
Flavin mononucleotide (FMN)444 (max)558 (max)4.27 to 4.6731, 32
Structural proteins
Collagen280 to 350370 to 4400.2 to 0.4, 0.4 to 2.532, 33
Elastin300 to 370420 to 4600.2 to 0.4, 0.4 to 2.532, 33
Retinol327 (max)510 (max)1.8, 5.0 (free), 0.7, 3.6, 12 (bound)26, 34
Riboflavin420 to 500520 to 7504.1232
Vitamin B6330 (max)420 (max)0.6 to 8.435, 36
Vitamin K335 (max)480 (max)26
Vitamin D390 (max)480 (max)26
Vitamin B12275 (max)305 (max)26
Melanin300 to 800440, 520, 5750.1 to 0.2, 0.5 to 1.8, 7.93233.34
Eumelanin3555200.058, 0.51, 2.9, 737, 38
Keratin277 (max)382 (max)1.439, 40
Protoporhorphyrin IX400 to 450630, 690, 7109.7 to 1626, 41
Lipofuscin340 to 395540, 430 to 4601.3432, 35
Bilirubin350 to 520480 to 6500.02 to 0.09, 1 to 242, 43
Amino acids
Phenylalanine258 (max)280 (max)7.532
Tryptophan280 (max)250 to 3103.0332
Tyrosine275 (max)300 (max)2.532


FLIM of NAD(P)H and FAD for metabolic imaging

Nicotinamide adenine dinucleotide (NAD) and flavin adenine dinucleotide (FAD) are two metabolic coenzymes that play a myriad of roles in cellular oxidation and reduction reactions. The reduced form NADH and oxidized form NAD+ are involved in mitochondrial function, energy metabolism, calcium homeostasis, gene expression, oxidative stress, aging, and apoptosis. The reduced NAD phosphate (NADPH) is associated with reductive fatty acid biosynthesis, steroid biosynthesis, oxidative stress, and antioxidation, while the oxidized form (NADP+) is involved with calcium homeostasis.9 Real-time monitoring of cellular metabolism during pathophysiological changes is possible by measuring the redox ratio (NADH/NAD+). NADH is the principal electron acceptor in glycolysis, which results in two NADH molecules per glucose molecule. The Krebs cycle also reduces NAD+ to NADH in three of its reactions. During oxidative phosphorylation, NADH is oxidized to NAD+ by donating electrons to the electron transport chain, and these electrons are ultimately accepted by oxygen.8,9 In the case of anaerobic glycolysis, NAD+ is converted to NADH and oxidative phosphorylation is diminished, which creates an overall increase in NADH abundance. Thus, the reduction–oxidation pair NADH/NAD+ serves as an indicator of balance between oxidative phosphorylation and glycolysis. Flavins such as FAD are also involved in cellular oxidation–reduction reactions. The reduced form (FADH2) is oxidized to FAD in complex II of the electron transport chain, while FAD is reduced to FADH2 in pyruvate decarboxylation and the Krebs cycle.

NADH and FAD are fluorescent while NAD+ and FADH2 are not. The fluorescence of NADH and NADPH are difficult to distinguish, and their combined fluorescence is referred to as NAD(P)H. Due to the pivotal role of NADH, NADPH, and FAD in cell biology and metabolism, these endogenous fluorophores have been used to monitor cellular redox reactions, energy metabolism, and mitochondrial anomalies under different pathophysiological conditions. Chance and others in the 1980s established NAD(P)H and FAD fluorescence for metabolic imaging.4448 The use of endogenous fluorescence enables noninvasive metabolic imaging of cells and tissue in their native physiological environment without perturbations associated with contrast agents. After the development of FLIM instrumentation, biophysicists began to relate the fluorescence lifetimes of NAD(P)H and FAD to cellular metabolism.24,25,29 The fluorescence lifetime of NAD(P)H is significantly shorter in the free state (400  ps) compared with the protein-bound state (1 to 5 ns) of the molecule.24,25,27 This is due to quenching in the free state as the NAD(P)H molecule folds and diminished quenching in the protein-bound state as the NAD(P)H molecule extends. Conversely, FAD has a longer lifetime in its free state (2.3 to 2.9 ns) compared with its protein-bound state (<0.1  ns).29,30,49,50 Bird et al. used FLIM to demonstrate a correlation between the redox ratio (NADH/NAD+) and the relative amounts of free to protein-bound NAD(P)H.51


FLIM-FRET Microscopy

The fluorescence lifetime of a donor fluorophore changes when it undergoes FRET with an acceptor molecule. As a result, FLIM can visualize changes in the proximity of FRET pairs.17,18,52 Specifically, the quenching of the donor emission by FRET leads to a decrease in its lifetime. FRET has been used to detect conformational changes within proteins, receptor/ligand interactions between proteins, hybridization or splitting of nucleic acid strands, membrane lipid interactions and distributions,16 the activity of proteases, chromatin architecture,53 and many other phenomena. Genetically engineered FRET pairs can be strategically expressed in biological systems for any application in which distances between proteins or protein subdomains are of interest (Fig. 2). A detailed review of FRET can be found elsewhere.54

Fig. 2

FLIM-FRET concept and applications. (a) At the Förster distance (R0, defined by the specific donor–acceptor pair), the efficiency of energy transfer between donor and acceptor is 50%, such that large distances exhibit low efficiencies of energy transfer. At these large distances, the fluorescence lifetime of the donor (τ) is not affected by FRET. (b) As the distance between donor and acceptor decreases, FRET can occur, quenching the emission of the donor and shortening the donor lifetime decay (τ). (c) Photobleaching of the acceptor can confirm that the change in donor emission lifetime was due to a FRET interaction. (d), (e) Donor and acceptor pairs within the same molecule can be used to detect changes in (d) protein confirmation or (e) cleavage of proteins by proteases. (f), (g) Donor and acceptor pairs in separate molecules can be used to detect (f) receptor/ligand binding and (g) hybridization or splitting of nucleic acid strands.


FLIM-FRET has a number of advantages over intensity-based FRET. In addition to the advantages of FLIM over intensity imaging that will be discussed in Sec. 1.4, there are also benefits specific to FRET interactions. Most importantly, FLIM-FRET only requires the measurement of the donor lifetime, so direct excitation of the acceptor is not needed and acceptors with poor quantum efficiencies can be used. In addition, less excitation intensity is required for FLIM-FRET because wider emission filters can be used, allowing for FLIM-FRET pairs that are less photostable. 55 A portion of the donor fluorophores can fail to excite in some FRET experiments, which introduces additional errors in intensity-based FRET that are avoided in FLIM-FRET. Finally, FLIM-FRET and multiexponential fitting can be used to quantify the proportion of quenched and unquenched donor molecules.18 FRET events can be confirmed by photobleaching the acceptor, which should result in a donor lifetime at pre-FRET levels. The drawbacks of FLIM-FRET versus intensity-based FRET mirror those of using FLIM in general with the additional stipulation that a carefully measured reference lifetime value for the donor alone (without acceptor present) is needed for accurate calibration.

A number of considerations need to be made to select FRET pairs (FPs) that are specifically suited for use with FLIM-FRET (Table 2). Some pairs that are not useful for intensity FRET are quite useful for FLIM-FRET. This is because intensity FRET requires spectral overlap between donor and acceptor, while FLIM-FRET is optimized using pairs with well-separated emission spectra.55,69 A donor fluorophore will ideally have a long, monoexponential decay. With simple decay kinetics, it is easier to determine the distinct lifetime of the quenched donor using multiexponential fitting of the decay data.55 The unquenched single lifetime of the donor should also be as long as possible to optimize the dynamic range of the FP. Donors should also have high photostability and should not photoconvert, which could lead to an overestimation of quenching. Acceptors should have a high absorbance coefficient but an extremely low quantum yield to avoid acceptor emission in the donor channel. By choosing an optimal acceptor, donor emission can be collected using a wider spectral window to increase signal while reducing excitation power. In addition, another probe could be added to the acceptor spectral window to correlate FRET interactions with the behavior of another labeled protein.55

Table 2

Examples of FRET pairs for FLIM-FRET imaging.

mTurquoiseYFP56, 57
mEGFPmCherry62, 63

New approaches to improve the FLIM-FRET continue to emerge. One goal of recent efforts is fast FLIM-FRET to capture rapid cellular events, increase imaging throughput, and quickly acquire volumes of three-dimensional (3-D) biological interactions. For example, Poland et al.70 developed a multifocal multiphoton system that simultaneously acquires multiple planes of FLIM-FRET using an array of beamlets produced by a spatial light modulator. Other techniques employ many parallel detectors to rapidly image protein–protein interactions in live cells.71 Recent developments in FLIM-FRET analysis techniques provide highly localized information on molecular interactions. Phasor analysis (see Sec. 3.1.2) of FLIM-FRET data has quantified chromatin organization at the nucleosome level, which is below the diffraction limit of most imaging modalities.72 Here, increased FRET between fluorescent histones signals an increase in nucleosome proximity.


Advantages of FLIM Over Intensity Imaging

FLIM offers many unique advantages over intensity-based fluorescence microscopy. Fluorescence intensity imaging provides information on the spatial distribution of fluorophores and can discriminate between fluorophores with distinct spectral properties. However, intensity alone cannot distinguish fluorophores with similar spectra or distinguish unique molecular environments around the same fluorophore. FLIM can frequently discriminate spectrally overlapping fluorophores using the fluorescence lifetime. For example, NAD(P)H often appears indistinguishable in different cellular environments based on fluorescence intensity and spectral information, but it can be distinguished easily using FLIM (Fig. 3). Overall, FLIM is advantageous in its ability to detect changes in the molecular environments of fluorophores to provide information about fluorophore function and behavior that could not be elucidated with intensity measurements alone.74

Fig. 3

FLIM provides metabolic contrast in 3-D tumor organoids treated with chemotherapy. The τm of endogenous NAD(P)H is sensitive to the metabolic response to chemotherapy in patient-derived pancreatic cancer organoids. NAD(P)H intensity measurements alone did not distinguish treatment. Here, τm is calculated from a two-exponential decay of the free and protein-bound lifetimes of NAD(P)H. Scale bar=50  μm. Adapted with permission from Ref. 73.


Unlike intensity-based measurements, FLIM is largely independent of fluorophore concentration. This means that FLIM can determine whether a change in fluorescence intensity is due to changes in quantum yield (e.g., fluorescence quenching), a variation in the overall concentration of the fluorophore, or both. FLIM measurements are also less vulnerable than intensity measurements to inner filter effects, which are absorption and scattering events that modulate the detected fluorescence intensity. Therefore, FLIM is well-suited for accurate measurements of quenching dynamics.75,76 Multiple configurations or states of a fluorophore can be detected with FLIM at a single location or pixel. For example, both bound and unbound fluorophores, as well as proteins with distinct folding states, will have different molecular environments that coexist within the same pixel.

FLIM is a self-referenced measurement (i.e., independent of absolute detected intensity), so FLIM experiments do not require the throughput calibration steps that are needed for intensity-based experiments. Lifetime is an absolute measurement that can be repeated across numerous device configurations (e.g., excitation intensity, detector sensitivity, and path length) after accounting for the instrument response function (IRF) of that device. Thus, artifacts caused by nonuniform illumination, which would greatly affect intensity measurements, are mitigated by measuring the lifetime. This internal calibration has the added advantage of making FLIM experiments more reproducible and comparable between different instrumentation configurations. In addition, lifetime measurements are independent of excitation and emission light scattering in cells and tissues, provided that any time delay is smaller than the resolution of the timing electronics. In addition, confounding scattering profiles can be modeled in a fitting routine by assuming a Gaussian spread of the IRF function. This is conventionally used in commercial time correlated single photon counting (TCSPC) packages such as SPCImage77 (Becker & Hickl), allowing FLIM to be performed accurately at deeper penetration depths.

The key downsides to time-resolved measurements include a long acquisition time that may prevent visualizing fast events, requirements for time-resolved electronics and accurate IRF measurements, and sensitivity to changes in temperature, pH, and viscosity that complicate data interpretation. Promising new techniques have increased the speed of FLIM, which should enable visualization of fast dynamics in the future.7880



FLIM measures the fluorescence decay rate of a fluorophore on the timescale of subnanoseconds to hundreds of nanoseconds. For reference, light travels at a speed of 3×108  m/s or 1 foot in 1 ns. Fast electronics coupled with efficient photon detectors have been integral tools for FLIM and other fast temporal measurements.

Time-domain and frequency-domain FLIM measurements are overviewed in Fig. 4, with detailed descriptions below. Briefly, time-domain fluorescence lifetime measurements use a short pulse of light for excitation (short relative to the lifetime of the sample) and then record the exponential decay of fluorescent molecules either directly (i.e., by gated detection or pulse sampling) or using time-resolved electronics that bin photons by their arrival times [Figs. 4(a) and 4(b)].74,8184 Alternatively, frequency-domain techniques can also measure fluorescence lifetimes [Figs. 4(c) and 4(d)].85,86 Here, the excitation is continuous with amplitude modulation over time as a sine wave. The fluorescence signal shifts in phase and amplitude with respect to the excitation wave. The phase delay and amplitude modulation for a fluorophore are visualized by plotting the phase changes over a range of modulation frequencies [Fig. 4(d)]. This resulting fluorescence sine signal can be demodulated in the frequency-domain to quantify the delay induced by the exponential decay of the fluorescence intensity.

Fig. 4

Schematic of time-domain (TCSPC) and frequency-domain FLIM. (a) TCSPC FLIM acquisition includes the short excitation pulse, single exponential fluorescence decay curve, and lifetime (τ) defined at the 1/e value. Inset shows detected single fluorescence photons (red circles) at different time periods within multiple excitation pulses. (b) Photon time of arrival histogram built from the detection time of multiple fluorescent photons (red circles); green line represents the IRF, and dotted red line represents the fit function. (c) Schematic diagram of frequency-domain measurement with sinusoidally modulated excitation (exc) and the resulting phase shifted emission (em) signal. The AC and DC components of each signal are also indicated. (d) Modulation and phase versus frequencies for different lifetimes. TM, modulation lifetime; TP, phase lifetime.


The most common implementation of FLIM is with a fast electronic method called TCSPC [Fig. 4(a)]. In TCSPC, a fast stop-watch measures the time between an excitation photon and emission photon. This time defines each emission photon’s time-of-arrival. The fast clock time is experimentally measured with a time–amplitude converter circuit (TAC), which converts the photon time-of-arrival to an analog voltage that can be recorded. In conventional TCSPC, at high photon count rates, most of the incoming photons will not be measured due to the instrument dead time. This will lead to the pile-up effect where only the photons with shorter arrival times will be recorded per excitation pulse. This loss of photons with longer arrival times will create an incorrect photon histogram, leading to overall shortening of the measured fluorescence lifetime. To avoid these effects, a low photon count at the detector is desirable, ideally <10% of the excitation repetition rate. Thus, in general, time-domain methods detect one fluorescence photon across several excitation pulses, so many excitation pulses are required to build a histogram [Figs. 4(a) and 4(b)]. The signal-to-noise ratio (SNR) of FLIM measured by photon counting, (SNRFLIM), depends on the number of photons detected per pixel (N), such that it changes with the square root of N:87

Eq. (4)

Therefore, to improve the SNRFLIM, the photon detection process is repeated thousands of times to generate a distribution of time-of-arrivals of fluorescence photons, which is the measured exponential fluorescence decay. This improves accuracy of FLIM data analysis by curve fitting, which will be discussed in Sec. 3.

Practically, TCSPC employs an efficient method known as reverse-TAC to measure the time between the emission photon and the next excitation photon. After the histogram is made with reverse-TAC, the time axis is inverted. Reverse-TAC mode is advantageous over the forward-TAC mode in systems with high repetition rates, such as laser sources in the range of 50 to 100 MHz, but with low photon count rates. In reverse-TAC mode, TAC is reset only when a photon is detected, using the reset signal from the consequent laser pulse, thus avoiding the requirement of additional TAC reset circuits. The measured decay is a convolution of the excitation pulse and fluorescence decay [Figs. 4(a) and 4(b)]. For an ideal delta excitation pulse, the measured fluorescence decay would equal the actual fluorescence decay. The alternative to TAC is a time–digital converter (TDC), which converts the time to a digital value of delay. Experimentally, both TAC and TDC are realized using a programmable logic gate array called a field programmable gate array (FPGA), but TAC and TDC use different electronic means to calculate the time delay between the excitation photon and detected photon. TAC systems are often credited with higher quality timing due to reduced timing jitters (error in timing estimation). Photon detection is usually achieved with a constant fraction discriminator (CFD) circuit to read the analog voltage output from the detector. This discriminator circuit determines the photon counts and triggers the stop clock for TCSPC. These methods are discussed in detail in previous publications.74,82

One of the main limitations of TCSPC is slow acquisition speed, which motivates new techniques to expedite FLIM imaging. TCSPC and other FPGA-based architectures have long dead times (tens of nanoseconds) between photon detection events.74,88 Therefore, recent adoption of high-speed digitizers (subnanosecond sampling) aim to decrease dead times to two or less nanoseconds between photon detection events, so fluorescent decays can be generated more rapidly. These digitizers have been used for FLIM of endogenous and exogenous fluorophores in cells and tissues7880 and offer great promise for shorter FLIM acquisition times. Other strategies to increase speed focus on minimizing dead times, parallelizing TCSPC, and implementing TCSPC with multifocal excitation.8894 A detailed discussion on challenges and current approaches to improve FLIM can be found in a previous review.95

Other time-domain methods include time-gating (TG) and pulse sampling. In TG FLIM, following a short excitation pulse, the fluorescence decay is directly sampled at two or more time gates that are sequentially delayed from the excitation pulse.9699 For a single exponential decay, the lifetime (τ) can be calculated using two equal time gates at Δt time separation83,100

Eq. (5)

where I1 and I2 are the intensities measured at the two gates, respectively. For multiple fluorophores, however, two time gates would yield only a mean lifetime. Thus, multiple precisely synchronized gates, in combination with decay analysis techniques such as exponential fitting and phasor approach, are employed for multiexponential lifetime calculation.83,98,101,102 Following the first demonstration of multiphoton laser scanning TG FLIM by Sytsma et al., the technique has been employed in multiple studies.103105 The TG approach has been more widely adopted for time-domain wide-field FLIM and will be discussed in Sec. 2.1.1. TG is advantageous over TCSPC because the electronics dead time is significantly lower (about 1 ns) and the high count rates result in faster acquisition times.98 However, TG lacks the sensitivity and time resolution of the TCSPC approach.106,107 Moreover, limitations in the number of gates and counters might result in undersampling of the decay.87 Another consideration is the convolution of the IRF with the intensity measured at the first time gate. This can be solved by either correcting for the IRF or delaying the first time gate after the IRF. However, this delay causes additional loss in photon efficiency.87 For a detailed discussion on TG FLIM and comparison with TCSPC, please refer to prior publications.98

Similar to TG, the time-domain pulse sampling approach involves direct measurement of the decay signal after a short excitation pulse; it was first demonstrated by Steingraber and Berlman in 1963.108 The fluorescence signal is continuously measured by a fast response detector and sampled by fast digitizers.84,106 The decay is reconstructed from all detected photons for a single excitation. The time-domain approach was traditionally employed for nonimaging point measurements of time-resolved fluorescence spectroscopy.109112 Recent developments include pulse sampling coupled with optical fiber-based multispectral fluorescence lifetime imaging (FLIm).113117 The fast acquisition speeds make this approach attractive for spectroscopy in clinical applications. Furthermore, the measured fluorescence signal is not affected by background light, allowing data acquisition in a clinical setting including operating rooms without the need to dim or turn off the room light.11,116 It has also been implemented with other imaging techniques such as optical coherence tomography118 or ultrasound119 for bimodal imaging. Current implementation of FLIm includes microchannel plate (MCP) detectors and a high-speed (12.5  GS/s) digitizer.116 In a recent development, a second pulsed laser was multiplexed in time with the excitation laser traditionally used for endogenous fluorophores. This allowed additional exogenous fluorophores to be imaged.117 MCPs have a high response time but low gain; thus they perform better with samples with high quantum efficiencies.74,106 The pulse sampling technique also suffers from uncertain accuracies in the lifetime decay curve reconstruction since the instrument noise characteristics are unknown.106 Pulse sampling techniques are discussed in detail in prior publications.84

Frequency-domain lifetime measurements date back to 1927 when, for the first time, lifetime was measured by Gaviola.120 Venetta in 1959 demonstrated lifetime measurements by coupling a phase fluorometer to a microscope, one of the seminal works leading to present day FLIM.121 In frequency-domain FLIM, the sample is excited by an amplitude modulated light source at high frequencies (MHz), and the harmonic response of the system is measured [Fig. 4(c)].122 The equation of this modulated excitation signal for a given frequency of modulation can be written as123

Eq. (6)

where E(t) is intensity at time t and E(0) is intensity at time=0. ME is the excitation modulation factor, and ω is the angular frequency and is given by ω=2πf, where f is the linear modulation frequency. With a sinusoidal excitation, the emission signal will also be modulated sinusoidally.122,124 However, the emission signal will be phase shifted with respect to the excitation due to delay between the absorption and emission. This can be written as

Eq. (7)

where F(t) is the fluorescence intensity at time t and F(0) is that at time t=0.124 MF is the emission modulation factor, and ϕ is the phase delay between excitation and emission. The modulation and phase shift of the emission is dependent on the relative values of the frequency of modulation, f, and lifetime τ [Fig. 4(d)]. In the case of single exponential decays, the phase lifetime (τP) is

Eq. (8)

The modulation factors can be expressed as the ratio of AC to DC components of the respective excitation (EX) and emission (EM) signals

Eq. (9)

From this, we can estimate the relative modulation M as

Eq. (10)

The relative modulation (M) and modulation lifetime (τM) are related by

Eq. (11)

In the case of a single exponential decay, lifetimes from phase and relative modulation are equal, (i.e., τP=τM for all ω). However, for multiexponential decays, τP<τM and absolute values will depend on the modulation frequency.123 The angular modulation frequency should be set to roughly the inverse of the lifetime (i.e., ωτ=1) to give maximum sensitivity. Linear modulation frequencies of 100  MHz to 1 GHz give picosecond temporal resolution, which is appropriate for fluorescence lifetime measurements. Given these constraints on modulation lifetimes in frequency-domain measurements, phase lifetimes are preferred to modulation lifetimes.

One major advantage of frequency-domain FLIM over time-domain FLIM techniques, such as TCSPC, is acquisition speed, making frequency-domain an ideal technique for measuring rapid cellular events. The slower processing electronics used in TCSPC can also limit the ability to accurately measure lifetime in very bright samples with high photon count rates. Previously described TG FLIM, pulse sampling techniques, and new faster TAC/TDC electronics (2 to 100 ns dead time) have improved current time-domain FLIM acquisition times, bringing them closer to frequency-domain (0 ns electronics dead time). One of the latest advancements includes implementation of frequency-domain FLIM in a multiphoton microscope capable of imaging deeper than conventional systems.125 Finally, frequency-domain FLIM can be implemented without the use of costly pulsed lasers.

On the other hand, TCSPC can provide better timing resolution and higher SNR for weakly fluorescent samples due to its ability to time individual photons. Thus, frequency-domain may be more advantageous for brighter, more dynamic samples, while TCSPC may be beneficial for weakly fluorescent, static samples. These dynamic range limitation errors are pronounced for fitting routines that use spatial binning for increasing accuracy, when the fitting would automatically be biased by the larger number of photons from neighboring pixels. In addition, since individual photons are timed, TCSPC can distinguish between individual components of a multiexponential decay with high accuracy. To resolve multiple components in the frequency-domain, the signal must either be recorded using multiple modulation frequencies123,126 or digital heterodyning aided with phasor analysis techniques.127,128 This is summarized in Table 3.

Table 3

Advantages and limitations of time and frequency-domain FLIM techniques.

FLIM techniqueAdvantagesDisadvantagesReferences
TCSPC• High accuracy of lifetime estimation• Requires costly pulsed lasers85, 87
• Poor performance with high photon count rates or dynamic samples
• Provides better SNR for weakly fluorescent samples
• Estimates multiple lifetime components
Time gating• Lower electronic dead time than TCSPC electronics• Low sensitivity and time resolution87, 98, 106
• High count rates (no count rate limitation)• Poor performance with low photon counts
• Fast acquisition speed• Undersampling of decay signal due to limited gates and counters
Pulse sampling• Fast acquisition speed• Requires costly electronics and pulsed laser11, 74, 84, 106, 116
• Relatively unaffected by background light
• Low gain of detector (MCP)
• Instrumentation noise characteristics unknown
Frequency domain• Fast acquisition speed• Poor performance with low photon counts85, 87
• Performs well with bright samples
• Can be implemented without pulsed lasers



Two microscopy imaging schemes are used for FLIM: laser scanning microscopy (LSM) and wide-field illumination (WFI) microscopy. LSM and WFI are compared in the FLIM imaging schematic shown in Fig. 5. LSM systems are subdivided based on their excitation-detection method as either confocal (CLSM) or multiphoton (MP-LSM). All of these microscopy techniques offer optical sectioning, which allows 3-D FLIM imaging. However, it is important to note that certain clinical applications such as surgical guidance, endoscopy, ophthalmoscopy, and others do not always require sectioning and can work in a topological imaging modality or a lensless imaging scheme.116,129 FLIM imaging optics require a light source for illumination and sensitive detection to distinguish the photons of interest from background photons. However, the key component of FLIM is the electronics used to estimate the timing of detected photons within each pixel. Numerous strategies are in development for more sensitive detection methods.130134

Fig. 5

Schematic showing FLIM implementation in scanning and wide-field configurations. Two imaging modalities are compared side by side: scanning TCSPC FLIM and wide-field TG FLIM. The scanning beam of laser light from a galvanometric mirror is projected onto the back focal plane of the objective lens (OL) using the scan lens (SL) and tube lens (TL). The size of the beam and scan angle is often adjusted by varying the SL-TL pair. The light from the back focal plane is then focused by the objective lens into an excitation light cone. The emission from the same light cone is retraced back through a dichroic beamsplitter into a photon detector unit such as a photomultiplier tube (PMT). The electrical current from the PMT is amplified and read by an electronic board to calculate photon arrival times. These photon times are linked to the pixel of illumination by the computer (PC) that controls the scanner position in the image and thereby produces a histogram of photon arrival times for each pixel as shown in the inset of the left. LSM FLIM typically achieves 4 to 10 frames per second (fps) acquisition speed, which is usually limited by the scanner speed. (Right) Wide-field FLIM requires a focusing lens (FL) to achieve a field of illumination. The fluorescence from the focus of the objective lens is magnified by a tube lens and then imaged onto a camera sensor. FLIM in wide-field systems is achieved using a short frame exposure time (ns per frame). However, wide-field FLIM requires repeated frame acquisitions over a total time of milliseconds to seconds to collect sufficient photons for a complete histogram of fluorescence decay, as shown in the inset on the right side. Recent FLIM cameras intelligently select modulation-demodulation waveforms to achieve faster FLIM frame rates of 15  fps.



Wide-field FLIM

WFI uses a parallel illumination field at the focus of the objective lens and collects fluorescence from the focal plane onto a camera (Fig. 5). Wide-field FLIM is often used for rapidly imaging large sample areas since light from the entire field of view is collected using a camera-based detection. Wide-field FLIM uses either time-domain techniques, such as TCSPC135 or TG, in which a series of fluorescence images are collected by shifting a timing window (order of nanoseconds) through the emission decay,136138 or frequency-domain methods of demodulating the fluorescence signal from the excitation frequency.139

Wide-field FLIM has the advantages of higher frame rates and less photodamage when compared with LSM. However, camera sensitivity and SNR are not as high as that of LSM detectors, which results in poorer axial resolution. In wide-field collection, every camera pixel simultaneously detects scattered light from all other pixels of the illuminated area, which intermixes the timing-spatial coordinates.133 Assuming a fixed photon emission rate from the sample, image optimization is a trade-off between either spatial or temporal resolution.133 Wide-field techniques such as structured illumination and spinning-disk confocal can achieve higher spatial resolution without compromising imaging speed.136,140143 Other advantages of wide-field systems are their simpler implementation and the low computational cost to assign photon detection times in each pixel. Some of the advancements in wide-field FLIM include its implementation with Nipkow disc microscopy for fast 3-D FLIM imaging,142 wide-field coupled with single plane illumination microscopy for high-resolution 3-D FLIM,144 TG single photon avalanche diode (SPAD) cameras for phasor-based high speed wide-field FLIM,145 multifrequency widefield,146,147 and image gating by pockel cells.148 Current wide-field FLIM systems are discussed in detail by Suhling and Hirvonen et al.149


Laser Scanning Microscope FLIM

Laser scanning microscopes have out-of-focus rejection methods that enable higher contrast and spatial resolution than wide-field systems. In comparison with wide-field FLIM, LSM-FLIM modalities are often coupled to faster electronics to generate precise photon detection times per pixel. As discussed previously, these methods fall under two modalities of timing calculation: TDC-based or TAC-based.


Laser sources

FLIM uses modulated laser sources for illumination. This can be achieved by many modern pulsed laser diodes that are modulated with an internal trigger or the high power density, ultrafast laser sources developed in the 1990s. These lasers have a remarkably short pulse duration (hundreds of femtoseconds), durable repetition frequency (in the order of 0.1 GHz), and tunable wavelengths in the near-infrared region. These lasers are currently used extensively for in vivo and small animal imaging due to their use as a multiphoton excitation (MPE) source (explained in the section below). Pulsed light sources are popular because of numerous applications in digital communications and remote sensing. Nonlinear light sources, such as supercontinuum sources, are also popular because they achieve near continuum tunability over a large wavelength range.


Confocal and multiphoton microscopes

Both CLSM-FLIM and MP-LSM FLIM are broadly used in applied sciences to study biology and materials. Confocal imaging methods use a pinhole (small aperture) to reject out-of-focus light. Most biomedical confocal systems use a low power laser for excitation and focus the light at one point in space using a pair of galvanometric scanners (XY scanner). Precise movement of the objective controls for the Z position. The fluorescence emission from the 3-D focal volume retraces back through the XY scanner (thus descanned) and reaches the detector. The focused spot is scanned across the sample to detect photon density pixel-by-pixel. A computer records the photon density (i.e., fluorescence intensity) along with the location of the XY scanner and Z position to generate a CLSM image. The difference between LSM and CLSM is the use of the pinhole in CLSM that enables axial (Z-plane) selection. Comprehensive reviews of CLSM are available elsewhere.150152

Multiphoton FLIM uses MPE, generally two-photon (2P) or three-photon (3P) excitation, which relies on high photon density to achieve nonlinear excitation of fluorescence. This high density is achieved by lower energy, higher wavelength photons in the near-infrared region. In 2P excitation, two photons of half the energy spontaneously come together to excite the molecule to a higher electronic energy level, which then follows its regular radiative decay (fluorescence) route to relax back to the ground state. Multiphoton FLIM is widely used for tissue imaging because near-infrared wavelengths achieve deeper penetration depths in tissues compared to the visible wavelengths that are commonly used in single photon excitation. This is due to reduced scattering and absorption in tissues within the near-infrared wavelength window. The nonlinear excitation scheme of MPE limits the fluorescence excitation to a small focal volume comparable to the confocal detection volume, but without a pinhole. This allows MP-LSM detectors to be placed in the transmission mode (or nondescanned mode), instead of descanning through the scanning optics. This nondescanned geometry enables higher detection efficiencies. Multiphoton systems use tunable, mode-locked lasers that provide ultrashort, high intensity pulses. A popular source is titanium–sapphire crystal lasers with tunability between 680 and 1100 nm. Most MP-LSM systems include pulsed sources to achieve high photon density, so additional FLIM capabilities only require the timing electronics to estimate photon arrival times. Therefore, many MP-LSM systems likely include FLIM, unlike CLSM systems that conventionally use a continuous-wave excitation source. MP-LSM systems often collect 3-D image tomograms over deeper depths than CLSM. Reviews of MP-LSM are available from Refs. 152 to 156.

Simultaneous excitation of multiple fluorophores is advantageous over sequential imaging because it minimizes FLIM acquisition times. Simultaneous FLIM of three endogenous fluorophores in addition to second harmonic generation (SHG) signals have been achieved by multiphoton wavelength mixing.6 Furthermore, one wavelength has been used to excite the two intrinsic fluorophores, NAD(P)H and FAD.157 In addition, two single-photon wavelengths have been temporally interleaved to alternately excite NAD(P)H and FAD.158



Detectors in LSM are often characterized by their sensitivity, reproducibility, quantum efficiency, photon-counting capability, narrow temporal responses, relatively fast transit time, low dark counts, and high SNR. Most LSM detectors are photomultiplier tubes (PMTs) that can be used in a photon-counting mode, which uses discriminators as described above. MCPs, avalanche photodiode, SPAD, hybrid PMTs, and SPAD arrays are also used for FLIM detection, each with merits and challenges. For example, SPAD arrays are capable of 256×256  pixels including a TDC in each pixel, so an entire FLIM image can be acquired with <100  picosecond resolution.159,160 However, SPAD arrays suffer from lower quantum efficiency at 460 nm (<35%) and a low fill factor (additional microlens arrays can help to effectively guide light), which results in lower photon collection efficiency.161 A detailed discussion on current detectors is given by Bruschini et al.160

Cameras with integrated FLIM capabilities have recently gained popularity. Chen et al. demonstrated wide-field FLIM using a frequency-domain CMOS FLIM camera,162 while Mitchell et al. implemented a frequency-domain CMOS FLIM camera in a lightsheet microscopy system.163 Raspe et al. developed single-image fluorescence lifetime imaging microscopy (siFLIM) with a modulated electron-multiplied-charged couple device FLIM camera capable of simultaneously recording phase-shifted images.164


Analysis of FLIM Data

Quantitative analysis of FLIM images provides insights into cell function, structure–function relationships, and spatial heterogeneity that are not apparent with qualitative observations of the images. Analytical tools for fluorescence lifetime estimations, image segmentation, and heterogeneity analysis are introduced in this section. Many popular methods are highlighted and summarized in Table 4, but this is not an exhaustive list of tools. Innovations in these areas are ongoing, and new FLIM analysis tools are frequently adopted from disparate disciplines.

Table 4

Comparison of commonly-used fluorescence lifetime estimation methods.

Lifetime estimation methodsAdvantagesDisadvantagesReferences
Curve fitting methods
Least squares fitting• High accuracy for high SNR data• Poor accuracy for low SNR data165, 166
• Easy to implement• Assumes Gaussian-distributed noise
• Requires pixel binning
Maximum likelihood estimation• High accuracy for low and high SNR data• Assumes Poisson-distributed noise165, 167
• Easy to implement• Requires pixel binning
• Provides flexible bin widths
Global analysis• Applicable to time-domain and frequency-domain data• Long computation times168171
• Incorporates spatial information to improve accuracy
Bayesian analysis• High accuracy for low and high SNR data• Susceptible to error from initial assumption of decay parameters172174
• Adaptable for fit-free analysis
• Long computation times
• Avoids pixel binning
Phasor methods
Phasor analysis• Fit-free• Poor accuracy for low SNR data172, 175, 176
• Intuitive representation of lifetime estimates
• Susceptible to error from instrument response
• Useful for large time-domain or frequency-domain datasets
• Visualizes lifetime heterogeneity
Deconvolution methods
Stretched exponential analysis• Rapid fitting procedure• Poor fitting for decays greater than three components177, 178
• No prior assumption of the decay distribution• Requires prior knowledge of fluorescent species present
• Visualizes lifetime heterogeneity
Lifetime moment analysis• No prior assumption of IRF• Poor accuracy for low SNR data81, 179
• Poor accuracy for exceptionally long lifetimes
• Visualizes lifetime heterogeneity
Transformation (e.g., Fourier, Laplace)• Rapid fitting procedure• Poor accuracy for low SNR data170
• Requires increased sampling of decay curve
Laguerre deconvolution• Fit-free• Poor accuracy compared with other fit-free techniques176, 180
• Rapid fitting procedure
• Precise (low variation) lifetime estimates• Poorly suited for estimating exceptionally short or long lifetimes


Fluorescence Lifetime Estimation


Curve fitting

Fluorescence lifetime measurements often capture multiple fluorophores within each pixel, resulting in multiexponential decays. In TCSPC, each photon is assigned to a time bin within a lifetime histogram.74 Histograms are fit to multiexponential decay functions described in Eq. (1).76 An IRF is measured from a sample with an instantaneous lifetime (e.g., SHG signal from a urea crystal for two-photon microscopy), which accounts for the temporal response of the optical system.74 Curve fitting analysis requires some prior assumptions, including the number of lifetime components, temporal offset of detected signals, and sources of background fluorescence.74 This method is also highly dependent on the number of photons detected per pixel as higher photon counts will improve the accuracy of the fit. The multicomponent exponential decay estimate is then convolved with the IRF and compared with the experimentally measured lifetime decay curve [Fig. 6(a)].74 The chi-squared (χ2) goodness-of-fit test is used to evaluate agreement between the fit and the measured data. Parameters of the model (ai, τi) are iterated to achieve a chi-squared value closest to 1, indicating the best model fit to the experimental data.74

Fig. 6

Examples of fluorescence lifetime estimation methods. (a) Curve fitting analysis determines lifetime decay variables (αi,τi) by fitting an estimated decay function and estimated or measured IRF to experimental data. This process is iterated with the measured data to optimize goodness-of-fit parameters (χ2). (b) Fit-free methods for estimating lifetime parameters of time-domain or frequency-domain lifetime data also exist. The phasor approach is one such technique frequently used for intuitive representation of fit-free lifetime estimation. Here, measured lifetime data can be transformed into phasor space to visualize pixels with similar lifetime values [Eqs. (12)–(15)]. Universal circle shown as blue semicircle. Example phasors for single exponential species (τ1,τ2) and a two-component mixture of τ1 and τ2 illustrate the rule of linear addition. (c) Neural networks can be trained with simulated or experimental FLIM data for fast generation of fluorescence lifetime maps. Adapted under CC BY-4.0 with permission from Ref. 181.


These parameter estimates and fit quality measurements can be determined from analytical approaches, such as least squares fitting, maximum likelihood estimation, and Bayesian analysis.165,168,172 These methods describe the likelihood of detecting specific photon counts within each time bin from the experimental decay, based on statistical assumptions unique to each method. For example, least squares fitting minimizes the squared difference between measured fluorescence and estimated signal and assumes Gaussian noise, whereas maximal likelihood methods assume Poisson-distributed noise.165167 Both approaches provide comparably accurate estimates of fluorescence decay parameters for lifetime histograms with high photon counts, though maximum likelihood analysis performs better for low photon counts.165 Maximum likelihood analysis also allows for varied bin widths.172,182

Global analysis is another approach to estimating fluorescence lifetimes from low SNR images.168170 One implementation of global analysis assumes that all fluorescent species are present within each pixel.171 This assumption improves estimation accuracy for data with low photon counts or high background signal. Lifetime parameters of the frequency-domain or transformed time-domain data are estimated by fitting for phase shift and modulation values at multiple frequencies, as shown in Fig. 4(d).169,171 Alternatively, simultaneous analysis of per-pixel lifetimes assumes fixed lifetime values across all pixels and iterates lifetime parameter (τi,αi) estimates to improve a whole-image goodness-of-fit measure.183185 This approach provides better parameter estimates by conserving the spatial information typically lost from averaging photon counts across pixels.169 Lifetime estimates can also improve with segmentation prior to global fitting, which is important when shorter acquisition times are needed to capture fast dynamics.170,186

Bayesian analysis has also been used to improve lifetime estimations. This method empirically determines both the prior distribution of the fluorescence decay (not limited to Gaussian distributions) and the likelihood function173 to establish the posterior distribution of parameters.172,187 Parameter estimates are iterated to maximize the posterior distribution and provide reliable lifetime estimates.187 In general, this method yields optimal lifetime fits even with high noise and low total photon counts, but careful selection of the prior distribution is critical to ensuring accurate estimates.173,187 Recent developments in Bayesian fluorescence lifetime estimation can bypass fits to the measured data and therefore bypass assumptions of the prior distribution of lifetime parameters, which can bias the lifetime estimates.172,174,188


Phasor analysis

Phasor analysis is a fit-free technique in which the fluorescence decay from each pixel is transformed into a point in two-dimensional (2-D) phasor space. Phasor representation provides a visual distribution of the molecular species in an image by clustering pixels with similar lifetimes [Fig. 6(b)]. Phasor analysis is instantaneous because it does not require an iterative fit procedure, and visualization in phasor space is especially advantageous for large FLIM datasets.175,189 In phasor analysis, phasor distributions corresponding to similar lifetimes (decays) can be selected to locate the corresponding pixels in the image with similar lifetimes, even if they are spatially separated.175

Phasor analysis can be applied to both time-domain and frequency-domain FLIM measurements. If P(i,j) represents a pixel in the FLIM image with coordinates (i,j) and Ii,j(t) is the fluorescence intensity decay at that pixel, the corresponding coordinates in the phasor plot (g,s) for time-domain measurements are given as190

Eq. (12)


Eq. (13)

where ω=2πf and f=1/T is the laser repetition rate. Notably, variations in background signal or the temporal response of the optical system may introduce error into time-domain lifetime measurements transformed into phasor space, which should be considered when performing phasor transformations.176

In the case of frequency-domain measurements, the coordinates are given as

Eq. (14)


Eq. (15)

where Mi,j is the modulation and ϕi,j is the phase shift of the emission signal with respect to the excitation. The phasor coordinates can also be expressed in terms of lifetime and angular laser repetition frequency (ω). In the case of a single exponential decay, the g and s coordinates are given as

Eq. (16)

gi,j(ω)=  11+(ωτ)2,

Eq. (17)

From the equations of phasor coordinates, (g,s), [i.e., Eqs. (16) and (17)], the following can be derived:123

Eq. (18)

From this equation, it can deduced that all single exponential lifetimes will fall on a semicircle of radius 1/2 and center (1/2,0). Since all possible single exponentials fall on this circle, it is referred to as the universal circle [Fig. 6(b)]. A short lifetime having a smaller phase will lie near the point (1,0), which corresponds to τ=0, while a long lifetime will fall near the universal circle coordinates (0,0), which corresponds to τ=.

Phasors follow the rule of linear addition. For example, the phasor location of a mixture of two species falls on a straight line joining the phasor location of the two individual species on the universal circle [Fig. 6(b)].189,191 The position on this line is determined by the relative fractional contributions of each species. Similarly, the phasor distribution of a three-exponential species will fall in the triangle formed by the three individual phasor locations and similarly for higher order exponentials.192,193 Hence, the phasor distribution of a heterogeneous sample will have a position inside the universal circle. This representation also provides a straightforward interpretation of the biological significance of lifetime values compared with other lifetime estimation analyses.102


Deconvolution analysis

Deconvolution methods recover the lifetime decay from the measured fluorescence signal by deconvolving the estimated optical system response. Deconvolution includes variations on the least squares approach discussed above, stretched exponential/lifetime moment analysis, and methods of transformation (e.g., Fourier and Laplace transforms).177,179,194198 Stretched exponential and lifetime moment methods are fundamentally similar because the lifetimes of individual species are estimated from the total measured decay distribution.177,179 Stretched exponential and lifetime moment methods have specific implementations for time-domain or frequency-domain data. Furthermore, Fourier and Laplace techniques transform the measured decay curve to be proportional to the product of the transformed source excitation pulse and system response.196,197 The lifetime and contribution of each species are then recovered from the transformed system response.196,197,199,200 Deconvolution methods avoid assumptions about the instrument response but still share limitations associated with standard model fitting methods.

Laguerre deconvolution is an alternative to model fitting for fluorescence lifetime estimation. Here, the measured fluorescence decay is transformed and represented in the form of Laguerre polynomials, resulting in a series expansion of the decay and convolved IRF.176,180,195,201 This Laguerre transformation produces linearly independent functions that enable expansion of decays with low susceptibility to noise and proportionality between fluorescence intensity and lifetime decay.176,180 The pixelwise linear combination of Laguerre coefficients provides the per-pixel fluorescence decay.180,201 The Laguerre method is less accurate but more precise in lifetime estimates compared with the similarly fit-free phasor method, especially for values at either extreme.176


Machine learning analysis

Machine learning techniques are another alternative to time-intensive curve fitting procedures. A few key algorithms commonly used for this purpose are highlighted here. Simple neural networks can estimate fluorescence lifetimes directly from TCSPC data by learning weights from examples of curve fitted pixels.181 Variations of convolutional neural networks (CNNs) can also rapidly calculate fluorescence lifetimes. Briefly, CNNs downsample recorded images through kernel convolution and window pooling steps, resulting in a low-resolution image on which predictions of pixel class membership are made (i.e., contraction).202 Pixel positions from the initial pooling steps are recalled to assign class predictions to pixels in upsampled images (i.e., expansion).202 This computational structure has been used to analyze hyperspectral fluorescence lifetime images and to dynamically monitor fluorescence lifetimes in vitro and in vivo [Fig. 6(c)].203,204 Neural networks continue to improve, including simultaneous prediction of fluorescence lifetimes and object segmentation masks, which will be discussed below.


Fluorescence Lifetime Heterogeneity Analysis


Pixel-level analysis

Pixel-level analysis of fluorescence lifetimes can inform on subcellular and cell-level heterogeneity within a sample. Lifetime histograms provide a useful quality check of curve fitting from TCSPC pixels, confirm the presence of distinct fluorescence lifetimes, and/or confirm expected changes in lifetime values from an experimental condition or FRET interaction.205 Distributions of pixels within phasor space provide complementary information on the identity of fluorophores in the sample and lifetime changes throughout an experiment.175,205 Pixel-level FLIM analysis has been previously used to quantify lipid membrane integrity and heterogeneity, immune cell heterogeneity, cell development, protein conformation and organization, and other phenomena.206211


Object-level analysis

Object-level analysis provides a biological context for interpreting FLIM images by averaging lifetime values across all pixels within a single object of interest (e.g., cells, organelles, and bacteria).32 This approach quantifies diversity across cells, organelles, and other features in a similar manner to established techniques such as flow cytometry or colony counting. Segmentation is required for object-level analysis. Automatic segmentation can be achieved with computational approaches such as multiresolution community detection and morphological filtering and thresholding.212214 In addition, unsupervised clustering techniques (e.g., K-means clustering) can segment single cells and intracellular compartments for phasor-based lifetime data.215 Open source packages such as FLIMfit and FLIM-FRET analyzer have been developed for multiple functionalities, including automatic segmentation, lifetime decay fitting, and data visualization.216,217

Machine learning is also popular for image segmentation due to its high accuracy and generalization across imaging formats. Neural networks can learn features of pixels within objects to generate segmentation masks.202 Several architectures have been designed to improve segmentation performance, primarily for intensity-based segmentation of cellular compartments. Variations of CNNs have been developed for object segmentation. These variations include UNets, feature pyramid networks, and Mask-RCNN.202 UNet employs the standard CNN framework, described earlier, but maintains symmetrical contraction and expansion branches, and concatenates contracted layers with expansion layers to better preserve the initial structure of the image.218 Feature pyramid networks follow this scheme, but they sum convolutional layers from contraction and expansion to inform image upsampling.219 Similarly, mask-RCNN uses feature pyramid networks to generate feature maps for regions of interest (ROIs) before convolutional steps that classify pixels within an ROI for object masks.220 Modified versions of these algorithms have been used to segment numerous cell types and intracellular features across imaging platforms,202 and a combination of these algorithms has enabled nuclei segmentation across image types [e.g., hematoxylin and eosin (H&E) and immunofluorescence].221 This is not an exhaustive list of machine learning approaches for fluorescence image segmentation, but it provides an overview of some object-level segmentation and classification tools.

Segmentation of autofluorescence images is challenging due to low SNR and poor spatial specificity of the fluorescence signal. CellProfiler has been used to isolate cells, nuclei, and cytoplasms from two-photon NAD(P)H autofluorescence images.222 This approach identifies nuclei within a specified size range by thresholding pixels above the background fluorescence but below cytoplasmic fluorescence values. Whole cell masks are defined by propagating outward from nuclear masks. Nuclear masks are subtracted from cell masks to isolate cytoplasmic areas. Similarly, a collection of ImageJ plugins have been developed for integrated lifetime decay fitting (SLIM curve) and object segmentation (Trainable Weka Segmentation).223 This approach first requires input of a small subset of fluorescence images in which the user annotates the objects of interest. A number of ImageJ-based filters and transformations are applied to the images to extract features specific to annotated objects.224 A suite of machine learning algorithms (Weka) are then applied to annotated inputs and extracted features to classify pixels in unlabeled images.225 CellProfiler has developed a similar annotation-based method, Ilastik, optimized for fluorescent proteins and dyes but applicable to autofluorescence images.226

Histograms of lifetime values can be plotted across numerous objects for population distribution analysis. This approach visualizes heterogeneity within an object class (e.g., cells or mitochondria) under basal conditions or in response to perturbations. Histograms can be fit with population density models to summarize the distribution of objects and to identify distinct subpopulations of objects. Gaussian mixture modeling is a common population density modeling approach in which multiple Gaussian probability density functions are iteratively fit to each frequency histogram [Fig. 7(a)].227,230 Goodness of fit is assessed by the minimum Akaike information criterion.230,231 However, this approach is limited by assumptions about the number of populations within the data and the Gaussian distribution of the data.232,233 Approaches have been developed to circumvent these assumptions, including density-based clustering [Fig. 8(a)].235 Density-based clustering defines subpopulations within data such that the highest density datapoints define the cluster for the nearest remaining points.235 Previous studies have shown that population distribution analysis of autofluorescence lifetimes can classify cell types, drug response, and disease states [Fig. 7(b)].32,227229,236238 In addition, lifetime distributions can identify objects without prior segmentation. Variations in lifetime distributions have identified other molecular features such as tagged neurons in C. elegans and metabolic activities within tumors.239,240 Overall, population distribution analysis provides unique insights into sample heterogeneity.

Fig. 7

Heterogeneity analysis of fluorescence lifetime data. (a) Histograms of lifetimes per object are fit to distribution models to describe subpopulations and variability in the data. The H-index and wH-index are derived from these fits [Eqs. (19) and (20), respectively]. Here, pi is the proportion of each subpopulation, di is the distance between subpopulation median and global median, and σi is the subpopulation standard deviation. (b) Distribution density models fit to cell-level NAD(P)H mean lifetimes can accurately identify distinct breast cancer cell lines (MDA-MB-231 and SKBr3) in mixed cocultures (proportion of mixtures indicated above plots). Errors (Err) in the model predictions for mean (x¯) and proportion (p) of each population are given within each plot. Adapted with permission from Ref. 227. (c) H-index of in vivo FaDu tumor cell NAD(P)H mean lifetime (right) correlates with in vivo treatment response (left). Adapted with permission from Ref. 228. (d) Cell autofluorescence wH-index is similar for (left) in vitro organoids derived from primary PyVmT tumors and (right) in vivo PyVmT tumors with vehicle and combination treatment. Adapted with permission from Ref. 229.


Fig. 8

Spatial analysis of fluorescence lifetime distribution. Spatial statistical analyses can quantify spatial heterogeneity in fluorescence lifetimes. Here, spatial heterogeneity in cell-level autofluorescence lifetimes is used as an example. (a) Density-based clustering defines cell subpopulations that are mapped back onto lifetime images. Relative proximity measurements define spatial distributions within (intrapopulation proximity) and between (interpopulation proximity) cell subpopulations. (b) Multivariate spatial heterogeneity is quantified with spatial autocorrelation and spatial principal components analysis. Adapted under CC BY-4.0 with permission from Ref. 234.


Heterogeneity in lifetime measurements is commonly quantified from coefficients of variation (CV).179,241245 The CV is the standard deviation divided by the mean of a measurement, which enables comparisons of variability between samples. However, the CV does not define whether distinct subpopulations exist within the lifetime data. Alternatively, quantitative metrics of heterogeneity can be derived from population density models, so the behavior of subpopulations can be compared between conditions. A heterogeneity index (H-index) was previously defined to quantify cell subpopulations of fluorescence lifetimes using population density models of cells in vivo in head and neck cancer.228 This H-index is based on the Shannon diversity index, widely used in ecological studies, and is defined as

Eq. (19)

where di is the distance between the medians of each subpopulation and the median of the overall distribution and pi is the proportion of each subpopulation [Fig. 7(a)]. Here, increases in the H-index reflect increases in the number of subpopulations within a sample, increases in the separation between subpopulations, and equality of population proportions. This H-index continues to be adapted for different applications. Specifically, a weighted heterogeneity index (wH-index) was developed to assess metabolic heterogeneity across in vitro and in vivo breast cancer models.229 The wH-index includes subpopulation standard deviations and is defined as

Eq. (20)

where σi is the standard deviation of each subpopulation and p and d denote their proportions and distances from the overall median of the distribution, respectively [Fig. 7(a)].

These heterogeneity metrics have provided valuable insight into diversity within biological systems. The H-index of cell autofluorescence in vivo in tumors shows more homogeneous activity across cells (lower H-index) for treatments that significantly reduce tumor volume (combination treatment), and more heterogeneous activity across cells (higher H-index) for treatments that do not change tumor volume (cetuximab or cisplatin alone) [Fig. 7(c)]. In addition, wH-index values of cell autofluorescence for control and treated conditions are similar between in vitro tumor organoids and in vivo tumors, indicating similar treatment-induced changes in metabolic heterogeneity in vivo and in vitro for the same tumor model229 [Fig. 7(d)]. Collectively, these studies show that quantitative metrics of fluorescence lifetime heterogeneity provide powerful tools to study diversity in biological systems.


Analysis of the Spatial Distributions of Fluorescence Lifetimes

False-colored fluorescence lifetime images can be generated from curve fitting parameters (e.g., τ1, τ2, α1, and α2) or phasor values (e.g., g and s).175,205 These images are used for qualitative assessments of molecular distribution in biological samples. For example, FLIM images can map lifetime sensors of intracellular molecules (e.g., magnesium, calcium, chromatin, myoglobin, and antigens), pH, oxygen, or temperature.246248 FLIM images of endogenous fluorophores provide qualitative information on the distribution of subcellular and cellular metabolism, biogenesis, and structure.51,80,205,209,243,249252 Furthermore, tissue-level lifetime images can distinguish cellular compartments across diverse tissue types (e.g., stroma, endothelium, epithelium, and cancerous tissue).11,208,241,253,254

Quantitative metrics of spatial heterogeneity have also been developed for FLIM. Spatial statistical analyses have quantified cell-level spatial heterogeneity in autofluorescence lifetimes across in vitro and in vivo tumor models. This approach uses density-based clustering to identify populations with distinct lifetimes, map them back to image space, and then extract proximity measurements to assess spatial distributions within a population and between populations [Fig. 8(a)].234,235 Multivariate spatial autocorrelation and spatial principal components analysis can further define patterns based on multiple fluorescence lifetime fit parameters and multiple fluorophores [e.g., NAD(P)H and FAD] [Fig. 8(b)]. Additional quantitative methods have been developed to evaluate spatial variations in intracellular fluorescence, though these have not been translated for lifetime data. For example, QuantEv measures the localization of fluorescently tagged proteins as a function of the global structure of a cell,255 and a similar approach was designed for spatial analysis of GFP-expressing plant Golgi proteins.256 Quantitative methods to assess fluorescence lifetime spatial distributions will be critical to exploiting the wealth of information in FLIM images.


Multiparametric Analysis of Fluorescence Lifetime Data

Fluorescence lifetime images usually have multiple variables per pixel (e.g., curve fit parameters, fluorescence intensity, and phasor values) that can be used in multivariate classification models to identify distinct cell subsets or functions. For example, partial least squares–discriminant analysis of autofluorescence lifetimes has been used to classify cell-cycle state in heterogeneous samples.257 Specifically, this model included NAD(P)H and FAD fluorescence lifetimes and intensities to separate apoptotic, proliferating, and quiescent tumor cells in FLIM images. Other studies used discriminant analysis of fluorescence intensity, lifetime, and morphological parameters to classify cell types (keratinocytes, adipocytes, myoblasts, cardiomyocytes, and stem cells) in response to metabolic perturbations (growth factor and nutrient starvation/supplementation and environmental stimuli).258 Multivariate FLIM analysis can also use more complex models including nonlinear classifiers (e.g., logistic regression and random forests) and CNNs. These models also achieve high accuracy for multigroup classification based on autofluorescence lifetimes, specifically for T-cell subtypes and activation states (e.g., quiescent/activated, CD3/CD4/CD8 coexpression).259,260 These studies illustrate the strength of multivariate classification models based on fluorescence lifetime data, which provide robust separation of cell types and cell function.


Examples of FLIM in Biology and Medicine

Numerous studies have used FLIM to understand molecular features of biological systems and changes due to disease progression or drug treatment. Below are a few examples of autofluorescence FLIM, FLIM of exogenous molecular probes, and FLIM-FRET.


Autofluorescence FLIM Applications


In vivo autofluorescence FLIM

Numerous sources of molecular contrast make FLIM attractive for in vivo imaging. One of the earliest in vivo FLIM studies was performed with intrinsic sources of contrast in human skin. Koenig et al. investigated changes in autofluorescence and SHG that occur with human skin disease in vivo.33 SHG is a frequency doubling process with instantaneous lifetime (for a review of SHG, see Ref. 261). Collagen fibers create strong SHG signals in skin. Koenig et al. showed that FLIM resolves changes in the autofluorescence lifetime of skin cells within different layers of the skin. Furthermore, FLIM detected cells that had become diseased due to melanoma or fungal infections.33 This study pioneered the use of FLIM for clinical applications.

Similarly, the first in vivo FLIM studies in animal models also focused on autofluorescence. Skala et al. identified differences in the autofluorescence lifetime between normal, low-grade, and high-grade precancerous epithelia in the hamster cheek pouch in vivo.241,262 Later, in vivo FLIM was used to predict treatment response in mouse tumor models. Specifically, NAD(P)H lifetime changes were found to directly correlate to standard tumor response measurements (i.e., tumor volume).228 Importantly, FLIM detected treatment-induced changes in tumors in vivo only 2 days post-treatment, which is earlier than detectable changes in tumor volume [6 days post-treatment, Fig. 7(c)]. A recent study also demonstrated that in vivo FLIM can measure the efficacy of chemotherapy agents in a mouse model of colorectal cancer.263 Furthermore, autofluorescence FLIM can capture metabolic features of specific cell types in vivo. Work by Szulczewski et al. indicated that macrophages have a fluorescence lifetime that differs from mammary tumor cells such that macrophages can be identified and monitored in vivo without labels264 [Fig. 9(a)]. Other in vivo applications of autofluorescence FLIM focus on metabolism in the mouse brain. For example, NAD(P)H lifetimes reveal metabolic preferences in the brain using a well-defined set of inhibitors that target-specific metabolic reactions7,266 [Fig. 9(b)].

Fig. 9

Autofluorescence FLIM applications. (a) NAD(P)H FLIM of a mammary mouse tumor (heatmap) overlaid on an SHG image of collagen (grayscale). Scale bar=100  μm. Adapted with permission from Ref. 264 (b) Mean NAD(P)H lifetimes in solution and in the rat cortex in vivo after metabolic inhibition. [2DG, 2-deoxy-d-glucose; IAA, iodoacetic acid; KCNm, potassium cyanide; FCCP, carbonyl cyanide-4-(trifluoromethoxy)phenylhydrazone; BMI, bicuculline methiodide; ETC, electron transport chain.] * indicates significantly different from baseline in vivo measurement; Error bars indicate standard error across all pixels over all measurements. Reproduced with permission from Ref. 7 (c) Optical redox ratio [NAD(P)H/FAD; first row], NAD(P)H τm (second row), and FAD τm (third row) images of organoids generated from primary human breast tumors obtained from resection surgeries. TNBC, triple negative breast cancer; ER, estrogen receptor. Scale bar=100  μm. Adapted with permission from Ref. 265.


In addition, autofluorescence FLIM has been performed in numerous non-mammalian in vivo models to study organs and whole-body processes that are not easily visualized in mammals. For example, the metabolic gradient along the germline of C. elegans was visualized with autofluorescence FLIM,6 which provided new insights into metabolic changes with germline differentiation. FLIM has also been performed in plants such as Arabidopsis, where FLIM estimated vacuolar pH inside intact plant cells with the lifetime of anthocyanin.267


Three-dimensional in vitro autofluorescence FLIM

3-D in vitro cultures, including organoids and cell constructs within microdevices, have also been assessed with FLIM. Optical sectioning techniques such as CLSM and MP-LSM are especially attractive for FLIM of 3-D cultures due to their high spatial resolution and volumetric imaging capabilities. Numerous cancer studies have focused on predicting in vivo drug response using primary tumor organoids. These organoids retain all of the cells of the original tumor in a 3-D matrix so that in vivo cell–cell interactions and relevant gradients of oxygen, nutrients, and drugs are preserved.268271 For example, MP-LSM FLIM indicates that autofluorescence lifetimes in primary tumor organoids can predict in vivo response in mouse models across a range of treatment conditions in breast265 and head and neck cancer.272 Furthermore, FLIM can detect differences in the metabolism of primary patient-derived tumor organoids based upon their surface marker expression [Fig. 9(c)]. FLIM has been used to investigate treatment response in patient-derived tumor organoids across multiple cancer types including breast,265 pancreatic,273 and colorectal cancer.274 In addition, FLIM of colorectal cancer organoids was used to inform a patient treatment regimen.274 Organoids provide important 3-D architecture for in vitro studies, but microdevices improve the relevance of 3-D cultures by mimicking in vivo structures. Specifically, FLIM monitored changes in the metabolism of ductal carcinoma in situ cells during invasion in a lumen microdevice model. FLIM captured changes in metabolism based on the position of a cell within the lumen or invading branch.275

Tissues ex vivo have also been imaged to determine whether FLIM can guide surgical resection of tumors.276 First, Lukina et al. compared NAD(P)H FLIM of in vivo and ex vivo samples using mouse models of colorectal cancer, lung carcinoma, and melanoma to determine optimal tissue maintenance protocols to preserve in vivo signals within ex vivo samples. Then, Lukina et al. used these protocols to perform NAD(P)H FLIM in postoperative samples obtained from colorectal cancer patients and found significant differences in NAD(P)H lifetimes between normal and malignant specimens.


Autofluorescence FLIM in two-dimensional samples

Autofluorescence FLIM of 2-D in vitro cultures provides simple and repeatable systems to test perturbations of autofluorescence lifetime properties. For example, Walsh et al. showed that NAD(P)H lifetimes can detect metabolic differences due to breast cancer subtype.236 In addition, the fluorescence lifetime of NAD(P)H correlates with the differentiation potential of neural progenitor and stem cells.190 Similarly, changes in the relative fluorescence lifetimes of NAD(P)H and lipid droplet associated granules discriminated differentiated and undifferentiated human embryonic stem cells,277 as well as human induced pluripotent stem cell-derived cardiomyocytes under oxidative stress.278 Further autofluorescence FLIM studies in 2-D culture discriminated activation states in multiple types of immune cells including macrophages209 and T cells.259

Finally, autofluorescence FLIM can resolve subcellular features to study intracellular dynamics, including communication between organelles, subcellular features of whole cell processes such as cell division, and bioenergetic demands of different cell types. Mitochondrial organization is often altered to accommodate cellular bioenergetics and biosynthetic demands. Changes in metabolism are also a hallmark of many diseases including cancer. Therefore, mitochondrial imaging has been especially popular for subcellular FLIM applications. Fluorescent dyes such as TMRE can measure mitochondrial membrane potential, which is closely related to cell health.279 However, mitochondrial dyes can alter cellular respiration,280 and therefore label-free methods are in development. NAD(P)H and FAD fluorescence signals are brightest in the mitochondria, which enables label-free visualization of mitochondria. Pouli et al. showed that FLIM of NAD(P)H and FAD can capture rapid changes in mitochondrial spatial dynamics and metabolism using high-resolution imaging of individual mitochondria within cells.281


FLIM of Exogenous Molecular Probes


Exogenous molecular probes for in vivo applications

Numerous optical probes have been developed for both in vivo and in vitro applications to capitalize on the sensitivity of FLIM to physical conditions, including viscosity,282 temperature,283 acidity,284 and oxygenation.104,247,285 Additional molecular probes have been generated that allow for FLIM-based monitoring of drug delivery.

Mouse models are widely used for in vivo FLIM studies of exogenous molecular probes. Ardeshirpour et al. detected mouse tumors in vivo that express human epidermal growth factor receptor (HER2) with FLIM of a fluorescent anti-HER2 antibody.286 Similarly, FLIM showed that the near-infrared fluorescence dye cypate localizes to mouse tumors in vivo [Fig. 10(a)]. FLIM of two fluorophores, cypate and bacteriochlorophyll, can identify the unique distribution of each fluorophore in vivo,239,287,290 which is difficult with intensity-based imaging alone. FLIM has also evaluated renal function in mice using the fluorescence lifetime reporter LS-288, which has a distinct lifetime when free in solution vs. bound to proteins. This approach provides contrast between the protein-rich viscera and the mostly protein-free bladder in mice in vivo [Fig. 10(b)].288 Furthermore, pH-sensitive fluorescence lifetime probes that provide a nonterminal method to quickly determine the acidity of a region in vivo have been developed.291 Overall, FLIM in conjunction with the development of these sophisticated probes is promising in cancer detection and other in vivo applications.

Fig. 10

Molecular probes for FLIM and FLIM-FRET. (a) Near-infrared fluorescence lifetime image of Cyp-GRD distribution (heatmap) in an A549-tumor-bearing mouse at 24-h postinjection. Adapted with permission from Ref. 287 (b) Fluorescence lifetime (heatmap) of mouse abdomen acquired 90 min after intravenous injection of LS-288. The low fluorescence lifetime region in the center of the abdomen is the filled urinary bladder. Adapted with permission from Ref. 288. (c) FLIM maps of the weighted mean fluorescence lifetime of T2AMPKAR-NES, a sensor for AMPK activation, in HEK293 spheroids. The blue end of the colormap indicates increased AMPK activation. Scale bar=100  μm. Adapted with permission from Ref. 289.


In non-mammalian in vivo models, fluorescence lifetime probes that change with both temperature292,293 and concentration of ions have been developed. For example, Zhang et al. generated a phosphorescent lifetime probe that is temperature dependent and demonstrated this temperature dependence in vivo in a zebrafish model.283 Another example of a non-mammalian application of fluorescence lifetime probes includes imaging chloride ion concentrations in cockroach salivary glands done by Hille et al.294


In vitro molecular probe FLIM

Many fluorescence lifetime probes exist for in vitro applications to measure whole cell changes and localize molecular trafficking within a cell. For example, a fluorescence lifetime probe was developed to track the location and use of Zn2+ within a cell. These probes can be localized to understand ion use within individual organelles. Other fluorescence lifetime probes have been developed to detect intracellular prodrug trafficking,295 as well as pH296 and oxygenation changes. Oxygen sensing via phosphorescent lifetime imaging has become a well-established method to monitor intracellular oxygen tension. Furthermore, simultaneous measurement of NAD(P)H FLIM and oxygen sensing by phosphorescence lifetime imaging of Ruthenium tris-(2,2′-bipyridyl) has also been demonstrated in 2-D cell cultures.297


FLIM-FRET Applications


FLIM-FRET for in vivo applications

Finally, FLIM can be used to better capture extracellular and subcellular interactions on the nanoscale both in vivo and in vitro via FLIM-FRET. FLIM-FRET interactions can be used to measure protein activity, gene regulation, and subcellular dynamics. For example, an activatable FRET probe has been developed with a donor–acceptor pair that can be cleaved by matrix metalloproteinases (MMP). This probe was used in a mouse model of breast cancer to monitor MMP activity.298 FLIM-FRET has also identified patterns in RhoA activity in vivo using a GFP-RFP Raichu-RhoA reporter. These studies found that active RhoA, which is associated with cellular cytoskeleton organization, has subcellular localization to the leading edge of invasion in a pancreatic cancer mouse model.299

FLIM-FRET probes have also been used in non-mammalian in vivo models including plants and zebrafish. In Arabidopsis roots, FLIM-FRET probes have been developed to investigate the role of transcription factors that regulate plant cell fates.300 In zebrafish, FLIM imaged a time-course of apoptosis after radiation treatment in 3-D over the entire zebrafish body using a FRET sensor,301 which provided an important whole-body context for the apoptotic process. These are just a few of the many non-mammalian in vivo applications of FLIM-FRET.


FLIM-FRET for in vitro applications

Subcellular dynamics can also be monitored with FLIM-FRET in vitro. T2AMPKAR-NES is a FLIM-FRET sensor for AMPK activation that was shown to detect spatial changes in the activity of AMPK between the inner and outer layers of human embryonic kidney spheroids in a 3-D culture [Fig. 10(c)].289 FLIM-FRET has also been used in a 2-D culture of mouse pituitary cells to detect dimerization between the transcription factor CAATT and the enhancer binding protein alpha. This dimerization corresponds with increased gene expression and adipogenesis.302 Finally, autofluorescence FLIM-FRET can detect molecular interactions within live cells as well, specifically between NAD(P)H and tryptophan. In these studies, tryptophan is the FRET donor and NAD(P)H is the FRET acceptor [i.e., NAD(P)H quenches tryptophan fluorescence]. These studies found that doxorubicin increases the abundance of FRET interactions between tryptophan and NAD(P)H.303 New developments in super-resolution FLIM can localize molecular features within smaller structures and is growing in popularity along with other super-resolution techniques.133,304


Challenges and Solutions in FLIM Applications

FLIM provides insight into molecular features of living systems, yet challenges remain. FLIM instrumentation is more costly than intensity-based measurements. In addition, FLIM analysis requires specific expertise, and the computational cost of FLIM analysis is often higher than intensity-based imaging. Furthermore, FLIM acquisition is generally slower than intensity imaging because more photons are needed to accurately estimate a fluorescence lifetime for each pixel. This can be limiting for biological processes that occur rapidly.305 Additionally, fluorescence lifetimes are affected by numerous factors (e.g., molecular interactions/binding activity, environmental factors such as pH, temperature, and viscosity), which can complicate the interpretation of fluorescence lifetime changes in biological systems. Like all light microscopy, FLIM suffers from scattering that limits SNR and resolution at deep imaging depths and/or highly scattering samples. The effects of scattering can be improved with MP-LSM, but imaging depths are still limited to <2  mm in most tissues.106 For in vivo applications, motion artifacts from animal breathing and heartbeat require specific sample preparation and/or image gating to maintain quality during FLIM acquisition,306 also in a similar fashion to other light microscopy techniques.

Endogenous fluorophores have quantum yields that are orders of magnitude lower than traditional dyes,76 which presents challenges for autofluorescence FLIM. Furthermore, disentangling the contributions from multiple endogenous fluorophores can be difficult when lifetime values overlap, such as NADH and NADPH or FAD and FMN.49,307310 Numerous drugs also naturally fluoresce, and these properties must be known when measuring fluorescence lifetime changes due to drug treatment. For example, Mohammed et al. could not separate the overlapping lifetimes of NAD(P)H and zinc-oxide nanoparticles, so their contribution was combined into one lifetime component.311 Other experimental parameters including cell density, media conditions, 2-D versus 3-D culture, and coatings on culture dishes can have significant effects on autofluorescence lifetime values.51,205,237 To minimize these challenges, care should be taken to keep parameters consistent across samples and a daily control sample should be imaged. Furthermore, drugs that can directly modulate autofluorescent molecules, such as duroquinone, which has been shown to modulate the ratio of NAD+ to NADH, might aid in interpretation.312

FLIM of molecular probes is often challenged by the nonspecific binding of the probe. Molecular probes in vivo also suffer from high background due to autofluorescence. FLIM probes with long lifetimes (>8  ns) can remove autofluorescence background, but acquisition must be optimized for these longer lifetimes. FLIM-FRET probes must also be designed to accurately measure protein interactions without interfering with the interaction, which is especially difficult when the labeled molecules are overexpressed and likely disrupt the normal state of the cells.313 Therefore, careful management of the experimental conditions and appropriate controls are needed for robust FLIM studies. Overall, FLIM is an enabling technique that requires specific training for reproducible results and appropriate data interpretation.



FLIM is a widely used tool for biomedical imaging and has advanced the field of microscopy in the past few decades. In this review, we discussed FLIM as a technique to measure biophysical changes at the molecular scale. FLIM coupled with fluorescence lifetime probes can quantify chemical and physical changes to molecules including changes in temperature, viscosity, pH, and others. Unlike intensity-based measurements, lifetime is self-referenced and independent of the absolute number of photons. Therefore, FLIM is not corrupted by variations in intensity between pixels.

FLIM instrumentation can be deployed in either the time-domain or frequency-domain, which generates either a photon timing histogram or a phase-frequency plot to measure the exponential decay rate of the fluorescence. Analytical equations describe the decay rate and fluorescence lifetime. FLIM images can be acquired either in a WFI scheme using a camera-based detection or by raster-scanning a focused point of excitation across a sample using a single channel detector. Instrumentation schemes are flexible and can be optimized for the desired field of view, spatial–temporal resolution, imaging speed, and other considerations. Multiphoton FLIM provides a unique tool for 3-D optical sectioning and deeper imaging depth into tissues, which is especially advantageous for in vivo and in situ imaging. Current fast-FLIM systems use electronics with short dead times to increase frame rates for medical applications in surgery. Algorithms for FLIM analysis are under rapid development to improve image segmentation, quantify multidimensional heterogeneity, and perform multiparametric analysis. These computational tools unravel spatial and molecular features of cellular physiology that are not apparent with qualitative observation of FLIM images.

Numerous biomedical applications were discussed including autofluorescence FLIM as a label-free method to monitor metabolism and protein–enzyme interactions with the endogenous fluorophores NAD(P)H and FAD. Autofluorescence FLIM has provided insight into metabolism in cancer, stem cells, immune cells, and the brain across diverse sample types including 3-D organoids, microfluidic physiological systems, in vivo mouse models, and human skin. FLIM-FRET sensors have also quantified molecular interactions related to cellular signaling, cellular proliferation, and cytokinesis. In the future, FLIM technologies, analysis, and applications will continue to develop toward advancements in biological research and clinical assessments.


The authors have no relevant financial interests in this manuscript and no potential conflicts of interest.


We would like to thank Dr. Alba Alfonso-Garcia, Dr. Suman Ranjit, and Dr. Jenu Varghese Chacko for their valuable inputs and suggestions. The authors are supported by grants from the NSF (CBET-1642287), Stand Up to Cancer (SU2C-AACR-IG-08-16, SU2C-AACR-PS-18), the NIH (R01 CA185747, R01 CA205101, R01 CA211082, R21 CA224280, U01 TR002383, R37 CA226526), and the University of Wisconsin Carbone Cancer Center (Support Grant P30 CA014520 and the UWCCC Pancreatic Cancer Taskforce).



G. G. Stokes, “On the change of refrangibility of light,” Philos. Trans. R. Soc. Lond., 142 463 –562 (1852). Google Scholar


Perspectives on Fluorescence: A Tribute to Gregorio Weber, Springer International Publishing, Switzerland (2016). Google Scholar


M. Chalfie et al., “Green fluorescent protein as a marker for gene expression,” Science, 263 (5148), 802 –805 (1994). SCIEAS 0036-8075 Google Scholar


S. Inouye and F. I. Tsuji, “Aequorea green fluorescent protein: expression of the gene and fluorescence characteristics of the recombinant protein,” FEBS Lett., 341 (2), 277 –280 (1994). FEBLAL 0014-5793 Google Scholar


C. Stringari et al., “In vivo single-cell detection of metabolic oscillations in stem cells,” Cell Rep., 10 (1), 1 –7 (2015). Google Scholar


C. Stringari et al., “Multicolor two-photon imaging of endogenous fluorophores in living tissues by wavelength mixing,” Sci. Rep., 7 (1), 3792 (2017). SRCEC3 2045-2322 Google Scholar


M. A. Yaseen et al., “Fluorescence lifetime microscopy of NADH distinguishes alterations in cerebral metabolism in vivo,” Biomed. Opt. Express, 8 (5), 2368 –2385 (2017). BOEICL 2156-7085 Google Scholar


I. Georgakoudi and K. P. Quinn, “Optical imaging using endogenous contrast to assess metabolic state,” Annu. Rev. Biomed. Eng., 14 351 –367 (2012). ARBEF7 1523-9829 Google Scholar


A. A. Heikal, “Intracellular coenzymes as natural biomarkers for metabolic activities and mitochondrial anomalies,” Biomark. Med., 4 (2), 241 –263 (2010). Google Scholar


Y. H. Sun et al., “Fluorescence lifetime imaging microscopy for brain tumor image-guided surgery,” J. Biomed. Opt., 15 (5), 056022 (2010). JBOPFO 1083-3668 Google Scholar


D. Gorpas et al., “Autofluorescence lifetime augmented reality as a means for real-time robotic surgery guidance in human patients,” Sci. Rep., 9 (1), 1187 (2019). SRCEC3 2045-2322 Google Scholar


P. M. Schaefer et al., “NADH autofluorescence-a marker on its way to boost bioenergetic research,” Cytom. Part J. Int. Soc. Anal. Cytol., 95 (1), 34 –46 (2019). Google Scholar


K. Okabe et al., “Intracellular temperature mapping with a fluorescent polymeric thermometer and fluorescence lifetime imaging microscopy,” Nat. Commun., 3 (1), 705 (2012). NCAOBW 2041-1723 Google Scholar


S. Ogikubo et al., “Intracellular pH sensing using autofluorescence lifetime microscopy,” J. Phys. Chem. B, 115 (34), 10385 –10390 (2011). JPCBFK 1520-6106 Google Scholar


K. Suhling et al., “Fluorescence lifetime imaging for viscosity and diffusion measurements,” Proc. SPIE, 10882 1088210 (2019). PSISDG 0277-786X Google Scholar


L. Loura, “FRET in membrane biophysics: an overview,” Front. Physiol., 2 82 (2011). FROPBK 0301-536X Google Scholar


D. W. Piston and G.-J. Kremers, “Fluorescent protein FRET: the good, the bad and the ugly,” Trends Biochem. Sci., 32 (9), 407 –414 (2007). TBSCDB 0968-0004 Google Scholar


B. Bajar et al., “A Guide to fluorescent protein FRET pairs,” Sensors, 16 (9), 1488 (2016). SNSRES 0746-9462 Google Scholar


W. R. Algar et al., “FRET as a biomolecular research tool: understanding its potential while avoiding pitfalls,” Nat. Methods, 16 (9), 815 –829 (2019). 1548-7091 Google Scholar


T. K. Christopoulos, E. P. Diamandis, “14 – Fluorescence immunoassays,” Immunoassay, 309 –335 Academic Press, San Diego (1996). Google Scholar


K. P. Wall, R. Dillon and M. K. Knowles, “Fluorescence quantum yield measurements of fluorescent proteins: a laboratory experiment for a biochemistry or molecular biophysics laboratory course,” Biochem. Mol. Biol. Educ., 43 (1), 52 –59 (2015). 1470-8175 Google Scholar


R. F. Chen, “Fluorescence quantum yields of tryptophan and tyrosine,” Anal. Lett., 1 (1), 35 –42 (1967). ANALBP 0003-2719 Google Scholar


A. U. Rehman et al., “Fluorescence quenching of free and bound NADH in HeLa cells determined by hyperspectral imaging and unmixing of cell autofluorescence,” Biomed. Opt. Express, 8 (3), 1488 –1498 (2017). BOEICL 2156-7085 Google Scholar


T. G. Scott et al., “Synthetic spectroscopic models related to coenzymes and base pairs. V. Emission properties of NADH. Studies of fluorescence lifetimes and quantum efficiencies of NADH, AcPyADH, [reduced acetylpyridineadenine dinucleotide] and simplified synthetic models,” J. Am. Chem. Soc., 92 (3), 687 –695 (1970). JACSAT 0002-7863 Google Scholar


J. R. Lakowicz et al., “Fluorescence lifetime imaging,” Anal. Biochem., 202 (2), 316 –330 (1992). ANBCA2 0003-2697 Google Scholar


N. Ramanujam, “Fluorescence spectroscopy of neoplastic and non-neoplastic tissues,” Neoplasia, 2 (1–2), 89 –117 (2000). Google Scholar


Q. Yu and A. A. Heikal, “Two-photon autofluorescence dynamics imaging reveals sensitivity of intracellular NADH concentration and conformation to cell physiology at the single-cell level,” J. Photochem. Photobiol. B, 95 (1), 46 –57 (2009). JPPBEG 1011-1344 Google Scholar


J. T. Sharick et al., “Protein-bound NAD(P)H lifetime is sensitive to multiple fates of glucose carbon,” Sci. Rep., 8 (1), 5456 (2018). SRCEC3 2045-2322 Google Scholar


N. Nakashima et al., “Picosecond fluorescence lifetime of the coenzyme of D-amino acid oxidase,” J. Biol. Chem., 255 (11), 5261 –5263 (1980). JBCHA3 0021-9258 Google Scholar


B. Koziol et al., “Riboflavin as a source of autofluorescence in Eisenia fetida coelomocytes,” Photochem. Photobiol., 82 (2), 570 –573 (2006). PHCBAP 0031-8655 Google Scholar


R. Leenders et al., “Flavin dynamics in reduced flavodoxins,” Eur. J. Biochem., 211 (1–2), 37 –45 (1993). EJBCAI 0014-2956 Google Scholar


M. Y. Berezin and S. Achilefu, “Fluorescence lifetime measurements and biological imaging,” Chem. Rev., 110 (5), 2641 –2684 (2010). Google Scholar


K. Koenig and I. Riemann, “High-resolution multiphoton tomography of human skin with subcellular spatial resolution and picosecond time resolution,” J. Biomed. Opt., 8 (3), 432 –439 (2003). JBOPFO 1083-3668 Google Scholar


Y. Dancik et al., “Use of multiphoton tomography and fluorescence lifetime imaging to investigate skin pigmentation in vivo,” J. Biomed. Opt., 18 (2), 026022 (2013). JBOPFO 1083-3668 Google Scholar


R. Richards-Kortum and E. Sevick-Muraca, “Quantitative optical spectroscopy for tissue diagnosis,” Annu. Rev. Phys. Chem., 47 (1), 555 –606 (1996). ARPLAP 0066-426X Google Scholar


C. Bueno et al., “Photophysics and photochemical studies of the vitamin B6 group and related derivatives,” Photochem. Photobiol., 86 (1), 39 –46 (2010). PHCBAP 0031-8655 Google Scholar


A. Chorvatova, D. Chorvat, “Tissue fluorophores and their spectroscopic characteristics,” Fluorescence Lifetime Spectroscopy and Imaging, 47 –84 CRC Press, Boca Raton, Florida (2015). Google Scholar


Y. Liu and J. D. Simon, “Isolation and biophysical studies of natural eumelanins: applications of imaging technologies and ultrafast spectroscopy,” Pigment Cell Res., 16 (6), 606 –618 (2003). PCREEA 0893-5785 Google Scholar


A.-M. Pena et al., “Spectroscopic analysis of keratin endogenous signal for skin multiphoton microscopy,” Opt. Express, 13 (16), 6268 –6274 (2005). OPEXFF 1094-4087 Google Scholar


A. Ehlers et al., “Multiphoton fluorescence lifetime imaging of human hair,” Microsc. Res. Tech., 70 (2), 154 –161 (2007). MRTEEO 1059-910X Google Scholar


L. Brancaleon et al., “Characterization of the photoproducts of protoporphyrin IX bound to human serum albumin and immunoglobulin G,” Biophys. Chem., 109 (3), 351 –360 (2004). BICIAZ 0301-4622 Google Scholar


S. D. Zucker et al., “Localization of bilirubin in phospholipid bilayers by parallax analysis of fluorescence quenching,” J. Lipid Res., 42 (9), 1377 –1388 (2001). JLPRAW 0022-2275 Google Scholar


C. D. Tran and G. S. Beddard, ““Excited state properties of bilirubin and its photoproducts using picosecond flourescence and ciruclarly polarized luminescence spectroscopy,” Biochim. Biophys. Acta BBA – Gen. Subj., 678 (3), 497 –504 (1981). Google Scholar


B. Chance et al., “Oxidation-reduction ratio studies of mitochondria in freeze-trapped samples. NADH and flavoprotein fluorescence signals,” J. Biol. Chem., 254 (11), 4764 –4771 (1979). JBCHA3 0021-9258 Google Scholar


K. Harbig et al., “In vivo measurement of pyridine nucleotide fluorescence from cat brain cortex,” J. Appl. Physiol., 41 (4), 480 –488 (1976). Google Scholar


S. Ji et al., “Two-dimensional analysis of the redox state of the rat cerebral cortex in vivo by NADH fluorescence photography,” Brain Res., 119 (2), 357 –373 (1977). BRREAP 0006-8993 Google Scholar


A. Mayevsky and B. Chance, “Intracellular oxidation-reduction state measured in situ by a multichannel fiber-optic surface fluorometer,” Science, 217 (4559), 537 –540 (1982). SCIEAS 0036-8075 Google Scholar


R. Scholz et al., “Flavin and pyridine nucleotide oxidation-reduction changes in perfused rat liver. I. Anoxia and subcellular localization of fluorescent flavoproteins,” J. Biol. Chem., 244 (9), 2317 –2324 (1969). JBCHA3 0021-9258 Google Scholar


R. D. Spencer and G. Weber, “Measurements of subnanosecond fluorescence lifetimes with a cross-correlation phase fluorometer,” Ann. N. Y. Acad. Sci., 158 (1), 361 –376 (1969). ANYAA9 0077-8923 Google Scholar


S. D. M. Islam et al., “Fluorescence quenching of flavin adenine dinucleotide in aqueous solution by pH dependent isomerisation and photo-induced electron transfer,” Chem. Phys., 295 (2), 137 –149 (2003). CMPHC2 0301-0104 Google Scholar


D. K. Bird et al., “Metabolic mapping of MCF10A human breast cells via multiphoton fluorescence lifetime imaging of the coenzyme NADH,” Cancer Res., 65 (19), 8766 –8773 (2005). CNREA8 0008-5472 Google Scholar


E. A. Jares-Erijman and T. M. Jovin, “FRET imaging,” Nat. Biotechnol., 21 (11), 1387 –1395 (2003). NABIF9 1087-0156 Google Scholar


D. Llères et al., “Quantitative FLIM-FRET microscopy to monitor nanoscale chromatin compaction in vivo reveals structural roles of condensin complexes,” Cell Rep., 18 (7), 1791 –1803 (2017). Google Scholar


S. S. Vogel, C. Thaler and S. V. Koushik, “Fanciful FRET,” Sci. STKE, 2006 (331), re2 (2006). Google Scholar


R. N. Day and M. W. Davidson, “Fluorescent proteins for FRET microscopy: monitoring protein interactions in living cells,” BioEssays, 34 (5), 341 –350 (2012). BIOEEJ 0265-9247 Google Scholar


F. Mérola et al., “Newly engineered cyan fluorescent proteins with enhanced performances for live cell FRET imaging,” Biotechnol. J., 9 (2), 180 –191 (2013). Google Scholar


J. Goedhart et al., “Structure-guided evolution of cyan fluorescent proteins towards a quantum yield of 93%,” Nat. Commun., 3 (1), 751 (2012). NCAOBW 2041-1723 Google Scholar


B. G. Abraham et al., “Fluorescent protein based FRET pairs with improved dynamic range for fluorescence lifetime measurements,” PLoS One, 10 (8), e0134436 (2015). POLNCL 1932-6203 Google Scholar


A. J. Lam et al., “Improving FRET dynamic range with bright green and red fluorescent proteins,” Nat. Methods, 9 (10), 1005 –1012 (2012). 1548-7091 Google Scholar


D. M. Grant et al., “Multiplexed FRET to image multiple signaling events in live cells,” Biophys. J., 95 (10), L69 –L71 (2008). BIOJAU 0006-3495 Google Scholar


H. Murakoshi et al., “A dark green fluorescent protein as an acceptor for measurement of Förster resonance energy transfer,” Sci. Rep., 5 (1), 15334 (2015). SRCEC3 2045-2322 Google Scholar


L. Albertazzi et al., “Quantitative FRET analysis with the E0GFP-mCherry fluorescent protein pair,” Photochem. Photobiol., 85 (1), 287 –297 (2009). PHCBAP 0031-8655 Google Scholar


M. Tramier et al., “Sensitivity of CFP/YFP and GFP/mCherry pairs to donor photobleaching on FRET determination by fluorescence lifetime imaging microscopy in living cells,” Microsc. Res. Tech., 69 (11), 933 –939 (2006). MRTEEO 1059-910X Google Scholar


M. Peter et al., “Multiphoton-FLIM quantification of the EGFP-mRFP1 FRET pair for localization of membrane receptor-kinase interactions,” Biophys. J., 88 (2), 1224 –1237 (2005). BIOJAU 0006-3495 Google Scholar


S. Padilla-Parra et al., “Quantitative comparison of different fluorescent protein couples for fast FRET-FLIM acquisition,” Biophys. J., 97 (8), 2368 –2376 (2009). BIOJAU 0006-3495 Google Scholar


S. Ganesan et al., “A dark yellow fluorescent protein (YFP)-based resonance energy-accepting chromoprotein (REACh) for Forster resonance energy transfer with GFP,” Proc. Natl. Acad. Sci. U. S. A., 103 (11), 4089 –4094 (2006). Google Scholar


H. Murakoshi, S.-J. Lee and R. Yasuda, “Highly sensitive and quantitative FRET–FLIM imaging in single dendritic spines using improved non-radiative YFP,” Brain Cell Biol., 36 (1–4), 31 –42 (2008). Google Scholar


H. Murakoshi and A. C. E. Shibata, “ShadowY: a dark yellow fluorescent protein for FLIM-based FRET measurement,” Sci. Rep., 7 (1), 6791 (2017). SRCEC3 2045-2322 Google Scholar


A. Pietraszewska-Bogiel and T. W. J. Gadella, “FRET microscopy: from principle to routine technology in cell biology,” J. Microsc., 241 (2), 111 –118 (2011). JMICAR 0022-2720 Google Scholar


S. P. Poland et al., “Multifocal multiphoton volumetric imaging approach for high-speed time-resolved Förster resonance energy transfer imaging in vivo,” Opt. Lett., 43 (24), 6057 –6060 (2018). OPLEDP 0146-9592 Google Scholar


E. J. Osterlund et al., Rapid imaging of BCL-2 family interactions in live cells using FLIM-FRETBCL-2 Family Proteins: Methods and Protocols, 305 –335 Springer, New York (2019). Google Scholar


J. Lou et al., “Phasor histone FLIM-FRET microscopy quantifies spatiotemporal rearrangement of chromatin architecture during the DNA damage response,” Proc. Natl. Acad. Sci. U. S. A., 116 (15), 7323 –7332 (2019). Google Scholar


J. T. Sharick et al., “Metabolic heterogeneity in patient tumor-derived organoids by primary site and drug treatment,” Front. Oncol., 10 553 (2020). Google Scholar


W. Becker, The BH TCSPC Handbook, 8th ed.Becker et Hickl, Berlin, Germany (2019). Google Scholar


B. J. Glasgow, “Fluorescence lifetime imaging microscopy reveals quenching of fluorescein within corneal epithelium,” Exp. Eye Res., 147 12 –19 (2016). EXERA6 0014-4835 Google Scholar


J. R. Lakowicz, Principles of Fluorescence Spectroscopy, 3rd ed.Springer Science & Business Media, Boston, Massachusetts (2013). Google Scholar


SPCImage Next Generation FLIM Data Analysis Software, Becker & Hickl GmbH, Berlin, Germany (2019). Google Scholar


M. G. Giacomelli et al., “Rapid imaging of surgical breast excisions using direct temporal sampling two photon fluorescent lifetime imaging,” Biomed. Opt. Express, 6 (11), 4317 (2015). BOEICL 2156-7085 Google Scholar


J. Ryu et al., “Real-time visualization of two-photon fluorescence lifetime imaging microscopy using a wavelength-tunable femtosecond pulsed laser,” Biomed. Opt. Express, 9 (7), 3449 (2018). BOEICL 2156-7085 Google Scholar


A. J. Bower et al., “High-speed imaging of transient metabolic dynamics using two-photon fluorescence lifetime imaging microscopy,” Optica, 5 (10), 1290 –1296 (2018). Google Scholar


W. Becker, “Fluorescence lifetime imaging: techniques and applications,” J. Microsc., 247 (2), 119 –136 (2012). JMICAR 0022-2720 Google Scholar


Advanced Time-Correlated Single Photon Counting Applications, Springer International Publishing, Cham (2015). Google Scholar


C. J. de Grauw and H. C. Gerritsen, “Multiple time-gate module for fluorescence lifetime imaging,” Appl. Spectrosc., 55 (6), 670 –678 (2001). APSPA4 0003-7028 Google Scholar


D. R. Yankelevich et al., “Pulse sampling technique,” Fluorescence Lifetime Spectroscopy and Imaging, 47 –84 CRC Press, Boca Raton, Florida (2015). Google Scholar


E. Gratton et al., “Fluorescence lifetime imaging for the two-photon microscope: time-domain and frequency-domain methods,” J. Biomed. Opt., 8 (3), 381 –390 (2003). JBOPFO 1083-3668 Google Scholar


R. M. Clegg, O. Holub and C. Gohlke, “Fluorescence lifetime-resolved imaging: measuring lifetimes in an image,” Methods Enzymol., 360 509 –542 (2003). MENZAU 0076-6879 Google Scholar


X. Liu et al., “Fast fluorescence lifetime imaging techniques: a review on challenge and development,” J. Innov. Opt. Health Sci., 12 (05), 1930003 (2019). Google Scholar


S. Isbaner et al., “Dead-time correction of fluorescence lifetime measurements and fluorescence lifetime imaging,” Opt. Express, 24 (9), 9429 –9445 (2016). OPEXFF 1094-4087 Google Scholar


P. Peronio et al., “Improving the counting efficiency in time-correlated single photon counting experiments by dead-time optimization,” Rev. Sci. Instrum., 86 (11), 113101 (2015). RSINAK 0034-6748 Google Scholar


A. Cominelli et al., “High-speed and low-distortion solution for time-correlated single photon counting measurements: a theoretical analysis,” Rev. Sci. Instrum., 88 (12), 123701 (2017). RSINAK 0034-6748 Google Scholar


G. Acconcia et al., “Fast and compact time-correlated single photon counting system for high-speed measurement with low distortion,” Proc. SPIE, 11246 1124608 (2020). PSISDG 0277-786X Google Scholar


J. L. Rinnenthal et al., “Parallelized TCSPC for dynamic intravital fluorescence lifetime imaging: quantifying neuronal dysfunction in neuroinflammation,” PLoS One, 8 (4), e60100 (2013). POLNCL 1932-6203 Google Scholar


S. P. Poland et al., “A high speed multifocal multiphoton fluorescence lifetime imaging microscope for live-cell FRET imaging,” Biomed. Opt. Express, 6 (2), 277 –296 (2015). BOEICL 2156-7085 Google Scholar


S. Kumar et al., “Multifocal multiphoton excitation and time correlated single photon counting detection for 3-D fluorescence lifetime imaging,” Opt. Express, 15 (20), 12548 –12561 (2007). OPEXFF 1094-4087 Google Scholar


A. L. Trinh et al., “Fast single-cell biochemistry: theory, open source microscopy and applications,” Methods Appl. Fluoresc., 7 (4), 044001 (2019). Google Scholar


E. P. Buurman et al., “Fluorescence lifetime imaging using a confocal laser scanning microscope,” Scanning, 14 (3), 155 –159 (1992). SCNNDF 0161-0457 Google Scholar


J. Sytsma et al., “Time-gated fluorescence lifetime imaging and microvolume spectroscopy using two-photon excitation,” J. Microsc., 191 (1), 39 –51 (1998). JMICAR 0022-2720 Google Scholar


H. C. Gerritsen et al., “Time domain FLIM: theory, instrumentation, and data analysis,” Lab. Tech. Biochem. Mol. Biol., 33 95 –132 (2009). LTBBDT 0075-7535 Google Scholar


R. G. Bennett, “Instrument to measure fluorescence lifetimes in the millimicrosecond region,” Rev. Sci. Instrum., 31 (12), 1275 –1279 (1960). Google Scholar


T. Ni and L. A. Melton, “Fluorescence lifetime imaging: an approach for fuel equivalence ratio imaging,” Appl. Spectrosc., 45 (6), 938 –943 (1991). APSPA4 0003-7028 Google Scholar


H. C. Gerritsen et al., “Fluorescence lifetime imaging in scanning microscopes: acquisition speed, photon economy and lifetime resolution,” J. Microsc., 206 (3), 218 –224 (2002). JMICAR 0022-2720 Google Scholar


F. Fereidouni et al., “A modified phasor approach for analyzing time-gated fluorescence lifetime images,” J. Microsc., 244 (3), 248 –258 (2011). JMICAR 0022-2720 Google Scholar


M. A. M. J. van Zandvoort et al., “Discrimination of DNA and RNA in cells by a vital fluorescent probe: lifetime imaging of SYTO13 in healthy and apoptotic cells,” Cytometry, 47 (4), 226 –235 (2002). CYTODQ 0196-4763 Google Scholar


H. C. Gerritsen et al., “Fluorescence lifetime imaging of oxygen in living cells,” J. Fluoresc., 7 (1), 11 –15 (1997). JOFLEN 1053-0509 Google Scholar


F. Fereidouni, K. Reitsma and H. C. Gerritsen, “High speed multispectral fluorescence lifetime imaging,” Opt. Express, 21 (10), 11769 –11782 (2013). OPEXFF 1094-4087 Google Scholar


L. Marcu, “Fluorescence lifetime techniques in medical applications,” Ann. Biomed. Eng., 40 (2), 304 –331 (2012). ABMECF 0090-6964 Google Scholar


C. Liu et al., “Timing and operating mode design for time-gated fluorescence lifetime imaging microscopy,” Sci. World J., 2013 (2013). Google Scholar


O. J. Steingraber and I. B. Berlman, “Versatile technique for measuring fluorescence decay times in the nanosecond region,” Rev. Sci. Instrum., 34 (5), 524 –529 (1963). RSINAK 0034-6748 Google Scholar


P. V. Butte et al., “Intraoperative delineation of primary brain tumors using time-resolved fluorescence spectroscopy,” J. Biomed. Opt., 15 (2), 027008 (2010). JBOPFO 1083-3668 Google Scholar


L. Marcu et al., “Detection of rupture-prone atherosclerotic plaques by time-resolved laser-induced fluorescence spectroscopy,” Atherosclerosis, 204 (1), 156 –164 (2009). ATHSBL 0021-9150 Google Scholar


M.-A. Mycek, K. T. Schomacker and N. S. Nishioka, “Colonic polyp differentiation using time-resolved autofluorescence spectroscopy,” Gastrointest. Endosc., 48 (4), 390 –394 (1998). Google Scholar


T. J. Pfefer et al., “Temporally and spectrally resolved fluorescence spectroscopy for the detection of high grade dysplasia in Barrett’s esophagus,” Lasers Surg. Med., 32 (1), 10 –16 (2003). LSMEDI 0196-8092 Google Scholar


H. Xie et al., “Multispectral scanning time-resolved fluorescence spectroscopy (TRFS) technique for intravascular diagnosis,” Biomed. Opt. Express, 3 (7), 1521 –1533 (2012). BOEICL 2156-7085 Google Scholar


S. Shrestha et al., “High-speed multispectral fluorescence lifetime imaging implementation for in vivo applications,” Opt. Lett., 35 (15), 2558 –2560 (2010). OPLEDP 0146-9592 Google Scholar


D. R. Yankelevich et al., “Design and evaluation of a device for fast multispectral time-resolved fluorescence spectroscopy and imaging,” Rev. Sci. Instrum., 85 (3), (2014). RSINAK 0034-6748 Google Scholar


A. Alfonso-Garcia et al., “Real-time augmented reality for delineation of surgical margins during neurosurgery using autofluorescence lifetime contrast,” J. Biophotonics, 13 e201900108 (2019). Google Scholar


A. Alfonso-Garcia et al., “Fiber-based platform for synchronous imaging of endogenous and exogenous fluorescence of biological tissue,” Opt. Lett., 44 (13), 3350 –3353 (2019). OPLEDP 0146-9592 Google Scholar


J. Park et al., “A dual-modality optical coherence tomography and fluorescence lifetime imaging microscopy system for simultaneous morphological and biochemical tissue characterization,” Biomed. Opt. Express, 1 (1), 186 –200 (2010). BOEICL 2156-7085 Google Scholar


H. Fatakdawala et al., “Fluorescence lifetime imaging combined with conventional intravascular ultrasound for enhanced assessment of atherosclerotic plaques: an ex vivo study in human coronary arteries,” J Cardiovasc. Transl. Res., 8 (4), 253 –263 (2015). Google Scholar


E. Gaviola, ““Ein Fluorometer. Apparat zur Messung von Fluoreszenzabklingungszeiten,” Z. Für Phys., 42 (11), 853 –861 (1927). Google Scholar


B. D. Venetta, “Microscope phase fluorometer for determining the fluorescence lifetimes of fluorochromes,” Rev. Sci. Instrum., 30 (6), 450 –457 (1959). RSINAK 0034-6748 Google Scholar


J. R. Lakowicz, “Frequency-domain lifetime measurements,” Principles of Fluorescence Spectroscopy, 157 –204 Springer US, Boston, Massachusetts (2010). Google Scholar


D. M. Jameson, E. Gratton and R. D. Hall, “The measurement and analysis of heterogeneous emissions by multifrequency phase and modulation fluorometry,” Appl. Spectrosc. Rev., 20 (1), 55 –106 (1984). APSRBB 0570-4928 Google Scholar


J. A. Ross et al., “Time-resolved methods in biophysics. 8. Frequency domain fluorometry: applications to intrinsic protein fluorescence,” Photochem. Photobiol. Sci., 7 (11), 1301 –1312 (2008). PPSHCB 1474-905X Google Scholar


A. Dvornikov, L. Malacrida and E. Gratton, “The DIVER microscope for imaging in scattering media,” Methods Protoc., 2 (2), 53 (2019). Google Scholar


J. R. Lakowicz et al., “Frequency-domain phase-modulation fluorometry: resolution of complex decays of fluorescence intensity and anisotropy,” J. Lumin., 31–32 699 –702 (1984). JLUMA8 0022-2313 Google Scholar


R. A. Colyer, C. Lee and E. Gratton, “A novel fluorescence lifetime imaging system that optimizes photon efficiency,” Microsc. Res. Tech., 71 (3), 201 –213 (2008). MRTEEO 1059-910X Google Scholar


B. A. Feddersen, D. W. Piston and E. Gratton, “Digital parallel acquisition in frequency domain fluorimetry,” Rev. Sci. Instrum., 60 (9), 2929 –2936 (1989). RSINAK 0034-6748 Google Scholar


J. E. Phipps, J. Bec, L. Marcu, “Dual-modality fluorescence lifetime and intravascular ultrasound for label-free intravascular coronary imaging,” Multimodality Imaging: For Intravascular Application, 153 –171 Springer, Singapore (2020). Google Scholar


W. Becker et al., “FLIM and FCS detection in laser-scanning microscopes: increased efficiency by GaAsP hybrid detectors,” Microsc. Res. Tech., 74 (9), 804 –811 (2011). MRTEEO 1059-910X Google Scholar


W. Becker, C. Wetzker and A. Benda, “Ultra-fast HPM detectors improve NAD(P)H FLIM,” Proc. SPIE, 10498 1049806 (2018). PSISDG 0277-786X Google Scholar


M. Gersbach et al., “High frame-rate TCSPC-FLIM using a novel SPAD-based image sensor,” Proc. SPIE, 7780 77801H (2010). PSISDG 0277-786X Google Scholar


M. Castello et al., “A robust and versatile platform for image scanning microscopy enabling super-resolution FLIM,” Nat. Methods, 16 (2), 175 (2019). 1548-7091 Google Scholar


A. S. Tremsin, J. V. Vallerga and O. H. W. Siegmund, “Overview of spatial and timing resolution of event counting detectors with microchannel plates,” Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip., 949 162768 (2020). Google Scholar


L. M. Hirvonen and K. Suhling, “Wide-field TCSPC: methods and applications,” Meas. Sci. Technol., 28 (1), 012003 (2016). MSTCEP 0957-0233 Google Scholar


S. E. D. Webb et al., “A wide-field time-domain fluorescence lifetime imaging microscope with optical sectioning,” Rev. Sci. Instrum., 73 (4), 1898 –1907 (2002). RSINAK 0034-6748 Google Scholar


K. Dowling et al., “High resolution time-domain fluorescence lifetime imaging for biomedical applications,” J. Mod. Opt., 46 (2), 199 –209 (1999). JMOPEW 0950-0340 Google Scholar


A. V. Agronskaia, L. Tertoolen and H. C. Gerritsen, “High frame rate fluorescence lifetime imaging,” J. Phys. Appl. Phys., 36 (14), 1655 –1662 (2003). Google Scholar


T. W. J. Gadella, T. M. Jovin and R. M. Clegg, “Fluorescence lifetime imaging microscopy (FLIM): spatial resolution of microstructures on the nanosecond time scale,” Biophys. Chem., 48 (2), 221 –239 (1993). BICIAZ 0301-4622 Google Scholar


M. J. Cole et al., “Time-domain whole-field fluorescence lifetime imaging with optical sectioning,” J. Microsc., 203 (3), 246 –257 (2001). JMICAR 0022-2720 Google Scholar


T. Hinsdale et al., “Optically sectioned wide-field fluorescence lifetime imaging microscopy enabled by structured illumination,” Biomed. Opt. Express, 8 (3), 1455 –1465 (2017). BOEICL 2156-7085 Google Scholar


D. M. Grant et al., “High speed optically sectioned fluorescence lifetime imaging permits study of live cell signaling events,” Opt. Express, 15 (24), 15656 –15673 (2007). OPEXFF 1094-4087 Google Scholar


D. M. Grant et al., “Optically sectioned fluorescence lifetime imaging using a Nipkow disk microscope and a tunable ultrafast continuum excitation source,” Opt. Lett., 30 (24), 3353 –3355 (2005). OPLEDP 0146-9592 Google Scholar


K. Greger et al., “Three-dimensional fluorescence lifetime imaging with a single plane illumination microscope provides an improved signal to noise ratio,” Opt. Express, 19 (21), 20743 –20750 (2011). OPEXFF 1094-4087 Google Scholar


A. Ulku et al., “Wide-field time-gated SPAD imager for phasor-based FLIM applications,” Methods Appl. Fluoresc., 8 (2), 024002 (2020). Google Scholar


H. Chen and E. Gratton, “A practical implementation of multifrequency widefield frequency-domain fluorescence lifetime imaging microscopy,” Microsc. Res. Tech., 76 (3), 282 –289 (2013). MRTEEO 1059-910X Google Scholar


H. Chen et al., “Widefield multifrequency fluorescence lifetime imaging using a two-tap complementary metal-oxide semiconductor camera with lateral electric field charge modulators,” J. Biophotonics, 12 (5), e201800223 (2019). Google Scholar


A. J. Bowman et al., “Electro-optic imaging enables efficient wide-field fluorescence lifetime microscopy,” Nat. Commun., 10 (1), 4561 (2019). NCAOBW 2041-1723 Google Scholar


K. Suhling et al., “Wide-field time-correlated single photon counting-based fluorescence lifetime imaging microscopy,” Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip., 942 162365 (2019). Google Scholar


C. L. Smith, “Basic confocal microscopy,” Curr. Protoc. Mol. Biol., 81 (1), 14.11.1 –14.11.18 (2008). Google Scholar


J. Jonkman and C. M. Brown, “Any way you slice it: a comparison of confocal microscopy techniques,” J. Biomol. Tech. JBT, 26 (2), 54 –65 (2015). Google Scholar


A. Diaspro et al., “Fluorescence microscopy,” Springer Handbook of Microscopy, 2 –2 Springer International Publishing, Cham (2019). Google Scholar


W. Denk, J. H. Strickler and W. W. Webb, “Two-photon laser scanning fluorescence microscopy,” Science, 248 (4951), 73 –76 (1990). SCIEAS 0036-8075 Google Scholar


D. W. Piston, “Imaging living cells and tissues by two-photon excitation microscopy,” Trends Cell Biol., 9 (2), 66 –69 (1999). TCBIEK 0962-8924 Google Scholar


A. Ustione and D. W. Piston, “A simple introduction to multiphoton microscopy,” J. Microsc., 243 (3), 221 –226 (2011). JMICAR 0022-2720 Google Scholar


B. R. Masters et al., “Multiphoton excitation microscopy, confocal microscopy, and spectroscopy of living cells and tissues; functional metabolic imaging of human skin in vivo,” Methods Enzymol., 307 513 –536 (1999). MENZAU 0076-6879 Google Scholar


R. Cao, H. K. Wallrabe and A. Periasamy, “Multiphoton FLIM imaging of NAD(P)H and FAD with one excitation wavelength,” J. Biomed. Opt., 25 (1), 014510 (2020). JBOPFO 1083-3668 Google Scholar


W. Becker et al., “Metabolic imaging by simultaneous FLIM of NAD(P)H and FAD,” Proc. SPIE, 10882 108820B (2019). PSISDG 0277-786X Google Scholar


P. Luca et al., “A 256 x 256 SPAD array with in-pixel time to amplitude conversion for fluorescence lifetime imaging microscopy,” (2015). Google Scholar


C. Bruschini et al., “Single-photon avalanche diode imagers in biophotonics: review and outlook,” Light Sci. Appl., 8 (1), 1 –28 (2019). Google Scholar


M. Popleteeva et al., “Fast and simple spectral FLIM for biochemical and medical imaging,” Opt. Express, 23 (18), 23511 (2015). OPEXFF 1094-4087 Google Scholar


H. Chen, G. Holst and E. Gratton, “Using of a modulated CMOS camera for fluorescence lifetime microscopy,” Microsc. Res. Tech., 78 (12), 1075 –1081 (2015). MRTEEO 1059-910X Google Scholar


C. A. Mitchell et al., “Functional in vivo imaging using fluorescence lifetime light-sheet microscopy,” Opt. Lett., 42 (7), 1269 –1272 (2017). OPLEDP 0146-9592 Google Scholar


M. Raspe et al., “siFLIM: single-image frequency-domain FLIM provides fast and photon-efficient lifetime data,” Nat. Methods, 13 (6), 501 –504 (2016). 1548-7091 Google Scholar


M. Maus et al., “An experimental comparison of the maximum likelihood estimation and nonlinear least-squares fluorescence lifetime analysis of single molecules,” Anal. Chem., 73 (9), 2078 –2086 (2001). ANCHAM 0003-2700 Google Scholar


M. Straume, S. G. Frasier-Cadoret, M. L. Johnson, “Least-squares analysis of fluorescence data,” Topics in Fluorescence Spectroscopy: Principles, 177 –240 Springer US, Boston, Massachusetts (2002). Google Scholar


Ž. Bajzer et al., “Maximum likelihood method for the analysis of time-resolved fluorescence decay curves,” Eur. Biophys. J., 20 (5), 247 –262 (1991). EBJOE8 0175-7571 Google Scholar


L. Marcu, P. M. W. French and D. S. Elson, Fluorescence Lifetime Spectroscopy and Imaging: Principles and Applications in Biomedical Diagnostics, CRC Press, Boca Raton, Florida (2014). Google Scholar


P. J. Verveer, A. Squire and P. I. H. Bastiaens, “Global analysis of fluorescence lifetime imaging microscopy data,” Biophys. J., 78 (4), 2127 –2137 (2000). BIOJAU 0006-3495 Google Scholar


S. Pelet et al., “A fast global fitting algorithm for fluorescence lifetime imaging microscopy based on image segmentation,” Biophys. J., 87 (4), 2807 –2817 (2004). BIOJAU 0006-3495 Google Scholar


A. Squire, P. J. Verveer and P. I. Bastiaens, “Multiple frequency fluorescence lifetime imaging microscopy,” J. Microsc., 197 (Pt 2), 136 –149 (2000). JMICAR 0022-2720 Google Scholar


M. I. Rowley et al., “Bayesian analysis of fluorescence lifetime imaging data,” Proc. SPIE, 7903 790325 (2011). PSISDG 0277-786X Google Scholar


B. Kaye et al., “Developing and testing a Bayesian analysis of fluorescence lifetime measurements,” PLoS One, 12 (1), e0169337 (2017). POLNCL 1932-6203 Google Scholar


S. Wang et al., “Nonparametric empirical Bayesian framework for fluorescence-lifetime imaging microscopy,” Biomed. Opt. Express, 10 (11), 5497 –5517 (2019). BOEICL 2156-7085 Google Scholar


S. Ranjit et al., “Fit-free analysis of fluorescence lifetime imaging data using the phasor approach,” Nat. Protoc., 13 (9), 1979 –2004 (2018). 1754-2189 Google Scholar


F. Fereidouni et al., “Rapid fluorescence lifetime estimation with modified phasor approach and Laguerre deconvolution: a comparative study,” Methods Appl. Fluoresc., 5 (3), 035003 (2017). Google Scholar


K. C. B. Lee et al., “Application of the stretched exponential function to fluorescence lifetime imaging,” Biophys. J., 81 (3), 1265 –1274 (2001). BIOJAU 0006-3495 Google Scholar


Advances in Intravital Microscopy: From Basic to Clinical Research, Springer Netherlands, Amsterdam (2014). Google Scholar


A. Esposito, H. C. Gerritsen and F. S. Wouters, “Fluorescence lifetime heterogeneity resolution in the frequency domain by lifetime moments analysis,” Biophys. J., 89 (6), 4286 –4299 (2005). BIOJAU 0006-3495 Google Scholar


P. Pande, C. A. Trivedi and J. A. Jo, “Analysis of fluorescence lifetime imaging microscopy (FLIM) data,” Methods Inf. Med., 49 (5), 531 –536 (2010). Google Scholar


G. Wu et al., “Artificial neural network approaches for fluorescence lifetime imaging techniques,” Opt. Lett., 41 (11), 2561 –2564 (2016). OPLEDP 0146-9592 Google Scholar


G. Nishimura and M. Tamura, “Artefacts in the analysis of temporal response functions measured by photon counting,” Phys. Med. Biol., 50 (6), 1327 –1342 (2005). PHMBA7 0031-9155 Google Scholar


J. R. Knutson, J. M. Beechem and L. Brand, “Simultaneous analysis of multiple fluorescence decay curves: a global approach,” Chem. Phys. Lett., 102 (6), 501 –507 (1983). CHPLBC 0009-2614 Google Scholar


J. M. Beechem, “Global analysis of biochemical and biophysical data,” Methods Enzymol., 210 37 –54 (1992). MENZAU 0076-6879 Google Scholar


E. Gratton et al., “Resolution of mixtures of fluorophores using variable-frequency phase and modulation data,” Biophys. J., 46 (4), 479 –486 (1984). BIOJAU 0006-3495 Google Scholar


T. Sanchez et al., “Metabolic imaging via fluorescence lifetime imaging microscopy for egg and embryo assessment,” Fertil. Steril., 111 (2), 212 –218 (2019). FESTAS 0015-0282 Google Scholar


M. I. Rowley et al., “Robust Bayesian fluorescence lifetime estimation, decay model selection and instrument response determination for low-intensity FLIM imaging,” PLoS One, 11 (6), e0158404 (2016). POLNCL 1932-6203 Google Scholar


K. Santra et al., “A Bayesian approach for extracting fluorescence lifetimes from sparse data sets and its significance for imaging experiments,” Photochem. Photobiol., 95 (3), 773 –779 (2019). PHCBAP 0031-8655 Google Scholar


M. A. Digman et al., “The phasor approach to fluorescence lifetime imaging analysis,” Biophys. J., 94 (2), L14 –L16 (2008). BIOJAU 0006-3495 Google Scholar


C. Stringari et al., “Phasor fluorescence lifetime microscopy of free and protein-bound NADH reveals neural stem cell differentiation potential,” PLoS One, 7 (11), e48014 (2012). POLNCL 1932-6203 Google Scholar


G. I. Redford and R. M. Clegg, “Polar plot representation for frequency-domain analysis of fluorescence lifetimes,” J. Fluoresc., 15 (5), 805 –815 (2005). JOFLEN 1053-0509 Google Scholar


R. Datta et al., “Fluorescence lifetime imaging of endogenous biomarker of oxidative stress,” Sci. Rep., 5 9848 (2015). SRCEC3 2045-2322 Google Scholar


S. Ranjit et al., “Multicomponent analysis of phasor plot in a single pixel to calculate changes of metabolic trajectory in biological systems,” J. Phys. Chem. A, 123 (45), 9865 –9873 (2019). JPCAFH 1089-5639 Google Scholar


W. R. Ware, L. J. Doemeny and T. L. Nemzek, “Deconvolution of fluorescence and phosphorescence decay curves. Least-squares method,” J. Phys. Chem., 77 (17), 2038 –2048 (1973). JPCHAX 0022-3654 Google Scholar


J. A. Jo et al., “Fast model-free deconvolution of fluorescence decay for analysis of biological systems,” J. Biomed. Opt., 9 (4), 743 –753 (2004). JBOPFO 1083-3668 Google Scholar


M. Ameloot, “Laplace deconvolution of fluorescence decay surfaces,” Methods Enzymol., 210 279 –304 (1992). MENZAU 0076-6879 Google Scholar


J. C. Andre et al., “Applications of fast Fourier transform to deconvolution in single photon counting,” J. Phys. Chem., 83 (17), 2285 –2294 (1979). JPCHAX 0022-3654 Google Scholar


D. U. Campos-Delgado et al., “Deconvolution of fluorescence lifetime imaging microscopy by a library of exponentials,” Opt. Express, 23 (18), 23748 –23767 (2015). OPEXFF 1094-4087 Google Scholar


R. Cubeddu et al., “Time-resolved fluorescence imaging in biology and medicine,” J. Phys. Appl. Phys., 35 (9), R61 –R76 (2002). Google Scholar


A. Gafni, R. L. Modlin and L. Brand, “Analysis of fluorescence decay curves by means of the Laplace transformation,” Biophys. J., 15 (3), 263 –280 (1975). BIOJAU 0006-3495 Google Scholar


J. A. Jo et al., “Laguerre nonparametric deconvolution technique of time-resolved fluorescence data: application to the prediction of concentrations in a mixture of biochemical components,” Proc. SPIE, 5326 8 –16 (2004). PSISDG 0277-786X Google Scholar


H. Ajmal et al., “Convolutional neural network based image segmentation: a review,” Proc. SPIE, 10649 106490N (2018). PSISDG 0277-786X Google Scholar


J. T. Smith et al., “Fast fit-free analysis of fluorescence lifetime imaging via deep learning,” Proc. Natl. Acad. Sci. U. S. A., 116 (48), 24019 –24030 (2019). Google Scholar


R. Yao et al., “Net-FLICS: fast quantitative wide-field fluorescence lifetime imaging with compressed sensing–a deep learning approach,” Light Sci. Appl., 8 (1), 26 (2019). Google Scholar


J. V. Chacko and K. W. Eliceiri, “Autofluorescence lifetime imaging of cellular metabolism: sensitivity toward cell density, pH, intracellular, and intercellular heterogeneity,” Cytometry A, 95 (1), 56 –69 (2019). Google Scholar


E. Gratton and T. Parasassi, “Fluorescence lifetime distributions in membrane systems,” J. Fluoresc., 5 (1), 51 –57 (1995). JOFLEN 1053-0509 Google Scholar


M. Stöckl et al., “Detection of lipid domains in model and cell membranes by fluorescence lifetime imaging microscopy of fluorescent lipid analogues,” J. Biol. Chem., 283 (45), 30828 –30837 (2008). JBCHA3 0021-9258 Google Scholar


J. J. Rico-Jimenez et al., “Automated detection of superficial macrophages in atherosclerotic plaques using autofluorescence lifetime imaging,” Atherosclerosis, 285 120 –127 (2019). ATHSBL 0021-9150 Google Scholar


A. Alfonso-García et al., “Label-free identification of macrophage phenotype by fluorescence lifetime imaging microscopy.,” J. Biomed. Opt., 21 (4), 046005 (2016). JBOPFO 1083-3668 Google Scholar


K. W. Y. Cho et al., “Label-free assessment of preimplantation embryo quality by the fluorescence lifetime imaging microscopy (FLIM)-phasor approach,” Reprod. Biomed. Online, 38 e34 –e35 (2019). Google Scholar


G. Krishnamoorthy, “Fluorescence lifetime distribution brings out mechanisms involving biomolecules while quantifying population heterogeneity,” Reviews in Fluorescence 2017, 75 –98 Springer, Cham, Switzerland (2018). Google Scholar


D. Hu et al., “Automatic segmentation of fluorescence lifetime microscopy images of cells using multiresolution community detection—a first study,” J. Microsc., 253 (1), 54 –64 (2014). JMICAR 0022-2720 Google Scholar


L. Baiazitova et al., “Three-dimensional fluorescence lifetime imaging in confocal microscopy of living cells,” in 25th Eur. Signal Process. Conf. (EUSIPCO), 439 –443 (2017). Google Scholar


W. Guo et al., “Automated fluorescence lifetime imaging high-content analysis of Förster resonance energy transfer between endogenously labeled kinetochore proteins in live budding yeast cells,” SLAS Technol. Transl. Life Sci. Innov., 24 (3), 308 –320 (2019). Google Scholar


Y. Zhang et al., “Automatic segmentation of intravital fluorescence microscopy images by K-means clustering of FLIM phasors,” Opt. Lett., 44 (16), 3928 –3931 (2019). OPLEDP 0146-9592 Google Scholar


S. C. Warren et al., “Rapid global fitting of large fluorescence lifetime imaging microscopy datasets,” PLoS One, 8 (8), e70687 –e70687 (2013). POLNCL 1932-6203 Google Scholar


J. Kim et al., “FLIM-FRET analyzer: open source software for automation of lifetime-based FRET analysis,” Source Code Biol. Med., 12 (1), 7 (2017). 1751-0473 Google Scholar


O. Ronneberger, P. Fischer and T. Brox, “U-net: convolutional networks for biomedical image segmentation,” Lect. Notes Comput. Sci., 9351 234 –241 (2015). LNCSD9 0302-9743 Google Scholar


T.-Y. Lin et al., “Feature pyramid networks for object detection,” in Proc. IEEE Conf. Comput. Vision Pattern Recognit., 2117 –2125 (2017). Google Scholar


K. He et al., “Mask R-CNN,” in Proc. IEEE Int. Conf. Comput. Vision, 2961 –2969 (2017). Google Scholar


J. C. Caicedo et al., “Nucleus segmentation across imaging experiments: the 2018 Data Science Bowl,” Nat. Methods, 16 (12), 1247 –1253 (2019). 1548-7091 Google Scholar


A. J. Walsh and M. C. Skala, “An automated image processing routine for segmentation of cell cytoplasms in high-resolution autofluorescence images,” Proc. SPIE, 8948 89481M (2014). PSISDG 0277-786X Google Scholar


E. T. Arena et al., “Quantitating the cell: turning images into numbers with ImageJ,” Wiley Interdiscip. Rev. Dev. Biol., 6 (2), e260 (2017). Google Scholar


I. Arganda-Carreras et al., “Trainable Weka segmentation: a machine learning tool for microscopy pixel classification,” Bioinformatics, 33 (15), 2424 –2426 (2017). BOINFP 1367-4803 Google Scholar


M. Hall et al., “The WEKA data mining software: an update,” ACM SIGKDD Explor. Newsl., 11 (1), 10 –18 (2009). Google Scholar


C. Sommer et al., “Ilastik: interactive learning and segmentation toolkit,” in IEEE Int. Symp. Biomed. Imaging: From Nano to Macro, 230 –233 (2011). Google Scholar


A. J. Walsh and M. C. Skala, “Optical metabolic imaging quantifies heterogeneous cell populations,” Biomed. Opt. Express, 6 (2), 559 –573 (2015). BOEICL 2156-7085 Google Scholar


A. T. Shah et al., “In vivo autofluorescence imaging of tumor heterogeneity in response to treatment,” Neoplasia, 17 (12), 862 –870 (2015). Google Scholar


J. T. Sharick et al., “Cellular metabolic heterogeneity in vivo is recapitulated in tumor organoids,” Neoplasia, 21 (6), 615 –626 (2019). Google Scholar


W. Pan, J. Lin and C. T. Le, “Model-based cluster analysis of microarray gene-expression data,” Genome Biol., 3 (2), RESEARCH0009 (2002). GNBLFW 1465-6906 Google Scholar


H. Akaike, “A new look at the statistical model identification,” Selected Papers of Hirotugu Akaike, 215 –222 Springer, New York (1974). Google Scholar


S. A. Haney, “Rapid assessment and visualization of normality in high-content and other cell-level data and its impact on the interpretation of experimental results,” J. Biomol. Screen., 19 (5), 672 –684 (2014). JBISF3 1087-0571 Google Scholar


S. Huang, “Non-genetic heterogeneity of cells in development: more than just noise,” Development, 136 (23), 3853 –3862 (2009). Google Scholar


T. M. Heaster, B. A. Landman and M. C. Skala, “Quantitative spatial analysis of metabolic heterogeneity across in vivo and in vitro tumor models,” Front. Oncol., 9 1144 (2019). FRTOA7 0071-9676 Google Scholar


A. Rodriguez and A. Laio, “Clustering by fast search and find of density peaks,” Science, 344 (6191), 1492 –1496 (2014). SCIEAS 0036-8075 Google Scholar


A. J. Walsh et al., “Optical metabolic imaging identifies glycolytic levels, sub-types and early treatment response in breast cancer,” Cancer Res., 73 (20), 6164 –6174 (2013). CNREA8 0008-5472 Google Scholar


T. M. Cannon, A. T. Shah and M. C. Skala, “Autofluorescence imaging captures heterogeneous drug response differences between 2D and 3D breast cancer cultures,” Biomed. Opt. Express, 8 (3), 1911 –1925 (2017). BOEICL 2156-7085 Google Scholar


T. M. Heaster et al., “Density-based clustering analyses to identify heterogeneous cellular sub-populations,” Proc. SPIE, 10043 100430X (2017). PSISDG 0277-786X Google Scholar


W. Akers et al., “In vivo resolution of multiexponential decays of multiple near-infrared molecular probes by fluorescence lifetime-gated whole-body time-resolved diffuse optical imaging,” Mol. Imaging, 6 (4), 2007 –7290 (2007). Google Scholar


V. K. Ramanujan et al., “Multiphoton fluorescence lifetime contrast in deep tissue imaging: prospects in redox imaging and disease diagnosis,” J. Biomed. Opt., 10 (5), 051407 (2005). JBOPFO 1083-3668 Google Scholar


M. C. Skala et al., “In vivo multiphoton microscopy of NADH and FAD redox states, fluorescence lifetimes, and cellular morphology in precancerous epithelia,” Proc. Natl. Acad. Sci. U. S. A., 104 (49), 19494 –19499 (2007). Google Scholar


E. S. Berge et al., “Fluorescence lifetime imaging of unstained, fixed FaDu tumor slides,” Proc. SPIE, 10890 108901D (2019). PSISDG 0277-786X Google Scholar


C. Dysli, S. Wolf and M. S. Zinkernagel, “Autofluorescence lifetimes in geographic atrophy in patients with age-related macular degeneration,” Invest. Ophthalmol. Vis. Sci., 57 (6), 2479 –2487 (2016). IOVSDA 0146-0404 Google Scholar


K. M. Katika et al., “In vivo time-resolved autofluorescence measurements on human skin,” Proc. SPIE, 6078 60780L (2006). PSISDG 0277-786X Google Scholar


A. Esposito et al., “Unsupervised fluorescence lifetime imaging microscopy for high content and high throughput screening,” Mol. Cell. Proteomics, 6 (8), 1446 –1454 (2007). Google Scholar


A. Sargenti et al., “Fluorescence lifetime imaging of intracellular magnesium content in live cells,” Analyst, 144 (6), 1876 –1880 (2019). ANLYAG 0365-4885 Google Scholar


R. Penjweini et al., “Intracellular oxygen mapping using a myoglobin-mCherry probe with fluorescence lifetime imaging,” J. Biomed. Opt., 23 (10), 107001 (2018). JBOPFO 1083-3668 Google Scholar


A. Pliss et al., “Fluorescence lifetime of fluorescent proteins as an intracellular environment probe sensing the cell cycle progression,” ACS Chem. Biol., 7 (8), 1385 –1392 (2012). Google Scholar


C.-W. Chang and M.-A. Mycek, “Precise fluorophore lifetime mapping in live-cell, multi-photon excitation microscopy,” Opt. Express, 18 (8), 8688 –8696 (2010). OPEXFF 1094-4087 Google Scholar


K. König et al., “Two-photon excited lifetime imaging of autofluorescence in cells during UV A and NIR photostress,” J. Microsc., 183 (3), 197 –204 (1996). JMICAR 0022-2720 Google Scholar


E. Stuntz et al., “Endogenous two-photon excited fluorescence imaging characterizes neuron and astrocyte metabolic responses to manganese toxicity,” Sci. Rep., 7 (1), 1041 (2017). SRCEC3 2045-2322 Google Scholar


T. S. Blacker and M. R. Duchen, “Investigating mitochondrial redox state using NADH and NADPH autofluorescence,” Free Radic. Biol. Med., 100 53 –65 (2016). FRBMEH 0891-5849 Google Scholar


T. Hato et al., “Two-photon intravital fluorescence lifetime imaging of the kidney reveals cell-type specific metabolic signatures,” J. Am. Soc. Nephrol., 28 (8), 2420 –2430 (2017). JASNEU 1046-6673 Google Scholar


J. D. Jones et al., “In vivo multiphoton microscopy detects longitudinal metabolic changes associated with delayed skin wound healing,” Commun. Biol., 1 (1), 198 (2018). Google Scholar


T. Pécot et al., “A quantitative approach for analyzing the spatio-temporal distribution of 3D intracellular events in fluorescence microscopy,” eLife, 7 e32311 (2018). Google Scholar


E. Biot et al., “Strategy and software for the statistical spatial analysis of 3D intracellular distributions,” Plant J., 87 (2), 230 –242 (2016). Google Scholar


T. M. Heaster et al., “Autofluorescence imaging identifies tumor cell-cycle status on a single-cell level,” J. Biophotonics, 11 (1), e201600276 (2018). Google Scholar


Z. Liu et al., “Mapping metabolic changes by noninvasive, multiparametric, high-resolution imaging using endogenous contrast,” Sci. Adv., 4 (3), eaap9302 (2018). STAMCV 1468-6996 Google Scholar


A. Walsh et al., “Label-free method for classification of T cell activation,” (2019). Google Scholar


Z. J. Wang et al., “Classifying T cell activity in autofluorescence intensity images with convolutional neural networks,” J. Biophot., 13 (3), e201960050 (2020). Google Scholar


X. Chen et al., “Second harmonic generation microscopy for quantitative analysis of collagen fibrillar structure,” Nat. Protoc., 7 (4), 654 –669 (2012). 1754-2189 Google Scholar


M. C. Skala et al., “In vivo multiphoton fluorescence lifetime imaging of protein-bound and free nicotinamide adenine dinucleotide in normal and precancerous epithelia,” J. Biomed. Opt., 12 (2), 024014 (2007). JBOPFO 1083-3668 Google Scholar


M. M. Lukina et al., “In vivo metabolic and SHG imaging for monitoring of tumor response to chemotherapy,” Cytometry A, 95 (1), 47 –55 (2019). Google Scholar


J. M. Szulczewski et al., “In vivo visualization of stromal macrophages via label-free FLIM-based metabolite imaging,” Sci. Rep., 6 (1), 25086 (2016). SRCEC3 2045-2322 Google Scholar


A. J. Walsh et al., “Quantitative optical imaging of primary tumor organoid metabolism predicts drug response in breast cancer,” Cancer Res., 74 (18), 5184 –5194 (2014). CNREA8 0008-5472 Google Scholar


M. A. Yaseen et al., “In vivo imaging of cerebral energy metabolism with two-photon fluorescence lifetime microscopy of NADH,” Biomed. Opt. Express, 4 (2), 307 (2013). BOEICL 2156-7085 Google Scholar


A. Chanoca et al., “Using fluorescence lifetime microscopy to study the subcellular localization of anthocyanins,” Plant J., 88 (5), 895 –903 (2016). Google Scholar


C. G. Hubert et al., “A three-dimensional organoid culture system derived from human glioblastomas recapitulates the hypoxic gradients and cancer stem cell heterogeneity of tumors found in vivo,” Cancer Res., 76 (8), 2465 –2477 (2016). CNREA8 0008-5472 Google Scholar


L. M. Langan et al., “Direct measurements of oxygen gradients in spheroid culture system using electron parametric resonance oximetry,” PLoS One, 11 (2), e0149492 (2016). POLNCL 1932-6203 Google Scholar


W. Lim and S. Park, “A microfluidic spheroid culture device with a concentration gradient generator for high-throughput screening of drug efficacy,” Mol. Basel Switz., 23 (12), 3355 (2018). Google Scholar


M. M. Lukina et al., ““Metabolic cofactors NAD(P)H and FAD as potential indicators of cancer cell response to chemotherapy with paclitaxel,” Biochim. Biophys. Acta BBA – Gen. Subj., 1862 (8), 1693 –1700 (2018). Google Scholar


A. T. Shah, T. M. Heaster and M. C. Skala, “Metabolic imaging of head and neck cancer organoids,” PLoS One, 12 (1), e0170415 (2017). POLNCL 1932-6203 Google Scholar


A. J. Walsh et al., “Optical imaging of drug-induced metabolism changes in murine and human pancreatic cancer organoids reveals heterogeneous drug response,” Pancreas, 45 (6), 863 –869 (2016). PANCE4 0885-3177 Google Scholar


C. A. Pasch et al., “Patient-derived cancer organoid cultures to predict sensitivity to chemotherapy and radiation,” Clin. Cancer Res., 25 (17), 5376 –5387 (2019). Google Scholar


J. M. Ayuso et al., “Organotypic microfluidic breast cancer model reveals starvation-induced spatial-temporal metabolic adaptations,” EBioMedicine, 37 144 –157 (2018). Google Scholar


M. M. Lukina et al., “Interrogation of tumor metabolism in tissue samples ex vivo using fluorescence lifetime imaging of NAD(P)H,” Methods Appl. Fluoresc., 8 (1), 014002 (2019). 10.1088/2050-6120/ab4ed8 Google Scholar


C. Stringari et al., “Label-free separation of human embryonic stem cells and their differentiating progenies by phasor fluorescence lifetime microscopy,” J. Biomed. Opt., 17 (4), 046012 (2012). JBOPFO 1083-3668 Google Scholar


R. Datta et al., “Label-free imaging of metabolism and oxidative stress in human induced pluripotent stem cell-derived cardiomyocytes,” Biomed. Opt. Express, 7 (5), 1690 –1701 (2016). BOEICL 2156-7085 Google Scholar


I. A. Okkelman, D. B. Papkovsky and R. I. Dmitriev, “Estimation of the mitochondrial membrane potential using fluorescence lifetime imaging microscopy,” Cytometry A, (2019). Google Scholar


R. C. Scaduto and L. W. Grotyohann, “Measurement of mitochondrial membrane potential using fluorescent rhodamine derivatives,” Biophys. J., 76 (1), 469 –477 (1999). BIOJAU 0006-3495 Google Scholar


D. Pouli et al., “Imaging mitochondrial dynamics in human skin reveals depth-dependent hypoxia and malignant potential for diagnosis,” Sci. Transl. Med., 8 (367), 367ra169 (2016). STMCBQ 1946-6234 Google Scholar


L. Hao et al., “Monitoring mitochondrial viscosity with anticancer phosphorescent Ir(III) complexes via two-photon lifetime imaging,” Chem. Sci., 10 (5), 1285 –1293 (2019). 1478-6524 Google Scholar


H. Zhang et al., “Dual-emissive phosphorescent polymer probe for accurate temperature sensing in living cells and Zebrafish using ratiometric and phosphorescence lifetime imaging microscopy,” ACS Appl. Mater. Interfaces, 10 (21), 17542 –17550 (2018). AAMICK 1944-8244 Google Scholar


A. Almutairi et al., “Biodegradable pH-sensing dendritic nanoprobes for near-infrared fluorescence lifetime and intensity imaging,” J. Am. Chem. Soc., 130 (2), 444 –445 (2008). JACSAT 0002-7863 Google Scholar


A. V. Zhdanov et al., “Imaging of oxygen gradients in giant umbrella cells: an ex vivo PLIM study,” Am. J. Physiol. Cell Physiol., 309 (7), C501 –C509 (2015). 1522-1563 Google Scholar


Y. Ardeshirpour et al., “In vivo fluorescence lifetime imaging monitors binding of specific probes to cancer biomarkers,” PLoS One, 7 (2), e31881 (2012). POLNCL 1932-6203 Google Scholar


S. R. Bloch et al., “Whole-body fluorescence lifetime imaging of a tumor-targeted near-infrared molecular probe in mice,” J. Biomed. Opt., 10 (5), 054003 (2005). JBOPFO 1083-3668 Google Scholar


R. J. Goiffon et al., “Dynamic noninvasive monitoring of renal function in vivo by fluorescence lifetime imaging,” J. Biomed. Opt., 14 (2), 020501 (2009). JBOPFO 1083-3668 Google Scholar


G. Chennell et al., “Imaging of metabolic status in 3D cultures with an improved AMPK FRET biosensor for FLIM,” Sensors, 16 (8), (2016). SNSRES 0746-9462 Google Scholar


S. Bloch et al., “Targeting beta-3 integrin using a linear hexapeptide labeled with a near-infrared fluorescent molecular probe,” Mol. Pharm., 3 (5), 539 –549 (2006). Google Scholar


M. Y. Berezin et al., “Near-infrared fluorescence lifetime pH-sensitive probes,” Biophys. J., 100 (8), 2063 –2072 (2011). BIOJAU 0006-3495 Google Scholar


N. Inada et al., “Temperature imaging using a cationic linear fluorescent polymeric thermometer and fluorescence lifetime imaging microscopy,” Nat. Protoc., 14 (4), 1293 –1321 (2019). 1754-2189 Google Scholar


H. Itoh et al., “Direct organelle thermometry with fluorescence lifetime imaging microscopy in single myotubes,” Chem. Commun., 52 (24), 4458 –4461 (2016). Google Scholar


C. Hille et al., “Two-photon fluorescence lifetime imaging of intracellular chloride in cockroach salivary glands,” Photochem. Photobiol. Sci., 8 (3), 319 –327 (2009). PPSHCB 1474-905X Google Scholar


D. Maji et al., “Cellular trafficking of Sn-2 phosphatidylcholine prodrugs studied with fluorescence lifetime imaging and super-resolution microscopy,” Precis. Nanomed., 1 (2), 128 –145 (2018). Google Scholar


Y. Ma et al., “Phosphorescent soft salt for ratiometric and lifetime imaging of intracellular pH variations,” Chem. Sci., 7 (5), 3338 –3346 (2016). 1478-6524 Google Scholar


S. Kalinina et al., “Correlative NAD(P)H-FLIM and oxygen sensing-PLIM for metabolic mapping,” J. Biophotonics, 9 (8), 800 –811 (2016). Google Scholar


M. Solomon et al., “Detection of enzyme activity in orthotopic murine breast cancer by fluorescence lifetime imaging using a fluorescence resonance energy transfer-based molecular probe,” J. Biomed. Opt., 16 (6), 066019 (2011). JBOPFO 1083-3668 Google Scholar


P. Timpson et al., “Spatial regulation of RhoA activity during pancreatic cancer cell invasion driven by mutant p53,” Cancer Res., 71 (3), 747 –757 (2011). CNREA8 0008-5472 Google Scholar


Y. Long et al., “In vivo FRET–FLIM reveals cell-type-specific protein interactions in Arabidopsis roots,” Nature, 548 (7665), 97 –102 (2017). Google Scholar


N. Andrews et al., “Visualising apoptosis in live zebrafish using fluorescence lifetime imaging with optical projection tomography to map FRET biosensor activity in space and time,” J. Biophotonics, 9 (4), 414 –424 (2016). Google Scholar


M. Elangovan, R. N. Day and A. Periasamy, “Nanosecond fluorescence resonance energy transfer-fluorescence lifetime imaging microscopy to localize the protein interactions in a single living cell,” J. Microsc., 205 (1), 3 –14 (2002). JMICAR 0022-2720 Google Scholar


R. Cao et al., “Single-cell redox states analyzed by fluorescence lifetime metrics and tryptophan FRET interaction with NAD (P) H,” Cytometry A, 95 (1), 110 –121 (2019). Google Scholar


Y. Zhang et al., “Super-resolution multiphoton frequency-domain fluorescence lifetime imaging microscopy by generalized stepwise optical saturation (GSOS),” Proc. SPIE, 10884 108840C (2019). PSISDG 0277-786X Google Scholar


P. Sarder, D. Maji and S. Achilefu, “Molecular probes for fluorescence lifetime imaging,” Bioconjug. Chem., 26 (6), 963 –974 (2015). BCCHES 1043-1802 Google Scholar


S. C. Warren et al., “Removing physiological motion from intravital and clinical functional imaging data,” eLife, 7 (2018). Google Scholar


R. D. Spencer, G. Weber, “Thermodynamics and kinetics of the intramolecular complex in flavin-adenine dinucleotide,” Structure and Function of Oxidation–Reduction Enzymes, 393 –399 Pergamon, Oxford (1972). Google Scholar


G. A. Wagnieres, W. M. Star and B. C. Wilson, “In vivo fluorescence spectroscopy and imaging for oncological applications,” Photochem. Photobiol., 68 (5), 603 –632 (1998). PHCBAP 0031-8655 Google Scholar


T. S. Blacker et al., “Separating NADH and NADPH fluorescence in live cells and tissues using FLIM,” Nat. Commun., 5 (1), 1 –9 (2014). NCAOBW 2041-1723 Google Scholar


P. Wahl and J. C. Auchet, “Time resolved fluorescence of flavin adenine dinucleotide,” FEBS Lett., 44 (1), 67 –70 (1974). FEBLAL 0014-5793 Google Scholar


Y. H. Mohammed et al., “Noninvasive in vivo human multiphoton microscopy: a key method in proving nanoparticulate zinc oxide sunscreen safety,” J. Biomed. Opt., 25 (1), 014509 (2020). JBOPFO 1083-3668 Google Scholar


D. Y. Gui et al., “Environment dictates dependence on mitochondrial complex i for NAD+ and aspartate production and determines cancer cell sensitivity to metformin,” Cell Metab., 24 (5), 716 –727 (2016). 1550-4131 Google Scholar


E. J. Osterlund, Q. Liu, D. W. Andrews, “The use of FLIM-FRET for the detection of mitochondria-associated protein interactions,” Mitochondrial Medicine: Volume I, Probing Mitochondrial Function, 395 –419 Springer, New York (2015). Google Scholar


Rupsa Datta is an assistant scientist in the Skala Lab at Morgridge Institute for Research, Madison. She graduated from the University of California, Irvine, with a PhD in biomedical engineering (BME). Trained in fluorescence lifetime imaging microscopy, her research entails employing optical imaging techniques to study cancer, specifically, to decipher the influence of stroma in development and progression of pancreatic cancer and thereby devise new approaches to cancer therapy.

Tiffany M. Heaster received her BS degree in biological engineering from Mississippi State University and her MS degree in BME from Vanderbilt University. She is currently completing her PhD in BME at the University of Wisconsin–Madison with Dr. Melissa Skala. Her current research involves developing metabolic autofluorescence imaging and analyses for assessing tumor microenvironment heterogeneity.

Joe T. Sharick received his BSE degree in BME from Duke University and his MS degree and PhD in BME from Vanderbilt University. His work focused on the development of optical metabolic imaging to detect lethal drug-resistant subpopulations of cells hidden within a patient’s tumor. He is currently a postdoctoral researcher at Ohio State University. He is interested in engineering new in vitro and in vivo models of cancer for personalized medicine and drug development.

Amani A. Gillette received her BS degree in BME from Michigan Technological University. She is currently completing her PhD in BME at the University of Wisconsin–Madison with Dr. Melissa Skala. Her research focuses on developing metabolic autofluorescence imaging to study the effects of single protein disruptions to metabolic pathways and analysis of mitochondria with label-free imaging.

Melissa C. Skala received her BS degree in physics from Washington State University in 2002, her MS degree in BME from the University of Wisconsin–Madison in 2004, and her PhD in BME from Duke University in 2007. Her postdoctoral training was also in BME at Duke University from 2007 to 2010. She started her lab at Vanderbilt University. She is now an investigator at the Morgridge Institute for Research and an associate professor of BME at the University of Wisconsin–Madison.

© The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Rupsa Datta, Tiffany M. Heaster, Joe T. Sharick, Amani A. Gillette, and Melissa C. Skala "Fluorescence lifetime imaging microscopy: fundamentals and advances in instrumentation, analysis, and applications," Journal of Biomedical Optics 25(7), 071203 (13 May 2020).
Received: 6 January 2020; Accepted: 24 April 2020; Published: 13 May 2020

Back to Top