Polarization-resolved Second Harmonic microscopy in anisotropic thick tissues

: We thoroughly analyze the linear propagation effects that affect polarization-resolved Second Harmonic Generation imaging of thick anisotropic tissues such as collagenous tissues. We develop a theoretical model that fully accounts for birefringence and diattenuation along the excitation propagation, and polarization scrambling upon scattering of the harmonic signal. We obtain an excellent agreement with polarization-resolved SHG images at increasing depth within a rat-tail tendon for both polarizations of the forward SHG signal. Most notably, we observe interference fringes due to birefringence in the SHG depth proﬁle when excited at π / 4 angle from the tendon axis. We also measure artifactual decrease of ρ = χ xxx / χ xyy with depth due to diattenuation of the excitation. We therefore derive a method that proves reliable to determine both ρ and the tendon birefringence and diattenuation.


Introduction
Second Harmonic Generation (SHG) microscopy is an efficient imaging technique to visualize the three-dimensional (3D) distribution of fibrillar collagen in biological tissues [1,2]. Incident circular polarization is usually used since it enables imaging of fibrils independently of their orientation in the focal plane. This approach is particularly relevant for quantitative biomedical studies such as fibrosis scoring [3,4]. However, fibrillar collagen exhibits a structural anisotropy that may be interesting to characterize for many applications. For that purpose, polarization-resolved SHG provides complementary information about the 3D distribution of nonlinear dipoles within the focal volume [5,6,7,8]. The usual approach is to measure the ratio ρ of the two main tensorial components of the nonlinear response, considering a cylindrical symmetry for fibrillar collagen [9, 10,11,12,13,14]. This ratio depends on the orientational distribution of the collagen triple helices and on the orientation of the nonlinear dipoles (along the peptidic bonds) within the triple helix. ρ therefore provides information about the ordering of collagen molecules or fibrils within the focal volume. It has been measured in various collagenous tissues and showed a wide dispersion of values (1.2 to 2.6). Comparison of healthy and pathological tissue however requires to characterize the accuracy of ρ measure- ments and their sensitivity to various artifacts. In particular, the anisotropy of many collagenous tissues may affect polarization-resolved SHG experiments as reported in tendons [10,15,16]. Polarization-resolved image processing improved when taking into account diattenuation [15] and birefringence [10,16]. Furthermore, polarization scrambling due to scattering has been shown in tendons and other biological tissues [17]. The aim of this paper is therefore to thoroughly characterize the linear propagation effects that affect polarization-resolved SHG experiments in anisotropic tissues and to develop a reliable method to determine ρ in collagenous tissues. For this purpose, we first develop a theoretical model that fully accounts for birefringence, diattenuation and polarization scrambling upon scattering. We then record polarization-resolved SHG images in rat-tail tendons that is a model tissue composed of aligned collagen fibrils. We finally compare simulated and experimental data that show an excellent agreement and give insight into the linear and nonlinear optical properties of the tendon. Most notably, we observe interference fringes in the SHG depth profiles when excited at π/4 angle from the tendon axis, and artifactual decrease of ρ with depth due to diattenuation of the excitation. We conclude that the correction for these effects enables a reliable determination of ρ.

Tendon preparation
Tendons were extracted from Sprague-Dawley rat-tails (female, ≈ 300g), centrifugated at 4700 rpm and stored at 4 • in phosphate buffer saline (PBS). Imaging was performed within a few days using a water-immersion objective. Tendons were first labelled with fluorescent latex beads to enable precise location of the tissue surface (1µm diameter, L1030, Sigma-Aldrich). After rinsing, they were fixed at both extremities and streched to get rid of the crimps and facilitate alignment of the fibrillar pattern along the x direction of the microscope stage. They were kept in PBS on a glass coverslip that maintained them in the focal plane of the microscope.

