Second harmonic generation light quantifies the ratio of type III to total (I + III) collagen in a bundle of collagen fiber

The ratio of type III to type I collagen is important for properly maintaining functions of organs and cells. We propose a method to quantify the ratio of type III to total (type I + III) collagen (λIII) in a given collagen fiber bundle using second harmonic generation (SHG) light. First, the relationship between SHG light intensity and the λIII of collagen gels was examined, and the slope (k1) and SHG light intensity at 0% type III collagen (k2) were determined. Second, the SHG light intensity of a 100% type I collagen fiber bundle and its diameter (D) were measured, and the slope (k3) of the relationship was determined. The λIII in a collagen fiber bundle was estimated from these constants (k1-3) and SHG light intensity. We applied this method to collagen fiber bundles isolated from the media and adventitia of porcine thoracic aortas, and obtained λIII = 84.7% ± 13.8% and λIII = 17.5% ± 15.2%, respectively. These values concurred with those obtained with a typical quantification method using sodium dodecyl sulfate–polyacrylamide gel electrophoresis. The findings demonstrated that the method proposed is useful to quantify the ratio of type III to total collagen in a collagen fiber bundle.

www.nature.com/scientificreports/ Against this background, this study has presented a novel method to quantify the ratio of type III to total (type I + type III) collagen in a bundle of collagen fiber using second harmonic generation (SHG) light. The method was applied to bundles of collagen fiber isolated from the media and adventitia of porcine thoracic aortas. We verified the method's accuracy by comparing the ratio of type III collagen/total collagen estimated using our proposed method with that obtained using SDS-PAGE.

