Translator Disclaimer
29 November 2019 Probing ECM remodeling in idiopathic pulmonary fibrosis via second harmonic generation microscopy analysis of macro/supramolecular collagen structure
Author Affiliations +

Idiopathic pulmonary fibrosis (IPF) is a progressive disease with poor prognosis with short lifespan following diagnosis as patients have limited effective treatment options. A fundamental limitation is a lack of knowledge of the underlying collagen alterations in the disease, as this could lead to better diagnostics, prognostics, and measures of treatment efficacy. While the fibroses is the primary presentation of the disease, the collagen architecture has not been well studied beyond standard histology. Here, we used several metrics based on second harmonic generation (SHG) microscopy and optical scattering measurements to characterize the subresolution collagen assembly in human IPF and normal lung tissues. Using SHG directional analysis, we found that while collagen synthesis is increased in IPF, the resulting average fibril architecture is more disordered than in normal tissue. Wavelength-dependent optical scattering measurements lead to the same conclusion, and both optical approaches are consistent with ultrastructural analysis. SHG circular dichroism revealed significant differences in the net chirality between the fibrotic and normal collagen, where the former has a more randomized helical structure. Collectively, the measurements reveal significant changes in the collagen macro/supramolecular structure in the abnormal fibrotic collagen, and we suggest these alterations can serve as biomarkers for IPF diagnosis and progression.



Idiopathic pulmonary fibrosis (IPF) is a form of interstitial lung disease that is a progressive disorder characterized by unrelenting scarring and stiffening of the lungs that leads to 34,000 deaths each year. The interstitial fibrosis is characterized by spatial and temporal heterogeneity; for example, areas of dense collagen accumulation (old scar) are juxtaposed with fibroblastic foci (new scar formation).1 As scarring increases, efficient lung function decreases, compromising the uptake of oxygen by blood, resulting in shortness of breath, limited exercise capability, and daily cough. The median survival for patients diagnosed with IPF is typically 3 to 5 years post-diagnosis.

Currently, IPF is clinically diagnosed using a combination of tests including chest x-ray, high-resolution computed tomography (HRCT), lung biopsy, and lung function tests. These techniques often lack the sensitivity and specificity needed to examine the extracellular matrix (ECM) changes and their role in IPF progression. For example, when the classic radiographic pattern is present, IPF can be diagnosed via HRCT.2 However, this pattern is not always uniquely defined due to the heterogeneity of the disease. Moreover, no radiographic method can probe the collagen morphology changes that are an integral part of the pathology. Biopsy followed by histology remains the gold-standard for IPF diagnosis; however, there is a significant risk of associated morbidity and longitudinal samples cannot be taken.3

There remains a clear need for better diagnostics as well as prognostic indicators. We suggest that probing the underlying macro/supramolecular changes in collagen presents such a new direction. Perhaps counterintuitively, the specific collagen changes in this fibrosis have not received significant attention beyond conventional H&E staining, showing increased deposition. As an alternative, second harmonic generation (SHG) microscopy offers considerable opportunities to study collagen alterations and provide insight into both disease etiology and progression. We previously showed that SHG combined with machine learning based on two-dimensional wavelet transforms/principal component analysis of the fiber morphology classified normal and IPF tissues with near 100% accuracy.4 In addition, we found that the collagen/elastin balance increased, consistent with increased clinical stiffness and difficulty breathing.

SHG microscopy can also be used to probe subresolution aspects of collagen architecture. For example, SHG polarization-resolved methods [e.g., linear SHG polarization analysis (P-SHG) and circular dichroism (SHG-CD)] can extract information on the helical structure of collagen (e.g., α-pitch angle and triple helical chirality).5 For example, we previously showed that collagen molecular aspects were different in human ovarian cancer from normal ovarian tissues.6 This is also relevant for IPF as there have been earlier reports using immunostaining that showed the Col I and Col III isoform balance changes during disease progression and is also different from that of normal tissue.79 Moreover, these isoforms have different pitch angles, and we have shown that these can be delineated by linear polarization analysis and application of the single-axis molecular model.10,11 Additionally, analysis of the SHG directional response, i.e., the forward-backward ratio (F/B), yields data on the subresolution size and packing. In this study, we used these SHG metrics along with bulk optical property measurements to characterize the macro/supramolecular structure of IPF and normal human lung tissues. Collectively, these optical readouts provide new insight into the collagen alterations in IPF that have been previously unattainable by other methods.


Materials and Methods


Human Tissues

