Translator Disclaimer
Open Access
23 October 2019 Characterization of a multicore fiber image guide for nonlinear endoscopic imaging using two-photon fluorescence and second-harmonic generation
Anikitos Garofalakis, Sergei G. Kruglik, Tigran Mansuryan, André Gillibert, Luc Thiberville, Frédéric Louradour, Christine Vever-Bizet, Geneviève Bourg-Heckly
Author Affiliations +

Multiphoton microscopy (MPM) has the capacity to record second-harmonic generation (SHG) and endogenous two-photon excitation fluorescence (2PEF) signals emitted from biological tissues. The development of fiber-based miniaturized endomicroscopes delivering pulses in the femtosecond range will allow the transfer of MPM to clinical endoscopy. We present real-time SHG and 2PEF ex vivo images using an endomicroscope, which totally complies with clinical endoscopy regulations. This system is based on the proximal scanning of a commercial multicore image guide (IG). For understanding the inhomogeneities of the recorded images, we quantitatively characterize the IG at the single-core level during nonlinear excitation. The obtained results suggest that these inhomogeneities originate from the variable core geometries that, therefore, exhibit variable nonlinear and dispersive properties. Finally, we propose a method based on modulation of dispersion precompensation to address the image inhomogeneity issue and, as a proof of concept, we demonstrate its capability to improve the nonlinear image quality.



Multiphoton microscopy (MPM) is now a well-established method for tissue structural and functional imaging.14 Studies have been reported in a variety of tissues, from animal models to human samples, and the technique is also used in vivo for dermatological applications.57

MPM is of particular value for imaging extracellular matrix (ECM) architecture in label-free tissues since the combination of second-harmonic generation (SHG) and endogenous two-photon excitation fluorescence (2PEF) allows the simultaneous visualization of collagen and elastin, respectively, the main ECM proteins. Several major lung pathologies such as chronic bronchial inflammations, chronic obstructive pulmonary disease, fibrosis, asthma, emphysema, and cancer are characterized by early modifications of the ECM fibrillary collagen and elastin network. A number of investigations on ex vivo animal811 and human lung samples1218 already demonstrated the potential of the method to assess the pathological changes occurring in ECM organization.

To date, fiber-optic, single-photon confocal autofluorescence endomicroscopy allows in vivo microscopic real-time imaging of the human respiratory tract, during bronchoscopy, from proximal bronchi down to the alveoli. In previous studies,19,20 the authors showed that, under the 488-nm single-photon excitation mode, the technique images the elastin network of the proximal and distal bronchial tree but does not give access to the collagen component of the ECM. Translation of MPM to endoscopy potentially offers one way to overcome this limitation: imaging both collagen (through SHG) and elastin (through 2PEF) networks would provide the capability to follow structural remodeling of the ECM in vivo during an endoscopic procedure, which would be of great clinical value.

There have been significant research efforts21,22 to transfer table-top MPM to clinical endoscopy by developing miniaturized endomicroscopes using optical fibers both to deliver light to the tissue and to collect the backward-emitted multiphoton tissular signals. Two endomicroscopy techniques have been developed, based either on the distal scanning of a single fiber or on the proximal scanning of a fiber bundle. In the first case, the scanning mechanism is placed at the tip of a single fiber, whereas, in the second case, the scanning mechanism, placed at the fiber entrance, sweeps the laser beam across a multicore image guide (IG), ensuring a sequential illumination of the fiber cores. In recent years, research has been focused mainly on the development of endomicroscopes based on the distal scanning of different types of single fibers.2330 These studies have proven the potential of providing uniform images, but they do not fully comply with the U.S. Food and Drug Administration (FDA) and European Commission (CE) conformity regulations required for a clinical use, namely because of the ability to accurately sterilize this kind of miniaturized device has not been demonstrated yet.

On the other hand, fiber-bundle-based endomicroscopes present a number of advantages: as the scanning mechanism and its electrical control are not included in the probe to be incorporated in the operating channel, there is no restriction on scanner size; the probe is compact, mechanically flexible, and robust, and a real-time acquisition rate (10  Hz) can be achieved. Moreover, it can be chemically sterilized against bacteria and viruses, in contrast to distal scanning endomicroscopes, and it allows efficient signal collection using a large surface, i.e., the total bundle area. All these advantageous features increase the clinical applicability of the technique. On the other hand, a limitation of this approach is the limited lateral resolution of multifiber IG due to several micrometers core-to-core distance and optical cross talk between the adjacent fiber cores.

Several studies have assessed the potential of multicore fiber-bundle-based systems for 2PEF and/or SHG imaging. Using a 140-cm-long fiber bundle connected distally to a gradient-index lens and proximally to a pair of diffraction gratings to precompensate the group velocity dispersion occurring in the fiber-bundle cores, images from fluorescein-labeled blood vessels in anesthetized rats have been obtained.31 The 2PEF images from fresh cresyl-violet-stained tissue canine bladder were also reported32 using a 50-cm-long fiber bundle, without compensating for dispersion and nonlinear effects in the fibers. Some of the present authors developed several systems allowing the precompensation of both linear and nonlinear distortions arising in long single-mode fibers (SMF),33,34 and using a specially designed multicore fiber bundle from Fujikura, 2PEF images of human colon cells stained with Rhodamine B were obtained.33 Recently, a system employing a fiber bundle has been developed enabling the simultaneous imaging of coherent anti-Stokes Raman scattering, SHG, and 2PEF in a human dermis sample.35,36 The bundle, used only to conduct the excitation wavelengths, is an IG consisting of 10,000 light-guiding cores with a relatively small length of 1 m, and thus, no compensation for dispersion in the fibers was implemented.

However, to the best of our knowledge, no attempt has been reported to perform joint imaging of a tissue by SHG and endogenous (label-free) 2PEF using a fiber-bundle-based endomicroscope directly operable in a routine endoscopic procedure. This means that the instrument should provide real-time imaging to the practitioner and avoid motion artifacts induced by heartbeat and breathing; the fibered probe should be flexible, small-sized, and long enough to be easily inserted in the operating channel of an endoscope, guided along curved lumen of hollow organs and put into contact with the tissue. Moreover, the probe should be fully biocompatible, resistant to decontamination processes and eligible for both 510(k) clearance from the U.S. FDA and the European CE conformity marking. All these features can be found today in single-photon endomicroscopes, but, to our knowledge, these requirements are not met by any multiphoton endomicroscope in development.