Strategy
For this study, we assumed that the collagen fiber was made up of type I and type III collagen only, a common assumption adopted in similar studies [4][5][6] . Hereafter, λ III refers to the type III collagen/total collagen ratio, calculated from the equation below: where V I and V III are the volumes of type I and III collagen, respectively, in a sample.
A schematic of the strategy to estimate the ratio of type III to total collagen, reflected by λ III , in a bundle of collagen fiber is shown in Fig. 1. SHG light intensity is a function of λ III and the volume of a bundle of collagen fiber. We use these two functions to estimate λ III in each collagen fiber bundle.
The first relationship is that SHG light intensity (I SHG-λIII ) decreases linearly with an increasing λ III (Fig. 1a), which was described by Lutz et al. 16 . Mathematically, this relationship is expressed as follows: where k 1 and k 2 are constants to be determined from the experiments. Constant k 1 means the slope of the fitted line and k 2 is the SHG light intensity at λ III = 0 (i.e., when the sample is made of pure type I collagen). Thus, the ratio of SHG light intensities (δ) of samples including type III collagen to pure type I collagen samples is expressed by: where I SHG-TypeI-gel is SHG intensity of samples constituted from pure type I collagen.
(1) III = V III V I + V III (2) I SHG− III = k 1 III + k 2 (3) δ = I SHG− III I SHG−TypeI−gel = I SHG− III k 2 , Figure 1. Two equations about SHG intensities are required to estimate λ III . Two relationships of type III ratio-SHG intensity (λ III -I SHG-λIII ) (a) and fiber bundle diameter-SHG intensity of type I collagen (D-I SHG-TypeI ) (b) were obtained from separately performed experiments to determine k 1 , k 2 , and k 3 . For a tested collagen fiber bundle where λ III is unknown, D and SHG light intensity (I SHG ) of the fiber bundle were measured (picture in the right bottom in (b). I SHG-TypeI was calculated from the D in the graph (b). Since the ratio of I SHG :I SHG-TypeI equals to the I SHG-λIII :k 2 , λ III can be determined. This study used the diameter of fiber bundle because the diameter should be proportional to the amount of collagen fiber bundle. This figure was created using Adobe Illustrator (Adobe, CS6).  16 is that the SHG light intensity in pure type I collagen (I SHG-TypeI ) increases linearly with an increasing total collagen volume (Fig. 1b). In our study, collagen fiber diameter (D) was used as a substitute to reflect total collagen volume (Fig. 1b). The relationship between SHG light intensity and collagen diameter was calculated as follows: where k 3 is a constant expressing the slope of the Eq. (4).
When attempting to estimate λ III of a given sample of collagen fiber bundle with a diameter D, we first compare the measured SHG light intensity of the actual collagen fiber bundle I SHG with that predicted if the bundle was made up of pure type I collagen (I SHG-TypeI ), calculated from Eq. (4). The I SHG value should be lower than the I SHG-TypeI value because the bundle comprises type III collagen. From this difference, we can estimate the ratio of SHG light intensity δ (Eq. 3) caused by collagen type III fiber, expressed as follows: Hence, to estimate the λ III of a given collagen fiber bundle, we can combine Eqs. (2)-(5) to obtain the following equation: Equation (6) tells that the λ III of the collagen fiber bundle can be estimated from the diameter D and SHG light intensity I SHG if the constants k 1 , k 2 , and k 3 are readily determined.
In this study, to determine the relationship between fiber bundle diameter and SHG light intensity for type I collagen (D-I SHG-TypeI ), tail tendon samples from rats were used owing to the high concentration of type I collagen (95%) in these samples 17 . This relationship (D-I SHG-Tendon ) can be described as follows: where k 3 ′ is a constant determined from experiments and denotes the slope of the Eq. (7). Type III/total collagen ratio, reflected by λ III , of the tendon is estimated as λ III-Tendon . The ratio of SHG light intensities in the tendon (δ Tendon ) including type III collagen to pure type I collagen can be determined by combining Eqs. (2) and (3) as follows: where I SHG-TypeI is calculated as follows (adapted from Eq. 4): Equation (9) shows that the constant k 3 is determined from the constants k 1 , k 2 , k 3 ʹ, and λ III in the tendon fiber λ III-Tendon , which are determined experimentally.

Results
Relationship between SHG light intensity and type III/total collagen ratio. The relationship between I SHG-λIII and λ III was determined using collagen gels made by mixing type I and III collagen. Figure 2a-l shows typical SHG light images of gels mixed with type I and III collagen. The SHG intensity I SHG-λIII of collagen gel appears to decrease with λ III increase in both the backward and forward directions of photomultiplier tubes. In the backward direction, SHG intensity was lower than in the forward direction. Our quantitative analysis confirmed these results (Fig. 2m). The SHG intensity I SHG of the collagen gel was significantly correlated with λ III in both the backward and forward directions (R 2 = 0.85, P < 0.001). From this estimation, we were able to determine constants k 1 = − 1.39 and k 2 = 171.8 in the forward and k 1 = − 0.89 and k 2 = 110.3 in the backward directions as described in Eq. (2). In further analyses, the SHG signal was collected in the backward direction to minimize the effects of fiber orientation and organization on the image 18 . Figure 2n shows the forward/backward ratio of SHG light intensity, which was almost constant in all samples examined, with no significant correlation observed (R 2 = 0.34, P = 0.98). This is in contrast to a previous study 19 , which reported that the forward-backward ratio of SHG light intensity decreased with an increase in type III collagen concentration.
Relationship between SHG light intensity and diameter of fiber bundles for type I collagen. Constant k 3 ′ in the relationship between D and I SHG-TypeI of a type I collagen fiber bundle was determined in rat tail tendons. Figure 3a,b show typical images of collagen fiber bundles obtained from the rat tail tendons in this study. The bundle diameter appears to be constant for each sample. SHG light intensity I SHG-Tendon is signifi- www.nature.com/scientificreports/

λ III of collagen fibers measured with SDS-PAGE. SDS-PAGE was conducted for two different sam-
ples. First, SDS-PAGE was applied to collagen fiber bundles of rat tail tendons to determine k 3 in Eq. (9) by obtaining λ III-Tendon . Second, SDS-PAGE was applied to collagen fiber bundles of porcine aorta to check the validity of the proposed method. Figure 4 shows a typical image obtained using SDS-PAGE. Although the molecular weight of type I collagen (283.3 ± 1.3 kDa) is similar to that of type III (288.3 ± 1.3 kDa) 20 , including urea in this method successfully distinguished these bands (lanes 2 and 3 in Fig. 4). Collagen fiber bundles from rat tail tendon (lane 4) and adventitia (lane 6) have a higher density of type I collagen (α 1 (I), α 2 (I), and β(I)) compared to type III collagen (α 1 (III) and β(III)), while collagen from media (lane 5) shows a lower density of type I collagen compared to type III collagen.
According to SDS-PAGE, the λ III in rat tail tendons (λ III-Tendon ) was 13.1% ± 3.8%, lane 4 in Fig. 4). Thus, we calculated the gradient of the relationship between SHG light intensity I SHG-TypeI and diameter D in a hypothetical 100% type I collagen fiber bundle, then determined constant k 3 = 3.79 as shown in Eq. (4) (Fig. 3c). Our quantitative analysis showed that the λ III of the aortic media obtained with SDS-PAGE was 73.9% ± 1.5% and that of the adventitia was 15.5% ± 5.1%. λ III of collagen fiber bundles of thoracic aorta using SHG light intensity. The method proposed was applied to collagen fiber bundles isolated from the media and adventitia of porcine thoracic aortas. Figure 5 shows typical images of collagen fiber bundles obtained from the media and the adventitia of porcine thoracic aortas with 5 × optical zoom (see Supplementary Information S1 online, image without optical zoom). When comparing bundles of similar diameters, samples obtained from the media (Fig. 5a) showed a lower SHG intensity than those obtained from the adventitia (Fig. 5b). This tendency was also observed in thin bundles (Fig. 5c,d); in both the media and the adventitia, thick bundles ( Fig. 5a,b) had higher intensities than thin bundles (Fig. 5c,d).
The λ III of the medial and the adventitial collagen fiber bundles were calculated using Eq. (6). Figure 6 shows the λ III value obtained from the SHG method along with that of SDS-PAGE. The λ III was significantly higher in the medial collagen (84.7% ± 13.8%, n = 16) than in the adventitial (17.5% ± 15.2%, n = 17). These results were similar to those obtained using SDS-PAGE.

