Why is VO2 max after altitude acclimatization still reduced despite normalization of arterial O2 content?

[EN] Acute hypoxia (AH) reduces maximal O2 consumption (VO2 max), but after acclimatization, and despite increases in both hemoglobin concentration and arterial O2 saturation that can normalize arterial O2 concentration ([O2]), VO2 max remains low. To determine why, seven lowlanders were studied at VO2 max (cycle ergometry) at sea level (SL), after 9-10 wk at 5,260 m [chronic hypoxia (CH)], and 6 mo later at SL in AH (FiO2 = 0.105) equivalent to 5,260 m. Pulmonary and leg indexes of O2 transport were measured in each condition. Both cardiac output and leg blood flow were reduced by approximately 15% in both AH and CH (P < 0.05). At maximal exercise, arterial [O2] in AH was 31% lower than at SL (P < 0.05), whereas in CH it was the same as at SL due to both polycythemia and hyperventilation. O2 extraction by the legs, however, remained at SL values in both AH and CH. Although at both SL and in AH, 76% of the cardiac output perfused the legs, in CH the legs received only 67%. Pulmonary VO2 max (4.1 +/- 0.3 l/min at SL) fell to 2.2 +/- 0.1 l/min in AH (P < 0.05) and was only 2.4 +/- 0.2 l/min in CH (P < 0.05). These data suggest that the failure to recover VO2 max after acclimatization despite normalization of arterial [O2] is explained by two circulatory effects of altitude: 1) failure of cardiac output to normalize and 2) preferential redistribution of cardiac output to nonexercising tissues. Oxygen transport from blood to muscle mitochondria, on the other hand, appears unaffected by CH.

Determinants of maximal oxygen uptake in severe acute hypoxia

[EN] To unravel the mechanisms by which maximal oxygen uptake (VO2 max) is reduced with severe acute hypoxia in humans, nine Danish lowlanders performed incremental cycle ergometer exercise to exhaustion, while breathing room air (normoxia) or 10.5% O2 in N2 (hypoxia, approximately 5,300 m above sea level). With hypoxia, exercise PaO2 dropped to 31-34 mmHg and arterial O2 content (CaO2) was reduced by 35% (P < 0.001). Forty-one percent of the reduction in CaO2 was explained by the lower inspired O2 pressure (PiO2) in hypoxia, whereas the rest was due to the impairment of the pulmonary gas exchange, as reflected by the higher alveolar-arterial O2 difference in hypoxia (P < 0.05). Hypoxia caused a 47% decrease in VO2 max (a greater fall than accountable by reduced CaO2). Peak cardiac output decreased by 17% (P < 0.01), due to equal reductions in both peak heart rate and stroke VOlume (P < 0.05). Peak leg blood flow was also lower (by 22%, P < 0.01). Consequently, systemic and leg O2 delivery were reduced by 43 and 47%, respectively, with hypoxia (P < 0.001) correlating closely with VO2 max (r = 0.98, P < 0.001). Therefore, three main mechanisms account for the reduction of VO2 max in severe acute hypoxia: 1) reduction of PiO2, 2) impairment of pulmonary gas exchange, and 3) reduction of maximal cardiac output and peak leg blood flow, each explaining about one-third of the loss in VO2 max.

Parasympathetic neural activity accounts for the lowering of exercise heart rate at high altitude

[EN] BACKGROUND: In chronic hypoxia, both heart rate (HR) and cardiac output (Q) are reduced during exercise. The role of parasympathetic neural activity in lowering HR is unresolved, and its influence on Q and oxygen transport at high altitude has never been studied. METHODS AND RESULTS: HR, Q, oxygen uptake, mean arterial pressure, and leg blood flow were determined at rest and during cycle exercise with and without vagal blockade with glycopyrrolate in 7 healthy lowlanders after 9 weeks' residence at >/=5260 m (ALT). At ALT, glycopyrrolate increased resting HR by 80 bpm (73+/-4 to 153+/-4 bpm) compared with 53 bpm (61+/-3 to 114+/-6 bpm) at sea level (SL). During exercise at ALT, glycopyrrolate increased HR by approximately 40 bpm both at submaximal (127+/-4 to 170+/-3 bpm; 118 W) and maximal (141+/-6 to 180+/-2 bpm) exercise, whereas at SL, the increase was only by 16 bpm (137+/-6 to 153+/-4 bpm) at 118 W, with no effect at maximal exercise (181+/-2 bpm). Despite restoration of maximal HR to SL values, glycopyrrolate had no influence on Q, which was reduced at ALT. Breathing FIO(2)=0.55 at peak exercise restored Q and power output to SL values. CONCLUSIONS: Enhanced parasympathetic neural activity accounts for the lowering of HR during exercise at ALT without influencing Q. The abrupt restoration of peak exercise Q in chronic hypoxia to maximal SL values when arterial PO(2) and SO(2) are similarly increased suggests hypoxia-mediated attenuation of Q.