The aim of this study is to evaluate the possibility of developing a multiphoton fiber-bundle-based endomicroscope, meeting all the requirements for clinical use, and thus immediately transferable into endoscopy at the patient bedside. This endomicroscope would be intended primarily for imaging ECM of elastin and collagen networks. For this purpose, we built an endomicroscope based on a commercial fiber bundle and present the acquired endogenous 2PEF and SHG images from ex vivo tissue samples. The main focus has been put into quantitative characterization of the properties of this fiber bundle as a nonlinear IG.


Materials and Methods


Multicore Fiber Image Guide

The IG used in our study is the Gastroflex™ UHD confocal miniprobe™ commercialized by Mauna Kea Technologies for use with the single-photon fluorescence endomicroscope Cellvizio®. This miniprobe is based on a Fujikura multicore fiber-optic bundle of the FIGH series. The bundle is composed of 30,000 separate optical cores, each one having a diameter no greater than 2  μm, with center-to-center intercore distance of about 3.5  μm. The effective working field of view of this miniprobe has a diameter of 600  μm and its total length is 3 m. Therefore, it can be used in an operating room, regardless the tissular target to be reached. The outer diameter of the IG is 2 mm, which makes it compatible with any endoscope-operating channel having a minimal diameter of 2.8 mm. The distal end of the bundle is connected to a proprietary micro-objective (MO) with working distance of 60  μm and magnification of 2.48, resulting in a lateral resolution of 1.4  μm and an axial resolution (optical section thickness) of 9  μm (specifications at 488 nm).


Tissue Sample Preparation

The human lung samples used in this investigation were prepared by Rouen University Hospital and were obtained from one patient undergoing surgery for lung cancer. The specimens were removed from a lobectomy sample in a healthy zone, remote from the cancer, and fixed in 4% paraformaldehyde. After fixation, samples were cut into 1-mm-thick tissue slices. Image acquisitions were performed on a healthy peripheral alveolar territory sample. Each tissue block was placed on a microscope slide.

Rat tail tendons were dissected from thawed rat tails and stored in phosphate-buffered saline solution (to remove hemoglobin traces) for several minutes before image acquisitions. A part of the tendon was placed on a microscope slide.


Experimental Setups

Two variations of the experimental setup were built: one for 2PEF and SHG imaging of tissue using fiber-bundle IG, and another one for quantitative characterization of the IG at the single-core level.

The setup is shown in Fig. 1. A tunable (690 to 1040 nm) Ti:sapphire laser (L) having 80-MHz repetition rate and 70-fs pulse duration (Spectra-Physics Mai Tai/DeepSee) is used as the excitation light source. A half-wavelength plate (WP1, Thorlabs AHWP05M-980) together with the Faraday isolator (FI, Newport ISO-05-800-BB-G) are placed at the laser output to vary the pulse energy and also to prevent the parasitic feedback by reflected light. Then, via a pair of reflecting mirrors (M1 and M2, Thorlabs PF10-03-P01), the light enters a home-built pulse-shaping device,34,37 which permits simultaneous second-order dispersion (SOD), third-order dispersion (TOD) dispersion, and nonlinear effect precompensation. This device was initially developed for delivering ultrashort pulses through a long SMF for various endoscopic applications.34,3740 It consists of (i) a short SMF (Thorlabs PM780-HP) with its focusing optics (L1 and L2, Thorlabs A260TM-B) and polarization rotation management (WP2, Thorlabs AHWP05M-980), (ii) mirrors for light delivery (M3 to M6, Thorlabs PF10-03-P01), and (iii) a grism-based pulse stretcher (G1 and G2) consisting of a pair of reflection diffraction gratings in close contact with two prisms.

Fig. 1

Experimental setup employing an endoscope-compatible multicore IG for ex vivo 2PEF and SHG imaging of tissular samples and characterization of the same IG at the single-core level. (a) The part of the setup that is common for the two type of measurements; L, a femtosecond laser; WP1 and WP2, half-waveplates; FI, Faraday isolator; M1 to M7, mirrors; L1 to L4, lenses; SMF, single-mode fiber; and G1 and G2, pair of grisms. (b) Module dedicated to 2PEF and SHG imaging; DM, dichroic mirror; F1, dichroic short-pass or bandpass filter; SM, set of two scanning mirrors; OBJ, objective; IG, multicore image guide; and MO, micro-objective lens. The tissue sample, S, is placed on a 3-D stage. The dashed line frame depicts the components composing the pulse-shaping device. (c) Module dedicated to the IG characterization at the single-core level; WP3, quarter-waveplate; IG, multicore image guide; OBJ1 and OBJ2, objectives on 3-D stages; P, polarizer; PM, power meter; BS, beam splitter; L8, lens; F2, neutral density filter; CMOS, video camera; WP4, half-waveplate; AC, autocorrelator; and M10 to M12, mirrors.


The principle of operation of this precompensation system has been previously described in detail.33,34 Briefly, a short piece of SMF introduces spectral and temporal pulse broadening under the combined effects of self-phase modulation and positive (i.e., normal) group-velocity dispersion, effectively decreasing the light intensity entering the first grism (G1). Then, an anomalous dispersive stretcher consisting of two grisms (G1 and G2) produces, in two passes, largely negatively chirped pulses. As a result, after propagation inside the long delivery fiber with positive dispersion and nonlinearity, the pulse impinging on the tissue can be optimally compressed to a value close to or even shorter than the initial laser pulse duration. This is achieved by proper adjustment of the anomalous stretcher, i.e., by changing the distance d between the two grisms and the angle of incidence θ, one can control the desired values of SOD and TOD, respectively. Although a pulse duration of about 10 picosecond (ps) is measured at the stretcher’s exit (between mirrors M5 and M7), a value in the range 50 to 100 fs can be routinely obtained after passing through a single-mode 3-m-long fiber. Note that our grisms-based precompensation device is tunable from 790 to 900 nm.39 An excitation wavelength of 850 nm, optimal for suppressing the 2PEF signal from collagen,39 was used in this study.

After the pulse-shaping device, the beam passes via an optical delivery system consisting of a mirror (M7, Thorlabs PF10-03-P01) and two lenses (L3, Thorlabs AC254-050-B and L4, Thorlabs AC254-100-B) that act as a telescope for expanding the laser beam and adjusting its divergence. Then, depending on the type of experiment, the beam follows different trajectories.


