Infection in the intensive care unit alters physiological networks
© Grossman et al; licensee BioMed Central Ltd. 2009
Published: 17 September 2009
Physicians use clinical and physiological data to treat patients every day, and it is essential for treating a patient appropriately. However, medical sources of clinical physiological data are only now starting to find use in bioinformatics research.
We collected 29 types of physiological and clinical data on a minute-by-minute basis from trauma patients in the intensive care unit along with whether they contracted an infection during their stay. Dividing the patients into two groups based on this criterion, we determined that the correlational network amongst pairs of physiological variables changes based on whether the patient contracted an infection.
Examining the variable pairs with the largest change in correlation across groups reveals potential changes in the way our treatments affect the patient's physiology and in how our bodies react to physiological insults. These findings highlight the usefulness of physiological informatics and suggest new relationships to study while also validating previously reported relationships.
While physicians have been making use of physiologic data for decades, these measurements have not always been collected and stored electronically for research use. Largely due to the advent of electronic medical records and improved methods of extracting data from patient monitoring equipment and paper charts, we can now obtain significant amounts of clinical data both during and after treatment [1–4]. This points to a relatively new source of data that has been made available to both clinicians and informatics researchers in ways that were not possible previously.
One of the main thrusts of informatics research using clinical data has been to provide greater context for defining disease using physiological and treatment data [5, 6]. Use of clinical data for basic biological research has been successful in linking physiological measurements to already known genetic markers of disease , and has yielded biomarkers related to human aging .
One area where there seems to be less research is in the physiological modelling of humans undergoing intensive care in a hospital setting and the various conditions these patients encounter. It is this gap that we seek to fill by providing data showing that even basic physiological models of patients in intensive care are changed under conditions of infection. We also show how treatments administered during intensive care can affect physiology in different ways depending on the state of the patient.
Our data were obtained from 19 patients who were admitted to the trauma intensive care unit (ICU) at San Francisco General Hospital from May 2004 to June 2005. These patients were generally heavily injured, with average stays in the ICU of 21 days and an average hospital stay of 36 days total. Upon admission to the ICU, patients were connected to standard ICU monitors and ventilators as necessary, and were also monitored with an experimental muscle microdialysis catheter capable of measuring metabolic parameters such as glucose, lactate, and pyruvate.
Types of data collected and their abbreviations
Arterial O2/CO2 partial pressure
Mean arterial blood pressure
Muscle oxygen level (temperature)
Oxygen saturation percentage
Fraction of inspired oxygen
Muscle lactate concentration
Muscle glucose concentration
Muscle glutamate concentration
Muscle pyruvate concentration
Muscle lactate/pyruvate ratio
Bedside glucose reading
Mechanical lung compliance
Positive end expiratory pressure
Volume of air per minute
Central venous pressure
Blood urea nitrogen
Prior to data analysis, the muscle microdialysis data required correction to reflect its semi-continuous sample collection method. As the microdialysis catheter would continuously collect extracellular fluid from the muscle for approximately 30 minutes between readings, each reading was an average of the metabolites collected over that span. To fill in the missing data from this data source, we performed linear interpolation between successive readings of the dialysate composition.
Data were analyzed by calculating Pearson correlation coefficients for each pair of variables under different physiological conditions (406 pairs in total). We grouped each patient according to whether he or she contracted an infection during his or her stay in the hospital. In our cohort 11/19 patients (58%) contracted an infection. The data used to calculate the correlation for each variable pair were all those rows containing data for both variables of that pair, i.e. "pairwise complete observations". While this can yield a correlation matrix that is not positive definite, we are concerned with the coefficients themselves rather than manipulating the resulting matrix.
Correlations were considered significant if the p-value derived from Student's t-test was less than 0.05 after correcting for multiple comparisons using the Holm-Bonferroni method. Taking those correlations that were significant in both the infection and non-infection groups and intersecting them resulted in the candidate set of variable pairs to further investigate. We then take the absolute value of the difference between the candidate correlations in our set and only consider correlations after applying a cut off value for the difference of 0.4.
Correlation coefficients of the variables with the largest change in correlation between conditions
ρ with infection
ρ without infection
We then consider the correlation pairs for which each correlation was significant. Four of the top five changes in correlation involved a switch in the direction of correlation in the two conditions – i.e. they are positively correlated in one condition and negatively correlated in the other.
Upon splitting the physiological data into patients who did and did not contract infections during their hospital stay, we have revealed that several physiological variables change the direction of their correlations. Here, we will explain these changes via physiology and general standards and methods of treatment in the ICU. It should also be noted that we are trying to explain relationships that do not indicate a direction of causality. As a tool, however, we will hypothesize causality in some of the following discussion.
Fraction of inspired O2 vs. muscle O2 content
In the non-infection case, as the fraction of oxygen inspired increases so too does the resulting partial pressure of oxygen throughout the remainder of the body's tissues, as expected. We would further expect the partial pressure to saturate very quickly with increasing oxygen inspiration, as in healthy individuals oxygen saturation is above 95% .
In the case of infection, though, there are many potential physiological changes, ranging from hypotension to impaired oxygen transfer, and other problems that could also lead to hypoperfusion. Sair et al. showed that in a rat model of endotoxemia the effect of increasing FiO2 on muscle oxygen saturation is significantly reduced in the infection group . Because there were no observable differences in tissue perfusion between the two groups in that experiment, it is possible that our data showing a reversal in the correlation are due to reduced tissue perfusion combined with the reduced oxygen delivery capability despite normal perfusion as shown by Sair.
Ikossi et al. showed that lower PmO2 was strongly associated with worse outcomes in the ICU, including infection and multiple organ failure . In patients who appeared to be well resuscitated by normal measures (MAP, heart rate, BD, and PaO2) reduced muscle oxygen content predicted worse outcomes and increased complications. Because this study used oxygen sensors placed in the deltoid muscle, this indicates a likely decrease in perfusion to distal organs using muscle as a proxy. This is despite the clinicians attempting to increase perfusion by increasing FiO2.
It is also relevant that inspired oxygen fraction is controlled by the ICU staff. Fraction of inspired O2 is primarily increased to compensate for decreases in tissue oxygen saturation. So the prior expectation is that patients in the ICU, especially those with impaired lung function (as those with pneumonia often have), will have their oxygen intake increased as their tissue oxygen pressure decreases, leading to the negative correlation seen in our data. The most likely explanation for the change seen in the two conditions is that the infection reduces the coupling of inspired oxygen fraction to oxygen present in tissue, and then clinical intervention results in increasing FiO2 in the face of decreasing PmO2.
Muscle pyruvate vs. muscle O2 content
In both the infection and non-infection groups, muscle pyruvate and muscle oxygen pressure are anti-correlated. In the non-infection group, the correlation is stronger, indicating a tighter connection between oxygen delivered to the tissue and the removal of pyruvate form the tissue via aerobic respiration. This could indicate consequences of infection similar to those proposed in the previous section – i.e. patients with infections have reduced ability to deliver oxygen to the tissue and so the rate of metabolism is less tightly coupled to the delivery of oxygen than in trauma patients without infection. This may point to mitochondrial dysfunction as a culprit in this relationship.
Minute volume vs. muscle glucose
It has been reported that glucose metabolism is altered in patients with sepsis, and furthermore that in critically ill patients reduced glucose can be seen largely as a marker for disease severity . We can therefore postulate the following. Patients with infections are clearly more ill than those without infections, suggesting that they will have reduced levels of glucose in their tissue. That these patients are sicker may also indicate a need for increased ventilation, either to expedite removal of CO2 from the blood or increase the ability of oxygen to enter the bloodstream. So in patients with infections, we would therefore expect to see the minute volume to be negatively correlated with tissue glucose.
This result may also not be explainable by traditional physiology. Glucose is tightly controlled by ICU physicians, so these results are also clouded by that intervention. We have shown here a relationship that does not yet exist in the literature, so this defines a new physiological relationship between these two variables. The difference in the relationship between these variables, then, may be an indicator of severity of sickness or a predictor of outcome rather than something to be explained.
Lung compliance vs. muscle pyruvate and PEEP
Lung compliance is a variable measured by mechanical ventilators. It is the ratio of the change in lung volume in response to a change in air pressure. Several factors affect lung compliance, which include the amount of surfactant present in the alveoli and pathologies, such as fibrosis.
In patients with infection, PEEP and compliance co-vary, while the opposite is true for non-infection patients. It has been shown that in rats given a dose of lipopolysaccharide (an antigen produced by the that bacteria cause sepsis) who also received increased constant airway pressure had increased lung damage . Experience dictates that sicker patients (like those with pneumonia) have stiffer and less functional lungs, which would necessitate the use of a higher PEEP to recruit otherwise collapsed alveoli. This increases the apparent lung volume for a given pressure in the sicker patients resulting in a measured increase in lung compliance. It is unclear why this should be the reverse in less sick patients, but in looking at the scatter plot it is clear that the majority of the data have a PEEP of approximately 5, with fewer samples at higher values.
Again, there is no simple reason why lung compliance should be at all related to the concentration of pyruvate in muscle tissue. It is not a problem that this relationship is not readily explainable, since this is now a new physiological relationship that we've defined. To be complete, however, a hypothesis is as follows. Lung compliance measures, in part, the ability of the alveoli to open and accept new air in response to increasing inspiratory pressure. In patients with infection, as compliance increases, muscle pyruvate decreases. This results in a greater volume of air being available to supply oxygen, thereby increasing the potential delivery rate of oxygen to the blood and peripheral tissues. Increased oxygen allows aerobic respiration to proceed faster, decreasing the amount of pyruvate present in the system.
Application as a decision tool
It should be noted that these results could find potential use as the basis for a decision tool to warn clinicians of impending infection. While our data set did not allow us to test this application, future work could include testing our results on novel data. This highlights the overall utility of pursuing physiological informatics and modelling using clinical data.
We have shown that physiological networks can be constructed from clinical measurements from an intensive care unit, and the topology of these networks can change as a patient contracts an infection. Between these two conditions many variables changed the strength or direction of their correlation. The five pairs of variables with the largest magnitude of change were:
FiO2 vs. muscle oxygen content
Muscle pyruvate vs. muscle oxygen content
Minute volume vs. muscle glucose
Lung compliance vs. PEEP
Lung compliance vs. muscle pyruvate
We presented potential physiological interpretations for these results as well. All but one of these variable pairs directly involve parameters for treatment/ventilation chosen by physicians in the ICU. Therefore, we can also conclude that the effects of treatment on physiological measurements changes when infections are contracted in the intensive care unit. Some of these differences have already been reported in the literature while others are novel. This shows that our technique can be used to discover previously unknown relationships between physiologic variables. This work points to opportunities to study the changes that have not yet been reported, potentially opening new doors to discover how our best efforts to heal patients can both alter and be informed by their physiology.
Supported by the Lucile Packard Foundation for Children's Health (AJB), and the Stanford Bio-X Graduate Fellowship (ADG), NIH R01 NS050173 (GTM), and UC Discovery Grant GCP06 10226 (GTM). The authors wish to thank Diane Morabito for her expertise in manually curating and cleaning these data.
This article has been published as part of BMC Bioinformatics Volume 10 Supplement 9, 2009: Proceedings of the 2009 AMIA Summit on Translational Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/10?issue=S9.
- Butte AJ: The ultimate model organism. Science 2008, 320: 325–327. 10.1126/science.1158343PubMed CentralView ArticlePubMedGoogle Scholar
- Chen L, McKenna TM, Reisner AT, Gribok A, Reifman J: Decision tool for the early diagnosis of trauma patient hypovolemia. Journal of biomedical informatics 2008, 41: 469–478. 10.1016/j.jbi.2007.12.002View ArticlePubMedGoogle Scholar
- Zhang Y, Szolovits P: Patient-specific learning in real time for adaptive monitoring in critical care. Journal of biomedical informatics 2008, 41: 452–460. 10.1016/j.jbi.2008.03.011PubMed CentralView ArticlePubMedGoogle Scholar
- Sorani MD, Hemphill JC 3rd, Morabito D, Rosenthal G, Manley GT: New approaches to physiological informatics in neurocritical care. Neurocrit Care 2007, 7: 45–52. 10.1007/s12028-007-0043-7View ArticlePubMedGoogle Scholar
- Sung NS, Crowley WF Jr, Genel M, Salber P, Sandy L, Sherwood LM, Johnson SB, Catanese V, Tilson H, Getz K, et al.: Central challenges facing the national clinical research enterprise. JAMA 2003, 289: 1278–1287. 10.1001/jama.289.10.1278View ArticlePubMedGoogle Scholar
- Payne PR, Johnson SB, Starren JB, Tilson HH, Dowdy D: Breaking the translational barriers: the value of integrating biomedical informatics and translational research. J Investig Med 2005, 53: 192–200. 10.2310/6650.2005.00402View ArticlePubMedGoogle Scholar
- Stoll M, Cowley AW Jr, Tonellato PJ, Greene AS, Kaldunski ML, Roman RJ, Dumas P, Schork NJ, Wang Z, Jacob HJ: A genomic-systems biology map for cardiovascular function. Science 2001, 294: 1723–1726. 10.1126/science.1062117View ArticlePubMedGoogle Scholar
- Chen DP, Weber SC, Constantinou PS, Ferris TA, Lowe HJ, Butte AJ: Novel integration of hospital electronic medical records and gene expression measurements to identify genetic markers of maturation. Pac Symp Biocomput 2008, 243–254.Google Scholar
- Despopoulos A, Silbernagl S: Color Atlas of Physiology. 5th edition. New York: Thieme; 2003.Google Scholar
- Sair M, Etherington PJ, Curzen NP, Winlove CP, Evans TW: Tissue oxygenation and perfusion in endotoxemia. Am J Physiol 1996, 271: H1620–1625.PubMedGoogle Scholar
- Ikossi DG, Knudson MM, Morabito DJ, Cohen MJ, Wan JJ, Khaw L, Stewart CJ, Hemphill C, Manley GT: Continuous muscle tissue oxygenation in critically injured patients: a prospective observational study. J Trauma 2006, 61: 780–788. discussion 788–790. discussion 788–790. 10.1097/01.ta.0000239500.71419.58View ArticlePubMedGoogle Scholar
- Van Cromphaut SJ, Vanhorebeek I, Berghe G: Glucose metabolism and insulin resistance in sepsis. Curr Pharm Des 2008, 14: 1887–1899. 10.2174/138161208784980563View ArticlePubMedGoogle Scholar
- Tsuchida S, Engelberts D, Roth M, McKerlie C, Post M, Kavanagh BP: Continuous positive airway pressure causes lung injury in a model of sepsis. Am J Physiol Lung Cell Mol Physiol 2005, 289: L554–564. 10.1152/ajplung.00143.2005View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.