Cardiovascular responses to dynamic exercise with acute anemia in humans

[EN] We hypothesized that reducing arterial O2 content (CaO2) by lowering the hemoglobin concentration ([Hb]) would result in a higher blood flow, as observed with a low PO2, and maintenance of O2 delivery. Seven young healthy men were studied twice, at rest and during two-legged submaximal and peak dynamic knee extensor exercise in a control condition (mean control [Hb] 144 g/l) and after 1-1.5 liters of whole blood had been withdrawn and replaced with albumin [mean drop in [Hb] 29 g/l (range 19-38 g/l); low [Hb]]. Limb blood flow (LBF) was higher (P < 0.01) with low [Hb] during submaximal exercise (i.e., at 30 W, LBF was 2.5 +/- 0.1 and 3.0 +/- 0.1 l/min for control [Hb] and low [Hb], respectively; P < 0.01), resulting in a maintained O2 delivery and O2 uptake for a given workload. However, at peak exercise, LBF was unaltered (6.5 +/- 0.4 and 6.6 +/- 0.6 l/min for control [Hb] and low [Hb], respectively), which resulted in an 18% reduction in O2 delivery (P < 0.01). This occurred despite peak cardiac output in neither condition reaching >75% of maximal cardiac output (approximately 26 l/min). It is concluded that a low CaO2 induces an elevation in submaximal muscle blood flow and that O2 delivery to contracting muscles is tightly regulated.

Hypoxia and the cardiovascular response to dynamic knee-extensor exercise

[EN] Hypoxia affects O2 transport and aerobic exercise capacity. In two previous studies, conflicting results have been reported regarding whether O2 delivery to the muscle is increased with hypoxia or whether there is a more efficient O2 extraction to allow for compensation of the decreased O2 availability at submaximal and maximal exercise. To reconcile this discrepancy, we measured limb blood flow (LBF), cardiac output, and O2 uptake during two-legged knee-extensor exercise in eight healthy young men. They completed studies at rest, at two submaximal workloads, and at peak effort under normoxia (inspired O2 fraction 0.21) and two levels of hypoxia (inspired O2 fractions 0.16 and 0.11). During submaximal exercise, LBF increased in hypoxia and compensated for the decrement in arterial O2 content. At peak effort, however, our subjects did not achieve a higher cardiac output or LBF. Thus O2 delivery was not maintained and peak power output and leg O2 uptake were reduced proportionately. These data are consistent then with the findings of an increased LBF to compensate for hypoxemia at submaximal exercise, but no such increase occurs at peak effort despite substantial cardiac capacity for an elevation in LBF.

A subpixel edge detector applied to aortic dissection detection

[EN] The aortic dissection is a disease that can cause a deadly situation, even with a correct treatment. It consists in a rupture of a layer of the aortic artery wall, causing a blood flow inside this rupture, called dissection. The aim of this paper is to contribute to its diagnosis, detecting the dissection edges inside the aorta. A subpixel accuracy edge detector based on the hypothesis of partial volume effect is used, where the intensity of an edge pixel is the sum of the contribution of each color weighted by its relative area inside the pixel. The method uses a floating window centred on the edge pixel and computes the edge features. The accuracy of our method is evaluated on synthetic images of different hickness and noise levels, obtaining an edge detection with a maximal mean error lower than 16 percent of a pixel.

Angiogenesis and angioregression gene expression analyses in swine corpus luteum