All tissue samples (de-identified) were obtained from lung transplant recipients at UW Hospital Madison, Wisconsin, under a current IRB-approved protocol. The normal tissues were from pathologist-defined normal adjacent tissue from biopsies of patients without fibrotic lung disease. All IPF samples were from patients with advanced or highly progressed IPF. Tissues were fixed in 4% formalin and sectioned using a vibratome (Leica VT1200) to 200  μm thickness. We have shown that fixation does not significantly alter the fibril structure.12 After sectioning, the tissues were stored at 4°C in phosphate-buffered saline (PBS) for conventional SHG imaging or optically cleared by immersion in 50% glycerol overnight to reduce scattering-induced de-polarization effects for SHG polarization-resolved imaging. For imaging, samples were mounted on glass slides in PBS or glycerol with #1.5 coverslips and nail polish to seal the slides. A total of three normal and four IPF-independent patient samples were imaged and also used for optical property measurements.


Collagen Concentration Assay and α-SMA Staining

Using a Sirius Red Collagen Detection Kit (catalog no. 9062, Chondrex, Redmond, Washington), collagen concentration of collagen type I standards (8, 16, 31.5, 63, 125, 250, and 500  μg/ml solutions), blanks (acetic acid only), and our test samples (lung tissues) were extracted, in accordance with the manufacturer’s instructions. Each of the lung tissue samples was homogenized in 1  mg/mL pepsin in 0.05-M acetic acid and incubated for 10 days at 4°C. After collagen digestion and Sirius Red staining, the supernatant was collected, and the total collagen concentration was detected. A Tecan Infinite M1000 Plate Reader was used to measure the optical density (OD) at 530 nm to obtain absorbance readings from the standards, blanks, and lung samples. We subtracted the blank OD values from the standards and test samples. Then, we plotted the OD values of the standard curve using linear regression analysis and then calculated the collagen concentration (μg/ml) of the lung tissues. Three tissues per group, each run in duplicate, were analyzed.

α-smooth muscle actin (SMA) expression was imaged to identify fibrotic regions to be imaged by SHG. Samples were fixed in 4% formaldehyde, blocked in 1% BSA in PBS, permeabilized with 0.1% Triton X-100 for 30 min, and then incubated overnight at 4°C with the primary antibody, anti-α-SMA (monoclonal, mouse). Secondary staining was performed using goat anti-mouse IgG2a Alexa Fluor 594 (ThermoFisher, 1:1000 dilution) and then incubated overnight, then mounted in PBS and imaged by two-photon excited fluorescence (TPEF). The two-photon excitation laser wavelength was 780 nm, and the emission collection was centered at 632 nm with a bandpass filter.


SHG Microscope System

The details of the SHG microscope have been described elsewhere13,14 and are only briefly described here. The system consists of a laser scanning unit (FluoView 300; Olympus, Melville, New York) mounted on an upright microscope (BX61; Olympus, Tokyo, Japan), where the excitation source is a mode-locked titanium sapphire laser (Mira; Coherent, Santa Clara, California). Imaging was performed with a fundamental laser wavelength of 890 nm for SHG forward-backward (F/B) and P-SHG analysis and 780 nm for SHG-CD; the shorter wavelength for the latter provides greater sensitivity.5 Average powers at the focus were 30 to 50 mW using a 40× 0.8 NA water immersion lens (LUMPlanFL; Olympus, Tokyo, Japan) and a 0.9 NA condenser. The resulting lateral and axial resolutions were 0.7 and 2.5  μm, respectively. Forward- and backward-directed SHG emission was collected using matched photon-counting detectors (7421 GaAsP; Hamamatsu, Hamamatsu City, Japan), where the collection efficiencies were calibrated as before using fluorescent beads.13 The SHG wavelengths (445 and 390 nm) were isolated with the respective 10-nm-wide bandpass filters (Semrock, Rochester, New York). The excitation wavelength was confirmed using a fiber-optic spectrometer (Ocean Optics, Dunedin, Florida). Fields of view were 170×170  μm for SHG F/B and 85×85  μm for both SHG-CD and P-SHG and were acquired with scanning speeds of 2.71  s/frame with three-frame Kalman averaging. The power was controlled by an electro-optic modulator (ConOptics, Danbury, Connecticut) run by a custom LabVIEW program (National Instruments, Austin, Texas), interfaced with the FluoView scanning system using a data acquisition card (PCI-6024E; National Instruments).

Linear polarization was obtained using a half-wave plate to define the state entering the microscope, and the desired linear rotation at the focal plane was achieved using a liquid-crystal rotator (LCR; Meadowlark Optics, Frederick, Colorado) mounted in the infinity space.14 Circular polarization is achieved with a quarter-wave plate after the LCR, where left- and right-handed states are achieved with 90 deg of linear rotation by the LCR.14 The linear and circular polarization states were validated as previously described by imaging cylindrically symmetric giant vesicles.5,14 The polarization control was also run by a custom LabVIEW program interfaced to the FluoView scanning system.


