Difference between arterial and end-tidal carbon dioxide and adverse events after non-cardiac surgery: a historical cohort study.

PURPOSE: The difference between arterial and end-tidal partial pressure of carbon dioxide (ΔCO2) is a measure of alveolar dead space, commonly evaluated intraoperatively. Given its relationship to ventilation and perfusion, ΔCO2 may provide prognostic information and guide clinical decisions. We hypothesized that higher ΔCO2 values are associated with occurrence of a composite outcome of re-intubation, postoperative mechanical ventilation, or 30-day mortality in patients undergoing non-cardiac surgery.
METHODS: We conducted a historical cohort study of adult patients undergoing non-cardiac surgery with an arterial line at a single tertiary care medical centre. The composite outcome, identified from electronic health records, was re-intubation, postoperative mechanical ventilation, or 30-day mortality. Student's t test and Chi-squared test were used for univariable analysis. Logistic regression was used for multivariable analysis of the relationship of ΔCO2 with the composite outcome.
RESULTS: A total of 19,425 patients were included in the final study population. Univariable analysis showed an association between higher mean (standard deviation [SD]) intraoperative ΔCO2 values and the composite outcome (6.1 [5.3] vs 5.7 [4.5] mm Hg; P = 0.002). After adjusting for baseline subject characteristics, every 5-mm Hg increase in the ΔCO2 was associated with a nearly 20% increased odds of the composite outcome (odds ratio, 1.20; 95% confidence interval, 1.12 to 1.28; P < 0.001).
CONCLUSIONS: In this patient population, increased intraoperative ΔCO2 was associated with an increased odds of the composite outcome of postoperative mechanical ventilation, re-intubation, or 30-day mortality that was independent of its relationship with pre-existing pulmonary disease. Future studies are needed to determine if ΔCO2 can be used to guide patient management and improve patient outcomes.

Postoperative respiratory complications are common and significant postoperative adverse events. Patients affected by complications such as pulmonary edema, postoperative respiratory failure, and pneumonia tend to have longer hospital stays, higher costs, and increased 30-day mortality.[1-4] Recent studies indicate that ventilator settings contribute to postoperative pulmonary complications, even in healthy patients.[5-8] These studies have focused primarily on the impact of inspiratory pressures and tidal volume (TV), with little data collected on ventilation efficiency. Dead space ventilation (VD) is affected by both lung function and perfusion,[9] and its intraoperative measurement may provide prognostic information and aid in optimizing ventilation management during surgery.

Dead space is the portion of TV that does not participate in gas exchange. Total dead space (VDphys) is composed of anatomic (VDana) and alveolar dead space (VDalv), represented by the equation VDphys = VDana + VDalv. The division between VDana and VDalv is not a strictly anatomic division, but depends on the interface between inspired air and alveolar gas. Dead space associated with conducting airways (endotracheal tube, trachea, proximal bronchi) is described as VDana and dead space associated with non-perfused alveoli (distal airways) is described as VDalv.[10] While VDana is relatively fixed in an intubated patient undergoing general anesthesia, VDalv can be influenced by perfusion (heart rate [HR], mean arterial pressure [MAP], cardiac output [CO], volume status, patient positioning), ventilation (respiratory rate [RR], TV, positive end-expiratory pressure [PEEP]), and metabolic (CO2 production)[11-13] variables, many of which are monitored and modified under general anesthesia.

Evaluation of VDphys is useful in assessment and management in some clinical settings. An elevated VDphys predicts an increased risk for mortality in acute respiratory disease syndrome (ARDS).[14-16] Through its relationship to atelectasis and ventilation/perfusion (V/Q) matching, VDphys is also useful for choosing the optimal level of PEEP in mechanically ventilated patients.[17-19] As VDana is often a relatively fixed component, evaluation and management based on VDphys or VDalv, might be useful in other patient populations undergoing mechanical ventilation.

