
==== Front
Intensive Care Med Exp
Intensive Care Med Exp
Intensive Care Medicine Experimental
2197-425X
Springer International Publishing Cham

39225817
663
10.1186/s40635-024-00663-2
Research Articles
Regional distribution of mechanical strain and macrophage-associated lung inflammation after ventilator-induced lung injury: an experimental study
Liggieri Francesco 12
Chiodaroli Elena 13
Pellegrini Mariangela 15
Puuvuori Emmi 4
Sigfridsson Jonathan 6
Velikyan Irina 4
Chiumello Davide 378
Ball Lorenzo 2
Pelosi Paolo 2
Stramaglia Sebastiano 9
Antoni Gunnar 46
Eriksson Olof 4
http://orcid.org/0000-0001-6834-6399
Perchiazzi Gaetano gaetano.perchiazzi@surgsci.uu.se

15
1 https://ror.org/048a87296 grid.8993.b 0000 0004 1936 9457 The Hedenstierna Laboratory, Department of Surgical Sciences, Uppsala University, Akademiska Sjukhuset-Ing. 40, Tr. 3, 75185 Uppsala, Sweden
2 https://ror.org/0107c5v14 grid.5606.5 0000 0001 2151 3065 Dipartimento di Scienze Diagnostiche e Chirurgiche Integrate, Università di Genova, Genoa, Italy
3 https://ror.org/03dpchx26 0000 0004 5373 4585 Department of Anesthesia and Intensive Care, ASST Santi Paolo e Carlo, San Paolo University Hospital, Milan, Italy
4 grid.8993.b 0000 0004 1936 9457 Science for Life Laboratory, Department of Medicinal Chemistry, Uppsala University, Uppsala, Sweden
5 https://ror.org/01apvbh93 grid.412354.5 0000 0001 2351 3333 Department of Anesthesia and Intensive Care Medicine, Uppsala University Hospital, Uppsala, Sweden
6 https://ror.org/01apvbh93 grid.412354.5 0000 0001 2351 3333 PET Center, Center for Medical Imaging, Uppsala University Hospital, Uppsala, Sweden
7 https://ror.org/00wjc7c48 grid.4708.b 0000 0004 1757 2822 Department of Health Sciences, University of Milan, Milan, Italy
8 https://ror.org/00wjc7c48 grid.4708.b 0000 0004 1757 2822 Coordinated Research Center on Respiratory Failure, University of Milan, Milan, Italy
9 grid.7644.1 0000 0001 0120 3326 Department of Physics, National Institute for Nuclear Physics, University of Bari Aldo Moro, Bari, Italy
3 9 2024
3 9 2024
12 2024
12 774 6 2024
23 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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/.
Background

Alveolar macrophages activation to the pro-inflammatory phenotype M1 is pivotal in the pathophysiology of Ventilator-Induced Lung Injury (VILI). Increased lung strain is a known determinant of VILI, but a direct correspondence between regional lung strain and macrophagic activation remains unestablished. [68Ga]Ga-DOTA-TATE is a Positron Emission Tomography (PET) radiopharmaceutical with a high affinity for somatostatin receptor subtype 2 (SSTR2), which is overexpressed by pro-inflammatory-activated macrophages. Aim of the study was to determine, in a porcine model of VILI, whether mechanical strain correlates topographically with distribution of activated macrophages detected by [68Ga]Ga-DOTA-TATE uptake.

Methods

Seven anesthetized pigs underwent VILI, while three served as control. Lung CT scans were acquired at incremental tidal volumes, simultaneously recording lung mechanics. [68Ga]Ga-DOTA-TATE was administered, followed by dynamic PET scans. Custom MatLab scripts generated voxel-by-voxel gas volume and strain maps from CT slices at para-diaphragmatic (Para-D) and mid-thoracic (Mid-T) levels. Analysis of regional Voxel-associated Normal Strain (VoStrain) and [68Ga]Ga-DOTA-TATE uptake was performed and a measure of the statistical correlation between these two variables was quantified using the linear mutual information (LMI) method.

Results

Compared to controls, the VILI group exhibited statistically significant higher VoStrain and Standardized Uptake Value Ratios (SUVR) both at Para-D and Mid-T levels. Both VoStrain and SUVR increased along the gravitational axis with an increment described by statistically different regression lines between VILI and healthy controls and reaching the peak in the dependent regions of the lung (for strain in VILI vs. control was at Para-D: 760 ± 210 vs. 449 ± 106; at Mid-T level 497 ± 373 vs. 193 ± 160; for SUVR, in VILI vs. control was at Para-D: 2.2 ± 1.3 vs. 1.3 ± 0.1; at Mid-T level 1.3 ± 1.0 vs. 0.6 ± 0.03). LMI in both Para-D and Mid-T was statistically significantly higher in VILI than in controls.

Conclusions

In this porcine model of VILI, we found a topographical correlation between lung strain and [68Ga]Ga-DOTA-TATE uptake at voxel level, suggesting that mechanical alteration and specific activation of inflammatory cells are strongly linked in VILI. This study represents the first voxel-by-voxel examination of this relationship in a multi-modal imaging analysis.

http://dx.doi.org/10.13039/501100003793 Hjärt-Lungfonden 20220536 20200841 20200877 20200825 20220681 20230767 Pellegrini Mariangela Perchiazzi Gaetano http://dx.doi.org/10.13039/501100003748 Svenska Sällskapet för Medicinsk Forskning 463402221 Pellegrini Mariangela http://dx.doi.org/10.13039/501100007687 Svenska Läkaresällskapet SLS-959793 Pellegrini Mariangela Swedish ALF Research FundALF-977974 ALF-977586 Pellegrini Mariangela http://dx.doi.org/10.13039/501100004359 Vetenskapsrådet 2018-02438 2020-02312 Eriksson Olof Perchiazzi Gaetano Alvar Gullstrand research grantALF-938050 Perchiazzi Gaetano Swedish Science for Life Laboratoryhttp://dx.doi.org/10.13039/501100006415 Stiftelsen Familjen Ernfors Fond Nils Erik Holmstens foundation for diabetes researchExoDiabEFSD/Lilly European Diabetes Research Programmehttp://dx.doi.org/10.13039/501100008316 P. O. Zetterlings stiftelse Uppsala DiabetescentrumHelmsley foundationhttp://dx.doi.org/10.13039/501100004973 Barndiabetesfonden http://dx.doi.org/10.13039/501100008550 Diabetesfonden Diabetes Wellness2409-PG Eriksson Olof Uppsala UniversityOpen access funding provided by Uppsala University.

issue-copyright-statement© European Society of Intensive Care Medicine and Springer Nature Switzerland AG 2024
==== Body
pmcIntroduction

In 1744 for the first time, it was observed that mechanical forces generated by uncontrolled manual ventilation could lead to injury and biological damage [1]. However, it is necessary to arrive at the implementation of mechanical ventilation on a routine basis to hypothesize the role of barotrauma as a complication [2]. Later, the focus moved from the high pressures exerted against the alveolar walls to the volumes that the lungs had to bear [3]. Still, it was necessary to come to the 2000s to convey several laboratory observations [4] into a unifying clinical hypothesis [5]. This, in turn led to two different research paths: the attempt to define safe limits of strain during mechanical ventilation [6] and the effort of measuring the elastic forces on a regional basis using image analysis of computer tomograms [7].

Different studies have reported increased inflammatory markers in animal models exposed to excessive lung strain during mechanical ventilation [8, 9]. However, these temporal associations never arrived to demonstrate the presence of inflammatory cells where a numerically estimated local strain of the lung was high: the studies of Borges and coworkers, based on Positron emission tomography (PET) with 2-deoxy-2-[18F]fluoro-d-glucose ([18F]FDG), although observed a concentration of the marker in specific areas of the lung, did not compute the local strain of these [10, 11].

A possible further complication has been the lack of a specific PET tracer for the cells intervening in the inflammatory process in the lung. The mentioned [18F]FDG is a glucose analog that labels all the cells using a glycolytic pathway: for these reasons, it requires a sophisticated data treatment for attenuating the effects of multicompartmental kinetics of the substance [12, 13]. A step forward has been performed by Bitker et al. [14, 15] using the tracer [N-methyl-11C]-1-(2-chlorophenyl)-N-(1-methylpropyl)-3-isoquinoline-carboxamide (from here on shortened as [11C]PK11195) which can mark the presence of macrophages. However, the generic presence of macrophages does not warrant the existence of an ongoing inflammatory process.

In fact, macrophages responding to environmental stimuli undergo a process of differentiation known as macrophagic polarization which results in a spectrum of diversely differentiated phenotypes, between the classically activated (M1) and alternatively activated (M2) macrophages. The M1 phenotype plays a pro-inflammatory role, producing cytokines that promote the acute phase of the inflammatory cascade, while the M2 macrophages play predominantly an immunoregulatory role, being involved in later stages of inflammation [16].

The discovery that activated, proinflammatory macrophages in phase M1 express on their surface the somatostatin receptors type 2 (SSTR2) [17] opened new perspectives in this field of research. The clinically available PET-tracer [68Ga] Ga-DOTA-TATE can avidly bind with high affinity to SSTR2 [18] and constitute a tool to analyze the spatial distribution of the inflammatory process related to lung injury at voxel level, as recently demonstrated by Puuvuori et al. [19].

We profited by a combined scanner performing positron emission tomography and computed tomography (PET–CT), using [68Ga]Ga-DOTA-TATE for activated macrophages and a computation method to estimate lung strain at the voxel level, to reveal whether activated macrophages are distributed in lung regions where high strain is located.

We hypothesized that high mechanical strain is topographically associated with the activation of the inflammatory cascade, as mirrored by the presence of polarized M1 macrophages. To assess this hypothesis, one group of animals experienced a dual-hit lung injury and was compared with a group of healthy controls. Both groups were ventilated using the same strategy, except for the limited duration of injurious ventilation in the VILI group. We analyzed the strain distribution to which the animal was exposed during injurious mechanical ventilation. Being computed tomography and positron emission tomography different modes of imaging, this purpose was performed using analysis methods able to observe whether information portrayed by one imaging mode was mirrored by the other.

Methods

The animal experiments were authorized by the Animal Research Ethics Committee (decision number 5.8.18-15648/2019) of Uppsala Region and carried out according to institutional guidelines (“Uppsala university guidelines on animal experimentation”, UFV 2007/724) and reported according to the ARRIVE guidelines [20]. The timeline of the experiment is presented in Fig. 1 and the outline of the experimental protocol is illustrated in Fig. 2.Fig. 1 Timeline of the experiment. After the arrival at the laboratory and the surgical preparation, one group of animals underwent the induction of lung injury by repeated lung lavages and injurious ventilation, while the healthy controls were ventilated with a tidal volume of 6 ml/kg and PEEP adjusted according to the ARDSnet FIO2/PEEP table (in the graph labeled as “protective ventilation”: see details and references in the main text). After this phase, for the remaining part of the experiment, all the animals were ventilated with 6 ml/kg and PEEP according to the ARDSnet table. T1, T2, and T3 denote the timing for the coordinated measurement of arterial blood gas analysis and hemodynamic parameters. T1 corresponds to baseline after animal preparation, T2 after VILI induction, and T3 just before the piglet was euthanized. Please note that the total time of ventilation in the two groups was the same: each animal arrived at the laboratory at the same time of the day and reached the PET facility at a fixed time in the afternoon. (Graph created with BioRender.com)

Fig. 2 Outline of the experiment. Phase 1. Two groups of animals were studied: one underwent a dual-hit lung injury composed of lung lavage followed by injurious mechanical ventilation; the other group consisted of healthy controls. Phase 2. All the animals underwent the computed tomography of the entire lung at different lung volumes between ZEEP and inspiratory capacity (IC). The map of lung strain distribution was obtained by subtracting voxel-by-voxel the gas content at ZEEP (here labeled as VFRC, volume at functional residual capacity—FRC, assuming that the end-expiratory lung volume at ZEEP can be considered practically equivalent to FRC in static conditions and in absence of intrinsic PEEP) from the gas content at IC. The difference was divided by VFRC: more details in the text. In phase 3, with the animal still in the PET–CT scanner, we administered the tracer [68Ga]Ga-DOTA-TATE, which has a high affinity for the SSTR2 receptors expressed on the surface of activated macrophages during inflammatory processes. The standardized uptake value ratios of the mentioned tracer allowed us to assess the regional distribution of activated macrophages. In phase 4, using information theory, we compared whether the activated macrophages foregather in the areas where the distribution of strain is higher (picture created with BioRender.com)