SHG Polarization Analyses


Helical pitch angle analysis

Polarization-dependent measurements were performed as previously described,10 where images were taken every 10 deg through 180 deg of rotation. Here, the method was applied to both optically cleared normal and IPF tissues, where clearing is required as scattering rapidly scrambles the excitation polarization.15 We had also showed that using thin tissues (10  μm) this process does not affect the polarization response. The α-helical pitch angle is extracted10 by combining the pixel-based generic model16 with the single-axis molecular model.11 As previously, we determined the pitch angle, θp, through analysis of the symmetry reduced tensor elements:

Eq. (1)



Second harmonic generation-circular dichroism

SHG-CD analysis was used to interrogate the overall chirality of the human lung tissues; the method has been described previously.5 In brief, images of IPF and normal lung tissues were taken using left- and right-handed (LH and RH) circularly polarized (CP) laser excitation and were obtained 30-μm deep into the optically cleared tissues to avoid boundary effects. To account for variations in intensity in the different samples, we report the normalized SHG-CD response defined as

Eq. (2)

where I2(ω)LHCP and I2(ω)RHCP represent the integrated pixel intensities of the SHG images of LHCP and RHCP, respectively. This is calculated on a pixel basis, where we first set a threshold mask above the noise background to identify nonzero pixel values. Absolute values were summed across the entire field of view as the sign of CD response will depend on fiber orientation.5,17,18


Bulk Property Measurements

The spectral dependence of the single scattering anisotropy, g, and the scattering coefficient, μs, was determined using the setup previously reported by Hall et al.19 The scattering coefficient, μs, is determined by on-axis attenuation measurements through the Beer–Lambert law:

Eq. (3)

where I is the transmission with the sample, I0 is the transmission without the sample, α is the factor for losses due to refractive index mismatches, d is the tissue thickness, μs is the scattering coefficient, and μa is the absorption coefficient.19 Since lung is a collagen-rich tissue, μaμs and can be considered negligible.20

The scattering anisotropy, g, is first determined by goniometry, where the angular distribution of increasing scattering light is measured from 0 deg to 180 deg through the use of a rotating photodiode detector on a motorized stage and fit to the Henyey–Greenstein phase function (valid for most tissues):

Eq. (4)

Multiple scattering broadens the angular response and yields a lower measured anisotropy, gef. We therefore developed a Monte Carlo simulation framework to extract the true anisotropy gsingle from the measured gef.19 gsingle is associated with the scattering directionality and structural organization of the tissue on a scale from 0 to 1, with the following limits: g=0, isotropic scattering, corresponding to randomly organized structure; and g=1, all forward-directed scattering, and associated with a highly organized structure relative to the excitation wavelength. The independent determination of g and μs also yields the reduced scattering coefficient, μs (which is used as an input parameter for the Monte Carlo simulations for SHG directionality):

Eq. (5)



Determination of SHG Emission Directionality

The emission or creation directionality (FSHG/BSHG) is reflective of the fibril assembly, i.e., size and packing in the axial direction21 and is determined by measuring the depth-dependent measured F/B through image stacks of 100  μm of thickness. The depth-dependent response then results from a convolution of the FSHG/BSHG with scattering (μs) of the SHG signal at λSHG. As previously described, a Monte Carlo simulation framework22 (modified for SHG from MCML)23 using the measured bulk optical properties (Sec. 2.5) extracts the emission directionality FSHG/BSHG.20,22 The SHG emission directionality was determined both across the whole field of view as well as locally to investigate the extent of heterogeneity in the tissues.24 Here, a 15×15  pixel patch area (corresponding to 5.3×5.3  μm) was used as it optimally captures the fiber and fiber bundle structures of the tissue and is adequately robust to Poisson noise. The number of stacks for this analysis was 34 and 75 for the IPF and normal tissues, respectively.


Statistical Analysis

One-way analysis of variance (ANOVA) with Student’s t-tests (polarization-resolved SHG, SHG-CD) and Tukey’s honest significant difference post-hoc tests (bulk optical property measurements, all other SHG methods) were performed in Origin 9.1 (OriginLab, Northampton, Massachusetts). p-values less than α=0.05 were considered statistically significant.




Assessment of Overall Collagen Assembly