The corpus luteum (CL) lifespan is characterized by a rapid growth, differentiation and controlled regression of the luteal tissue, accompanied by an intense angiogenesis and angioregression. Indeed, the CL is one of the most highly vascularised tissue in the body with a proliferation rate of the endothelial cells 4- to 20-fold more intense than in some of the most malignant human tumours. This angiogenic process should be rigorously controlled to allow the repeated opportunities of fertilization. After a first period of rapid growth, the tissue becomes stably organized and prepares itself to switch to the phenotype required for its next apoptotic regression. In pregnant swine, the lifespan of the CLs must be extended to support embryonic and foetal development and vascularisation is necessary for the maintenance of luteal function. Among the molecules involved in the angiogenesis, Vascular Endothelial Growth Factor (VEGF) is the main regulator, promoting endothelial cells proliferation, differentiation and survival as well as vascular permeability and vessel lumen formation. During vascular invasion and apoptosis process, the remodelling of the extracellular matrix is essential for the correct evolution of the CL, particularly by the action of specific class of proteolytic enzymes known as matrix metalloproteinases (MMPs). Another important factor that plays a role in the processes of angiogenesis and angioregression during the CL formation and luteolysis is the isopeptide Endothelin-1 (ET-1), which is well-known to be a potent vasoconstrictor and mitogen for endothelial cells. The goal of the present thesis was to study the role and regulation of vascularisation in an adult vascular bed. For this purpose, using a precisely controlled in vivo model of swine CL development and regression, we determined the levels of expression of the members of VEGF system (VEGF total and specific isoforms; VEGF receptor-1, VEGFR-1; VEGF receptor-2, VEGFR-2) and ET- 1 system (ET-1; endothelin converting enzyme-1, ECE-1; endothelin receptor type A, ET-A) as well as the activity of the Ca++/Mg++-dependent endonucleases and gelatinases (MMP-2 and MMP-9). Three experiments were conducted to reach such objectives in CLs isolated from ovaries of cyclic, pregnant or fasted gilts. In the Experiment I, we evaluated the influence of acute fasting on VEGF production and VEGF, VEGFR-2, ET-1, ECE-1 and ET-A mRNA expressions in CLs collected on day 6 after ovulation (midluteal phase). The results indicated a down-regulation of VEGF, VEGFR-2, ET-1 and ECE-1 mRNA expression, although no change was observed for VEGF protein. Furthermore, we observed that fasting stimulated steroidogenesis by luteal cells. On the basis of the main effects of VEGF (stimulation of vessel growth and endothelial permeability) and ET-1 (stimulation of endothelial cell proliferation and vasoconstriction, as well as VEGF stimulation), we concluded that feed restriction possibly inhibited luteal vessel development. This could be, at least in part, compensated by a decrease of vasal tone due to a diminution of ET-1, thus ensuring an adequate blood flow and the production of steroids by the luteal cells. In the Experiment II, we investigated the relationship between VEGF, gelatinases and Ca++/Mg++-dependent endonucleases activities with the functional CL stage throughout the oestrous cycle and at pregnancy. The results demonstrated differential patterns of expression of those molecules in correspondence to the different phases of the oestrous cycle. Immediately after ovulation, VEGF mRNA/protein levels and MMP-9 activity are maximal. On days 5–14 after ovulation, VEGF expression and MMP-2 and -9 activities are at basal levels, while Ca++/Mg++-dependent endonuclease levels increased significantly in relation to day 1. Only at luteolysis (day 17), Ca++/Mg++-dependent endonuclease and MMP-2 spontaneous activity increased significantly. At pregnancy, high levels of MMP-9 and VEGF were observed. These results suggested that during the very early luteal phase, high MMPs activities coupled with high VEGF levels drive the tissue to an angiogenic phenotype, allowing CL growth under LH (Luteinising Hormone) stimulus, while during the late luteal phase, low VEGF and elevate MMPs levels may play a role in the apoptotic tissue and extracellular matrix remodelling during structural luteolysis. In the Experiment III, we described the expression patterns of all distinct VEGF isoforms throughout the oestrous cycle. Furthermore, the mRNA expression and protein levels of both VEGF receptors were also evaluated. Four novel VEGF isoforms (VEGF144, VEGF147, VEGF182, and VEGF164b) were found for the first time in swine and the seven identified isoforms presented four different patterns of expression. All isoforms showed their highest mRNA levels in newly formed CLs (day 1), followed by a decrease during mid-late luteal phase (days 10–17), except for VEGF182, VEGF188 and VEGF144 that showed a differential regulation during late luteal phase (day 14) or at luteolysis (day 17). VEGF protein levels paralleled the most expressed and secreted VEGF120 and VEGF164 isoforms. The VEGF receptors mRNAs showed a different pattern of expression in relation to their ligands, increasing between day 1 and 3 and gradually decreasing during the mid-late luteal phase. The differential regulation of some VEGF isoforms principally during the late luteal phase and luteolysis suggested a specific role of VEGF during tissue remodelling process that occurs either for CL maintenance in case of pregnancy or for noncapillary vessel development essential for tissue removal during structural luteolysis. In summary, our findings allow us to determine relationships among factors involved in the angiogenesis and angioregression mechanisms that take place during the formation and regression of the CL. Thus, CL provides a very interesting model for studying such factors in different fields of the basic research.

