Lysyl oxidases regulate fibrillar collagen remodelling in idiopathic pulmonary fibrosis

Abstract

Idiopathic pulmonary fibrosis (IPF) is a progressive scarring disease of the lung with few effective therapeutic options. Structural remodelling of the extracellular matrix [i.e. collagen cross-linking mediated by the lysyl oxidase (LO) family of enzymes (LOX, LOXL1-4)] might contribute to disease pathogenesis and represent a therapeutic target. This study aimed to further our understanding of the mechanisms by which LO inhibitors might improve lung fibrosis. Lung tissues from IPF and non-IPF subjects were examined for collagen structure (second harmonic generation imaging) and LO gene (microarray analysis) and protein (immunohistochemistry and western blotting) levels. Functional effects (collagen structure and tissue stiffness using atomic force microscopy) of LO inhibitors on collagen remodelling were examined in two models, collagen hydrogels and decellularized human lung matrices. LOXL1/LOXL2 gene expression and protein levels were increased in IPF versus non-IPF. Increased collagen fibril thickness in IPF versus non-IPF lung tissues correlated with increased LOXL1/LOXL2, and decreased LOX, protein expression. β-Aminoproprionitrile (β-APN; pan-LO inhibitor) but not Compound A (LOXL2-specific inhibitor) interfered with transforming growth factor-β-induced collagen remodelling in both models. The β-APN treatment group was tested further, and β-APN was found to interfere with stiffening in the decellularized matrix model. LOXL1 activity might drive collagen remodelling in IPF lungs. The interrelationship between collagen structural remodelling and LOs is disrupted in IPF lungs. Inhibition of LO activity alleviates fibrosis by limiting fibrillar collagen cross-linking, thereby potentially impeding the formation of a pathological microenvironment in IPF.

INTRODUCTION

Idiopathic pulmonary fibrosis (IPF) is a specific form of progressive fibrosing lung disease of unknown cause that leads to chronic respiratory failure and death (Raghu et al., 2011), with a median survival of ∼3 years after diagnosis (Demedts et al., 2001). IPF is a subcategory under the general disease category of interstitial lung disease (ILD). Common traits of ILDs are the presence of inflammation and pulmonary fibrosis. IPF patients have a poor prognosis, with 44% mortality rate for IPF patients expected within 5 years, compared with the mortality rates in two other ILDs, ILD-associated connective tissue disease (33%) and pulmonary sarcoidosis (2%) (Demedts et al., 2001).

Like many fibrotic disorders, IPF is characterized by enhanced deposition and remodelling of the extracellular matrix (ECM). Although data are emerging about the altered composition of the ECM in IPF (Booth et al., 2012; Estany et al., 2014; Kuhn et al., 1989; Kuhn and McDonald, 1991; Ward and Hunninghake, 1998), less is known about the structural changes that occur in collagens in this disease. Collagens undergo extensive post-translational modifications, including formation of cross-links between fibres that stabilize the molecules and generate a strong collagen fibril, which provides tensile strength. Aberrant cross-linking is known to contribute to increased tissue stiffness, a pathological characteristic of IPF (Clarke et al., 2013).

Fibrillar collagen structures (predominantly collagens I and III) can be visualized and quantified using second harmonic generation (SHG) microscopy (Abraham et al., 2012; Tjin et al., 2014). SHG imaging has been used to analyse ECM structures in biological systems in vitro, ex vivo and in vivo (Brown et al., 2003; Caetano-Lopes et al., 2010; Cox et al., 2003; Guo et al., 1999; Nadiarnykh et al., 2007; Sivaguru et al., 2010; Suzuki et al., 2012; Williams et al., 2001) and can semi-quantitatively measure the ratio of mature organized collagen to immature disorganized collagen fibrils (Williams et al., 2005). This technique, therefore, enables comparisons of collagen rearrangement in tissues, analysis of collagen structure (Caetano-Lopes et al., 2010; Nadiarnykh et al., 2007; Sivaguru et al., 2010; Suzuki et al., 2012) and even measurement of the thickness of the fibres (Chu et al., 2007). We have recently reported the robustness of this technique to quantify fibrillar collagen I remodelling in airway tissues in chronic obstructive pulmonary disease (Tjin et al., 2014).

Maturity of collagen fibres and thus the mechanical properties of the fibres are dependent on the level of intermolecular cross-linking. Although the intermolecular fibrillar organization and certain cross-linking processes can happen spontaneously, enzymes such as lysyl oxidases (LOs) are necessary for further cross-linking. The LO family of enzymes consists of five known paralogues: lysyl oxidase (LOX) and LOX like 1-4 (LOXL1-4). During normal development, the expression of LOs is tightly regulated, with aberrant expression and activity being associated with pathological manifestation in various diseases, including cancer (Erler et al., 2006; Payne et al., 2005) and, as recently discovered, IPF (Aumiller et al., 2017; Barry-Hamilton et al., 2010).

Collagen cross-linking contributes to the stiffness of the ECM of tissues and is tightly regulated in normal tissue homeostasis. However, in diseased tissue, aberrant cross-linking could lead to the development of a pathological microenvironment. The stiffness of the ECM influences the behaviour and function of cells within its vicinity, including an increase in fibroblast proliferation and contraction (Balestrini et al., 2012; Marinkovic et al., 2013). In addition, a stiffer ECM leads to an increase in latent transforming growth factor-β (TGF-β) activation (Burgess et al., 2016; Shi et al., 2011; Wipff et al., 2007). In IPF, the lung tissue has a stiffness of 16.52±2.25 kPa, as opposed to 1.96±0.13 kPa in healthy lung tissue (Booth et al., 2012). A stiffer matrix also promotes IPF-like phenotypes, such as enhanced differentiation, proliferation and resistance to apoptosis in non-IPF fibroblasts (Liu et al., 2010; Marinkovic et al., 2012; Paszek et al., 2005; Wipff et al., 2007). Despite the availability of two recently approved drugs for slowing the progress of IPF [nintedanib (Mazzei et al., 2015; Richeldi et al., 2014a,b) and pirfenidone (King et al., 2014; Valeyre et al., 2014)], their associated incidence of side effects and modest clinical benefit warrants further investigation into possible new therapeutic directions for IPF (Moodley et al., 2015). Selective inhibition of LOXL2 activity was recently investigated as a therapeutic approach for IPF [Barry-Hamilton et al., 2010; ClinicalTrials.gov NCT01769196 (Raghu et al., 2017)]. Although this trial was terminated early because of lack of efficacy, there is still potential for the development of effective LO inhibitors for targeting fibrosis. To enhance understanding of the mechanisms by which LOX and LOXL enzyme inhibition might be beneficial in IPF, this study aimed to investigate fibrillar collagen structural remodelling in IPF lung tissues and to determine the expression profile of LOs in IPF. It also investigated the effectiveness of pan-LO and selective LOXL2 inhibition on fibrillar collagen structural remodelling and tissue stiffness in the context of IPF.

RESULTS

Increased fibrillar collagen maturity/organization in IPF lung tissues

To investigate the collagen structural remodelling in IPF, SHG microscopy was performed on lung tissue sections from non-IPF and IPF subjects. The SHG forward (F) signal is predominantly from mature/organized fibrillar collagen, whereas the backward (B) signal is from immature/disorganized fibrillar collagen. Thus, increased F/B area corresponds to increased amounts of mature/immature collagen, whereas increased F/B intensity ratios correspond to an increased degree of collagen maturity/thickness and vice versa (Tjin et al., 2014). Quantification of SHG images showed that IPF lung parenchymal tissue had greater levels of fibrillar collagen maturity/organization compared with the non-IPF tissues (Fig. 1A,B). By contrast, the collagen structural organization was not different between IPF and non-IPF ILD tissues (Fig. S1).

Fig. 1. Increased level of collagen organization/maturity in IPF lung parenchyma. (A) Representative SHG images of lung parenchyma from non-IPF (n=7; open circles) and IPF (n=8; filled diamonds) subjects (average of three or four samples per subject; yellow indicates backward immature/disorganized collagen; cyan indicates mature/organized collagen; scale bar: 100 [mu]m). Brightness and contrast have been enhanced for display purposes only and were applied equally for all images. (B) Image analysis quantification of SHG intensity forward/backward signal ratio (intensity F/B; each point represents the mean of three or four lung parenchymal samples per subject); a higher ratio indicates higher proportion of organized/mature fibrillar collagen. Data were analysed with Student's two-tailed parametric unpaired t-test. Data are presented as means±s.d. Data were obtained from one experimental replicate per subject. *P[ less than ]0.05.

Lysyl oxidase genes were differentially expressed in IPF and healthy control lung tissues

Fibrillar collagen maturity is catalysed by the LO family of enzymes. To investigate the levels of LO enzymes in lung tissues, lung tissue gene expression profiles of IPF and healthy controls were analysed from two publically available data sets (Table 1). LOXL1 gene expression was increased in IPF lung tissue compared with healthy controls in two independent studies [false discovery rate (FDR) adjusted P-value < 0.25], whereas the LOXL2 gene expression was increased in IPF lung tissue in a single data set. No differences were found in LOX gene expression levels.