Discussion
In this study, we proposed a new method to quantify the ratio of type III to total collagen (λ III ) of a given fiber bundle from its diameter and SHG light intensity. We also confirmed the validity of our method against a known procedure (SDS-PAGE). Our method is based on two relationships: the first between type I collagen diameter and SHG light intensity; the second between λ III and SHG light intensity 16 . Combining these enables us to quantify λ III within a collagen fiber bundle. As the measurement of λ III using previous methods is timeconsuming (taking at least several hours), our method is beneficial. We applied this method to collagen fiber bundles obtained from the media and the adventitia of porcine thoracic aortas and found that λ III in the media and adventitia concurred with the data obtained using SDS-PAGE (Fig. 6). Moreover, these results also concur www.nature.com/scientificreports/ with other previous studies: the ratio of type III and I to total collagen in the human aortic media are 70% and 30%, respectively 21 and adventitia is almost exclusively composed of type I collagen 22 . Thus, we can say that our proposed method is suitable for the evaluation of λ III . An increase in the ratio of type III resulted in a decrease in the SHG intensity (Fig. 2). Tilbury 19 reported that relative increase in the type III ratio within fibrils induced changes in the structure of collagen fibrils: decrease in fibrillary organization shown by shorter and more randomly arranged collagen fibers: decrease in the dipole moment shown by SHG signal anisotropy within fibers: an increase in the pitch angle of α-helices shown by pixel-based polarization-resolved approach. Such changes in the structure of fibrils with increasing type III collagen would have decreased SHG intensity with an increase in the type III ratio.
Our SHG light images contained heterogeneous intensity distributions within the collagen fiber bundles observed. The intensity at the center of the fiber bundle image was higher than in the peripheral area (Fig. 3a,b). This appears to be caused by the higher laser power in the central region of the images by the optical devices used in this study. To avoid this heterogeneity, we used the same setting throughout the experiment. However, heterogeneous distribution can be seen even within a collagen fiber bundle on a scale of only a few micrometers (e.g., in Fig. 5). It is not clear whether the heterogeneous distribution reflects the actual spatial difference in the local density of type III collagen, as there is no method to verify this hypothesis. Since the level of the standard deviation in SHG intensity shown in Fig. 5 is much larger than that of Fig. 2, this high frequency is not caused by noise. Previous studies demonstrated that glycation of collagen decreased SHG intensity 23 and a lower level  www.nature.com/scientificreports/ of the crosslinking increased SHG intensity 24 . It is thus considered that the heterogeneity of SHG light intensity is possibly produced by the local non-enzymatic glycation and lysyl oxidase-mediated crosslinking of collagen.
Our proposed method has several advantages: first, it allows researchers to quantify λ III more quickly than existing methods such as SDS-PAGE, gene expression 25 , and immunofluorescence 26 , and picrosirius red staining 27 . Our method only requires a collagen fiber bundle to be placed and observed under a microscope, which is much quicker and handier than previous methods. Potential applications of the proposed method include cosmetic, healthcare, and clinical research. Second, this method measures λ III without the need for invasive procedures such as fixation and homogenization. Therefore, we can simultaneously measure multiple properties (e.g., mechanical properties) from a fiber bundle.
Although the method proposed in this study has not been fully established, it has wide potential for practical application. For instance, this method can be applied to tissues in vivo as it uses SHG light in the backward direction (i.e., the light does not have to penetrate through the whole thickness of the specimen). Since the SHG light of human skin in vivo can be observed 28,29 , λ III of the skin can be estimated. Although our method requires the isolation of a collagen fiber bundle from a given tissue, techniques to quantify the volume of collagen in tissues enable us to measure the λ III of skin, albeit invasively. Moreover, the spatial distribution of λ III in a specimen can be obtained. Although we have to investigate whether the distribution of intensity is accurate, measuring the local distribution of type I/III ratio in tissues can be used as an indicator for skin wound healing 12,13 .
Users have to bear in mind that this method uses the intensity. When the SHG light of a collagen fiber bundle is repeatedly captured, the intensity is slightly decreased after every image captured (see Supplementary Information S2 online). Although this decrease was small (0.4% after each image capture) and did not affect the results of this study, this decay might affect results in future measurements and should thus be taken into account. Furthermore, the light intensity varies depending on the devices used, such as lasers, quantum efficiency and sensitivity of photon multiplier tubes, as well as optical devices (e.g., objectives and filters). To use our proposed method, users will need to adjust the constant parameters k 1 , k 2 , and k 3 in their system. Preparing samples of the various concentrations of type III collagen and several diameters of collagen fiber bundles can help establish a reference standard.
In this study, we assumed that the collagen fiber was made up of type I and type III collagen only, as commonly assumed in past studies [4][5][6] . Reportedly, collagen compositions of human arteries were type I of 60%, III of 30% and others (type V and minor collagens) of 10% 30 , indicating that type I and III were dominant compared to other types. These compositions are obtained from human tissues, not from the porcine ones used in this study. Carrasco et al 31 implemented comparative studies on the collagen of vertebrate arteries, and reported that vertebrates shared morphologic and histochemical features of collagens in arteries. These reports suggest that it could be reasonable to assume that the porcine aorta is made mainly of type I and type III collagens like the human aorta, although they are not identical.
Care for other types of collagen and non-collagenous tissues is needed in the application of the proposed method to actual biological tissues. The potential effects of non-collagenous proteins, for example, collagenbound small leucine-rich proteoglycans that regulate collagen fibril formation and possibly intermolecular space, on SHG cannot be entirely excluded. SHG signals thus may include signal generated from non-collagenous proteins. Type II collagen is the major component in articular cartilages 32 . Because type II collagen generates strong SHG light 33 , direct application of the proposed method to such tissues might not be adequate. Type IV and type V collagen are relatively common in tissues, but they do not generate the SHG light 33 . Elastin and cells also do not affect SHG light, because SHG light generates from noncentrosymmetric molecular assemblies 34 . In support of it, high SHG light intensity areas do not colocalize with non-collagenous tissue regions 18,32 . These indicate that the method can be applied to most biological tissues except type II collagen-rich ones such as articular cartilage.
In SDS-PAGE, λ III was quantified with the β-component and α chains to obtain a more accurate measurement. To consider the effect of the β-component on our result, we also measured the λ III using only the bands of α chains. The values obtained from the rat tail tendon and the aortic media and adventitia were 10.9% ± 3.0%, 75.3% ± 2.5%, and 11.6% ± 3.5%, respectively. Compared with the results shown in Fig. 6, the λ III value showed a < 3% difference in the three kinds of specimens, suggesting that the effect of β-component is not large enough to change our conclusion. More accurate results would be obtained if γ-component was included although it would have minor effects on SHG.
Reportedly, the ratio of type III collagen in the tendon was less than 5% 35 , which is smaller than our result (13.1% ± 3.8%). One possible cause for this might be the way to select the area of type III band in SDS-PAGE gel. Because the band of the rat tail-tendon (lane 4 in Fig. 4) was not clearly recognized, we defined the type III region in reference to the purified type III band (lane 3 in Fig. 4). In this process, we might have overestimated the type III band region. Other possible factor is the background noise of SDS-PAGE staining. As described in the method section, the intensity of the whole image was subtracted by the background value. It however remains unclear if the background noise of SDS-PAGE staining is totally canceled.
For SDS-PAGE analysis, samples were treated with pepsin, the solubilized collagen was precipitated and subjected to SDS-PAGE analysis. The data of the λ III obtained with the SDS-PAGE was compared with the λ III obtained by the proposed method that uses SHG light intensity, and they were quantitatively in good agreement. These λ III were also congruent with Mccullagh and Balian 21 that showed the ratio of type III and I to total collagen in the human aortic media are 70% and 30%, respectively. Only soluble collagens were however analyzed in Mccullagh and Balian 21 . According to Menashi et al 36 , type I/III collagen ratio is 2:1 if the whole collagen including insoluble ones is analyzed. Since λ III obtained with the proposed method is close to the type I/III ratio of soluble collagens, the proposed method may quantify the type I/III ratio of soluble collagens.
At the beginning of this study, we attempted to quantify λ III from the difference in SHG light intensity recorded on the photon multiplier tubes between the backward and forward directions, since the forward-backward ratio of SHG light intensity was reported to decrease with increasing type III collagen concentration 19  www.nature.com/scientificreports/ forward-backward ratio was almost constant regardless of λ III gels (Fig. 2n). Since the strength of SHG intensity is affected by the width of fiber in the forward direction, thinner fibers in the gels might have led to the constant value of the forward-backward ratio. Therefore, when using SHG light intensity to measure λ III , the backward direction is more suitable.
In conclusion, we propose a novel method to quantify the ratio of type III to total collagen in various tissues, based on SHG light intensity and the diameter of collagen fiber bundles. We applied this method to collagen fiber bundles in the media and adventitia of porcine thoracic aortas and confirmed its measurement accuracy and suitability by comparing data obtained with SDS-PAGE.