Mechanosensitivity in the myenteric plexus of the guinea pig ileum

The enteric nervous system regulates autonomously from the central nervous system all the reflex pathways that control blood flow, motility, water and electrolyte transport and acid secretion. The ability of the gut to function in isolation is one of the most intriguing phenomenons in neurogastroenterology. This requires coding of sensory stimuli by cells in the gut wall. Enteric neurons are prominent candidates to relay mechanosensitivity. Surprisingly, the identity of mechanosensitive neurons in the enteric nervous system as well as the appropriate stimulus modality is unknown despite the evidence that enteric neurons respond to sustained distension. Objectives: The aim of our study was to record from mechanosensitive neurons using physiological stimulus modalities. Identification of sensory neurons is of central importance to understand sensory transmission under normal conditions and in gut diseases associated with sensorimotor dysfunctions, such as Irritable Bowel Syndrome. Only then it will be possible to identify novel targets that help to normalise sensory functions. Methods: We used guinea-pig ileum myenteric plexus preparations and recorded responses of all neurons in a given ganglion with a fast neuroimaging technique based on voltage sensitive dyes. To evoke a mechanical response we used two different kinds of stimuli: firstly we applied a local mechanical distortion of the ganglion surface with von Frey hair. Secondarily we mimic the ganglia deformation during physiological movements of myenteric ganglia in a freely contracting ileal preparation. We were able to reliably and reproducibly mimic this distortion by intraganglionic injections of small volumes of oxygenated and buffered Krebs solution using stimulus parameters that correspond to single contractions. We also performed in every ganglion tested, electrical stimulations to evoke fast excitatory postsynaptic potentials. Immunohistochemistry reactions were done with antibodies against Calbindin and NeuN, considered markers for sensory neurons. Results: Recordings were performed in 46 ganglia from 31 guinea pigs. In every ganglion tested we found from 1 to 21 (from 3% to 62%) responding cells with a median value of 7 (24% of the total number of neurons). The response consisted of an almost instantaneous spike discharge that showed adaptation. The median value of the action potential frequency in the responding neurons was 2.0 Hz, with a recording time of 1255 ms. The spike discharge lasted for 302 ± 231 ms and occurred only during the initial deformation phase. During sustained deformation no spike discharge was observed. The response was reproducible and was a direct activation of the enteric neurons since it remained after synaptic blockade with hexamethonium or ω-conotoxin and after long time perfusion with capsaicin. Muscle tone appears not to be required for activation of mechanosensory neurons. Mechanosensory neurons showed a response to mechanical stimulation related to the stimulus strength. All mechanosensory neurons received fast synaptic inputs. There was no correlation between mechanosensitivity and Calbindin-IR and NeuN-IR (44% of mechanosensitive neurones Calb-IR-/NeuN-IR-). Conclusions: We identified mechanosensitive neurons in the myenteric plexus of the guinea pig ileum which responded to brief deformation. These cells appear to be rapidly accommodating neurons which respond to dynamic change. All mechanosensitive neurons received fast synaptic input suggesting that their activity can be highly modulated by other neurons and hence there is a low stimulus fidelity which allows adjusting the gain in a sensory network. Mechanosensitivity appears to be a common feature of many enteric neurons belonging to different functional classes. This supports the existence of multifunctional enteric neurons which may fulfil sensory, integrative and motor functions.

Development of Clinical Decision Support Systems based on Mathematical Models of Physiological Systems