Imaging setup
Multiphoton imaging was performed using a custom-built laser scanning microscope [18] as depicted in fig. 1.a. Excitation was provided by a femtosecond Titanium-sapphire laser tuned at 860 nm and scanned in the xy directions using galvanometric mirrors. It was focused in the tendon using a water-immersion 20x, 0.95 NA objective that resulted in typically 0.4 µm lateral and 1.6 µm axial resolution near the sample surface. Power was adjusted to 10-15 mW at the focus using a rotating half waveplate and a Glan prism that filtered out y-polarization before entering the microscope setup. Nonlinear optical signals were detected using photon-counting photomultiplier tubes and appropriate spectral filters to reject the laser excitation (FF01-680SP, Semrock) and select two-photon excited fluorescence (2PEF) signal (GG455 high-pass filter, Schott) or SHG signal (HQ430/20 interferential filter, Chroma). 2PEF was detected in the backward direction and SHG in the forward direction. Multimodal images were recorded using 100 to 200 kHz pixel rate, 0.4 to 0.8 µm pixel size and 1 to 2 µm z-step. They were combined using Matlab and ImageJ softwares as shown in fig. 1.c. Note that there is a slight index mismatch between the tendon (n = 1.5) [20] and water (n = 1.33), so that the depth z within the tendon is related to the lens displacement dz lens by: z = dz lens n tendon /n water .

Polarization-resolved measurements
The laser polarization at the back pupil of the objective was approximately along the x direction but exhibited a 14% ellipticity due to the optical components within the microscope (9% ellipticity without the dichroic mirror for epi-detection). We therefore inserted an infrared polarizer and achieved an ellipcity less than 1% with small scanning angles. This linear polarization was tuned from −2π/3 to 2π/3 (usually with π/12 steps) by rotating an achromatic half waveplate (MWPAA2-22-700-1000, CVI-Melles Griot) placed just before the objective (see fig.  1.a). Forward SHG signals were analyzed using a polarizing beamsplitter cube (BBPC-550, CVI-Melles Griot). To improve the extinction ratio of the x-and y-polarized detection channels, linear polarizers (03FPG021, CVI-Melles Griot) were set in front of the detectors. The relative transmission of these two channels was calibrated using a fluorescent slide (Chroma): we took advantage of the isotropy of the setup within the xy plane and compared the x-polarized signal to the y-polarized one excited with polarizations shifted by π/2 angle (see fig.1.b). The ratio of both channels was typically 1.1. Calibration was performed before every experiment and enabled quantitative comparison between x-and y-polarized signals.

Polarization-resolved SHG
Let's consider the nonlinear optical interaction of the incident laser beam with a rat-tail tendon. The polarization induced in the medium by the elecric field E is given by: where we considered only the first-and second-order electric susceptibility tensors χ (1),(2) . Rattail tendon is commonly assumed to have a cylindrical symmetry (C 6v symmetry) [9, 10,7], which reduces the number of independent nonvanishing χ (2) components. Moreover, we assume that the Kleinman symmetry is valid as usually considered because of the nonresonant character of the interaction. Within these approximations, there are only 2 independent nonvanishing χ (2) tensorial components: χ xxx and χ xyy = χ xzz = χ yxy = χ zxz = χ yyx = χ zzx , where x represents the main axis of the tendon [19]. Considering a laser beam propagating in the z direction with a linear polarization at angle α to the tendon lying in the xy plane (see fig. 1.a), the electric fields reads: E ω x = E 0 cos α and E ω y = E 0 sin α near the focus and induces a SH polarization: which radiates at frequency 2ω. In this paper, we are interested in ratiometric measurement of the second-order response and we consider the ratio of the two independent tensorial components of the second-order susceptibility: ρ = χ xxx /χ xyy . This ratio provides insight into the orientational distribution of the collagen molecules within the focal volume. It is taken to be real in the Kleiman approximation. The SH intensity detected for each polarisation then reads: where K is a constant merging various parameters such as setup geometry and squared incident beam intensity I 2 0 . Polar diagrams of eq. (3) are displayed in fig. 2.a and b using ρ = 1.4. A common method to determine ρ from polarization-resolved SH experiments is to fit I 2ω x using eq. (3a).
In the following, we propose a more general method. Since I 2ω x and I 2ω y are even functions of α and contain only even powers of trigonometric functions, they can be represented as a sum of cos 2nα functions: In that framework, one can notice that: This expression enables the determination of ρ from polarization-resolved SH data by fitting I 2ω x using eq. (4a). The advantage of this method is that it can take into account various physical effects that appear to modify the expression of I 2ω x but not its representation as a sum of cos 2nα functions. These physical effects will be presented in the following sections and the new expressions of the A, B and C parameters will be derived.  2)). (c) I 2ω x and (d) I 2ω y at 24µm depth considering birefringence (∆n = 0.0066), diattenuation (∆l a = 175µm) and polarisation cross-talk due to scattering (η xy = 0.13) (see Eq. (11)). ρ = 1.40 in all simulations.