Methods
Quantification of the relationship between SHG light intensity and λ III . Type III collagen solutions of various λ III were prepared by mixing 3 mg/mL of type I (PSC-1-100-20, Nippi, Tokyo, Japan) and type III (PSC-3-100-20, Nippi) collagen solutions in ratios of 5:0, 4:1, 3:2, 2:3, 1:4, and 0:5. A silicone rubber container (internal dimensions: 15 × 24 × 3 mm) was placed on a glass slide and cleaned with ethanol. Then, 1-5 µL of 1 M NaOH filtered with a φ0.20-µm filter was added to each collagen solution to obtain a final concentration of 0.1%, and 1 mL of each mixed solution was poured into the container and incubated at 37 °C for 24 h to obtain a gel. Thereafter, the container was removed and a coverslip was placed on the collagen gel to prepare the observation sample. A control sample was prepared using 1 mL of incubated phosphate-buffered saline [PBS(−)].
SHG light intensity of the collagen was imaged under a multiphoton microscope (FV1200MPE, Olympus, Tokyo, Japan), as reported previously 18,37 . Briefly, a laser (wavelength, 800 nm; pulse width, 100 fs; repeated frequency, 80 MHz) was applied to the prepared gels through a 25 × objective lens (NA = 1.05, XLPLN25XW, Olympus). The generated SHG light was observed using a bandpass filter (400 ± 5 nm). Data were recorded on photomultiplier tubes in both forward (transmitted) and backward (reflected) directions. Observation conditions were the same for all specimens (HV 600 V, laser power = 18%, dwell time = 10 µs/pixel, Karman filter = 2). To reduce the influence of structural variations of collagen fiber such as direction and dimeter, we imaged 51 slices with an interval of 2.0 µm at 5 positions per specimen. Five specimens were tested for each experimental condition, and the average was used.
Image analyses were performed using ImageJ (v. 1.51i, National Institute of Health, Bethesda, MD, USA). From the image stack, 40 slices, which cover the whole thickness, were selected and the average intensity in the image stack I gel was measured. I SHG-λIII was determined by subtracting the background intensity in the SHG image I BG , as follows: Finally, fitting constants k 1 and k 2 in Eq. (2) were determined using the least-squares regression from the relationship between I SHG-λIII and λ III .