Anesthesia induction and animal preparation

Piglets (n = 10; 7 cases, 3 controls; weight 24.6 ± 1.6 kg) were pre-medicated through intramuscular injection of tiletamine 2.2 mg/kg and Zolazepam 6 mg/kg and lifted on the operatory table once unconscious. Subsequently, a peripheral intravenous line was established while peripheral blood saturation and 3 derivation ECG monitoring was applied. Fentanyl (1.5–2 mcg/kg) was administered intravenously before performing surgical tracheostomy and airway was secured with endotracheal cuffed tube (ETT) number 6 (6.0 Hi-Contour; Mallinckrodt Medical, Athlone, Ireland).

General anesthesia was maintained with intravenous Ketamine 32 mg/kg/h, Fentanyl 4 μg/kg/h, and Midazolam 0.12 mg/kg/h in Glucose solution 25 mg/mL. Muscle relaxation was achieved with intravenous Rocuronium 2.5 mg/kg/h. Mechanical ventilation was started with a tidal volume (VT) 8 ml/kg, respiratory rate (RR) 30 breath/min; positive end-expiratory pressure (PEEP) 5 cmH2O, FIO2 0.5. A high-flow central venous catheter, an arterial catheter, a pulmonary artery catheter, and a suprapubic urinary catheter were positioned via open surgical techniques. An esophageal catheter was positioned in the distal third of the esophagus to obtain a continuous measurement of esophageal pressure (Pes). Another balloon catheter was positioned in the gastric lumen, and continuous measurement of gastric pressure (Pga) was obtained. The correct position of the esophageal and gastric catheters was ascertained using a C-arm portable radiological machine. Hemodynamic advanced monitoring consisted of invasive arterial pressure, central venous pressure, pulmonary arterial pressure, wedge pressure, and cardiac output (this last by the technique of the cold bolus of saline executed in triplicate, injected randomly during the respiratory cycle). The animals were then ventilated for 30 min without performing any procedure to allow for stabilization.

Instrumental setup for lung mechanics measurements

Airway flow was measured continuously with a Fleisch pneumotachograph (Laminar Flow Element type PT, Special Instruments GmbH, Nördlingen, Germany) installed between the ETT and the ventilator, and connected to a differential pressure transducer (Diff-Cap Pressure Transducer, Special Instruments GmbH, Nördlingen, Germany). Pressure at airway opening, Pes, and Pga were measured continuously via three pressure transducers (Digimaclic Pressure Transducers, Special Instruments GmbH, Nördlingen, Germany). All respiratory signals were acquired via a data acquisition system (PowerLab PL 16/35, ADInstruments, Dunedin, New Zealand) controlled by the LabChart Software (ver. 7.0, ADInstruments, Dunedin, New Zealand). The inspiratory and expiratory airway volumes were obtained via flow integration.

Lung injury

Acute lung injury was induced in 7 of 10 animals using a “two-hit” procedure, combining lung lavages with warm saline (30 ml/kg sodium chloride 9 mg/ml at 37 °C), followed by injurious mechanical ventilation with high tidal volume and low PEEP [10]. Lung lavages were repeated to reach a PaO2/FIO2 ratio (P/F) ≤ 200 mmHg. Subsequently, injurious mechanical ventilation was started and maintained for a minimum of 30 min and a maximum of 60 min until a P/F ratio < 150 mmHg was reached. It consisted of pressure controlled mechanical ventilation with applied airway pressure 40 cmH2O; ZEEP (zero end-expiratory pressure); RR 30 breath/min; FIO2 = 1.

After the induction of VILI, protective ventilation was initiated (baseline ventilation) with a VT of 6 ml/kg. PEEP was adjusted according to the ARDSnet FIO2/PEEP table [21], aiming to a SpO2 target of 89–95%. The respiratory rate was set aiming to normocapnia. This was the ventilation mode throughout the entire experiment, with the temporary exception relative to the CT exposure at different inflation volumes. Animals belonging to the control group underwent directly protective ventilation without lung injury. The total time of mechanical ventilation was equal for both the control and VILI group of animals.

CT-imaging

To estimate the regional distribution of lung strain at different inflation volumes, we used an already-established procedure [7] that we summarize here following.

The first step was to measure spirometrically the inspired volume by applying a constant pressure of 40 cmH2O for 20 s, which was denoted as inspiratory capacity (IC). Based on this, CT scans of the lungs were performed at 12 different inflation volumes, corresponding to incremental steps of 1/12 of the IC. Subsequently, this value was divided into twelve volumetric steps (VT = (IC/12) × 1, (IC/12) × 2, (IC/12) × 3… up to IC) and noted down.

The animals were moved to the PET Center on a specially designed trolley, ventilated by the same ventilator and the baseline ventilation described before. After positioning on the PET–CT scanner bed in a supine position, they underwent RM at 40 cmH2O for 40 s to eventually restore the history of volumes. The twelve incremental VT steps were then delivered sequentially, each followed by an inspiratory hold maneuver (IHM) of approximately 10 s, during which an attenuation CT scan of the entire lung was acquired. Each IHM was followed by 2–3 min of tidal breathing to allow for the restoration of mechanical steady-state conditions. After IC was reached, two more CT scans were acquired during end-expiratory hold maneuvers (EHMs) at PEEP and ZEEP conditions (Fig. 2). The 64-slice CT scanner used a digital 4-ring system and had a 198-mm axial field of view (FOV) and the following parameters: 100 kV, 80–400 mA, noise index 10, rotation 0.5″, full spiral, slice thickness 3.75 mm, pitch 0.98:1, recon diameter 50 mm.

Measurements of hemodynamics

A further coordinated set of measurements was performed three times during the experiment: at baseline after animal preparation (T1), after VILI induction (T2), and just before the piglet was euthanized (T3). These measurements consisted of arterial blood gas (ABG) analysis and the above-mentioned hemodynamic measurements. At the end of the protocol, the animals were euthanized under deep anesthesia with an intravenous bolus of potassium chloride.

Lung strain regional distribution maps

For each animal, two transverse planes per inflation volume were analyzed. One corresponded to a para-diaphragmatic plane, located 2 cm cranially to the diaphragm; and one mid thoracic transverse plane, equidistant from lung apex and diaphragmatic dome.

Utilizing purposely written MatLab scripts by the authors (GP, MP), the perimeter of the lung parenchyma was manually contoured for each analyzed CT slice. This was done to avoid potential errors that automatic segmentation software may introduce, especially when analyzing consolidated lung areas [22]. To preserve the correspondence between voxels deriving from images shot at different inflation volumes, a lung image registration [23] procedure was applied before comparing images. It was based on purposely written scripts for the Image Processing Toolbox of MatLab (MatLab R2022b, The MathWorks, Natick, USA) and already validated in previous research by our group [22, 24].

For each image, a corresponding volume map was created by transforming the HU densities into gas volume at a voxel level, according to the concept proposed by Gattinoni et al. in 1988 [25]:1 Vvox,gas=Vvox·-HU1000,

where Vvox is the volume of a voxel, Vvox,gas is the volume of gas contained in the voxel, HU is the measure of the density in Hounsfield units. This method operates by approximating the density of a fully aerated area of lung parenchyma to air and the density of a fully consolidated area of lung parenchyma to water. Following this assumption, each voxel is assigned a gas content percentage that is directly proportional to its HU number in a spectrum between the values of HU for gas and water.

The Hookean strain is a measure of displacement of an elastic body from its equilibrium point: if the lung volume at Functional Residual Capacity (VFRC) is considered the elastic equilibrium point of the lung, the proportion of strain acting on the lung, at a given inflation volume (V1) can be written as [26]2 LungstrainV1=V1-VFRCVFRC.

We applied Eq. 2 to all the voxels corresponding to lung parenchyma. By hypothesizing that strain is (a) homogenously distributed inside each voxel and (b) all deformation corresponds to normal strain (i.e., perpendicular to the faces of the voxels), it is possible to quantify the voxel-associated Normal Strain (VoStrain) [27]:3 VoStraini,V1=13Vi,1-Vi,FRCVi,FRC.

This allowed us to represent the VoStrain associated with the different parts of the lung (field of strain, FOS) as computed from the CT scans.

For voxels belonging to completely atelectatic areas of the lung at FRC (i.e. presenting a null content of gas) an arbitrary gas content of 0.01% of voxel volume was assigned, to avoid values of VoStrain → ∞ in the computation of the field of strain. The strain that was analyzed during the successive part of the study was the one that the animals of the VILI group faced during the injurious mechanical ventilation, i.e., the strain developed ventilating between ZEEP and 40 cmH2O, while for the healthy controls, it was the strain corresponding to their end-inspiratory volume.

PET imaging

Immediately after having completed the acquisitions of CT images, the pigs were injected intravenously with a target dose of 2 MBq/kg of [68Ga]Ga-DOTA-TATE and underwent, in the same scanner (Discovery MI, GE Healthcare) a dynamic PET scan of 60 min (30 frames: 12 × 10 s″, 6 × 30 s″, 5 × 2 min′, 5 × 5 min′, 2 × 10 min′) while mechanically ventilated with the same pattern of breathing described before. The PET and CT acquisitions were automatically co-registered, as the pigs were not moved between the scans.

The PET images were generated by computing the standardized uptake value ratio (SUVR) [28] of the lung tissue, computing the ratio of radioactivity concentration (expressed as MBq/mL) between each lung voxel (cvoxel) at the coordinates x, y and the average concentration inside the left ventricle of the heart (cref) within the same PET image, at 3596 s (≈ 1 h) from the injection of [68Ga]Ga-DOTA-TATE:SUVRx,y=cvoxelx,ycref.

This allowed us to draw the regional distribution of SUVR at voxel level. The images depicting the SUVR distribution corresponding to the two chosen transversal planes of the CT-derived strain maps were further analyzed.

Image analysis and statistics

The field of strain and the SUVR of corresponding lung sections were further assessed. The images were analyzed using descriptive statistics (as mean ± standard deviation and median with interquartile range of all the animals and separately for VILI and control animals) and inferential statistics of FOS and SUVR. We could not exclude a non-normal distribution of the studied variables: for this reason, we used the nonparametric Wilcoxon signed rank test (α = 0.05) throughout the present manuscript, except where explicitly stated. The following question were addressed:Q1. Whether the SUVR of the PET images pertaining to controls and VILI have the same median;

Q2. Whether the FOS derived from CT images pertaining to controls and VILI have the same median;

After having divided all the PET and the FOS images into 100 consecutive iso-gravitational layers and having computed the means at each gravitational level from the four groups of animals (four conditions: mid-thorax or para-diaphragmatic; VILI or healthy controls) we tested whether:Q3. The gravitational course of the SUVR from PET images pertaining to controls and VILI are equal

Q4. The gravitational course of the FOS images pertaining to controls and VILI are equal

by performing the linear regression of the progressively leaning part of the curves (corresponding to the two central quartiles of the distribution, i.e. from slice 25 to 75) and evaluating whether the regression lines subtended different models by applying the F test [29].

To ascertain:Q5. Whether the local amount of activated macrophages is determined only by the distribution of blood due to gravitational forces,

We measured the relation between the volume of nongaseous component per each voxel (corresponding to the volume of blood or tissue) and SUVR, separately for diaphragm and for mid-thorax level, in healthy and VILI. This calculation of nongaseous content was performed using the method described by Pelosi et al. [30]. We computed the linear regressions between SUVR and nongaseous content deriving from the leaning portion of the two distribution profiles in the four conditions mentioned above. We compared the angular coefficients during VILI with the corresponding healthy conditions by computing the Z-statistic, according to Paternoster et al. [31].

Then, we compared the informative content between each pair of corresponding images, obtained with a different modality (e.g. FOS vs. SUVR pertaining to the same level, the same animal and separately for controls and VILI animals). For each couple of images, we computed the Linear Mutual Information (LMI) according to the method described by Lange and Grubmüller [32], in orderQ6. To verify that the information contained in one of the modalities was also present in the other.