Table 1. LOX expression levels in lung tissue

Lysyl oxidase protein levels were differentially expressed between IPF and non-diseased lung tissues

To confirm the findings from the gene expression study, tissue sections from the same samples used for the examination of collagen structural organization were examined for LOX, LOXL1 and LOXL2 protein levels. Quantification of the captured images of the stained lung tissue sections showed that LOX percentage surface area and density were decreased in IPF compared with non-IPF, whereas the percentage surface area and density of LOXL1 was increased in IPF compared with non-IPF tissues (Table 2; representative images in Fig. 2A). In agreement with prior studies (Aumiller et al., 2017; Barry-Hamilton et al., 2010), LOXL2 density was also increased in IPF compared with non-IPF tissues. By contrast, only LOXL1, but not LOX or LOXL2, density increased in IPF compared with non-IPF ILD (Table S1). In vitro, primary human lung parenchymal fibroblasts from IPF patients produced significantly greater levels of LOXL2 but not LOX or LOXL1 in the cell lysate compared with non-IPF controls (Fig. S2A,B).

Table 2. Lysyl oxidase enzymes are differentially expressed in IPF lung tissues

Fig. 2. Collagen organization and lysyl oxidase family of enzymes are differentially expressed and interrelate differently in IPF. (A) Representative IHC images of lysyl oxidase enzymes in the lungs. Images are representative of lung tissues from non-IPF (n=10) and IPF (n=8) subjects (brown, protein of interest; blue, nucleus; pink, cytoplasm; scale bar: 200 [mu]m). (B) The ratio of LOX and LOXL1 percentage surface area obtained by image analysis of IHC-stained lung tissues of non-IPF (n=10; open circles) compared with IPF (n=8; filled diamonds) subjects; data were analysed with Student's two-tailed parametric unpaired t-test. (C) Correlation of the ratio of LOX and LOXL1 with F/B ratio obtained through image analysis of SHG microscopy images of collagen structure in non-IPF (n=7) and IPF (n=8) subjects (P[ less than ]0.0001; R2=0.7527; average of three or four samples per subject). (D) 3D scatter plot of collagen organization (intensity F/B, i.e. SHG intensity forward/backward signal ratio), LOXL1 and LOXL2 expression (open circles, non-IPF, n=7; filled circles, IPF, n=8). (E) Scree plot showing the number of relationship groups needed to explain the associations within measured variables in non-IPF and IPF tissues (eigenvalue threshold≥1). (F,G) Weighted average values of the variables contributing to each relationship group; comparing non-IPF with IPF subjects (average of three or four samples per subject); data were analysed with Student's two-tailed parametric unpaired t-test. Data are presented as means±s.d. Data were obtained from one experimental replicate per subject. *P[ less than ]0.05, **P[ less than ]0.01, ***P[ less than ]0.005.

Expression of LOs was correlated with collagen remodelling in lung tissues

To begin to understand the relationship between LOs and collagen structural arrangement in lung tissues, we used correlation analysis. The data from both non-IPF and IPF subjects were combined in the correlation studies to examine the whole spectrum of collagen remodelling in both cohorts. There was a positive correlation between both LOXL1 (Table 3; Fig. S3) and LOXL2 (Table 3) with collagen remodelling (SHG F/B ratios), whereas there was a negative correlation between LOX and SHG forward to backward signal (F/B) ratios (Table 3). LOXL1 positively correlated with LOXL2, whereas LOX negatively correlated with both LOXL1 and LOXL2 (Table 3).

Table 3. Expression of lysyl oxidase enzymes correlates with the maturity/organization of the fibrillar collagen and other lysyl oxidase enzymes in lung tissues

Lysyl oxidase family members were differentially expressed in IPF

The ratio of LOX/LOXL1 (Fig. 2B) and LOX/LOXL2 (Fig. S4) densities and LOX/LOXL1 percentage surface areas were significantly decreased in IPF lung tissue samples compared with non-IPF subjects (Fig. 2B; Fig. S4). The ratio of LOX/LOXL1 percentage surface area inversely correlated with collagen maturity/organization (F/B; Fig. 2C). There was no significant difference in the ratios of LO family enzymes in the tissues of IPF and non-IPF ILD subjects (Fig. S5; P=NS).

Collagen organization and LOX enzymes interrelate differently in IPF

Given that LOX, LOXL1 and LOXL2 were differentially expressed in IPF compared with non-IPF (Fig. 2B,C) and the expression of these enzymes also correlated strongly with collagen remodelling (Table 3), it was important to understand the interrelationship between these parameters in IPF. Visualization of the three-dimensional (3D) association among these variables showed spatial separations between the variables in IPF compared with non-IPF (Fig. 2D). We next performed factor analysis to investigate the number of relationship groups (factors) accounting for the structure of these associations. A minimum of two groups were required to explain the associations within the data set (Fig. 2E), and Table S2 shows the strength of the contribution of the variables [LOX, LOXL1, LOXL2, immunohistochemistry (IHC) data and SHG data] to the two relationship groups. The strengths of contributions of the variables were used to determine the weighted average values for each relationship group (Fig. 2F,G). The values for relationship Groups 1 and 2 were significantly different between non-IPF and IPF, indicating a difference in the relationships between the variables (LOX, LOXL1, LOXL2 IHC data and SHG data) in disease.

Inhibition of LO activity reduced TGF-β-induced collagen remodelling in collagen I hydrogels

To confirm the role of LO enzymes in fibrillar collagen remodelling, 3D in vitro cell culture experiments using inhibitors of LO activity were performed. β-APN is a non-selective inhibitor of LO activity, whereas Compound A is a specific inhibitor of LOXL2 activity (Schilter et al., 2014). Pan-inhibition of LO activity, but not inhibition of LOXL2-specific activity, reduced TGF-β-induced collagen remodelling in collagen hydrogel fibroblast cultures (Fig. 3). TGF-β increased the deposition of mature fibrillar collagen by non-IPF fibroblasts but not IPF fibroblasts, which was alleviated by the pan-LO inhibitor β-APN (Fig. 3D), but TGF-β did not affect collagen fibre maturity/thickness (Fig. 3E). By contrast, Compound A did not alleviate the increased deposition of mature collagen compared with immature collagen, from both fibroblast groups, induced by TGF-β, but on the contrary, increased it (Fig. 3F). Compound A did reduce the thickness of the fibres generated by the IPF fibroblasts in the presence and absence of TGF-β stimulation (Fig. 3G).

Fig. 3. Inhibition of lysyl oxidase activity reduces TGF-[beta]-induced collagen remodelling. Collagen I hydrogels seeded with primary human lung parenchymal fibroblasts (250,000 cells ml−1) were stimulated with TGF-[beta] (10 ng ml−1) to induce a fibrotic phenotype in the presence or absence of the pan-lysyl oxidase inhibitor, [beta]-APN (100 [mu]M), or the LOXL2-specific inhibitor, Compound A (Cmpd A; 300 nM). Samples were treated for 7 days. At the end of the treatment period, the samples were formalin fixed and paraffin embedded. (A) Haematoxylin and Eosin stain of fibroblasts in collagen gel (blue, nucleus; pink, cytoplasm and general tissue stain; scale bar: 100 [mu]m) at day 7. (B) SHG image of collagen gel with fibroblasts at day 7 (cyan, forward mature/organized collagen; magenta, cell autofluorescence; scale bar=100 [mu]m). (C) Representative SHG images of collagen gels with primary human lung parenchymal fibroblasts treated with or without TGF-[beta] (10 ng ml−1) in the presence or absence of Compound A (300 nM) (yellow, backward immature/disorganized collagen; cyan, forward mature/organized collagen; scale bar: 20 [mu]m). (D-G) Quantification of SHG microscopy on collagen I hydrogels embedded with fibroblasts from non-IPF (n=5; open circles) and IPF (n=5; filled diamonds) subjects treated with pan-lysyl oxidase inhibitor ([beta]-APN; D,E) or LOXL2-specific inhibitor (Compound A; F,G). Data are presented as a percentage of the respective no-cell control (means±s.d.). (D,F) SHG forward/backward signal surface area ratio (F/B area); higher ratio indicates higher amount of mature compared with immature collagens. (E,G) SHG forward/backward signal intensity ratio (F/B intensity); higher ratio indicates greater maturity/thickness level of fibrillar collagen fibres. Statistical analysis was by two-way ANOVA with matching and multiple comparisons with Tukey's correction. Data are presented as means±s.d. Data were obtained from one experimental replicate per subject. Brightness and contrast have been enhanced for display purposes only and were applied equally for all images. *Significant difference compared with BSA control; [delta]significant difference compared with BSA+Compound A control; [phi]significant difference compared with TGF-[beta] alone; and #significant difference between groups.

