The stiffness of the extracellular matrix strongly influences multiple behaviors of adherent cells. Matrix stiffness varies spatially throughout a tissue, and undergoes modification in various disease conditions. Here we develop methods to characterize spatial variations in stiffness in normal and fibrotic mouse lung tissue using atomic force microscopy microindentation.
Matrix stiffness strongly influences growth, differentiation and function of adherent cells1-3. On the macro scale the stiffness of tissues and organs within the human body span several orders of magnitude4. Much less is known about how stiffness varies spatially within tissues, and what the scope and spatial scale of stiffness changes are in disease processes that result in tissue remodeling. To better understand how changes in matrix stiffness contribute to cellular physiology in health and disease, measurements of tissue stiffness obtained at a spatial scale relevant to resident cells are needed. This is particularly true for the lung, a highly compliant and elastic tissue in which matrix remodeling is a prominent feature in diseases such as asthma, emphysema, hypertension and fibrosis. To characterize the local mechanical environment of lung parenchyma at a spatial scale relevant to resident cells, we have developed methods to directly measure the local elastic properties of fresh murine lung tissue using atomic force microscopy (AFM) microindentation. With appropriate choice of AFM indentor, cantilever, and indentation depth, these methods allow measurements of local tissue shear modulus in parallel with phase contrast and fluorescence imaging of the region of interest. Systematic sampling of tissue strips provides maps of tissue mechanical properties that reveal local spatial variations in shear modulus. Correlations between mechanical properties and underlying anatomical and pathological features illustrate how stiffness varies with matrix deposition in fibrosis. These methods can be extended to other soft tissues and disease processes to reveal how local tissue mechanical properties vary across space and disease progression.
1. Lung Tissue Strip Preparation
2. AFM Microindentation and Fluorescence Imaging
We recommend using a silicon nitride triangle cantilever with a 5-μm diameter borosilicate spherical tip (Novascan). Using AFM probes with a spring constant of 0.06 N/m, we have mechanically characterized soft materials with shear moduli spanning 100 to 50,000 Pa .
The indentation velocity should be selected to be sufficiently slow to explore elastic rather viscoelastic properties of the soft sample7 . A velocity range of 2-20 μm per second is suitable for lung tissue.
3. AFM Data Extraction
Figure 1 (A) Schematic of AFM microindentation of a soft sample on a rigid substrate using a spherical probe (Reproduced with permission from Dimitriadis E., et al. 2002 Biophys J8). (B) Representative force-displacement curve collected from a clean glass slide in PBS to determine cantilever deflection sensitivity. (C) Representative force-displacement curves collected from soft (blue) and stiff (red) samples showing the corresponding parameters from (A).
4. Representative Results:
Figure 2A shows a lung parenchyma strip attached on poly-l-lysine coated coverslip in PBS with an AFM probe in place directly above region of interest and ready for microindentation mapping. Figure 2B shows that in properly cut and stained tissue, the alveolar microarchitecture of lung parenchyma is well preserved, as observed by immunofluorescent staining for basement membrane component laminin without fixation or permeabilization. Figure 2 C-D show immunofluorescent staining for collagen I in fresh unfixed lungs harvested from mice previously treated with PBS (Fig. 2C), or treated with bleomycin (Fig. 2D) to induce fibrosis10.
Figure 3A shows sample force-displacement plots obtained from AFM microindentation. With the same applied force, AFM tip generates a large indentation on a soft region (blue line), resulting in a relatively “flat” force-displacement curve versus a small indentation and a “steep” force-displacement curve for a stiffer region (red line). In order to obtain clean force curves, after each indentation the AFM tip needs to be retracted completely from the sample surface and free of contact before the next indentation. This contact free state corresponds to the flat region of the curve in Figure 3A, where the tip translates without cantilever deflection. Figure 3B shows a typical false curve without a flat region which occurs when the tip is not completely retracted from the sample surface (e.g. soft sample attaches to tip) making it impossible to determine the contact point. If the tip is completely trapped in soft sample, the curve may look as shown in Figure 3C, without clear deflection but small noise.
Figure 4 shows shear modulus data extracted from force-displacement curves displayed spatially as stiffness maps (elastographs), where color scales with shear modulus. Elastographs demonstrate striking differences in the range and distribution of tissue shear modulus in normal (Fig. 4A) and fibrotic (Fig. 4B) lung parenchyma, and large spatial variations in shear modulus, particularly within the fibrotic lung sample.
Figure 2. (A) Phase contrast micrograph showing an AFM probe over a mouse lung parenchyma strip. Scale bar, 50 μm. (B) Immunostaining of laminin in normal mouse lung tissue showing the alveolar microarchitecture of normal lung parenchyma. White box indicates typical 80-by-80μm elasticity mapping area. Scale bar: 20 μm. (C and D) Immunostaining of collagen I in saline (C) and bleomycin (D) treated mouse lung parenchyma. Scale bars, 100 μm.
Figure 3. (A) Representative interpretable force-displacement curves collected from saline (blue) and bleomycin-treated (red) mouse lung parenchyma, respectively. Black lines are the best fit resulting from the spherical Hertz model and are labeled with their corresponding calculated Young’s moduli. (B, C) Representative un-interpretable force-displacement curves.
Figure 4. Representative elastographs collected from saline (A) and bleomycin-treated (B) mouse lung parenchyma. The color bars indicate shear modulus in kilopascals. Axis labels indicate spatial scale in micrometers. Scale bar: 20 μm
Mechanical characterization of lung tissue using AFM microindentation offers unprecedented spatial resolution (Fig. 4), providing a unique perspective on microscale variations in tissue stiffness. As an example of its utility, previous macro-scale measurements in normal and fibrotic lung tissue strips indicated an approximate 2-3-fold increase in elastance with fibrosis11,12. In contrast, AFM microindentation reveals that tissue stiffening is highly localized, with some regions exhibiting up to ~30-fold increases in shear modulus above the median observed in normal lung tissue10. As matrix stiffness is now known to critically influence cell function, these local measurements provide invaluable parameters to enhance the biofidelity of cell culture study of lung cells.
Several practical issues arise with the use of thin strips of lung tissue. The surfaces of the strips are not perfectly flat, as the tissue profile follows the architecture of the underlying alveoli. The AFM system automatically adjusts tip position in the Z-direction during indentation when the sample surface height variation is smaller than 15-μm to help overcome this challenge. Measurements are made at room temperature, not 37°C, so deviations in mechanical properties caused by this variation in temperature cannot be evaluated, though they would be expected to be minor. The influence of the underlying alveolar wall architecture on observed mechanical properties is difficult to determine with the current light microscopy setup. For instance, it would be desirable to determine if alveolar walls exhibit anisotropy and different mechanical properties when indented in directions aligned with or transverse to the plane of the wall. However, the current specimens are thick and the imaging system does not have 3D capabilities, hence it is not possible to determine the local alveolar wall orientation at each point of contact. Finally, the influence of cellular constituents on measured mechanical properties remains to be fully elucidated. In the methods detailed here, no efforts are made to specifically remove cellular constituents of the tissue. However, the cells present at the surface available for indentation are unlikely to be viable given the time elapsed since tissue harvest and the necessity of cutting the tissue to gain access to the alveoli. Specific experiments to remove cells, or repopulate matrices with viable cells, and evaluate the resulting changes in tissue stiffness would appear warranted.
Because fresh unfixed tissue is needed for these measurements, the time elapsed from tissue harvest to measurement should be minimized and samples should be stored at 4 °C to avoid changes in mechanical properties. Particular attention should be paid whenever tissue strips are transferred between containers during washing or staining so that minimum distortion or damage is generated. For AFM application in liquid, a crucial step is to cut the tissue as flat as possible and immobilize the sample on the supporting coverslip. If available, an automated sectioning machine such as a vibratome or tissue slicer can be used to cut slices of highly uniform thickness. It’s important to attach tissue strips immediately before AFM measurements and minimize the time elapsed for AFM measurements as the sample will eventually detach from the coverslips. One useful observation is that larger strips appear to attach more stably to cover slips and remain in place for longer durations in PBS than smaller strips.
AFM microindentation can characterize samples spanning a broad range from 100 Pa to 50 kPa (shear modulus) when using a standard 0.06 N/m cantilever with a 5 μm diameter spherical tip. This range can be expanded using probes with different spring constants; AFM probes with spherical glass tips ranging from 0.6 to 12 μm in diameter and spring constants ranging from 0.01 to 0.58 N/m are commercially available (e.g. Novascan) and commonly used3. With a 5 μm spherical tip, the theoretical contact area between the tip and tissue is about 5-9 μm2 for 400-700 nm indentation (Fig. 1A). Smaller or larger tips can be used to provide smaller or larger scales of spatial resolution. Pyramidal tips have also been used in AFM microindentation13-16, providing smaller contact areas and thus increasing spatial resolution in mapping, though data fitting is more complex for this tip geometry.
Several limitations to this method should be noted. The lung has traditionally been mechanically characterized non-invasively, for instance using pressure-volume analysis17 or punch-indentation of whole lungs19,20. Invasive methods such as the one described here alter the lung architecture in important ways through the loss of the air-liquid interface that normally exists in the air-filled lung and the loss of pre-stress that maintains lung partial inflation upon relaxation of respiratory muscles. These limitations are common to all measurements made in lung tissue strips18. Notably, however, the median stiffness measured in the parenchyma of normal lung tissue (shear modulus ~0.5kPa) does not differ substantially from estimates based on punch-indentation of intact lungs at resting volumes19,20. While lung tissue is known to exhibit non-linear stiffening with increasing deformation, it is not possible to test in a rigorous fashion whether this property persists down to the micro-scale with the methods employed here. The Hertz model assumes homogeneity of the sample. However, most biological materials, including lung parenchyma, are increasingly heterogeneous at decreasing spatial scales. Heterogeneity of the sample can result in artifacts like variation of the Young’s modulus depending on indentation depth, i.e. depending on the layer or component that is deforming. The heterogeneity in the xy-plane can be limited by carefully choosing the appropriate spherical tip size depending on the microstructure of the biomaterial as proposed by Dimitriadis EK et al.8 It is much more difficult to predict or correct the Hertz model error due to material heterogeneity in the z-direction. Azeloglu et al. recently proposed a hybrid computational model to characterize the elastic properties of heterogeneous substrate with discrete embedded inclusions21. Their new technique provides a potential means to calculate inclusion properties of heterogeneous materials overcoming the limitations of Hertzian analysis.
The Hertz model also assumes absolute elastic behavior, while biological materials typically display time-dependent viscoelastic behaviors. A full viscoelastic characterization of tissue can be obtained by varying the indentation velocities used. Importantly, previous macro-scale mechanical testing of normal and fibrotic lung tissue demonstrates weak frequency dependence of lung tissue mechanical properties, and preservation of the differences between normal and fibrotic tissue mechanical properties across all frequencies tested11. These findings strongly suggest that the measurement of mechanical properties using a single indentation velocity with AFM captures an essential aspect of the changes in tissue mechanical properties that accompany fibrosis.
The Poisson’s ratio of 0.4 for lung tissue used in this work is from a macroscopic measurement9. Unfortunately, the Poisson’s ratio at micro-scale and any changes under disease condition are not available in the literature. As alternatives to E, E/(1-υ2) or (1-υ2)/πE (denoted the elastic constant k)22 can be calculated from AFM microindentation and reported when the Poisson’s ratio is unknown. For most biomaterials Poisson’s ratio is in the range of 0.4 to 0.5 due to their high water content. Within the range 0.3-0.5, the factor 1/(1-υ2) varies only from 1.10-1.33, such that reasonable variations in the Poisson’s ratio exert only modest effects on the reported modulus. The increase in shear modulus that we report for fibrotic tissue relative to normal tissue is several fold in magnitude, implying that errors associated with variations in Poisson’s ratio are minor relative to the changes in mechanical properties observed.
The actual algorithm and code that can be used for the analysis of force-displacement data is subject to the specific application condition and the subsequent characteristics of the diverse populations of force-displacement curves. If more sophisticated analysis is of interest, one may consult the work of Lin et al.23. The authors compiled a series of synergistic strategies into an algorithm that overcomes many of the complications that have previously impeded efforts to automate the fitting of Hertz models of indentation data.
Several other areas are available for further development and exploitation of this method. In cases where one is interested in visualizing the alveolar walls without antibody labeling, both elastin and collagen can be visualized from their autofluorescent signal in the green spectrum. On the other hand, better imaging, using either thinner tissue sections, 3D imaging techniques, or both, could enhance the ability to correlate tissue architecture with underlying mechanical properties. While the current methods allow staining and visualization of extracellular matrix components such as collagen and laminin, additional efforts could be aimed at staining cell surface markers to identify specific cell populations and to characterize the mechanical microenvironments in the vicinity of such populations. Alternatively, tissue could be harvested from mice expressing fluorescently-tagged lineage markers or cell specific proteins to pursue the same goal. Finally, the method detailed here appears well suited to characterizing other anatomical features in the lung, such as vessels which remodel in hypertension, and airways which remodel in asthma. Based on its current state of development and potential for further enhancement, AFM microindentation appears poised to yield valuable insights into the changes in tissue stiffness that accompany disease progression in the lung, and will no doubt be of value in characterizing spatial and temporal changes in the stiffness of a variety of other soft tissues.
The authors have nothing to disclose.
We thank A. Tager and B. Shea for their ongoing collaboration, and for providing the lung tissue used for demonstration purposes here. This work was supported by National Institutes of Health grant HL-092961. This work was performed in part at the Center for Nanoscale Systems (CNS), a member of the National Nanotechnology Infrastructure Network, which is supported by the National Science Foundation under NSF award ECS-0335765. CNS is part of the Faculty of Arts and Sciences at Harvard University.
Name of the reagent | Company | Catalogue number | Comments |
---|---|---|---|
AFM tip | Novascan | PT.GS | 5 micron Borosilicate bead, 0.06 N/m |
Poly-L-Lysine coverslips | VWR | 354085 | BD BioCoat 12 mm round No.1 glass |
Agarose, Low Gelling Temperature | Sigma | A0701 | |
Rabbit anti-Collagen I (Rb pAb) antibody | Abcam | ab34710-100 | |
Alexa Fluor 546 goat anti-rabbit antibody | Invitrogen | A11010 | |
Rabbit anti-Laminin | Sigma | L9393 |