There are multiple ways to calculate or estimate VDphys, VDana, and VDalv, but these are often not available in many clinical settings. The difference between arterial partial pressure of carbon dioxide (PaCO2) and partial pressure of end-tidal CO2 (PetCO2), ΔCO2, is one of the most readily available methods to estimate VDalv as the difference between these two values is attributed primarily to VDalv.[13,20-28] Increased ΔCO2 has been associated with poor outcomes in trauma[29] and critically ill patients requiring major surgery.[13]

The objective of this study was to evaluate if ΔCO2 is independently associated with worse outcomes in patients undergoing non-cardiac surgery. We hypothesized that, in patients undergoing general anesthesia, greater ΔCO2 is associated with an increased odds of the postoperative composite outcome of postoperative mechanical ventilation, re-intubation, or 30-day mortality.

Methods

Following approval from the Institutional Review Board, we conducted a historical cohort study at University of Michigan Medical Center (Ann Arbor, MI, USA); the requirement for patient consent was waived. We queried the electronic anesthesia records for all adult cases from 23 July 2009 to 31 January 2019 for which general endotracheal anesthesia was performed and at least one arterial blood gas (ABG) was recorded during the procedure. Exclusion criteria were an American Society of Anesthesiology (ASA) Physical Status classification of VI, preoperative mechanical ventilation, sodium bicarbonate administration within 15 min prior to ABG measurement, cardiopulmonary bypass or extracorporeal membrane oxygenation, absent PetCO2 or PaCO2 data, or intraoperative one-lung ventilation.

Patient and perioperative data

Demographic information was collected from the preoperative history and physical documentation. This comprised age, sex, height, weight (actual and ideal), body mass index (BMI), ASA Physical Status score, smoking status (current or former), chronic obstructive pulmonary disease (COPD), and asthma. Preoperative values for albumin, blood urea nitrogen, creatinine, and peripheral oxygen saturation (SpO2) were obtained. Creatinine clearance was calculated using the Cockcroft–Gault formula. Procedure variables included (i) type, (ii) duration, and (iii) emergent vs elective procedure. Other intraoperative variables were the total transfusion requirement including volume of packed red blood cells, fresh frozen plasma, platelets, and cryoprecipitate; estimated blood loss; total crystalloid administration; and grams of albumin administered.

Blood gas, ventilator, and hemodynamic data

For all patients, the first ABG obtained in the operating room after induction of general anesthesia was utilized in our analysis and pH, PaCO2, arterial partial pressure of oxygen (PaO2), and arterial oxygen saturation (SaO2) were collected. Ventilator parameters (TV, RR, minute ventilation, pressure support, peak airway pressure, PEEP, and fraction of inspired oxygen), hemodynamic data (HR and MAP), and PetCO2 were collected for the ten-minute epoch surrounding the ABG time and the median values were utilized in the analysis.

Outcome

Our study outcome was a composite of the postoperative complications, postoperative mechanical ventilation, re-intubation, or mortality, as identified from the electronic health record (EHR). Postoperative mechanical ventilation was defined as any patient who required mechanical ventilation after the procedure for any duration not associated with another procedure. We defined re-intubation as a period of time in the EHR without receipt of mechanical ventilation subsequently followed by the documented presence of mechanical ventilation. Mortality was included as an outcome since it would act to censor future pulmonary complications. Patients were assumed to be alive at 30 days unless noted to be deceased in the medical record.

Statistics

We calculated descriptive statistics for each variable. Continuous variables that were normally distributed were presented as mean (standard deviation [SD]). Continuous variables that were non-normally distributed were summarized as median [interquartile range (IQR)]. Categorical variables were summarized using counts and proportions.