Inhibition of LO activity reduced TGF-β-induced collagen remodelling in reseeded decellularized matrices

To model the in vivo cell ECM interactions more closely, we used decellularized human IPF and non-IPF lung scaffolds. There were significantly more primary human lung fibroblasts remaining in the supernatant after reseeding of the IPF matrices compared with the non-IPF matrices (Fig. 4A), potentially indicating higher cell attachment to the non-IPF matrix. This was supported by histological staining of the reseeded matrices, which showed that fewer cells were detectable in the IPF matrix post-reseeding, compared with the non-IPF matrix (Fig. 4B). Given the sparsity of cellular attachment in the IPF matrices, these matrices were not examined further.

Fig. 4. Inhibition of lysyl oxidase activity reduces TGF-[beta]-induced collagen remodelling. Decellularized non-IPF matrices reseeded with primary human lung parenchymal fibroblasts (1×106 cells per matrix) were stimulated with TGF-[beta] (10 ng ml−1) to induce a fibrotic phenotype in the presence or absence of the pan-lysyl oxidase inhibitor, [beta]-APN (100 [mu]M), or the LOXL2-specific inhibitor, Compound A (Cmpd A; 300 nM). Cells were seeded in the matrix for 24 h before being transferred into separate tissue culture wells with treatment media for 7 days. Samples were collected at the end of the treatment period and formalin fixed and paraffin embedded. (A) Manual cell counts of cells remaining in supernatant after reseeding in matrices [n=9; combined data from 4 non-IPF and 5 IPF fibroblasts (500,000 cells ml−1; 2 ml per matrix) into matrices from one non-IPF and one IPF subject]. (B) Masson's Trichrome stain of fibroblasts reseeded into decellularized matrices for 7 days (blue, collagen; pink, cell cytoplasm; scale bar: 100 [mu]m). (C) Representative SHG images of collagen gels with fibroblasts treated with or without TGF-[beta] (10 ng ml−1) in the presence of or absence of Compound A (300 nM) (yellow, backward immature/disorganized collagen; cyan, forward mature/organized collagen; scale bar: 100 [mu]m). (D) Quantification of SHG microscopy on decellularized non-IPF matrices reseeded for 7 days with fibroblasts from control (n=5; open circles) and IPF (n=5; filled diamonds) subjects treated with [beta]-APN (100 [mu]M) or Compound A (300 nM) in the presence or absence of TGF-[beta] (10 ng ml−1). Data are presented as a percentage of the respective no-cell control (means±s.d.); SHG forward/backward signal intensity ratio (F/B intensity); higher ratio indicates greater maturity/thickness level of fibrillar collagen fibres. Statistical analysis was by two-way ANOVA with matching and multiple comparisons with Tukey's correction. Data are presented as means±s.d. Data were obtained from one experimental replicate per subject. Brightness and contrast have been enhanced for display purposes only and were applied equally for all images. *Significant difference compared with BSA control; and #significant difference between groups.

The non-IPF cells increased the maturity/thickness of the collagen fibres in the non-IPF matrix in the presence of TGF-β (Fig. 4C,D) but did not increase the amount of mature collagen compared with immature collagen (Fig. S6). The pan-inhibitor of LOs, β-APN, returned TGF-β-induced collagen remodelling to baseline. In addition, the LOXL2-specific inhibitor blocked TGF-β-induced remodelling of fibrillar collagen; remodelling levels in the presence of either inhibitor without TGF-β were not different compared with control (Fig. 4D). Primary IPF lung fibroblasts seeded into decellularized normal matrices did not remodel the fibrillar collagen in the presence of TGF-β (Fig. 4D).

Inhibition of LO activity reduced TGF-β-induced stiffening in reseeded decellularized matrices

To confirm the functional significance of TGF-β-induced collagen remodelling in the decellularized matrices reseeded with non-IPF cells, tissue stiffness (elastic modulus) was measured by atomic force microscopy (AFM; micro-indentation). TGF-β induced a significant increase in tissue stiffness that was inhibited by the pan-inhibitor of LOs, β-APN (Fig. 5). Individual stiffness measurements for each patient segregated by treatment can be found in Fig. S7, and representative force curves used to measure stiffness can be found in Fig. S8C.

Fig. 5. Inhibition of lysyl oxidase activity reduces TGF-[beta]-induced stiffening of reseeded decellularized matrices. Decellularized non-IPF matrices reseeded with primary human lung parenchymal fibroblasts (106 cells per matrix) were stimulated with TGF-[beta] (10 ng ml−1) in the presence or absence of the pan-lysyl oxidase inhibitor, [beta]-APN (100 [mu]M). Cells were seeded in the matrix for 24 h before being transferred into separate tissue culture wells with treatment media for 7 days. Samples were collected at the end of the treatment period and fresh frozen in OCT compound (10 [mu]m slices). Box and whiskers plot of the stiffness (elastic modulus). Stiffness was measured with AFM, and > 100 measurements spread over three areas were taken per sample. Each n is a single AFM measurement. Data were obtained from one experimental replicate per subject per treatment. Statistical analysis was by the Kruskal–Wallis test with Dunn's correction. **P[ less than ]0.01, ****P[ less than ]0.0001.

DISCUSSION

Although increased collagen deposition is well appreciated in fibrotic lung tissues, this study identifies, for the first time, alterations in the fibrillar organization of the collagen fibres in IPF lung tissue. In addition, this study reports that LOXL1 and LOXL2 gene and protein levels are increased in the lungs of IPF compared with non-IPF subjects, with a positive correlation between increased protein expression and increased fibrillar collagen organization.

These findings show the complex interrelationship between LOs and collagen structure and provide further support for focusing on LO inhibition as a potential therapeutic approach in IPF. Indeed, LOXL2 inhibition prevented and even reversed bleomycin-induced pulmonary fibrosis in mice (Barry-Hamilton et al., 2010). In addition, LOX and LOXL2 were also found to be elevated in liver fibrosis (Barry-Hamilton et al., 2010; Vadasz et al., 2005) and were associated with increased deposition of collagen around the hepatocytes (Vadasz et al., 2005). Although the recent clinical trial targeting LOXL2 inhibition (ClinicalTrials.gov NCT01769196) was halted early because of a lack of efficacy (Raghu et al., 2017), the present study suggests that targeting LOXL1 in IPF might be as important as targeting LOXL2.

In the present study, IPF alveoli had higher collagen maturity/thickness (increased F/B SHG) compared with non-diseased subjects. It can be postulated that the increased collagen fibre maturity/thickness contributes to the increased stiffness characteristic of IPF lungs (Booth et al., 2012). This is supported by the increased stiffnesses measured in this study in the decellularized lung scaffolds reseeded with fibroblasts exposed to TGF-β. In IPF there is augmented TGF-β (Broekelmann et al., 1991), which induces IPF fibroblasts to produce enhanced ECM proteins (Westergren-Thorsson et al., 2004), thereby creating a positive feedback loop. In fact, collagen cross-linking might account for organ stiffening in a variety of fibrotic diseases. The results of the present study are also concordant with those of Kottmann et al. (2015), who demonstrated an altered F/B SHG signal in IPF compared with non-diseased lungs. Interestingly, the study by Kottmann et al. (2015) showed a decrease in F/B SHG between healthy and IPF lungs; this difference might be attributable to different sampling and analysis protocols in comparison to the present study. For example, in the study by Kottmann et al. (2015), the majority of the F SHG signal (numerator in the F/B SHG signal ratio) from healthy subjects originated from muscular arteries, whereas the present study excluded such regions from analysis. That said, neither study identified differences in collagen remodelling between IPF and other fibrotic interstitial lung diseases.

A potential role of other LO family members, not only LOXL2, in collagen remodelling in IPF was also highlighted in the present study. A significant increase in LOXL1 gene and protein expression in IPF lungs, compared with non-IPF tissues, was observed, with a positive correlation of protein levels with collagen organization. By contrast, LOX gene expression was unchanged, whereas the LOX protein level was decreased in IPF lung, highlighting differential LO family expression in IPF lung tissue, which is similar to that reported during tumour progression and metastasis (Erler et al., 2006; Payne et al., 2005). In concordance with our findings, a recent study by Aumiller et al. (2017) also reported increased LOXL1 and LOXL2 protein levels in IPF lungs compared with healthy lungs. By contrast, however, they also reported an increased LOX protein level in IPF (Aumiller et al., 2017). These differential results might be attributable to the different analysis approaches. Aumiller et al. (2017) scored the intensity of LO staining in discrete structures of the lung, whereas in our study we used digitized image analysis to measure the staining intensity of the whole lung tissue section whilst controlling for the increased tissue mass in IPF lungs.

