Skip to main content
Download PDF

Examining lung mechanical strains as influenced by breathing volumes and rates using experimental digital image correlation

Respiratory Research volume 23, Article number: 92 (2022) Cite this article

Abstract

Background

Mechanical ventilation is often employed to facilitate breathing in patients suffering from respiratory illnesses and disabilities. Despite the benefits, there are risks associated with ventilator-induced lung injuries and death, driving investigations for alternative ventilation techniques to improve mechanical ventilation, such as multi-oscillatory and high-frequency ventilation; however, few studies have evaluated fundamental lung mechanical local deformations under variable loading.

Methods

Porcine whole lung samples were analyzed using a novel application of digital image correlation interfaced with an electromechanical ventilation system to associate the local behavior to the global volume and pressure loading in response to various inflation volumes and breathing rates. Strains, anisotropy, tissue compliance, and the evolutionary response of the inflating lung were analyzed.

Results

Experiments demonstrated a direct and near one-to-one linear relationship between applied lung volumes and resulting local mean strain, and a nonlinear relationship between lung pressures and strains. As the applied air delivery volume was doubled, the tissue surface mean strains approximately increased from 20 to 40%, and average maximum strains measured 70–110%. The tissue strain anisotropic ratio ranged from 0.81 to 0.86 and decreased with greater inflation volumes. Local tissue compliance during the inflation cycle, associating evolutionary strains in response to inflation pressures, was also quantified.

Conclusion

Ventilation frequencies were not found to influence the local stretch response. Strain measures significantly increased and the anisotropic ratio decreased between the smallest and greatest tidal volumes. Tissue compliance did not exhibit a unifying trend. The insights provided by the real-time continuous measures, and the kinetics to kinematics pulmonary linkage established by this study offers valuable characterizations for computational models and establishes a framework for future studies to compare healthy and diseased lung mechanics to further consider alternatives for effective ventilation strategies.

Background

The global pandemic has pushed pulmonary mechanics research to the forefront. Prior to COVID-19, respiratory illnesses were already one of the leading causes of death worldwide [1, 2]. Burdensome respiratory ailments such as chronic obstructive pulmonary disease (COPD) are the leading causes of disability amongst an expansive list of lung diseases which includes tuberculosis and acute lower respiratory tract infection, all detrimental to lung function [3]. Other examples of respiratory illnesses affecting countless individuals which lead to irreversible tissue damage are pneumonia, asthma, and emphysema [4,5,6].

Severe manifestations of such respiratory diseases can often lead to the need for additional breathing support in the form of mechanical ventilation (MV). Preexisting lung conditions play a major role in inflicting tissue damage when mechanically ventilated, and the risk of ventilator induced lung injuries (VILI), such as barotrauma and edema, are greater [7]. Another possibility of MV is ventilator-associated pneumonia which affects 8–28% of patients with high chances of obtaining this complication after just 48 h on ventilation, demonstrative of the quick onset of complications due to ventilation [8, 9].

Factors that contribute to these injuries are heavily related to the ventilator driving pressure as well as the mechanical elasticity, integrity, and behavior of the lung tissue when ventilated, which is compromised in diseased states [10,11,12,13]. The behavioral mechanics of lung tissue has been explored through investigations of isolated parenchymal strips and tracheal and bronchial airway segments via tensile [14,15,16], biaxial [17, 18], and indentation tests [19,20,21]. However, there are limited studies considering the local stretch behavior of tissue in situ, where the tissue still resides on the intact lung organ. The local tissue behavior in response to globally applied variable inflation volumes and breathing rates is particularly key to understanding strain heterogeneities and concentrations [22, 23], and can support studies focusing on alternative ventilation strategies to minimize VILI (e.g. multi-oscillatory and high-frequency ventilation) [24, 25].

As such, the objective of this current study is to characterize ventilation strains. We do so by utilizing digital image correlation (DIC) as a novel technique for analyzing topological deformations of the inflating whole porcine lung organ. Building upon previously introduced methods, we interface two systems to simultaneously assess local strain and global pulmonary measurements throughout the lung expansion cycle [26, 27]. Existing strain measurement studies, such as digital volume correlation (DVC), optical coherence tomography (OCT) [28,29,30], computerized tomography (CT) scans [31, 32] and doppler elastography and speckle tracking for biological tissues [33, 34] are limited because they are at discrete time points, retrieve 2D contour images with projected strains, or are time intensive, which affects the resulting measured lung behavior.

In contrast, in this study we employ DIC as a real-time continuous imaging approach to examine lung deformations instantaneously and divulge regional pulmonary behaviors, such as tissue strain heterogeneity and isotropy, and further investigate volume and rate ventilation effects on these lung behaviors. The resulting comprehensive characterization of local–global lung mechanics in pig specimens offers critical pulmonary insights for new oscillatory ventilation schemes. Given the similarity between porcine and human lung size and anatomy, findings can help support studies aimed at lowering the risk of ventilation injuries and complications [35]. Furthermore, the variable inflation volume and breathing rate data collected in this work can inform finite element computational models to enable predictive ventilation techniques [36].

Methods

Sample preparation

Pulmonary samples were obtained from four pigs (A, B, C, and D) weighing approximately 200–250 lbs, ranging 6–8 months in age (local abattoir, IACUC approval not required). These specimens were shipped with wet ice and stored in a plastic container, enveloped in a barrier of bubble wrap to protect the tissue from direct contact with the ice. Testing was conducted immediately upon arrival and within 36 h postmortem.

Each lung specimen was weighed, and the initial lung volume was measured based on liquid submersion [37]. The lung was inflated by a laboratory air-line pressure outlet (< 1 psi) to open collapsed airways and to expand the surface for the speckling process and digital image correlation (DIC) technique previously detailed in Mariano et al. (2020) [27]. Briefly, a thin layer of quick-drying white enamel paint (rust-oleum) and an exfoliator pad was used to apply a randomized speckling pattern onto the specimen surface. The speckled lung was then placed inside our custom electromechanical pressure–volume (PV) system within an air-tight transparent tank for imaging and DIC data collection (Trilion ARAMIS Adjustable 12M system) as seen in Fig. 1 [26].