Preoperative characteristics and ΔCO2 were compared between those with and without the composite outcome using Chi square tests for categorical variables and t tests or Mann–Whitney U tests for continuous variables. Multivariable logistic regression models were generated to evaluate whether ΔCO2 was associated with the composite outcome while adjusting for other patient and surgical factors based on clinical and statistical significance. For entry into the model, we chose clinically relevant variables that were selected by two anesthesiologists based on the literature. Additionally, we considered variables with standardized differences greater than 0.2 using Cohen’s d. We also collapsed several variables into one variable, or chose a representative variable from a group, when the groups of variables were similar or highly correlated (ρ > 0.8). Upon constructing the multivariable model, we removed, collapsed, or grouped variables with variance inflation factors > 2.5.

For each of the continuous variables included in the model (age, creatinine clearance, preoperative albumin, total albumin, case duration, total volume blood products, percent inspired oxygen, TV, RR, MAP, HR, and ΔCO2), the following potential transformations were considered: X-2, X-1, X-0.5, ln(X), X0.5, X1, and X2. The continuous variables were binned and plotted as the proportion within each bin with the outcome. Variables with values of zero were shifted by one or one plus their minimum to use all the listed transformations. Univariate regressions with the composite outcome as the outcome and each transformation of each continuous variable were then used to predict values at the median of the bins and plot the predicted outcome for the transformation. These transformations were inspected visually and used to determine that no non-linear terms were necessary.

Additionally, all variables that were significant predictors in the initial multivariable model were also considered for interactions with the primary predictor, ΔCO2. Each significantly associated variable was centred and an interaction between it and ΔCO2 was added individually to the initial full regression. All interactions that were significant individually were used together in a full regression. Interactions that were not significant at a level of P < 0.05 in this new full regression were removed. Even though multiple interactions remaining in the model had P < 0.05, most interaction terms had a small magnitude of association with the outcome, none changed the direction of the relationship between ΔCO2 and the primary outcome, and none of the variables were conceptually thought to have a multiplicative effect for any of the variables with the ΔCO2. Based on this, none of the interactions were added to the final, clinically relevant model.

Akaike Information Criteria, Bayes Information Criteria, -2 log likelihood, and area under the receiver operating characteristic curve (95% confidence interval [CI]) were used to evaluate the model. P values < 0.05 and CIs excluding 1 were considered statistically significant. Analyses were performed in R version 3.6.2 (R Core Team, R Foundation for Statistical Computing, Vienna, Austria).

Results

After excluding patients with missing data or data concerning for errant manual entry, 19,425 patients remained for analysis (Fig. 1). The mean (SD) time from anesthesia start to ABG was 73 (26) min. In our patient population, 15% (n = 2,830) had COPD, 10% (n = 1,951) had asthma, and 7% (n = 1,274) were on supplemental oxygen preoperatively (Table 1). The mean values for PaCO2 and PetCO2 were 40.2 (5.8) mm Hg and 34.3 (4.0) mm Hg, respectively, with a mean ΔCO2 of 5.8 (4.6) mm Hg. A total of 3,006 patients (16%) had the composite outcome: 85 (0.4%) required postoperative re-intubation, 2,769 (14%) required postoperative mechanical ventilation, and 452 (2%) had 30-day mortality.

Fig. 1

Exclusion criteria: Using initial query parameters, information for 39,002 patients were included. Patients were excluded based on pre-determined exclusion criteria. After initial exclusions, box plots of patient data were reviewed, and outliers were eliminated if data were determined to be entered in error (venous blood gas charted as arterial) or unlikely based on known physiologic variables. Ultimately, 19,425 patients were included in the final analysis.

Table 1

Patient and demographic characteristics

		Total N	n	(%)
Male		19,355	10,981	(56.7%)
Comorbidities
Supplemental oxygen Current or former smoker COPD Asthma Dialysis		19,355	1,274	(6.6%)
		19,355	10,669	(55.2%)
		19,355	2,830	(14.6%)
		19,335	1,951	(10.1%)
		19,354	554	(2.9%)
ASA Physical Status score III, IV, or V		19,334	15,627	(81%)
BMI	Normal	18,724	4,905	(26.2%)
	Underweight	18,724	498	(2.7%)
	Overweight	18,724	5,848	(31.2%)
	Obesity class I	18,724	3,943	(21.1%)
	Obesity class II	18,724	1,979	(10.6%)
	Obesity class III	18,724	1,551	(8.3%)