By testing whether LMI was statistically different between the mentioned groups (Wilcoxon rank sum test, α = 0.05). The above-described image analysis and statistical testing were conducted using scripts specially developed by the authors GP and SS, for the Image Processing and the Statistics and Machine Learning Toolbox of MatLab (MatLab R2022b, The MathWorks, Natick, USA).

Results

All the animals survived the experiment. They presented at baseline a P/F ratio of 531.0 ± 43.4 mmHg. The pigs that underwent lung injury reached a P/F ratio of 79.0 ± 30.8 mmHg (p = 0.01) while the controls remained stable at P/F ratio 464.0 ± 44.5 mmHg (p = 0.25). The main cardiorespiratory data are presented in Table 1. The animals belonging to the VILI group were ventilated (during the phase of injurious ventilation) with a tidal volume of 436.6 ± 134 [ml], which corresponded to 18.0 ± 5 [ml/kg]. Table 1 Cardiorespiratory variables during the three main phases of the experiment in the two groups of healthy controls and animals exposed to ventilator-induced lung injury (VILI)

Time	T1	T2	T3	
Healthy	VILI	Healthy	VILI	Healthy	VILI	
Heart rate [bpm]	94.3 ± 13.1	98.9 ± 14.7	84.3 ± 5.3	88.6 ± 14.7	86.0 ± 3.7	95.0 ± 13.1	
Systemic systolic pressure [mmHg]	102.3 ± 13.6	88.7 ± 7.8	104.7 ± 10.9	86.1 ± 5.2*	103.0 ± 4.2	103.4 ± 11.2	
Systemic diastolic pressure [mmHg]	71.3 ± 11.6	57.3 ± 6.1	71.0 ± 7.8	56.0 ± 6.3*	63.3 ± 1.2	65.0 ± 7.3	
Pulmonary systolic pressure [mmHg]	22.7 ± 2.4	22.9 ± 3.8	28.0 ± 0.0	31.6 ± 10.3	28.7 ± 3.1	36.6 ± 8.0	
Pulmonary diastolic pressure [mmHg]	10.7 ± 3.7	10.7 ± 2.2	11.3 ± 0.9	17.6 ± 7.1	11.7 ± 3.4	24.1 ± 10.5	
Wedge pressure [mmHg]	9.0 ± 1.4	7.7 ± 1.8	9.3 ± 1.9	8.7 ± 2.3	9.7 ± 1.2	12.0 ± 2.4	
Central venous pressure [mmHg]	8.3 ± 1.2	6.0 ± 1.8	8.0 ± 1.4	7.7 ± 2.2	9.7 ± 2.9	9.7 ± 2.8	
Cardiac output [l/min]	2.9 ± 0.9	2.5 ± 0.5	2.7 ± 0.3	2.6 ± 0.6	3.1 ± 0.5	3.3 ± 0.4	
Tidal volume [ml]	243.3 ± 20.5	227.1 ± 8.8	190.0 ± 8.2	202.9 ± 16.0	203.3 ± 4.7	221.7 ± 22.7	
Respiratory rate [bpm]	31.7 ± 2.4	31.4 ± 3.5	35.0 ± 0.0	35.7 ± 3.2	35.0 ± 0.0	36.7 ± 2.4	
PEEP [cmH2O]	5.0 ± 0.0	5.0 ± 0.0	5.0 ± 0.0	8.7 ± 2.1*	5.0 ± 0.0	9.8 ± 1.5*	
Minute volume [l/min]	7.6 ± 0.8	7.5 ± 1.0	6.9 ± 0.2	7.4 ± 1.2	7.5 ± 0.3	8.3 ± 0.9	
Mechanical power [J/min]	11.1 ± 2.3	11.1 ± 2.2	9.8 ± 0.2	17.1 ± 5.0	10.5 ± 0.4	21.4 ± 4.0*	
The three times of data sampling were: at baseline after animal preparation (T1), after VILI induction (T2), and just before the piglet was euthanized (T3). The asterisk marks a statistically significant difference (Student t test, α = 0.05) between healthy controls and VILI

Their total lung compliance went from 21.6 ± 4.9 to 12.2 ± 2.6 [ml/cmH2O] in the animals subjected to lung injury (p = 0.01) and from 25.6 ± 6.5 to 25.3 ± 6.0 (p = 0.75) for the controls.

The analysis of respiratory mechanics showed that when all the animals were examined at the same reference volume of (IC/12) × 6, corresponding to 273.6 ± 87.8 [ml], the healthy controls showed a Pes 14.7 ± 5.0 [cmH2O], a transpulmonary pressure (PTP) of 4.0 ± 5.9 [cmH2O] and a driving pressure (ΔP) of 13.7 ± 1.6 [cmH2O], while the animals subjected to VILI presented Pes 13.1 ± 2.7 [cmH2O], PTP of 12.9 ± 4.3 [cmH2O] and ΔP of 18.0 ± 4.0 [cmH2O]. Statistical analysis comparing healthy controls with VILI group reported a significant difference in terms of PTP (p = 0.03), while the difference in Pes and ΔP did not reach statistical significance (respectively p = 0.58 and p = 0.13).

The computation of strain was performed between lung volume at ZEEP and the highest inflation volume to which the animals were exposed during the induction of VILI, corresponding to the inspiratory capacity. By this way the experiment produced 21,300 CT scans, we selected 240 CT slices and brought to the final analysis 40 images (10 animals × 2 slice levels × 2 inflation volumes), see Fig. 3.Fig. 3 Representative picture of standardized uptake value ratio (SUVR), voxel-associated normal strain (VoStrain) distribution and computed tomography (CT) in one animal of the VILI group and one healthy control, at para-diaphragmatic and mid-thoracic levels. SUVR and strain are dimensionless measures. HU: hounsfield number. Please note that the color scale of the VoStrain is different between VILI and healthy controls: the regional voxel-associated normal strain in VILI can reach levels three orders of magnitude higher than in healthy controls

Strain

We computed the VoStrain to which the animal was exposed during the induction of VILI, using the CT images obtained at inspiratory capacity and the CT images taken at ZEEP. This was calculated at para-diaphragmatic level and at mid-thoracic level. Data are reported as mean ± standard deviation, followed by (median; interquartile range).

For the VILI group, average VoStrain over the entire slice was 315 ± 79 (284; 181) at para-diaphragmatic level and 185 ± 165 (122; 155) at mid-thoracic level. For the control group, it was respectively: 108 ± 6 (104; 52) and 24 ± 4 (21; 22).

Plotting the profile of VoStrain along the gravitational axis (x-axis in Fig. 4) it is possible to observe that it reaches a maximum of 760 ± 210 (812; 268) at para-diaphragmatic level, while at mid-thoracic level max value was 497 ± 373 (668; 582); in the controls it was respectively: 449 ± 106 (485; 131) and 193 ± 160 (255; 101).Fig. 4 VoStrain in the different groups of analyzed images: in VILI and healthy controls, at para-diaphragmatic and mid-thoracic levels. The graphs represent the mean ± standard error obtained by averaging the 100 iso-gravitational slices of the lung across the different animals. The dashed line depicts the linear regression of the leaning part of the profiles, for homogeneity identified as the central quartiles (from slice 25 to 75) in all the graphs. The graph also reports the respective regression lines and regression coefficients. ε represents VoStrain and it is a dimensionless variable

The gravitational course of strain, analyzed using regression lines, showed significant differences between healthy and VILI conditions at both mid-thorax (p < 0.01) and para-diaphragmatic levels (p < 0.01).

PET

The animals received an average dose of 2.3 ± 0.9 MBq/kg of [68Ga]Ga-DOTA-TATE. After the PET study, we gathered a total of 2130 images (10 animals × 30 times × 71 sections). The following data are reported as mean ± standard deviation, followed by (median; interquartile range). The mean SUVR for entire images at para-diaphragmatic level in healthy conditions was 0.6 ± 0.0 (0.6; 0.0); after lung injury it was 1.1 ± 0.6 (0.9; 0.7); at thorax level, it was, respectively, 0.3 ± 0.0 (0.3; 0.1) in healthy and 0.5 ± 0.3 (0.3; 0.4) after injury.

The topographic distribution of SUVR (see Fig. 5) of [68GA]GA-DOTA-TATE was higher in the dependent regions and reached a maximum of 1.3 ± 0.1 (1.3; 0.1) at para-diaphragmatic level in controls; at the same level, it was 2.2 ± 1.3 (1.8; 1.6) during VILI; at mid-thorax it was 0.6 ± 0.03 (0.7; 0.0) in controls and 1.3 ± 1.0 (0.8; 1.3) in VILI.Fig. 5 Standardized uptake value ratios (SUVR) in the different groups of analyzed images: in VILI and healthy controls, at para-diaphragmatic and mid-thoracic levels. The graphs represent the mean ± standard error obtained by averaging the 100 iso-gravitational slices of the lung across the different animals. The dashed line depicts the linear regression of the leaning part of the profiles, for homogeneity identified as the central quartiles (from slice 25 to 75) in all the graphs. The graph also reports the respective regression lines and regression coefficients. SUVR is a dimensionless variable

The profiles of the regression lines describing the increase of SUVR along the gravitational axis showed that at both mid-thoracic and para-diaphragmatic levels, the curves followed statistically different regression lines between VILI and healthy conditions at both mid-thoracic and para-diaphragmatic level, both expressed by p < 0.01.

SUVR vs. nongaseous content of the images

The course of the nongaseous content is presented in Fig. 6 and, as expected, increases progressively in the most dependent regions of the lung. The regression lines (see Fig. 7) of SUVR (y) vs. amount of nongaseous component (x) yielded for para-diaphragmatic level in VILI: y = 15.68x − 0.61 (R2 = 0.90); for para-diaphragmatic level in healthy: y = 11.70x − 0.17 (R2 = 0.88); for mid-thorax in VILI: y = 12.09x − 0.09 (R2 = 0.97); for mid-thorax in healthy: y = 9.86x − 0.02 (R2 = 0.98). All the regression equations were statistically significant as well as there was a statistically significant difference (by applying the Z test on the angular coefficients) between healthy and VILI at both diaphragm level (p = 1.6 × 10–5) and mid-thorax level (p = 6.8 × 10–9).Fig. 6 Nongaseous content in the different groups of analyzed images: in VILI and healthy controls, at para-diaphragmatic and mid-thoracic levels. The graphs represent the mean ± standard error obtained by averaging the 100 iso-gravitational slices of the lung across the different animals. Nongaseous content is represented in microliters [µl]

Fig. 7 Scatter plot of SUVR vs. nongaseous content, deriving from the same iso-gravitational slices, sampled in the leaning part of both distribution profiles. The two panels compare VILI vs. healthy controls at para-diaphragmatic and mid-thorax levels. The respective linear regression equations are also represented (regression lines for VILI are represented in red and for healthy in green). The rate of increase of SUVR and non-gaseous content along the gravitational axis is represented by the angular coefficients of the regressions. For clarity, on the right of the scatter plots, it is reported the relation between the same regression lines, imposing an intercept of zero, to facilitate the comparison of the angular coefficients. The Z tests performed on both levels (VILI vs. healthy at para-diaphragmatic and mid-thoracic levels) yielded a statistically significant difference between angular coefficients (marked with “*”), demonstrating that the rate of increase of SUVR is higher in VILI than in healthy. This means that at the same amount of nongaseous content in the lung (that in the VILI group may derive from excess tissue as in atelectasis or from a higher local perfusion due to vasodilation and inflammation), the concentration of macrophages that are present on-site is higher than the one that can be considered “passively carried” by blood flow

Mutual information

Computation of the linear mutual information between FOS and SUVR images yielded a global 0.734 ± 0.144 during VILI and for controls 0.420 ± 0.162. Their difference was statistically significant with p = 6·10–4. In particular, at the para-diaphragmatic level, it presented 0.793 ± 0.108 in VILI and 0.564 ± 0.046 in healthy conditions (statistically different, p = 0.017); at the mid-thorax level, it was: during VILI 0.676 ± 0.159 and in healthy conditions: 0.275 ± 0.017 (statistically different, p = 0.017).

Discussion

This study demonstrated that in an animal model of VILI induced by surfactant depletion and injurious ventilation, activated macrophages plausibly in phase M1 accumulate in the areas that were subjected to high local strain. These areas can be located according to the classical ventilation distribution of the lung or in proximity of closed/atelectatic parenchyma because exposed to high traction force [33]. This is the first time this relation has been examined on a voxel-by-voxel basis in a multi-modal image analysis.