Module for ex vivo two-photon excitation fluorescence and second-harmonic generation imaging of tissular samples

For the 2PEF and SHG imaging of tissue, the M8 and M9 mirrors guide the beam to the module presented in Fig. 1(b). The collimated light passes through a dichroic mirror (DM, Semrock FF665-Di02-25×36) toward a proprietary scanning module (Mauna Kea CellVizio®) composed of two scanning mirrors (SMs), a telescope (L5 and L6), and a coupling objective (OBJ). The laser beam is swept across the proximal end of the multicore IG (Gastroflex™ UHD miniprobe™) ensuring sequential illumination of all fiber cores. Horizontal line scanning is performed using a 4-kHz resonance mirror while a galvanometric mirror performs frame scanning. The dwell time of the laser light per core is <1  μsec and varies depending on the area scanned. The minimal frame rate of 9 Hz is set for compatibility with in vivo in-clinic endoscopic application, which requires real-time imaging during intervention on patients. After propagating through the IG, the beam is focused onto the sample tissue (S) by a proprietary MO, which is an integral part of the employed IG. The sample is placed on a three-dimensional (3-D) differential micrometer translation stage; the horizontal in-plane movement is used to manually select the sample’s region of interest, whereas the vertical movement enables probing of different tissue depths.

The total average power on the sample is deliberately limited to 40 mW, corresponding to a pulse energy of 0.5 nJ. The 2PEF and SHG signals are collected by the whole IG in epi-detection configuration, delivered back to the proximal end of the IG, descanned by SMs (as for usual linear fluorescence imaging but without pinhole filtering), reflected by the dichroic beamsplitter (DM) and focused on a fast photomultiplier (PMT, Hamamatsu H10723-110MOD) by a lens (L7) with f=50  mm. Dichroic filter is placed before the PMT for spectral selection of the informative signal: a low-pass filter (Semrock FF01-720/SP-25) allows simultaneous recording of 2PEF and SHG signals, whereas a bandpass filter centered at 427 nm with 10-nm FWHM (Semrock FF01-427/10-25) allows acquisition of the SHG signal exclusively. The PMT signal is then synchronously processed using a 14-bit digitizer, and two-dimensional images are reconstructed using a proprietary software provided by Mauna Kea Technologies.

The recorded videos are then subjected to postprocessing treatment frame by frame. First, a 2×2 binning is applied for reducing the read noise. Then a circular averaging filter with a radius of 7.5  μm is applied for smoothing the image.


Module for image guide characterization at single-core level

For the quantitative characterization of the IG at the single-core level, the M8 mirror is removed and the beam exiting the telescope (L3 and L4) is focused onto the proximal end of the fiber-bundle IG by a 20× objective (OBJ1, Nikon Plan, NA 0.4). The proximal end of the fiber bundle is positioned on a 3-D stepper-motor translation stage with a minimal incremental motion of 50 nm and unidirectional repeatability of 100 nm (three linear stages M-111.12S and stepper motor controller C-600, Physik Instrumente). The choice of an arbitrary single core is performed by moving the stage along the XY plane. Then, translation along the optical (Z) axis is used to optimize the coupling into this core. Various cores can be selected by moving the 3-D stage in the XY plane without changing the Z axis position. Light exiting the distal end of the IG is collimated by a 50× objective (OBJ2, Olympus MS Plan, NA 0.80) and then split by a beam sampler (BS, Thorlabs BSF10-B) into two measurement channels containing relevant optical elements: a half-waveplate (WP4, Thorlabs AHWP05M-980) before an autocorrelator (AC), a lens, L8, and a neutral density filter (F2 Thorlabs NE10A-B) before a CMOS camera. Both XY and Z positions of the stepper motor are fine-tuned to optimize the signal of the cores on the CMOS camera. The process of this optimization consists in looking at the spatial distribution of one core on the CMOS camera. Then the stepper motor moves to maximize the core intensity and minimize the cross talk of the targeted core with respect to the adjacent ones. Optionally, a polarizer (P, Thorlabs LPNIR050) can be placed in front of the mirror, M11, and a quarter-waveplate (WP3, Thorlabs AQWP05M-980) before the OBJ1.

The above setup was used for measuring the following parameters.

  • a. Average power of transmitted light with a power meter (PM, Newport 1918C with the 918D-SL-OD3R measuring head) from which the pulse energy was calculated. The power was usually limited to 10 mW, corresponding to pulse energy of 125 pJ.

  • b. Pulse duration, by AC (APE PulseCheck, long-pulse version: 50 fs to 50 ps). Signal is maximized by rotating the half-waveplate WP4.

  • c. The degree of polarization (defined as Pol=[ImaxImin]/[Imax+Imin]), by inserting and rotating the polarizer P in front of the PM. For measurements with circularly polarized input light, the quarter-waveplate (WP3) was placed before the first objective, OBJ1.

  • d. Beam spatial profile, by CMOS camera (Thorlabs DCC1545M). The lens L8 with f=50  cm is used for focusing the parallel beam into the camera, providing real-time images of the cores with a field of view of 25×25  μm2. The neutral density filter F2 prevents signal saturation on the camera.



Spearman’s correlation coefficients were computed to assess the correlations of pulse energy and degree of polarization with pulse duration. The p-value and 95% confidence intervals (CIs) were estimated by bias-corrected accelerated (BCa) bootstrap according to Efron’s method.41 A correlation coefficient was considered negligible if it was below 0.30, weak if it was between 0.31 and 0.70, strong if it was between 0.71 and 0.90, and very strong if it was above 0.91. Correlation coefficients were compared to the value 0 (null hypothesis). A p-value below 0.05 was considered statistically significant. A semilogarithmic scale was used for pulse duration in scatter plots since the pulse duration distribution was highly skewed. All statistics were computed in R statistical software (version 3.5.0, The R Foundation for Statistical Computing, Vienna, Austria).


Results and Discussion


Unlabeled Tissues Imaging

Two types of dissected samples were imaged with the first configuration of the setup; (a) 2PEF of human alveolar territory sample without any staining, and (b) SHG from a rat tail tendon (Fig. 2).

Fig. 2