Tendon birefringence
In the former derivation of I 2ω , we have assumed that the electric field propagation takes place in an isotropic medium. Nevertheless, the tendon is considered to exhibit a C 6v symmetry, so that the permittivity tensor ε 0 1 + χ (1) i j has two different components and the tendon is analogous to a uniaxial birefringent crystal. Indeed, it has been reported that birefringence ∆n = n e − n o attains 5 · 10 −3 [20,21] in tendon. As a consequence, the x (extraordinary wave) and y (ordinary wave) excitation field components undergoes a relative phase retardation when propagating within the tendon. SH intensity then reads: accounts for the birefringence in the excitation propagation (SHG intensity is unaffected by birefringence in the harmonic propagation). Note that z represents the depth in tendon at which SH takes place: due to index mismatch, z = dz lens n tendon /n water , where dz lens is the microscope lens displacement from the tendon upper surface.
These expressions show that I 2ω y is unaffected by birefringence, while I 2ω x has a supplemen- tary term compared to eq. (4). Decomposing I 2ω x in the cos 2nα basis, we get: One observes that B biref = B and A biref +C biref = A +C, so that ρ can be determined using eq. 6 as in section 3.1. To get better insight into the effect of birefringence, we have plotted in fig.  3 the variation of I 2ω x and I 2ω x as a function of the incident polarization angle α for increasing depth within the tendon. Comparison of model calculations with and without birefringence shows that birefringence results in interference fringes in the x-polarized SH intensity I 2ω x when excited with a linear polarization at π/4 from the tendon axis, as shown by the cos ∆φ term in eq. 7a (see z-profiles in 3.k and 3.n).

Polarization cross-talk
We have assumed yet that propagation occurs in an optically perfect medium. However, tendon is a scattering medium like most biological tissues. It results in the decrease of the excitation intensity since scattered light is not intense enough to induce a nonlinear response. This effect, that may be different for x-and y-polarized incident light, will be considered in the next section. Here, we are interested in the effect of light scattering on the harmonic beam polarization. When propagating in an anisotropic scattering medium, a wave with well-defined polarization accumulates scrambling over its polarization direction since scattering processes slightly rotate the polarization direction. As a consequence, a small amount of SH light that is initially polarized along y is detected in the x-polarized channel and vice versa. Let η XY (resp. η Y X ) be the amount of that polarization "cross-talk" from y detection channel to x detection channel (resp. from x to y), so I 2ω reads: Decomposing I 2ω x in the cos 2nα basis, we obtain: Once again, B pol = B and A pol +C pol = A+C, and ρ can be determined as before. Fig. 3.b and g show that the angular profiles of I 2ω x and I 2ω y are slightly modified by polarization cross-talk.

Diattenuation in tendon
Multiphoton microscopy uses near-infrared excitation which minimizes absorbance and scattering in biological tissues. Nevertheless, light scattering and spherical aberrations cause focal volume deterioration upon propagation within the tendon, which results in smaller excitation intensity and weaker SH signal. For simplicity, we assume that the excitation intensity undergoes exponential attenuation with depth. Due to the tendon anisotropy, the effective attenuation lengths for I ω x and I ω y are different and the excitation light experiences diattenuation. These attenuation lengths will be noted as l e a and l o a to account for the uniaxial symmetry of the tendon. We assume that diattenuation of the SH signal is negligible because signal forward-detection in multiphoton setups is not much sensitive to light scattering and aberrations. I 2ω x,y then reads: The expression of the A, B and C parameters are the same as in eq. 5 with the modified parameters: l o a and ρ → ρe − z ∆la . Consequently, once ∆l a is known, ρ can be determined in a similar way as in section 3.1. Note that this approach can be generalized to an arbitrary diattenuation profile f a (z) = f x a (z)/ f y a (z) by using the transformation ρ → ρ f a (z) in eq. 6. Fig. 3.c, h and m show the modifications of I 2ω x and I 2ω y induced by diattenuation. In particular, z-profiles for an incident polarization parallel (resp. perpendicular) to the tendon exhibit a l e a (resp. l o a ) exponential decay.