Emergent		19,334	2,624	(13.6%)
Postoperative re-intubation		19,355	85	(0.4%)
30-day mortality		19,355	452	(2.3%)
Postoperative ventilator		19,355	2,769	(14.3%)
Procedure type
General surgery Gynecologic surgery		19,355 19,355	3,145 329	(16.2%) (1.7%)
Head and neck surgery		19,355	2,339	(12.1%)
Neurosurgery Orthopedic surgery Plastic surgery Thoracic surgery Transplant surgery Urologic surgery Vascular surgery Other surgery		19,355 19,355 19,355 19,355 19,355 19,355 19,355 19,355	4,676 1,441 231 800 1,099 1,546 2,315 1,434	(24.2%) (7.4%) (1.2%) (4.1%) (5.7%) (8.0%) (12.0%) (7.4%)

Body mass index (BMI) was subdivided into World Health Organization categories: underweight (BMI < 18.5 kg·m-2), normal (BMI 18.5–24.9 kg·m-2) overweight (BMI 25–29.9 kg·m-2), obesity class I (BMI 30–34.9 kg·m-2), obesity class II (BMI 35.0–39.9 kg·m-2), and obesity class III (BMI ≥ 40 kg·m-2). Tidal volume calculated by dividing tidal volume by ideal body weight.

*The median [interquartile range] time of arterial blood gas from anesthesia start was 71 [54–91] min. †As most patients did not receive transfusions or colloid, numbers reflect only those patients who received a particular product with the corresponding means (standard deviations)/medians [interquartile ranges].

ASA = American Society of Anesthesiologists; BMI = body mass index; COPD = chronic obstructive pulmonary disease; FFP = fresh frozen plasma; IQR = interquartile range; PaCO2 = arterial partial pressure of carbon dioxide; PaO2 = arterial partial pressure of oxygen; PEEP = positive end expiratory pressure; PetCO2 = end-tidal partial pressure of carbon dioxide; pRBC = packed red blood cells; SD = standard deviation; SpO2 = peripheral oxygen saturation.

Case duration (min)

Time of ABG from anesthesia start (min)*

19,354

357

73

(171)

(26)

Preoperative and intraoperative variables and outcomes

Many of the evaluated variables were significantly different between groups of patients who did and did not experience the composite outcome (Table 2). Notably, patients who had the composite outcome were more likely to use preoperative supplemental oxygen (19% vs 4%; P < 0.001) and have COPD (17% vs 14%; P < 0.001), and less likely to have asthma (8% vs 10%; P < 0.001). Similar smoking rates (55% vs 55%; P = 0.386) were present in each cohort. Patients who had the composite outcome had a significantly higher mean (SD) ΔCO2 (6.1 [5.3] vs 5.7 [4.5] mm Hg; P = 0.002).

Table 2

Univariable analysis for study variables

		No composite outcome			Composite outcome			P value	Cohen’s d
		n	Total N	(%)	n	Total N	(%)
Male		9,184	16,349	(56.2%)	1,797	3,006	(59.8%)	0.001	0.07
Preoperative supplemental O2		712	16,349	(4.4%)	562	3,006	(18.7%)	< 0.001	0.46
Comorbidities:
Current or former smoker		9,035	16,349	(55.3%)	1,634	3,006	(54.4%)	0.36	0.02
Diagnosis of COPD		2,333	16,349	(14.3%)	497	3,006	(16.5%)	0.002	0.06
Diagnosis of asthma		1,708	16,349	(10.4%)	243	3,006	(8.1%)	< 0.001	0.08
Preoperative dialysis		355	16,349	(2.2%)	199	3,005	(6.6%)	< 0.001	0.22
BMI	Normal	4,154	15,926	(26.1%)	751	2,798	(26.8%)	< 0.001	0.02
	Underweight	387	15,926	(2.4%)	111	2,798	(4.0%)		0.09
	Overweight	5,031	15,926	(31.6%)	817	2,798	(29.2%)		0.05
	Obesity class I	3,396	15,926	(21.3%)	547	2,798	(19.5%)		0.04
	Obesity class II	1,681	15,926	(10.6%)	298	2,798	(10.7%)		< 0.01
	Obesity class III	1,277	15,926	(8.0%)	274	2,798	(9.8%)		0.06