We first present an overall comparison of the collagen content in normal and IPF lung tissues. The top row of Fig. 1 shows representative SHG images of normal and IPF lung tissues. We note that the IPF tissue has greater coverage across the field and denser collagen accumulation in comparison to normal. IPF also appears to have thinner, wavier fiber structures, whereas the fibers are more linear in normal lung tissue. We point out that collagen morphologies tend to vary significantly across sampling areas, where regions of diseased tissue may resemble that of normal/healthy lung tissue. Still, we previously were able to differentiate these tissues with high accuracy using machine learning analysis of the SHG images.4 In order to validate the apparent difference in coverage, we calculated a packing efficiency for each group, where this is quantified by creating a binary mask over a lower threshold of 15 counts (on a 12-bit image stack) and then calculating the fraction of the resulting nonvanishing pixels. We found that IPF has a significantly higher packing (normal: 0.56±0.03, IPF: 0.69±0.04; p<0.01).

Fig. 1

Combined collagen and α-SMA staining in normal [(a) and (c)] and IPF [(b) and (d)] tissues. The top row gives SHG images only and the bottom row is an overlap of SHG (grayscale) and TPEF (red) for α-SMA, identifying fibrotic regions. Scale bar=30  μm.


We next verified that we were interrogating fibrotic regions, where this is important given the heterogeneity of the collagen morphology in both types of tissues. A standard marker for fibrosis is increased α-SMA expression, as this is associated with a fibroblast-to-myofibroblast transition, where myofibroblasts secrete significantly larger amounts of collagen and other ECM proteins than undifferentiated fibroblasts.25,26 We thus overlapped SHG and TPEF images from an anti α-SMA antibody, where this is shown in the bottom row of Fig. 1. The normal and IPF lungs show little and extensive α-SMA expression, respectively.

We further measured the collagen concentrations for IPF and normal lung tissues using a Sirius Red detection kit. The Sirius Red dye binds to the repeating Gly-X-Y pattern that forms the majority sequence in the fibrillar collagen triple helix. The measured concentrations are shown in Fig. 2, and we found approximately a two-fold increase in collagen concentration in IPF over normal tissues, in agreement with the SHG imaging data showing increased coverage.

Fig. 2

Collagen concentration data of normal (blue) and IPF (red) lung tissues. Standard error bars are shown. Number of samples were six for both IPF and normal with three separate slices in each case. Note: * indicates p<0.05.



Bulk Optical Parameter Measurements

We further characterized the respective normal and IPF tissue structures through the wavelength dependence of the scattering coefficient, μs, the single scattering anisotropy, g, and the reduced scattering coefficient, μs. We will also use these values in the Monte Carlo simulations in Sec. 3.3. The scattering coefficient, μs, is associated with the density and polarizability of a tissue and is equal to the reciprocal of the mean free path of a photon before scattering. We used a set of wavelengths (390, 495, 545, 780, 990, and 1070 nm) that are near those of the SHG (445 nm) and fundamental (890 nm), and as examples, the latter are given in Table 1. We found the IPF tissues had consistently higher values, corresponding to increased densities, consistent with the analyses in Sec. 3.1. The anisotropies were not significantly different between the tissues but increased at longer wavelengths, in good accordance with Mie theory.

Table 1

Bulk optical parameters for normal and IPF tissues tabulated as mean±standard error.

Groupλ (nm)μs (cm−1)gsingleμs′ (cm−1)
Normal (N=20)445327±230.914±0.00426.5±1.6
IPF (N=18)445405±230.908±0.00741.9±4.2

For quantitation, we compare the spectral dependence of the reduced scattering coefficient, μs, which arises from the spatial distribution of the refractive index due to structural differences on size scales smaller than the diffraction limit. We follow the treatment of Backman using the Whittle–Matérn correlation function, where the spectral dependence is given by μs(λ)λ(2m4). Here, the output is the shape factor m, which corresponds to one half of the fractal dimension,27,28 where higher m values are associated with larger, more ordered structures on the approximate size scale of 50 nm to 1  μm. The resulting μs(λ) values are plotted in Fig. 3, where the best fits for the normal and IPF tissues were m=1.61 and 1.34, respectively. In this analysis, this is a large difference that corresponds to very different tissue structures, specifically indicating IPF tissues have a larger distribution of scatter sizes that contribute to the response. We note that we observed similar behavior in comparing normal and malignant ovarian tissues, where the latter had higher values of μs due to increased collagen deposition but stronger wavelength dependence, i.e., lower m due to the reduced regularity of the fibril structure.29

Fig. 3

Spectral dependence of μs over UV/Vis and NIR wavelengths for normal and IPF tissues where the fit is to the Whittle–Matérn correlation function. The IPF tissues are more highly scattering but have a stronger spectral slope (lower m), indicting a broader range of scatter sizes. There were 20 independent measurements at each wavelength using the different tissues.



Characterization of Fibril Assembly by Local SHG Emission Directionality