Of the LO family enzymes, only the LOXL1 protein level was increased in IPF compared with non-IPF ILD lungs, indicating that changes in LOXL1 might be specific to IPF rather than to the process of fibrosis itself. To our knowledge, this study is the first to show a difference in the interrelationships of the organization of collagen, LOX, LOXL1 and LOXL2 enzymes between the lungs of IPF and non-IPF subjects. This is particularly important, because a Phase II clinical trial using simtuzumab, a LOXL2-specific inhibitor [ClinicalTrials.gov NCT01769196 (Raghu et al., 2017)] was recently terminated early owing to a lack of efficacy. The data presented in the present study support the continued development of new therapeutic strategies targeting LOs for IPF, especially focusing on exploration of the potential of LOXL1-specific inhibitors or combination approaches.

To understand the role of LOs in the remodelling of collagen fibres in IPF, two 3D in vitro culture models were used to mimic the complexity of the matrix in vivo. Interestingly, the decellularized IPF matrices supported lower primary fibroblast infiltration and viability compared with non-IPF matrices. Although prior works (Booth et al., 2012; Parker et al., 2014) are in line with the present data that decellularized IPF matrices influence cellular phenotype, previous evaluations had not examined differences in cell attachment between non-IPF and IPF matrices. Thus, these observations highlight another difference between IPF and non-IPF lung matrices that should be investigated as a contributor to the pathogenesis of fibrosis. These findings further highlight the challenge ahead in finding solutions for repairing the disrupted ECM in IPF.

In collagen gels, TGF-β induced fibroblast production of mature/thick collagen fibres, but did not affect the overall maturity/thickness of the fibres. By contrast, TGF-β did not significantly alter the amount of mature collagen fibres in decellularized human lung matrices, whereas overall maturity/thickness was increased. This apparent paradox might simply be attributable to the model compositions, in that collagen hydrogels lack the complex ECM components found in decellularized matrices, components that might influence cellular responsiveness. Alternatively, the artificiality of the collagen hydrogel, where collagen fibres are not in a native structural configuration, leads cells to respond by depositing and remodelling collagen fibres to attain a native environment. By contrast, decellularized matrices are composed of collagen fibres and ECM already in a native configuration, thereby driving cellular driven ECM remodelling rather than production. This latter conclusion is supported by the presented data (Figs 4 and 5) showing collagen gel F/B intensity > 150% of the no cell control even in unstimulated conditions, whereas decellularized matrices were recorded at ∼100% in the same conditions. The differences in the presence of artificial matrix (collagen hydrogel) and native matrix (decellularized matrix) support previous findings that the ECM contributes a very important regulatory control on cellular behaviour (Booth et al., 2012; Parker et al., 2014).

Pan-inhibition of LO activity blocked TGF-β-induced collagen remodelling in both collagen hydrogels and decellularized matrix systems. By contrast, LOXL2-specific inhibition reduced collagen fibril maturity/thickness only in IPF fibroblasts and in decellularized matrices. As this effect was not seen with the pan-inhibitor of LOs, it is likely that this effect was independent of LOXL2. However, prior data demonstrate that fibroblasts grown in two-dimensional (2D) culture express a 50 kDa proteolytically processed form of LOXL2 (Hollosi et al., 2009), rather than the intact 95 kDa LOXL2 enzyme that is inhibited by Compound A. Although the functional activity of the LOXL2 50 kDa isoform is unknown, it seems unlikely to account for the present observations, because Compound A was able to block other TGF-β-induced effects in 2D and 3D culture. Additionally, pan-inhibition of LO activity blocked TGF-β-induced tissue stiffening in the decellularized matrix system (Fig. 5). Both collagen fibril maturity/thickness and tissue stiffness were increased by TGF-β treatment, with both being interrupted in the presence of β-APN. These data indicate that LO enzymes regulate the TGF-β-induced increase in collagen fibril maturity/thickness that results in increased tissue stiffness. The pro-fibrotic features of stiff matrices contribute to the TGF-β-induced fibrosis with potential feedback mechanisms through the promotion of latent TGF-β activation (Wipff et al., 2007), thus contributing to the fibrotic process in IPF.

It is worth noting that the stiffness measured in the reseeded decellularized matrices (mean 40-60 kPa) was higher than the values measured in native human lung parenchymal tissue (1.34±0.36 kPa; Liu et al., 2010) and IPF tissue (16.52±2.25 kPa; Booth et al., 2012). It has been reported in the literature that the decellularization method affects the stiffness of acellular tissue (Melo et al., 2014) and, more generally, that decellularized tissue is stiffer than normal lung tissue (Giménez et al., 2017), although the origins of these differences in decellularized versus normal tissue stiffness remain poorly understood. The reseeded decellularized matrices did not form alveolar structures in the time span of this study and were densely packed (Fig. 4B), possibly as a result of lack of inflation of the tissues, which is not possible within the present experimental set-up. This tissue structure might more closely mimic scar tissue or remodelled IPF lung tissue (characterized by progressive scarring of the lungs) rather than non-diseased lung tissues. In such a case, this would strengthen the significance of reseeded decellularized matrices in the study of IPF tissue remodelling. It should also be noted that β-APN might have non-cell-induced effects on the stiffness, as the stiffness appears lower in β-APN-treated non-reseeded matrices than those without β-APN treatment (∼40 versus ∼60 kPa, respectively; Fig. S7), but the present experiment was not designed to address this hypothesis.

Overall, this study highlights a potential mechanism for preventing ongoing fibrotic processes, with LO inhibition preventing lung fibrogenesis in IPF by blocking fibrillar collagen organization and subsequent tissue stiffening. This, in turn, might leave deposited collagen more susceptible to proteolytic degradation and clearance, although further research is needed to investigate this possibility. The present findings also highlight the role of multiple members of the LO family in regulating collagen cross-linking and the disrupted balance of LOs in the IPF lung. These findings indicate that selective inhibition of LO enzymes, such as LOXL1, either alone or in combination with LOXL2 inhibition, might be a potential therapeutic strategy for targeting fibrosis in IPF.

MATERIALS AND METHODS

Human tissues

Human lungs were procured with written and informed consent from the patients or next of kin from non-transplantable donors and IPF patients undergoing lung transplantation as previously described (Booth et al., 2012). Approval for all experiments with human lung tissue was provided by the Ethics Review Committees of the South West Sydney Area Health Service, St Vincent's Hospital Sydney, Royal Prince Alfred Hospital (RPAH), and the University of Sydney Human Research Ethics Committee. As tissues were de-identified, the University of Michigan Institutional Review Board deemed this work exempt from oversight. Primary human lung fibroblast samples used as non-diseased controls (non-IPF) were from the macroscopically normal regions of lungs from patients with non-small cell carcinoma (NSCCa). Human lung tissue samples were classified by pathologist diagnosis. Patient demographics are described in Table S3.

Tissue and cell isolation

Tissue isolation was performed on fresh explanted lung tissues. Human lung parenchymal tissues were isolated and cut into pieces no bigger than 1 cm³. These pieces were fixed in 4% formaldehyde (Fronine; Thermo Fisher Scientific, Waltham, MA, USA). After formaldehyde fixation, isolated lung tissues were embedded in paraffin.

Primary human lung fibroblasts were isolated as previously described (Krimmer et al., 2012; White et al., 2003) and were grown in growth media [Dulbecco's modified Eagle's medium (DMEM) low glucose (Gibco, Thermo Fisher Scientific) with 0.025 M HEPES (Gibco), 0.375% sodium hydrogen carbonate (Ajax Finechem, Thermo Fisher Scientific), 10% (v/v) FBS (JRH Biosciences, Brooklyn, Victoria, Australia) and 1% (v/v) antibiotics (Invivogen, San Diego, CA, USA), pH 7.1-7.2]. All fibroblast cell lines were tested for mycoplasma contamination after cell isolation and culture. For experiments, mycoplasma-free fibroblasts were used between the third and ninth passages.

In vitro cell culture

Primary human lung fibroblasts were grown in growth media [DMEM low glucose (Gibco) with 0.025 M HEPES (Gibco) and 0.375% sodium hydrogen carbonate (Ajax Finechem, Thermo Fisher Scientific)] incubated in a humidified CO₂ incubator (5% CO₂ in air) at 37°C. Growth media were aspirated from cells every 4 days and replaced with fresh growth media until confluence.

Protein extraction from primary human lung fibroblasts

Whole-cell lysates were prepared using protein extraction buffer [20 mM Tris (Amresco, Dallas, TX, USA) pH 7.4 with 150 mM sodium chloride (PanReac Applichem, Darmstadt, Germany), 1 mM EGTA (Sigma, St Louis, MO, USA), 1 mM EDTA (Amresco), 1 mM sodium fluoride (NaF; Sigma), 20 mM sodium pyrophosphate (Sigma), 2 mM sodium orthovanadate (Na₃VO₄; Sigma), 1% (v/v) Triton X-100 (Sigma), 10% (v/v) glycerol (Sigma), 0.1% (w/v) SDS (Amresco), 0.5% (w/v) sodium deoxycholate (Sigma), 1 mM PMSF (Life Technologies, Solon, OH, USA), 1% (v/v) Protease Inhibitor Cocktail Set III (Calbiochem, San Diego, CA, USA)] on ice, scraped into a tube, and centrifuged at 150×g for 5 min. The supernatants were collected and stored at −80°C for further experiments. Total sample protein was measured using a standard BCA assay (Sigma).