Combined effects of birefringence, polarization cross-talk and diattenuation
When taking into account all the effects introduced above, we obtain the following expression: l o a ρe − z ∆la cos 2 αe i∆φ + sin 2 α 2 + η XY e − z ∆la |sin2α| 2 (12) Fig. 2.c and d display the polar diagrams of this expression along with the similar expression for I 2ω y . They clearly look very different from the polar diagrams a and b plotted from the simplified model (eq. 3). It shows that determination of ρ using eq. 3a is hazardous. Conversely, I 2ω x can still be decomposed as a sum of cos 2nα functions as in eq. 4a: the parameters A eff , B eff and C eff are obtained by combining the effects of birefringence (eq. 8) and polarization crosstalk (eq. 10) only, and considering the modified parameters: l o a , ρ → ρe − z ∆la and η xy → η xy e − z ∆la . In that framework, ρ can be determined as follows: given that z-profiles of I 2ω x for α = 0 and π 2 can be used to determine l e a and l o a (or f a (z) in a more general case).
The right hand side of fig. 3 displays the x-and y-polarized SH signal given by eq. 12: diattenuation is revealed by differences in z attenuation for incident polarization parallel and perpendicular to the tendon axis, birefringence by oscillations in the z-profile at π/4 polarization angle from the tendon axis and polarization cross-talk by deformation of the angular profiles. In order to better evidence the latter two effects, we introduce the parameter ∆ that sorts out the contributions of ∆φ and η XY to A eff and C eff : At a given depth, contributions of birefringence and polarization cross-talk are intrinsically mixed, but information on these effects can be obtained by analysis of ∆ depth variation.