ASA Physical Status score III, IV, or V		12,946	16,335	(79.3%)	2,681	2,999	(89.4%)	< 0.001	0.28
Emergent		1,295	16,335	(7.9%)	1,329	2,999	(44.3%)	< 0.001	0.91
Procedure type
General surgery		2,470	16,349	(15.1%)	675	3,006	(2.5%)	< 0.001	0.19
Gynecologic surgery		284	16,349	(1.7%)	45	3,006	(1.5%)	0.39	0.02
Head and neck surgery		1,944	16,349	(11.9%)	395	3,006	(13.1%)	0.06	0.04
Neurosurgery		4,267	16,349	(26.1%)	409	3,006	(13.6%)	< 0.001	0.32
Orthopedic surgery		1,294	16,349	(7.9%)	147	3,006	(4.9%)	< 0.001	0.12
Plastic surgery		201	16,349	(1.2%)	30	3,006	(1.0%)	0.32	0.02
Thoracic surgery		726	16,349	(4.4%)	74	3,006	(2.5%)	< 0.001	0.11
Transplant surgery		626	16,349	(3.8%)	473	3,006	(15.7%)	< 0.001	0.41
Urologic surgery		1,432	16,349	(8.8%)	114	3,006	(3.8%)	< 0.001	0.21
Vascular surgery		2,020	16,349	(12.4%)	295	3,006	(9.8%)	< 0.001	0.08
Other surgery		1,085	16,349	(6.6%)	349	3,006	(11.6%)	< 0.001	0.17

Body mass index (BMI) was subdivided into World Health Organization categories: underweight (BMI < 18.5 kg·m-2), normal (BMI 18.5–24.9 kg·m-2), overweight (BMI 25–29.9 kg·m-2), obesity class I (BMI 30–34.9 kg·m-2), obesity class II (BMI 35.0–39.9 kg·m-2), and obesity class III (BMI ≥ 40 kg·m-2).

*t test P value unless otherwise indicated; † Median values from the 10 min epoch surrounding arterial blood gas analysis were utilized in calculation; ‡ Mann–Whitney U test P value. § As most patients did not receive transfusions or colloid, numbers reflect only those patients who received a particular product with the corresponding means (standard deviations)/medians [interquartile ranges].

ABG = arterial blood gas; ASA = American Society of Anesthesiologist ; BMI = body mass index; BUN = blood urea nitrogen; COPD = chronic obstructive pulmonary disease; Cr = creatinine; FFP = fresh frozen plasma; IQR = interquartile range; PaCO2 = arterial partial pressure of carbon dioxide; PaO2 = arterial partial pressure of oxygen; PEEP = positive end expiratory pressure; PetCO2 = end-tidal partial pressure of carbon dioxide; pRBC = packed red blood cells; SD = standard deviation; SpO2 = peripheral oxygen saturation.

Multivariable regressions

After adjusting for other factors, ΔCO2 was associated with an increased occurrence of re-intubation, postoperative mechanical ventilation, or 30-day mortality (Table 3). In our regression, every 5-mm Hg increase in the ΔCO2 was associated with a nearly 20% increased odds of the composite outcome (odds ratio [OR], 1.20; 95% CI, 1.12 to 1.28; P < 0.001). For our multivariable regression, at ΔCO2 = 0, the median [IQR] risk of re-intubation, postoperative mechanical ventilation, or 30-day mortality was 9.2 [7.9–10.5]%, at ΔCO2 = 5 mm Hg, the risk was 10.8 [9.5–12.2]%, and at ΔCO2 = 10 mm Hg, the risk was 12.7 [11.0–14.4]%. For the rare patient with ΔCO2 = 20 mm Hg, the risk was 17.2 [13.8–20.6]% (Fig. 2).