The 2PEF and SHG ex vivo real-time imaging with the use of endoscopic multicore IG. (a) Still image of the 2PEF signal taken from a cine loop of 9 Hz of a human alveolus territory from a biopsy sample (see Video 1, MP4, 8.88 MB [URL:]). (b) Still image taken from a cine loop of 9 Hz showing the SHG signal of the tendons of rat tail excised sample (see Video 2, MP4, 4.58 MB [URL:]). Additional 2PEF experiments using fluorescein-labeled tissue were performed and the signals obtained were stronger (see Video 3, MP4, 1.00 MB [URL:]).


These images were recorded in selected regions of interest with the size 190×140  μm (see also Videos 1 and 2). Elastin, which is found in abundance in human alveoli, is considered as the main source of 2PEF signal; the image in Fig. 2(a) provides information on the structure of elastin network. In Fig. 2(b), SHG signal originates from tendons rich in collagen.

Then, we sought to explore the possibility of parallel assessment of 2PEF and SHG signals within the same region of a sample (Fig. 3). A dissected portion of a rat tail was continuously imaged and manually scanned until a tendon structure appeared in the SHG spectral window [see Fig. 3(a)]. Then without moving the sample, a low-pass filter was positioned (see Sec. 2.3.1), and the SHG filter was removed in front of the PMT. From the recorded image containing both 2PEF and SHG signals (data not shown), we performed a pixel-to-pixel subtraction of the image shown in Fig. 3(a) to obtain the signal of 2PEF [Fig. 3(b)].

Fig. 3

(a) SHG and (b) 2PEF in-vitro imaging of a tendon from a dissected rat tail sample. The module of Fig. 1(b) has been used for this acquisition. Both images correspond to the same field of view of the sample. (c) Fusion image of (a) and (b).


Our approach aimed to be conformed with the medical regulations related to standard endoscopic procedures, meaning that the fiber guide had to be approved for clinical endoscopic imaging. Moreover, real-time imaging is indispensable to endoscopic practices. The above results presented in Figs. 2 and 3 are, to the best of our knowledge, the first images of 2PEF of elastin and SHG of collagen from ex vivo tissue without any external labeling, using a multicore IG and at endoscopy-relevant frame rate. The images shown here were subjected to postprocessing by applying a combination of numerical filters, as described above. However, the time of image treatment was significantly shorter than the time interval between two successive frame acquisitions at 9 Hz. Therefore, such a treatment can be directly integrated into real-time monitoring of 2PEF and SHG signals during in-clinic endoscopy intervention.

Although these first results are promising, the recorded images are far from being optimal. First, the signals, especially in the case of SHG, were relatively weak and, second, the images exhibited high spatial inhomogeneity contrary to images provided by distal imaging schemes.

Regarding signal intensity, a number of factors contribute to the observed problem. First, we chose an endoscopically relevant fast frame rate of 9 Hz which translates in a very short exposure time of laser excitation per fiber in one frame. As a result the non-linear signals from nonlabeled tissue are weak. We performed additional 2PEF experiments using fluorescein-labeled tissue and we obtained much stronger signals with better quality of images (Video 3). On the other side, the total average power of femtosecond radiation at the distal end of an IG was limited at a certain level (40 mW) and cannot be substantially increased to avoid tissue damage. Importantly, in this study, we used the available scanning module from Mauna Kea Technologies, the optics of which were optimized for linear excitation in the visible range and not for nonlinear regime using near-infrared (NIR) light at 850 nm. Moreover, we would like to stress that the signal intensity is not the most important parameter per se, as useful information is derived from the morphology of elastin and collagen networks. To accurately resolve the network, a few dozen photons per pixel above the dark current threshold is sufficient. Therefore, we believe that gradual improvement of both scanning and registration modules will allow one to obtain signals sufficiently strong for real-time in-clinic endoscopic imaging of unlabeled tissue.

The second major issue, which is image spatial inhomogeneity, seems to be much more important, and, as described in the following sections, further experiments were performed for the better understanding of the underlying problems.


Uniform Sample Imaging

The postprocessing treatment improves image quality in Figs. 2 and 3; however, the image inhomogeneity is still present. To investigate this problem in detail, experiments were performed on homogeneous samples to examine the signal intensity variation between the cores. The results are shown in Fig. 4.

Fig. 4

The distribution of the intensities of the cores of the IG. (a) An image of 2PEF measurement of a fluorescein solution. (b) The corresponding image for the case of a SHG measurement of a BBO crystal. (c), (d) The distributions of the cores based on their intensities following their segmentation as shown in (a) and (b). The range of values for the (a) image was from 811 to 1760 a.u. Similarly, the range of values for the (b) image was from 1030 to 5180 a.u.


Fluorescein water solution (2.5  μM) and barium borate (α-BBO) crystal were used for obtaining reference images of 2PEF [Fig. 4(a)] and SHG [Fig. 4(b)], respectively. In both cases, the recorded images are surprisingly highly heterogeneous; the corresponding intensity histograms are shown in panels of Figs. 4(c) and 4(d). The core intensities exhibit big variability. In the linear regime with continuous-wave (cw) excitation, the variability in transmission of individual cores also exists. However, it is much lower and, therefore, can be compensated by proper image scaling procedure based on calibration using a fluorophore solution as a standard. However, in the case of nonlinear imaging, an analogous scaling procedure is not immediately possible as it leads to very large correction errors due to very large spread of intensities between the cores.


Multicore Fiber Image Guide Characterization at Single-Core Level

We examined the origin of nonlinear image heterogeneity by studying the IG at the single-core level using the module described in Sec. 2.3.2. In total, 107 different individual cores randomly located within the IG were examined. The distributions of pulse duration, pulse energy, and degree of polarization at the distal exit are presented in Fig. 5. Pulse duration of a substantial part of examined cores lies between 70 and 500 fs (n=30, 28%), and, within this group, 53% of cores exhibit pulse durations between 100 and 200 fs [insert in Fig. 5(a)]. However, for the majority of cores (n=64, 60%), pulse duration is measured to be longer than 1 ps and may even reach 7.4 ps. This means that the precompensation stretcher device (G1 and G2) does not properly compensate the SOD and TOD in the multicore IG for the ensemble of cores but only for some of them. These longer pulse durations translate into almost zero-intensity image points as both 2PEF and SHG signals depend inversely on the pulse duration.

Fig. 5