Fig. 4 displays typical
x-and y-polarized SHG images obtained at various depths within a tendon for an incident beam polarized at π/4 from the tendon axis. Transverse profiles show that the SHG signal is higher at the tendon edges as expected because of the cylindrical shape of the tendon: attenuation of the excitation beam is more effective in the center of the section when the incident light propagates along a larger distance within the tendon. However, the xand y-polarized SHG transverse profiles exhibit different behaviour as a function of the depth from the tendon upper surface. In fig. 4.a and b, the x-polarized SHG image is more attenuated than the y-polarized SHG image at the tendon center, whereas it is not the case in fig.  4.c and d. Moreover, the x-polarized SHG image at 57µm ( fig. 4.c) displays more signal at the center than at intermediate positions along the transverse profile. This complex behaviour is characteristic for birefringence effects. The x and y components of the π/4-polarized incident beam experience different optical indices when propagating within the tendon and accumulate a relative phase shift. The tensorial components of the nonlinear response excited by the x and y components of the π/4-polarized incident beam then interfere in a constructive or destructive way, depending on the propagation distance within the tendon. It results in interference fringes in the depth profile of the SHG signal for incident excitation mixing up x and y-polarization components. These fringes appear only in x-polarized SHG images because they result from the coherent summation of two tensorial components excited along x and y polarizations respectively (see eq. 2). Conversely, y-polarized SHG images are not sensitive to birefringence because they probe only one tensorial component. To get a better insight into the physics of polarization-resolved SHG images, we carried out SHG experiments at increasing depth within the tendon. Fig. 5.a and b display typical x-and y-polarized SHG signals in 28 x 20µm regions of interest such as the one depicted in fig. 1.a. The depth profile of the 2PEF signal from the latex beads is also displayed in fig. 5.c. It enables the location of the tendon surface z 0 that corresponds to the maximum of the 2PEF peak. The specific patterns observed in the SHG images were consistently observed in all our samples provided that the ellipticity of the incident polarization was negligible. Practically, experiments performed without a linear polarizer at the back pupil of the objective lens exhibited a non symmetrical pattern relative to the tendon axis (α = 0). Fig. 5.f displays z-profiles of the x-polarized SHG signal for various incident polarization angles. Interference fringes characteristic for birefringence in the excitation propagation are clearly observed for an incident polarization at π/4 angle to the tendon axis. The distance between the two maxima is δ z ≈ 47µm. In a first approximation, it is related to the birefringence phase shift by ∆n = λ /2δ z. It gives: ∆n ≈ 0.008 (after index mixmatch correction), in qualitative agreement with reported birefringence in tendon [20,21]. It will be determined more precisely later on. Diattenuation of the laser excitation is suggested by the z-profiles of the xpolarized SHG signal at 0 and π/2 angles. Exponential fitting gives the following attenuation lengths for x and y incident polarizations : l e a = 91µm and l o a = 190µm (see fig. 6.a). Finally, fig. 5.d and e display I 2ω x and I 2ω y at the surface of the tendon (z = z 0 ) where birefringence and diattenuation effects do not apply. We therefore fit our data using eq. 9 that accounts only for polarization cross-talk due to light scattering (or equivalently, using eq. 12 with z = 0 and ∆φ = 0). We observe that the experimental data are better fitted using eq. 9 than in the simplified approach (eq. 3). It shows that polarization cross-talk due to light scattering significantly distorts the polarization pattern at the tendon surface.     x (α) and (e) I 2ω y (α) at the tendon surface (black dots), along with fits using the simplified approach (eq. 3, red line) or accounting for polarization cross-talk (eq. 9, blue line). (f) experimental z-profiles of I 2ω x for α = 0 (green dotted line), α = π/4 (blue dashed line) and α = π/2 (red dash-dot line).
Altogether, these experimental observations show that the usual approach with eq. 3 does not satisfactorily account for polarization-resolved SHG imaging in thick tendons. Birefringence, diattenuation and polarization cross-talk due to scattering must all be taken into account to explain experimental data. Accordingly, we used eq. 4 with parameters A eff , B eff and C eff to fit our data. Note that we introduced an angular shift α 0 to account for possible slight misalignment of the tendon axis from x direction (−2 • < α 0 < 2 • for all our data). We then calculated ρe −z ∆la by use of eq. 13. This value was subsequently corrected for diattenuation using l e a = 91µm and l o a = 190µm obtained from exponential fitting of I 2ω x (α = 0) and I 2ω x (α = π/2). Both the raw and corrected values of ρ are depicted in fig. 6.b. Diattenuation correction successfully removes the artifactual variation of ρ with depth and enables the reliable determination of this parameter: ρ = 1.40 ± 0.03. Moreover, as already stated, this approach is applicable to any attenuation profile f a (z). To fully characterize the optical response of our sample, we also calculated the parameter ∆(z) from eq. 14a and plotted its depth profile in fig. 6.c. Eq. 14b shows that the values at the maxima (∆φ = 2πk with k integer) reduce to η xy . The first maximum is obtained at the tendon surface: η xy = 0.13. The second maximum is zero which shows that ∆la η xy is negligible ≈ 60µm deep within the tendon. We then consider that η xy diminishes exponentionally from 0.13 to 0 within the first 40µm from the surface and we get ∆n = 0.0066 by fitting our data using eq. 14b (∆n = 0.0074 without any index mismatch correction).
Finally, to verify the consistency of our data processing, let's compare our experimental data ( fig. 5.a and b) to the simulations using our experimentally determined parameters for birefringence, diattenuation and polarization cross-talk ( fig. 3.e and j). The theoretical and experimental intensity maps are in excellent agreement which proves that our theoretical approach satisfactorily reproduces polarization-resolved SHG experiments in tendons. ρ determined from polarization-resolved SHG measurements with (black) and without (green) correction for diattenuation. (c) ∆ parameter evidencing birefringence (oscillations), polarization cross-talk due to scattering (non-vanishing value at the tendon surface) and diattenuation (exponential attenuation with depth). The solid line represents fitting with the following parameters: ∆n = 1.40, η xy = 0.13 near the surface and ∆l a = 134µm.