In the last years of research, I focused my studies on different physiological problems. Together with my supervisors, I developed/improved different mathematical models in order to create valid tools useful for a better understanding of important clinical issues. The aim of all this work is to develop tools for learning and understanding cardiac and cerebrovascular physiology as well as pathology, generating research questions and developing clinical decision support systems useful for intensive care unit patients. I. ICP-model Designed for Medical Education We developed a comprehensive cerebral blood flow and intracranial pressure model to simulate and study the complex interactions in cerebrovascular dynamics caused by multiple simultaneous alterations, including normal and abnormal functional states of auto-regulation of the brain. Individual published equations (derived from prior animal and human studies) were implemented into a comprehensive simulation program. Included in the normal physiological modelling was: intracranial pressure, cerebral blood flow, blood pressure, and carbon dioxide (CO2) partial pressure. We also added external and pathological perturbations, such as head up position and intracranial haemorrhage. The model performed clinically realistically given inputs of published traumatized patients, and cases encountered by clinicians. The pulsatile nature of the output graphics was easy for clinicians to interpret. The manoeuvres simulated include changes of basic physiological inputs (e.g. blood pressure, central venous pressure, CO2 tension, head up position, and respiratory effects on vascular pressures) as well as pathological inputs (e.g. acute intracranial bleeding, and obstruction of cerebrospinal outflow). Based on the results, we believe the model would be useful to teach complex relationships of brain haemodynamics and study clinical research questions such as the optimal head-up position, the effects of intracranial haemorrhage on cerebral haemodynamics, as well as the best CO2 concentration to reach the optimal compromise between intracranial pressure and perfusion. We believe this model would be useful for both beginners and advanced learners. It could be used by practicing clinicians to model individual patients (entering the effects of needed clinical manipulations, and then running the model to test for optimal combinations of therapeutic manoeuvres). II. A Heterogeneous Cerebrovascular Mathematical Model Cerebrovascular pathologies are extremely complex, due to the multitude of factors acting simultaneously on cerebral haemodynamics. In this work, the mathematical model of cerebral haemodynamics and intracranial pressure dynamics, described in the point I, is extended to account for heterogeneity in cerebral blood flow. The model includes the Circle of Willis, six regional districts independently regulated by autoregulation and CO2 reactivity, distal cortical anastomoses, venous circulation, the cerebrospinal fluid circulation, and the intracranial pressure-volume relationship. Results agree with data in the literature and highlight the existence of a monotonic relationship between transient hyperemic response and the autoregulation gain. During unilateral internal carotid artery stenosis, local blood flow regulation is progressively lost in the ipsilateral territory with the presence of a steal phenomenon, while the anterior communicating artery plays the major role to redistribute the available blood flow. Conversely, distal collateral circulation plays a major role during unilateral occlusion of the middle cerebral artery. In conclusion, the model is able to reproduce several different pathological conditions characterized by heterogeneity in cerebrovascular haemodynamics and can not only explain generalized results in terms of physiological mechanisms involved, but also, by individualizing parameters, may represent a valuable tool to help with difficult clinical decisions. III. Effect of Cushing Response on Systemic Arterial Pressure. During cerebral hypoxic conditions, the sympathetic system causes an increase in arterial pressure (Cushing response), creating a link between the cerebral and the systemic circulation. This work investigates the complex relationships among cerebrovascular dynamics, intracranial pressure, Cushing response, and short-term systemic regulation, during plateau waves, by means of an original mathematical model. The model incorporates the pulsating heart, the pulmonary circulation and the systemic circulation, with an accurate description of the cerebral circulation and the intracranial pressure dynamics (same model as in the first paragraph). Various regulatory mechanisms are included: cerebral autoregulation, local blood flow control by oxygen (O2) and/or CO2 changes, sympathetic and vagal regulation of cardiovascular parameters by several reflex mechanisms (chemoreceptors, lung-stretch receptors, baroreceptors). The Cushing response has been described assuming a dramatic increase in sympathetic activity to vessels during a fall in brain O2 delivery. With this assumption, the model is able to simulate the cardiovascular effects experimentally observed when intracranial pressure is artificially elevated and maintained at constant level (arterial pressure increase and bradicardia). According to the model, these effects arise from the interaction between the Cushing response and the baroreflex response (secondary to arterial pressure increase). Then, patients with severe head injury have been simulated by reducing intracranial compliance and cerebrospinal fluid reabsorption. With these changes, oscillations with plateau waves developed. In these conditions, model results indicate that the Cushing response may have both positive effects, reducing the duration of the plateau phase via an increase in cerebral perfusion pressure, and negative effects, increasing the intracranial pressure plateau level, with a risk of greater compression of the cerebral vessels. This model may be of value to assist clinicians in finding the balance between clinical benefits of the Cushing response and its shortcomings. IV. Comprehensive Cardiopulmonary Simulation Model for the Analysis of Hypercapnic Respiratory Failure We developed a new comprehensive cardiopulmonary model that takes into account the mutual interactions between the cardiovascular and the respiratory systems along with their short-term regulatory mechanisms. The model includes the heart, systemic and pulmonary circulations, lung mechanics, gas exchange and transport equations, and cardio-ventilatory control. Results show good agreement with published patient data in case of normoxic and hyperoxic hypercapnia simulations. In particular, simulations predict a moderate increase in mean systemic arterial pressure and heart rate, with almost no change in cardiac output, paralleled by a relevant increase in minute ventilation, tidal volume and respiratory rate. The model can represent a valid tool for clinical practice and medical research, providing an alternative way to experience-based clinical decisions. In conclusion, models are not only capable of summarizing current knowledge, but also identifying missing knowledge. In the former case they can serve as training aids for teaching the operation of complex systems, especially if the model can be used to demonstrate the outcome of experiments. In the latter case they generate experiments to be performed to gather the missing data.