SHG in tissues is characterized by nonideal phase-matching, i.e., Δk=k2ω2kω0, where ω and 2ω correspond to the fundamental and SHG angular frequencies. As a consequence, to conserve momentum, this results in a distribution of forward (FSHG) and backward (BSHG) emitted components where we have dubbed the quantity FSHG/BSHG as the creation ratio or emission directionality.21 In this treatment, lower values correspond to greater phase-mismatch and more disorganized structures relative to the size scale of λSHG. The SHG emission becomes coupled with scattering during tissue propagation and Monte Carlo simulations using μs at λSHG are used to extract the creation ratio as previously described.20 As we are interested in the heterogeneity within the tissues, the analysis is performed on 15×15  pixel patches instead of the whole field of view. In previous work, we found this size range to be optimal in examining heterogeneity.24

The measured F/B versus depth and best simulation for the creation ratio, FSHG/BSHG, are shown in Fig. 4(a) and summarized in Table 2. We note that at a few depths, there are similar experimental F/B values between the normal and IPF tissues. This is likely due to nonideal F/B behavior for some of the IPF samples due to intrinsic heterogeneity of these tissues relative to the normal lung. Still, the extracted best-fit values of the SHG creation ratio, FSHG/BSHG and associated reduced chi-squared values are negligibly affected. The locally extracted FSHG/BSHG values for the patches of representative images for normal and IPF are shown in Figs. 4(b) and 4(c), respectively. Compared to normal, IPF has a lower FSHG/BSHG creation ratio, which suggests smaller and less organized collagen fibrils in the axial direction relative to λSHG.

Fig. 4

Local analysis of the SHG emission directionality. (a) Measured F/B as a function of depth into normal (blue) and IPF (red) tissues; solid and open symbols correspond to measured and simulated responses, respectively. The resulting FSHG/BSHG in 15×15  pixel patches for (b) normal and (c) IPF tissues. Number of stacks were 34 and 75 for IPF and normal, respectively. Scalebar=30  μm.


Table 2

Creation ratios for normal and IPF lung, tabulated as mean±SE.

GroupFSHG/BSHGFSHG/BSHG std. dev.a
Normal (N=75)2.67±0.141.16±0.07
IPF (N=34)2.02±0.120.74±0.07

aCalculated per individual tissue.

We propose that the lower organization of collagen in IPF is either due to fibrosis (i.e., fibrotic collagen is intrinsically less organized) or other changes in the ECM that affect the collagen morphology. This finding is in agreement with previous SEM studies that showed that the collagen fibrils were more irregular in IPF than in normal tissue.29 We also examined the heterogeneity within these tissues. Interestingly, we found a lower FSHG/BSHG standard deviation for IPF, which suggests that Δk values are more uniform within the obtained stacks. Thus, while overall the IPF tissues have more heterogeneity, containing both normal and fibrotic regions as well, the fibrotic regions themselves are more uniform than normal tissues.


Polarization-Resolved Second Harmonic Generation

There are previous accounts using immunostaining that showed the relative proportion of Col I and Col III is different in IPF relative to normal lung. Specifically, Col III is relatively increased in early stage disease, and then Col I is more prevalent in later stages, corresponding to an older scar.25,30,31 Immunostaining is not always very quantitative in general and less so for the current case as Col III antibodies have crosstalk with those for Col I as the epitope is the same.

To examine this proposed change in relative collagen isoform abundance, we employed SHG polarization analysis based on the single-axis molecular model, which is sensitive to the α-helical pitch angle.11 Based on structural biology, the pitch angle (angle of coil relative to long molecular axis) for Col III is about 2 deg higher than that of Col I. Previously, we showed that SHG could discriminate between the fibrillar morphology of varying collagen I/collagen III concentrations in mixed gels.10 We also successfully differentiated these gels based on the extracted pitch angles, where the results were consistent with known difference of Col I and III from the protein database.32

Applying this same polarization-resolved SHG technique to image human lung tissues [Eq. (1)], we obtained pitch angles (Fig. 5) of 48.25 and 48.2 for normal and IPF, respectively, suggesting there is no measurable collagen isoform change in IPF. This could occur as there is no significant increase in Col III or it is not identifiable due to the spatial heterogeneity in IPF.

Fig. 5

Linear polarization analysis of normal (blue) and IPF (red) tissues. (a) The reconstructed pixel-based response; (b) the extracted pitch angles. The data were similar to each other, inconsistent with an increase in Col III abundance in IPF.


We further analyzed the helical properties of normal and IPF tissues via the SHG-CD protocol as described previously [Eq. (2)].5 Representative SHG-CD images are shown in Fig. 6(a) for both normal and IPF, where red and blue correspond to the sign of the response and corresponds to the polarity of the fiber and the magnitude arises from the alignment of the triple helices within the focal volume.5 The response is calculated using the absolute value [Eq. (2)]; we found the average SHG-CD was significantly higher (almost two-fold) in normal versus IPF [Fig. 6(b)]. This decreased chirality in IPF suggests either improper collagen fibril formation or changes in crosslinking. In principle, this would also be consistent with an increase in Col III; however, that is not consistent with the pitch angle analysis (Fig. 5).