Fig. 1
figure 1

Experimental set-up of the electromechanical pressure–volume ventilation apparatus (right) interfaced with the digital image correlation system (left). The Trilion ARAMIS Adjustable 12M two-camera system is positioned above a transparent, air-tight tank containing the lung specimen which is controlled by the dual-piston apparatus to apply inflation volumes and measure resulting lung volumes and pressures. The combination of these two systems enables the simultaneous collection of global pressures and volumes and local lung topological strain measurements

Full size image

The volume of trapped air within each lung was calculated by subtracting the initial lung volume from its total volume computed by using the weight density equation, where the lung tissue density of 1.06 g/ml was used [38]. The peak lung volume and pressure were measured at the highest applied inflation volume’s (1350 ml) maximum values (Table 1).

Table 1 Specimen characteristics for varying breathing rates
Full size table

Experimental protocol and data collection

The sample was hydrated with 1XPBS solution atop a frictionless platform. Using our custom-built PV system depicted in Fig. 1 (design described in detail in Sattari et al. 2020), a preload of 0.05psi ensured common initial states between lung specimens [26]. The specimens underwent three preconditioning cycles inflated to specified volumes before the fourth cycle was analyzed [39]. The sample was inflated within the total lung capacity of porcine lung to three maximum volumes of 675, 900, and 1350 ml corresponding to tidal volumes of 6, 8, and 12 ml/kg in clinical practice and at rates of 10, 15, and 20 breaths per minute (BPM) [40, 41].

DIC data was collected at frequencies corresponding to the rate of inflation as follows: 10 Hz for 10 BPM, 15 Hz for 15 BPM, and 20 Hz for 20 BPM. The total dimensional DIC measuring volume was 375 × 295 × 295 mm. A facet size of 30 pixels and point distance of 10 pixels were used across all samples. Surface components were created at the maximum inflation stage of the lung to analyze the entire surface at full expansion as previously described [27]. The principal strains and stretches, major and minor, were calculated based on the reference stage of the uninflated, preconditioned state of the specimen for each volume and rate combination, which were then averaged across all lungs for each set of testing parameters.

Data analysis

Strain data collected from commercial GOM software (Trillion, King of Prussia, Pennsylvania, USA) was plotted as topological color contour maps to observe qualitative lung inflation patterns and visually assess the localized strains (Figs. 2, 7B). Histograms of the strain values at the maximum inflation stage were quantified for the various inflation volumes and multiple breathing rates to assess the heterogenous and isotropic tissue behaviors relative to their corresponding fraction of the lung surface (Figs. 3 and 5).

Fig. 2
figure 2

Topological strain contour maps of representative porcine lung at varying breathing rates of 10, 15, and 20 breaths per minute and inflation volumes of 675 ml, 900 ml, and 1350 ml at the maximum inflation stage. As the applied volume increased, surface strain heterogeneity became more prominent as indicated by appearance of regional high strain concentrations in the upper caudal left lung lobe and middle lobe of the right lung. The influence of breathing rate on local strains was not observed to be unidirectional

Full size image
Fig. 3
figure 3

Averaged quantified strains across the fraction of the specimen’s topological surface at the maximum inflation stage divulged the effects of varying breathing rate and inflation volumes. Surface strains increased with increasing inflation volumes. Lower inflation volumes exhibited right skewed distribution and spread to a centered distribution at higher volumes, indicative of expanding strain ranges and heterogeneity. No statistical rate dependency was found, however the strains did depend on the applied volume (see Fig. 4)

Full size image

Table 2 documents individual specimen’s mean, median, maximum, and range of major principal strains at maximum inflation. The anisotropic ratio was calculated as the fraction of minor to major stretch which resulted in a ratio range of 0–1, where unity represents isotropic stretch. Outlying strain values were marked using the interquartile range analysis [42, 43].

Table 2 Individual and averaged specimen strain data measures and anisotropic ratio
Full size table

Associations between local and global measurements of mean strain, applied volumes, and lung pressures were collected and analyzed throughout inflation. The local compliance, the relationship between local tissue strains and lung pressures was bilinear (Fig. 6). To calculate the initial (slope 1) and final slopes (slope 2) indicative of alveolar recruitment and optimal compliance respectively [44], the curves were fitted with R2 = 0.9, where R2 is the coefficient of determination [45]. Additional local to global behavioral insights were assessed by analyzing intermediate inflation stages of a representative pig at 15BPM and 1350 ml applied volume (Fig. 7A), through topological contour maps (Fig. 7B) and histograms (Fig. 7C).

Data was analyzed by a two-way repeated measurements ANOVA test (GraphPad Prism) with Bonferroni correction for strain means, medians, ranges, and maximums at various volumes and breathing rates, in addition to anisotropic ratio comparisons (Table 2, Fig. 4). Significance was set at p < 0.05. p-values less than 0.05, 0.01, and 0.001 were represented by *, **, and ***, respectively.

Fig. 4
figure 4

Statistical analysis of mean (A), median (B), range (C), and maximum strains (D), and anisotropic ratios (E) for varying rates and inflation volumes. Significant volume dependency was observed for all data sets between 675 and 1350 ml, and no significant rate effects were found

Full size image

Results

Specimens’ weights, initial lung volumes, peak lung volumes, and peak pressures for three breathing rates were displayed in Table 1. The mean weight of the lung samples was 903 g. Based on the averages of the four specimens, the amount of air trapped within the lung is 46% of the initial lung volume. Peak pressures experienced for all rates were in the range of 0.36–0.52psi and the peak lung volume inflation due to compressed air ranged 1023.8–1215.3 ml, less than the applied 1350 ml. At the largest applied inflation volume of 1350 ml, the peak pressures increased with increasing breathing rate and alternatively, the lung volume response decreased (10 to 20BPM: peak pressure p = 0.009 and peak volume p  <0.001; 15 to 20BPM: peak pressure p = 0.015 and peak volume p = 0.003).