Outcome respiratorio a lungo termine nei soggetti affetti da cardiopatia congenita sottoposti ad intervento cardiochirurgico

Introduction: In the last years cardiac surgery for congenital heart disease (CHD) reduced dramatically mortality modifying prognosis, but, at the same time, increased morbidity in this patient population. Respiratory and cardiovascular systems are strictly anatomically and functionally connected, so that alterations of pulmonary hemodynamic conditions modify respiratory function. While very short-term alterations of respiratory mechanics after surgery were investigated by many authors, not as much works focused on long-term changes. In these subjects rest respiratory function may be limited by several factor: CHD itself (fetal pulmonary perfusion influences vascular and alveolar development), extracorporeal circulation (CEC), thoracotomy and/or sternotomy, rib and sternal contusions, pleural adhesions and pleural fibrosis, secondary to surgical injury. Moreover inflammatory cascade, triggered by CEC, can cause endothelial damage and compromise gas exchange. Aims: The project was conceived to 1) determine severity of respiratory functional impairement in different CHD undergone to surgical correction/palliation; 2) identify the most and the least CHD involved by pulmonary impairement; 3) find a correlation between a specific hemodynamic condition and functional anomaly, and 4) between rest respiratory function and cardiopulmonary exercise test. Materials and methods: We studied 113 subjects with CHD undergone to surgery, and distinguished by group in accord to pulmonary blood flow (group 0: 28 pts with normal pulmonary flow; group 1: 22 pts with increased flow; group 2: 43 pts with decreased flow; group 3: 20 pts with total cavo-pulmonary anastomosis-TCPC) followed by the Pediatric Cardiology and Cardiac Surgery Unit, and we compare them to 37 age- and sex-matched healthy subjects. In Pediatric Pulmonology Unit all pts performed respiratory function tests (static and dynamic volumes, flow/volume curve, airway resistances-raw- and conductance-gaw-, lung diffusion of CO-DLCO- and DLCO/alveolar volume), and CHD pts the same day had cardiopulmonary test. They all were examined and had allergological tests, and respiratory medical history. Results: restrictive pattern (measured on total lung capacity-TLC- and vital capacity-VC) was in all CHD groups, and up to 45% in group 2 and 3. Comparing all groups, we found a significant difference in TLC between healthy and group 2 (p=0.001) and 3 (p=0.004), and in VC between group 2 and healthy (p=0.001) and group 1(p=0.034). Inspiratory capacity (IC) was decreased in group 2 related to healthy (p<0.001) and group 1 (p=0.037). We showed a direct correlation between TLC and VC with age at surgery (p=0.01) and inverse with number of surgical interventions (p=0.03). Reduced FEV1/FVC ratio, Gaw and increased Raw were mostly present in group 3. DLCO was impaired in all groups, but up to 80% in group 3 and 50% in group 2; when corrected for alveolar volume (DLCO/VA) reduction persisted in group 3 (20%), 2 (6.2%) and 0 (7.1%). Exercise test was impaired in all groups: VO2max and VE markedly reduced in all but especially in group 3, and VE/VCO2 slope, marker of ventilatory response to exercise, is increased (<36) in 62.5% of group 3, where other pts had anyway value>32. Comparing group 3 and 2, the most involved categories, we found difference in VO2max and VE/VCO2 slope (respectively p=0.02 and p<0.0001). We evidenced correlation between rest and exercise tests, especially in group 0 (between VO2max and FVC, FEV1, VC, IC; inverse relation between VE/VCO2slope and FVC, FEV1 and VC), but also in group 1 (VO2max and IC), group 2 (VO2max and FVC and FEV1); never in group 3. Discussion: According with literature, we found a frequent impairment of rest pulmonary function in all groups, but especially in group 2 and 3. Restrictive pattern was the most frequent alteration probably due to compromised pulmonary (vascular and alveolar) development secondary to hypoperfusion in fetal and pre-surgery (and pre-TCPC)life. Parenchymal fibrosis, pleural adhesions and thoracic deformities can add further limitation, as showed by the correlation between group 3 and number of surgical intervention. Exercise tests were limited, particularly in group 3 (complex anatomy and lost of chronotropic response), and we found correlations between rest and exercise tests in all but group 3. We speculate that in this patients hemodynamic exceeds respiratory contribution, though markedly decreased.

Fluidodinamica biomedica computazionale del flusso sanguigno in presenza di affezioni strumentali e/o patologiche

