Full-Field Optical Coherence Tomography as a Diagnosis Tool : Recent Progress with Multimodal Imaging

Full-field optical coherence tomography (FF-OCT) is a variant of OCT that is able to register 2D en face views of scattering samples at a given depth. Thanks to its superior resolution, it can quickly reveal information similar to histology without the need to physically section the sample. Sensitivity and specificity levels of diagnosis performed with FF-OCT are 80% to 95% of the equivalent histological diagnosis performances and could therefore benefit from improvement. Therefore, multimodal systems have been designed to increase the diagnostic performance of FF-OCT. In this paper, we will discuss which contrasts can be measured with such multimodal systems in the context of ex vivo biological tissue examination. We will particularly emphasize three multimodal combinations to measure the tissue mechanics, dynamics, and molecular content respectively.


Introduction
Over the past 25 years, optical coherence tomography (OCT) has rapidly developed and spread into many fields of biomedical research thanks to its ability to capture label-free structural contrast deep inside a tissue [1].Over the years, OCT has been developed into numerous different modalities that emphasize the 3D interplay between tissue microstructures and various specific complementary contrasts [2][3][4][5].
Full-field optical coherence tomography (FF-OCT) is a less familiar OCT-derived technology [6,7].FF-OCT differs from standard OCT by acquiring 2D en face images in a widefield configuration at a given depth in a tissue.To record a 3D tomogram, the sample is further scanned in depth, instead of acquiring series of A lines as in regular OCT.Such a configuration offers the possibility to use high numerical aperture (NA) microscope objectives, without limiting the exploration depth, and can therefore access submicrometric and subcellular resolution.
FF-OCT has been successfully applied for non-invasive label-free ex vivo pathology assessment in various tissues [8][9][10][11][12][13][14][15].Thanks to its isotropic micrometric resolution, FF-OCT can often be compared to histology slices.In contrast to histology, FF-OCT does not require time for sample fixation, slicing, or labeling, and an histology-like image can be produced within a few minutes after tissue extraction.It also helps to reduce fixation and slicing artifacts, as FF-OCT does not necessarily require physical slicing, and can image in depth, eventually far from the incision edges for a biopsy.
However, FF-OCT may suffer from its lack of specificity, and sometimes fails to equal the diagnostic capacity of gold standard histology.Yet, FF-OCT detection scores are above 80% for most pathologies in both sensitivity and specificity [14,16], and Manu Jain and colleagues could even report a 100% correct identification of tumors in ex vivo kidneys [17], obtaining similar results to histology, even though the tumor subtyping is still not 100%.In order to increase diagnosis scores, FF-OCT has therefore expanded to multimodal systems combining FF-OCT structural contrast with other pathology biomarkers inside a single microscope.The interest of such multimodal systems is to reduce the ambiguity in structure identification that can arise from sequential imaging with two different microscopes, and to reduce operation time.Some FF-OCT multimodal systems have been reviewed in the recently published Handbook of Full-Field Optical Coherence Microscopy, edited by Arnaud Dubois [16].We can non-exhaustively mention the multimodal combinations of FF-OCT with fluorescence [18], photothermal imaging [19], elastography [20], polarization [21], multispectral and hyperspectral imaging [22][23][24].Most of these multimodal systems have been inspired by similar systems developed in standard OCT, but they needed to be adapted and redesigned to fit FF-OCT's particular configuration.
In this paper, our goal was not to exhaustively describe these multimodal systems, or to systematically review all the potential fields in which multimodal FF-OCT systems can be of interest.Instead, we will start from the need to increase the FF-OCT's diagnostic ability up to equal that of gold standard histology to describe three different multimodal systems that we are currently developing or improving in our lab.After a brief description of FF-OCT, and the optical systems that we commonly use, we describe and emphasize the differences between FF-OCT and histology contrasts, to better apprehend what contrasts would be optimal to add to FF-OCT.Then, we will describe our recent development related to multimodal systems combining FF-OCT with elastography, with cellular metabolism, and with fluorescence, and how these multimodal systems could increase the diagnosis performance of FF-OCT.

Scanning OCT versus Full-Field OCT and Technological Description of the Setups Description
The main difference between classical OCT systems and FF-OCT lies in the absence of transverse scanning.In OCT, or in optical coherence microscopy (OCM), an axial line (A-line) is captured with a single point [1,25,26] or with a line detection [27,28], and further scanned in the transverse plane, to successively obtain 2D planes (B-scan) and 3D tomograms (C-scan).On the contrary, the full-field configuration captures a 2D image at a single depth with a 2D camera, and an axial scan enables the recording of a 3D tomogram.In both configurations, the optical section is mostly controlled by the temporal coherence length, given by the inverse of the spectral bandwidth.Scanning in one dimension compared to two dimensions has several major consequences:

•
Because it looks at a single depth at a time, FF-OCT can be performed with a small depth of field, and therefore with high numerical aperture (NA) objectives.This enables high transverse resolution to be obtained, below 1 µm for visible light.Other elegant OCT systems have greatly increased their transverse resolution [1,[29][30][31][32], however these techniques remain more complex and often have to sacrifice their temporal resolution.

•
Moreover, scanning in two directions imposes a trade-off between exposure time and the field of view.In scanning OCT, each pixel is only illuminated for a few microseconds, while in FF-OCT, the entire field of view is usually evenly illuminated during a few milliseconds.It can make a difference in the light dose delivered to the sample, and the most rapid speckle fluctuations are often averaged out [33].

•
The time required to get information from the sample is much longer in FF-OCT, since it is limited by the camera frame rate, which is hardly below 1 ms, while some scanning OCT systems can now acquire lines at more than 100,000 lines per second [34].The consequence is that FF-OCT is much more sensitive to sample motion, and is difficult to perform in vivo, other than on anesthetized animal or by pressing the sample against an imaging window [35].
• A last important difference when looking at en face images either performed in scanning OCT or FF-OCT is that all the pixels have been acquired at different times in OCT, while they are acquired at the same time in FF-OCT (if the camera is operated in global shutter mode), which can introduce some artifacts when looking at moving objects.On the other hand, all axial pixels are acquired simultaneously in OCT, while they are largely separated in time for FF-OCT.An important consequence is that the optical phase can be directly retrieved in OCT by subtracting two adjacent axial pixels [36] while FF-OCT requires multiple phase-shifted measurements to retrieve the optical phase.
In practice, we mainly use two different FF-OCT layouts presented in Figure 1.The first setup is the LLTech commercial system [37], based on a Linnik interferometer using two identical 10X, 0.3 NA microscope objectives, and a thermal source to produce an isotropic resolution of 1 µm.We also developed custom-built versions of this system, exhibiting higher mechanical stability, or different cameras for example, but the general layout stays identical.We used these systems to add elastography [20,38] and motility contrasts [39,40].The second setup is a multimodal setup combining FF-OCT and structured illumination fluorescence microscopy (SIM) [41], in a relatively similar layout to the one developed by H.Makhlouf and colleagues [42].To enhance the photon collection, we used higher NA objectives, either 40X, 0.8 NA in water, or 30X, 1.05 NA in silicone oil, achieving transverse resolution below 500 nm, and depth of focus between 1 and 2 µm.Interestingly, in this system, the optical sectioning ability comes from the lateral spatial coherence [43,44] controlled by the depth of focus, and therefore the objective NA rather than the temporal coherence, meaning that simple LEDs can be used without sacrificing the axial resolution.Both systems are illuminated by spatially incoherent sources, which eliminate the cross talks [45], and preserve the transverse resolution from degradation due to low order aberrations [46].The illumination sources are mounted in a Kohler illumination to provide an illumination as homogeneous as possible.
Both systems are now operated with a new CMOS camera developed by Adimec during an ERC project (CAReIOCA) led by LLTech [47].This camera can handle a Full Well Capacity of about 2Me − on 2 million pixels at frame rates up to 700 Hz, which typically enhanced our signal-to-noise ratio by a three-fold factor, and acquisition speed by a five-fold factor, with respect to previously used cameras (PhotonFocus, MV1-D1024E-160CL, Lachen, Switzerland).

FF-OCT and Histology Contrasts
In order to further compare histology and FF-OCT, we describe in this section histology contrast and procedures.
It is commonly said that Marcello Malpighi is the precursor of modern histology, as he was the first scientist to carefully describe micro-anatomic features of plants and animal tissues by the end of the 17th century with the first microscopes [48,49].Since then, histology has continuously benefited from advances in chemical labeling and fixatives, and has benefited from a better understanding of tissue nanostructures and organization thanks to electron microscopy [50].It now benefits from an extensive literature, and from digital atlases [51][52][53] enabling direct comparisons with normal or pathological slices of most tissues.However, the basic microscopy principles have mainly remained unchanged since the beginning of histology.They consist of imaging a thin slice of a sample with a simple transmission optical microscope.Since light absorption and scattering by biological structures exhibit a low contrast, histology commonly uses chemical agents to reveal specific contrasts of interest [54][55][56].The most common staining uses two associated stains, Hematoxylin and Eosin (H & E), that respectively label chromatin and chromosomes in blue, and cytoplasm in red.Hemoglobin usually appears orange, and collagen fibers appear pink.This staining is particularly versatile, as it can be used with a variety of fixatives.A wide range of different staining can be applied to a tissue slice.We can non-exhaustively evoke Golgi staining that labels 1% to 2% of neurons with silver that absorbs light to appear dark, Nissl staining that labels nucleic acids, or tetraoxide osmium that labels unsaturated lipids and lipoproteins.If many subcellular structures can be selectively labeled, it is often difficult to combine different stains in a single slide.However, using specific labels imposes a prior knowledge of the tissue composition, and of the pathology, and some unlabeled biomarkers could be missed.
Moreover, histology requires tissue fixation, performed either chemically or by freezing the sample, and physical slicing of the tissue, usually performed by a microtome.The fixation first goal is to harden the sample to be able to cut it in thin sections.However, fixation can introduce strong artifacts, as it often removes lipidic structures, and can increase the tissue volume.The physical slicing can induce physical damage, notably in neuronal tissues, in which many axial neurites are usually cut, leading to neuronal death.Finally, these two steps are usually time consuming, and can lead to the physical destruction of the samples if not adequately performed.
In the case of histopathology, the tissue is often fixated by replacing water by paraffin, and is then sliced in 4 to 5 µm slices.For each slice, the dyes have then to diffuse through the paraffin scaffold.In total, these processes can last for a few days.
In the case of intra-operative histology, the time allowed to the pathologist is usually not more than 15 to 30 min.To meet these time requirements, the sample is frozen, and directly sectioned and colored.However, many tissues, such as fat, brain, or retina, are not well suited for this process, and the number of available labeling is limited.
On the other hand, OCT was only invented 25 years ago, and its full-field version is 18 years old, so comparatively little literature and atlases can be found [57].However, OCT and FF-OCT provide a label-free, and optically-sectioned imaging.This means that imaging can be performed right after tissue extraction, which can provide a very fast diagnosis, and could even ideally provide an intra-operative diagnosis to efficiently define a tumor margin for example.
The main advantage of OCT is its virtual sectioning ability, enabling imaging at a few hundred microns below the tissue surface, hence reasonably far from the sectioned surface no matter which section is required [58].To completely remove slicing artifacts and reach large depths, OCT systems have successfully been coupled to a vibratome, in a procedure where an entire volume is imaged before a physical slice smaller than the imaged volume is performed [59].
Additionally, OCT systems capture the light naturally backscattered by the sample, without using any label.Working label-free can be advantageous, as staining artifacts can be avoided, and as it can importantly reduce the time prior to imaging.On the other hand, the recorded signal often lacks specificity.The obtained contrast is complicated to predict, as it is described by multi-scale complex theories, notably combining Rayleigh and Mie theories, depending on the scatterer size [60].However, for a given tissue, the obtained contrast is highly reproducible, so the diagnosis performance could be increased only by simply better characterizing the typical contrast obtained on each tissue, and each condition.
Quite generally, we can nevertheless emphasize the strong ability of FF-OCT to detect high refractive index (high lipid or high fat content) fibers, such as collagen fibers [8], myelin fibers [61], and even some large unmyelinated axons in the eye [62].Moreover, in terms of contrast, looking at the backscattered light can differ substantially from transmission microscopes in thick samples.However, one can assume that for thin slices, the contrast in reflection is mainly the inverted contrast for scattering structures and the same contrast from absorbing structures compared to a transmission microscope.
The label-free ability of OCT has a final practical advantage, since the imaging can be performed prior to any tissue modification.It means that the same tissue can be further used for another procedure, including a standard histology procedure [8], or fluorescence labeling.It is usually preferred to perform the OCT imaging before any sample modification, as OCT signal alteration can occur after fixation, depending on the tissue content [62] (especially since fixatives often remove lipidic structures of the cells).Keeping this in mind, it is nevertheless possible to perform histology-like measurements in fixed tissues [58].
Scanning OCT has also been used successfully to perform histology-like measurements [58,[63][64][65], and to detect tumor regions [66,67] through changes in the collagen content [68], the tissue microangiogram and backscattering [69], or in the tissue fractal dimension [70].However, we believe that the high resolution of FF-OCT can significantly improve the diagnostic abilities of OCT [8,17].With higher numerical apertures, the light collection is improved which helps to reveal more structures and increase the signal-to-noise ratio for the calculation of backscaterring coefficients or fractal dimensions [40,71].Signal from cell contours and from intracellular structures can also be detected.This will be described in Section 5 on Dynamic FF-OCT.Finally, the en face view of the sample is quite natural for histologists, as it looks similar to histology slides.FF-OCT is also faster than scanning OCT to record a single en face image.
Figure 2 shows a comparison between FF-OCT and a regular H & E histology in a breast sample, the organization of which is altered by a ductal carcinoma in situ (DCIS).FF-OCT notably reveals normal adipocytes, normal and abnormal collagen fibers (based on their scattering strength), and enlarged ducts.These indicators are characteristics of DCIS, that could be identified with 90% sensitivity and 75 % specificity in this case [14].

FF-OCT and Biomechanics
Since the early beginning of medicine, palpation has been used to qualitatively assess abnormal tissue stiffness associated with several diseases.Palpation-based diagnosis of heart, abdomen, and wounds diseases or infections has been reported in many ancient cultures, such as Greek, Egyptian, or Chinese medicine [72,73].The mechanical properties of tissues and cells are mostly related to their structure and function, and changes in a tissue's biomechanical properties can reveal its pathological state [74].Notably, biomechanical changes can include the stiffness increase due to an extracellular matrix accumulation, including collagen accumulation, associated with liver fibrosis, and ultimately cirrhosis [75,76].Moreover, tumor development significantly and progressively impacts the tissue biomechanics [77].More importantly, the tumor stiffness actively increases its progression, and thus its malignancy, partly through a mechanically-induced molecular circuit, leading to the reduction in the number of tumor suppressors [78].Therefore, not only tissue stiffness, but also stiffness anisotropy can be used as relevant diagnosis parameters.
In the 1980s, the thorough and quantitative measurement of tissue elasticity was developed, giving access to stiffness in the approximation of purely elastic deformation of tissues.Elastography has been introduced in biomedical imaging with ultrasound and MRI elastography [79,80], and has been proven to be a precious diagnosis biomarker for physicians since then [81,82].In this paper, we will only discuss tissue biomechanics in comparison to palpation-based diagnosis, and therefore at macroscopic and mesoscopic scales.However, we can add that cellular and molecular biomechanics can also be significant biomarkers of diseases, but they stay difficult to measure in tissues.To our knowledge, the most promising technique to study microscopic biomechanics inside a tissue is fluorescent FRET sensors [83].
In OCT, mechanical contrast has been introduced by Schmitt in 1998 [84].Since this pioneering work, many methods have been developed with success such as 2D, 3D, static or dynamic elastography [84][85][86][87][88].In FF-OCT, our group adapted two similar approaches to measure mechanical contrast.The first approach is a static method based on a finite element 3D correlation algorithm that gives access to the 3D local strain tensor inside the samples [20], giving access to the local shear modulus, which is itself proportional to the tissue elasticity, or stiffness for soft tissues.The second approach is a quantitative dynamic method based on shear wave speed imaging, inspired by similar experiments performed in ultrasounds developed in our lab.Using a fast FF-OCT system working at up to 30,000 frames per second, we could measure the local shear wave speed which is directly related to the local shear modulus [38].
In this section, we present a brief description of both methods and some results obtained on biological samples.

Description of the Different Mechanical Parameters and Approximations
Before describing how it is possible to optically measure tissue biomechanics, we briefly summarize here the numerous different mechanical parameters that can be found in the extensive literature on tissue biomechanics.We will emphasize the relationship between these parameters, and their relevance for diagnosis.
Most studies in MRI, as well as in ultrasound or optics, focus on the elastic properties of tissues, meaning the reversible linear deformation ε caused by a small mechanical stress σ.In the elastic regime, and neglecting viscosity, the generalized Hooke's law gives: where C is a fourth order tensor, named stiffness tensor, that depends on the mechanical parameters of the sample.All the other parameters that we will describe are elements of this fourth order tensor, and are usually obtained when the dimensions of the generalized problem are reduced.For example, Young's modulus E corresponds to the linear coefficient that relates a purely axial stress, and its associated deformation along the stress axis.In this case, the problem is reduced to a single dimension: We can add that Young's modulus is qualitatively the main measurement that is performed with manual palpation.Basically, the first information that one can get with manual palpation is how the tissue resists to an axial deformation or stress, which is given by the Young's modulus.
From Equation (1), it is easy to realize that a given sample is not only deformed in the direction of the stress, but can also be deformed in the transverse direction, with a preferential axis for anisotropic material.We will describe in Section 4.2 an example metric to quantify the local anisotropy of a sample.
In addition to these parameters, many mechanical measurements rely on shear wave speed measurement that we will describe more carefully in Section 4.3.The mechanical parameters on which the propagation of mechanical waves depend are simply defined as a result of the wave equation in a medium of a given stiffness.Inside biological media, two waves can propagate [89]: A longitudinal wave that corresponds to compression and traction of the medium with a mechanical displacement parallel to the direction of propagation.Additionally, shear waves can propagate, which correspond to perpendicular mechanical displacements with respect to the direction of propagation.The two waves propagate at speeds c L and c shear respectively that are controlled by the mechanical parameters of the medium.We can therefore define the bulk modulus K, and the shear modulus µ with: where ρ is the medium local density.These parameters are related to Young's modulus [89]: Finally, for soft tissues, which represent most biological tissues, except bones, the bulk modulus is much larger than the shear modulus.Bulk moduli of biological samples range from 1 GPa to 10 GPa, while shear moduli range from 0.1 kPa in the brain to 0.1 GPa in epidermis and cartilage, and up to 10 GPa in bones [90].When considering Equation ( 5) for soft tissues, the Young's modulus becomes dependent on the shear modulus only: With such approximations, we have shown that it was possible to locally measure the Young's modulus by measuring the propagation speed of shear waves.

Static Elastography Based on Digital Volume Correlation (DVC)
The principle of the static approach is to register a volumetric image before and after mechanical solicitation of the sample [20].From these two sets of images, we can estimate the 3D local displacement field, u, inside the sample using a finite element 3D digital image correlation algorithm developed by S. Roux and his team [91,92].From the optically measured 3D displacement field, we can calculate the 3D strain tensor ε: However, biological samples can be very complex and mechanically heterogeneous.This complexity usually prevents direct access to the local stress tensor and so, in most cases, methods based on digital image correlation are limited to strain contrast and do not provide quantitative stiffness information.This is the main limitation of this type of approach but already the strain information recreates a sort of "virtual local palpation" that can highlight mechanical heterogeneity at the micrometer scale given by FF-OCT or OCT .From the local strain tensor, it is possible to extract different parameters that give access to mechanical information such as strain magnitude or strain anisotropy.

Strain magnitude
In order to quantify the strain magnitude, we chose to use an equivalent Von Mises strain defined as: This parameter is independent of the base in which the strain tensor is expressed and allows quantification of the strain amplitude in the case of incompressible media such as biological samples.To illustrate this, we present strain amplitude maps calculated for a human breast tissue in Figure 3A,B, which highlight mechanical differences.

Strain anisotropy
The full 3D strain tensor also allows calculation of eigenvalues and eigenvectors of the strain tensor that gives access to information on the mechanical anisotropy of the sample.One should keep in mind that, since the compression of the sample is generally performed along the optical axis, the local stress field is inherently anisotropic.In order to display and to quantify strain anisotropy, we mainly use two parameters.The first parameter is the principal strain direction.For each pixel, the eigenvector corresponding to the largest eigenvalue of the strain tensor is calculated.This vector gives the principal strain direction.For example, Figure 3C is a cross-sectional slice from the volumetric FF-OCT image of a rat heart.In the FF-OCT image, it is difficult to evaluate the fiber direction, and the two different fiber orientations are not visible.However, if the projection angle of the principal strain direction vector on the plane perpendicular to the compression (see Figure 3D) is plotted, it is clear that there are two different fiber orientations.
Another parameter used to quantify the strain anisotropy is the fractional anisotropy FA.
In a similar manner to classical diffusion tensor imaging (DTI) in MRI [93,94] or ultrasound [95], we define the fractional anisotropy as: where ε 1 , ε 2 and ε 3 are the strain tensor eigenvalues and ε is the mean of the absolute value of the strain tensor eigenvalues.Using this parameter, it is possible to differentiate isotropic from anisotropic strain.In Figure 3E,G, we show that computing the fractional anisotropy in isotropic and anisotropic polyvinyl alcohol (PVA) polymer gels [96] enables discrimination of the two samples, in a case where FF-OCT fails to detect a difference.
In parallel to our work on 3D digital image correlation, we adapt the "tactile" approaches developed by Kennedy et al [97] for standard OCT.The principle is to use a thin layer of polymer deposited on the imaging windows to mimic manual palpation at the micrometer scale [98].Knowing the geometry and the mechanical properties of the polymer layer, it is possible to directly measure the stress at the surface of the sample by measuring the displacement at the surface of the polymer layer.With FF-OCT, the 2D en face displacement can be measured by taking the difference between two phase images taken at the surface of the polymer before and after the compression.This method allows a fast stress measurement at the surface of the sample simultaneously with FF-OCT imaging.This method is particularly interesting in the case of micro-biopsy where different tissues are usually juxtaposed.Figure 3I,J show a measurement on a biopsy of bovine kidney.At the micrometer scale, the stress map shows different mechanical behaviors for the muscle fibers and the kidney at the same stress level.The higher pressure at the muscle fiber level shows that this part of the sample is stiffer than the kidney area because the polymer is more deformed by this tissue.According to an experiment that consists of the observation of the tissue behavior during the compression, the pressure is negative at the kidney level because the tissue tends to move sideways.The negative pressure where there is no tissue can be explained by the fact that the silicone is incompressible and needs to expand somewhere to conserve its total volume.Finally, this picture shows the high resolution of this method because each muscle fiber with 20 µm diameter can be distinguished by comparing the FF-OCT image and the stress map.FF-OCT static elastography measurements can be performed with a standard FF-OCT system, such as the LLTech commercial system [37], or can be associated with other multimodal FF-OCT systems to diversify the types of contrast available.However, this technique can be computationally intensive, and requires the calculation of two 3D stacks and a physical compression of the sample.In total, the acquisition time can last for one to a few minutes, and the processing time for several hours.FF-OCT static elastography measurements only require a static FF-OCT signal to correlate, so penetration depth is only limited by FF-OCT capabilities.

Transient Elastography Based on Shear Wave Imaging
Transient elastography based on shear wave imaging is the second approach used to combine elastography with FF-OCT [38].Transient elastography based on shear wave imaging was introduced in OCT by Razani et al. [99].As shown in Equation ( 6), a direct measurement of the local shear wave speed gives direct access to a quantitative value of Young's modulus.Dynamic elastography methods involve two steps: first, a shear wave is locally sent onto the sample; then, the shear wave propagation is recorded.We were inspired by a technique named supersonic shear wave imaging, partly developed in our lab, which uses the ultrasound radiation force to create short shear impulses and ultrafast ultrasound imaging to image the shear wave propagation [100,101].We adapted this method to FF-OCT by generating shear waves using an ultrasound system and by recording the generated wave using a FF-OCT system operated with a high speed camera that can run at up to 30 kHz for a 512 × 512 pixel field of view [38].Shear waves were induced using a programmable ultrasound scanner (Aixplorer R , Supersonic Imaging, Aix-en-Provence, France), and a 15-MHz ultrasound linear array (128 elements, pitch 100 µm, elevation focus at 12 mm, Vermon, France), as described in A. Nahas et al. [38].This system focuses ultrasonic beams for tens of microseconds in the sample, which create a local radiation force, localized at the ultrasound focal spot and oriented along the beam axis.In response to this force, the tissue relaxation generates a shear wave, polarized along the ultrasonic beam and propagating transversely [90,102].To increase the amplitude of the radiated shear wave, an original solution consisting of successively focusing the ultrasonic beams at different depths along the ultrasonic path is performed [38,102].If the shear source is moved at a higher speed than the speed of the radiated shear wave, the source induces constructive interference along a Mach cone between all the generated shear waves, and the resulting wave can propagate with minimized diffraction over a large area.For short ultrasonic bursts (shorter than one millisecond), the spatial distribution of large tissue displacements initially corresponds to the acoustic intensity distribution, i.e., a cylinder with a lateral extension of 200 µm and an axial extension of about 1 mm in this configuration.The shear wave location can further be swept transversely by electronically focusing the ultrasound beam at different lateral locations, without any mechanical translation involved.The pressure levels involved are within the limits set by the Food and Drug Administration (FDA) for the application of diagnostic ultrasound in several organs.
By taking advantage of the interferometric sensitivity of FF-OCT to small axial displacement, we can further detect the propagation of the shear wave with high sensitivity.Unfortunately, high speed cameras do not offer large full well capacities (and it would be hard to generate hundreds of thousands of photoelectrons per pixel at 30 KHz), so our axial sensitivity is about 10 times lower than in our standard systems.Nevertheless, we could successfully follow shear waves in a sample and calculate their local speeds, by computing the temporal cross-correlation between neighboring pixels.One such propagation in an ex vivo brain tissue is shown in Figure 4. To better capture the propagation in time, we did not modulate the reference piezo, but only recorded direct images.However, we get rid of the incoherent intensity captured by the camera by subtracting the direct images to a few initial frames captured before the generation of the shear wave.The signal of each pixel is therefore directly proportional to the cosine of the phase: where a δ is linked to the backscattered wave amplitude interfering with the reference signal; δ is the path length difference at the considered pixel; λ 0 the central wavelength; ϕ(t) the relative phase that depends on time; and ϕ ini the initial phase.In Figure 4, we can see the wavefront of the shear wave crossing the field of view.With this method, we measured a shear wave speed of 0, 64 m•s −1 , corresponding to a Young's modulus of 1.25 kPa, which agrees well with the literature [103,104].We should add that although this dynamic measurement seems quite promising, it is performed here in an ideal case, since the brain is one of the softest tissues of an organism, and is quite reflective.In other cases, dynamic elastography FF-OCT can be challenging mainly due to current technological limitations.Indeed, we could not find any high sensitivity and high speed camera with parameters which would be sufficient to accurately follow shear waves in most tissues.However, since the shear wave source is controlled here, it should be possible to average several measurements to increase the sensitivity, or to perform a stroboscopic measurement by shifting the time delay between the shear wave generation and therefore measure the wave propagation with lower speed cameras.Ultimately, a shear speed can be calculated from only three measurements, so mechanical maps can be computed almost at video rate with dynamic elastography FF-OCT, and has been demonstrated at several thousand frames per second here.Here again, the penetration depth is only limited by FF-OCT characteristics.However, due to the low sensitivity provided by fast cameras, the penetration depth is, in practice, limited to 50 µm.To conclude this section, we have shown that FF-OCT can be coupled with elastography to measure the local mechanical parameters of a tissue.In terms of coupling with biomechanical contrast, FF-OCT and standard OCT measurements are not very different, since they rely on similar techniques to generate the mechanical contrast.Although FF-OCT-based mechanical measurements should theoretically have an isotropic 1 µm resolution, we could never measure any relevant differences at this scale.Additionally, the mechanical models used to link local deformations with mechanical stiffness are mesoscopic models, whose validity can become questionable when decreasing the resolution below 10 µm.
We should add that as the deformation is applied along the optical axis, it seems more natural to look at its propagation in the transverse plane, especially when looking at shear waves.Therefore, FF-OCT, or any other en face OCT techniques should be well adapted for such dynamic mechanical measurements.However, mainly due to current technological limitations, these measurements are still difficult to perform in FF-OCT.
Finally, a promising technique, named passive elastography [105,106], allows one to measure dynamic mechanical parameters in the case where the sampling rate that is lower than the shear wave propagation has recently been demonstrated in spectral domain OCT in our lab [107].Passive elastography is well adapted to FF-OCT measurements, and is currently under development in our group.

Dynamic FF-OCT
In this section, we will briefly review D-FF-OCT, including details of available contrasts and structures that can be revealed, and we present a few potential applications.
We developed dynamic full-field OCT (D-FF-OCT), inspired by the work on quasi elastic dynamic light scattering and other dynamic measurements to measure biological activity [33,[108][109][110][111][112][113].We therefore wondered whether the intracellular movements and vibrations can be important enough to modulate the optical path length of our interferometer by themselves, even without reference piezo modulation.Interestingly, we found out that in fresh biological samples, calculating the standard deviation of the interferometric signal over time gives a substantially different contrast in comparison with the regular FF-OCT signal.We discovered later that this dynamic signal was related to cellular metabolism [39].Interestingly, D-FF-OCT can access motility contrast in 3D, taking advantage of the FF-OCT low coherence to only select the phase variations at a given plane.
Now, what we generally name D-FF-OCT is a set of mathematical techniques to analyze the temporal fluctuations of the FF-OCT signal.Trying to understand and quantify these fluctuations, we were inspired and influenced by similar measurements performed with more standard OCT [4,[114][115][116], other interferometric techniques [33,111], or general dynamic techniques [117].
The new contrast revealed by D-FF-OCT originates from a different interpretation of the temporal fluctuations of the signal.In regular FF-OCT, the temporal fluctuations are averaged out to increase the signal-to-noise ratio of strongly backscattering structures.On the contrary, in the first versions of D-FF-OCT [39], we emphasized the amplitude of fluctuations of the signal by calculating the standard deviation, as described previously.This amplitude of fluctuations not only depends on the backscattering strength, but also on the scatterers dynamics inside a given voxel.An interesting layout of D-FF-OCT is to remove the piezo modulation, since the incoherent light varies less, and more slowly than the interference term.In such a configuration, the backscattering amplitude and the phase fluctuation amplitude are coupled in the final image.However, we found out that D-FF-OCT images calculated from four phase signals are often less contrasted, so we often stop the piezo movement and only acquire direct images.
To extract more information from intracellular dynamics, we now try to emphasize not only the amplitude of fluctuations, but also the timescales of variations.We measured the interferometric signal autocorrelation to extract two characteristic times, and could measure the metabolic-related effects of a drug that depolymerizes F-actin [40].Similarly, we also emphasize fluctuation timescales by calculating the Fourier transform of the time signal at a given pixel, or set of pixels.Then, we integrate the signals between different frequency bands to measure the dynamic signal originating from different timescales.We displayed the values by using an HSV color space, with the central frequency as the hue, the spectral width as the saturation and the amplitude of fluctuations as the value for each pixel of an image.This method is illustrated in Figure 5B,E, and enables segregation of cell populations based not only on their shapes but also on their fluctuation dynamics.We can distinguish cells with highly dynamic membranes, such as red blood cells [118] or lymphocytes [119] from slowly varying cells.
Finally, both our group and H. Ammari et al. [120] also used singular value decomposition (SVD) to better separate large structures that are slowly varying from smaller structures that are varying more rapidly.SVD can also eliminate slow movements from all the samples, similarly to what was obtained in ultrasound [121].
Although the location of cell bodies can often be inferred from the FF-OCT image in regions with low signal, D-FF-OCT can distinguish between a cell and a local oedema (which appears black in D-FF-OCT) for example, and gives a better precision of the cell contour.For example, in Figure 6 and in the cortex region of coronal brain slices, the position of neuron somas can be inferred from the FF-OCT image, but the dynamic image can better reveal individual neuron shapes and sizes, from which the neuron type might be inferred (pyramidal or stellate neurons).Ultimately, the identification of cortical layers can be performed in situ.In addition, D-FF-OCT may play an important role as a diagnostic tool as it can detect nuclei in large cells (Figures 5E,F and 6).Nuclear size and shape is an important biomarker of cancer, with higher nuclear to cytoplasmic ratios in tumor cells [122,123].The nuclei mainly appear black in D-FF-OCT, even though the inner dynamics of the nuclei are significant and usually faster than cytoplasm dynamics [110].Our hypothesis for the dark appearance of the nuclei is that nuclei are more densely packed than cytoplasm, and appear somehow homogeneous at a 1 µm resolution.Inside a homogenous nucleus, there should be no backscattering except on its contour and no dynamics to image, which can explain why nuclei appear black in both FF-OCT and D-FF-OCT.Interestingly, imaging cells with higher NA objectives, and with a transverse resolution below 500 nm, enables us to recover some dynamic signal inside nuclei (See Figure 5F).Another interesting D-FF-OCT feature for diagnosis is its ability to detect red blood cells and immune cells (See Figure 5B,D,E).Indeed, in ex vivo tissues, the red blood cells are trapped inside the capillaries, and can be imaged even at low frame rates.Due to their strong membrane fluctuations [118,124], they produce an intense and fast dynamic signal (See Figure 5E,G) on their edges.Similarly, immune cells are known for their intense motility as they are constantly exploring their environment and changing their shapes, as their protrusions are looking for abnormal markers to suppress [119].Interestingly, the nature of the immune cell might be inferred from their motility coefficient [119].Firstly, D-FF-OCT can therefore detect immune cells infiltration inside the tumors, which seems to be correlated with the tumor grade [125,126].Secondly, FF-OCT can usually detect capillaries, but cannot make the difference between a blood vessel and a lymphatic vessel.Thanks to the additional information from D-FF-OCT, and by detecting either red blood cells, or lymphocytes in a vessel, a more accurate angiogram and lymphogram could be established in cancerous tissues.Interestingly, correlations between an increase of the blood vessel density and tumor grades have been well established for a long time ( [127,128] among many others).Similar results with the lymphangiogram have been obtained [129].Quite generally, denser angiograms and lymphangiograms increase the risk of metastasis.D-FF-OCT is promising, when compared with digital holography techniques [33], dynamic light scattering techniques [110], or quantitative phase imaging [124,130,131], as it provides optical sectioning, and can access a motility-based contrast in 3D in scattering tissues with a subcellular resolution.Additionally, compared with dynamic OCT techniques [116,132], we take advantage of the higher spatial resolution to capture subcellular variations, which further enable the cells to be revealed.To our knowledge, the contrast revealed by D-FF-OCT is therefore quite unique.Similar dynamic measurements in OCT mainly focus on the blood flow velocity [116,133,134], even with high transverse resolution [135] similar to ours, but, to our knowledge, they never directly show dynamic signals coming from intracellular dynamics.We are wondering whether it comes from a loss in the signal collection due to the use of lower NA illumination, from a difference in the algorithms used (if only blood flow was expected), or from a difference in the exposure time per pixel.
In terms of multimodality, D-FF-OCT is attractive, as it does not require any modification from a regular FF-OCT.It can be associated with any other dual modality system to provide three different contrasts [41], as will be pointed out in the next section.D-FF-OCT is highly complementary to standard FF-OCT especially in fibrous environments, as fibers are highly backscattering structures and often hide cells entangled inside the matrix.However, fibers are usually still, while cells have intracellular activity that produces a highly dynamic signal.Hence, D-FF-OCT can easily reveal these cells, and suppress the stationary signal from the fibers.Nevertheless, one should keep in mind that D-FF-OCT significantly lowers the acquisition rates of FF-OCT microscopy.Indeed, even though each image is acquired at the camera frame rate, significant cellular dynamics occur in the second timescale, which is not very compressible because it is mainly sample dependent.In practice, our fastest D-FF-OCT measurements were at one to two seconds per plane that therefore expand to several minutes for a 3D volume.As for penetration depth, D-FF-OCT shows similar features to FF-OCT, except that the D-FF-OCT signal is often lost before the FF-OCT signal when going in depth, mainly because D-FF-OCT is generally weaker.Nevertheless, we succeeded in imaging D-FF-OCT signals at 100 µm in brain slices, or even through the entire retina, going through up to 200 µm.

Full-Field OCT and Fluorescence
In principle, usual histology staining and fluorescence techniques are quite similar, except that fluorescence techniques use fluorescence probes instead of absorbing probes in histology.Some fluorescence techniques have been part of the histology toolkit for decades [136], such as immunofluorescence, which aims at detecting the interaction of a labeled antibody with a structure of interest [137], or fluorescence in situ hybridization (FISH) ( [138,139], Chapter 6), which aims at detecting the interaction of a labeled nucleotide sequence with its specific complementary sequence in a cell.More recently, confocal microscopy [140] and two photon microscopy [141] have been used in order to detect fluorescent probes in depth in tissues.Interestingly, thanks to the optical sectioning provided by these techniques, it has become possible to follow some fluorescent structures in ex vivo tissues, and in vivo [142,143].
However, fluorescence techniques still require fluorescent labels, that need to be injected in the tissue.It can take a long time before the label diffuses in depth inside a tissue.Moreover, and similarly to classic histology stains, these labels only tag a few structures of interest, and it might be hard to understand the organization of these structures inside their microenvironment.Multiple labeling can be used to tag several structures, but it increases the complexity, and the cost of such experiments.
We can also add that usually, confocal and two photon microscopy can only access shallow depth inside a tissue, down to a few hundred µm at best.OCT techniques on the other hand can access deeper regions down to one millimeter, since the signal is often enhanced by the reference arm.Recently, clearing techniques [144,145] with low tissue distortions have been developed in order to increase tissue transparency, and to image through the entire tissue.OCT can also benefit from these approaches [146], even if we would expect the intrinsic signal to be reduced in OCT, since clearing tends to remove lipid structures and to average the refractive index fluctuations.However, it can help to detect exogenous OCT contrast agent at important depths, such as metallic nanoparticles.Here again, clearing processes are often long, and can induce tissue distortions.
Fluorescence and OCT techniques are quite complementary, and multimodal systems combining the two can help with diagnosis, as we will illustrate throughout this section.Many such multimodal setups have been developed with both scanning OCT approaches [5,147] and full-field approaches, using either a flip mirror to switch from one modality to the other [18,148], or a simultaneous approach [41,42].If all these multimodal approaches share a common interest, scanning OCT approaches are more easily combined with scanning fluorescence approaches, such as confocal or two photon microscopies, whereas FF-OCT systems have been developed in combination with widefield fluorescence techniques, such as structured illumination microscopy (SIM).SIM acquisition speed intrinsically depends on the fluorescence camera frame rate divided by the number of grid modulation steps.The SIM acquisition rate can go up to 100 frames per second in our measurements, even though the usual frame rate limit comes from the fluorophores' brightness, and from the minimal exposure time required to produce a decent image.Interestingly, simultaneous combination of FF-OCT and SIM [41,42] enable the parallel acquisition of both modalities, so images with both modalities can be acquired at more than 100 frames per second.However, SIM is not very robust to aberrations and multiple scattering, and can hardly penetrate inside scattering samples after 30 to 50 µm.Nevertheless, in mostly transparent samples, such as cornea or retina, with only a few structures labeled, we could perform multimodal FF-OCT and SIM images through the entire sample.
The first interest of combining OCT with fluorescence is to identify structures that are captured in OCT, in order to better understand and characterize OCT contrast [62].The second interest is to use fluorescence to reveal some structures, proteins, or ultimately any structure that cannot be captured by OCT, therefore combining fluorescence molecular contrast with structural contrast from OCT.It can help to understand how a structure of interest integrates and evolves inside its micro environment, such as neurons in the intestine [149], or to find the location of stem cells in the cornea [12].Figure 7 shows the example of imaging a tumorous intestine tissue, in which cell nuclei are labeled with fluorescence.This image shows the complementary contrast that can be obtained with both modalities.Additionally, with both sources of information, we can localize the cells on the contours of organized structures, named intestine villi, indicating that this field of view has been acquired in a healthy region of the cancerous tissue.A final interest, that is often neglected, especially in the context of histology, is that fluorescence enables dynamic measurements that can probe the tissue biochemistry and its evolution, using specific probes.For examples, fluorescence measurements can follow the electrical activity of cells by quantitatively measuring their intracellular calcium concentration [150,151], or directly measure the membrane potential [152].It can be useful to assess the viability of a cortical or retinal explant from its spontaneous activity, or can assess the endocrine pancreatic function, as insulin production is triggered by calcium uptake.Fluorescence can also directly measure cell function by measuring the intracellular ATP concentration [153], and can specifically detect apoptotic cells, with the ClickIt TUNEL assay (ThermoFisher).This can be of major importance to study the efficiency of pharmaceuticals on a cancerous tissue, as the cellular recovery of apoptosis within the tumor is usually a promising sign of a treatment's success.
Figure 8 shows the result of the ClickIt TUNEL assay (from ThermoFisher Scientific [154]) on a fresh liver tissue.The tissue is first purposely excised with a corner shape to elicit and follow the cell death propagation from the margin to the center of the tissue.After several minutes, the tissue is imaged with FF-OCT and D-FF-OCT, and an inactive region (dark region from the tissue margin) can be detected with D-FF-OCT (panels A to C). Directly after the imaging, the tissue is fixed, in order to start the TUNEL assay and detect the apoptotic regions of this tissue.Unfortunately, this assay requires fixed tissues, so we cannot overlay a D-FF-OCT, and a fluorescence image.However, we could detect that almost all cells in a 200 µm region from the tissue border were in apoptosis (panels D to F).The fluorescence signal started decreasing as we went further from the tissue margin, and is almost equal to zero at the center of the tissue.Figure 8 tends to show that apoptotic cells do not produce a dynamic contrast, at least for cells in apoptosis for a long time.Surprisingly, this indirect result is different from what was obtained with OCT measurements [155], where the authors showed a decrease of the speckle decorrelation time (i.e., an increase of the intracellular activity) in apoptotic cells.Explanations for this result may be the following: In contrast to G.Farhat et al. [155], we used a brute force method to generate apoptosis, as we physically cut the sample, and imaged close to the tissue edge.It is likely that the cells in the dark D-FF-OCT region are in late stages of apoptosis.Moreover, since cells in apoptosis tend to dilute and homogenize their intracellular material [156], it should lower the backscattered signal from the cell to a level that may be beneath our system detection limit, thus creating a dark D-FF-OCT image, even if the dynamics have increased.

Comparison with Other Popular Novel Microscopies
As previously pointed out, although fixation and labeling protocols have been greatly improved these past years in histology, this field could still benefit from the important technological developments that optical imaging has known over the past twenty years.So far, we compared FF-OCT and its multimodal versions to standard histology procedures.However, to be fair, we should also compare them to other recently developed optical techniques, such as multiphoton microscopies, and coherent anti-Stokes Raman scattering (CARS) microscopy, that also are label-free and enable optical sectioning.Similarly to OCT, these techniques rely on the intrinsic properties of tissues to generate contrast.They both can provide optical sectioning, since they rely on the combination of two or three photons, that can statistically only happen at the focal plane where the power density is high enough.Second harmonic generation (SHG) is measured in dense non-centrosymmetric structures, such as microtubule bundles, and collagen I fibers [157].Third harmonic generation (THG) is measured at interfaces displaying high index mismatch, such as cell or nuclear contours, and lipid droplets [158].Both techniques, especially when combined, are therefore useful to study 3D tissue organization, e.g., in the cornea [159], or in developing embryos [160].The association of the two techniques gives a contrast similar to FF-OCT, but non-linear microscopes are a little more sensitive and specific than the FF-OCT signal, and can be naturally combined with two-photon fluorescence microscopy [161].The transverse resolution is often higher than FF-OCT, as high NA objectives are required to generate a high power density at the focal point.We can add that SHG has already been extensively used to assess tissue disorganization, especially in the context of liver fibrosis [162].A commercial SHG microscope for histology is already available at HistoIndex [163].
CARS [164], and its hyperspectral version [165], are other promising non-linear techniques for histology, which reveal a high density of either a single or multiple molecular vibrations inside a sample.Interestingly, CARS is one of the only techniques that is both label-free, and specific at the molecular level, as vibration spectra of given structures can be measured, and computationally extracted [166].It is therefore one of the most promising label-free techniques today for quantification of tissue structure, and molecular content.Unfortunately, to our knowledge, there is no commercial CARS system available, and CARS systems often require careful and time consuming alignment steps.
Finally, we can add that all these techniques are based on non-linear processes that require similar kinds of sources, i.e., femtosecond lasers, and can be ingeniously combined all together from a single laser source [167].However, such non-linear processes are quite weak.They typically display a low signal-to-noise ratio, especially in dense tissues, and often require long exposure times to generate a measurable signal.Combined with high NA objectives, the measurable field of view within one hour of imaging is fairly limited.Additionally, femtosecond lasers are expensive, especially compared with FF-OCT sources that are simple thermal light or LEDs.Finally, as non-linear processes, SHG, THG, and CARS signals are mainly forward scattered, so they often have to be used within thin tissue sections (a few hundred microns to 1 mm).We can nevertheless note that multiple scattering happening in highly scattering tissues allows the observation of non-linear processes in an epi configuration [164].
FF-OCT appears better suited than competing techniques for fast imaging of large tissue volumes in the operating room, to complement frozen-section histology procedures.

Conclusions and Perspectives
Full-field optical coherence tomography is a modified version of OCT that can access submicrometric resolution deep inside a sample.Similarly to OCT, it can provide label-free 3D imaging of ex vivo tissues, and can perform measurements similar to histology in terms of contrast and resolution.The interest of FF-OCT lies in its ability to obtain 3D images similar to histology with almost no tissue preparation, which can reduce slicing, staining, and fixation artifacts.FF-OCT can also significantly shorten the biopsy-to-diagnosis time, down to a few minutes.Unfortunately, FF-OCT diagnosis scores are not yet 100% when compared with histology.In this review, we highlighted three promising recent multimodal systems combining FF-OCT structural contrast with other specific contrasts such as mechanical, motility, and molecular contrasts, providing hope that FF-OCT diagnosis scores will increase in the future, and might even ultimately perform better than histology.We need to add that despite our choice here, there are other promising multimodal systems that could be of interest for diagnosis improvement.For example, polarization FF-OCT [21] can be used to increase the sensitivity to fibers, and should be able to detect fiber disorganizations, similarly to OCT measurements, e.g., that could detect polarization changes caused by age-related macular degeneration (AMD) or glaucoma [168].Multicolor FF-OCT [22] might also help the detection of subtle wavelength-dependent scattering changes, related to tumor malignancy [69].Despite all these exciting developments, and the possibility to combine many different modalities in a single optical system, we have to emphasize that a compromise will have to be found between speed, processing time, compactness, penetration depth and the number of contrasts available, especially when FF-OCT is used during surgery.We believe that a better characterization of FF-OCT (and D-FF-OCT) contrasts that would be performed specifically for each tissue would help to increase FF-OCT diagnosis ability, without even changing the optical design of the microscope.
Beyond the scope of histology, we hope that the development of such multimodal systems will increase the potential interest in FF-OCT, not only as an ex vivo diagnosis tool, but also as a tool to answer complex biological questions.We showed in this article that FF-OCT reveals a label-free original contrast, and can be used either in ex vivo tissues, or in vivo in small animals.If scanning OCT techniques were now to be extensively used in vivo and in clinics, FF-OCT would suffer from the slow speed of 2D detectors that do not allow a dynamic correction of the sample position, which could cancel the sample motion; in vivo implementation is therefore challenging.However, a few in vivo experiments have been successfully performed in paralyzed or anesthetized small animals [35,61], and in human skins pressed against an imaging window.We recently successfully performed a simultaneous FF-OCT, D-FF-OCT, and fluorescence measurement in living zebrafish larvae.We illustrate our ability to obtain a motility contrast in vivo, even in the presence of blood flow in Figure 9.The dynamic signal is saturated in the blood vessels, and below, but cells can be detected in between the vessels.Additionally, an important effort is being made to enable FF-OCT measurements in patients, with the recent development of two adaptive optics FF-OCT systems [169,170] for retinal imaging, a FF-OCT system coupled with an OCT system to dynamically compensate for eye motion for cornea and retina assessment, and a dark field FF-OCT system [171] to record internal fingerprints.Additionally, our group recently developed an en face matrix based approach of OCT, which allows to increase the penetration depth by a factor of two [172], and could be advantageous for in vivo experiments.

Cells
Blood Vessels Finally, we are now interested in the development of functional imaging in FF-OCT, in fresh tissues, or in living animals, either by measuring changes in the blood flow, similarly to what has been achieved in OCT [132,173], or by directly measuring changes in the tissue biomechanics, or cellular activity.We believe that diagnosis capabilities will further benefit from this research, since it might soon be possible not only to detect static structural changes in a tissue, but also to detect how it reacts to different stimuli, or different drugs.

Figure 1 .
Figure 1.FF-OCT systems layouts.(A) Standard layout, as in the LLTech commercial system.It is based on a Linnik interferometer, with identical microscope objectives in both arms.The system can be illuminated by a simple thermal source, or a standard LED, and the image is recorded with a CMOS camera; (B) Principle of 3D acquisition in FF-OCT.A 2D plane can be recorded in a single camera acquisition, then the field is scanned in the axial direction; (C) Another FF-OCT layout developed in the lab, offering the possibility to work with higher NA objectives, and to add a fluorescence channel simultaneously to the FF-OCT image.

Figure 2 .
Figure 2. Comparison between FF-OCT and Histology.(A,B) Large field image of a breast sample affected by ductal carcinoma in situ (DCIS) acquired with FF-OCT (A) and histology (H & E staining; (B-E) Zoom in from boxes represented in panel A, respectively showing a normal fibrous tissue, enlarged ducts, and the collagen organization around the carcinomatous cells in a tumorous region [14].These images have been acquired by LLTech, and are displayed here with their authorization.Scale bars in A and B represent 1 mm, and 200 µm in C to E.

Figure 3 .
Figure 3. FF-OCT and biomechanical contrast in various tissues obtained with static elastography.Panels A and B show a standard four-phases FF-OCT of an ex vivo human breast tissue and the corresponding strain map (equivalent von Mises strain) after sample compression.As expected, a greater deformation can be observed around the adipose inclusions (black domains in panel A) when at the surrounding fibrous breast tissue.The color bar corresponds to a strain in percent from 0 for blue to 10 for red.Panel C corresponds to a FF-OCT cross-sectional image of an ex vivo rat heart.In Panel D, we show the corresponding angle of the principal strain direction.The color bar corresponds to the angle in degree from 0 for blue to 90 for red.Panels E and G respectively correspond to 3D volumetric FF-OCT images of a mechanically isotropic and anisotropic polyvinyl alcohol (PVA) polymers, and panels F and H correspond to their respective fractional anisotropy maps.The color bar corresponds to the fractional anisotropy from 0 for blue to 0.7 for red.Panels I and J finally show a standard FF-OCT image of a biopsy of bovine kidney, and its associated stress map obtained with the "tactile" approach.The color coded stress map varies from −0.5 kPa in blue to 0.4 kPa in red.All the orange scale bars correspond to 100 µm.

Figure 4 .
Figure 4. Dynamic elastography performed with FF-OCT by measuring shear wave propagation.This figure displays different FF-OCT images of an ex vivo rat brain recorded at 30,000 Hz with a field of view of 471 × 235 µm 2 , and their evolution during shear wave propagation.Panel A shows a standard four-phase FF-OCT amplitude image of the sample.Panels B to D show the evolution of direct images after subtraction of the incoherent intensity, during shear wave propagation.Panels B to D have been respectively acquired at t = 0 ms, t = 4 ms, and t = 5 ms after the shear wave emission.Finally, Panel E shows a plot of the temporal evolution of the average image intensity along the y axis.The horizontal scale bars represent 50 µm and the vertical scale bar 1 ms.

Figure 5 .
Figure 5. Dynamic FF-OCT accessible contrasts.(A,B) Same field of view of an ex vivo subcutaneous pancreatic tumor respectively with FF-OCT and D-FF-OCT.Panel A shows a fibrous environment dense in collagen fibers without specific cues of cancerous cells.On the other hand, panel B identifies slowly varying cell bodies (in blue), and rapidly varying immune cells (bright white cells); (C,D) Zoom corresponding to the white areas of panels A and B respectively; (E) D-FF-OCT image of a fresh rat liver with 10X 0.3 NA objectives, based on a frequency analysis to compute the color.D-FF-OCT can reveal hepatocyte cytoplasms, and nuclei (black region at the center of cells), organized along stationary capillaries, characterized by the presence of rapidly varying red blood cells trapped inside;(F) D-FF-OCT image of a similar fresh rat liver with 40X 0.8 NA objectives.The increased resolution allows to better capture intracellular dynamics, and to capture some dynamic signal inside the nucleus; (G) High resolution image of red blood cells trapped in a capillary of a fresh rat retina, acquired with the 40X 0.8 objectives.The better resolution allows to capture the enhanced dynamics localized at the red blood cells membrane.

Figure 6 .
Figure 6.D-FF-OCT as a histology tool for the brain.This figure presents a combination of FF-OCT (Panel A) and D-FF-OCT (Panel B) in a thick coronal section of fresh ex vivo rat brain.This image has been acquired at a few dozen microns below the section surface and presents the first cellular layers of the cortex.It reveals both myelinated axons in FF-OCT, as well as the neuronal cell bodies in D-FF-OCT.Panel C shows a frequency analysis performed in D-FF-OCT, which reveals highly dynamic nuclei in red, neuronal cytoplasms with intermediate dynamics (mostly in green) and mostly static myelin fibers in blue.The scale bar represents 60 µm.

Figure 7 .
Figure 7. Revealing cells in tumors with FF-OCT/SIM.A fixed tumorous tissue of intestine is imaged in FF-OCT (A); in fluorescence (B); and the overlay is displayed in Panel C. Fluorescence reveals cell nuclei in the tissue.These images show that the cells are located within a well defined structure captured by FF-OCT, indicating that we are imaging in a healthy part of the tissue.The scale bars represent 40 µm.

Figure 8 .
Figure 8. Imaging apoptosis with multimodal FF-OCT systems.This figure presents the imaging of the edge of an ex vivo liver.Here, the liver has been excised with a triangular edge to illustrate the propagation of cellular death from the excised edge.First, a large field FF-OCT and D-FF-OCT images have been acquired with the commercial LLTech system (A, B and overlay in C) 30 min after the excision.Then, the tissue has been fixed, and labeled with ClickIt Plus TUNEL Alexa 488 assay (ThermoFisher) to reveal apoptotic cells; Panels D, E and F respectively show the FF-OCT and fluorescence signals (and the overlay) acquired from the part of the tissue in the yellow box in panel A. In this region close to the excised edge, almost all cells are apoptotic, as suggested by the low level of dynamic signal in panel B. Scale bars for panels A, B, and C represent 125 µm, and scale bars in D to F represent 40 µm.

Figure 9 .
Figure 9.In vivo dynamic FF-OCT in a 2 dpf zebrafish larva.Despite strong phase variations caused by the blood flow, which generates a strong dynamic signal, some cells can still be revealed in between the capillaries.The scale bar represents 40 µm.