Histogram plots from parameters measured by the module of the setup of Fig. 1(c) designed for the evaluation of the IG at the single-core level. (a) Distribution of pulse duration, (b) transmitted pulse energy, and (c) degree of polarization over the single cores of multicore IG. The insert in panel (a) shows the extended distribution of pulse duration between 70 and 500 fs.


The pulse energy distribution is shown in Fig. 5(b). The energy of incident pulses was set to 130±5  pJ at the proximal entrance of the IG. At the distal exit, the pulse energy is measured to be within the range from 10 to 35 pJ. About 35% of the cores (n=36) exhibited output energies between 25 and 30 pJ, and the maximum value is measured at 35 pJ, giving a coefficient of transmission of about 27%. Note that in this study we do not focus on energy transmission optimization but rather on a robust comparative characterization of different individual cores. Clearly, the pulse energy distribution is rather compact and does not exhibit big variations as compared to those of pulse duration. Finally, when pulse energies are plotted against pulse durations [Fig. 6(a)] of the corresponding cores, the correlation is zero or negligible (rho=0.01, 95% CI=0.19 to +0.16).

Fig. 6

Scatter plots of (a) pulse energy and (b) degree of polarization against pulse duration. In each case, the Spearman’s correlation coefficient (rho) and the 95% CI are given.


To study the birefringence due to light propagation in a single core of the IG, we measured the degree of polarization of transmitted pulses and its distribution, as shown in Fig. 5(c). It varies between 3% and 55%, with 59% of the cores (n=63) exhibiting a value below 25%. This result confirms that cores are optically anisotropic and alter the linear polarization of an incident light. At the same time, no correlation or a negligible correlation (rho=+0.07, 95% CI=0.14 to +0.28) between the degree of polarization and pulse duration was found [Fig. 6(b)], suggesting that the induced birefringence is not the main cause of pulse duration variations.

To verify the above conclusion, a modified version of the precedent experiment was performed. The polarization of an incident beam was changed from linear to circular using the quarter-waveplate [WP3, Fig. 1(c)] inserted just before the proximal side of the IG. Durations of the output pulses with and without the quarter-waveplate were measured for 54 individual cores, and the result is shown in Fig. 7. The change from linear to circular polarization of the incident beam does not substantially alter the output pulse durations’ distribution, which remains very broad, i.e., from 70 fs to 7.4 ps. This result confirms that the state of polarization during propagation along the fiber core does not influence in a major way the transmitted pulse duration.

Fig. 7

The effect of circular versus linear input polarization on pulse duration of transmitted light. Gray bars represent distribution of pulse duration in the case of linearly polarized incident light, and black bars represent the case of circular polarized incident light.


Note that the difference between the two distributions of pulse duration for linearly polarized incident light [Fig. 5(a) versus the gray bars of Fig. 7] is due to different individual cores studied in each experiment. One can clearly see that, in spite of substantial local differences, the general tendency remains the same: both distributions are very broad, extending 2 orders of magnitude, from 70 fs to about 7 ps.

Finally, we examined the cross talk between adjacent cores. For this study, CMOS camera images of spatial distribution of transmitted core intensities were used; two representative examples are shown in Fig. 8. In both images, the strong peak from the studied core is located in the middle of the camera, whereas adjacent cores with a much lower intensity are also visible around the targeted core. Figure 8(a) shows the particular case of strong cross talk where the intensities from adjacent cores amount up to 44% of the targeted core peak. However, the latter case is rare and for the large majority of studied cores, the cross talk effect is minimal, with adjacent peak intensities below 10% of the main peak; a typical example is shown in Fig. 8(b). Intensity profiles along the white lines in Figs. 8(a) and 8(b) are shown in Figs. 8(c) and 8(d), respectively. These data reveal insignificant cross talk between the adjacent cores for Gastroflex™ UHD confocal miniprobe™ IG for NIR light at 850 nm.

Fig. 8

The cross talk effect of the IG. (a) Image showing a single core of the IG as measured by the CMOS camera of the setup of the Fig. 1(c) as laser light passes through the cores. In this image, a cross talk signal from the adjacent cores is shown. (b) Image showing a single core with typical cross talk signal from the adjacent cores. (c), (d) The distributions of the signal of the cores along the white lines of (a) and (b), respectively.


There is in fact a trade-off between core uniformity and cross talk; if the transversal core geometries were exactly the same (i.e., cores with identical phase velocities), the cross talk between the adjacent cores would be very strong because of intercore coupling. Fujikura fiber-optic bundle of the FIGH series was specifically designed to reduce the cross talk, thus increasing the resolution of the IG. This is achieved by manufacturing every adjacent core with statistically different transversal profile in size and shape, as depicted in the onset of Fig. 1(b). This intended and well-controlled inhomogeneity is not detrimental for intensity transmission in the linear imaging mode because it can be compensated through image processing. In the case of nonlinear excitation using 2PEF and SHG with femtosecond pulses, the situation is different. The interaction between the propagating ultrafast light pulse and the core medium depends on the square of intensity and is detrimentally affected by core geometry inhomogeneity.

For a given single core, by fine-tuning the precompensation stretcher, pulse durations in the order of 100 fs can always be achieved at the exit of the core. Indeed, as mentioned in the description of the grism-based stretcher, the value of the introduced anomalous dispersion is adjusted in order to correctly compensate the positive dispersion in the cores of an IG. Modifying the distance between the two grisms (d) and/or the incidence angle (θ) onto the stretcher [see Fig. 1(a)] permits independent tuning of SOD and TOD to the desired values. The pulses that are anomalously dispersed at the exit of the stretcher can be temporally compressed during their propagation through any chosen core of the bundle.

In the case of our multicore IG, the stretcher should be tuned in a way for precompensating the dispersion of all cores. In our approach, we chose to keep the incidence angle θ constant and tune the distance d between the grisms [Fig. 1(a)]. The criterion for the optimal distance d was chosen based on mean intensity of the whole IG. However, each core responds individually to the distance d, i.e., values of d that maximize the intensity for some cores are not optimal for the others. This observation shows that each core has its own dispersion and nonlinearity. We conclude that optimal precompensation cannot be achieved simultaneously for the ensemble of cores in the bundle.


Effect of Modulating the Intergrism Distance on the Nonlinear Image Homogeneity