A representative topological strain colormap was illustrated in Fig. 2. All specimens were observed to experience higher strains at higher volumes as well as greater strains predominantly located where the cranial and caudal lobe meet on the left lung and within the middle lobe of the right lung [46]. Localized high strains often manifested as stripes extending downward towards the lower end of each caudal lobe.

Figure 3 demonstrated the corresponding histograms of the major strain data at maximum inflation; there was a right skewed distribution populating ~ 20% mean strain at lower volumes of 675 ml demonstrating surface homogeneity in comparison to the expanded distribution at increased inflation volume of 1350 ml averaging ~ 40% strain.

Table 2 listed the strain means, medians, maximums, and ranges for each specimen. Statistical significance was found for each rate of inflation as the averaged mean strain values increased with greater inflation volumes (Fig. 4A). Median strain values had statistical significance between the volumes of 675 ml with 1350 ml (10BPM: p < 0.0001, 15BPM: p < 0.0001, 20BPM: p < 0.0001), and 900 ml with 1350 ml (10BPM: p = 0.002, 15BPM: p = 0.003, 20BPM: p = 0.004) for all three rates (Fig. 4B). Average range and maximum strains were also higher with increasing volumes (Fig. 4C and D) and showed statistically significant differences between the lowest and highest inflation volumes of 675 ml and 1350 ml (range—10BPM: p = 0.008, 15BPM: p = 0.029, 20BPM: p = 0.021; strain—10BPM: p < 0.001, 15BPM: p < 0.001, 20BPM: p < 0.001). Breathing rate did not affect the strain means, medians, maximums, and ranges, however those values increased with increasing applied volumes.

In contrast, the anisotropic ratio decreased for increasing applied volumes, averaging 0.86, 0.84, and 0.81 for 675 ml, 900 ml, and 1350 ml, respectively for 20BPM (Table 2). Breathing rate had a negligible effect on the anisotropic ratio. With increasing volume, the tissue exhibited greater directional expansion preference as defined by the decreasing anisotropic ratio, which was found to be statistically significant between the 675 and 1350 ml volumes (10BPM: p = 0.005, 15BPM: p = 0.007, 20BPM: p = 0.009) (Fig. 4E). The histogram of the anisotropic ratio as a fraction of the lung surface found the lung tissue is near isotropic as the averaged means of the ratios stayed within the range of 0.81–0.86, closer to the value of 1 (Fig. 5). In agreement with previous insights, the specimens primarily expanded in the lateral-medial and ventral-dorsal direction (across the ribcage) [27].

Fig. 5
figure 5

Anisotropic ratio for the fractional surface of lung specimens across varying inflation rates and volumes. An averaged mean ranging 0.81–0.86 was observed. Rate had no effect on the anisotropic ratio however increased inflation volumes decreased the anisotropic ratio (Fig. 4E)

Full size image

Associations between local and global behavior for each specimen was highlighted in Fig. 6. The trend for average slope 1 appears to increase for increasing volumes, however, is only statistically significant for 20BPM when comparing 675 ml and 900 ml to 1350 ml (675 ml: p = 0.011, 900 ml: p = 0.048). Slope 2 also trends higher with increasing volumes however is significant only for volume comparisons of 675 ml and 1350 ml for 10BPM (p = 0.005) and between the rates of 10BPM and 20BPM for 1350 ml inflation volume (p = 0.009). Breathing rates did not generally affect the slopes but decreased slope 2 trends were noted at faster rates for 900 ml and 1350 ml, (significant between 10 and 20BPM for 1350 ml). The mean strain evolution as a function of applied volume plot (not shown) was also analyzed and was found to not be dependent on the rate nor the amount of volume applied.

Fig. 6
figure 6

Associated global pressures and local mean strain. The slopes of these bilinear graphs and the calculated inflection points were analyzed for statistical significance to assess local tissue compliancy

Full size image

Figure 7A demonstrated the inflation trajectory for an applied volume of 1350 ml of a representative pig at four nominal stages, which linked lung pressures and volumes to the spatial distribution of local strains as exhibited on the contour map (Fig. 7B) and was quantified in the corresponding histograms (Fig. 7C). The evolution of the contour map showed local strains (including the mean, range, and maximums) increased with greater global inflation volume and pressure. The strains were observed to continuously evolve in real-time and resulted in regions of high strain concentration at the middle and upper caudal lobes of the right and left lung respectively. In the select stages corresponding to a global lung inflation volume increase of 322–1350 ml, the global lung pressure rises from 0.30 to 0.40psi, and the maximum strain grows from 29.5 to 79.8%. Histograms illustrated right skewed strain distributions initially with lower mean strain of 9.4% and 0–30% range which transitioned to higher mean strains of 37.4% coupled with increased strain ranges of 0–80% and greater surface heterogeneity.

Fig. 7
figure 7

Evolution of the local strain as related to global inflation pressures and volumes. Four intermediate stages of a representative pig lung at 15 BPM inflated to 1350 ml (A) was presented as qualitative strain contours (B) and quantitative histograms (C). Increasing global lung pressure and volumes corresponded to increasing local strain maximum and mean values as well as a larger range of experienced strains. The transition from a right skewed distribution at stage 1 to a more distributed spread of strain values for stages 2–4 was observed and illustrated greater surface heterogeneity

Full size image

Discussion