La tesi di Dottorato studia il flusso sanguigno tramite un codice agli elementi finiti (COMSOL Multiphysics). Nell’arteria è presente un catetere Doppler (in posizione concentrica o decentrata rispetto all’asse di simmetria) o di stenosi di varia forma ed estensione. Le arterie sono solidi cilindrici rigidi, elastici o iperelastici. Le arterie hanno diametri di 6 mm, 5 mm, 4 mm e 2 mm. Il flusso ematico è in regime laminare stazionario e transitorio, ed il sangue è un fluido non-Newtoniano di Casson, modificato secondo la formulazione di Gonzales & Moraga. Le analisi numeriche sono realizzate in domini tridimensionali e bidimensionali, in quest’ultimo caso analizzando l’interazione fluido-strutturale. Nei casi tridimensionali, le arterie (simulazioni fluidodinamiche) sono infinitamente rigide: ricavato il campo di pressione si procede quindi all’analisi strutturale, per determinare le variazioni di sezione e la permanenza del disturbo sul flusso. La portata sanguigna è determinata nei casi tridimensionali con catetere individuando tre valori (massimo, minimo e medio); mentre per i casi 2D e tridimensionali con arterie stenotiche la legge di pressione riproduce l’impulso ematico. La mesh è triangolare (2D) o tetraedrica (3D), infittita alla parete ed a valle dell’ostacolo, per catturare le ricircolazioni. Alla tesi sono allegate due appendici, che studiano con codici CFD la trasmissione del calore in microcanali e l’ evaporazione di gocce d’acqua in sistemi non confinati. La fluidodinamica nei microcanali è analoga all’emodinamica nei capillari. Il metodo Euleriano-Lagrangiano (simulazioni dell’evaporazione) schematizza la natura mista del sangue. La parte inerente ai microcanali analizza il transitorio a seguito dell’applicazione di un flusso termico variabile nel tempo, variando velocità in ingresso e dimensioni del microcanale. L’indagine sull’evaporazione di gocce è un’analisi parametrica in 3D, che esamina il peso del singolo parametro (temperatura esterna, diametro iniziale, umidità relativa, velocità iniziale, coefficiente di diffusione) per individuare quello che influenza maggiormente il fenomeno.

Inverse problems in classical and quantum physics

The subject of this thesis is in the area of Applied Mathematics known as Inverse Problems. Inverse problems are those where a set of measured data is analysed in order to get as much information as possible on a model which is assumed to represent a system in the real world. We study two inverse problems in the fields of classical and quantum physics: QCD condensates from tau-decay data and the inverse conductivity problem. Despite a concentrated effort by physicists extending over many years, an understanding of QCD from first principles continues to be elusive. Fortunately, data continues to appear which provide a rather direct probe of the inner workings of the strong interactions. We use a functional method which allows us to extract within rather general assumptions phenomenological parameters of QCD (the condensates) from a comparison of the time-like experimental data with asymptotic space-like results from theory. The price to be paid for the generality of assumptions is relatively large errors in the values of the extracted parameters. Although we do not claim that our method is superior to other approaches, we hope that our results lend additional confidence to the numerical results obtained with the help of methods based on QCD sum rules. EIT is a technology developed to image the electrical conductivity distribution of a conductive medium. The technique works by performing simultaneous measurements of direct or alternating electric currents and voltages on the boundary of an object. These are the data used by an image reconstruction algorithm to determine the electrical conductivity distribution within the object. In this thesis, two approaches of EIT image reconstruction are proposed. The first is based on reformulating the inverse problem in terms of integral equations. This method uses only a single set of measurements for the reconstruction. The second approach is an algorithm based on linearisation which uses more then one set of measurements. A promising result is that one can qualitatively reconstruct the conductivity inside the cross-section of a human chest. Even though the human volunteer is neither two-dimensional nor circular, such reconstructions can be useful in medical applications: monitoring for lung problems such as accumulating fluid or a collapsed lung and noninvasive monitoring of heart function and blood flow.

Entwicklungen zur Absolutquantifizierung der Lungenperfusion mittels kontrastmittelverstärkter Magnetresonanztomographie