Following the above findings, we propose a method for improving the homogeneity of the acquired nonlinear images using a modulation of the dispersion precompensation. This approach is based on the observation that different d values (hence precompensations of SOD) are optimal for different groups of cores. As a proof of concept, an experiment was performed where the distance d was changed in five steps (Δd=0.5  mm) from its initial value. Images of 2PEF from a fluorescein aqueous solution were recorded at all six precompensation settings. Figure 9(a) shows the whole recorded image at initial d position. A region of interest of 30×20  μm was chosen [frame in Fig. 9(a)] and shown in Fig. 9(b) for the six different distances d between the grisms. At each d value, the signal of core intensities varies, as can be seen in Fig. 9(d), for the three chosen cores marked by orange circles in the upper panel of Fig. 9(b). By superposing and averaging the six images of Fig. 9(b) pixel by pixel, a synthetic mean image can be created [Fig. 9(c)], which is more homogeneous than any of the constituent images. Indeed the aggregated image in Fig. 9(c) exhibited less variance in the core intensities (σ2=1.3×103) than any of the six images in Fig. 9(b) (σ2=2.7×103 to 5.3×103).

Fig. 9

A method for improving the homogeneity of nonlinear images based on modulation of the distance d between the grisms, hence on modulation of precompensation of SOD. (a) A 2PEF image of a fluorescein solution (2.5  μM). (b) Six zoomed images corresponding to six different distances d between the grisms. The parameter d changes in steps of 0.5 mm. (c) Composed image by the mean intensity projection of the six images of (b). (d) The variations of the signals for the three encircled cores of (b). The different points correspond to the different distances between the grisms.



Future Perspectives for Improving the Image Homogeneity

We demonstrated that modulation of the dispersion precompensation is a valuable method for improving the image quality. For implementing this method in our setup, an accessible technical solution would consist in modulating the dispersion at high frame rate with a variability amount defined by the results in Fig. 9. In this case, several successive images, having complementary inhomogeneities, would have to be properly associated in order to form a combined final image with optimal homogeneity. One possible discrete dispersion modulation arrangement could be achieved using a relatively simple optomechanical device in combination with the already-existing pulse stretcher. This additional device could be composed of a set of parallel glass plates of different thicknesses positioned in front of a mirror in the image space of a converging lens. A one-dimensional galvanometric mirror, synchronized with the imaging system frame rate and positioned in the focal plane of the lens, could allow quick sequential illumination of various plates, one image frame corresponding to one glass thickness. Due to the above mirror, the double-pass configuration would allow avoiding the transverse displacement of the beam during dispersion modulation.

Other realizations of precompensation modulation, for example, using acousto- or electro-optical beam deflection, are also possible and will be the subject of further studies.

In the experiment depicted in Fig. 9, we note also the existence of cores with very low intensity, regardless of the value d. This phenomenon is attributed to nonoptimal coupling of NIR light into the IG since both the scanning module and the IG were initially designed for cw excitation in the visible range. Indeed, from Fig. 7(b) one can see that, in the conditions of optimal static single-core coupling, the variance of the transmitted energy is rather low; thus, the NIR pulsed laser radiation can in principle be efficiently coupled into any core of the entire IG, provided the optimal matching conditions are satisfied during the proximal beam scanning. Therefore, a new scanning module specifically optimized for NIR laser light and synchronized with an ad hoc precompensation modulation system is to be designed and built to achieve the optimal image quality during real-time in vivo endoscopy.

It is noteworthy that, for an IG with high intrinsic inhomogeneity, the proposed precompensation modulation method leads either to substantial lowering of the useful frame rate or to the reduction of the exposure time in one core. Therefore, in the future, instead of working with a commercially available imaging guide being fabricated and distributed for decades by Fujikura Ltd, one may consider developing a new bundle possessing an optimal equilibrium between cores homogeneity, cores distance, and cross talk between cores that effectively reduces the spatial resolution of the resulting nonlinear image. For such an “optimal” IG, the proposed precompensation modulation procedure is expected to be much faster or eventually may even not be required.

Note that several research groups across the world42,43 are now actively working on the development of new imaging multicore bundles with increased resolution and homogeneity. Such IGs should exhibit diameter sizes and mechanical properties that would make them good candidates to be inserted in commercial endomicroscopic probes currently manufactured. Future development of new optical fiber technologies, in particular one based on air–silica microstructures,42 represents, in our opinion, one promising route toward nonlinear endomicroscopy and its effective application in the clinic.



We report a nonlinear endomicroscopy apparatus employing multicore fiber IG compatible with real-time in-clinic endoscopy procedure. The system allowed us to obtain the first real-time cine loops of 2PEF of elastin and SHG images of collagen from dissected tissues of human alveoli and rat tail tendon, respectively, without any labeling. A frame rate of 9 Hz was achieved with a total excitation power of 40 mW at 850 nm. These results suggest good potential for future in-clinic application of nonlinear fiber-bundle-based endomicroscopy, although the obtained image quality is far from optimal, being highly inhomogeneous. To examine this, we performed a rigorous study of the multicore IG employed (Fujikura FIGH series) at a single-core level. We concluded that the nonlinear image inhomogeneity is primarily due to strong variation across the bundle of core transversal geometries, the various cores having variable dispersions. Finally, we proposed a correction method based on modulation of dispersion precompensation to address the image inhomogeneity issue and demonstrated its capability to improve the nonlinear image quality.


The authors declare that there are no conflicts of interest related to this article.


We are grateful to Mauna Kea Technologies for providing the scanning module and the image acquisition software in the frame of a former project funded by the Agence Nationale de la Recherche (ANR).The major part of the experimental equipment used in this study was funded by the ANR (project ONL-in vivo ANR-08-TECS-0006-01).



M. Strupler et al., “Second harmonic imaging and scoring of collagen in fibrotic tissues,” Opt. Express, 15 (7), 4054 –4065 (2007). OPEXFF 1094-4087 Google Scholar


J. Chen et al., “Multiphoton microscopy study of the morphological and quantity changes of collagen and elastic fiber components in keloid disease,” J. Biomed. Opt., 16 (5), 051305 (2011). JBOPFO 1083-3668 Google Scholar


F. Aptel et al., “Multimodal nonlinear imaging of the human cornea,” Invest. Ophthalmol. Visual Sci., 51 (5), 2459 –2465 (2010). IOVSDA 0146-0404 Google Scholar