This work builds upon our previous methods introducing DIC as a novel technique for analyzing topological strains of the whole porcine lung organ and associates quantifications of global pulmonary loading and local strain measures to further characterize porcine surface heterogeneity, anisotropy, and the effect of volume and breathing rate on lung deformations for the first time. Strains as high as 111% are reported with an average of 40% at the greatest inflation tidal volume of 1350 ml. We observe isotropic-leaning organ distension tendencies, and significant local strain dependencies on the global pressures and volumes imposed, especially at the greatest volumes and breathing rates; on the other hand, we observe frequency dependent trends only for a single comparison of local compliancy. Interestingly, we find breathing rates play no significant role in local lung strain values. Findings are specifically translatable to clinical research investigating the role of pulmonary frequency dependencies, such as multi-oscillatory ventilation and inspiratory rate induced lung injury [24, 47].

Lung surface heterogeneity

Upper regions of the caudal lobe in the left, and middle lobe in the right lung experience higher strains as seen in Fig. 2 [46]. This difference in behavior at the upper and lower regions of the lung have been previously seen in the literature where inflation tests on dog and mice lung found similar regional behaviors of nonuniform air distribution; greater distension is seen in the upper airways of mice where the large bronchi of the left and right lung reside [48, 49]. This regional behavior may be attributable to the bronchi network, which is heterogenous in geometry and material properties [50]; passageways narrow in diameter distally and therefore, air expansion when approaching the distal regions [51]. Similarly, in humans the alveolar size is four times greater in the upper regions compared to the lower regions and can more readily conduct the higher observed lung stretches [52].

Lung surface heterogeneity can also be qualitatively observed in the strain distributions and strain topological maps at various applied volumes and rates (Figs. 2 and 3). The high striped areas of strain can be due to the presence of underlying airways, which alter the deformation behavior of the enveloping parenchymal lung tissue [53]. Post-experimental dissection demonstrated that the left caudal lobe and right middle and upper caudal lobes were heavily populated by large underlying airways corresponding to the striped areas. The greater regional surface strains observed in our experiments may well be caused by this underlying bronchial network, especially as these strain contour hotspots extended distally across the surface as the values increased.

Volume effects

Our study finds statistically significant evidence of volume dependence, where lung surface strains as well as pressure increases with increasing inflation volume. In a study employing CT scans where porcine lungs were inflated using our same low and high tidal volumes, this increase in local strain in response to increasing applied lung volume was also observed within the parenchyma [54]. Volume-strain interdependence was also seen in cat [55] and mice [56] specimens. This real-time ex-vivo DIC technique additionally offers insights regarding the temporal evolution of strains as the volume is incrementally applied; as such, we observe the spatial expansion of the lower lungs as more air is delivered (Fig. 2). Compared to previous studies finding the highest strain regions in pig specimens to be central-dorsally located [54], we find the greatest concentration of strains to be dorsal but to be more upper than central; this may be attributed to the ex-vivo nature of our experiments compared to lung specimens with an intact chest cavity.

Rate effects

Our experiments observe that the breathing rate did not significantly affect the local strains (Fig. 4). In comparable DIC experiments conducted on mice specimens, similar rate insignificance was concluded [56]. We do however see significance in rate variations where peak pressure increases and peak lung volume decreases with greater respiratory frequency between the two smaller rates of 10 and 15 BPM compared to 20 BPM (Table 1). This may be attributable to the viscoelastic nature of lung tissue [57], which causes pressure losses at slow breathing rates [55, 58]. The faster rates may not allow the tissue to adequately react and the air to fully permeate throughout the airways, causing higher pressures and lower volumes in the inflated lung. The limited lung studies exploring this phenomenon on isolated tissue strips in tension or compression have found greater stiffness at higher rates [19, 59], which is a similar trend for our whole organ inflation tests: the final slope tends to stiffen with faster breathing rates (not significant, Fig. 6). Despite our study concluding insignificant rate effects, the additional exploration of breathing frequencies is merited as inspiratory rate is one of the main causes of concern in VILI [47].

Strain anisotropy

Lung surface strain is found to be anisotropic with a ratio ranging 0.81–0.86. Previous studies in the literature consisting of parenchymal tissue samples [60, 61] have concluded lung isotropy. This can partly be due to the presence of isotropically expanding alveoli [62], however, our observed anisotropy, in agreement with our previous conclusions finding reduced expansion of the lung in the cranial-caudal direction, may be attributed to the underlying airway network previously demonstrated to be twice as stiff along the bronchial length compared to circumferentially [15, 27]. While past experiments inflating whole dog lungs filled with saline [63] have concluded isotropy, our air-filled lungs are known to exhibit alternate mechanics not comparable to liquid inflation [64, 65]. Historical studies on air-filled orthogonally tracked markers on rabbit lungs found near equal directional deformations but were technologically limiting and may not capture pulmonary anisotropy at the resolution presented here [66]. A comparison of both air-filled and saline-filled lung tissue observed higher degrees of anisotropy in air, in agreement with our experiments [53].

Global to local

Utilizing our system, we are able to simultaneously analyze the evolutionary relationship between the local and global lung response of the lungs. Associated mean strains and pressures in Fig. 6 shows how the tissue exhibits a bilinear behavior, similar to existing lung studies on whole mice lungs [27, 56]. This bilinearity is thought to be an effect of alveoli opening and recruitment at low pressures and expansion at high pressures [62], and could also be attributed to the bilinear response observed in pulmonary constituents (i.e. airways, parenchyma, visceral pleura) [14, 67,68,69,70].

Our study facilitates insights regarding the spatial distribution of strains on the lung surface at varying pressure values throughout inflation stages (Fig. 7B), which cannot be observed with classical pressure–volume plots. From our histograms (Fig. 7C), we see that low pressures and volumes exhibit predominantly low strains but as global loading increases, a shift occurs where the lung inflates more heterogeneously and strains extend to a wider range of simultaneously existing low and high values. Such insights can support investigations of oscillatory ventilation schemes to better understand concentrated regions of strain and acute respiratory distress syndrome [25].