Table 3

Multivariable regression for composite outcome

	Odds ratio	(95% CI)	P value
(Intercept)	0.01	(0.004 to 0.024)	< 0.001
Preoperative variables
Emergent procedure	5.45	(4.69 to 6.33)	< 0.001
Age (yr)	1.01	(1.00 to 1.01)	0.007
Male	1.03	(0.90 to 1.17)	0.71
Supplemental O2	2.37	(1.97 to 2.84)	< 0.001
Current or former smoker	0.91	(0.81 to 1.03)	0.15
COPD	1.16	(0.99 to 1.36)	0.06
ASA Physical Status III–IV or V	1.21	(1.01 to 1.45)	0.04
Asthma	0.75	(0.61 to 0.93)	0.007
Preoperative dialysis	1.45	(1.08 to 1.95)	0.01
BMI category
Normal	1
Underweight	1.03	(0.74 to 1.44)	0.85
Overweight	1.01	(0.87 to 1.18)	0.89
Obesity class I	1.14	(0.96 to 1.35)	0.15
Obesity class II	1.12	(0.90 to 1.39)	0.32
Obesity class III	1.34	(1.06 to 1.71)	0.02
Procedure type
Head and neck + plastic surgery	1.32	(1.09 to 1.60)	0.004
Neurosurgery	0.57	(0.48 to 0.68)	< 0.001
Orthopedic surgery	0.58	(0.45 to 0.75)	< 0.001
Urologic surgery	0.62	(0.48 to 0.80)	< 0.001
General/gynecologic/thoracic/transplant surgery	0.69	(0.57 to 0.85)	< 0.001
Other surgery	1.56	(1.21 to 2.03)	0.001
Creatinine clearance (mL·min-1)	1.00	(1.00 to1.00)	0.001
Preoperative albumin (g·dL-1)	0.53	(0.48 to 0.58)	< 0.001
Intraoperative variables
Albumin transfused (g)	1.01	(1.01 to 1.01)	< 0.001
Case duration (min)	1.29	(1.24 to 1.34)	< 0.001
Total volume blood products (mL)	1.81	(1.68 to 1.94)	< 0.001
Inspired oxygen (%)	1.01	(1.01 to 1.02)	< 0.001
Tidal volume (mL·kg-1)	1.05	(1.00 to 1.09)	0.046
Respiratory rate (breaths·min-1)	1.07	(1.05 to 1.10)	< 0.001
Mean arterial pressure (mm Hg)	1.0	(0.99 to 1.00)	0.243
Heart rate (beats·min-1)	1.02	(1.01 to 1.02)	< 0.001
PaCO2 - PetCO2 (5 mm Hg)	1.20	(1.12 to 1.28)	< 0.001

BMI was subdivided into World Health Organization categories: underweight (BMI < 18.5 kg·m-2), normal (BMI 18.5–24.9 kg·m-2), overweight (BMI 25–29.9 kg·m-2), obesity class I (BMI 30–34.9 kg2), obesity class II (BMI 35.0–39.9 kg·m-2), and obesity class III (BMI ≥ 40 kg·m-2). Akaike Information Criteria = 8356.444, Bayes Information Criteria = 8599.658, -2logLikelihood = 928292.444, and C statistic = 0.87 (95% CI, 0.86 to 0.88).

ASA = American Society of Anesthesiologists; BMI = body mass index; CI = confidence interval; COPD = chronic obstructive pulmonary disease; PaCO2 = arterial partial pressure of carbon dioxide; PetCO2 = end-tidal partial pressure of carbon dioxide.

Fig. 2