First, we have separately analyzed the distribution of strain and SUVR for [68Ga]Ga-DOTA-TATE. We have observed that for both the techniques, the images pertaining to healthy controls and VILI had different medians, presenting the images in VILI a higher strain and a higher SUVR than healthy lungs (hypothesis Q1 and Q2).

Then, we studied separately the distribution of both strain and SUVR along the gravitational axis, by dividing the images in 100 isogravitational levels. With both methods (SUVR and FOS) there was a statistically significant difference between the regression lines subtending the behavior in healthy and VILI conditions (hypothesis Q3 and Q4).

The meaning of these findings resides in the fact that the strain and the SUVR increase in the most dependent areas of the lung, growing differently in healthy and VILI conditions. The gravitational increase of strain confirms the well-established concept that ventilation is higher in the most dependent areas of the lung [34], as well as the hypo- or non-ventilated areas, so consequently, there is the potential for a higher strain [7]. Moreover, the difference between the inclinations of the two regression lines of SUVR in healthy and VILI conditions (Fig. 7), suggests that the higher presence of macrophages in the dependent areas is not a simple consequence of gravitational forces. To confirm this, we have measured the regression lines between SUVR and the amount of nongaseous component of each voxel (composed of fluids, like blood and of tissues, like atelectatic lung). If the distribution of macrophages (mirrored by the SUVR) followed only a gravitational distribution, the ratio between SUVR and nongaseous component should have followed the same equation line, independently if the overall amount of SUVR was high or low. However, the comparison has demonstrated that the equations subtending SUVR vs. nongaseous component are different between healthy and VILI condition, supporting the idea that the accumulation of macrophages is not a simple consequence of gravitational forces.

We hypothesize that the SUVR of [68Ga]Ga-DOTA-TATE portrays activated macrophages polarizing toward the M1 phenotype and reflects an ongoing inflammatory process. It is still matter of debate the proportion between resident (intra-alveolar and interstitial) macrophages and the latterly arrived on the site of inflammation by the recruitment of circulating monocytes [35]. A classical view is that the overall amount of blood is higher in the gravitational dependent areas [36], however, is equally clear that the regional distribution of blood depends also on local factors [37]. This means that the arrival of recruited monocytes–macrophages can follow also other pathways than the gravity-paved one.

Strain derived from CT images is a measure of the displacement caused by the force field, and PET depicts the distribution of a tracer inside the lung.

In a second phase of our study, to make a voxelwise comparison between the distribution of strain and of activated macrophages, we applied Shannon’s mutual information laws, included in Lange and Grubmüller method [32]. To compare these couples of images (strain map derived by CT and macrophage map obtained by SUVR and [68Ga]Ga-DOTA-TATE) it is necessary not only to see whether there is a generic association between them but also to measure whether they express the same informative content.

Many methods are available [38] to perform this task. Shannon’s mutual information is one of the most efficient and most widely used methods for calculating the similarity between images obtained with different modalities. Based on Shannon’s theorem of information [39] and qualitatively described, it yields the probability that the intensity of corresponding pixels (in two images obtained with a different technique) shows a coordinated variation, thus depicting the same objects. Measurement of mutual information in our images demonstrated that activated macrophages converge in high numbers where the strain is high (Q6). In fact, during VILI, PET images detect many more activated macrophages that gather in the same position where there is high strain. For this reason, these “flashes of light” that PET images show during VILI inform the viewer that a higher intensity is present in the corresponding strain maps. This information is mutual: analyzing strain map intensity reveals the location of activated macrophages.

Computation of voxel-associated strain

Differently from the “spirometrical” strain computed at airways opening [40], the voxel-associated strain (VoStrain) is computed from a change of density determined by a CT [41, 42]

Assuming that each voxel subtends anatomical structures that have the same behavior, the voxel-associated strain yields the amount of strain that is distributed on each of the three axes of the voxel and is defined by classical physics of elastic materials [27]. Assigning an amount of strain to each axis makes it possible to deduce the magnitude of deformation to which molecules composing the so-called stress-bearing elements are exposed [43] and eventually brought to mechanical failure [44].

The structures encompassed in one voxel include alveoli, gas, vessels, fluids, and interstitial tissues. If we assign to the alveoli the role of loading the inspiratory gas, we should affirm that the computed strain must be attributed to their walls. However, inside one typical CT voxel, there are hundreds of alveoli [45] and their inflation process is not linear. There are sequential processes of recruitment and inflation [46], and moreover, after the alveolus starts its inflation, the amount of strain born by the alveolar walls depends on the morphological alteration of its shape [47] rendering the stress/strain relation of the single groups of alveoli as nonlinear [48]. In fact, the process of sequential recruitment increases the number of available alveoli and make it possible to subdivide the inflation-related strain between newly opened units.

Distribution pro-inflammatory M1 macrophages using PET

[68Ga] Ga-DOTA-TATE is a clinically available radiopharmaceutical binding to SSTR2 and widely used for neuroendocrine tumors imaging. The SSTR are normally expressed on endocrine organs, spleen, thyroid, adrenal, and pituitary glands, while they show little to no expression in healthy lungs.

Among somatostatin receptors, [68Ga]Ga-DOTA-TATE has the highest affinity for somatostatin receptors type 2 [18], widely expressed specifically by pro-inflammatory M1 macrophages [17], and for this reason, has already been used successfully to detect macrophagic activation in atherosclerosis-related inflammatory processes [49–51].

The macrophages, both resident in the alveoli and recruited from the bloodstream, have a pivotal role during the initial phase of lung inflammation and ARDS [52, 53]. For this reason, the [68Ga]Ga-DOTA-TATE is an ideal candidate for investigating the presence of activated M1 alveolar macrophage during ARDS.

In healthy lungs, long-lived resident alveolar macrophages are the predominant population in the alveoli and are mainly expressing immunosuppressive M2 characteristics. In response to noxious stimuli of the lung, however, circulating monocytes are recruited into the alveolar space, undergo an M1 polarization and thus secrete cytokines and mediators that promote a pro-inflammatory response. This process happens in the early stages of the pathogenesis of ARDS, also known as exudative phase, when mediators produced by M1 macrophages, attract neutrophils from the intravascular space, contributing to the elimination of noxious factors, but also to the tissue damage seen in ARDS. In the later stages of ARDS, after the pathogenic stimuli are eliminated, a gradual shift of resident and recruited macrophages from the M1 to the M2 phenotype occurs. At this stage, M2 macrophages play a key role in mitigating the immune response, producing anti-inflammatory cytokines, and participating in lung tissue repair by phagocytizing debris of necrotic and apoptotic cells [54, 55].

Our group has already tested the efficacy of [68Ga]Ga-DOTA-TATE to signal the presence of inflammation during lung injury in a porcine model of ARDS, comparing it with [18F]FDG–PET and with histology [19]. It is worth noting that SSTR2 receptors are actively expressed only in M1 pro-inflammatory macrophages [17], while other SSTR subtypes are not upregulated on the M1 polarized macrophages.

This makes the [68Ga]Ga-DOTA-TATE highly specific for detecting the presence of activated macrophages and undoubtedly gives it an advantage over the FDG in detecting cell-mediated early inflammation in acute lung injury.

18F-FDG is a glucose analog and lacks specificity, although it is acknowledged that 18F-FDG may accumulate in activated macrophages [56] than in other subtypes because of the higher glycolytic activity in the pro-inflammatory population. However, it should be mentioned that the increase of glycolytic activity also varies depending on the type of stimulation received by macrophages, with stimulation via the alternative pathway resulting in fewer metabolic effects [57].

Bitker and colleagues have studied lung inflammation by measuring the presence of macrophages using the different PET-tracer named [11C]PK11195 which is a ligand of the so-called translocator protein (TSPO). TSPO is a mitochondrial receptor found in macrophages and neutrophils [58]. Although associated with the state of lung inflammation by tissue examination by Bitker [15], this protein is expressed by macrophages in both conditions M1 (proinflammatory) and M2 (anti-inflammatory). Moreover, a recent paper reports that the TSPO is downregulated on pro-inflammatory ‘M1’ human macrophages [59]. This makes [11C]PK11195 not genuinely specific for inflammation per se, but surely able to generally signal the presence of macrophages.

Critical appraisal of the methods

The lung injury model

We used a well-established dual-hit VILI model [10, 60, 61]. The advantage of applying in sequence the two-hits injury resides mainly in the steadiness of the alterations in lung mechanics and gas exchange during the course of the experiment [60]. To prevent differences associated with varying mechanical ventilation times, the total ventilation time was kept equal in both animal groups.

Computation of voxel-associated strain is surely one of the strengths of the present paper. However, the assessment of strain with such a high granularity warrants important reflections. CT scans portray voxels subtending alveoli in a broad range of inflation volumes and not necessarily “functionally connected” to the main airways. In particular, the denominator of Eqs. 2 and 3 can reach very low numbers (i.e. when at FRC the corresponding voxel subtends hypo-inflated or even non-inflated alveoli) and the consequently computed strain reaches very high values. The problem posed by non-ventilated areas in the computation of voxel-associated strain on the lung during mechanical ventilation has been already encountered in other studies. Already in 2010, Caironi and others [62] hypothesized a relation between the potentially recruitable lung and the end-inspiratory alveolar strain, by adjusting the denominator of strain equation for the recruited volumes by assuming their CT density on the basis of preceding observations. The same problem was acknowledged by Paula et al. [63] who adjusted strain calculation using the experimental values deriving from average measurements of regional aeration. However, the use of the recruited volume for adjusting the denominator of the strain equations (see Eqs. 2 and 3) is also prone to generate flaws in strain computation. In our study, we decided not to correct for the recruited volumes and to present the original, untreated data. This has a scientific rationale behind it. There is a general consensus that strain should be computed using as a reference the unstressed volume of lung structures [5]. At the microscopic level an atelectatic structure is certainly in an unstressed state; however, when gas start entering a group of alveoli (that would allow it to be classified as recruited area), it does not guarantee that the equilibrium point between inward and outward forces is reached. In fact, as pointed out by the studies of micromechanics [48], although groups of alveoli can start “loading” gas (and becoming “recruited” by definition), they are still in a process of unfolding their scaffolding and still unable to generate any elastic recoil [64]. This means that the use of the recruited volume for adjusting the denominator of Eqs. 2 and 3 is also by itself prone to generate flaws in strain computation. In contrast, when considering the translational applicability of strain calculation in a clinical context, it is important to recognize that the modification of end-expiratory lung volume through PEEP application can alter the strain exerted on lung structures.

In this context, it is correct to take into account the possibility that animals in the VILI group, ventilated after VILI induction with a tidal volume of 6 ml/kg and PEEP adjusted according to the ARDSnet FIO2/PEEP table [21], may have experienced an additional strain-related inflammation due to higher levels of PEEP applied during ventilation (see Table 1). However, it is worth noting that several translational observations [65] suggest that high dynamic strain (applied during injurious ventilation in the VILI group) causes more evident damage than the steady application of a high PEEP.

Due to the study design, which involved two defined groups (healthy controls vs. animals exposed to VILI), it was not possible to quantitatively analyze a dose–response relationship between strain and inflammation. This type of analysis would have required a larger sample size and potentially the induction of varying degrees of damage at the beginning of the experiment.

In our model, during VILI induction, the animals were exposed to remarkable strain levels (corresponding to repeated step changes from ZEEP to 40 cmH2O) that can be considered very high when compared with the pattern of ventilation customarily used in clinical routine. With the premises deriving from the choice of the model and the computation method described before, the calculation of voxel-associated strain using CT imaging surmounts the value of strain measured at airway opening [6].

PET-tracer

The non-specificity of [18F]FDG obliges complex calculations of the tracer uptake to assess the contribution of the active use of glucose by the cells involved the inflammation cascade and differentiate from other cells that, although using glucose, are transient bystanders. Instead, for [68Ga]Ga-DOTA-TATE, we could apply a plain SUVR calculation. Being in our lung PET slice the activated macrophages and the [68Ga]Ga-DOTA-TATE tracer either in the lung structure or in the blood, the SUVR computation [28] allowed us to know the proportion of [68Ga]Ga-DOTA-TATE in relation to its concentration inside the heart chambers representing the blood pools (there is no uptake of [68Ga]Ga-DOTA-TATE in the myocardium walls which may confound the segmentation of the ventricle [19], as can happen with [18F]FDG). We shot a series of PET images per each transverse plane and we chose a timing for their analysis corresponding to 3596 s (≈ 60 min) in line with the available literature and our previous tests [19].