Discussion
In this paper, we developed a method to account for linear propagation effects affecting the polarization when determining ρ by polarization-resolved SHG experiments. Previous works studied the effects of diattenuation [15] and birefringence [10,16] separately, but they did not propose an approach accounting for both effects and enabling the determination of ρ at any depth. Our model considers birefringence and diattenuation along the excitation beam propagation and polarization cross-talk due to scattering of the SHG signal. Our experimental results show that this model perfectly fits our data on rat-tail tendons and enables the reliable determination of ρ. We also observed interference fringes in the SHG depth profiles which show unambiguously that birefringence affects polarization-resolved SHG microscopy. Our model may be refined by numerical simulations of the field distribution within the focal volume using a vectorial approach to properly describe the polarizations [22]. Such a calculation would account for polarization mixing through the high numerical aperture objective lens [22,23] and for deformation of the focal volume due to birefringence. The description of the polarization may also be further refined by using Mueller matrices and tracking Stokes vectors through the tendon [24]. The advantage of our model however is to include all the characteristic features of the biological sample in a phenomenological approach, including light scattering or depolarization. It presumably explains why it shows an excellent agreement with our experimental data and proves relevant to determine optical parameters of rat-tail tendon.
Our data gives insight into both the linear and nonlinear optical response of the tendon. It first enables the determination of the tendon birefringence at the excitation wavelength. We obtain ∆n = 6.6 10 −3 in good agreement with recently reported values measured in a more  [20,21]. We also measure the diattenuation at 860nm, ∆l a = 175µm, and the polarization cross-talk due to scattering, η xy = 13%. Our measured values seem reasonable considering usual optical properties in biological tissues. Previous work also reported polarization scrambling in tendon that could be reduced by optical clearing [17]. Our phenomenological parameter η xy = 13% amounts to an average SHG polarization rotation of about 20 • . It shows that the SHG polarization is reasonably well preserved in the forward direction although propagation to the detector takes place over several scattering lengths. Indeed, forward detection picks up mainly quasi-ballistic photons and the polarization of these photons is not strongly modified. Similarly, the measured attenuation lengths (l e a = 91µm and l o a = 190µm) are in good agreement with our previous measurements in rat-artery [26] and the diattenuation shows qualitative agreement with previously reported measurements in rat-tail tendon using OCT (∆l a = 125µm [20,21]) or in horse flexor tendon using SHG microscopy (∆l a = 480µm [15]). Note that the precise determination of these parameters depends on the tendon index considered for index mismatch correction. Finally, as intended from these polarization-resolved SHG experiments, they determine the ratio ρ of the two main tensorial components of the tendon second-order susceptibility. We obtain ρ = 1.40 ± 0.03 from the experimental data displayed in fig. 5. Similar values in the range 1.3 to 1.5 were consistently obtained in all our samples. The reproducibility of these measurements is excellent considering the dispersion of tissue properties in biological samples. Slight differences from one tendon to the other may moreover be attributed to differences in the tendon preparation, particularly in the tendon stretching, that could affect the 3D distribution of fibrils at a submicrometric scale.
Most importantly, our measurements show that the raw value of ρ varies a lot with depth as depicted in fig. 6.b. It decreases from 1.40 at the surface of the tendon to 0.8 at ≈ 90µm depth within the tendon. These data proves that correction for diattenuation is essential for reliable measurements of ρ in thick anisotropic tissues and that our model provides an efficient correction along the full stack of data. We expect that the determination of ρ deeper in the tissue would be hampered because of the low signal to noise ratio. When possible, the most reliable method is to perform polarization-resolved measurements at the sample surface to get rid of diattenuation and birefringence effects.
Finally, note that this method is also applicable to epidetected SHG signals. In that configuration, the infrared polarizer set before the half waveplate at the back pupil of the objective has to be replaced by a quarter waveplate set at a suitable angle to correct for incident ellipticity without rejecting some epi-SHG signal (see fig. 1.a). Polarization analysis of this epidetected SHG signal is a complex task because of the presence of the waveplates. Furthermore, backwarddetected SHG is highly affected by scattering and related polarization distortion [25,17]. However, SHG polarization analysis was required here for properly characterizing the different linear optical effects that may affect the determination of ρ. Once this study is completed, the total SHG intensity can be fitted like the x-polarized component as a sum of cos 2nα functions, using an expression similar to eq. 12. Favourably, this configuration does not require the alignment of the tendon axis relative to polarization analysis.

Conclusion
In this paper, we combined model calculation and experiments in rat-tail tendon to characterize linear optical effects that affect the polarization in polarization-resolved SHG experiments. We evidenced that birefringence, diattenuation and polarization scrambling upon scattering significantly distort the SHG response for tunable incident linear polarizations. In particular, birefringence results in interference fringes in the depth profile of the SHG signal excited with polarization at π/4 angle from the tendon axis. Most importantly, diattenuation results in an ar- tifactual decrease of ρ with the depth within the tendon. To address these problems, we derived a method to process polarization-resolved SHG data and we successfully retrieved relevant optical parameters in rat-tail tendon. This method is applicable to any anisotropic sample that exhibits SHG signals, including other collagenous tissues and presumably skeletal muscles. It may also be generalized to other nonlinear optical processes, that is 2PEF, THG, CARS or SRS polarization-resolved microscopies. Altogether, our work proves unambiguously that uncorrecting for polarization distortion results in misleading determination of the tensorial nonlinear response in anisotropic thick tissues.