Ventilator efficacy and pulmonary disease are also dependent on tissue compliancy [25, 71, 72]; high compliance is found in lungs experiencing obstructive diseases such as emphysema and reduced compliance is found in restrictive lung diseases such as fibrosis [73,74,75]. Previous studies analyzing these global organ compliances from pressure–volume loops have observed lower lung compliance as higher lung volumes were reached due to incremental lung volume gains at higher pressures [74]; however, our local compliance measures linking mean surface strains to global pressures finds that the regional tissue compliance increases at higher inflation volumes as well as at slower breathing rates. This discrepancy is centered on the definition of local tissue versus global organ compliance; the former has not been widely measured but ex-vivo DIC real-time continuous evolutionary measurements between local strains and global loading makes this now possible.

Limitations

Despite this study being the first of its kind to systematically characterize global and local lung tissue behavior under varied applied volumes and breathing rates in porcine specimens, there are practical limitations relating to replicating physiological loading conditions as seen in past works [76]. First, the lung is explanted and removed from the restraining ribcage. While we anticipate that absolute deformation values may differ from intact thoracic cavity analyses, the relative local strain alterations due to sequential increases in applied loading or faster breathing frequencies is unlikely to be affected. Nonetheless, the ex vivo nature of these experiments must be considered when considering physiological translation, particularly as alveoli surfactants degrade ex vivo and may impact airway recruitment as seen in the inflection points of Fig. 6 [77].

Another limitation in using DIC is that the measurements are topological. Our study lacks the ability to assess internal lung behavior, such as those afforded by digital volume correlation (DVC) or CT scans; however, DVC and CT scans are time intensive and limit real-time observations of continuous strain evolution during a respiratory cycle where qualitative local strain assessments can be regionally observed as instantaneously interfaced with our global pressure volume data, as is the advantage of DIC techniques. Additionally, to maintain the focus of this study, we have limited our analyses to the inflation limb of the respiratory cycle; we note that the deflation portion of lung mechanics may exhibit alternate mechanics and intend to consider such behavior in future works.

Conclusion

The present study relates real-time topological strain data with global lung organ measurements and finds the local tissue response exhibits volume and pressure loading dependence but rate independence. The first time-continuous measures and associations between the tissue and organ level kinetics and kinematics of the porcine lung provided here inaugurates new insights which can be used in future computational models and subsequent studies seeking to contrast healthy and diseased pulmonary mechanics to improve ventilation techniques.

Availability of data and materials

The data from this study are available upon reasonable request.

Abbreviations

COPD:

Chronic obstructive pulmonary disease

MV:

Mechanical ventilation

VILI:

Ventilator-induced lung injuries

DIC:

Digital image correlation

DVC:

Digital volume correlation

OCT:

Optical coherence tomography

CT:

Computerized tomography

PV:

Pressure–volume

BPM:

Breaths per minute