Fig. 6

(a), (b) Normalized SHG-CD data of cleared normal (blue) and IPF (red) lung tissues measured at 780-nm excitation wavelength. (a) The red and blue correspond to positive and negative SHG-CD values, respectively, which are determined by the fiber polarity. (b) Standard error bars are shown. Number of unique images were 134 and 121 for IPF and normal, respectively. Field size=85×85  μm. Note: **** represents p<0.00001.




While seemingly paradoxical, the collagen macro/supramolecular changes in IPF are not well studied. SHG can directly visualize the collagen assembly, is sensitive to morphological changes within the assembly, and has significant potential, especially when coupled with polarization-resolved approaches and further with optical scattering measurements.33 For example, the SHG directional analysis of the measured F/B versus depth response combined with Monte Carlo simulations based on optical properties yields data on the relative fibril size and assembly that is consistent with the (limited) available SEM measurements.29 Moreover, the spectral analysis of the reduced scattering coefficient revealed that IPF is both more dense (higher μs, packing efficiency) and less organized (lower m) than normal lung tissue, further consistent with the lower SHG creation ratio, FSHG/BSHG. Thus, these optical measurements are consistent with both increased collagen deposition and decreased order than that of normal tissues. Importantly, both the SHG and optical properties can be performed on intact tissues, where thin sections are required for any high-resolution electron microscopy work.

We did not find any differences in the respective α-helix pitch angle and thus no difference in relative Col I/III abundance. It is possible that due to the heterogeneity within the IPF tissues, our imaging regions were not optimized. However, the SHG-CD response was very different, so there were clear changes in chirality in these locations. Moreover, there was clear enhanced α-SMA expression in these tissues (Fig. 1), consistent with fibrotic regions. It is now known there is marked difference in crosslinking in IPF relative to normal tissues,34 where this could affect the net chirality. There are other possible underlying ECM changes as well, e.g., increased fibronectin deposition,35 to which SHG is not directly sensitive.

Interestingly, we found a similar trend in our work on ovarian cancer, where Col III was also reported to be increased.36 Similarly, this was not borne out by extracted pitch angles while the SHG-CD was significantly less for ovarian cancer than normal stroma.6 Moreover, the SHG directional analysis and optical properties analysis all trended in the same direction (i.e., denser and more disordered) for IPF and high grade ovarian cancer relative to the corresponding normal tissues, suggesting commonalities in the misformed collagen. Intriguingly, several of the same pathways are altered in both diseases, e.g., upregulation of proteases.3638 Moreover, the microenvironments of fibroses and cancers have many similarities including fibroblast activation, increased collagen synthesis, and stiffness.37 SHG is thus sensitive to a range of physical changes associated with increased collagen deposition accompanying different diseases.



IPF prognosis is poor due to the lack of effective treatment options, limited knowledge of the disease etiology and underlying molecular and temporal changes associated with disease progression.4 To help solve this problem, we have used SHG imaging in combination with optical property measurements to examine macro/supramolecular and fibril changes in the fibrotic collagen. These metrics indicated significant differences in collagen assembly between the normal and IPF tissues, with the latter being characterized by increased disorder, where this is consistent with the limited available structural biology data.29 As SHG can be done on whole tissues, the ability to obtain subresolution structural data without the constraints of historical methods offers great promise for this imaging modality as a diagnostic tool. For example, a laser scanning microendoscope39 could be developed to monitor disease progression as well as response to treatment.


The authors declare no competing financial interests.


PJC and NKS gratefully acknowledge support under 1R21HL126190-01A1; PJC gratefully acknowledges NSF CBET – 1402757; and DSJ gratefully acknowledges support under NSF DGE-1747503.



A. L. Katzenstein and J. L. Myers, “Idiopathic pulmonary fibrosis: clinical relevance of pathologic classification,” Am. J. Respir. Crit. Care Med., 157 1301 –1315 (1998). AJCMED 1073-449X Google Scholar


G. Raghu et al., “An official ATS/ERS/JRS/ALAT statement: idiopathic pulmonary fibrosis: evidence-based guidelines for diagnosis and management,” Am. J. Respir. Crit. Care Med., 183 788 –824 (2011). AJCMED 1073-449X Google Scholar


C. J. Lettieri et al., “Outcomes and safety of surgical lung biopsy for interstitial lung disease,” Chest, 127 1600 –1605 (2005). CHETBF 0012-3692 Google Scholar


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