Die regionale Bestimmung der Durchblutung (Perfusion) ermöglicht differenzierte Aussagen über den Gesundheitszustand und die Funktionalität der Lunge. Durch neue Messverfahren ermöglicht die Magnetresonanztomographie (MRT) eine nicht-invasive und strahlungsfreie Untersuchung der Perfusion. Obwohl die Machbarkeit qualitativer MRT-Durchblutungsmessungen bereits gezeigt wurde, fehlt bisher eine validierte quantitative Methode. Ziel dieser Arbeit war eine Optimierung der bestehenden Messprotokolle und mathematischen Modelle zur Absolutquantifizierung der Lungenperfusion mit Magnetresonanztomographie. Weiterhin sollte die Methodik durch Vergleich mit einem etablierten Referenzverfahren validiert werden. Durch Simulationen und Phantommessungen konnten optimale MRT-Messparameter und ein standardisiertes Protokoll festgelegt werden. Des Weiteren wurde eine verallgemeinerte Bestimmung der Kontrastmittelkonzentration aus den gemessenen Signalintensitäten vorgestellt, diskutiert und durch Probandenmessungen validiert. Auf der Basis dieser Entwicklungen wurde die MRT-Durchblutungsmessung der Lunge tierexperimentell mit der Positronenemissionstomographie (PET) intraindividuell verglichen und validiert. Die Ergebnisse zeigten nur kleine Abweichungen und eine statistisch hochsignifikante, stark lineare Korrelation. Zusammenfassend war es durch die Entwicklungen der vorgestellten Arbeit möglich, die kontrastmittelgestützte MRT-Durchblutungsmessung der Lunge zu optimieren und erstmals zu validieren.

Lesion detection and classification in breast cancer: evaluation of approaches based on morphological features, tracer kinetic modelling and semi-quantitative parameters in MR functional imaging (DCE-MRI)

The diagnosis, grading and classification of tumours has benefited considerably from the development of DCE-MRI which is now essential to the adequate clinical management of many tumour types due to its capability in detecting active angiogenesis. Several strategies have been proposed for DCE-MRI evaluation. Visual inspection of contrast agent concentration curves vs time is a very simple yet operator dependent procedure, therefore more objective approaches have been developed in order to facilitate comparison between studies. In so called model free approaches, descriptive or heuristic information extracted from time series raw data have been used for tissue classification. The main issue concerning these schemes is that they have not a direct interpretation in terms of physiological properties of the tissues. On the other hand, model based investigations typically involve compartmental tracer kinetic modelling and pixel-by-pixel estimation of kinetic parameters via non-linear regression applied on region of interests opportunely selected by the physician. This approach has the advantage to provide parameters directly related to the pathophysiological properties of the tissue such as vessel permeability, local regional blood flow, extraction fraction, concentration gradient between plasma and extravascular-extracellular space. Anyway, nonlinear modelling is computational demanding and the accuracy of the estimates can be affected by the signal-to-noise ratio and by the initial solutions. The principal aim of this thesis is investigate the use of semi-quantitative and quantitative parameters for segmentation and classification of breast lesion. The objectives can be subdivided as follow: describe the principal techniques to evaluate time intensity curve in DCE-MRI with focus on kinetic model proposed in literature; to evaluate the influence in parametrization choice for a classic bi-compartmental kinetic models; to evaluate the performance of a method for simultaneous tracer kinetic modelling and pixel classification; to evaluate performance of machine learning techniques training for segmentation and classification of breast lesion.

CT perfusion image processing: analysis of liver tumors

Perfusion CT imaging of the liver has potential to improve evaluation of tumour angiogenesis. Quantitative parameters can be obtained applying mathematical models to Time Attenuation Curve (TAC). However, there are still some difficulties for an accurate quantification of perfusion parameters due, for example, to algorithms employed, to mathematical model, to patient’s weight and cardiac output and to the acquisition system. In this thesis, new parameters and alternative methodologies about liver perfusion CT are presented in order to investigate the cause of variability of this technique. Firstly analysis were made to assess the variability related to the mathematical model used to compute arterial Blood Flow (BFa) values. Results were obtained implementing algorithms based on “ maximum slope method” and “Dual input one compartment model” . Statistical analysis on simulated data demonstrated that the two methods are not interchangeable. Anyway slope method is always applicable in clinical context. Then variability related to TAC processing in the application of slope method is analyzed. Results compared with manual selection allow to identify the best automatic algorithm to compute BFa. The consistency of a Standardized Perfusion Index (SPV) was evaluated and a simplified calibration procedure was proposed. At the end the quantitative value of perfusion map was analyzed. ROI approach and map approach provide related values of BFa and this means that pixel by pixel algorithm give reliable quantitative results. Also in pixel by pixel approach slope method give better results. In conclusion the development of new automatic algorithms for a consistent computation of BFa and the analysis and definition of simplified technique to compute SPV parameter, represent an improvement in the field of liver perfusion CT analysis.