J. Lin et al., “Assessment of liver steatosis and fibrosis in rats using integrated coherent anti-Stokes Raman scattering and multiphoton imaging technique,” J. Biomed. Opt., 16 (11), 116024 (2011). JBOPFO 1083-3668 Google Scholar


K. König et al., “Clinical two-photon microendoscopy,” Microsc. Res. Tech., 70 (5), 398 –402 (2007). MRTEEO 1059-910X Google Scholar


H. G. Breunig et al., “Multiphoton excitation characteristics of cellular fluorophores of human skin in vivo,” Opt. Express, 18 (8), 7857 –7871 (2010). OPEXFF 1094-4087 Google Scholar


K. Koenig, “Hybrid multiphoton multimodal tomography of in vivo human skin,” IntraVital, 1 (1), 11 –26 (2012). Google Scholar


A.-M. Pena et al., “Three-dimensional investigation and scoring of extracellular matrix remodeling during lung fibrosis using multiphoton microscopy,” Microsc. Res. Tech., 70 (2), 162 –170 (2007). MRTEEO 1059-910X Google Scholar


S. Zhuo et al., “Nonlinear optical microscopy of the bronchus,” J. Biomed. Opt., 13 (5), 054024 (2008). JBOPFO 1083-3668 Google Scholar


T. Abraham et al., “Minimally invasive multiphoton and harmonic generation imaging of extracellular matrix structures in lung airway and related diseases,” Pulm. Pharmacol. Ther., 24 (5), 487 –496 (2011). Google Scholar


I. Pavlova et al., “Multiphoton microscopy and microspectroscopy for diagnostics of inflammatory and neoplastic lung,” J. Biomed. Opt., 17 (3), 036014 (2012). JBOPFO 1083-3668 Google Scholar


C.-C. Wang et al., “Differentiation of normal and cancerous lung tissues by multiphoton imaging,” J. Biomed. Opt., 14 (4), 044034 (2009). JBOPFO 1083-3668 Google Scholar


T. Abraham and J. Hogg, “Extracellular matrix remodeling of lung alveolar walls in three dimensional space identified using second harmonic generation and multiphoton excitation fluorescence,” J. Struct. Biol., 171 (2), 189 –196 (2010). JSBIEM 1047-8477 Google Scholar


X. Xu et al., “Multimodal non-linear optical imaging for label-free differentiation of lung cancerous lesions from normal and desmoplastic tissues,” Biomed. Opt. Express, 4 (12), 2855 –2868 (2013). BOEICL 2156-7085 Google Scholar


K. Tilbury et al., “Second harmonic generation microscopy analysis of extracellular matrix changes in human idiopathic pulmonary fibrosis,” J. Biomed. Opt., 19 (8), 086014 (2014). JBOPFO 1083-3668 Google Scholar


G. Tjin et al., “Quantification of collagen I in airway tissues using second harmonic generation,” J. Biomed. Opt., 19 (3), 036005 (2014). JBOPFO 1083-3668 Google Scholar


M. Jain et al., “Multiphoton microscopy: a potential ‘optical biopsy’ tool for real-time evaluation of lung tumors without the need for exogenous contrast agents,” Arch. Pathol. Lab. Med., 138 (8), 1037 –1047 (2014). APLMAS 0003-9985 Google Scholar


R. M. Kottmann et al., “Second harmonic generation microscopy reveals altered collagen microstructure in usual interstitial pneumonia versus healthy lung,” Respir. Res., 16 61 (2015). Google Scholar


L. Thiberville et al., “In vivo imaging of the bronchial wall microstructure using fibered confocal fluorescence microscopy,” Am. J. Respir. Crit. Care Med., 175 (1), 22 –31 (2007). AJCMED 1073-449X Google Scholar


L. Thiberville et al., “Human in vivo fluorescence microimaging of the alveolar ducts and sacs during bronchoscopy,” Eur. Respir. J., 33 (5), 974 –985 (2009). Google Scholar


B. A. Flusberg et al., “Fiber-optic fluorescence imaging,” Nat. Methods, 2 (12), 941 –950 (2005). 1548-7091 Google Scholar


L. Fu and M. Gu, “Fibre-optic nonlinear optical microscopy and endoscopy,” J. Microsc., 226 (Pt. 3), 195 –206 (2007). Google Scholar


H. Bao et al., “Fast handheld two-photon fluorescence microendoscope with a 475 microm x 475 microm field of view for in vivo imaging,” Opt. Lett., 33 (12), 1333 –1335 (2008). OPLEDP 0146-9592 Google Scholar


C. M. Brown et al., “In vivo imaging of unstained tissues using a compact and flexible multiphoton microendoscope,” J. Biomed. Opt., 17 (4), 040505 (2012). JBOPFO 1083-3668 Google Scholar


G. Ducourthial et al., “Development of a real-time flexible multiphoton microendoscope for label-free imaging in a live animal,” Sci. Rep., 5 18303 (2015). SRCEC3 2045-2322 Google Scholar


D. G. Ouzounov et al., “Dual modality endomicroscope with optical zoom capability,” Biomed. Opt. Express, 4 (9), 1494 –1503 (2013). BOEICL 2156-7085 Google Scholar


A. Lombardini et al., “High-resolution multimodal flexible coherent Raman endoscope,” Light Sci. Appl., 7 (1), 10 (2018). Google Scholar


W. Liang et al., “Nonlinear optical endomicroscopy for label-free functional histology in vivo,” Light Sci. Appl., 6 (11), e17082 (2017). Google Scholar


F. Akhoundi et al., “Compact fiber-based multi-photon endoscope working at 1700 nm,” Biomed. Opt. Express, 9 (5), 2326 –2335 (2018). BOEICL 2156-7085 Google Scholar


C. H. Hage et al., “Towards two-photon excited endogenous fluorescence lifetime imaging microendoscopy,” Biomed. Opt. Express, 9 (1), 142 –156 (2018). BOEICL 2156-7085 Google Scholar


W. Göbel et al., “Miniaturized two-photon microscope based on a flexible coherent fiber bundle and a gradient-index lens objective,” Opt. Lett., 29 (21), 2521 –2523 (2004). OPLEDP 0146-9592 Google Scholar


K. Thomsen et al., “Development of a confocal and two-photon endomicroscope—preliminary results of qualitative evaluation,” Med. Laser Appl., 25 (3), 166 –172 (2010). Google Scholar