Quantification of the relationship between SHG light intensity and the diameter of type I collagen fiber bundles. All animal experiments were approved by the Institutional Review Board for
Animal Care of the Nagoya Institute of Technology, following recommendations from their Guide for Animal Experimentation and this study is reported under ARRIVE guidelines (https:// arriv eguid elines. org). Three male Wistar/ST rats (15-week-old, Chubu Kagaku Shizai, Nagoya, Japan) were used in this experiment. After the rats were euthanized using CO 2 , their tail tendons were excised and collagen fiber bundles (D = 5-100 µm) were extracted using forceps. Both ends of the fiber bundle were glued on a polyethylene terephthalate (PET) film sheet with a cyanoacrylate adhesive. The samples were then sent to a laboratory-made tensile tester 38 to stretch the specimen. When the tester recorded 2 mN of force, the specimen was deemed to have been sufficiently stretched. We selected a small sample size to reduce the animal experiments but multiple animals to consider the individual difference (three rats).
SHG light intensity from a collagen fiber bundle was imaged under a multiphoton microscope (FV1200MPE, Olympus) using a 5 × objective lens (MPLN5X, NA = 0.1, Olympus). In this condition, image depth was 109.3 µm. Thus, a single slice was considered to include all depth-wise SHG light generated from a collagen fiber bundle with a diameter several 10 µm. However, because the long axis of the collagen fiber bundle might not be parallel to the image plane, several slices were imaged with 10-µm intervals. Observation conditions were the same for all specimens (HV, 250 V; laser power = 18%; dwell time = 10 µs/pixel; Karman filter = 2). After measuring the SHG light intensity, the fiber bundle was imaged again using a 5 × optical zoom function in addition to the 5 × objective lens to measure its exact diameter.
The maximum intensity projection image was obtained from the stacked image of the fiber bundle. The fiber bundle region was defined and the average SHG light intensity in the region I F was recorded from triplicate measurements. Similarly, average SHG light intensity in the background area I BG was measured in triplicates, and the SHG light intensity generated from the tendon collagen fibers I SHG-Tendon was determined as follows: From the fiber image obtained using the 5 × optical zoom function, the diameter was manually measured at any three locations and its average was determined as fiber diameter D. From the relationship between I SHG-Tendon and D, constant k 3 ′ was determined in Eq. (7) using a least-squares regression analysis.
Measurement of λ III in a collagen fiber bundle obtained from porcine thoracic aorta using SHG light intensity. Porcine thoracic aortas were purchased from a local slaughterhouse (Handa meat market,  www.nature.com/scientificreports/ Handa, Japan) and transported to our laboratory in PBS( −) at 4 °C. The connective tissues were removed with forceps, and the ventral side of the aorta between the first and the second intercostal arteries from the proximal side was dissected using a surgical knife. The adventitia was separated from the aorta with forceps, while the medial specimens were cut into blocks (5 mm in the longitudinal and 10 mm in the circumferential directions), immersed in a 6% w/w agar solution, and kept at 4 °C for 30 min. The specimens were then sliced into 100-µm sections using a micro slicer (DTK-1000, Dosaka-EM, Kyoto, Japan). Adventitial and medial samples were separated into in microtubes and immersed in 68. 5 20 . Therefore, we added urea in the gel to distinguish the band positions of type I and III collagens, as described previously 39 .
Collagen fiber bundles from rat tail tendons and porcine aortic media and adventitia samples were obtained, as described above, and specimen solutions were obtained as described previously 40 . Specimens were homogenized in 20% w/v NaCl in 0.05 M Tris buffer (pH 7.5) using an ultrasonic homogenizer (VP-050 N, TAITEC, Koshigaya, Japan). The samples were centrifuged and the sedimentation was resuspended in 0.6 mg/mL pepsin (P7012, Sigma-Aldrich, St. Louis, MO, USA) and 0.5 M acidic solution for 24 h. The solution was centrifuged and its supernatant was dialyzed overnight in 0.02 M disodium hydrogen phosphate (pH 9) using a 20-kDa cassette (87,735, Thermo Fisher Scientific, Waltham, MA, USA). The dialyzed solution was centrifuged, and the sediment was obtained as a sample for SDS-PAGE.
The sample buffer used for SDS-PAGE comprised 62.5 mM Tris buffer (pH 6.8), 2% SDS, 25% glycerol, and 0.5% bromophenol blue. The prepared collagen samples were dissolved in an equal volume of sample buffer and 0.5 M Tris buffer (pH 6.8) to obtain a final concentration of 50 mg/mL. Additionally, 3 mg/mL of type I (PSC-1-100-20, Nippi, Tokyo, Japan) and type III collagen (PSC-3-100-20, Nippi, Tokyo, Japan) solutions were separately prepared for use as control samples for type I and III collagen. After adding 2-mercaptoethanol to obtain a 5% concentration, the solution was spun down, incubated at 100 °C for 5 min, and centrifuged.
The sample was electrophoresed at 200 V for 100 min and proteins were visualized with Coomassie Brilliant Blue R-250. The gel was imaged using ChemiDoc MP Imaging System (Universal Hood III, Bio-Rad). SDS-PAGE was performed in triplicates using different samples.
For image analysis, the intensity of the whole image was subtracted by the background value. Then, band regions were selected in reference to type I (lane 2 in Fig. 4) and type III samples (lane 3 in Fig. 4) and the intensity of each region was summed. Although dimers (β-component) and trimers (γ-component) have been reported previously 41 , regions of α 1 (I)-and α 2 (I)-chains and β(I)-components were considered type I collagen and those of α 1 (III)-chain and β(III)-component were considered type III collagen. The measurements were performed in triplicates for each gel.
Statistics. Pearson's correlation coefficient was calculated for the plots between the I SHG-λIII -λ III , forward/ backward ratio-λ III , I SHG-Tendon -D, and I SHG-TypeI -D. P values of < 0.05 were considered statistically significant. Data were averaged and presented as mean ± standard deviation (SD). All analyses were performed using Microsoft Excel (2016).

Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.