The visualization of macrophages during inflammation can also benefit from the use of tracers other than [68Ga]Ga-DOTA-TATE. Apart from the previously discussed [11C]PK11195, promising results have also been achieved using markers that target the folate receptor β in studies of myocarditis [66] or rheumatoid arthritis [67]. Studies about the specificity and sensitivity of [68Ga]Ga-DOTA-TATE as a diagnostic test for detecting pulmonary inflammation are lacking. More studies have been conducted to assess the diagnostic performance of [68Ga]Ga-DOTA-TATE in detecting other kinds of inflammation processes, like atherosclerosis. [68Ga]Ga-DOTA-TATE demonstrated a better power than [18F]FDG to discriminate high-risk from low-risk coronary atherosclerotic lesions, allowing for a clear interpretation of coronary signals [51]. Even if there are no studies that have directly compared [68Ga]Ga-DOTA-TATE and [18F]FDG as diagnostic methods to assess pulmonary inflammation, a possible hypothesis that is worth investigating has been enunciated by Li et al. [49]. They affirm that the main difference lies in tracer kinetics. [68Ga]Ga-DOTA-TATE binds specifically to SSTR2, and its uptake is limited by the number of saturated SSTR2 receptors. On the other hand, [18F]FDG is continuously metabolized by activated macrophages and is more susceptible to nonspecific uptake mechanisms, such as increased vascular permeability due to inflammation. Therefore, theoretically, [18F]FDG should be more sensitive than [68Ga]Ga-DOTA-TATE, but potentially less specific. As already said, this hypothesis needs to be confirmed by a purposely designed study.

Shannon’s mutual information

Mutual information measures the amount of information that one variable contains about another or the reduction in uncertainty of one variable when we know the other [68]. Although nowadays extensively used as a tool for image registration [69], mutual information was originally introduced as a similarity measure between images [70]. In particular, the rationale for using it in the present study is that no assumption is needed of a prior functional relationship between the images, it is robust towards outliers, and can be calculated efficiently [68]. The absence of assumptions about the relationship between the image intensities of the two studied modalities makes mutual information highly versatile and powerful without the need for prior segmentation [70].

In its general definition, computation of Shannon’s mutual information yields values that might vary notably in relation to the size and the content of the analyzed images [38]. This may render difficult the comparison of two series of images instead of a couple of them. For this reason, many attempts have been proposed to “normalize” the output of Shannon’s mutual entropy to a fixed range [38] to compare different pools of data and to simplify their interpretation. A mathematical solution to this problem has been published some years ago by Lange and Grubmüller [32] in the field of molecular dynamics, with the scope of avoiding the original range of the output from the original Shannon’s equations [between zero and infinite], which they defined “unfamiliar” and without obvious interpretation. To solve this, they proposed the method of calculation that yields results in the range of [0…1] that we have used in the present paper and allowed the reported comparisons.

The shear strain

Estimation of strain from density changes derived from CT does not allow to compute the so-called shear strain (the strain parallel to the surfaces of the voxel), although in terms of energy dissipation, this component can play a non-negligible role.

Clinical perspective

It is theoretically reasonable to use a PET tracer in combination with a PET–CT scanner to detect inflammation and diagnose active lung injury in patients admitted to intensive care. The use of PET–CT methods has been increasing in recent years [71], but more research is needed before it can become part of routine clinical practice. It is important to conduct extensive studies to test the specificity and sensitivity of new tracers for each lung inflammatory disease that requires a PET diagnosis. In addition, each new method should undergo a Health Technology Assessment evaluation before its introduction [57].

Conclusion

Using PET–CT image analysis, we have been able to compute the field of lung strain and to mark the position of activated macrophages using the selective SSTR2 PET-ligand [68Ga]Ga-DOTA-TATE. We have analyzed the information provided by these imaging methods first separately and then in a joint fashion, by applying a method based on the theory of information. Macrophages, among the first inflammatory cells to be activated in the place of inflammation, gather in the areas where strain induced by mechanical ventilation is high.

For the first time we show a strong regional coupling between high mechanical strain and cell-mediated inflammation. These findings may have future clinical relevance for both monitoring of VILI-associated inflammation and its potential treatments (e.g., for targeted anti-inflammatory therapies of ARDS).

Abbreviations

[18F]FDG 2-Deoxy-2-[18F]fluoro-d-glucose

ABG Arterial blood gas

clung Radioactivity concentration in voxel belonging to the lung

cref Radioactivity concentration in a reference group of voxels

CT Computed tomography

ΔP Driving pressure

ETT Endotracheal tube

FIO2 Fraction of inspired oxygen

FOS Field of strain

HU Hounsfield units

IC Inspiratory capacity

IHM Inspiratory hold maneuver

MBq Mega Becquerel

Mid-T Computed tomography slice at mid-thoracic levels

Para-D Computed tomography slice at para-diaphragmatic level

PEEP Positive end-expiratory pressure

Pes Esophageal pressure

PTP Transpulmonary pressure

PET Positron emission tomography

PET–CT Combined positron emission tomography and computed tomography

Pga Gastric pressure

[11C]PK11195 [N-Methyl-11C]-1-(2-chlorophenyl)-N-(1-methylpropyl)-3-isoquinoline-carboxamide

RM Recruitment maneuver

RR Respiratory rate

SpO2 Peripheral oxygen saturation

SSTR2 Somatostatin receptor subtype 2

SUVR Standardized uptake value ratios

VFRC Lung volume at functional residual capacity

VILI Ventilator-induced lung injury

VoStrain Voxel-associated normal strain

VT Tidal volume

Vvox,gas Volume of gas contained in a voxel

Vvox Volume of one voxel

ZEEP Zero end-expiratory pressure

Acknowledgements

The authors express their sincere gratitude to Kerstin Ahlgren, Agneta Roneus, Liselotte Pihl, Mariette Andersson and Maria Swälas for their assistance and support during the experiments at the Hedenstierna Laboratory. The authors would also like to acknowledge the grateful assistance of the staff at the PET center at the Uppsala University Hospital.

Author contributions

FL, EP, EC performed the experiments and preliminary data analysis; MP, DC, LB, PP, GA performed data analysis and drafted the manuscript; IV performed the radiochemistry. JS and EP performed the imaging. GP wrote the image analysis and statistical scripts. SS wrote the scripts for Shannon Mutual Information test and analyzed its results. OE and GP designed the study, performed the experiments, analyzed data, and wrote the first draft of the manuscript, which was then revised by all authors for structure and scientific content. Finally, all authors read and approved the final manuscript.

Funding