SHG microscopy

SHG imaging and analysis were performed on paraffin-embedded, formalin-fixed lung tissue sections as previously published (Tjin et al., 2014).

Thirty-micrometre-thick samples were used for SHG analysis in human tissue samples (cohort 1, non-IPF=7 versus IPF=8; and cohort 2, non-IPF ILD=8 versus IPF=26) and reseeded decellularized human lung tissue samples. Five-micrometre-thick tissue samples were used for the SHG in collagen gel samples, as the softness of the sample prevented the use of thicker samples because of the unavoidable compression of these samples during processing. Confocal images were taken at ×63 magnification at 16-bit resolution for both forward and backward propagated SHG signals, with identical detector configurations for all images taken. Laser power was measured at the objective before the start of each experiment performed on different days, and laser power settings were adjusted so that all experiments use the same laser power (25 mW).

Owing to the heterogeneous nature of the lung tissues, a 10×10 grid was applied on the human tissue samples using the tile function of the Leica LAS AF software (Leica, Wetzlar, Germany). Three regions for imaging were then selected using a random number generator. A 5×5 grid was used on decellularized tissue samples because of the smaller size of the samples compared with the size of human lung tissue samples. Collagen gel samples were more homogeneous, and three regions were randomly selected for SHG.

In the figures, adjustments of brightness and contrast as indicated in the figure legends were applied uniformly for every set of images and were used for display purposes only.

SHG image analysis

Image analysis was performed using Fiji (Schindelin et al., 2012) with images imported from the LAS AF software (Leica). The stack was deliberately oversampled by 20% in each direction to ensure that the entirety of the tissue section was imaged. The confocal 3D image was analysed using the histogram to determine the mean (µ) and s.d. (σ) of the signal for the image. A threshold of µ+σ was then applied to the image to exclude background emissions as previously published (Abraham et al., 2012).

The pixel area (the number of pixels with intensity above threshold) and pixel density (average signal intensity per pixel with intensity above threshold) was measured for every slice in the image using batch processing. The total signal intensity (total intensity for all pixels with intensity above threshold) was determined as the product of the pixel area and density. The average values for these measurements were calculated for every image in the stack. The above process was repeated for both forward and backward propagated data to generate area, intensity and density data for both forward and backward propagated signals. The ratio of forward to backward signal was then calculated for area, intensity and density measurements.

Lysyl oxidase gene expression in lung tissues

Two publically available microarray data sets from lung tissue samples of IPF patients and healthy controls were analysed (GSE48149, non-IPF=9 versus IPF=13; and GSE24206, non-IPF=6 versus IPF=10). Microarray analysis was conducted using R software version 3.02, using the Bioconductor-limma package, and normalized using Robust Multi-array Average (RMA). A Benjamini–Hochberg procedure (false discovery rate) was undertaken to adjust for multiple testing.

Lysyl oxidase protein expression in lung tissues

Immunohistochemistry was performed using a protocol adapted from Faiz et al. (2013) with antibodies against LOX (Abcam, Cambridge, UK; ab31238; www.abcam.com/lox-antibody-ab31238.html, accessed 2 May 2016), LOXL1 (Novus Biologicals, Littleton, CO, USA; NBP1-82827; www.novusbio.com/LOXL1-Antibody_NBP1-82827.html, accessed 2 May 2016) and LOXL2 (Novus Biologicals; NBP1-32954; www.novusbio.com/Lysyl-Oxidase-Homolog-2-LOXL2-Antibody_NBP1-32954.html, accessed 2 May 2016). Primary antibodies were as follows: LOX (1 µg ml⁻¹), LOXL1 (0.1 µg ml⁻¹) or LOXL2 (1 µg ml⁻¹) and isotype control antibody (same concentration as the respective primary antibodies). Following washing, samples were treated with EnVision+ anti-rabbit horseradish peroxidase (HRP)-conjugated secondary antibody (K4003, Dako, Santa Clara, CA, USA). After a second wash step, colour development was detected with 3′,3′-diaminobenzidine (DAB; Dako), and the tissues were then counterstained in ready-to-use Mayer's Haematoxylin solution (Sigma) for 5 min. Sections were then dehydrated, mounted and coverslipped. All sections for each primary antigen were stained in the same staining run.

IHC image analysis

Images were captured using a Wide-field FL and TL microscope ZEISS Axio Scan.Z1 Slide Scanner (Zeiss, Oberkochen, Germany). The immunostained tissue was scanned at ×20 magnification, and a whole brightfield tissue image was stitched together from the individual serial images using ZEN Software (Zeiss). All images for each set of immunostaining were captured on the same day with the same microscope settings to ensure uniformity and accuracy of quantified image analysis. Fiji ImageJ (Schindelin et al., 2012) was used to quantify the density and distribution of staining. Colour deconvolution (Ruifrok and Johnston, 2001) vectors in ImageJ were optimized to ensure accurate separation of Eosin and DAB. Overlaid thresholded images of Eosin and DAB were used to measure total tissue surface area. The image analysis calculated the number of pixels above the threshold within the image (‘area’) and the average intensity of the pixels above the threshold (‘average intensity’). Macros were used to batch process the images.

Data were represented as ‘percentage tissue surface area’ (DAB area/total tissue surface area) and ‘density’ (average intensity calculated for the whole tissue).

Lysyl oxidase protein expression in cell lysates

Western immunoblot was performed using antibodies against LOX (Abcam; ab31238; 1:2000), LOXL1 (Novus Biologicals; NBP1-82827; 1:2000) and LOXL2 (Novus Biologicals; NBP1-32954; 1:2000) using a protocol adapted from Burgess et al. (2003). Cell lysates were run on 10% polyacrylamide gels and transferred to polyvinyledene difluoride membranes (Merck Millipore, Darmstadt, Germany) and protein blocked for 1 h in 5% skimmed milk. After washing, samples were treated with HRP-conjugated secondary antibody (Dako). After a second wash step, detection was performed using chemiluminsecent reagent Immobilon Western (Merck Millipore) and bands were analysed using a Kodak imaging system and Carestream MI software (Molecular Bioimaging, Bend, OR, USA). Membranes were stripped and reprobed for GAPDH (Merck Millipore; 1:5000 dilution) as above as a loading control.

Adjustments of brightness and contrast as indicated in the figure legends were applied uniformly for every set of images and used for display purposes only.

Antifibrotic potential of LO activity inhibitors (collagen gel 3D model)

Collagen gel solution was prepared for a final collagen gel concentration of 2 mg ml⁻¹. Briefly, 10% (v/v) of 10× PBS and 10% (v/v) 0.1 M NaOH were added to stock rat-tail collagen I (BD Biosciences, Franklin Lakes, NJ, USA) gel solution and mixed while cooled on ice. DMEM media [DMEM low glucose (Gibco) with 0.025 M HEPES (Gibco) and 0.375% sodium hydrogen carbonate (Ajax Finechem)] containing resuspended cells was then added at higher cell concentration for dilution in the collagen gel solution (final concentration of 2.5×10⁵ cells ml⁻¹ in the collagen gels). Collagen gel mixed with cells (n=5 non-IPF, n=5 IPF cell lines and n=1 no-cell control) was pipetted into a 48-well plate and allowed to set in a humidified CO₂ incubator at 37°C for 2 h. Once the gel had set, the gels were then released from the wells by gently running a sterile syringe needle around between the gel and the wall of the wells. Treatment media [Set 1: 0.1% bovine serum albumin (BSA) ±10 ng ml⁻¹ TGF-β ±0.1 mM β-APN; Set 2: 0.01% DMSO +0.1% BSA ±10 ng ml⁻¹ TGF-β ±300 nM Compound A (Pharmaxis, Sydney, Australia)] of the same volume as the gel were pipetted gently on top of the gels in the wells. The collagen gels were collected at day 7, fixed in 4% paraformaldehyde solution and paraffin embedded.

Decellularization of human lung tissues

Decellularized tissues from a non-diseased subject and an IPF subject were prepared as previously described (Booth et al., 2012). Essentially, lung tissues were sequentially incubated in sterile water, deoxycholic acid (Sigma), sodium chloride (PanReac) and DNAse (Sigma) to clear cellular and nuclear materials from the tissues. Tissues were washed three times in sterile PBS between each step. The resultant acellular tissue was biopsied using a 12-mm punch and then sliced into 1-mm slices on a vibratome. The acellular tissues slices were sterilized in peracetic acid (Sigma)/ethanol and rinsed extensively in sterile PBS. Slices were stored in sterile PBS at 4°C until use.

Antifibrotic potential of LO activity inhibitors (decellularized matrix 3D model)