References

  1. Guarascio AJ, Ray SM, Finch CK, Self TH. The clinical and economic burden of chronic obstructive pulmonary disease in the USA. Clin Outcomes Res. 2013;5(1):235–45.

    Google Scholar 

  2. Sattari S, Eskandari M. Characterizing the viscoelasticity of extra- and intra-parenchymal lung bronchi. J Mech Behav Biomed Mater. 2020;110(103824):1–10. https://doi.org/10.1016/j.jmbbm.2020.103824.

    CAS  Article  Google Scholar 

  3. Lopez AD, Murray CCJL. The global burden of disease, 1990–2020. Am J Respir Crit Care Med. 1998;4(11):1241–3.

    CAS  Google Scholar 

  4. Braun J, Dalhoff K, Schaaf B, Wood WG, Wiessmann KJ. Characterization of protein-antiproteinase imbalance in bronchoalveolar lavage from patients with pneumonia. Eur Respir J. 1994;7(1):127–33.

    CAS  PubMed  Article  Google Scholar 

  5. Murdoch JR, Lloyd CM. Chronic inflammation and asthma. Mutat Res - Fundam Mol Mech Mutagen. 2010;690(1–2):24–39. https://doi.org/10.1016/j.mrfmmm.2009.09.005.

    CAS  Article  Google Scholar 

  6. Eskandari M, Pfaller MR, Kuhl E. On the role of mechanics in chronic lung disease. Materials (Basel). 2013;6(12):5639–58.

    Article  Google Scholar 

  7. Gupta VK, Alkandari BM, Mohammed W, Tobar AM, Abdelmohsen MA. Ventilator associated lung injury in severe COVID-19 pneumonia patients – Case Reports: Ventilator associated lung injury in COVID-19. Eur J Radiol Open. 2021;8(100310):1–6. https://doi.org/10.1016/j.ejro.2020.100310.

    Article  Google Scholar 

  8. Chastre J, Fagon J. State of the art ventilator-associated pneumonia. Am J Respir Crit Care Med. 2002;165(23):867–903.

    PubMed  Article  Google Scholar 

  9. Morehead RS, Pinto SJ. Ventilator-associated pneumonia. Arch Internal Med. 2000;160(13):1926–36.

    CAS  Article  Google Scholar 

  10. Das A, Camporota L, Hardman JG, Bates DG. What links ventilator driving pressure with survival in the acute respiratory distress syndrome? A computational study. Respir Res. 2019;20(29):1–10.

    CAS  Google Scholar 

  11. Parker JC, Hernandez LA, Peevy KJ. Mechanisms of ventilator induced lung injury. Crit Care Med. 1993;21(1):131–43.

    CAS  PubMed  Article  Google Scholar 

  12. Briscoe WA, Cournand A. Uneven ventilation of normal and diseased lungs studied by an open-circuit method ’. J Appl Physiol. 1959;14(23):284–90.

    CAS  PubMed  Article  Google Scholar 

  13. Nieman GF, Satalin J, Andrews P, Aiash H, Habashi NM, Gatto LA. Personalizing mechanical ventilation according to physiologic parameters to stabilize alveoli and minimize ventilator induced lung injury (VILI). Intensive Care Med Exp. 2017;5(1):1–21.

    Article  Google Scholar 

  14. Hoppin FG, Lee GC, Dawson SV. Properties of lung parenchyma in distortion. J Appl Physiol. 1975;39(5):742–51.

    PubMed  Article  Google Scholar 

  15. Codd SL, Lambert RK, Alley MR, Pack RJ. Tensile stiffness of ovine tracheal wall. J Appl Physiol. 1994;76(6):2627–35.

    CAS  PubMed  Article  Google Scholar 

  16. Roberts CR, Rains JK, Paré PD, Walker DC, Wiggs B, Bert JL. Ultrastructure and tensile properties of human tracheal cartilage. J Biomech. 1997;31(1):81–6.

    Article  Google Scholar 

  17. Dassow C, Wiechert L, Martin C, Schumann S, Müller-Newen G, Pack O, et al. Biaxial distension of precision-cut lung slices. J Appl Physiol. 2010;108(3):713–21.

    CAS  PubMed  Article  Google Scholar 

  18. Debes JC, Fung YC. Effect of temperature on the biaxial mechanics of excised lung parenchyma of the dog. J Appl Physiol. 1992;73(3):1171–80.

    CAS  PubMed  Article  Google Scholar 

  19. Andrikakou P, Vickraman K, Arora H. On the behaviour of lung tissue under tension and compression. Sci Rep. 2016;6(36642):1–10.

    Google Scholar 

  20. Hajji M, Wilson T, Lai-Fook J. Improved measurements of shear modulus and pleural membrane tension of the lung. J Appl Physiol. 1979;47(1):175–81.

    CAS  PubMed  Article  Google Scholar 

  21. Polio SR, Kundu AN, Dougan CE, Birch NP, Ezra Aurian-Blajeni D, Schiffman JD, et al. Cross-platform mechanical characterization of lung tissue. PLoS ONE. 2018;13(10):1–17.

    Article  CAS  Google Scholar 

  22. Gattinoni L, Carlesso E, Cadringher P, Valenza F, Vagginelli F, Chiumello D. Physical and biological triggers of ventilator-induced lung injury and its prevention. Eur Respir J. 2003;22(47):15–25.

    Article  Google Scholar 

  23. Protti A, Andreis DT, Monti M, Santini A, Sparacino CC, Langer T, et al. Lung stress and strain during mechanical ventilation: any difference between statics and dynamics? Crit Care Med. 2013;41(4):1046–55.

    PubMed  Article  Google Scholar 

  24. Kaczka DW, Herrmann J, Zonneveld CE, Tingay DG, Lavizzari A, Noble PB, et al. Multifrequency oscillatory ventilation in the premature lung: effects on gas exchange, mechanics, and ventilation distribution. Anesthesiology. 2015;123(6):1394–403.

    CAS  PubMed  Article  Google Scholar 

  25. Meyers M, Rodrigues N, Ari A. High-frequency oscillatory ventilation: a narrative review. Can J Respir Ther. 2019;55(1):40–6.

    PubMed  PubMed Central  Article  Google Scholar 

  26. Sattari S, Mariano CA, Vittalbabu S, Velazquez JV, Postma J, Horst C, et al. Introducing a custom-designed volume-pressure machine for novel measurements of whole lung organ viscoelasticity and direct comparisons between positive- and negative-pressure ventilation. Front Bioeng Biotechnol. 2020;8(578762):1–12.

    Google Scholar 

  27. Mariano CA, Sattari S, Maghsoudi-Ganjeh M, Tartibi M, Lo DD, Eskandari M. Novel mechanical strain characterization of ventilated ex vivo porcine and murine lung using digital image correlation. Front Physiol. 2020;11(600492):1–12.

    Google Scholar 

  28. Arora H, Mitchell RL, Johnston R, Manolesos M, Howells D, Sherwood JM, et al. Correlating local volumetric tissue strains with global lung mechanics measurements. Materials (Basel). 2021;14(2):1–17.

    Article  CAS  Google Scholar 

  29. Acosta Santamaría VA, Flechas García M, Molimard J, Avril S. Three-dimensional full-field strain measurements across a whole porcine aorta subjected to tensile loading using optical coherence tomography-digital volume correlation. Front Mech Eng. 2018;4(3):1–14.

    Google Scholar 

  30. Germaneau A, Doumalin P, Dupré JC. 3D strain field measurement by correlation of volume images using scattered light: Recording of images and choice of marks. Strain. 2007;43(3):207–18.

    Article  Google Scholar 

  31. Hurtado DE, Erranz B, Lillo F, Sarabia-Vallejos M, Iturrieta P, Morales F, et al. Progression of regional lung strain and heterogeneity in lung injury: assessing the evolution under spontaneous breathing and mechanical ventilation. Ann Intensive Care. 2020;10(107):1–10. https://doi.org/10.1186/s13613-020-00725-0.

    CAS  Article  Google Scholar 

  32. Fredman AJ, Frolik JL, Garra BS. Lung strain profiles using computed tomography elastography. Annu Int Conf IEEE Eng Med Biol - Proc. 2004;26:1545–8.

    Google Scholar 

  33. Wang RK, Ma Z, Kirkpatrick SJ. Tissue Doppler optical coherence elastography for real time strain rate and strain mapping of soft tissue. Appl Phys Lett. 2006;89(14):3–5.

    Google Scholar 

  34. Pavlopoulos H, Nihoyannopoulos P. Strain and strain rate deformation parameters: from tissue Doppler to 2D speckle tracking. Int J Cardiovasc Imaging. 2008;24(5):479–91.

    PubMed  Article  Google Scholar 

  35. Judge EP, Hughes JML, Egan JJ, Maguire M, Molloy EL, O’Dea S. Anatomy and bronchoscopy of the porcine lung: a model for translational respiratory medicine. Am J Respir Cell Mol Biol. 2014;51(3):334–43.

    CAS  PubMed  Article  Google Scholar 

  36. Maghsoudi-Ganjeh M, Mariano CA, Sattari S, Arora H, Eskandari M. Developing a lung model in the age of COVID-19: a digital image correlation and inverse finite element analysis framework. Front Bioeng Biotechnol. 2021;9(October):1–14.

    Google Scholar 

  37. Rush SA, Maddox T, Fisk AT, Woodrey MS, Cooper RJ. A precise water displacement method for estimating egg volume. J F Ornithol. 2009;80(2):193–7.

    Article  Google Scholar 

  38. Okazawa M, D’Yachkova Y, Paré PD. Mechanical properties of lung parenchyma during bronchoconstriction. J Appl Physiol. 1999;86(2):496–502.

    CAS  PubMed  Article  Google Scholar 

  39. Eskandari M, Sattari S, Mariano CA. Investigating the mechanics of positive- versus negative-pressure ventilation. 2021;A4671–A4671.

  40. Protti A, Andreis DT, Milesi M, Iapichino GE, Monti M, Comini B, et al. Lung anatomy, energy load, and ventilator-induced lung injury. Intensive Care Med Exp. 2015;3(34):1–18. https://doi.org/10.1186/s40635-015-0070-1.

    Article  Google Scholar 

  41. Chiumello D, Carlesso E, Cadringher P, Caironi P, Valenza F, Polli F, et al. Lung stress and strain during mechanical ventilation for acute respiratory distress syndrome. Am J Respir Crit Care Med. 2008;178(4):346–55.

    PubMed  Article  Google Scholar 

  42. Amelon R, Cao K, Ding K, Christensen GE, Reinhardt JM, Raghavan ML. Three-dimensional characterization of regional lung deformation. J Biomech. 2011;44(13):2489–95. https://doi.org/10.1016/j.jbiomech.2011.06.009.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Eskandari M, Arvayo AL, Levenston ME. Mechanical properties of the airway tree: heterogeneous and anisotropic pseudoelastic and viscoelastic tissue responses. J Appl Physiol. 2018;125:878–88.

    CAS  PubMed  Article  Google Scholar 

  44. Hickling KG. The pressure-volume curve is greatly modified by recruitment. A mathematical model of ARDS lungs. Pneumologie. 1999;53(1):57–8.

    Google Scholar 

  45. Svantesson C, Drefeldt B, Sigurdsson S, Larsson A, Brochard L, Jonson B. A single computer-controlled mechanical insufflation allows determination of the pressure–volume relationship of the respiratory system. J Clin Monit Comput. 1999;15:9–16.

    CAS  PubMed  Article  Google Scholar 

  46. Chevaleyre C, Riou M, Bréa D, Vandebrouck C, Barc C, Pezant J, et al. The pig: a relevant model for evaluating the neutrophil serine protease activities during acute Pseudomonas aeruginosa lung infection. PLoS ONE. 2016;11(12):1–17.

    Article  CAS  Google Scholar 

  47. Akoumianaki E, Vaporidi K, Georgopoulos D. The injurious effects of elevated or nonelevated respiratory rate during mechanical ventilation. Am J Respir Crit Care Med. 2019;199(2):149–57.

    CAS  PubMed  Article  Google Scholar 

  48. Faridy EE, Kidd R. Topological distribution of inpsired gas in excised lobes of dogs. J Appl Physiol. 1967;22(4):760–6.

    CAS  PubMed  Article  Google Scholar 

  49. Nickles HT, Sumkauskaite M, Wang X, Wegner I, Puderbach M, Kuebler WM. Mechanical ventilation causes airway distension with proinflammatory sequelae in mice. Am J Physiol - Lung Cell Mol Physiol. 2014;307(1):27–37.

    Article  CAS  Google Scholar 

  50. Eskandari M, Kuschner WG, Kuhl E. Patient-specific airway wall remodeling in chronic lung disease. Ann Biomed Eng. 2015;43(10):2538–51.

    PubMed  PubMed Central  Article  Google Scholar 

  51. Noble PB, McLaughlin RA, West AR, Becker S, Armstrong JJ, McFawn PK, et al. Distribution of airway narrowing responses across generations and at branching points, assessed in vitro by anatomical optical coherence tomography. Respir Res. 2010;11(9):1–12.

    Google Scholar 

  52. Casha A, Alexander M, Marilyn G, Andre N, Emanuel F. Common pathophysiological pathways for apical and upper lobe lung disease. J Infect Dis Epidemiol. 2018;4(3):1–9.

    Google Scholar 

  53. Tai RC, Lee GC. Isotropy and homogeneity of lung tissue deformation. J Biomech. 1981;14(4):243–52.

    CAS  PubMed  Article  Google Scholar 

  54. Paula LF, Wellman TJ, Winkler T, Spieth PM, Güldner A, Venegas JG, et al. Regional tidal lung strain in mechanically ventilated normal lungs. J Appl Physiol. 2016;121(6):1335–47.

    PubMed  PubMed Central  Article  Google Scholar 

  55. Kochi T, Okubo S, Zin WA, Milic-Emili J. Flow and volume dependence of pulmonary mechanics in anesthetized cats. J Appl Physiol. 1988;64(1):441–50.

    CAS  PubMed  Article  Google Scholar 

  56. Nelson T, Quiros K, Mariano CA, Sattari S, Ulu A, Dominguez E, Nordgren T, Eskandari M. Associating Local Strains to Global Pressure-Volume Mouse Lung Mechanics using Digital Image Correlation. Manuscript submitted for publication. 2022.

  57. Bayliss LE, Robertson GW. The visco-elastic properties of the lungs. Exp Physiol. 1938;29(1):27–47.

    Article  Google Scholar 

  58. Saibene F, Mead J. Frequency pulmonary dependence of hysteresis. J Appl Physiol. 1969;25(6):732–7.

    Article  Google Scholar 

  59. Birzle AM, Wall WA. A viscoelastic nonlinear compressible material model of lung parenchyma—experiments and numerical identification. J Mech Behav Biomed Mater. 2019;94:164–75. https://doi.org/10.1016/j.jmbbm.2019.02.024.

    CAS  Article  PubMed  Google Scholar 

  60. Lai-Fook SJ. Lung tissue mechanics. Respir Biomech. 1990;89:39–40.

    Article  Google Scholar 

  61. Weed B, Patnaik S, Rougeau-Browning M, Brazile B, Liao J, Prabhu R, et al. Experimental evidence of mechanical isotropy in porcine lung parenchyma. Materials (Basel). 2015;8(5):2454–66.

    Article  Google Scholar 

  62. Roan E, Waters CM. What do we know about mechanical strain in lung alveoli? Am J Physiol - Lung Cell Mol Physiol. 2011;301(5):625–35.

    Article  CAS  Google Scholar 

  63. Ligas JR, Primiano FP, Saidel GM. Static mechanics of excised whole lung: theoretical framework and experimental studies. Ann Biomed Eng. 1984;12(5):421–35.

    CAS  PubMed  Article  Google Scholar 

  64. Morgan TE. Pulmonary surfactant. N Engl J Med. 1977;284(21):1185–93.

    Article  Google Scholar 

  65. Gil J, Bachofen H, Gehr P, Weibel ER. Alveolar volume-surface area relation in air- and saline-filled lungs fixed by vascular perfusion. J Appl Physiol Respir Environ Exerc Physiol. 1979;47(5):990–1001.

    CAS  PubMed  Google Scholar 

  66. Ardila R, Horie T, Hildebrandt J. Macroscopic isotropy of lung expansion. Respir Physiol. 1974;20(2):105–15.

    CAS  PubMed  Article  Google Scholar 

  67. Eskandari M, Nordgren TM, O’Connell GD. Mechanics of pulmonary airways: Linking structure to function through constitutive modeling, biochemistry, and histology. Acta Biomater. 2019;97:513–23. https://doi.org/10.1016/j.actbio.2019.07.020.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  68. McFawn PK, Mitchell HW. Bronchial compliance and wall structure during development of the immature human and pig lung. Eur Respir J. 1997;10(1):27–34.

    CAS  PubMed  Article  Google Scholar 

  69. Humphrey JD, Vawter DL, Vitoi RP. Mechanical behavior of excised canine visceral pleura. Ann Biomed Eng. 1986;14:451–66.

    CAS  PubMed  Article  Google Scholar 

  70. Vawter DL, Fung YC, West JB. Elasticity of excised dog lung parenchyma. J Appl Physiol Respir Environ Exerc Physiol. 1978;45(2):261–9.

    CAS  PubMed  Google Scholar 

  71. Irvin CG, Bates JHT. Measuring the lung function in the mouse: the challenge of size. Respir Res. 2003;9:1–9.

    Article  Google Scholar 

  72. Eskandari M, Javili A, Kuhl E. Elastosis during airway wall remodeling explains multiple co-existing instability patterns. J Theor Biol. 2016;403:209–18. https://doi.org/10.1016/j.jtbi.2016.05.022.

    Article  PubMed  PubMed Central  Google Scholar 

  73. Maghsoudi-Ganjeh M, Sattari S, Eskandari M. Mechanical behavior of the airway wall in respiratory disease. Curr Opin Physiol. 2021;22:100445. https://doi.org/10.1016/j.cophys.2021.05.008.

    Article  Google Scholar 

  74. Papandrinopoulou D, Tzouda V, Tsoukalas G. Lung compliance and chronic obstructive pulmonary disease. Pulm Med. 2012;2012(542769):1–6.

    Article  Google Scholar 

  75. Zachary Edwards PA. Physiology, Lung Compliance . StatPearls. 2020. https://www.ncbi.nlm.nih.gov/books/NBK554517/

  76. Wong WD, Wang L, Paré PD, Seow CY. Bronchodilatory effect of deep inspiration in freshly isolated sheep lungs. Am J Physiol - Lung Cell Mol Physiol. 2017;312(2):L178–85.

    PubMed  Article  Google Scholar 

  77. Beretta E, Romanò F, Sancini G, Grotberg JB, Nieman GF, Miserocchi G. Pulmonary interstitial matrix and lung fluid balance from normal to the acutely injured lung. Front Physiol. 2021;12(December):1–18.

    Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Funding support from the Dassault Systèmes U.S. Foundation Grant and the Hellman Fellows Program is gratefully acknowledged.