Predicted risk of pulmonary complications for different ΔCO2 values overlaid on a histogram showing the distribution of ΔCO2 values included in this study. The ΔCO2 range spanned values of -10 mm Hg to 30 mm Hg with most being between 0 mm Hg and 10 mm Hg. After adjusting for all factors in the multivariable model, a curvilinear relationship can be seen between ΔCO2 values and postoperative pulmonary complications with higher ΔCO2 values predicting greater risk. ΔCO2 = difference between arterial and end-tidal partial pressure of carbon dioxide.

Several other variables were independently associated with the composite outcome (Table 3). The use of preoperative supplemental oxygen was associated with the composite outcome (OR, 2.37; 95% CI, 1.97 to 2.84; P < 0.001), while higher albumin levels were associated with a decrease in the incidence of the composite outcome (OR, 0.53; 95% CI, 0.48 to 0.58; P < 0.001 per g·dL-1). Chronic obstructive pulmonary disease was not associated with re-intubation, postoperative mechanical ventilation, or 30-day mortality (OR, 1.16; 95% CI, 0.99 to 1.36; P = 0.06). The model had good discrimination (C statistic, 0.87; 95% CI, 0.86 to 0.88).

Discussion

In this historical cohort study, we found that increased ΔCO2 was associated with an increase in the composite outcome of postoperative mechanical ventilation, re-intubation, or 30-day mortality in patients undergoing general anesthesia for a non-cardiac surgical procedure. This relationship remained significant even when we adjusted for other relevant pre-existing clinical factors, including COPD, preoperative oxygen use, dialysis, and age. Since the magnitude of association for a 5 mm Hg of change is nearly 20%, even relatively small changes in ΔCO2 produce clinically important changes to the risk of the composite outcome. For a patient with a ΔCO2 of 2–4 mm Hg, consistent with normal healthy levels, the median [IQR] predicted risk of pulmonary complications or death in this patient population was between 9.9 [9.6–10.1]% and 10.5 [10.1–11.2]%. For a ΔCO2 of 10, 15, and 20 mm Hg, the risk increased to 12.7 [ 11.0–14.4]%, 14.8 [ 12.9–17.1]%, and 17.2 [13.8–20.6]%, respectively.

The magnitude of ΔCO2 is affected by both ventilated but unperfused alveoli (true dead space) and alveoli with high ventilation/perfusion ratios (V/Q inequality). Shunt (perfusion without ventilation) will cause a small increase in ΔCO2 by raising the PaCO2 as mixed venous blood, which is higher in PCO2, crosses the pulmonary circulation and mixes with the systemic arterial circulation. True dead space and V/Q inequality are markers of disease severity and increase with COPD severity or low CO states. ΔCO2 is also influenced by dynamic, modifiable, and clinically important variables of ventilation (RR, TV, PEEP), perfusion (HR, MAP, CO, volume status, patient positioning), and metabolism (CO2 production).[11-13,30] Some of these are under the control of the anesthesiologist, and studies should be conducted to determine how anesthesiologists’ manipulation of these parameters affects both ΔCO2, and more importantly, patient outcome.

Clinically, ΔCO2 has been utilized as a surrogate for VDalv, primarily in the intensive care unit, where levels correlate with the amount of VDalv.[31] In addition, ΔCO2 has been shown to correlate with disease severity in ARDS patients, with a higher ΔCO2 being associated with more severe ARDS[32] and with death.[33] Similarly, elevated ΔCO2 in patients who have survived cardiac arrest is associated with high hospital mortality.[34] Nevertheless, use of ΔCO2 as an intraoperative prognostic measure is less common but is supported by our findings.