Decellularized lung slices from both non-diseased and IPF subjects were washed three times in sterile DMEM solution and were subsequently added to a suspension of 5.0×10⁵ fibroblasts ml⁻¹ in growth media at a ratio of 1 matrix:2 ml of cell solution (n=5 non-IPF, n=5 IPF cell lines and n=1 no-cell control) up to a maximum of 25 ml of solution in a 50 ml Falcon tube. The Falcon tubes containing the matrices and cell suspension were mounted on a suspension mixer (Selby, ThermoFisher Scientific, Waltham, MA, USA) and constantly rotated at low speed (approximately one rotation every 20 s) at 37°C for 24 h. The reseeded ECMs were placed in separate wells of a 48-well plate with fresh media [0.01% DMSO +0.1% BSA ±10 ng ml⁻¹ TGF-β (R&D Systems, MN, USA) ±300 nM Compound A or 0.1 mM β-APN (Sigma)] and incubated at 37°C and 5% CO₂. The reseeded decellularized matrices were collected at day 7; half were fixed in 4% paraformaldehyde solution for 24 h and paraffin embedded, whereas the other half were fresh frozen in optimal cutting temperature (OCT) compound and stored at −80°C.

Tissue stiffness measurements using atomic force microscopy

AFM measurements were performed on decellularized human lung matrices with or without reseeded non-IPF cells. Matrices were embedded in optimal OCT compound, and 10-µm-thick tissue slices were cryosectioned at −21°C and mounted on poly-L-lysine-coated glass slides. PBS solution was added on the tissue slice to avoid drying.

Parenchymal tissue areas were identified under an optical microscope (×200 magnification; Olympus, Tokyo, Japan) and mechanically analysed with a BioScope Catalyst AFM (Bruker, Billerica, MA, USA). Microindentations were performed using a 2.5-µm-radius sphere-tipped probe (Novascan, Ames, IA, USA) with a spring constant determined at ∼100 pN nm⁻¹ by the thermal fluctuation method (Thundat et al., 1994). For each sample, three randomly selected areas from three non-consecutive tissue slices were analysed in PBS at room temperature. Force curves were acquired with MIRO 2.0 (NanoScope 9.1; Bruker) at an indentation rate of 20 μm s⁻¹ and a ramp size of 10 μm at different points along the parenchymal tissue. More than 100 force curves were performed per sample (∼35 per area; Fig. S8).

The Young's modulus, E, was determined by fitting of the force curve by the Hertz sphere model (Dimitriadis et al., 2002; Hertz, 1881) using NanoScope Analysis software (Bruker) and considering Poisson's ratio of 0.4 (Butler et al., 1986).

Statistical analysis

Data were entered and sorted using Microsoft Excel (Microsoft, Redmond, WA, USA). Sample sizes for comparing IPF versus non-IPF were determined based on a previous study (Tjin et al., 2014). Statistical analysis was performed using GraphPad Prism Version 6 Software (GraphPad, La Jolla, CA, USA). Data were tested for normal distribution with the D'Agostino and Pearson omnibus normality test and analysed with Student's t-test (normally distributed) or the Mann–Whitney U-test (non-normally distributed) as appropriate. Kruskal–Wallis and two-way ANOVA statistical analyses were performed where appropriate. Factor analysis was performed using SPSS software (IBM, Armonk, NY, USA) using ‘principal axis factoring’ analysing ‘correlation matrix’ and extracting ‘2’ factors from the data set. Promax rotation was used for the factor analysis, using a regression method and excluding cases listwise where there were missing data points.

Acknowledgements

The authors would like to acknowledge the subjects of the study, the cardiopulmonary transplant team and pathologists at St Vincent's Hospital and the team at the Department of Tissue Pathology and Diagnostic Oncology, Royal Prince Alfred Hospital. The authors acknowledge the facilities and the scientific and technical assistance of the Australian Microscopy & Microanalysis Research Facility at the Australian Centre for Microscopy & Microanalysis at the University of Sydney and the University of Sydney Histopathology laboratory. Dr Wolfgang Jarolimek and Dr Heidi Schilter from Pharmaxis kindly provided Compound A. Mr Mario D'Souza, a research statistician employed by Sydney Local Health District Clinical Research Centre, provided assistance with the statistical analysis.

Footnotes

Author contributions metadata

Conceptualization: G.T., E.P.W.K., J.K.B.; Methodology: G.T., E.P.W.K.; Software: G.T.; Validation: G.T., A.F., D.S., D.J.T.; Formal analysis: G.T., E.S.W., A.F., D.S., D.J.T., A.M., E.P.W.K., J.K.B.; Investigation: G.T., A.F., D.S.; Resources: G.T., E.S.W., D.J.T., A.M., E.P.W.K., J.K.B.; Data curation: G.T., J.K.B.; Writing - original draft: G.T.; Writing - review & editing: G.T., E.S.W., A.F., D.S., D.J.T., A.M., E.P.W.K., J.K.B.; Visualization: G.T., A.F., D.S.; Supervision: E.P.W.K., J.K.B.; Project administration: G.T., E.S.W., D.J.T., E.P.W.K., J.K.B.; Funding acquisition: J.K.B.

Funding

J.K.B. was supported by a National Health and Medical Research Council Australia Fellowship (1032695) and a Rosalind Franklin Fellowship, co-funded by Rijksuniversiteit Groningen, Universitair Medisch Centrum Groningen and the European Union. A.F. was funded by a RESPIRE2 fellowship from the European Respiratory Society. D.J.T., D.S. and E.S.W. were funded by the National Institutes of Health (HL092961, HL133320 and U01HL111016; UL1TR00043 and UL1TR002240 to E.S.W.).

Data availability

Two publically available microarray data sets from lung tissue samples of IPF and healthy controls, GSE48149 and GSE24206, were used in this study.

Competing interests

The authors declare no competing or financial interests.

References

Abraham , T. , Kayra , D. , McManus , B. and Scott , A. ( 2012 ). Quantitative assessment of forward and backward second harmonic three dimensional images of collagen Type I matrix remodeling in a stimulated cellular environment . J. Struct. Biol. 180 , 17 - 25 . https://doi.org/10.1016/j.jsb.2012.05.004

Aumiller , V. , Strobel , B. , Romeike , M. , Schuler , M. , Stierstorfer , B. E. and Kreuz , S. ( 2017 ). Comparative analysis of lysyl oxidase (like) family members in pulmonary fibrosis . Sci. Rep. 7 , 149 . https://doi.org/10.1038/s41598-017-00270-0

Balestrini , J. L. , Chaudhry , S. , Sarrazy , V. , Koehler , A. and Hinz , B. ( 2012 ). The mechanical memory of lung myofibroblasts . Integr. Biol. 4 , 410 - 421 . https://doi.org/10.1039/c2ib00149g

Barry-Hamilton , V. , Spangler , R. , Marshall , D. , McCauley , S. , Rodriguez , H. M. , Oyasu , M. , Mikels , A. , Vaysberg , M. , Ghermazien , H. , Wai , C. , et al. ( 2010 ). Allosteric inhibition of lysyl oxidase-like-2 impedes the development of a pathologic microenvironment . Nat. Med. 16 , 1009 - 1017 . https://doi.org/10.1038/nm.2208

Booth , A. J. , Hadley , R. , Cornett , A. M. , Dreffs , A. A. , Matthes , S. A. , Tsui , J. L. , Weiss , K. , Horowitz , J. C. , Fiore , V. F. , Barker , T. H. , et al. ( 2012 ). Acellular normal and fibrotic human lung matrices as a culture system for in vitro investigation . Am. J. Respir. Crit. Care. Med. 186 , 866 - 876 . https://doi.org/10.1164/rccm.201204-0754OC

Broekelmann , T. J. , Limper , A. H. , Colby , T. V. and McDonald , J. A. ( 1991 ). Transforming growth factor beta 1 is present at sites of extracellular matrix gene expression in human pulmonary fibrosis . Proc. Natl. Acad. Sci. USA 88 , 6642 - 6646 . https://doi.org/10.1073/pnas.88.15.6642

Brown , E. , McKee , T. , diTomaso , E. , Pluen , A. , Seed , B. , Boucher , Y. and Jain , R. K. ( 2003 ). Dynamic imaging of collagen and its modulation in tumors in vivo using second-harmonic generation . Nat. Med. 9 , 796 - 801 . https://doi.org/10.1038/nm879

Burgess , J. K. , Johnson , P. R. A. , Ge , Q. , Au , W. W. , Poniris , M. H. , McParland , B. E. , King , G. , Roth , M. and Black , J. L. ( 2003 ). Expression of connective tissue growth factor in asthmatic airway smooth muscle cells . Am. J. Respir. Crit. Care. Med. 167 , 71 - 77 . https://doi.org/10.1164/rccm.200205-416OC

Burgess , J. K. , Mauad , T. , Tjin , G. , Karlsson , J. C. and Westergren-Thorsson , G. ( 2016 ). The extracellular matrix - the under-recognized element in lung disease? J. Pathol. 240 , 397 - 409 . https://doi.org/10.1002/path.4808

