Fast epi-detected broadband multiplex CARS and SHG imaging of mouse skull cells

: We present a bimodal imaging system able to obtain epi-detected mutiplex coherent anti-Stokes Raman scattering (M-CARS) and second harmonic generation (SHG) signals coming from biological samples. We studied a fragment of mouse parietal bone and could detect broadband anti-Stokes and SHG responses originating from bone cells and collagen respectively. In addition we compared two post-processing methods to retrieve the imaginary part of the third-order nonlinear susceptibility related to the spontaneous Raman scattering.


Introduction
The use of non-destructive, non-invasive and label-free methods is particularly interesting to realize medical diagnosis of diseases or to develop tissue engineering and regenerative medicine [1]. In bone tissue engineering, the discrimination of different cellular pools is crucial prior to the application of any clinical treatment [2,3]. In this context spontaneous Raman [1,[3][4][5] and coherent anti-Stokes Raman scattering (CARS) [3,6,7] spectroscopies have proven that they could distinguish progenitor cells by detecting specific molecular vibrations.
To obtain a CARS signal, the difference between the optical frequencies of a pump and a Stokes laser beam, ω p − ω S , has to match the vibrational frequency of a molecule. The use of a third beam as a probe wave can then, by four-wave mixing, lead to the generation of anti-Stokes signal both emitted in the forward and the backward directions. Usually, the role of the probe is given to the pump beam so that the frequency of the anti-Stokes wave is ω aS = 2ω p − ω S . In most cases, CARS imaging is implemented in forward configuration where the anti-Stokes photons are collected after going through the whole thickness of the sample. The use of microscope objectives with high numerical aperture permits the phase matching condition to be always fulfilled. On the other hand, in epi-configuration the wave vector mismatch is more important [8]. Therefore, the signal coming from objects larger than the anti-Stokes wavelength is less or not collected due to destructive interferences. Consequently, epi-detection seems to be more suitable for living cells monitoring because the solvent does not participate in the generation of backward signal [9]. Furthermore, in this case, the implementation of the experimental setup is facilitated and the study of in vivo samples become possible. Epi-CARS analysis of biological objects has been demonstrated by measuring single-frequency anti-Stokes wave [10,11], but the simultaneous acquisition of different anti-Stokes wavelengths has not been achieved to our knowledge. To achieve multicolor epi-CARS, a supercontinuum radiation can be used as a broadband Stokes wave. In this manner, a broad anti-Stokes spectrum covering the fingerprint and the CH vibration regions can be instantaneously acquired. This system is well established in the forward configuration and called multiplex CARS (M-CARS) [12][13][14]. One publication introduces its implementation in epi-detection but it was only applied to the study of a mineral sample with long pixel dwell time [15]. In addition to the vibrational signal, second harmonic generation (SHG) from the pump beam, due to the presence of non-centrosymmetric components in the sample such as type I collagen fibrils, can also be acquired [16][17][18]. Thereby, CARS and SHG detection can provide rich and complementary molecular information.
We propose here a bimodal system, with epi-detected M-CARS and SHG, applicable to biological tissues. To demonstrate the ability of our setup, we study biopsies from mouse skull bone. We measure CARS spectra from 700 cm −1 to 3200 cm −1 with a spectral resolution of 2 cm −1 and a short pixel dwell time of 5 ms. Maps of typical vibrations of skull bone cells and SHG signal coming from collagen fibrils are obtained. We then compare qualitatively two algorithms that can extract the resonant signal from raw CARS data and we evaluate their respective post-processing time. Finally we compare our results with two-photon fluorescence images obtained on similar biopsies with different staining methods in order to confirm the presence of skull bone cells.

Materials and methods
The schematic of the epi-detected M-CARS and SHG system is presented in Fig.1. We use a laser source (Spark Lasers, Sirius series, 1064 nm, 60 ps, 110 kHz, linearly polarized) to generate the pump radiation. By using a half-wave plate and a Glan-Taylor polarizer (GTP), a part of the pump is directed into a photonic crystal fiber (PCF) to generate an infrared supercontinuum Stokes wave from 600 nm to 1650 nm. To compensate for the difference of optical length between the pump and the Stokes beams and, thus, to assure their temporal overlapping, a delay-line is set on the path of the pump. A long-pass filter at 1050 nm (Thorlabs, FEL1050) is used on the Stokes beam to suppress wavelengths that could be confused with anti-Stokes signals. The two beams are then spatially combined by mean of a notch dichroic mirror (Semrock, NFD01-1064-25x36) which can reflect the pump in the direction of propagation of the Stokes. Then, the beams travel together toward a dichroic short-pass mirror at 1000 nm (Thorlabs, DMSP1000R) that reflects them toward a microscope objective (Olympus, UPlanSApo 60x, N.A. = 1.2, water immersion) that focuses them onto the sample. The backscattered signal, collimated by the objective, is transmitted by the short-pass filter in the direction of a spectrometer (Horiba, LabRam Evolution). A notch filter at 1064 nm (Thorlabs, NF1064-44) is set to suppress the remaining pump signal. The different anti-Stokes wavelengths are then diffracted by a grating (Horiba, 520 29, 600 gr.mm −1 ) and spanned on the CCD camera (Horiba, Synapse).
The use of animals has received the approval of the local Ethical and Animal Care Committee (registration number 15-2015-15, MoNaPulse and CRAMpON). All animal cares and experimental procedures were conducted in conformity with the French legislation n • 2013-118 1 st February 2013 in accordance with European Community guidelines (directive 2010/63/UE for the Care and Use of Laboratory Animals). 12-week-old wild-type (no label), SwissTg(CAG-EGFP)1Osb (ubiquitous GFP expression) and B6.CgTg(Col1a1*2.3GFP)1Rowe/J (GFP in the osteogenic lineage) mice were dissected immediately following euthanasia. Cranium bone was removed from the parietal region with a dental drill as illustrated in Fig. 2. One label-free fragment was counterstained with 50 µM Rhodamine B (R6626, Sigma-Aldricht, λ em 627 nm) for 10 min. For two-photon fluorescence measurements we used an upright multiphoton microscope (Olympus, BX61WI/FV1200MPE) coupled with a tunable femtosecond Ti:Sapphire laser (Coherent, Chameleon Ultra II) via 25x immersion objective (Olympus, XLPLN25XWMP, N.A. = 1.05). The bregma and the lambda are both classical anatomical landmarks located on the midline of the outer surface of the skull, that are often used as reference points for stereotaxic surgery in this region. The bregma and the lambda are indeed the meeting point respectively between frontal and parietal bones, and between parietal and occipital bones. Here the cranium bone was extracted from the parietal region, down left from the lambda, as highlighted by the red circle.

Results and discussion
Using our CARS setup, we mapped an area of 260 µm × 260 µm with a 1 µm-step of the wild mouse parietal bone. The lateral resolution was 540 nm (FWHM). Because the grating and the CCD camera used could not collect the entire anti-Stokes spectrum and SHG response without loss of spectral resolution, we split the measurement into three sets of acquisitions: from 700 cm −1 to 1200 cm −1 (5 ms per spectrum); from 2500 cm −1 to 3200 cm −1 (5 ms per spectrum); and from 510 nm to 555 nm (1 ms per spectrum). For each point of the mapping, two accumulations were made to avoid random spikes in spectra due to cosmic rays.
In Fig. 3 are presented typical raw CARS spectra obtained at two different acquisition points in the fingerprint (700-1200 cm −1 ) and CH (2500-3200 cm −1 ) vibrational regions. In the fingerprint region of Figs. 3(a) and 3(b), typical dispersive profiles with maxima at 950 cm −1 and 1061 cm −1 demonstrate the presence of two vibrational modes. In the CH region the dispersive lines with maxima at 2869 cm −1 and 2932 cm −1 , only present in Fig. 3(b), attest to the presence of two other vibrational modes. The frequency of all these maxima does not exactly match the frequency of the corresponding spontaneous Raman band. As a matter of fact, due to the dispersive aspect of the CARS spectra, the frequency of the peaks is shifted in the low-energy direction. In Table 1, we list the vibrational modes for the different CARS peaks of Fig. 3, including the corresponding Raman shifts and associated molecule types. It appears that PO 3− 4 symmetric stretching and CO in plane stretching of CO 2− 3 characterize the presence of apatite crystals with phosphate and carbonate content [5,7,19,20]. CH 2 asymmetric stretching indicates the presence of lipids [21]. CH 3 symmetric stretching, which besides is coupled to overtone of asymmetric deformation [22], points out proteins [23].
We used two types of post-treatment methods to suppress the non-resonant background (NRB) from the CARS spectra. We first used an algorithm based on the Maximum Entropy Method (MEM) [24] in order to retrieve the imaginary part of the third order nonlinear susceptibility χ (3) involved in the four-wave mixing. It is this physical variable that carries the vibrational information, for it is directly comparable to the first order susceptibility χ (1) of the spontaneous Raman scattering. The second method utilizes a Savitzky-Golay filter [25] (second order polynomial fit with a 21-point window) to obtain the first order derivative of the CARS spectra with respect to the Raman shift. We consider that the NRB, varying slowly with the Raman shift, would present a weak participation to the derivative in comparison with the sharp variations of CARS peaks. Furthermore we computed in Fig. 4 the difference ∆ω between the positions of theoretical spontaneous Raman and CARS peaks (blue curve) according to the spectral half-width of the Raman peak (Γ) from 1 cm −1 to 200 cm −1 . We also computed the difference between the positions of the Raman peak and the derivative of the CARS peak (red curve). We can see that, compared to the position of a CARS peak, the position of its derivative is closer to the one of the Raman peak with a small increase of ∆ω with Γ. This indicates that we can use the derivative of a CARS peak to retrieve the position of the corresponding Raman line with good precision.  Fig. 4. Evolution of the difference ∆ω between the positions of a Raman peak and a CARS peak maximum (blue curve) and between a Raman peak and the derivative of a CARS peak maximum (red curve) according to Γ, the spectral half-width of the Raman peak.
The measured CARS spectra presented before have been processed using the two methods. The resulting spectra are plotted in Fig. 5. We observe that in both cases the dispersive profile of the CARS peaks has disappeared but the spectral features carrying the vibrational information have remained. This correction of the dispersive aspect is highlighted by the shifting of each CARS peak towards the Raman shifts presented in the literature. The two methods provided almost the same maxima positions: 961 cm −1 , 1070 cm −1 , 2878 cm −1 and 2950 cm −1 for PO 3− 4 symmetric, CO in plane of CO 2− 3 stretching modes and the CH 2 asymmetric and CH 3 symmetric stretching modes respectively. We also noted that the post-processing time of the derivative algorithm was shorter than that of the MEM algorithm: 25 ms per spectrum instead of 1 s per spectrum.
We then produced normalized intensity maps integrated over a bandwidth of 20 cm −1 for these four vibrational modes. To do so, we used both the spectra extracted with the MEM algorithm and the derivative spectra. By computing the standard-deviation and mean of each map, we found that the contrast was the same for the MEM and derivative algorithms. However a better signal to noise ratio of the maps was observed when using the derivative method. The intensity maps of PO 3− 4 , CO 2− 3 , CH 2 and CH 3 stretching modes obtained with the derivative method are displayed in Fig. 6, as well as the SHG intensity map. We merged three of those maps (PO 3− 4 at 961 cm −1 in red, CH 3 at 2950 cm −1 in blue, SHG in green) into a single image, as depicted Fig. 7(a). In this image, the red and blue clusters represent hydroxyapatite and proteins signal respectively. Thus these clusters would indicate the presence of bone cells embedded in a fibrous collagen matrix, the latter coinciding with the SHG signal in green.
To confirm this hypothesis, two-photon fluorescence intensity mapping of the same type of sample with different labels was realized. Thanks to the resulting maps we can observe the characteristic features of indeterminate cells labelled with GFP and Rhodamine (Figs. 7(b) and 7(c)); and osteogenic cells labelled with GFP ( Fig. 7(d)). In Fig. 7(b) the SHG signal originating from collagen is also plotted. This allows us to compare the previous M-CARS/SHG map with a reference measurement. The two images show a similarity between the size, shape and density of hydroxyapatite/proteins clusters on one hand ( Fig. 7(a)), and of GFP labeled cells on the other hand ( Fig. 7(b)). As well, the clusters and the cells both appear as embedded objects in the collagen matrix. We can thus affirm that the hydroxyapatite and proteins M-CARS signals originate from mouse bone cells. This highlights the relevance of epi-detected M-CARS for tracking biological cells in a complex environment.
Finally we calculated the number of cells from each of the five maps, i.e., from the two PO 3− 4 / CH 3 vibrational intensity maps and from the three fluorescence intensity maps. The resulting cell counting is plotted in the histogram of Fig. 7(e). At this stage, the variability in population density between the different modalities cannot be assigned to any specific chemico-biological assumption, because it would require a full biology protocol for making quantitative cell counting. However the difference of population density between PO 3− 4 and CH 3 modalities relates to the only-partial colocalization of hydroxyapatite crystals and proteins, which is also confirmed by the colorimetric variations in Fig. 7(a). This confirms the potential of M-CARS as a rich, label-free and chemically-selective imaging method for studying biological systems.

Conclusion
We have demonstrated that the epi-detection of M-CARS and SHG signals coming from mouse bone cells was possible at low average power and short pixel dwell time (5 ms and 1 ms, respectively). The repartition maps of PO 3− 4 symmetric stretching and CH 3 stretching vibrational modes could clearly reveal the presence of mineralized hydroxyapatite and proteins, respectively, and permitted the identification of bone cells. This study shows the potential of the multiplex configuration in the promising domain of CARS endoscopy [26]. Furthermore, we have proved that the use of the derivative of CARS spectra was as effective as the MEM for retrieving Raman lines, with high vibrational contrast. The derivative method also demonstrated shorter postprocessing time. This represents an interesting track towards real-time processing of multiplex CARS data. Finally, we have illustrated that, because the width of the CARS spectrum is crucial for calculating the derivative, broadband M-CARS measurements are highly compatible with this post-processing method.