Tyburski et al. studied 501 trauma patients undergoing emergency surgery. In patients who died (n = 147 [29%]), ΔCO2 was higher at all three time points evaluated (initial, post-resuscitation, and final). In those who died, lower values were predictive of longer survival with a final mean (SD) ΔCO2 of 22.25 (14.32) mm Hg in patients who died intraoperatively, 19.96 (14.91) mm Hg in patients who died within 24 hr after surgery, and 10.90 (10.13) mm Hg in patients who died after 24 hr but before hospital discharge. Survivors had a ΔCO2 of 7.37 (6.30) mm Hg (all P < 0.006 compared with non-survivors).[29] In another study of critically ill or injured patients requiring major surgery with 41% mortality, Domsky et al. found highest mortality rates in patient with the highest ΔCO2, and highest estimated VDalv fraction.[13]Our study is in agreement with these two studies. While these studies are limited to trauma and critically ill surgery patients, our study is more generalizable as it includes a large variety of non-cardiac operations. To our knowledge, our study is the first large study showing that increased ΔCO2 is associated with our composite outcome in a broad patient population of adults undergoing general anesthesia with arterial lines. Our results indicate that high intraoperative ΔCO2 should not be ignored in patients undergoing general anesthesia. Additional studies are necessary to better understand the clinical utility of our results.

There are a few important limitations to this study. As a single-centre study, it might not be generalizable to other centres with different intraoperative ventilator strategies or practice surrounding PaCO2 measurement indications or timing. At our institution, 97.7% of patients who had an arterial line placed for their procedure had an ABG measured with a mean (SD) time from anesthesia start to ABG measurement of 73 (26) min. This limits our findings to patients and procedures for which arterial lines are utilized and likely explains the wider ΔCO2 found in this study compared with previous studies that excluded ASA III–V patients or patients with “respiratory or cardiac abnormalities.”[35,36] Second, PaCO2 and PetCO2 are dynamic as patients undergo hours of anesthesia. The first ABG was selected to facilitate comparability of patients undergoing anesthesia as a baseline since many patients had only one ABG at the beginning of the case. Selecting the first ABG allowed for analysis of a large patient population; it also minimized bias. Nevertheless, this measure was more representative of the initial clinical condition, rather than the final clinical condition, of patients undergoing anesthesia. Furthermore, we did not evaluate whether ΔCO2 changes over the course of surgery and how that is associated with outcomes. Evaluating these changes and whether they can predict outcome needs further investigation.

Previous studies have investigated risk factors for postoperative respiratory failure and have developed clinically predictive models.[37] Some of these clinically significant variables were not included in our analysis because they were not collected in our EHR. For example, New York Heart Association heart failure class has been shown to be a risk factor for postoperative pulmonary complications,[37] but this variable is not present in our anesthesia preoperative history. Where a previous study found that lower preoperative SpO2 levels were associated with postoperative pulmonary complications,[38] our patients with low values had been placed on supplemental oxygen. Additionally, risk scores developed on a more general surgical population may not be applicable to our unique patient population, specifically patients with an arterial catheter who underwent non-cardiac surgery.

Despite these limitations, our findings support the need for future research on the relationship between ΔCO2 and postoperative adverse events. This study was retrospective and must be evaluated prospectively in patients undergoing non-cardiac surgery with a study design to address some of the limitations encountered in this patient population. Next, additional patient populations, such as those undergoing anesthesia without mechanical ventilation, pediatric patients, or those undergoing one-lung ventilation should be evaluated. Patients who require one-lung ventilation might represent a specific population risk of ventilation efficiency and where ΔCO2 might have clinical predictive value.

In addition to evaluating the ΔCO2 in other patient populations, additional research is needed to determine if specific therapies are beneficial when elevated ΔCO2 is detected. For example, one animal study of ARDS found that the PEEP value that minimized capnography-measured VDalv also maximized CO and mixed venous oxygen saturation.[39] Further study is needed to determine if individualizing care, through ventilator adjustments or manipulation of hemodynamics, based on minimizing ΔCO2 reduces adverse events.

Conclusions

Our study shows an independent association between increases in ΔCO2 and a composite outcome of re-intubation, postoperative mechanical ventilation, or 30-day mortality in patients undergoing non-cardiac surgery. Future research is needed to determine if the ΔCO2 can be used for risk stratification or as a target for the optimization of intraoperative ventilator and hemodynamic management.