Butler , J. P. , Nakamura , M. , Sasaki , H. , Sasaki , T. and Takishima , T. ( 1986 ). Poissons’ ratio of lung parenchyma and parenchymal interaction with bronchi . Jpn. J. Physiol. 36 , 91 - 106 . https://doi.org/10.2170/jjphysiol.36.91

Caetano-Lopes , J. , Nery , A. M. , Canhão , H. , Duarte , J. , Cascão , R. , Rodrigues , A. , Perpétuo , I. P. , Abdulghani , S. , Amaral , P. M. , Sakaguchi , S. , et al. ( 2010 ). Chronic arthritis leads to disturbances in the bone collagen network . Arthritis Res. Ther. 12 , R9 . https://doi.org/10.1186/ar2908

Chu , S.-W. , Tai , S.-P. , Chan , M.-C. , Sun , C.-K. , Hsiao , I.-C. , Lin , C.-H. , Chen , Y.-C. and Lin , B.-L. ( 2007 ). Thickness dependence of optical second harmonic generation in collagen fibrils . Opt. Express 15 , 12005 - 12010 . https://doi.org/10.1364/OE.15.012005

Clarke , D. L. , Carruthers , A. M. , Mustelin , T. and Murray , L. A. ( 2013 ). Matrix regulation of idiopathic pulmonary fibrosis: the role of enzymes . Fibrogenesis Tissue Repair 6 , 20 . https://doi.org/10.1186/1755-1536-6-20

Cox , G. , Kable , E. , Jones , A. , Fraser , I. , Manconi , F. and Gorrell , M. D. ( 2003 ). 3-dimensional imaging of collagen using second harmonic generation . J. Struct. Biol. 141 , 53 - 62 . https://doi.org/10.1016/S1047-8477(02)00576-2

Demedts , M. , Wells , A. U. , Anto , J. M. , Costabel , U. , Hubbard , R. , Cullinan , P. , Slabbynck , H. , Rizzato , G. , Poletti , V. , Verbeken , E. K. , et al. ( 2001 ). Interstitial lung diseases: an epidemiological overview . Eur. Respir. J. Suppl. 32 , 2s - 16s .

Dimitriadis , E. K. , Horkay , F. , Maresca , J. , Kachar , B. and Chadwick , R. S. ( 2002 ). Determination of elastic moduli of thin layers of soft material using the atomic force microscope . Biophys. J. 82 , 2798 - 2810 . https://doi.org/10.1016/S0006-3495(02)75620-8

Erler , J. T. , Bennewith , K. L. , Nicolau , M. , Dornhöfer , N. , Kong , C. , Le , Q.-T. , Chi , J.-T. A. , Jeffrey , S. S. and Giaccia , A. J. ( 2006 ). Lysyl oxidase is essential for hypoxia-induced metastasis . Nature 440 , 1222 - 1226 . https://doi.org/10.1038/nature04695

Estany , S. , Vicens-Zygmunt , V. , Llatjós , R. , Montes , A. , Penín , R. , Escobar , I. , Xaubet , A. , Santos , S. , Manresa , F. , Dorca , J. , et al. ( 2014 ). Lung fibrotic tenascin-C upregulation is associated with other extracellular matrix proteins and induced by TGFbeta1 . BMC Pulm. Med. 14 , 120 . https://doi.org/10.1186/1471-2466-14-120

Faiz , A. , Tjin , G. , Harkness , L. , Weckmann , M. , Bao , S. , Black , J. L. , Oliver , B. G. G. and Burgess , J. K. ( 2013 ). The expression and activity of cathepsins D, H and K in asthmatic airways . PLoS ONE 8 , e57245 . https://doi.org/10.1371/journal.pone.0057245

Giménez , A. , Uriarte , J. J. , Vieyra , J. , Navajas , D. and Alcaraz , J. ( 2017 ). Elastic properties of hydrogels and decellularized tissue sections used in mechanobiology studies probed by atomic force microscopy . Microsc. Res. Tech. 80 , 85 - 96 . https://doi.org/10.1002/jemt.22740

Guo , Y. , Savage , H. E. , Liu , F. , Schantz , S. P. , Ho , P. P. and Alfano , R. R. ( 1999 ). Subsurface tumor progression investigated by noninvasive optical second harmonic tomography . Proc. Natl. Acad. Sci. USA 96 , 10854 - 10856 . https://doi.org/10.1073/pnas.96.19.10854

Hertz , H. ( 1881 ). Üeber die Berührung fester elastischer Körper . J. Reine Angew. Math. 92 , 156 - 171 .

Hollosi , P. , Yakushiji , J. K. , Fong , K. S. K. , Csiszar , K. and Fong , S. F. T. ( 2009 ). Lysyl oxidase-like 2 promotes migration in noninvasive breast cancer cells but not in normal breast epithelial cells . Int. J. Cancer 125 , 318 - 327 . https://doi.org/10.1002/ijc.24308

King , T. E. , Jr, Bradford , W. Z. , Castro-Bernardini , S. , Fagan , E. A. , Glaspole , I. , Glassberg , M. K. , Gorina , E. , Hopkins , P. M. , Kardatzke , D. , Lancaster , L. , et al. ( 2014 ). A phase 3 trial of pirfenidone in patients with idiopathic pulmonary fibrosis . N. Engl. J. Med. 370 , 2083 - 2092 . https://doi.org/10.1056/NEJMoa1402582

Kottmann , R. M. , Sharp , J. , Owens , K. , Salzman , P. , Xiao , G.-Q. , Phipps , R. P. , Sime , P. J. , Brown , E. B. and Perry , S. W. ( 2015 ). Second harmonic generation microscopy reveals altered collagen microstructure in usual interstitial pneumonia versus healthy lung . Respir. Res. 16 , 61 . https://doi.org/10.1186/s12931-015-0220-8

Krimmer , D. I. , Burgess , J. K. , Wooi , T. K. , Black , J. L. and Oliver , B. G. G. ( 2012 ). Matrix proteins from smoke-exposed fibroblasts are pro-proliferative . Am. J. Respir. Cell Mol. Biol. 46 , 34 - 39 . https://doi.org/10.1165/rcmb.2010-0426OC

Kuhn , C. and McDonald , J. A. ( 1991 ). The roles of the myofibroblast in idiopathic pulmonary fibrosis. Ultrastructural and immunohistochemical features of sites of active extracellular matrix synthesis . Am. J. Pathol. 138 , 1257 - 1265 .

Kuhn , C. , III, Boldt , J. , King , T. E. , Jr, Crouch , E. , Vartio , T. and McDonald , J. A. ( 1989 ). An immunohistochemical study of architectural remodeling and connective tissue synthesis in pulmonary fibrosis . Am. Rev. Respir. Dis. 140 , 1693 - 1703 . https://doi.org/10.1164/ajrccm/140.6.1693

Liu , F. , Mih , J. D. , Shea , B. S. , Kho , A. T. , Sharif , A. S. , Tager , A. M. and Tschumperlin , D. J. ( 2010 ). Feedback amplification of fibrosis through matrix stiffening and COX-2 suppression . J. Cell Biol. 190 , 693 - 706 . https://doi.org/10.1083/jcb.201004082

Marinkovic , A. , Mih , J. D. , Park , J. A. , Liu , F. and Tschumperlin , D. J. ( 2012 ). Improved throughput traction microscopy reveals pivotal role for matrix stiffness in fibroblast contractility and TGF-beta responsiveness . Am. J. Physiol. Lung Cell. Mol. Physiol. 303 , L169 - L180 . https://doi.org/10.1152/ajplung.00108.2012

Marinkovic , A. , Liu , F. and Tschumperlin , D. J. ( 2013 ). Matrices of physiologic stiffness potently inactivate idiopathic pulmonary fibrosis fibroblasts . Am. J. Respir. Cell Mol. Biol. 48 , 422 - 430 . https://doi.org/10.1165/rcmb.2012-0335OC

Mazzei , M. E. , Richeldi , L. and Collard , H. R. ( 2015 ). Nintedanib in the treatment of idiopathic pulmonary fibrosis . Ther Adv Respir Dis . 9 , 121 - 129 . https://doi.org/10.1177/1753465815579365

Melo , E. , Garreta , E. , Luque , T. , Cortiella , J. , Nichols , J. , Navajas , D. and Farré , R. ( 2014 ). Effects of the decellularization method on the local stiffness of acellular lungs . Tissue Eng. Part C Methods 20 , 412 - 422 . https://doi.org/10.1089/ten.tec.2013.0325

Moodley , Y. , Corte , T. , Richeldi , L. and King , T. E. , Jr. ( 2015 ). Do all patients with idiopathic pulmonary fibrosis warrant a trial of therapeutic intervention? A pro-con perspective . Respirology 20 , 389 - 394 . https://doi.org/10.1111/resp.12484

