- Dipartimento di Medicina e Chirurgia, Università Milano-Bicocca, Monza, Italy
The architecture of the air-blood barrier is effective in optimizing the gas exchange as long as it retains its specific feature of extreme thinness reflecting, in turn, a strict control on the extravascular water to be kept at minimum. Edemagenic conditions may perturb this equilibrium by increasing microvascular filtration; this characteristically occurs when cardiac output increases to balance the oxygen uptake with the oxygen requirement such as in exercise and hypoxia (either due to low ambient pressure or reflecting a pathological condition). In general, the lung is well equipped to counteract an increase in microvascular filtration rate. The loss of control on fluid balance is the consequence of disruption of the integrity of the macromolecular structure of lung tissue. This review, merging data from experimental approaches and evidence in humans, will explore how the heterogeneity in morphology, mechanical features and perfusion of the terminal respiratory units might impact on lung fluid balance and its control. Evidence is also provided that heterogeneities may be inborn and they could actually get worse as a consequence of a developing pathological process. Further, data are presented how in humans inter-individual heterogeneities in morphology of the terminal respiratory hinder the control of fluid balance and, in turn, hamper the efficiency of the oxygen diffusion-transport function.
Introduction
In general, the architecture of the air-blood barrier is considered as being very effective in optimizing the gas exchange (Swan and Tawhai, 2011) by preserving the specific feature to assure it, namely its extreme thinness, that, in turn, reflects a strict control to keep lung extravascular water at minimum (Miserocchi, 2009). This occurs in physiological conditions despite the existence of morphological and mechanical heterogeneities of the terminal lung units. This review will explore how the heterogeneity in the morphology and perfusion of the terminal respiratory units come into play when the lung is facing edemagenic conditions that notably impact on the control of extravascular water and, in parallel, on the efficiency of the oxygen diffusion-transport at the level of the alveolar-capillary unit (Miserocchi et al., 2022). Further, evidence is provided concerning inter-individual differences in humans justifying the proneness to develop lung edema.
The structure-function of the air-blood barrier
The air–blood barrier represents the interface for gas exchanges; it is made of the capillary endothelium, an interstitial compartment and the alveolar epithelium. The surface area of the air-blood barrier is of the order of 2000 cm2/g lung tissue. Its thin portion, accounting for almost 50% of the whole surface, is <0.5 μm thick (Figure 1A) and is simply made of endothelium, epithelium and an intervening fused basement membrane representing the interstitial compartment. The overall surface of the thin portion is estimated at ∼100 m2 in the human lung (Weibel, 2017).
FIGURE 1. (A) Microphotograph of the thin portion of the air-blood barrier separating the capillary blood (CAP) from the alveolar air (ALV); endothelial cell (EN), epithelial cell (EP) and basement membrane (BM) representing the interstitial compartment (from Botto et al., 2006). (B) Lower magnitude TEM image to include the thick and the thin portion of the air-blood barrier (abb); IC are the interstitial cells, ECM is extracellular matrix (from Conforti et al., 2002).
A valid representative significance of the structure-function of the air-blood barrier rests on the comparison between oxygen diffusion, that is in the range of 15 × 10−2 mL/(min cm2) at rest, while capillary microvascular filtration is at least 10,000 fold less (Miserocchi, 2009); thus, the air-blood barrier allows easy transfer of respiratory gases but it is minimally permeable to water. The diffusional oxygen uptake in the air-blood barrier (
being
The fibrillar molecular component is mostly present in the “thick portion” of the air blood barrier (Figure 1B) and includes collagen I and III and elastic fibers (Conforti et al., 2002). The other important extracellular matrix component includes the proteoglycans (PGs). PGs consist of a “core protein” with one or more covalently attached glycosaminoglycan (GAG) chains. GAGs are linear polysaccharides made of repeating disaccharide units; they are highly polar and attract water at their hydrophilic end. GAGs are widely present in the body acting as lubricants, in particular when pressures are being generated among adjacent structures. PGs family comprises the large chondroitin sulphate sub-family (>1,000 kDa) that fills the voids between the fibrillar component. PGs act as link proteins binding with several other molecules, as well as with cell surface. The low-energy ionic and/or non-covalent bonds of chondroitin sulphate PGs with the collagen-elastin network allow reciprocal movement/sliding by reducing shear stresses during movements but also assure mechanical stability of the lung parenchyma (Cavalcante et al., 2005; Ma and Bates, 2014). In the thin portion of the basement membrane the PG-heparan sulphate sub-family (300–500 kDa) play a key role in the control on extravascular lung water assuring low microvascular permeability of the capillary surface to water and proteins as well as low porosity of the interstitial matrix (Negrini et al., 1996).
The molecular architecture of the alveolar compartment and the presence of surfactant on the alveolar surface represent stabilizing factors that redistribute forces between adjacent alveoli, so as to avoid alveolar collapse or overdistension in physiological conditions (Knudsen and Ochs, 2018).
The Starling law governing fluid exchanges
The law governing fluid exchanges (
being P and Π the hydraulic the colloidosmotic pressures in the compartments. Factor σ, named proteins reflection coefficient, is the permselectivity of the membrane to plasma proteins. The value of σ varies between 0 and 1, reflecting the molecular sieve of membrane pores relative to the size of the proteins. The equation includes the filtration coefficient Kf = Lp∙A, being Lp the hydraulic conductance and A the surface area available for fluid exchange. It is of interest to recall that the endothelial cells of the microvascular district, when compared to conduit vessels, exhibited a ∼20 fold lower Lp (Parker et al., 2006) and σ approaching 1 (Parker, 2007). Accordingly, concerning fluid and protein exchanges, a functional division has been proposed, between the two compartments placed in series, the “extra-alveolar” and the “true alveolar” compartment (Figure 1A), the latter being specifically designed to shield the lung against edema in order to favour gas exchange. Fluid drainage in the extra-alveolar compartment is drained by an extended lymphatic (Levick and Michel, 2010), while the “true alveolar” compartment, despite its incredible surface, is not served by an extended lymphatic network (Schraufnagel, et al., 2003).
Table 1 reports the estimated trans-endothelial and trans-epithelial Starling pressure gradients based on interstitial and alveolar liquid hydraulic pressure (Pcap, Pi and Pliq alv, respectively) and of oncotic pressure in the same compartments (Πcap, Πint and Πliq alv, respectively). Protein reflection coefficients for endothelium and epithelium (σ endo and σ epi, respectively) are also reported. Alveolar surface tension (γ) is in dyne/cm. In bold are hydraulic (ΔP) and oncotic (σΔΠ) pressure gradients and total Starling pressure gradient. Positive values of the Starling gradient at endothelial level indicate filtration into interstitium; negative value at epithelial level indicate alveolar reabsorption (from Beretta et al., 2021). The water content of the lung is well defined by the wet-to-dry weight ratio (W/D, index of amount of extravascular water) that, in physiological conditions, is ∼5 (Miserocchi et al., 1993). The thinness of the air-blood barrier reflects a strict control of extravascular water volume that is maintained at minimum, thanks to the extremely low water permeability (Lp) and a relatively high value of σ of the endothelial and the epithelial barriers.
TABLE 1. Trans-endothelial and trans-epithelial Starling pressure gradients at end-expiration (FRC) in physiological condition.
Figure 2 presents a conceptual diagram to depict the fluid dynamic and mechanical features involved in the control of lung fluid balance in physiological conditions. The interstitial compliance, defined as the ratio of the change in interstitial volume to the corresponding change in Pi, is also indicated; its low value in physiological conditions, 0.2 mL mmHg-1 100 g-1, hinders the accumulation of interstitial fluid (Miserocchi et al., 2001).
FIGURE 2. Conceptual diagram depicting the fluid dynamic and mechanical features involved in the control of lung fluid balance in physiological conditions.
The intrinsic resistance of the lung to developing edema
Experimental models of edema were developed by slowly inducing a perturbation in fluid balance, this allowed to study the time course of changes in mechanical structure/properties of the interstitial compartment to affect the control on the volume of extravascular water. These kind of informations could not be provided by experimental approaches causing acute severe lesion.
It was found that in the initial phase of developing lung edema, water filtering in the interstitium is captured by the hyaluronan-versican complex to form gel, whose increase in steric hindrance raises the hydraulic interstitial pressure from −10 cmH2O up to ∼ +5 cmH2O, reflecting the low tissue compliance (Miserocchi et al., 1993; Negrini et al., 1996; Miserocchi et al., 2001). The corresponding increase in extravascular lung water, defined as interstitial edema, averaged ∼10%, from 5 to 5.5 W/D. On fluid dynamic basis, the increase in interstitial pressure due to low lung tissue compliance counters filtration, representing therefore a strong “safety factor” (Figures 3A, B). The condition of interstitial edema is likely quite effective, providing a satisfactory answer to the question posed by Hopkins (2010): “given the sporadic nature of exercise-induced edema … why do not more athletes develop alveolar flooding during exercise?”
FIGURE 3. TEM image of clusters of cross and longitudinally sectioned collagen fibrils in the interstitium of the alveolar septa in the less dependent region of a control lung (A) and in the more dependent region after inducing interstitial edema (B) Scale bar 450 nm (from Conforti et al., 2002).
A major question remains: how long can the lung resist in sustained edemagenic conditions?
The kern of severe edema formation is the progressive loss of function of PGs families causing the waning of the “safety factor”.
Figure 3 shows clusters of cross sectioned collagen fibrils in the interstitium of the alveolar septa in the less dependent region of a control lung (Figure 3A) and in the more dependent region after inducing interstitial edema (Figure 3B) (data from Conforti et al., 2002). The white voids surrounding the fibrillary component are occupied by PG-chondroitin sulphate sub-family. Clearly, interstitial edema caused an increase in inter-fibrillar distance of variable extent, roughly doubling from ∼50 to ∼100 nm (Conforti et al., 2002). Intermolecular forces are governed by the inverse-square law so that attraction force falls off with the square of the distance between the interacting molecules. Accordingly, were the estimates relative to interstitial hydration representative of the average increase in intermolecular distance, the force of the non-covalent bond would fall to 1/4 of its original strength.
Besides loosing force of bonds, fragmentation of PGs also occurs down a sequence of events including the production of reactive oxygen species (ROS) due to inflammation (bacterial, viral or of sterile type, e.g. hypoxia, surgery) and the activation of metallo-proteases (Negrini et al., 1996; Negrini et al., 1998; Miserocchi et al., 1999; Passi et al., 1999). Fragmentation of intermediate/small PGs leads to progressive increase in microvascular permeability to water (increase in Lp) and solutes (decrease in σ). Figure 4 summarizes the events leading to the development of severe lung edema; in particular, fragmentation of large PGs vanish the increase in interstitial pressure (increase in tissue compliance) giving rise to uncontrolled filtration (Figures 4B, C).
FIGURE 4. Conceptual diagram showing pulmonary interstitial pressure (Pi) as a function of W/D. Control conditions (A), interstitial edema (B); the role of fragmentation of proteoglycans is highlighted in the development of severe edema (C).
Total proteoglycan content of the lung was evaluated by hexuronate assay on tissue samples. The technique to estimate the integrity of PGs was based on their extractability from the interstitial matrix using a chemical agent (guanidinium hydrochloride). It was found that proteoglycan extractability increased with increasing time using different models of edemagenic exposure, suggesting a weakening of noncovalent bonds linking intact proteoglycans to other interstitial matrix components (Negrini et al., 1996; Negrini et al., 1998; Passi et al., 1999; Miserocchi et al., 2001). Data allowed to derive the decrease of integrity of PGs, relative to a control value in physiological conditions (118 ± 4 μg/g wet tissue). Figure 5 LEFT shows that the loss of integrity of PGs correlates with increase in W/D ratio (points A, B and C refer to the conditions presented in Figure 4). Molecular size distribution of extracted PGs was also determined revealing a progressive decrease of large towards smaller molecular size suggesting fragmentation of native molecules (Negrini et al., 1998).
FIGURE 5. LEFT. Relationship between loss of integrity of interstitial hyaluronan and increase in W/D ratio of the lung. Points A, B and C refer to the conditions presented in Figure 4. RIGHT. Amount of extravascular lung water in the interstitial (Wint) and alveolar (Walv) compartments on increasing W/D (data from Beretta et al., 2021).
Figure 5 (RIGHT) shows the distribution of water in the interstitial (Wint) and in the alveolar compartment (Walv) with increasing W/D. Note that at W/D of ∼6.5, the capacity of the interstitial compartment for fluid accumulation (Wint) is saturated while alveolar flooding (Walv) will continue to increase on increasing W/D. Accordingly, a W/D of ∼6.5, corresponding to a loss of integrity of hyaluronan from 100 down to ∼55%, represents a potential threshold for alveolar flooding and thus development of severe edema. Exposure to 3–6 h of strongly edemagenic conditions impact remarkably on the integrity of the macromolecular PGs structure. Accordingly, identifying valid markers of developing edema appears crucial in the clinical setting as the time constant for developing of the severe stage ranges 4–7 min (Parker and Townsely, 2004; Mazzuca et al., 2016).
Deactivation of surfactant is a further cause leading to lung edema; this occurs on cyclic recruitment/de-recruitment of collapsed and/or fluid-filled alveoli that generate inter-alveolar forces that favour epithelial lesion and increase in microvascular permeability (Hamlington et al., 2018).
Heterogeneity of alveolar morphology, mechanical compliance and perfusion
We will now explore how the heterogeneity of the terminal respiratory units come into play when the lung is facing edemagenic conditions.
Heterogeneity in alveolar size has been reported in the healthy experimental animal based on direct imaging of subpleural alveoli and by optical coherence tomography (Meissner et al., 2010).
Figure 6 shows that the distribution of cross sectional area (data taken at equal height relative to the bottom of the lung) is remarkably right skewed at a distending pressure of 4 cmH2O; smaller alveoli (<2000 μm2) are ∼6 times more abundant than larger ones >4,000 μm2) (data from the experimental animal, Salito et al., 2014; Mazzuca et al., 2014).
FIGURE 6. Frequency distribution of cross-sectional area (A4) at a transpulmonary pressure of 4 cm H2O (data from experimental animal; from Salito et al., 2014; Mazzuca et al., 2014).
Heterogeneity in alveolar morphology implied heterogeneity in alveolar compliance in the range of transpulmonary pressure from 4 up to 8 cmH2O (Mazzuca et al., 2014; Salito et al., 2014).
Figure 7A shows trans-pleural images of alveoli of different size at distending pressure of 4 and 8 cmH2O. Figure 7B shows that, despite some variability, an index of absolute alveolar compliance (
FIGURE 7. (A) Trans-pleural images of alveoli of different cross-sectional area at distending pressure of 4 and 8 cmH2O. (B) Absolute alveolar compliance Cabs (ratio of change in alveolar surface to the change in distending pressure,
Interestingly, in lesional lung models causing severe edema (bleomycin treatment or surfactant deactivation by Tween 20), remarkable heterogeneity was observed ranging from atelectasis to over-distention of adjacent alveoli. This caused alveolar instability due to increased forces via mechanical interdependence exacerbating the spreading of lung edema during cycling recruitment/derecruitment of the alveoli (Bates & Smith, 2018; Bates et al., 2020; Gaver et al., 2020). Also in emphysema, a pathology implying high risk of edema, the heterogeneity in alveolar size, a specific trait of the disease, was greatly enhanced (Mishima et al., 1999). Further, a model simulation of heterogeneous distribution of molecular fragmentation leading to emphysema, correlated with the development of remarkable heterogeneity in alveolar size (Suki et al., 2003).
Concerning lung perfusion, extensive heterogeneities are reported in physiological conditions. Analysis of results in the upright posture, with and without gravity, and in the inverted, prone, and supine postures reveal significant flow heterogeneity and a persistent decrease in flow in the cranial and caudal regions for all postures suggesting that vascular geometry makes a major contribution to regional flow with gravity having a lesser role. A computational model of blood flow through the human pulmonary arterial tree suggested that, besides gravity and posture, the geometry of the branching structure plays a significant role in flow heterogeneity (Burrowes and Tawhai, 2006). Further, fractal geometry analysis showed that small asymmetries along the branching system might account for heterogeneity in both airways flow and blood flow distribution in terminal units (Glenny and Robertson, 2011). The last paper also confirmed that the regional lung perfusion reflects the complex interaction between vascular and surrounding alveolar pressures, particularly when local vasoconstriction occurs in response to hypoxia.
In interstitial edema, the immediate response is precapillary vasoconstriction, involving arterioles ∼80 µm in diameter (Negrini et al., 2001). The functional significance of this reflex is to avoid an increase of capillary pressure relative to the physiological value (∼10 cmH2O), particularly when capillary recruitment occurs on cardiac output as a response to increased oxygen requirement. A comprehensive model of fluid exchange at capillary level (based on Eq. 1, suggests indeed that doubling capillary pressure (say from 10 up to 20 cmH2O) is, by and large, more edemagenic compared to an eight fold increase in water permeability Lp (Mazzuca et al., 2016). The reason is that the increase in pressure contributes to capillary recruitment, thus implying also an increase in Kf due to increase A (see Eq. 1).
The lung vasomotion in distribution vessels <100 µm in response to a strong edemagenic factor (hypoxia) is shown in Figure 8, actually revealing two different patterns: either no change, or remarkable decrease in diameter approaching complete closure after 30 min of hypoxia exposure (Mazzuca et al., 2019). Closure was followed by return towards control at about 80 min, remaining thereafter essentially steady. On average, considering all the data, a significant decrease in diameter down to about 50% of control was observed at 30 min.
FIGURE 8. Time course of vasomotion in arterial pulmonary distribution vessels with diameter <100 nm in response to hypoxia (12% oxygen) (from Mazzuca et al., 2019).
Greater vasoconstriction occurs in regions developing edema with redirection of blood flow from edematous to non edematous regions (Rivolta et al., 2011; Mazzuca et al., 2019).
An attempt was also done to model the alveolar capillary network after identifying corner vessels and assuming a random distribution of Voronoi points (Mazzuca et al., 2019); the number of points was chosen so as to match the ratio between capillary volume and interstitial volume. The functional significance of this ratio is that to balance the filtration surface of the capillary network with the capacity of the surrounding interstitial compartment to host microvascular filtration. An estimate of this ratio, based on TEM images (Conforti et al., 2002), was set at 1.54.
Figure 9 shows an example of the capillary perfusion pattern with a color-coded log-scale in alveolar districts after inducing hypoxic edema (Mazzuca et al., 2019). Data refer to two different networks indicated as 1 and 2. Progressive closure of microvessels (blue colour) was observed over time in alveolar regions where interstitial edema was more marked, as judged from the thickness of the perivascular interstitial space. Blood flow switching was described as reflecting the geometry of branching, as well as the transcapillary pressure (capillary pressure—interstitial pressure) along the longitudinal profile of the capillary. The model supports the heterogeneous distribution of edema further suggesting that capillary closure is favoured by the compressive action of the increase in interstitial pressure, as predicted by the model of Glenny and Robertson, 2011). The heterogeneity in edema development was found to correlate with uneven pulmonary blood flow distribution (Urano et al., 2005). Capillary de-recruitment due to perivascular edema may well occur also in humans, as administration of a vasodilator agent was unable to restore blood flow in edematous lung regions (Scherrer et al., 1996).
FIGURE 9. Perfusion analysis for alveolar-capillary units relative to two networks (1 and 2). Red and light blue dots identify respectively arteriolar accesses and venular exits. The color panels show the capillary blood flow at different time points (baseline, 30 min and 120 min of hypoxia exposure) with color-coded log-scale intensity (from Mazzuca et al., 2019).
A recent study (Buxton et al., 2022) provides data from a novel analysis on the spatial/temporal dynamics of blood flow in the human lung measured noninvasively, potentially useful to estimate the vascular control mechanisms involved.
Finally, semi-quantitative estimates of microvascular filtration were obtained from a theoretical model of the morphologically based–alveolar capillary unit (Mazzuca et al., 2016). Giving the heterogeneous distribution of blood flow, it was found that the filtration flow was also heterogeneously distributed. A key point of this analysis was that greater microvascular filtration occurred around alveoli of larger size. Although TEM images were not available, one may hypothesize that, on geometrical ground, the capillary surface available for filtration in large alveoli is actually oversized compared to the capacity of the interstitial compartment to host microvascular filtration.
Inter-individual difference in human alveolar-capillary phenotype as potential factor to develop lung edema
As long as the architecture of the air-blood barrier is preserved, it appears to be very effective in optimizing the gas exchange (Swan and Tawhai, 2011). It is a common experience however that, in edemagenic conditions, some subjects are more prone to develop lung edema compared to others. The characteristic example is the exposure to hypoxia at high altitude that allowed to distinguish between HAPE-S (High Altitude Pulmonary Edema-Sensitive) and HAPE-R (High Altitude Pulmonary Edema-Resistant) subjects (Busch et al., 2001; Dehnert et al., 2006; Richalet et al., 2012; 2021; Eichstaedt et al., 2020a; b). Further, studies based on the MIGET technique (Multiple Inert Gas Elimination Technique, Wagner, 2008) have reported remarkable inter-individual differences in the perturbation of the distribution of the VA/Q ratio that were attributed to regional development of lung edema in severe edemagenic conditions (exercise in hypoxia) (Hammond et al., 1986; Hopkins et al., 1996; Hopkins, 2005; Hopkins, 2020).
Since lung edema is primarily affecting the alveolar-capillary compartment, a line of research aimed at characterizing inter-individual differences concerning the oxygen diffusion/transport function. At morpho-functional level, subject were found to differ remarkably considering the ratio between the pulmonary capillary blood volume (Vc) and the alveolar membrane diffusive capacitance (Dm) (Miserocchi et al., 2008). At the functional respiratory capacity (FRC, corresponding to end-expiration), the Vc/Dm ratio was found to vary from ∼ 1 up to ∼ 6, with a normal distribution (Miserocchi et al., 2022). It was considered that subjects with low Vc/Dm have a less extended capillary network relative to the surface of the air-blood barrier; on the opposite, subjects with high Vc/Dm have a more extended capillary network compared to the surface of the air-blood barrier. At FRC people with low Vc/Dm ratio, had, on the average, a 3 times greater Dm value (Miserocchi et al., 2008) likely reflecting a greater alveolar surface. It was considered that a higher value of Dm would be justified by a greater number of alveoli with smaller radius (set as n and r, respectively), compared to the case of low Dm implying alveoli of smaller number but greater radius (N and R, respectively). For a ratio between high and low Dm set at 3, and assuming either a spherical shape or regular polyhedron for the alveoli, the ratio of alveolar surfaces can be numerically modelled, being proportional to the product
FIGURE 10. (A) Relative increase in the ratio of number of alveoli on decreasing
One can therefore hypothesize that a greater capacity of the peri-alveolar interstitial volume to accept fluid in edemagenic condition would help to prevent alveolar edema; this would be the case for a greater number of alveoli of smaller radius. This concept fits with the experimental finding of greater interstitial fluid accumulation around large alveoli (Mazzuca et al., 2019).
The potential advantage of small alveoli in favouring gas exchange has been considered based on a geometrical approach similar to that proposed for the interpretation of differences in Dm. A critical parameter optimizing gas-exchange would be the perimeter of the surface on a planar cut of the acinus; the conclusion by the authors is that “the best possible acinus” should be small so that the length of its perimeter is of the order of the length travelled by a gas molecule before being absorbed (Sapoval et al., 2002; Felici et al., 2003).
Considering oxygen uptake, model by Swan and Tawhai, 2011, suggested that, despite some morphological heterogeneity, it is very efficient at rest, with only a small degree of intra-acinar PO2 heterogeneity. Further, it was recognized that anisotropy due to deformation causing non-uniformity of the tissue stretch among alveoli, could become important to hinder oxygen uptake only at high rates of ventilation (Swan et al., 2012).
Inter-individual heterogeneity of the air-blood barrier in humans and its impact on oxygen uptake-transport in edemagenic conditions
It is easy to hypothesize that edema developing at the level of the air-blood barrier might affect the efficiency of oxygen diffusion/transport. Given the differences in human alveolo-capillary phenotype, as suggested by the Vc/Dm ratio (Miserocchi et al., 2008), and the potential correlation with the control of lung fluid balance, a line of research was dedicated to identify inter-individual differences in the response to edemagenic conditions. The study relied on a biophysical model defining the alveolar-capillary oxygen equilibration (Leq) in the air-blood barrier (Piiper and Scheid, 1981) as:
being DO2 the O2 diffusive capacity,
Leq defines the equilibrium reached in the blood exiting the pulmonary capillaries between the PO2 in the arterial blood relative to the alveolar PO2. Leq varies from 0 (the case where PO2 in the arterial blood attains alveolar PO2) to 1 (the case of arterial-venous shunt).
The model assumes a kinetics of equilibration of exponential type, accordingly Eq. (2) can be rewritten as:
being Tt the blood transit time in the pulmonary capillaries that can be estimated as
All variables appearing in Eqs. 2, 3 could be estimated allowing to describe the kinetics of alveolar-capillary equilibration. This approach, reflecting the combined changes of the variables involved, allowed to characterize at individual level the efficiency of the oxygen diffusion-transport in edemagenic condition. The analysis was carried out in subjects at rest and exposed to strongly edemagenic conditions (work at 3,840 m, PIO2 ∼90 mmHg) (Beretta et al., 2017; Beretta et al., 2019). Subjects with low Vc/Dm at rest still attained Leq ∼ 0. Conversely, subjects having a high Vc/Dm showed a remarkable pulmonary capillary de-recruitment (decrease in Vc) and a decrease in DO2. The latter could be attributed to a perturbation in lung fluid balance that was assessed by the changes in mechanical properties of the respiratory system estimated with the frequency oscillation technique (Dellacà et al., 2008; Bartesaghi et al., 2014). Figure 11A shows that the distribution of Leq among subjects in severe edemagenic conditions is normal. The blue point at the left refers to a subject with Vc/Dm ∼ 1 at FRC at rest attaining an optimal equilibration for O2 exchange (no perturbation in lung water). The red point to the right refers to a subjects with Vc/Dm ∼ 4.3 showing that equilibration only achieved 45%. Due to impairment of the gas exchange, subjects with high Vc/Dm were forced to increase cardiac output. This reflex response, adding to the de-recruitment of the capillary network, led to an increase in blood velocity and shortening of
FIGURE 11. (A) Distribution of Leq among subjects exposed to severe edemagenic conditions (work in hypoxia). The blue point at the left refers to a subject with Vc/Dm ∼ 1 at FRC at rest attaining an optimal equilibration for O2 exchange (with no perturbation in lung water), while the red point to right refers to a subjects with Vc/Dm ∼ 4.3 showing that equilibration was only achieved at ∼45%. (B) Plot of Leq vs. 1/Tt (an index of blood velocity). Red and blue dots refer to the same subjects as in (A) (from: Miserocchi et al., 2022).
Figure 11B supports the functional correlation between the increase in Leq and an index of increase in capillary blood velocity expressed by the reciprocal of transit time 1/Tt.
In general, one shall note that
Besides alveolar morphological heterogeneity one shall mention a further pathophysiological factor to possibly explain a greater proneness to develop lung edema, namely an inborn greater microvascular permeability. A basic, so far unanswered pathophysiological question, remains indeed at clinical level: why does pulmonary hypertension develop in intrinsic lung disease such as idiopatic pulmonary fibrosis, interstitial pulmonary fibrosis associated with pulmonary artery fibrosis, combined pulmonary fibrosis and emphysema (Munson, 2010)? A possible answer would come by hypothesizing a common pathophysiological interpretation for these conditions. One may comment that precapillary vasoconstriction, as a reflex protective response against lung edema, inevitably leads to pulmonary hypertension. One could hypothesize that the deposition of fibrotic tissue could also be regarded as a long-term functional adaptive response to protect against the risk of lung edema in case of high inborn microvascular permeability. In support of this interpretation is the finding that the relative ratio of capillary area to total area of alveolar walls (the case of high Vc/Dm) was significantly higher at low grades of fibrosis than in control healthy lungs (Ebina et al., 2004). The same study reports that vascular density gradually decreased with increasing the degree of fibrosis becoming definitely lower than that of control in the most extensive lung fibrotic lesions. Obviously, if decreasing capillary surface helps preventing the occurrence of lung edema, on the other side, the complication is that excessive deposition of fibrotic tissue along the arterial vascular segment occurs at the expense of increasing right-ventricular afterload (a maladaptive case of compensation). The transdifferentiation of lung fibroblasts into myofibroblasts proliferation and extensive deposition of extracellular matrix are currently the object of growing interest from a therapeutic point of view in the pathogenesis of pulmonary fibrosis (Gan et al., 2022; Li et al., 2022; Gu et al., 2023; Qi et al., 2022).
Author contributions
GM conceived and wrote the paper. The author confirms being the sole contributor of this work and has approved it for publication.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author GM declared that he was an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Bartesaghi, M., Beretta, E., Pollastri, L., Scotti, V., Mandolesi, G., Lanfranconi, F., et al. (2014). Inter-individual differences in control of alveolar capillary blood volume in exercise and hypoxia. Respir. Physiol. Neurobiol. 190, 96–104. doi:10.1016/j.resp.2013.08.021
Bates, J., Gaver, D. P., Habashi, N. M., and Nieman, G. F. (2020). Atelectrauma versus volutrauma: A tale of two time-constants. Crit. Care Explor. 2, e0299. doi:10.1097/CCE.0000000000000299
Bates, J., and Smith, B. J. (2018). Ventilator-induced lung injury and lung mechanics. Ann. Transl. Med. 6, 378. doi:10.21037/atm.2018.06.29
Beretta, E., Grasso, G. S., Forcaia, G., Sancini, G., and Miserocchi, G. (2019). Differences in alveolo-capillary equilibration in healthy subjects on facing O2demand. Sci. Rep. Nov. 13 (9), 16693. PMID: 31723148; PMCID: PMC6854051. doi:10.1038/s41598-019-52679-4
Beretta, E., Lanfranconi, F., Grasso, G. S., Bartesaghi, M., Alemayehu, H. K., Pratali, L., et al. (2017). Air blood barrier phenotype correlates with alveolo-capillary O2 equilibration in hypobaric hypoxia. Respir. Physiol. Neurobiol. 246, 53–58. doi:10.1016/j.resp.2017.08.006
Beretta, E., Romanò, F., Sancini, G., Grotberg, J. B., Nieman, G. F., and Miserocchi, G. (2021). Pulmonary interstitial matrix and lung fluid balance from normal to the acutely injured lung. Front. Physiol. 12, 781874. doi:10.3389/fphys.2021.781874
Botto, L., Beretta, E., Daffara, R., Miserocchi, G., and Palestini, P. (2006). Biochemical and morphological changes in endothelial cells in response to hypoxic interstitial edema. Respir. Res. 7, 7. doi:10.1186/1465-9921-7-7
Burrowes, K., and Tawhai, M. (2006). Computational predictions of pulmonary blood flow gradients: Gravity versus structure. Respir. Physiol. Neurobiol. 154 (3), 515–523. Epub 2005 Dec 28. PMID: 16386472. doi:10.1016/j.resp.2005.11.007
Busch, T., Bärtsch, P., Pappert, D., Grünig, E., Hildebrandt, W., Elser, H., et al. (2001). Hypoxia decreases exhaled nitric oxide in mountaineers susceptible to high-altitude pulmonary edema. Am. J. Respir. Crit. Care Med. 163 (2), 368–373. doi:10.1164/ajrccm.163.2.2001134
Buxton, R. B., Prisk, G. K., and Hopkins, S. R. (2022). A novel nonlinear analysis of blood flow dynamics applied to the human lung. J. Appl. Physiol. (1985) 132 (6), 1546–1559. Epub 2022 Apr 14. PMID: 35421317; PMCID: PMC9190737. doi:10.1152/japplphysiol.00715.2021
Cavalcante, F., Ito, S., Brewer, K., Sakai, H., Alencar, A., Almeida, M., et al. (2005). Mechanical interactions between collagen and proteoglycans: Implications for the stability of lung tissue. J. Appl. Physiol. 98, 672–679. doi:10.1152/japplphysiol.00619.2004
Chen, Q., and Liu, Y. (2020). Heterogeneous groups of alveolar type II cells in lung homeostasis and repair. Am. J. Physiol. Cell. Physiol. 1 319 (6), C991–C996. Epub 2020 Sep 9. PMID: 32903031; PMCID: PMC7768230. doi:10.1152/ajpcell.00341.2020
Conforti, E., Fenoglio, C., Bernocchi, G., Bruschi, O., and Miserocchi, G. A. (2002). Morpho-functional analysis of lung tissue in mild interstitial edema. Am. J. Physiol. Lung Cell. Mol. Physiol. 282, L766–L774. doi:10.1152/ajplung.00313.2001
Dehnert, C., Risse, F., Ley, S., Kuder, T. A., Buhmann, R., Puderbach, M., et al. (2006). Magnetic resonance imaging of uneven pulmonary perfusion in hypoxia in humans. Am. J. Respir. Crit. Care Med. 174 (10), 1132–1138. doi:10.1164/rccm.200606-780OC
Dellacà, R. L., Zannin, E., Sancini, G., Rivolta, I., Leone, B. E., Pedotti, A., et al. (2008). Changes in the mechanical properties of the respiratory system during the development of interstitial lung edema. Respir. Res. 9 (1), 51. doi:10.1186/1465-9921-9-51
Ebina, M., Shimizukawa, M., Shibata, N., Kimura, Y., Suzuki, T., Endo, M., et al. (2004). Heterogeneous increase in CD34-positive alveolar capillaries in idiopathic pulmonary fibrosis. Am. J. Respir. Crit. Care. Med. 169 (11), 1203–1208. Epub 2004 Jan 30. PMID: 14754760. doi:10.1164/rccm.200308-1111OC
Eichstaedt, C. A., Benjamin, N., and Grünig, E. (2020a). Genetics of pulmonary hypertension and high-altitude pulmonary edema. J. Appl. Physiol. (1985) 128 (5), 1432–1438. doi:10.1152/japplphysiol.00113.2020
Eichstaedt, C. A., Mairbäurl, H., Song, J., Benjamin, N., Fischer, C., Dehnert, C., et al. (2020b). Genetic predisposition to high-altitude pulmonary edema. High. Alt. Med. Biol. 21 (1), 28–36. doi:10.1089/ham.2019.0083
Felici, M., Filoche, M., and Sapoval, B. (2003). Diffusional screening in the human pulmonary acinus. J. Appl. Physiol. (1985) 94 (5), 2010–2016. doi:10.1152/japplphysiol.00913.2002
Gan, C., Zhang, Q., Liu, H., Wang, G., Wang, L., Li, Y., et al. (2022). Nifuroxazide ameliorates pulmonary fibrosis by blocking myofibroblast Genesis: A drug repurposing study. Respir. Res. 23 (1), 32. doi:10.1186/s12931-022-01946-6
Gaver, D. P., Nieman, G. F., Gatto, L. A., Cereda, M., Habashi, N. M., and Bates, J. (2020). The POOR get POORer: A hypothesis for the pathogenesis of ventilator-induced lung injury. Am. J. Respir. Crit. Care Med. 202, 1081–1087. doi:10.1164/rccm.202002-0453CP
Glenny, R., and Robertson, H. (2011). Distribution of perfusion. Compr. Physiol. Jan. 1 (1), 245–262. doi:10.1002/cphy.c100012
Gu, S., Liang, J., Zhang, J., Liu, Z., Miao, Y., Wei, Y., et al. (2023). Baricitinib attenuates bleomycin-induced pulmonary fibrosis in mice by inhibiting TGF-β1 signaling pathway. Molecules 28 (5), 2195. Published Feb 27. doi:10.3390/molecules28052195
Hamlington, K. L., Bates, J., Roy, G. S., Julianelle, A. J., Charlebois, C., Suki, B., et al. (2018). Alveolar leak develops by a rich-get-richer process in ventilator-induced lung injury. PLoS One 13 (3), e0193934. doi:10.1371/journal.pone.0193934
Hammond, M. D., Gale, G. E., Kapitan, K. S., Ries, A., and Wagner, P. D. (1986). Pulmonary gas exchange in humans during normobaric hypoxic exercise. J. Appl. Physiol. (1985) 61 (5), 1749–1757. doi:10.1152/jappl.1986.61.5.1749
Hopkins, S. R., Belzberg, A. S., Wiggs, B. R., and McKenzie, D. C. (1996). Pulmonary transit time and diffusion limitation during heavy exercise in athletes. Respir. Physiol. 103 (1), 67–73. doi:10.1016/0034-5687(95)00028-3
Hopkins, S. R. (2005). The lung at maximal exercise: Insights from comparative physiology. Clin. Chest Med. 26 (3), 459–468. doi:10.1016/j.ccm.2005.05.003
Hopkins, S. R. (2020). Ventilation/perfusion relationships and gas exchange: Measurement approaches. Compr. Physiol. 10 (3), 1155–1205. doi:10.1002/cphy.c180042
Hopkins, S. R. (2010). Point: Pulmonary edema does occur in human athletes performing heavy sea-level exercise. J. Appl. Physiol. (1985) 109 (4), 1270–1272. doi:10.1152/japplphysiol.01353.2009
Knudsen, L., and Ochs, M. (2018). The micromechanics of lung alveoli: Structure and function of surfactant and tissue components. Histochem. Cell. Biol. Dec 150 (6), 661–676. Epub 2018 Nov 2. PMID: 30390118; PMCID: PMC6267411. doi:10.1007/s00418-018-1747-9
Lettau, M., Timm, S., Dittmayer, C., Lopez-Rodriguez, E., and Ochs, M. (2022). The ultrastructural heterogeneity of lung surfactant revealed by serial section electron tomography: Insights into the 3-D architecture of human tubular myelin. Am. J. Physiol. Lung Cell. Mol. Physiol. 1 322 (6), L873–L881. Epub 2022 Apr 19. PMID: 35438000. doi:10.1152/ajplung.00020.2022
Levick, J. R., and Michel, C. C. (2010). Microvascular fluid exchange and the revised Starling principle. Cardiovasc. Res. 87, 198–210. doi:10.1093/cvr/cvq062
Li, Y., Sun, Y., Wu, N., and Ma, H. (2022). GRK2 promotes activation of lung fibroblast cells and contributes to pathogenesis of pulmonary fibrosis through increasing Smad3 expression. Am. J. Physiol. Cell Physiol. 322 (1), C63–C72. doi:10.1152/ajpcell.00347.2021
Lynch, T., Ievlev, V., and Parekh, K. (2019). Heterogeneity of pulmonary stem cells. Adv. Exp. Med. Biol. 1169, 95–117. PMID: 31487021; PMCID: PMC8421072. doi:10.1007/978-3-030-24108-7_6
Ma, B., and Bates, J. (2014). Mechanical interactions between adjacent airways in the lung. J. Appl. Physiol. 116 (6), 628–634. doi:10.1152/japplphysiol.01180.2013
Mazzuca, E., Aliverti, A., and Miserocchi, G. (2016). Computational micro-scale model of control of extravascular water and capillary perfusion in the air blood barrier. J. Theor. Biol. 400, 42–51. doi:10.1016/j.jtbi.2016.03.036
Mazzuca, E., Aliverti, A., and Miserocchi, G. (2019). Understanding vasomotion of lung microcirculation by in vivo imaging. Imaging 5 (2), 22. doi:10.3390/jimaging5020022
Mazzuca, E., Salito, C., Rivolta, I., Aliverti, A., and Miserocchi, G. (2014). From morphological heterogeneity at alveolar level to the overall mechanical lung behavior: An in vivo microscopic imaging study. Physiol. Rep. 7 (2), e00221. PMID: 24744890; PMCID: PMC3966245. doi:10.1002/phy2.221
Meissner, S., Tabuchi, A., Mertens, M., Kuebler, W., and Koch, E. (2010). Virtual four-dimensional imaging of lung parenchyma by optical coherence tomography in mice. J. Biomed. Opt. 15 (3), 036016. doi:10.1117/1.3425654
Miserocchi, G., Beretta, E., Rivolta, I., and Bartesaghi, M. (2022). Role of the air-blood barrier phenotype in lung oxygen uptake and control of extravascular water. Front. Physiol. 28 (13), 811129. PMID: 35418875; PMCID: PMC8996119. doi:10.3389/fphys.2022.811129
Miserocchi, G. (2009). Mechanisms controlling the volume of pleural fluid and extravascular lung water. Eur. Respir. Rev. 18 (114), 244–252. doi:10.1183/09059180.00002709
Miserocchi, G., Messinesi, G., Tana, F., Passoni, E., Adamo, S., Romano, R., et al. (2008). Mechanisms behind inter-individual differences in lung diffusing capacity. Eur. J. Appl. Physiol. 102 (5), 561–568. doi:10.1007/s00421-007-0625-2
Miserocchi, G., Negrini, D., Del Fabbro, M., and Venturoli, D. (1993). Pulmonary interstitial pressure in intact in situ lung: Transition to interstitial edema. J. Appl. Physiol. (1985) 74 (3), 1171–1177. doi:10.1152/jappl.1993.74.3.1171
Miserocchi, G., Passi, A., Albertini, R., Negrini, D., and De Luca, G. (1999). Interstitial pressure and proteoglycan degradation in hydraulic- and elastase-induced lung edema. Chest 116 (1), 31S. doi:10.1378/chest.116.suppl_1.31s
Miserocchi, G., Passi, A., Negrini, D., Del Fabbro, M., and De Luca, G. (2001). Pulmonary interstitial pressure and tissue matrix structure in acute hypoxia. Am. J. Physiol. Lung Cell Mol. Physiol. 280 (5), L881–L887. doi:10.1152/ajplung.2001.280.5.L881
Mishima, M., Hirai, T., Itoh, H., Nakano, Y., Sakai, H., Muro, S., et al. (1999). Complexity of terminal airspace geometry assessed by lung computed tomography in normal subjects and patients with chronic obstructive pulmonary disease. Proc. Natl. Acad. Sci. U. S. A. 96 (16), 8829–8834. doi:10.1073/pnas.96.16.8829
Munson, J. C. (2010). Combined pulmonary fibrosis and emphysema: A high-pressure situation. Eur. Respir. J. 35 (1), 9–11. doi:10.1183/09031936.00156909
Negrini, D., Passi, A., de Luca, G., and Miserocchi, G. (1996). Pulmonary interstitial pressure and proteoglycans during development of pulmonary edema. Am. J. Physiol. 270 (6), H2000–H2007. doi:10.1152/ajpheart.1996.270.6.H2000
Negrini, D., Passi, A., De Luca, G., and Miserochi, G. (1998). Proteoglycan involvement during development of lesional pulmonary edema. Am. J. Physiol. 274 (2), L203–L211. doi:10.1152/ajplung.1998.274.2.L203
Negrini, D., Candiani, A., Boschetti, F., Crisafulli, B., Del Fabbro, M., Bettinelli, D., et al. (2001). Pulmonary microvascular and perivascular interstitial geometry during development of mild hydraulic edema. Am. J. Physiol. Lung Cell. Mol. Physiol. 281 (6), L1464–L1471. doi:10.1152/ajplung.2001.281.6.L1464
Parker, J. C. (2007). Hydraulic conductance of lung endothelial phenotypes and Starling safety factors against edema. Am. J. Physiol. Lung Cell. Mol. Physiol. 292 (2), L378–L380. doi:10.1152/ajplung.00196.2006
Parker, J. C., Stevens, T., Randall, J., Weber, D. S., and King, J. A. (2006). Hydraulic conductance of pulmonary microvascular and macrovascular endothelial cell monolayers. Am. J. Physiol. Lung Cell Mol. Physiol. 291, L30–L37. doi:10.1152/ajplung.00317.2005
Parker, J. C., and Townsley, M. I. (2004). Evaluation of lung injury in rats and mice. Am. J. Physiol. Lung Cell. Mol. Physiol. 286, L231–L246. doi:10.1152/ajplung.00049.2003
Passi, A., Negrini, D., Albertini, R., Miserocchi, G., and De Luca, G. (1999). The sensitivity of versican from rabbit lung to gelatinase A (MMP-2) and B (MMP-9) and its involvement in the development of hydraulic lung edema. FEBS Lett. 456 (1), 93–96. doi:10.1016/s0014-5793(99)00929-1
Piiper, J., and Scheid, P. (1981). Model for capillary-alveolar equilibration with special reference to O2 uptake in hypoxia. Respir. Physiol. 46 (3), 193–208. doi:10.1016/0034-5687(81)90121-3
Qi, Y., Zhang, H., Fan, H., Wang, X., Zhao, A., Tian, Y., et al. (2022). PPARγ/LXRα axis mediated phenotypic plasticity of lung fibroblasts in silica-induced experimental silicosis. Environ. Pollut. 292, 118272. doi:10.1016/j.envpol.2021.118272
Richalet, J. P., Larmignat, P., Poitrine, E., Letournel, M., and Canouï-Poitrine, F. (2012). Physiological risk factors for severe high-altitude illness: A prospective cohort study. Am. J. Respir. Crit. Care Med. 185 (2), 192–198. doi:10.1164/rccm.201108-1396OC
Richalet, J. P., Pillard, F., Le Moal, , Rivière, D., Oriol, P., Poussel, M., et al. (2021). Validation of a score for the detection of subjects with high risk for severe high-altitude illness. Med. Sci. Sports Exerc. 53 (6), 1294–1302. doi:10.1249/MSS.0000000000002586
Rivolta, I., Lucchini, V., Rocchetti, M., Kolar, F., Palazzo, F., Zaza, A., et al. (2011). Interstitial pressure and lung oedema in chronic hypoxia. Eur. Respir. J. 37 (4), 943–949. doi:10.1183/09031936.00066710
Salito, C., Aliverti, A., Mazzuca, E., Rivolta, I., and Miserocchi, G. (2015). The effect of exogenous surfactant on alveolar interdependence. Respir. Physiol. Neurobiol. May 210, 7–13. Epub 2015 Jan 17. PMID: 25600053. doi:10.1016/j.resp.2015.01.009
Salito, C., Aliverti, A., Rivolta, I., Mazzuca, E., and Miserocchi, G. (2014). Alveolar mechanics studied by in vivo microscopy imaging through intact pleural space. Respir. Physiol. Neurobiol. 1, 44–49. Epub 2014 Jul 22. PMID: 25058162. doi:10.1016/j.resp.2014.07.012
Sapoval, B., Filoche, M., and Weibel, E. R. (2002). Smaller is better–but not too small: A physical scale for the design of the mammalian pulmonary acinus. Proc. Natl. Acad. Sci. U.S.A. 99, 10411–10416. doi:10.1073/pnas.122352499
Scherrer, U., Vollenweider, L., Delabays, A., Savcic, M., Eichenberger, U., Kleger, G. R., et al. (1996). Inhaled nitric oxide for high-altitude pulmonary edema. N. Engl. J. Med. 334 (10), 624–629. doi:10.1056/NEJM199603073341003
Schraufnagel, D. E., Agaram, N. P., Faruqui, , Jain, S., Jain, L., Ridge, K. M., et al. (2003). Pulmonary lymphatics and edema accumulation after brief lung injury. Am. J. Physiol. Lung Cell. Mol. Physiol. 284, L891–L897. doi:10.1152/ajplung.00333.2002
Suki, B., Lutchen, K. R., and Ingenito, E. P. (2003). On the progressive nature of emphysema: Roles of proteases, inflammation, and mechanical forces. Am. J. Respir. Crit. Care Med. 168 (5), 516–521. doi:10.1164/rccm.200208-908PP
Swan, A., Clark, A., and Tawhai, M. (2012). A computational model of the topographic distribution of ventilation in healthy human lungs. J. Theor. Biol. 7 (300), 222–231. doi:10.1016/j.jtbi.2012.01.042
Swan, A. J., and Tawhai, M. H. (2011). Evidence for minimal oxygen heterogeneity in the healthy human pulmonary acinus. J. Appl. Physiol. (1985) 110 (2), 528–537. Epub 2010 Nov 11. PMID: 21071589; PMCID: PMC3043789. doi:10.1152/japplphysiol.00888.2010
Urano, T., Kuwahira, I., Iwamoto, T., Kamiya, U., Ohta, Y., Wood, J. G., et al. (2005). Exposure to hypoxia results in uneven pulmonary blood flow distribution prior to pulmonary edema. Tokai J. Exp. Clin. Med. 30 (4), 193–202. PMID: 16482939.
Wagner, P. D. (2008). The multiple inert gas elimination technique (MIGET). Intensive Care Med. 34 (6), 994–1001. doi:10.1007/s00134-008-1108-6
Keywords: Starling law, edema, lung diffusion, lung perfusion, lung interstitial matrix
Citation: Miserocchi G (2023) The impact of heterogeneity of the air-blood barrier on control of lung extravascular water and alveolar gas exchange. Front. Netw. Physiol. 3:1142245. doi: 10.3389/fnetp.2023.1142245
Received: 11 January 2023; Accepted: 28 April 2023;
Published: 11 May 2023.
Edited by:
Akira Tsuda, Tsuda Lung Research, United StatesReviewed by:
Bindi Shah Brook, University of Nottingham, United KingdomLars Knudsen, Hannover Medical School, Germany
Copyright © 2023 Miserocchi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Giuseppe Miserocchi, Z2l1c2VwcGUubWlzZXJvY2NoaUB1bmltaWIuaXQ=