Author information

Affiliations

  1. Department of Mechanical Engineering, University of California at Riverside, Riverside, CA, USA

    C. A. Mariano, S. Sattari, K. A. M. Quiros, T. M. Nelson & M. Eskandari

  2. BREATHE Center, School of Medicine, University of California at Riverside, Riverside, CA, USA

    M. Eskandari

  3. Department of Bioengineering, University of California at Riverside, Riverside, CA, USA

    M. Eskandari

Authors
  1. C. A. Mariano
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. S. Sattari
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. K. A. M. Quiros
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. T. M. Nelson
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. M. Eskandari
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

ME conceptualized and supervised research. CM, ME, and SS designed the study. CM, KA, SS, TN conducted experiments. CM and ME interpreted results, analyzed data, drafted the figures, and wrote the manuscript. All authors read and approved the submitted version.

Corresponding author

Correspondence to M. Eskandari.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mariano, C.A., Sattari, S., Quiros, K.A.M. et al. Examining lung mechanical strains as influenced by breathing volumes and rates using experimental digital image correlation. Respir Res 23, 92 (2022). https://doi.org/10.1186/s12931-022-01999-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12931-022-01999-7

Keywords

  • Digital image correlation
  • Topological strain
  • Pulmonary mechanics
  • Anisotropy
  • Heterogeneity
  • Biomechanics
  • Ventilation
  • Inflation

Respiratory Research

ISSN: 1465-993X

Contact us