K. R. Campbell and P. J. Campagnola, “Wavelength-dependent second harmonic generation circular dichroism for differentiation of Col I and Col III isoforms in stromal models of ovarian cancer based on intrinsic chirality differences,” J. Phys. Chem. B, 121 1749 –1757 (2017). JPCBFK 1520-6106 Google Scholar


K. R. Campbell et al., “Polarization-resolved second harmonic generation imaging of human ovarian cancer,” J. Biomed. Opt., 23 066501 (2018). JBOPFO 1083-3668 Google Scholar


G. Raghu et al., “Extracellular matrix in normal and fibrotic human lungs,” Am. Rev. Respir. Dis., 131 281 –289 (1985). ARDSBL 0003-0805 Google Scholar


J. M. Kirk et al., “Biochemical evidence for an increased and progressive deposition of collagen in lungs of patients with pulmonary fibrosis,” Clin. Sci. Lond., 70 39 –45 (1986). Google Scholar


J. D. Fulmer et al., “Collagen concentration and rates of synthesis in idiopathic pulmonary fibrosis,” Am. Rev. Respir. Dis., 122 289 –301 (1980). ARDSBL 0003-0805 Google Scholar


K. Tilbury et al., “Differentiation of Col I and Col III isoforms in stromal models of ovarian cancer by analysis of second harmonic generation polarization and emission directionality,” Biophys. J., 106 354 –365 (2014). BIOJAU 0006-3495 Google Scholar


S. V. Plotnikov et al., “Characterization of the myosin-based source for second-harmonic generation from muscle sarcomeres,” Biophys. J., 90 693 –703 (2006). BIOJAU 0006-3495 Google Scholar


O. Nadiarnykh et al., “Second harmonic generation imaging microscopy studies of osteogenesis imperfect,” J. Biomed. Opt., 12 051805 (2007). JBOPFO 1083-3668 Google Scholar


X. Chen et al., “Second harmonic generation microscopy for quantitative analysis of collagen fibrillar structure,” Nat. Protoc., 7 654 –669 (2012). 1754-2189 Google Scholar


C. H. Lien et al., “Precise, motion-free polarization control in second harmonic generation microscopy using a liquid crystal modulator in the infinity space,” Biomed. Opt. Express, 4 1991 –2002 (2013). BOEICL 2156-7085 Google Scholar


O. Nadiarnykh and P. J. Campagnola, “Retention of polarization signatures in SHG microscopy of scattering tissues through optical clearing,” Opt. Express, 17 5794 –5806 (2009). OPEXFF 1094-4087 Google Scholar


J. Duboisset et al., “Generic model of the molecular orientational distribution probed by polarization-resolved second-harmonic generation,” Phys. Rev. A, 85 043829 (2012). Google Scholar


H. Lee et al., “Chiral imaging of collagen by second-harmonic generation circular dichroism,” Biomed. Opt. Express, 4 909 –916 (2013). BOEICL 2156-7085 Google Scholar


N. Mazumder et al., “Polarization resolved second harmonic microscopy,” Methods, 128 105 –118 (2017). MTHDE9 1046-2023 Google Scholar


G. Hall et al., “Goniometric measurements of thick tissue using Monte Carlo simulations to obtain the single scattering anisotropy coefficient,” Biomed. Opt. Express, 3 2707 –2719 (2012). BOEICL 2156-7085 Google Scholar


R. Lacomb, O. Nadiarnykh and P. J. Campagnola, “Quantitative second harmonic generation imaging of the diseased state osteogenesis imperfecta: experiment and simulation,” Biophys. J., 94 4504 –4514 (2008). BIOJAU 0006-3495 Google Scholar


R. Lacomb et al., “Phase matching considerations in second harmonic generation from tissues: effects on emission directionality, conversion efficiency and observed morphology,” Opt. Commun., 281 1823 –1832 (2008). OPCOB8 0030-4018 Google Scholar


G. Hall, K. W. Eliceiri and P. J. Campagnola, “Simultaneous determination of the second-harmonic generation emission directionality and reduced scattering coefficient from three-dimensional imaging of thick tissues,” J. Biomed. Opt., 18 116008 (2013). JBOPFO 1083-3668 Google Scholar


L. Wang, S. L. Jacques and L. Zheng, “MCML—Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed., 47 131 –146 (1995). CMPBEK 0169-2607 Google Scholar


K. R. Campbell and P. J. Campagnola, “Assessing local stromal alterations in human ovarian cancer subtypes via second harmonic generation microscopy and analysis,” J. Biomed. Opt., 22 116008 (2017). JBOPFO 1083-3668 Google Scholar