Open access funding provided by Uppsala University. MP received research funding from the Swedish Heart and Lung Foundation (20220536, and 20200841), the Swedish Society for Medical Research (463402221), the Swedish Society of Medicine (SLS-959793), and the Alvar Gullstrand research grant, and ALF research fund (ALF-977974, ALF-977586). GP received grants from the following: the Swedish Research Council (2018-02438), the Swedish Heart and Lung Foundation (20200877, 20200825, 20220681, 20230767), and the Alvar Gullstrand research grant (ALF-938050). OE was supported by grants from Science for Life Laboratory, the Swedish Research Council (2020-02312), the Ernfors Family Foundation, the Nils Erik Holmstens Foundation for diabetes research, ExoDiab, the EFSD/Lilly European Diabetes Research Programme, PO Zetterlings Foundation, Uppsala Diabetescentrum, the Helmsley foundation, Barndiabetesfonden, Diabetesfonden, and Diabetes Wellness (#2409-PG).

Availability of data and materials

The data set analyzed during the current study are available from the corresponding author on reasonable request.

Declarations

Ethics approval and consent to participate

The animal experiments were authorized by the Animal Research Ethics Committee (decision number 5.8.18-15648/2019) of Uppsala Region and carried out according to institutional guidelines (“Uppsala university guidelines on animal experimentation”, UFV 2007/724).

Consent for publication

Not applicable.

Competing interests

OE is an employee and minority shareholder of Antaros Tracer AB. However, the mentioned company had no role in the planning, execution, and reporting of the present manuscript and is not involved in the production of the tracer object of the present research. MP serves on the journal's editorial board. Otherwise, the remaining authors have no competing interests.

Publisher's Note

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

Francesco Liggieri and Elena Chiodaroli are first authors and contributed equally to this work.

Olof Eriksson and Gaetano Perchiazzi are senior authors and contributed equally to this work.

Paolo Pelosi—deceased.
==== Refs
References

1. Fothergill J Observations on a case published in the last volume of the Medical Essays & c., of recovering a man dead in appearance, by distending the lungs with air Philosophical transactions of the Royal Society of London 1744 London Royal Society of London 275 281
Fothergill J (1744) Observations on a case published in the last volume of the Medical Essays & c., of recovering a man dead in appearance, by distending the lungs with air. Philosophical transactions of the Royal Society of London. Royal Society of London, London, pp 275–281
2. Kumar A Pontoppidan H Falke KJ Wilson RS Laver MB Pulmonary barotrauma during mechanical ventilation Crit Care Med 1973 1 181 186 10.1097/00003246-197307000-00001 4587509
Kumar A, Pontoppidan H, Falke KJ, Wilson RS, Laver MB (1973) Pulmonary barotrauma during mechanical ventilation. Crit Care Med 1:181–186. 10.1097/00003246-197307000-000014587509 10.1097/00003246-197307000-00001
3. Dreyfuss D Saumon G Barotrauma is volutrauma, but which volume is the one responsible? Intensive Care Med 1992 18 139 141 10.1007/BF01709236 1644960
Dreyfuss D, Saumon G (1992) Barotrauma is volutrauma, but which volume is the one responsible? Intensive Care Med 18:139–141. 10.1007/BF017092361644960 10.1007/BF01709236
4. Liu M Tanswell AK Post M Mechanical force-induced signal transduction in lung cells Am J Physiol Cell Mol Physiol 1999 277 L667 L683 10.1152/ajplung.1999.277.4.L667
Liu M, Tanswell AK, Post M (1999) Mechanical force-induced signal transduction in lung cells. Am J Physiol Cell Mol Physiol 277:L667–L683. 10.1152/ajplung.1999.277.4.L66710.1152/ajplung.1999.277.4.L667
5. Chiumello D Carlesso E Cadringher P Caironi P Valenza F Polli F Tallarini F Cozzi P Cressoni M Colombo A Marini JJ Gattinoni L Lung stress and strain during mechanical ventilation for acute respiratory distress syndrome Am J Respir Crit Care Med 2008 178 346 355 10.1164/rccm.200710-1589OC 18451319
Chiumello D, Carlesso E, Cadringher P, Caironi P, Valenza F, Polli F, Tallarini F, Cozzi P, Cressoni M, Colombo A, Marini JJ, Gattinoni L (2008) Lung stress and strain during mechanical ventilation for acute respiratory distress syndrome. Am J Respir Crit Care Med 178:346–355. 10.1164/rccm.200710-1589OC18451319 10.1164/rccm.200710-1589OC
6. Protti A Cressoni M Santini A Langer T Mietto C Febres D Chierichetti M Coppola S Conte G Gatti S Leopardi O Masson S Lombardi L Lazzerini M Rampoldi E Cadringher P Gattinoni L Lung stress and strain during mechanical ventilation: any safe threshold? Am J Respir Crit Care Med 2011 183 1354 1362 10.1164/rccm.201010-1757OC 21297069
Protti A, Cressoni M, Santini A, Langer T, Mietto C, Febres D, Chierichetti M, Coppola S, Conte G, Gatti S, Leopardi O, Masson S, Lombardi L, Lazzerini M, Rampoldi E, Cadringher P, Gattinoni L (2011) Lung stress and strain during mechanical ventilation: any safe threshold? Am J Respir Crit Care Med 183:1354–1362. 10.1164/rccm.201010-1757OC21297069 10.1164/rccm.201010-1757OC
7. Perchiazzi G Rylander C Vena A Derosa S Polieri D Fiore T Giuliani R Hedenstierna G Lung regional stress and strain as a function of posture and ventilatory mode J Appl Physiol 2011 110 1374 1383 10.1152/japplphysiol.00439.2010 21393463
Perchiazzi G, Rylander C, Vena A, Derosa S, Polieri D, Fiore T, Giuliani R, Hedenstierna G (2011) Lung regional stress and strain as a function of posture and ventilatory mode. J Appl Physiol 110:1374–1383. 10.1152/japplphysiol.00439.201021393463 10.1152/japplphysiol.00439.2010
8. Vlahakis NE Schroeder MA Limper AH Hubmayr RD Stretch induces cytokine release by alveolar epithelial cells in vitro Am J Physiol Lung Cell Mol Physiol 1999 277 167 173 10.1152/ajplung.1999.277.1.l167
Vlahakis NE, Schroeder MA, Limper AH, Hubmayr RD (1999) Stretch induces cytokine release by alveolar epithelial cells in vitro. Am J Physiol Lung Cell Mol Physiol 277:167–173. 10.1152/ajplung.1999.277.1.l16710.1152/ajplung.1999.277.1.l167
9. Laffey JG Kavanagh BP Insight into acute respiratory distress syndrome: from models to patients Am J Respir Crit Care Med 2017 196 18 28 10.1164/rccm.201612-2415CI 28146637
Laffey JG, Kavanagh BP (2017) Insight into acute respiratory distress syndrome: from models to patients. Am J Respir Crit Care Med 196:18–28. 10.1164/rccm.201612-2415CI28146637 10.1164/rccm.201612-2415CI
10. Borges JB EL Costa V Suarez-Sipmann F Widström C Larsson A Amato M Hedenstierna G Early inflammation mainly affects normally and poorly aerated lung in experimental ventilator-induced lung injury Crit Care Med 2014 42 e279 e287 10.1097/CCM.0000000000000161 24448197
Borges JB, EL Costa V, Suarez-Sipmann F, Widström C, Larsson A, Amato M, Hedenstierna G (2014) Early inflammation mainly affects normally and poorly aerated lung in experimental ventilator-induced lung injury. Crit Care Med 42:e279–e287. 10.1097/CCM.000000000000016124448197 10.1097/CCM.0000000000000161
11. Borges JB Costa ELV Bergquist M Lucchetta L Widström C Maripuu E Suarez-Sipmann F Larsson A Amato MBP Hedenstierna G Lung inflammation persists after 27 hours of protective acute respiratory distress syndrome network strategy and is concentrated in the nondependent lung Crit Care Med 2015 43 e123 e132 10.1097/CCM.0000000000000926 25746507
Borges JB, Costa ELV, Bergquist M, Lucchetta L, Widström C, Maripuu E, Suarez-Sipmann F, Larsson A, Amato MBP, Hedenstierna G (2015) Lung inflammation persists after 27 hours of protective acute respiratory distress syndrome network strategy and is concentrated in the nondependent lung. Crit Care Med 43:e123–e132. 10.1097/CCM.000000000000092625746507 10.1097/CCM.0000000000000926
12. Cereda M Xin Y Goffi A Herrmann J Kaczka DW Kavanagh BP Perchiazzi G Yoshida T Rizi RR Imaging the injured lung: mechanisms of action and clinical use Anesthesiology 2019 131 716 749 10.1097/ALN.0000000000002583 30664057
Cereda M, Xin Y, Goffi A, Herrmann J, Kaczka DW, Kavanagh BP, Perchiazzi G, Yoshida T, Rizi RR (2019) Imaging the injured lung: mechanisms of action and clinical use. Anesthesiology 131:716–749. 10.1097/ALN.000000000000258330664057 10.1097/ALN.0000000000002583
13. Schroeder T Vidal Melo MF Musch G Harris RS Venegas JG Winkler T Modeling pulmonary kinetics of 2-deoxy-2-[18F]fluoro-d-glucose during acute lung injury Acad Radiol 2008 15 763 775 10.1016/j.acra.2007.12.016 18486012
Schroeder T, Vidal Melo MF, Musch G, Harris RS, Venegas JG, Winkler T (2008) Modeling pulmonary kinetics of 2-deoxy-2-[18F]fluoro-d-glucose during acute lung injury. Acad Radiol 15:763–775. 10.1016/j.acra.2007.12.01618486012 10.1016/j.acra.2007.12.016
14. Dhelft F Lancelot S Mouton W Le Bars D Costes N Roux E Orkisz M Benzerdjeb N Richard JC Bitker L Prone position decreases acute lung inflammation measured by [11C](R)-PK11195 positron emission tomography in experimental acute respiratory distress syndrome J Appl Physiol 2023 134 467 481 10.1152/japplphysiol.00234.2022 36633865
Dhelft F, Lancelot S, Mouton W, Le Bars D, Costes N, Roux E, Orkisz M, Benzerdjeb N, Richard JC, Bitker L (2023) Prone position decreases acute lung inflammation measured by [11C](R)-PK11195 positron emission tomography in experimental acute respiratory distress syndrome. J Appl Physiol 134:467–481. 10.1152/japplphysiol.00234.202236633865 10.1152/japplphysiol.00234.2022
15. Bitker L Costes N Le Bars D Lavenne F Orkisz M Hoyos MH Benzerdjeb N Devouassoux M Richard JC Noninvasive quantification of macrophagic lung recruitment during experimental ventilation-induced lung injury J Appl Physiol 2019 127 546 558 10.1152/japplphysiol.00825.2018 31169472
Bitker L, Costes N, Le Bars D, Lavenne F, Orkisz M, Hoyos MH, Benzerdjeb N, Devouassoux M, Richard JC (2019) Noninvasive quantification of macrophagic lung recruitment during experimental ventilation-induced lung injury. J Appl Physiol 127:546–558. 10.1152/japplphysiol.00825.201831169472 10.1152/japplphysiol.00825.2018
16. Sica A Mantovani A Macrophage plasticity and polarization: in vivo veritas J Clin Invest 2012 122 787 795 10.1172/JCI59643 22378047
Sica A, Mantovani A (2012) Macrophage plasticity and polarization: in vivo veritas. J Clin Invest 122:787–795. 10.1172/JCI5964322378047 10.1172/JCI59643
17. Dalm VASH Van Hagen PM Van Koetsveld PM Achilefu S Houtsmuller AB Pols DHJ Van Der Lely AJ Lamberts SWJ Hofland LJ Expression of somatostatin, cortistatin, and somatostatin receptors in human monocytes, macrophages, and dendritic cells Am J Physiol Endocrinol Metab 2003 285 344 353 10.1152/ajpendo.00048.2003
Dalm VASH, Van Hagen PM, Van Koetsveld PM, Achilefu S, Houtsmuller AB, Pols DHJ, Van Der Lely AJ, Lamberts SWJ, Hofland LJ (2003) Expression of somatostatin, cortistatin, and somatostatin receptors in human monocytes, macrophages, and dendritic cells. Am J Physiol Endocrinol Metab 285:344–353. 10.1152/ajpendo.00048.200310.1152/ajpendo.00048.2003
18. Reubi JC Schär JC Waser B Wenger S Heppeler A Schmitt JS Mäcke HR Affinity profiles for human somatostatin receptor subtypes SST1–SST5 of somatostatin radiotracers selected for scintigraphic and radiotherapeutic use Eur J Nucl Med 2000 27 273 282 10.1007/s002590050034 10774879
Reubi JC, Schär JC, Waser B, Wenger S, Heppeler A, Schmitt JS, Mäcke HR (2000) Affinity profiles for human somatostatin receptor subtypes SST1–SST5 of somatostatin radiotracers selected for scintigraphic and radiotherapeutic use. Eur J Nucl Med 27:273–282. 10.1007/s00259005003410774879 10.1007/s002590050034
19. Puuvuori E Liggieri F Velikyan I Chiodaroli E Sigfridsson J Romelin H Ingvast S Korsgren O Hulsart-Billström G Perchiazzi G Eriksson O PET–CT imaging of pulmonary inflammation using [68Ga]Ga-DOTA-TATE EJNMMI Res 2022 12 19 10.1186/s13550-022-00892-0 35394238
Puuvuori E, Liggieri F, Velikyan I, Chiodaroli E, Sigfridsson J, Romelin H, Ingvast S, Korsgren O, Hulsart-Billström G, Perchiazzi G, Eriksson O (2022) PET–CT imaging of pulmonary inflammation using [68Ga]Ga-DOTA-TATE. EJNMMI Res 12:19. 10.1186/s13550-022-00892-035394238 10.1186/s13550-022-00892-0
20. Percie du Sert N Hurst V Ahluwalia A Alam S Avey MT Baker M Browne WJ Clark A Cuthill IC Dirnagl U Emerson M Garner P Holgate ST Howells DW Karp NA Lazic SE Lidster K MacCallum CJ Macleod M Pearl EJ Petersen OH Rawle F Reynolds P Rooney K Sena ES Silberberg SD Steckler T Würbel H The ARRIVE guidelines 2.0: updated guidelines for reporting animal research PLoS Biol 2020 18 e3000410 10.1371/journal.pbio.3000410 32663219
Percie du Sert N, Hurst V, Ahluwalia A, Alam S, Avey MT, Baker M, Browne WJ, Clark A, Cuthill IC, Dirnagl U, Emerson M, Garner P, Holgate ST, Howells DW, Karp NA, Lazic SE, Lidster K, MacCallum CJ, Macleod M, Pearl EJ, Petersen OH, Rawle F, Reynolds P, Rooney K, Sena ES, Silberberg SD, Steckler T, Würbel H (2020) The ARRIVE guidelines 2.0: updated guidelines for reporting animal research. PLoS Biol 18:e3000410. 10.1371/journal.pbio.300041032663219 10.1371/journal.pbio.3000410
21. The Acute Respiratory Distress Syndrome Network Ventilation with lower tidal volumes as compared with traditional tidal volumes for acute lung injury and the acute respiratory distress syndrome N Engl J Med 2000 342 1301 1308 10.1056/NEJM200005043421801 10793162
The Acute Respiratory Distress Syndrome Network (2000) Ventilation with lower tidal volumes as compared with traditional tidal volumes for acute lung injury and the acute respiratory distress syndrome. N Engl J Med 342:1301–1308. 10.1056/NEJM20000504342180110793162 10.1056/NEJM200005043421801
22. Perchiazzi G Rylander C Derosa S Pellegrini M Pitagora L Polieri D Vena A Tannoia A Fiore T Hedenstierna G Regional distribution of lung compliance by image analysis of computed tomograms Respir Physiol Neurobiol 2014 10.1016/j.resp.2014.07.001 25026158
Perchiazzi G, Rylander C, Derosa S, Pellegrini M, Pitagora L, Polieri D, Vena A, Tannoia A, Fiore T, Hedenstierna G (2014) Regional distribution of lung compliance by image analysis of computed tomograms. Respir Physiol Neurobiol. 10.1016/j.resp.2014.07.00125026158 10.1016/j.resp.2014.07.001
23. Gu X Pan H Liang Y Castillo R Yang D Choi D Castillo E Majumdar A Guerrero T Jiang SB Implementation and evaluation of various demons deformable image registration algorithms on a GPU Phys Med Biol 2010 55 207 219 10.1088/0031-9155/55/1/012 20009197
Gu X, Pan H, Liang Y, Castillo R, Yang D, Choi D, Castillo E, Majumdar A, Guerrero T, Jiang SB (2010) Implementation and evaluation of various demons deformable image registration algorithms on a GPU. Phys Med Biol 55:207–219. 10.1088/0031-9155/55/1/01220009197 10.1088/0031-9155/55/1/012
24. Pellegrini M Hedenstierna G Larsson AS Perchiazzi G Inspiratory efforts, positive end-expiratory pressure, and external resistances influence intraparenchymal gas redistribution in mechanically ventilated injured lungs Front Physiol 2021 11 1 11 10.3389/fphys.2020.618640
Pellegrini M, Hedenstierna G, Larsson AS, Perchiazzi G (2021) Inspiratory efforts, positive end-expiratory pressure, and external resistances influence intraparenchymal gas redistribution in mechanically ventilated injured lungs. Front Physiol 11:1–11. 10.3389/fphys.2020.61864010.3389/fphys.2020.618640
25. Gattinoni L Pesenti A Bombino M Baglioni S Rivolta M Rossi F Rossi G Fumagalli R Marcolin R Mascheroni D Relationships between lung computed tomographic density, gas exchange, and PEEP in acute respiratory failure Anesthesiology 1988 69 824 832 10.1097/00000542-198812000-00005 3057937
Gattinoni L, Pesenti A, Bombino M, Baglioni S, Rivolta M, Rossi F, Rossi G, Fumagalli R, Marcolin R, Mascheroni D (1988) Relationships between lung computed tomographic density, gas exchange, and PEEP in acute respiratory failure. Anesthesiology 69:824–832. 10.1097/00000542-198812000-000053057937 10.1097/00000542-198812000-00005
26. 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 15s 25s 10.1183/09031936.03.00021303
Gattinoni L, Carlesso E, Cadringher P, Valenza F, Vagginelli F, Chiumello D (2003) Physical and biological triggers of ventilator-induced lung injury and its prevention. Eur Respir J 22:15s–25s10.1183/09031936.03.00021303
27. Megson THG Megson THG Chapter 7—Stress and strain Structural and stress analysis 2019 4 Oxford Butterworth-Heinemann 165 201
Megson THG (2019) Chapter 7—Stress and strain. In: Megson THG (ed) Structural and stress analysis, 4th edn. Butterworth-Heinemann, Oxford, pp 165–201
28. Graham MM Ross BD Gambhir SS Quantification of radiotracer uptake into tissue Molecular imaging: principles and practice 2021 2 London Academic Press 1613 1624
Graham MM (2021) Quantification of radiotracer uptake into tissue. In: Ross BD, Gambhir SS (eds) Molecular imaging: principles and practice, 2nd edn. Academic Press, London, pp 1613–1624
29. Motulsky H Christopoulos A Comparing models using the extra sum-of-squares F test Fitting models to biological data using linear and nonlinear regression 2003 San Diego GraphPad Software Inc. 138 142
Motulsky H, Christopoulos A (2003) Comparing models using the extra sum-of-squares F test. Fitting models to biological data using linear and nonlinear regression. GraphPad Software Inc., San Diego, pp 138–142
30. Pelosi P D’Andrea L Vitale G Pesenti A Gattinoni L Vertical gradient of regional lung inflation in adult respiratory distress syndrome Am J Respir Crit Care Med 1994 149 8 13 10.1164/ajrccm.149.1.8111603 8111603
Pelosi P, D’Andrea L, Vitale G, Pesenti A, Gattinoni L (1994) Vertical gradient of regional lung inflation in adult respiratory distress syndrome. Am J Respir Crit Care Med 149:8–138111603 10.1164/ajrccm.149.1.8111603
31. Paternoster R Brame R Mazerolle P Piquero A Using the correct statistical test for the equality of regression coefficients Criminology 1998 36 859 866 10.1111/j.1745-9125.1998.tb01268.x
Paternoster R, Brame R, Mazerolle P, Piquero A (1998) Using the correct statistical test for the equality of regression coefficients. Criminology 36:859–866. 10.1111/j.1745-9125.1998.tb01268.x10.1111/j.1745-9125.1998.tb01268.x
32. Lange OF Grubmüller H Generalized correlation for biomolecular dynamics Proteins Struct Funct Genet 2006 62 1053 1061 10.1002/prot.20784 16355416
Lange OF, Grubmüller H (2006) Generalized correlation for biomolecular dynamics. Proteins Struct Funct Genet 62:1053–1061. 10.1002/prot.2078416355416 10.1002/prot.20784
33. Mead J Takishima T Leith D Stress distribution in lungs: a model of pulmonary elasticity J Appl Physiol 1970 28 596 608 10.1152/jappl.1970.28.5.596 5442255
Mead J, Takishima T, Leith D (1970) Stress distribution in lungs: a model of pulmonary elasticity. J Appl Physiol 28:596–6085442255 10.1152/jappl.1970.28.5.596
34. Milic-Emili J Macklem PT Mead J Static distribution of lung volumes Handbook of physiology 1986 Bethesda American Physiological Society 561 574
Milic-Emili J (1986) Static distribution of lung volumes. In: Macklem PT, Mead J (eds) Handbook of physiology, vol 3. Mechanics of breathing (section 3). American Physiological Society, Bethesda, pp 561–574
35. Evren E Ringqvist E Willinger T Origin and ontogeny of lung macrophages: from mice to humans Immunology 2020 160 126 138 10.1111/imm.13154 31715003
Evren E, Ringqvist E, Willinger T (2020) Origin and ontogeny of lung macrophages: from mice to humans. Immunology 160:126–138. 10.1111/imm.1315431715003 10.1111/imm.13154
36. West JB Dollery CT Distribution of blood flow and ventilation-perfusion ratio in the lung, measured with radioactive carbon dioxide J Appl Physiol 1960 15 405 410 10.1152/jappl.1960.15.3.405 13844133
West JB, Dollery CT (1960) Distribution of blood flow and ventilation-perfusion ratio in the lung, measured with radioactive carbon dioxide. J Appl Physiol 15:405–41013844133 10.1152/jappl.1960.15.3.405
37. Galvin I Drummond GB Nirmalan M Distribution of blood flow and ventilation in the lung: gravity is not the only factor Br J Anaesth 2007 98 420 428 10.1093/bja/aem036 17347182
Galvin I, Drummond GB, Nirmalan M (2007) Distribution of blood flow and ventilation in the lung: gravity is not the only factor. Br J Anaesth 98:420–428. 10.1093/bja/aem03617347182 10.1093/bja/aem036
38. Goshtasby AA Similarity and dissimilarity measures Image registration: principles, tools and methods 2012 London Springer 7 66
Goshtasby AA (2012) Similarity and dissimilarity measures. Image registration: principles, tools and methods. Springer, London, pp 7–66
39. Shannon CE A mathematical theory of communication Bell Syst Tech J 1948 27 379 423 10.1002/j.1538-7305.1948.tb01338.x
Shannon CE (1948) A mathematical theory of communication. Bell Syst Tech J 27:379–423. 10.1002/j.1538-7305.1948.tb01338.x10.1002/j.1538-7305.1948.tb01338.x
40. Gattinoni L Carlesso E Caironi P Stress and strain within the lung Curr Opin Crit Care 2012 18 42 47 10.1097/MCC.0b013e32834f17d9 22157254
Gattinoni L, Carlesso E, Caironi P (2012) Stress and strain within the lung. Curr Opin Crit Care 18:42–4722157254 10.1097/MCC.0b013e32834f17d9
41. Pelosi P Crotti S Brazzi L Gattinoni L Computed tomography in adult respiratory distress syndrome: what has it taught us? Eur Respir J 1996 9 1055 1062 10.1183/09031936.96.09051055 8793469
Pelosi P, Crotti S, Brazzi L, Gattinoni L (1996) Computed tomography in adult respiratory distress syndrome: what has it taught us? Eur Respir J 9:1055–1062. 10.1183/09031936.96.090510558793469 10.1183/09031936.96.09051055
42. Simon BA Computed tomography studies of lung mechanics Proc Am Thorac Soc 2005 2 517 521 10.1513/pats.200507-076DS 16352757
Simon BA (2005) Computed tomography studies of lung mechanics. Proc Am Thorac Soc 2:517–521. 10.1513/pats.200507-076DS16352757 10.1513/pats.200507-076DS
43. Faffe DS Zin WA Lung parenchymal mechanics in health and disease Physiol Rev 2009 89 759 775 10.1152/physrev.00019.2007 19584312
Faffe DS, Zin WA (2009) Lung parenchymal mechanics in health and disease. Physiol Rev 89:759–775. 10.1152/physrev.00019.200719584312 10.1152/physrev.00019.2007
44. Ritter MC Jesudason R Majumdar A Stamenović D Buczek-Thomas JA Stone PJ Nugent MA Suki B A zipper network model of the failure mechanics of extracellular matrices Proc Natl Acad Sci USA 2009 106 1081 1086 10.1073/pnas.0808414106 19144920
Ritter MC, Jesudason R, Majumdar A, Stamenović D, Buczek-Thomas JA, Stone PJ, Nugent MA, Suki B (2009) A zipper network model of the failure mechanics of extracellular matrices. Proc Natl Acad Sci USA 106:1081–1086. 10.1073/pnas.080841410619144920 10.1073/pnas.0808414106
45. Lum H Mitzner W A species comparison of alveolar size and surface forces J Appl Physiol 1987 62 1865 1871 10.1152/jappl.1987.62.5.1865 3597260
Lum H, Mitzner W (1987) A species comparison of alveolar size and surface forces. J Appl Physiol 62:1865–1871. 10.1152/jappl.1987.62.5.18653597260 10.1152/jappl.1987.62.5.1865
46. Bates JHT A recruitment model of quasi-linear power-low stress adaptation in lung tissue Ann Biomed Eng 2007 35 1165 1174 10.1007/s10439-007-9291-0 17380389
Bates JHT (2007) A recruitment model of quasi-linear power-low stress adaptation in lung tissue. Ann Biomed Eng 35:1165–117417380389 10.1007/s10439-007-9291-0
47. 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 990 1001 10.1152/jappl.1979.47.5.990 511725
Gil J, Bachofen H, Gehr P, Weibel ER (1979) Alveolar volume-surface area relation in air- and saline-filled lungs fixed by vascular perfusion. J Appl Physiol Respir Environ Exerc Physiol 47:990–1001. 10.1152/jappl.1979.47.5.990511725 10.1152/jappl.1979.47.5.990
48. Knudsen L Ochs M The micromechanics of lung alveoli: structure and function of surfactant and tissue components Histochem Cell Biol 2018 150 661 676 10.1007/s00418-018-1747-9 30390118
Knudsen L, Ochs M (2018) The micromechanics of lung alveoli: structure and function of surfactant and tissue components. Histochem Cell Biol 150:661–676. 10.1007/s00418-018-1747-930390118 10.1007/s00418-018-1747-9
49. Li X Samnick S Lapa C Israel I Buck AK Kreissl MC Bauer W 68Ga-DOTATATE PET/CT for the detection of inflammation of large arteries: correlation with 18F-FDG, calcium burden and risk factors EJNMMI Res 2012 2 1 10 10.1186/2191-219X-2-52 22251281
Li X, Samnick S, Lapa C, Israel I, Buck AK, Kreissl MC, Bauer W (2012) 68Ga-DOTATATE PET/CT for the detection of inflammation of large arteries: correlation with 18F-FDG, calcium burden and risk factors. EJNMMI Res 2:1–10. 10.1186/2191-219X-2-5222251281 10.1186/2191-219X-2-52
50. Rominger A Saam T Vogl E Übleis C La FC Förster S Haug A Cumming P Reiser MF Nikolaou K Bartenstein P Hacker M In vivo imaging of macrophage activity in the coronary arteries using 68Ga-DOTATATE PET/CT: correlation with coronary calcium burden and risk factors J Nucl Med 2010 51 193 197 10.2967/jnumed.109.070672 20080898
Rominger A, Saam T, Vogl E, Übleis C, La FC, Förster S, Haug A, Cumming P, Reiser MF, Nikolaou K, Bartenstein P, Hacker M (2010) In vivo imaging of macrophage activity in the coronary arteries using 68Ga-DOTATATE PET/CT: correlation with coronary calcium burden and risk factors. J Nucl Med 51:193–197. 10.2967/jnumed.109.07067220080898 10.2967/jnumed.109.070672
51. Tarkin JM Joshi FR Evans NR Chowdhury MM Figg NL Shah AV Starks LT Martin-Garrido A Manavaki R Yu E Kuc RE Grassi L Kreuzhuber R Kostadima MA Frontini M Kirkpatrick PJ Coughlin PA Gopalan D Fryer TD Buscombe JR Groves AM Ouwehand WH Bennett MR Warburton EA Davenport AP Rudd JHF Detection of atherosclerotic inflammation by 68Ga-DOTATATE PET compared to [18F]FDG PET imaging J Am Coll Cardiol 2017 69 1774 1791 10.1016/j.jacc.2017.01.060 28385306
Tarkin JM, Joshi FR, Evans NR, Chowdhury MM, Figg NL, Shah AV, Starks LT, Martin-Garrido A, Manavaki R, Yu E, Kuc RE, Grassi L, Kreuzhuber R, Kostadima MA, Frontini M, Kirkpatrick PJ, Coughlin PA, Gopalan D, Fryer TD, Buscombe JR, Groves AM, Ouwehand WH, Bennett MR, Warburton EA, Davenport AP, Rudd JHF (2017) Detection of atherosclerotic inflammation by 68Ga-DOTATATE PET compared to [18F]FDG PET imaging. J Am Coll Cardiol 69:1774–1791. 10.1016/j.jacc.2017.01.06028385306 10.1016/j.jacc.2017.01.060
52. Huang X Xiu H Zhang S Zhang G The role of macrophages in the pathogenesis of ALI/ARDS Mediat Inflamm 2018 2018 1 8 10.1155/2018/1264913
Huang X, Xiu H, Zhang S, Zhang G (2018) The role of macrophages in the pathogenesis of ALI/ARDS. Mediat Inflamm 2018:1–8. 10.1155/2018/126491310.1155/2018/1264913
53. Frank JA Wray CM Mcauley DF Schwendener R Matthay MA Alveolar macrophages contribute to alveolar barrier dysfunction in ventilator-induced lung injury Am J Physiol Lung Cell Mol Physiol 2023 10.1152/ajplung.00055.2006 37366533
Frank JA, Wray CM, Mcauley DF, Schwendener R, Matthay MA (2023) Alveolar macrophages contribute to alveolar barrier dysfunction in ventilator-induced lung injury. Am J Physiol Lung Cell Mol Physiol. 10.1152/ajplung.00055.200637366533 10.1152/ajplung.00055.2006
54. Johnston LK Rims CR Gill SE McGuire JK Manicone AM Pulmonary macrophage subpopulations in the induction and resolution of acute lung injury Am J Respir Cell Mol Biol 2012 47 417 426 10.1165/rcmb.2012-0090OC 22721830
Johnston LK, Rims CR, Gill SE, McGuire JK, Manicone AM (2012) Pulmonary macrophage subpopulations in the induction and resolution of acute lung injury. Am J Respir Cell Mol Biol 47:417–426. 10.1165/rcmb.2012-0090OC22721830 10.1165/rcmb.2012-0090OC
55. Herold S Mayer K Lohmeyer J Acute lung injury: how macrophages orchestrate resolution of inflammation and tissue repair Front Immunol 2011 2 1 13 10.3389/fimmu.2011.00065 22566792
Herold S, Mayer K, Lohmeyer J (2011) Acute lung injury: how macrophages orchestrate resolution of inflammation and tissue repair. Front Immunol 2:1–13. 10.3389/fimmu.2011.0006522566792 10.3389/fimmu.2011.00065
56. Joseph P Tawakol A Imaging atherosclerosis with positron emission tomography Eur Heart J 2016 37 2974 2980b 10.1093/eurheartj/ehw147 27125951
Joseph P, Tawakol A (2016) Imaging atherosclerosis with positron emission tomography. Eur Heart J 37:2974–2980b. 10.1093/eurheartj/ehw14727125951 10.1093/eurheartj/ehw147
57. Dahabreh IJ Gatsonis C A flexible, multifaceted approach is needed in health technology assessment of PET J Nucl Med 2014 55 1225 1227 10.2967/jnumed.114.142331 25047328
Dahabreh IJ, Gatsonis C (2014) A flexible, multifaceted approach is needed in health technology assessment of PET. J Nucl Med 55:1225–1227. 10.2967/jnumed.114.14233125047328 10.2967/jnumed.114.142331
58. Hatori A Yui J Yamasaki T Xie L Kumata K Fujinaga M Yoshida Y Ogawa M Nengaki N Kawamura K Fukumura T Zhang MR PET imaging of lung inflammation with [18F]FEDAC, a radioligand for translocator protein (18 kDa) PLoS ONE 2012 10.1371/journal.pone.0045065 22984611
Hatori A, Yui J, Yamasaki T, Xie L, Kumata K, Fujinaga M, Yoshida Y, Ogawa M, Nengaki N, Kawamura K, Fukumura T, Zhang MR (2012) PET imaging of lung inflammation with [18F]FEDAC, a radioligand for translocator protein (18 kDa). PLoS ONE. 10.1371/journal.pone.004506522984611 10.1371/journal.pone.0045065
59. Narayan N Mandhair H Smyth E Dakin SG Kiriakidis S Wells L Owen D Sabokbar A Taylor P The macrophage marker translocator protein (TSPO) is down-regulated on pro-inflammatory ‘M1’ human macrophages PLoS ONE 2017 12 1 19 10.1371/journal.pone.0185767
Narayan N, Mandhair H, Smyth E, Dakin SG, Kiriakidis S, Wells L, Owen D, Sabokbar A, Taylor P (2017) The macrophage marker translocator protein (TSPO) is down-regulated on pro-inflammatory ‘M1’ human macrophages. PLoS ONE 12:1–19. 10.1371/journal.pone.018576710.1371/journal.pone.0185767
60. Matute-Bello G Frevert CW Martin TR Animal models of acute lung injury Am J Physiol Cell Mol Physiol 2008 295 L379 L399 10.1152/ajplung.00010.2008
Matute-Bello G, Frevert CW, Martin TR (2008) Animal models of acute lung injury. Am J Physiol Cell Mol Physiol 295:L379–L399. 10.1152/ajplung.00010.200810.1152/ajplung.00010.2008
61. Broche L Perchiazzi G Porra L Tannoia A Pellegrini M Derosa S Sindaco A Batista Borges J Degrugilliers L Larsson A Hedenstierna G Wexler AS Bravin A Verbanck S Smith BJ Bates JHT Bayat S Dynamic mechanical interactions between neighboring airspaces determine cyclic opening and closure in injured lung Crit Care Med 2017 45 687 694 10.1097/CCM.0000000000002234 28107207
Broche L, Perchiazzi G, Porra L, Tannoia A, Pellegrini M, Derosa S, Sindaco A, Batista Borges J, Degrugilliers L, Larsson A, Hedenstierna G, Wexler AS, Bravin A, Verbanck S, Smith BJ, Bates JHT, Bayat S (2017) Dynamic mechanical interactions between neighboring airspaces determine cyclic opening and closure in injured lung. Crit Care Med 45:687–694. 10.1097/CCM.000000000000223428107207 10.1097/CCM.0000000000002234
62. Caironi P Cressoni M Chiumello D Ranieri M Quintel M Russo SG Cornejo R Bugedo G Carlesso E Russo R Caspani L Gattinoni L Lung opening and closing during ventilation of acute respiratory distress syndrome Am J Respir Crit Care Med 2010 181 578 586 10.1164/rccm.200905-0787OC 19910610
Caironi P, Cressoni M, Chiumello D, Ranieri M, Quintel M, Russo SG, Cornejo R, Bugedo G, Carlesso E, Russo R, Caspani L, Gattinoni L (2010) Lung opening and closing during ventilation of acute respiratory distress syndrome. Am J Respir Crit Care Med 181:578–586. 10.1164/rccm.200905-0787OC19910610 10.1164/rccm.200905-0787OC
63. Paula LF Wellman TJ Winkler T Spieth PM Güldner A Venegas JG Gama de Abreu M Carvalho AR Vidal Melo MF Regional tidal lung strain in mechanically ventilated normal lungs J Appl Physiol 2016 121 1335 1347 10.1152/japplphysiol.00861.2015 27445302
Paula LF, Wellman TJ, Winkler T, Spieth PM, Güldner A, Venegas JG, Gama de Abreu M, Carvalho AR, Vidal Melo MF (2016) Regional tidal lung strain in mechanically ventilated normal lungs. J Appl Physiol 121:1335–1347. 10.1152/japplphysiol.00861.201527445302 10.1152/japplphysiol.00861.2015
64. Roan E Waters CM What do we know about mechanical strain in lung alveoli? Am J Physiol Lung Cell Mol Physiol 2011 301 L625 L635 10.1152/ajplung.00105.2011 21873445
Roan E, Waters CM (2011) What do we know about mechanical strain in lung alveoli? Am J Physiol Lung Cell Mol Physiol 301:L625–L635. 10.1152/ajplung.00105.201121873445 10.1152/ajplung.00105.2011
65. Protti A Votta E Gattinoni L Which is the most important strain in the pathogenesis of ventilator-induced lung injury: dynamic or static? Curr Opin Crit Care 2014 20 33 38 10.1097/MCC.0000000000000047 24247615
Protti A, Votta E, Gattinoni L (2014) Which is the most important strain in the pathogenesis of ventilator-induced lung injury: dynamic or static? Curr Opin Crit Care 20:33–38. 10.1097/MCC.000000000000004724247615 10.1097/MCC.0000000000000047
66. Jahandideh A Uotila S Ståhle M Virta J Li XG Kytö V Marjamäki P Liljenbäck H Taimen P Oikonen V Lehtonen J Mäyränpää MI Chen Q Low PS Knuuti J Roivainen A Saraste A Folate receptor b–targeted PET imaging of macrophages in autoimmune myocarditis J Nucl Med 2020 61 1643 1649 10.2967/jnumed.119.241356 32284397
Jahandideh A, Uotila S, Ståhle M, Virta J, Li XG, Kytö V, Marjamäki P, Liljenbäck H, Taimen P, Oikonen V, Lehtonen J, Mäyränpää MI, Chen Q, Low PS, Knuuti J, Roivainen A, Saraste A (2020) Folate receptor b–targeted PET imaging of macrophages in autoimmune myocarditis. J Nucl Med 61:1643–1649. 10.2967/jnumed.119.24135632284397 10.2967/jnumed.119.241356
67. Verweij NJF Yaqub M Bruijnen STG Pieplenbosch S ter Wee MM Jansen G Chen Q Low PS Windhorst AD Lammertsma AA Hoekstra OS Voskuyl AE van der Laken CJ First in man study of [18F]fluoro-PEG-folate PET: a novel macrophage imaging technique to visualize rheumatoid arthritis Sci Rep 2020 10 1 10 10.1038/s41598-020-57841-x 31913322
Verweij NJF, Yaqub M, Bruijnen STG, Pieplenbosch S, ter Wee MM, Jansen G, Chen Q, Low PS, Windhorst AD, Lammertsma AA, Hoekstra OS, Voskuyl AE, van der Laken CJ (2020) First in man study of [18F]fluoro-PEG-folate PET: a novel macrophage imaging technique to visualize rheumatoid arthritis. Sci Rep 10:1–10. 10.1038/s41598-020-57841-x31913322 10.1038/s41598-020-57841-x
68. Russakoff DB Tomasi C Rohlfing T Maurer CR Pajdla T Matas J Image similarity using mutual information of regions Computer vision—ECCV 2004 2004 Berlin Springer 596 607
Russakoff DB, Tomasi C, Rohlfing T, Maurer CR (2004) Image similarity using mutual information of regions. In: Pajdla T, Matas J (eds) Computer vision—ECCV 2004. Springer, Berlin, pp 596–607
69. Pluim JPW Maintz JBA Viergever MA Mutual-information-based registration of medical images: a survey IEEE Trans Med Imaging 2003 22 986 1004 10.1109/TMI.2003.815867 12906253
Pluim JPW, Maintz JBA, Viergever MA (2003) Mutual-information-based registration of medical images: a survey. IEEE Trans Med Imaging 22:986–1004. 10.1109/TMI.2003.81586712906253 10.1109/TMI.2003.815867
70. Maes F Collignon A Vandermeulen D Marchal G Suetens P Multimodality image registration by maximization of mutual information IEEE Trans Med Imaging 1997 16 187 198 10.1109/42.563664 9101328
Maes F, Collignon A, Vandermeulen D, Marchal G, Suetens P (1997) Multimodality image registration by maximization of mutual information. IEEE Trans Med Imaging 16:187–198. 10.1109/42.5636649101328 10.1109/42.563664
71. van Leer B van Rijsewijk ND Nijsten MWN Slart RHJA Pillay J Glaudemans AWJM Practice of 18F-FDG–PET/CT in ICU patients: a systematic review Semin Nucl Med 2023 53 809 819 10.1053/j.semnuclmed.2023.05.003 37258380
van Leer B, van Rijsewijk ND, Nijsten MWN, Slart RHJA, Pillay J, Glaudemans AWJM (2023) Practice of 18F-FDG–PET/CT in ICU patients: a systematic review. Semin Nucl Med 53:809–819. 10.1053/j.semnuclmed.2023.05.00337258380 10.1053/j.semnuclmed.2023.05.003
