Abstract
Feasible estimations of perioperative changes in oxygen consumption could enable larger studies of its role in postoperative outcomes. Current methods, by pulmonary artery catheterisation or breathing gas analysis, are often regarded as either invasive or technically requiring. In this pilot study, we investigated the relationship between estimations of oxygen consumption, based on minimal-invasive cardiac output and arterial-central venous blood gas sampling, and indirect calorimetry in the perioperative period using the data collected during a clinical trial on perioperative oxygen transport.
In 20 patients >65 years during epidural and general anaesthesia for open abdominal surgery, Fick-based estimations of oxygen consumption(EVO2), the product of cardiac output from LiDCO™plus(LiDCO Ltd, Cambridge, UK) and arterial-central venous oxygen content difference, were compared with indirect calorimetry(GVO2) using QuarkRMR(COSMED, srl. Italy). Eighty-five simultaneous intra- and postoperative measurements at different time-points were analysed for prediction, parallelity and by traditional agreement assessment. There was an overall association between GVO2 and EVO2, 73(95% CI 62 to 83) + 0.45(95% CI 0.29 to 0.61) EVO2 ml min-1m-2, P<0.0001. GVO2 and EVO2 changed in parallel intra- and postoperatively when normalised to their respective overall means. Unadjusted mean difference between GVO2 and EVO2 indexed for body surface area was 26(95% CI 20 to 32) with limits of agreement (1.96SD) of -32 to 85 ml min-1m-2 and did not change over time. There was low correlation for absolute agreement, ICC(A,1) 0.37(95% CI 0.34 to 0.65) [F(84,10.2)=3.07, P=0.0266].
Despite lack of absolute agreement, the estimated oxygen consumption changed in parallel to the metabolic measurements in the perioperative period. Prediction or trending of oxygen consumption by this or similar methods could be further evaluated in larger samples.
Introduction
A postoperative imbalance between oxygen consumption and delivery, leading to increased oxygen extraction, has been associated with increased morbidity and mortality after major surgery.(1) The focus of goal-directed haemodynamic therapy (GDHT) has traditionally been on oxygen delivery, which is often easier to assess and to develop measurable optimisation strategies for.(2) Recently, interest is growing to reassess perioperative oxygen consumption in current surgical populations using modern monitoring and analytic methodologies.(3-5) Feasible estimations could enable larger studies on the role of oxygen consumption in postoperative outcomes. Available techniques, by pulmonary artery catheterisation or indirect calorimetry, are either deemed too invasive or difficult to manage in a clinical study setting during non-cardiac surgery. Using oxygen uptake calculated from fractions of inspiratory and expiratory oxygen in the closed breathing circuit during low-flow anaesthesia(6) has not demonstrated agreement when compared to standard methods.(7) Importantly, it can not be used in awake patients in the postoperative period. Commonly used haemodynamic monitoring in major surgery, such as minimal-invasive cardiac output with arterial and central venous access, could offer a possibility not only to estimate intra- and postoperative oxygen consumption but also to follow changes over time. By substituting mixed with central venous oxygen content and using the cardiac output derived from a minimal-invasive monitor, an estimation of oxygen consumption can be calculated by the reverse Fick principle.(8) The lack of absolute agreement between calorimetric and Fick-based methods has been reported previously, the latter do not include pulmonary oxygen consumption and global oxygen consumption values are usually reported around 20-40 ml min-1m-2 lower compared to those obtained from breathing gas analysis.(9-11) Yet, if this bias remains unchanged in the intra- and postoperative period, such estimations could be studied in larger samples and related to other clinical parameters and outcomes.
Our aim of this pilot study was to investigate the relationship and temporal changes between estimations of oxygen consumption, based on minimal-invasive cardiac output monitoring and arterial-central venous blood gas samples (EVO2), and measurements by indirect calorimetry (GVO2) in the perioperative period using the data collected during our oxygen transport study in elderly undergoing major abdominal surgery.(12)
Materials and methods
The present analysis was a secondary objective of a prospective observational study on perioperative oxygen transport in elderly patients undergoing major upper abdominal surgery (clinicaltrials.gov NCT03355118). The Regional Ethics Review Board of the Stockholm Region (ID 2017/291-31/4) approved the study and written informed consent was obtained from all participants. The primary aim of the study, i.e. the perioperative oxygen transport changes, has been published.(12) Data collection and analysis for the present publication were pre-planned and conducted simultaneously. The original cohort data was prospectively collected 2017-2018.
Patients and settings
A detailed description of the selection criteria, patient characteristics, perioperative management and oxygen transport outcomes can be found in the previous publication.(12) As stated there, 20 ASA II-IV patients over 65 years undergoing open pancreatic or liver resection surgery in epidural and general anaesthesia were included. The study was conducted at the Karolinska University Hospital in Huddinge, a tertiary referral center for upper abdominal surgery.
Data extraction and time-points
Paired values of oxygen consumption by estimations based on minimal-invasive cardiac output monitoring and arterial-central venous blood gas samples (EVO2) and indirect calorimetry GVO2) from five perioperative time-points were analysed; T1: during anaesthesia, right before surgical skin incision; T2: early during surgery, directly after skin incision; T3: later during surgery, >2h after skin incision; T4: early postoperatively, <12h after extubation; T5: late postoperatively, on postoperative day 1. The mean values for GVO2 during the approximate 20-minute measurement periods were compared with simultaneous cardiac output measurements averaged for each minute exported from LiDCOviewPRO (LiDCO Ltd, Cambridge, UK). The blood gas parameters were calculated as means of two simultaneously drawn arterial and central venous samples at 5 and 15 minutes into the measurement period.
Measurements of VO2 by indirect calorimetry (GVO2)
Indirect calorimetry was performed by QuarkRMR (COSMED srl, Italy). This device applies a breath-by-breath technique to measure gas flow and concentrations that are synchronised by data processing algorithms. The Haldane transformation is used to calculate oxygen consumption.(13) During intraoperative measurements, the flow meter (Flow-REE, COSMED srl, Italy), gas sampling line and moist filter were placed between the endotracheal tube and the Y-piece of the ventilator. The ventilator was set to a fresh gas flow of 2 L min-1 and FiO2 of 0.5 during measurements to allow for gas sampling. All other ventilation settings were left unchanged. Postoperative measurements were made with a tight-fitting face mask connected to a bidirectional turbine flow meter and a gas sampling line. No supplemental oxygen was administered during the postoperative measurements. The calorimeter was calibrated before start of intraoperative measurements and before each postoperative measurement after a warm-up time of 20 minutes with a standardised gas mixture containing 16% oxygen and 5% carbon dioxide. The gas sampling line, Flow-REE and moist filter were changed before each measurement (except before T2, continuous to T1) and all flowmeters were calibrated with a 3L-syringe.
Estimation of VO2 by minimal-invasive cardiac output and arterial-central venous oxygen content difference (EVO2)
EVO2 was calculated by the reverse Fick’s principle with central venous instead of pulmonary artery blood using the following formulas:(14) [CO; cardiac output in L min-1, Ca-cvO2; oxygen content difference between arterial and central venous blood in ml dl-1, Hb; haemoglobin in g dl-1, SaO2; arterial oxygen saturation, ScvO2; central venous saturation, PaO2; partial pressure of oxygen in arterial blood, PcvO2; partial pressure of oxygen in central venous blood, constants 1.31 and 0.0225, referring to the Hüfner constant and the solubility coefficient of oxygen (ml O2 dl-1 kPa-1), and 10 as a conversion factor from dL to L.]
Cardiac output was obtained from LiDCO™plus (LiDCO Ltd, Cambridge, UK). The device was calibrated and recalibrated according to the manufacturer’s instructions. Missing values from CO measurements (averaged for each minute) were substituted by linear interpolation if no more than 3 data points and no major haemodynamic changes occurred. Blood gases were analysed immediately after sampling by ABL800 Flex or ABL90 Flex (Radiometer Medical ApS, Denmark). Cardiac output and measured oxygen consumption were indexed for body surface area using the DuBois formula yielding values of GVO2 and EVO2 in ml min-1m-2.(15)
Statistical analysis
The sample size calculation was performed for the primary study, in which 20 patients were expected to demonstrate a relevant change in oxygen consumption after induction of anaesthesia. This would yield a maximum of 100 paired measurements of EVO2 and GVO2 which was considered sufficient even in the presence of >10% data loss. Continuous data was tested for normality distribution and statistical tests applied accordingly. Statistical analyses were performed and constructed in R (version 3.5.3; R Foundation for Statistical Computing, Vienna, Austria, URL; https://www.R-project.org) and SAS (version 9.4; SAS Institute Inc, Cary, NC, U.S.). The statisticians conducting the analyses were not involved in the data collection. Mean difference between EVO2 and GVO2 with 95% confidence interval were calculated from the individual paired measurements and grouped by time point (T1-5). These changes over time were analysed by linear mixed models with Holm-adjusted Tukey post-hoc tests. To investigate the overall association between EVO2 and GVO2, a random coefficient model was used based on individual slopes and coefficients. Analyses of the perioperative changes over time of GVO2 compared to EVO2 and its input variables (CI; cardiac index and Ca-cvO2) were conducted by random effect mixed models with method or component and time as fixed effects. Adjustment for differences in variances of the methods or components was made. In these models, the relative changes were normalised to the patients’ individual baseline measurements (T1). In the model analysing changes of each method in awake and anaesthetised subjects, the changes were normalised to the respective overall mean. Traditional agreement assessment was also performed by intraclass correlation and Bland-Altman analysis. Single score intraclass correlation was used, a in a two-way model yielding ICC coefficients with 95% CI. Bias and limits of agreement with 95% CI was visualised in Bland-Altman plots. Both ICC and Bland-Altman analyses were performed separately for each time-point T1-T5. The overall ICC and Bland-Altman analyses were not adjusted for repeated measurements as these were performed under varying intra- and postoperative conditions. Normality and homoscedasticity were assessed in residual plots. An alpha of 0.05 was considered significant.
Results
A total of 85 paired measurements of EVO2 and GVO2 were obtained in 20 subjects; 58 were obtained intraoperatively and 27 in the postoperative period. Four paired intraoperative measurements were not performed due to early termination of surgery (unexpected metastatic spread of malignancy) in two patients. Thirteen paired measurements could not be performed in the postoperative period because of technical or arterial line failure (n=2), logistical reasons (n=2), patients’ decline (n=3), exclusion due to short postoperative stay (n=4) and need for supplemental oxygen (n=2). Correct positioning of the CVC was confirmed by postoperative chest x-ray in all patients.
There was a overall mean difference between EVO2 and GVO2 [F(1, 167) = 72.8, P<0.0001] estimated to -26 (95% CI -20 to -32; P<0.001) ml min-1m-2. This difference was largely unchanged at the different time-points [F(4, 168)=1.39, P=0.241]. The means of GVO2 and EVO2 at the different perioperative time-points (T1-T5) are presented in Fig 1.
GVO2 and EVO2 changed in parallel when separated to the anaesthetised intraoperative state [F(2, 49.9)=0.57, P=0.5669] and the awake postoperative state F(1, 22) = 0.00, P=0.9604), see Figs 2a-b. An overall association between GVO2 and EVO2 was demonstrated in a random coefficient model (Fig 3). The two patients with early termination of surgery were excluded from this analysis. The variances of EVO2 and its components, oxygen content difference in arterial and central venous blood (Ca-cvO2) and cardiac index (CI) were larger compared to GVO2 at all time-points, these analyses are presented in Suppl S1.
The traditional agreement analyses are presented as supplementary results. The overall unadjusted mean bias was 26 ml min-1 m-2 with limits of agreement (1.96SD) of -32 to 85 ml min-1m-2. Excluding one outlier in the late postoperative period (a patient with a large Ca-cvO2 difference) changed the unadjusted bias to 28 (LoA -20 to 75) ml min-1 m-2. Bland-Altman plots were constructed to illustrate the bias and limits of agreement at the different time-points Suppl S2. The overall correlation for absolute agreement was poor, with an intraclass coefficient ICC(A,1) of 0.37 (95% CI 0.34 to 0.65) [F(84,10.2)=3.07, P=0.0266], and did not improve much when adjusted for lower overall mean difference of EVO2, ICC(A,1)=0.51 (95% CI 0.34 to 0.65) [F(84, 84) = 3.07, P<0.001]. Graphs depicting the correlation between indexed GVO2 and EVO2 at the different time-points (T1-5) including the unadjusted overall correlation are presented in Suppl S3.
Discussion
Estimates of increased oxygen extraction, i.e low mixed or central venous saturation, are associated with poor surgical outcomes.(16, 17) However, cut-off levels remain unclear and the quality of evidence is low.(18) In order to further study and distinguish the role of oxygen consumption in the perioperative period, feasible estimations are needed. To the best of our knowledge, this is the first pilot study investigating how a Fick-derived estimation method based on minimal-invasive haemodynamic monitoring (LiDCO™plus and blood gas sampling from arterial and central venous lines) can be used intra- and postoperatively. As expected, this estimation method had poor agreement concerning absolute values when compared to indirect calorimetry, but approximate to previous studies using pulmonary artery catheters.(11, 19-21) Importantly, the estimations based on routine haemodynamic monitoring for major surgery, were shown to change in parallel with the metabolic measurements in the perioperative period. We suggest that this or similar methods could be evaluated in larger samples and related to clinical outcomes.
Most previous studies investigating methods for oxygen consumption monitoring perioperatively or in critically ill patients were performed decades ago using traditional method comparison analytical methods. Some of the earlier method comparison studies are summarised in Table 1. Newer studies using non-invasive cardiac output monitors have not shown agreement with oxygen consumption measurements from indirect calorimetry(22) or pulmonary artery catheters.(23) However, the monitors used were not calibrated by transpulmonary or indicator dilution such as the PiCCO™ or LiDCO™plus systems and did not analyse changes over time.
A feasible oxygen consumption estimation method does not necessarily have to demonstrate perfect agreement to reflect changes in the perioperative period. Time effects and repeated measurements in the same subject under changing conditions are important statistical challenges in studies involving perioperative patients. Previous studies have often used simple linear regression or correlation(9, 10, 22, 24) or Bland-Altman analysis(25) without correction for repeated measurements(11, 22, 26, 27) except for some.(21, 28) Only a few addressed the relationship between measurements over time.(19, 29) In the present study, we developed a prediction model for EVO2 and GVO2 by using a random coefficient model based on individual slopes and intercepts. A significant positive association was demonstrated here, but such prediction models should obviously be evaluated in larger samples. We also present analyses of relative changes of EVO2 and its components with GVO2. The parallelity that was demonstrated could indicate an ability of EVO2 to track changes in oxygen consumption. To address this further, multiple measurements during shorter periods of time would be required. Analytic models previously used for cardiac output monitors such as polar plot approaches could be used to assess the magnitude and direction of changes.(30) Intraclass correlation (ICC) was used as it better reflects reliability and agreement based on analysis of variance of the pooled data.(31) When adjusted for the consistently lower values of EVO2, the ICC estimates of the model improved but not so much (ICC coefficient 0.51 vs 0.37). Bland-Altman analysis has since long been the standard method for visualisation of agreement when comparing different methods of VO2 monitoring.(21) Myles and Cui further elaborated the methodological issues related to repeated measurements in the same subject already considered by Bland and Altman(32) and proposed different random effects models to adjust limits of agreement.(21) As measurements were performed under varying perioperative conditions, we present the time-points separately and did not adjust the overall limits of agreement for repeated measurements in the same patient.
Oxygen consumption calculated by the reverse Fick equation is consistently reported lower than simultaneous measurements by analysis of respiratory gas exchange.(7, 9-11, 19, 24, 26, 28, 29, 33-36) This difference or bias has been attributed to the pulmonary oxygen consumption.(35, 37, 38) However, variability of Fick-derived measurements(34, 36)and wide limits of agreement(21) has made it difficult to estimate a systematic methodological bias. Many previous studies have either been performed in thoracic or cardiovascular surgery(9-11, 19, 37) or in critically ill patients.(24, 26, 36, 39) Pulmonary VO2 can be expected to increase after thoracic surgery(11) and in intensive care patients with varying degrees of lung injury.(40) Some studies that involve patients undergoing predominately abdominal surgery have shown acceptable agreement between the methods.(27, 41)The age of the studies is also reflected by the frequent use of the Deltatrac Metabolic Monitor® (Datex Instrumentarium, Helsinki, Finland), a metabolic monitor using a mixing-chamber technique and which is no longer in production. Many metabolic monitors in modern clinical use are based on breath-by-breath technology such as the Es-COVX® (GE Healthcare, Helsinki, Finland) or the QuarkRMR® in our study. Although there is supporting evidence for some overestimation of VO2, the technology has shown clinically acceptable agreement when compared with mixing-chamber methods(42, 43) and it has been validated in the semi-closed circle absorber systems commonly used in anaesthesia.(44) Our results on GVO2 were comparable with studies using Deltatrac II when corrected for difference in units (Table 1). The estimations of oxygen consumption rely on accurate cardiac output determinations and oxygen content difference measurements. The LiDCO™plus has shown acceptable performance against the pulmonary artery catheter and other devices in cross-comparisons in cardiac output accuracy studies. (45, 46) During rapidly changing haemodynamic situations, concerns regarding trending ability and underestimation of cardiac output have been raised.(30, 47) The 20-minute data extraction periods in this study were specifically chosen to represent perioperative time points that usually are without considerable circulatory instability. Central and mixed venous oxygen saturation have not shown interchangeability(48-50) but some studies have suggested that trends in ScvO2 can replace SvO2.(51-53) During stable intraoperative conditions, oxygen content difference is not expected to vary to a large extent whereas cardiac output can show considerable in- and between patient variability.(23) In our study, oxygen content difference and cardiac output demonstrated similar coefficients of variation.
The present study has several major limitations in addition to those discussed above. The sample size and number of observations are small, although in the same range as many earlier studies comparing calculated and measured VO2, and this could explain the large variability of EVO2. There was a considerable loss of data in the postoperative period limiting the conclusions on changes over time.
In summary, we have evaluated an oxygen consumption method that requires no extra equipment if perioperative CO-monitoring and CVC are used. Our results indicate that the performance could be equivalent to pulmonary artery catheters, and we have demonstrated parallelity in changes over time. Compared to earlier studies, we have used modern and commonly used calibrated CO-monitoring. We have also addressed the relationship between measurements over time using updated statistical analysis. The results should be regarded as indicative and further studies on larger samples are needed to establish if this method of estimating perioperative VO2 can prove useful.
Data Availability
All relevant data are within the manuscript and its Supporting Information files
Statements and declarations
The authors have no conflicts of interest to declare.
Supporting information
S1 File
S2 File
S3 File
Acknowledgements
The authors would like to thank research assistant nurses at the Clinical Research Unit, Dept. of Upper Abdominal Surgery, Karolinska University Hospital Huddinge, Birgitta Holmberg and Sirje Laur, for their help in recruiting and informing the patients included in this study.