M. Lelek et al., “Coherent femtosecond pulse shaping for the optimization of a non-linear micro-endoscope,” Opt. Express, 15 (16), 10154 –10162 (2007). OPEXFF 1094-4087 Google Scholar


C. Lefort et al., “Pulse compression and fiber delivery of 45 fs Fourier transform limited pulses at 830 nm,” Opt. Lett., 36 (2), 292 –294 (2011). OPLEDP 0146-9592 Google Scholar


A. Lukić et al., “Fiber probe for nonlinear imaging applications,” J. Biophotonics, 9 (1–2), 138 –143 (2016). Google Scholar


A. Lukic et al., “Endoscopic fiber probe for nonlinear spectroscopic imaging,” Optica, 4 (5), 496 –501 (2017). Google Scholar


M. Kalashyan et al., “Ultrashort pulse fiber delivery with optimized dispersion control by reflection grisms at 800 nm,” Opt. Express, 20 (23), 25624 –25635 (2012). OPEXFF 1094-4087 Google Scholar


C. Lefort et al., “Characterization, comparison, and choice of a commercial double-clad fiber for nonlinear endomicroscopy,” J. Biomed. Opt., 19 (7), 076005 (2014). JBOPFO 1083-3668 Google Scholar


D. A. Peyrot et al., “Development of a nonlinear fiber-optic spectrometer for human lung tissue exploration,” Biomed. Opt. Express, 3 (5), 840 –853 (2012). BOEICL 2156-7085 Google Scholar


C. Lefort et al., “Sub-30-fs pulse compression and pulse shaping at the output of a 2-m-long optical fiber in the near-infrared range,” J. Opt. Soc. Am. B, 31 (10), 2317 –2324 (2014). Google Scholar


B. Efron, “Better bootstrap confidence intervals,” J. Am. Stat. Assoc., 82 (397), 171 –185 (1987). Google Scholar


H. A. C. Wood et al., “High-resolution air-clad imaging fibers,” Opt. Lett., 43 (21), 5311 –5314 (2018). OPLEDP 0146-9592 Google Scholar


B. Morova et al., “Fabrication and characterization of large numerical aperture, high-resolution optical fiber bundles based on high-contrast pairs of soft glasses for fluorescence imaging,” Opt. Express, 27 (7), 9502 –9515 (2019). OPEXFF 1094-4087 Google Scholar


Anikitos Garofalakis received his PhD in biomedical physics in 2007 from the University of Crete. Currently, he is the head of the Research and Development (R&D) Department of the BiospaceLab Company. His work focuses on the development of optical imaging systems for biomedical applications. Prior to joining BiospaceLab, he pursued a career in academic research where he worked on the fusion of various imaging techniques in hybrid multimodal schemes.

Sergei G. Kruglik defended his PhD thesis in physics and mathematics in 1995 at the Institute of Physics, Minsk, Belarus. Following 15 years of research at this institute, he continued his endeavors at scientific laboratories throughout the world, in Japan, France, The Netherlands, and Canada. He is now a research engineer at Sorbonne University, Paris, France. His scientific interests are at the crossroads of laser physics, biophotonics, and life sciences, combining advanced laser techniques with optical spectroscopies and medical imaging.

Tigran Mansuryan received his PhD in laser physics in 2008 from the Physics Department of Yerevan State University. From 2009 to 2014, he was a postdoctoral researcher at XLIM, Limoges University. From 2014 to 2018, he was a R&D engineer at KAMAX Innovative System. He is now a researcher at XLIM, Limoges University. His research is focused on nonlinear optical methods of microscopy, fiber endoscopy, fiber nonlinearity, ultrafast optics, and femtosecond technologies.

André Gillibert, after a bachelor’s degree in pure mathematics, studied medicine at Rouen and Caen universities (2004 to 2017). He holds a master’s degree in methodology biostatistics and clinical research (Paris-Sud, 2016). Currently, he is a hospital assistant lecturer in the Department of Biostatistics of Rouen University Hospital.

Luc Thiberville is a professor of respiratory medicine and head of the Respiratory, Thoracic Oncology and Respiratory Care Medicine Department at Rouen University Hospital, France. His research fields include precancerous lesions of the bronchial tree and endoscopic detection and treatment of lung cancer. He introduced and is developing the technique of in vivo endomicroscopy of the respiratory tract, focusing on microimaging of peripheral lung nodules, proximal preneoplastic epithelium, bronchial lung cancers, and interstitial lung disease.

Frédéric Louradour is a professor at XLIM Institute (CNRS 7252), Limoges University, France. His preferred topics are nonlinear optics, femtosecond pulse characterization and shaping, and propagation in optical fibers applied to biomedical instrumentation. He studied endomicroscopy working with polarization or multiphoton or fluorescence lifetime imaging and using an image guide made of a high-resolution bundle of thousands of optical fibers or using a single dual-cladding fiber associated to a miniature fiber scanner.

Christine Vever-Bizet has a PhD in biophysics. She joined, in 2004, the Pierre and Marie Curie University, now known as Sorbonne University, in the Jean Perrin Laboratory, CNRS UMR 8237. Her scientific interests focus on photodynamic therapy, cellular and tissular autofluorescence spectroscopy and imaging, as well as confocal laser linear and nonlinear endomicroscopy

Genevieve Bourg-Heckly holds a PhD in physics. After twelve years of professional experience in laser and medical imaging industries, she joined the Pierre and Marie Curie University (currently known as Sorbonne-University, Paris, France) to develop projects in biomedical optics. Her research interests are in the field of noninvasive diagnostic methods based on optical spectroscopy and imaging. Over the past several years, her studies have mainly focused on the development of single-photon and multiphoton endomicroscopies.

CC BY: © 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.
Anikitos Garofalakis, Sergei G. Kruglik, Tigran Mansuryan, André Gillibert, Luc Thiberville, Frédéric Louradour, Christine Vever-Bizet, and Geneviève Bourg-Heckly "Characterization of a multicore fiber image guide for nonlinear endoscopic imaging using two-photon fluorescence and second-harmonic generation," Journal of Biomedical Optics 24(10), 106004 (23 October 2019).
Received: 20 January 2019; Accepted: 25 July 2019; Published: 23 October 2019

Cited by 8 scholarly publications.
Second-harmonic generation





Image quality


Back to Top