W. McKleroy, T. H. Lee and K. Atabai, “Always cleave up your mess: targeting collagen degradation to treat tissue fibrosis,” Am. J. Physiol. Lung Cell Mol. Physiol., 304 L709 –L721 (2013). APLPE7 1040-0605 Google Scholar


K. Zhang et al., “Myofibroblasts and their role in lung collagen gene expression during pulmonary fibrosis. A combined immunohistochemical and in situ hybridization study,” Am. J. Pathol., 145 114 –125 (1994). AJPAA4 0002-9440 Google Scholar


J. D. Rogers, I. R. Capoglu and V. Backman, “Nonscalar elastic light scattering from continuous random media in the Born approximation,” Opt. Lett., 34 1891 –1893 (2009). OPLEDP 0146-9592 Google Scholar


V. Turzhitsky et al., “A predictive model of backscattering at subdiffusion length scales,” Biomed. Opt. Express, 1 1034 –1046 (2010). BOEICL 2156-7085 Google Scholar


P. H. Saldiva et al., “Histochemical evaluation of lung collagen content in acute and chronic interstitial diseases,” Chest, 95 953 –957 (1989). CHETBF 0012-3692 Google Scholar


A. J. Hance and R. G. Crystal, “The connective tissue of lung,” Am. Rev. Respir. Dis., 112 657 –711 (1975). ARDSBL 0003-0805 Google Scholar


G. F. Rozin et al., “Collagen and elastic system in the remodelling process of major types of idiopathic interstitial pneumonias (IIP),” Histopathology, 46 413 –421 (2005). HISTDD 1365-2559 Google Scholar


R. Z. Kramer et al., “Sequence dependent conformational variations of collagen triple-helical structure,” Nat. Struct. Biol., 6 454 –457 (1999). NSBIEW 1072-8368 Google Scholar


P. J. Campagnola and L. M. Loew, “Second-harmonic imaging microscopy for visualizing biomolecular arrays in cells, tissues and organisms,” Nat. Biotechnol., 21 1356 –1360 (2003). NABIF9 1087-0156 Google Scholar


M. G. Jones et al., “Nanoscale dysregulation of collagen structure-function disrupts mechano-homeostasis and mediates pulmonary fibrosis,” Elife, 7 e36354 (2018). Google Scholar


E. E. Torr et al., “Myofibroblasts exhibit enhanced fibronectin assembly that is intrinsic to their contractile phenotype,” J. Biol. Chem., 290 6951 –6961 (2015). JBCHA3 0021-9258 Google Scholar


C. Ricciardelli and R. J. Rodgers, “Extracellular matrix of ovarian tumors,” Semin. Reprod. Med., 24 270 –282 (2006). Google Scholar


F. J. Martinez et al., “Idiopathic pulmonary fibrosis,” Nat. Rev. Dis. Primers, 3 17074 (2017). Google Scholar


V. J. Craig et al., “Matrix metalloproteinases as therapeutic targets for idiopathic pulmonary fibrosis,” Am. J. Respir. Cell Mol Biol., 53 585 –600 (2015). Google Scholar


Y. C. Wu et al., “Scanning all-fiber-optic endomicroscopy system for 3-D nonlinear optical imaging of biological tissues,” Opt. Express, 17 7907 –7915 (2009). OPEXFF 1094-4087 Google Scholar


Darian S. James is currently a PhD student in the Biomedical Engineering Department at the University of Wisconsin.

Alexander N. Jambor is currently a research intern in the Biomedical Engineering Department at the University of Wisconsin.

Hsin-Yu Chang was a research intern in the Biomedical Engineering Department at the University of Wisconsin.

Zachary Alden was a BS student in the Biomedical Engineering Department at the University of Wisconsin.

Karissa B. Tilbury was a PhD student in the Biomedical Engineering Department at the University of Wisconsin. She is now an assistant professor at the University of Maine.

Nathan K. Sandbo is an associate professor of medicine at the University of Wisconsin. A medical doctor (MD), he specializes in research and clinical care of idiopathic pulmonary fibrosis.

Paul J. Campagnola is a professor in the Biomedical Engineering and Medical Physics Departments at the University of Wisconsin. He specializes in developing nonlinear optical approaches to study changes in the ECM in diseased states.

© The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Darian S. James, Alexander N. Jambor, Hsin-Yu Chang, Zachary Alden, Karissa B. Tilbury, Nathan K. Sandbo, and Paul J. Campagnola "Probing ECM remodeling in idiopathic pulmonary fibrosis via second harmonic generation microscopy analysis of macro/supramolecular collagen structure," Journal of Biomedical Optics 25(1), 014505 (29 November 2019).
Received: 12 September 2019; Accepted: 6 November 2019; Published: 29 November 2019

Back to Top