Nadiarnykh , O. , Plotnikov , S. , Mohler , W. A. , Kalajzic , I. , Redford-Badwal , D. and Campagnola , P. J. ( 2007 ). Second harmonic generation imaging microscopy studies of osteogenesis imperfecta . J. Biomed. Opt. 12 , 051805 . https://doi.org/10.1117/1.2799538

Parker , M. W. , Rossi , D. , Peterson , M. , Smith , K. , Sikström , K. , White , E. S. , Connett , J. E. , Henke , C. A. , Larsson , O. and Bitterman , P. B. ( 2014 ). Fibrotic extracellular matrix activates a profibrotic positive feedback loop . J. Clin. Invest. 124 , 1622 - 1635 . https://doi.org/10.1172/JCI71386

Paszek , M. J. , Zahir , N. , Johnson , K. R. , Lakins , J. N. , Rozenberg , G. I. , Gefen , A. , Reinhart-King , C. A. , Margulies , S. S. , Dembo , M. , Boettiger , D. , et al. ( 2005 ). Tensional homeostasis and the malignant phenotype . Cancer Cell 8 , 241 - 254 . https://doi.org/10.1016/j.ccr.2005.08.010

Payne , S. L. , Fogelgren , B. , Hess , A. R. , Seftor , E. A. , Wiley , E. L. , Fong , S. F. T. , Csiszar , K. , Hendrix , M. J. C. and Kirschmann , D. A. ( 2005 ). Lysyl oxidase regulates breast cancer cell migration and adhesion through a hydrogen peroxide-mediated mechanism . Cancer Res. 65 , 11429 - 11436 . https://doi.org/10.1158/0008-5472.CAN-05-1274

Raghu , G. , Collard , H. R. , Egan , J. J. , Martinez , F. J. , Behr , J. , Brown , K. K. , Colby , T. V. , Cordier , J.-F. , Flaherty , K. R. , Lasky , J. A. , et al. ( 2011 ). 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 . https://doi.org/10.1164/rccm.2009-040GL

Raghu , G. , Brown , K. K. , Collard , H. R. , Cottin , V. , Gibson , K. F. , Kaner , R. J. , Lederer , D. J. , Martinez , F. J. , Noble , P. W. , Song , J. W. , et al. ( 2017 ). Efficacy of simtuzumab versus placebo in patients with idiopathic pulmonary fibrosis: a randomised, double-blind, controlled, phase 2 trial . Lancet Respir. Med. 5 , 22 - 32 . https://doi.org/10.1016/S2213-2600(16)30421-0

Richeldi , L. , Cottin , V. , Flaherty , K. R. , Kolb , M. , Inoue , Y. , Raghu , G. , Taniguchi , H. , Hansell , D. M. , Nicholson , A. G. , Le Maulf , F. , et al. ( 2014a ). Design of the INPULSIS™ trials: two phase 3 trials of nintedanib in patients with idiopathic pulmonary fibrosis . Respir. Med. 108 , 1023 - 1030 . https://doi.org/10.1016/j.rmed.2014.04.011

Richeldi , L. , du Bois , R. M. , Raghu , G. , Azuma , A. , Brown , K. K. , Costabel , U. , Cottin , V. , Flaherty , K. R. , Hansell , D. M. , Inoue , Y. , et al. ( 2014b ). Efficacy and safety of nintedanib in idiopathic pulmonary fibrosis . N. Engl. J. Med. 370 , 2071 - 2082 . https://doi.org/10.1056/NEJMoa1402584

Ruifrok , A. C. and Johnston , D. A. ( 2001 ). Quantification of histochemical staining by color deconvolution . Anal. Quant. Cytol. Histol. 23 , 291 - 299 .

Schilter , H. , Buson , A. , Cock , T.-A. , Deodhar , M. , Findlay , A. D. , Foot , J. S. , Jarnicki , A. , Moses , J. , Turner , C. I. , Yow , T. T. , et al. ( 2014 ). Inhibition of LOXL2: Pharmaxis’ small molecule approach to treat fibrosis . In C30. Late Breaking Abstracts in Disease Treatment and Clinical Outcomes , pp. A6669-A6669 , American Thoracic Society 2014 International Conference. New York: ATS Journals .

Schindelin , J. , Arganda-Carreras , I. , Frise , E. , Kaynig , V. , Longair , M. , Pietzsch , T. , Preibisch , S. , Rueden , C. , Saalfeld , S. , Schmid , B. , et al. ( 2012 ). Fiji: an open-source platform for biological-image analysis . Nat. Methods 9 , 676 - 682 . https://doi.org/10.1038/nmeth.2019

Shi , M. , Zhu , J. , Wang , R. , Chen , X. , Mi , L. , Walz , T. and Springer , T. A. ( 2011 ). Latent TGF-beta structure and activation . Nature 474 , 343 - 349 . https://doi.org/10.1038/nature10152

Sivaguru , M. , Durgam , S. , Ambekar , R. , Luedtke , D. , Fried , G. , Stewart , A. and Toussaint , K. C. , Jr. ( 2010 ). Quantitative analysis of collagen fiber organization in injured tendons using Fourier transform-second harmonic generation imaging . Opt. Express 18 , 24983 - 24993 . https://doi.org/10.1364/OE.18.024983

Suzuki , M. , Kayra , D. , Elliott , W. M. , Hogg , J. C. and Abraham , T. ( 2012 ). Second harmonic generation microscopy differentiates collagen type I and type III in diseased lung tissues. Proc. SPIE 8226, Multiphoton Microscopy in the Biomedical Sciences X11 , 82263F . https://doi.org/10.1117/12.910815

Thundat , T. , Warmack , R. J. , Chen , G. Y. and Allison , D. P. ( 1994 ). Thermal and ambient-induced deflections of scanning force microscope cantilevers . Appl. Phys. Lett. 64 , 2894 - 2896 . https://doi.org/10.1063/1.111407

Tjin , G. , Xu , P. , Kable , S. H. , Kable , E. P. W. and Burgess , J. K. ( 2014 ). Quantification of collagen I in airway tissues using second harmonic generation . J. Biomed. Opt. 19 , 036005 . https://doi.org/10.1117/1.JBO.19.3.036005

Vadasz , Z. , Kessler , O. , Akiri , G. , Gengrinovitch , S. , Kagan , H. M. , Baruch , Y. , Izhak , O. B. and Neufeld , G. ( 2005 ). Abnormal deposition of collagen around hepatocytes in Wilson's disease is associated with hepatocyte specific expression of lysyl oxidase and lysyl oxidase like protein-2 . J. Hepatol. 43 , 499 - 507 . https://doi.org/10.1016/j.jhep.2005.02.052

Valeyre , D. , Albera , C. , Bradford , W. Z. , Costabel , U. , King , T. E. , Jr, Leff , J. A. , Noble , P. W. , Sahn , S. A. and du Bois , R. M. ( 2014 ). Comprehensive assessment of the long-term safety of pirfenidone in patients with idiopathic pulmonary fibrosis . Respirology 19 , 740 - 747 . https://doi.org/10.1111/resp.12297

Ward , P. A. and Hunninghake , G. W. ( 1998 ). Lung inflammation and fibrosis . Am. J. Respir. Crit. Care. Med. 157 , S123 - S129 . https://doi.org/10.1164/ajrccm.157.4.nhlbi-10

Westergren-Thorsson , G. , Sime , P. , Jordana , M. , Gauldie , J. , Särnstrand , B. and Malmström , A. ( 2004 ). Lung fibroblast clones from normal and fibrotic subjects differ in hyaluronan and decorin production and rate of proliferation . Int. J. Biochem. Cell Biol. 36 , 1573 - 1584 . https://doi.org/10.1016/j.biocel.2004.01.009

White , E. S. , Thannickal , V. J. , Carskadon , S. L. , Dickie , E. G. , Livant , D. L. , Markwart , S. , Toews , G. B. and Arenberg , D. A. ( 2003 ). Integrin alpha4beta1 regulates migration across basement membranes by lung fibroblasts: a role for phosphatase and tensin homologue deleted on chromosome 10 . Am. J. Respir. Crit. Care. Med. 168 , 436 - 442 . https://doi.org/10.1164/rccm.200301-041OC

Williams , R. M. , Zipfel , W. R. and Webb , W. W. ( 2001 ). Multiphoton microscopy in biological research . Curr. Opin. Chem. Biol. 5 , 603 - 608 . https://doi.org/10.1016/S1367-5931(00)00241-6

Williams , R. M. , Zipfel , W. R. and Webb , W. W. ( 2005 ). Interpreting second-harmonic generation images of collagen I fibrils . Biophys. J. 88 , 1377 - 1386 . https://doi.org/10.1529/biophysj.104.047308

Wipff , P. J. , Rifkin , D. B. , Meister , J. J. and Hinz , B. ( 2007 ). Myofibroblast contraction activates latent TGF-beta1 from the extracellular matrix . J. Cell Biol. 179 , 1311 - 1323 . https://doi.org/10.1083/jcb.200704042