Two-photon imaging of endogenous fluorescence can provide physiological and metabolic information from intact tissues. However, simultaneous imaging of multiple intrinsic fluorophores, such as nicotinamide adenine dinucleotide(phosphate) (NAD(P)H), flavin adenine dinucleotide (FAD) and retinoids in living systems is generally hampered by sequential multi-wavelength excitation resulting in motion artifacts. Here, we report on efficient and simultaneous multicolor two-photon excitation of endogenous fluorophores with absorption spectra spanning the 750-1040 nm range, using wavelength mixing. By using two synchronized pulse trains at 760 and 1041 nm, an additional equivalent two-photon excitation wavelength at 879 nm is generated, and achieves simultaneous excitation of blue, green and red intrinsic fluorophores. This method permits an efficient simultaneous imaging of the metabolic coenzymes NADH and FAD to be implemented with perfect image coregistration, overcoming the difficulties associated with differences in absorption spectra and disparity in concentration. We demonstrate ratiometric redox imaging free of motion artifacts and simultaneous two-photon fluorescence lifetime imaging (FLIM) of NADH and FAD in living tissues. The lifetime gradients of NADH and FAD associated with different cellular metabolic and differentiation states in reconstructed human skin and in the germline of live C. Elegans are thus simultaneously measured. Finally, we present multicolor imaging of endogenous fluorophores and second harmonic generation (SHG) signals during the early stages of Zebrafish embryo development, evidencing fluorescence spectral changes associated with development.Multiphoton microscopy is a powerful tool for label-free and non-invasive functional imaging in small organisms and tissues 1, 2 . Pulsed near infrared excitation light allows in-depth imaging based on contrasts such as endogenous fluorescence 2 , second harmonic generation (SHG) 3 and third harmonic generation (THG) 4 . Endogenous fluorescence in living tissues arises from several intrinsic biomarkers that play important roles in physiological processes 2 . The primary intracellular sources are NAD(P)H and FAD, the two major cofactors of redox reactions in the cell and central regulators of energy production and metabolism 5,6 . Their fluorescence reports on the metabolic activity of cells allowing tissue physiology and processes such as stem cell differentiation, cancer development and neurodegenerative diseases to being non-invasively monitored [7][8][9][10][11][12] . The fluorescence lifetimes of NADH and FAD are different upon binding to the protein during the electron transport chain. FLIM provides sensitive measurements of the free and protein-bound NAD(P)H ratio and of the redox states (NADH/NAD + ) of cells, and can be used to distinguish glycolytic and oxidative phosphorylation metabolic states [13][14][15][16][17] . Monitoring lifetime of free and protein-bound FAD has also been exploited to quantify redox ratio FAD/FADH 2 , and used as a b...

This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Characterizing melanins in situ and determining their 3D z-epidermal distribution is paramount for understanding physiological/pathological processes of melanin neosynthesis, transfer, degradation or modulation with external UV exposure or cosmetic/pharmaceutical products. Multiphoton fluorescence intensity- and lifetime-based approaches have been shown to afford melanin detection, but how can one quantify melanin in vivo in 3D from multiphoton fluorescence lifetime (FLIM) data, especially since FLIM imaging requires long image acquisition times not compatible with 3D imaging in a clinical setup? We propose an approach combining (i) multiphoton FLIM, (ii) fast image acquisition times, and (iii) a melanin detection method called Pseudo-FLIM, based on slope analysis of autofluorescence intensity decays from temporally binned data. We compare Pseudo-FLIM to FLIM bi-exponential and phasor analyses of synthetic melanin, melanocytes/keratinocytes coculture and in vivo human skin. Using parameters of global 3D epidermal melanin density and z-epidermal distribution profile, we provide first insights into the in vivo knowledge of 3D melanin modulations with constitutive pigmentation versus ethnicity, with seasonality over 1 year and with topical application of retinoic acid or retinol on human skin. Applications of Pseudo-FLIM based melanin detection encompass physiological, pathological, or environmental factors-induced pigmentation modulations up to whitening, anti-photoaging, or photoprotection products evaluation.

Quantifying skin aging changes and characterizing its 3D structure and function in a non-invasive way is still a challenging area of research, constantly evolving with the development of imaging methods and image analysis tools. In vivo multiphoton imaging offers means to assess skin constituents in 3D, however prior skin aging studies mostly focused on 2D analyses of dermal fibers through their signals’ intensities or densities. In this work, we designed and implemented multiphoton multiparametric 3D quantification tools for in vivo human skin pigmentation and aging characterization. We first demonstrated that despite the limited field of view of the technic, investigation of 2 regions of interest (ROIs) per zone per volunteer is a good compromise in assessing 3D skin constituents in both epidermis and superficial dermis. We then characterized skin aging on different UV exposed areas—ventral and dorsal forearms, face. The three major facts of aging that are epidermal atrophy, the dermal–epidermal junction (DEJ) flattening and dermal elastosis can be non-invasively quantified and compared. Epidermal morphological changes occur late and were only objectified between extreme age groups. Melanin accumulation in suprabasal layers with age and chronic exposure on ventral and dorsal forearms is less known and appears earlier. Superficial dermal aging changes are mainly elastin density increase, with no obvious change in collagen density, reflected by SHGto2PEF ratio and SAAID index decrease and ImbrN index increase on all skin areas. Analysis of the z-dermal distribution of these parameters highlighted the 2nd 20 µm thickness normalized dermal sub-layer, that follows the DEJ shape, as exhibiting the highest aging differences. Moreover, the 3D ImbrN index allows refining the share of photoaging in global aging on face and the 3D SAAID index on forearm, which elastin or fibrillar collagens densities alone do not allow. Photoaging of the temple area evolves as a function of chronic exposure with a more pronounced increase in elastin density, also structurally modified from thin and straight elastic fibers in young volunteers to dense and compact pattern in older ones. More generally, multiphoton multiparametric 3D skin quantification offers rich spatial information of interest in assessing normal human skin condition and its pathological, external environment or product induced changes.

This study shows that multiphoton microscopy associated to specific 3D image processing tools allows cutaneous effects induced by topical retinoids in this in vivo model to be non-invasively detected, quantified and followed over time. This innovative approach could be applied to the evaluation of other active compounds.