Assessment of need for intravascular volume resuscitation remains challenging for anesthesiologists. Dynamic waveform indices, including systolic and pulse pressure variation, are demonstrated as reliable measures of fluid responsiveness for mechanically ventilated patients. Despite widespread use, real-world reference distributions for systolic and pulse pressure variation values have not been established for euvolemic intraoperative patients. The authors sought to establish systolic and pulse pressure variation reference distributions and assess the impact of modifying factors.
The authors evaluated adult patients undergoing general anesthetics for elective noncardiac surgery. Median systolic and pulse pressure variations during a 50-min postinduction period were noted for each case. Modifying factors including body mass index, age, ventilator settings, positioning, and hemodynamic management were studied via univariate and multivariable analyses. For systolic pressure variation values, effects of data entry method (manually entered vs. automated recorded) were similarly studied.
Among 1,791 cases, per-case median systolic and pulse pressure variation values formed nonparametric distributions. For each distribution, median values, interquartile ranges, and reference intervals (2.5th to 97.5th percentile) were, respectively, noted: these included manually entered systolic pressure variation (6.0, 5.0 to 7.0, and 3.0 to 11.0 mmHg), automated systolic pressure variation (4.7, 3.9 to 6.0, and 2.2 to 10.4 mmHg), and automated pulse pressure variation (7.0, 5.0 to 9.0, and 2.0 to 16.0%). Nonsupine positioning and preoperative β blocker were independently associated with altered systolic and pulse pressure variations, whereas ventilator tidal volume more than 8 ml/kg ideal body weight and peak inspiratory pressure more than 16 cm H2O demonstrated independent associations for systolic pressure variation only.
This study establishes real-world systolic and pulse pressure variation reference distributions absent in the current literature. Through a consideration of reference distributions and modifying factors, the authors’ study provides further evidence for assessing intraoperative volume status and fluid management therapies.
Systolic and pulse pressure variations are well-studied measures of fluid responsiveness in mechanically ventilated patients
However, reference distributions and clinical modifying factors for these parameters in adult patients undergoing elective noncardiac surgery are poorly described
Systolic pressure variation (SPV) and pulse pressure variation (PPV) reference distributions were established. Nonsupine positioning and preoperative β blocker were independently associated with altered SPV and PPV, whereas ventilator tidal volume more than 8 ml/kg ideal body weight and peak inspiratory pressure more than 16 cm H2O demonstrated independent associations for SPV only
ASSESSMENT of a patient’s intravascular volume status is a vital component guiding an anesthesiologist’s perioperative fluid management. During the past decade, studies have demonstrated an association between fluid administration strategy and postoperative outcome.1,2 However, accurate assessment of intravascular volume remains challenging.
Many clinical signs of hypovolemia (e.g., skin turgor, capillary refill, nausea, and syncope) are unreliable and become masked during general anesthesia, historically leaving anesthesiologists with blood pressure, heart rate, and urine output as indicators of volume status.3,4 Static measures—including central venous pressure, pulmonary artery occlusion pressure, and “noninvasive” cardiac output—have been shown inadequate to guide volume resuscitation therapy.5–7 While echocardiography including inferior vena cava imaging may be useful,8 associated equipment and training costs pose barriers to widespread use. Recently, dynamic indices—derived from the respirophasic variation of arterial pressure waveforms—have been developed and are suggested for use in mechanically ventilated patients. These include systolic pressure variation (SPV),9,10 pulse pressure (PP) variation (PPV),11 and stroke volume variation (SVV).12
SPV and PPV have been studied in the perioperative setting to assess fluid responsiveness or the benefit of a fluid challenge. Compared to SVV, SPV and PPV have demonstrated similar or superior ability to predict fluid responsiveness without the need for specialized SVV measurement equipment13,14 and eclipse other static measures.5,15 Limitations of SPV and PPV include requiring a regular heart rhythm16,17 and standardized ventilator settings,17,18 as well as sensitivity to confounding factors including positioning,19 intrathoracic pressure,16 intra abdominal pressure,16,20 and heart failure.21 Despite limitations, SPV and PPV have endured as commonly used measures of fluid responsiveness, owing to literature support, clinical safety, lack of other cost-effective, reliable measures, and automatic near-continuous calculation by newer physiologic monitors.
Current studies examining SPV and PPV are confined to subpopulations with limited sample sizes. Meta-analyses and multicenter studies have been performed in attempts to overcome this limitation, including a “gray zone” analysis of PPV among cardiac, vascular, and abdominal surgeries.13,22 Through these studies, theoretical ranges of SPV and PPV values have been proposed for responders and nonresponders to volume expansion. In spite of advances, there is a persistent lack of data establishing baseline distributions of SPV and PPV values among a broad population of patients in a euvolemic state undergoing elective surgical procedures. Additionally, although univariate analyses have examined changes in SPV and PPV associated with specific factors—including patient positioning,23,24 obesity,25 ventilator settings,26,27 and hemodynamic management28 —studies continue to be limited in clinical scope and size and have been frequently performed outside the perioperative setting. Furthermore, no multivariable analysis has evaluated independent associations among these factors.
To characterize the distribution of SPV and PPV values during elective surgical procedures, we performed a retrospective observational study at an academic tertiary care hospital. We hypothesized that for adult patients undergoing a wide range of elective noncardiac surgical procedures, distributions of median SPV and PPV values can be characterized for the euvolemic state. Additionally, we hypothesized that changes in distributions are associated with perioperative clinical factors including age, body mass index (BMI), positioning, preoperative medications, ventilator settings, and hemodynamic management.
Materials and Methods
This retrospective observational study was approved by the Institutional Review Board (HUM00052066; Ann Arbor, Michigan). Patient consent was waived. Deidentified data were extracted from our local single-center Multicenter Perioperative Outcomes Group database, used for storage of our institution’s electronic health record (EHR; Centricity® General Electric Healthcare, USA). Methodology for observational data collection, storage, and quality assurance within the Multicenter Perioperative Outcomes Group is described elsewhere.29,30 Per departmental policy, a detailed study protocol, including patient population, primary outcome, and planned statistical analyses, was presented at our local anesthesia clinical research committee on September 23, 2015, and registered on the department’s internal research website.
At our hospital, physiologic monitoring is acquired via automated and validated interfaces at each anesthetizing location (CARESCAPE B850 or Solar 9500®; General Electric Healthcare). For cases including an arterial blood pressure monitor, SPV and PPV measurements are automatically recorded into the EHR from the CARESCAPE B850 monitor; anesthesiologists can also choose to manually calculate and record SPV values as a distinct EHR database concept. Both manually entered and automated recorded values were calculated as the difference between the maximum and minimum systolic arterial pressures (SAPmax and SAPmin, respectively) observed on the arterial waveform throughout phases of the respiratory cycle, expressed in mmHg31 :
PPV values, expressed as a percentage, were calculated in an automated fashion using the established formula:
PPmax and PPmin represent the maximum and minimum PPs observed on the arterial waveform throughout phases of the respiratory cycle.31
All adult patients (greater than or equal to 18 yr) undergoing elective surgical procedures with general anesthesia, tracheal intubation, and invasive arterial blood pressure monitoring from January 1, 2009, to June 10, 2016, were included for analysis. To minimize the effect of prolonged preoperative fasting on intravascular volume status, cases were restricted to first-case morning start times, determined by documentation of anesthetic induction end between 7:30 and 8:30 AM. Additional exclusion criteria for the study were (1) cardiac or cardiothoracic operating room procedures, (2) open-chest procedures, (3) laparoscopic procedures, (4) procedures utilizing a double-lumen endotracheal tube or one-lung ventilation, (5) patients receiving regional or spinal anesthetics to supplement general anesthesia pre-induction, (6) inpatient admission status before the date of procedure, (7) patients with American Society of Anesthesiologists physical status class 5 and 6, (8) blood product transfusion during SPV/PPV measurement, (9) vasopressor or inotropic infusion during SPV/PPV measurement, and (10) patients with significant cardiopulmonary comorbidities. Cardiopulmonary comorbidity exclusion criteria were heart failure, chronic obstructive pulmonary disease, acute respiratory distress syndrome, sleep apnea, low functional capacity, or dysrhythmia other than sinus tachycardia/bradycardia, occasional premature ventricular or atrial contractions, or first-degree heart block. Inclusion/exclusion criteria are illustrated in figure 1; cardiopulmonary comorbidity diagnoses were established by the pick-list variables within the preoperative history and physical electronic documentation described in the appendix.
To account for a shift in measurement practices at our institution during the study period, patients were stratified by data entry method: those with manually entered SPV values only (PPV data unavailable) and those with automated recorded SPV and PPV values. This shift in practice was attributable to updated physiologic monitors capable of automated SPV/PPV recording (CARESCAPE B850). For cases in which manually entered and automated recorded SPV values were both available, automated recorded values were used. Validation of arterial blood pressure values—and subsequently derived SPV/PPV values—was performed by a manual review of a simple random sample of 125 patient intraoperative records. Through this manual review, an algorithm for removing possible artifactual arterial line values (e.g., dampened waveform and arterial line flushing for blood draws) was developed and subsequently applied to all cases. These included values for which the systolic blood pressure (SBP) was more than 200 mmHg and PP less than 50 mmHg, SBP more than 150 mmHg and PP less than 30 mmHg, SBP more than or equal to 100 mmHg and PP less than 20 mmHg, or SBP less than 100 mmHg and PP less than 10 mmHg. Automated SPV and PPV values concurrently recorded at the time of possible arterial line artifact were also removed.
For all cases, SPV and PPV measurements were obtained starting 10 min after anesthetic induction end and concluding 60 min after anesthetic induction end, or at surgical procedure end if earlier. This measurement period was selected to minimize the effects of confounders on SPV and PPV values (e.g., anesthetic induction agents before measurement period and blood loss, insensible loss, and fluid shifts after the measurement period). For automated recorded values, measurements were obtained once per minute; for manually entered values, measurements were obtained as frequently as recorded by the anesthesiologist—most commonly once every 15 min, per institutional practice. For each case meeting inclusion criteria, median SPV values (and PPV values when available) were calculated within the measurement period. To derive a standard distribution of arterial pressure variation values, histograms were developed for SPV values (manually entered and automated recorded) and PPV values.
After characterizing SPV and PPV distributions, associations between clinical factors and altered SPV and PPV measurements were analyzed. These factors included patient age and BMI, preoperative medications (preoperative angiotensin-converting enzyme inhibitor/angiotensin receptor blocker, β blocker, or diuretic), ventilator settings, decreased mean arterial pressure (MAP) more than 20% below baseline, and total fluid administration at the end of SPV/PPV measurement period. For purposes of statistical analysis, continuous study variables were collapsed to binary variables. BMI, tidal volume, peak inspiratory pressure (PIP), positive end-expiratory pressure (PEEP), and MAP cutoffs were chosen as used in previous literature.12,26,32–35 A patient age cutoff was chosen based upon the median age of the study population. Fluid administration was converted to crystalloid equivalents (colloids multiplied by a factor of 1.5),36 with a 1,000-ml binary cutoff. MAP values were obtained concurrently for all SPV/PPV measurements; a majority of MAP values more than 20% below initial intraoperative baseline defined hypotension as a binary variable for each case.
Analyses were performed using SAS version 9.3 (SAS Institute, USA) and SPSS version 21.0 (IBM, USA). Basic descriptive statistics were calculated for demographic and relevant intraoperative data. Pearson chi-square or Fisher exact tests (for categorical variables) and independent two-tailed Student’s t tests or Mann–Whitney U tests (for continuous variables) were used to assess baseline univariate clinical differences between patients with manually entered SPV values versus those with automated recorded SPV (and PPV) values (table 1). Each distribution of median SPV/PPV values was tested against the normal, lognormal, γ, and Weibull distributions using the Kolmogorov–Smirnov goodness-of-fit test within the SAS univariate procedure. A P < 0.05 denoted statistical significance.
Univariate clinical differences between median SPV/PPV value distributions were next assessed using Mann–Whitney U tests for all modifying factors studied. To assess for independent predictors of SPV and PPV values, separate full-fit multivariable linear regression models were derived for three study cohorts: cases with manually entered SPV values, cases with automated recorded SPV values, and cases with automated recorded PPV values. Before developing a prediction model to determine independent risk factors, all variables were tested for collinearity by investigating the correlations. Pairwise Pearson and Spearman correlation matrices were constructed to determine high correlation between variables, with a correlation coefficient of 0.70 as the threshold for high collinearity. If the correlation coefficient was less than 0.70, then no collinearity was detected, and all variables were eligible for multivariable model entry. Residual plots were tested for homoscedasticity and the presence of a nonlinear relationship. Finally, to assess for an independent association between SPV and data entry method (manually entered vs. automated recorded), a multivariable model including all univariate analysis variables that met collinearity criteria, plus a variable for the method of SPV data entry, was constructed for all patients.
Given the lack of previous studies describing multivariable-independent associations between shifts in SPV and PPV value distributions and clinical modifying factors, our sample size was selected based on the availability of reliable clinical data. Based upon study inclusion/exclusion criteria and procedural volume at our tertiary care center, we projected 1,000 to 2,000 patients included in analysis. To assess for robustness of study results, three sensitivity analyses were performed: (1) a study population excluding cases with epidurals dosed intraoperatively to supplement general anesthesia (spinal and regional supplemental anesthetics previously excluded), (2) a study population restricting SPV/PPV measurements to 15 min after anesthetic induction end and before surgical incision, and (3) a separate population of cases with an anesthetic induction end time documented between 9:00 AM and 1:00 PM.
We studied 1,791 cases, 1,323 with manually entered SPV values and 468 with automated SPV and PPV values (fig. 1). Patients’ perioperative characteristics are described in table 1. Compared to patients with manually entered SPV values, patients with automated recorded values were older, taller, less likely to be positioned supine, ventilated with lower tidal volumes and higher PEEP levels, and less likely to have general surgery. Before artifact removal, arterial line data quality assessment demonstrated that on a per-case basis, a mean of 95% and minimum of 90% of values were nonartifactual.
The 2.5 to 97.5 percentile reference interval as recommended by the Clinical and Laboratory Standards Institute37,38 for manually entered SPV was 3.0 to 11.0 mmHg, automated recorded SPV 2.2 to 10.4 mmHg, and PPV 2.0 to 16.0% (table 2). Interquartile ranges (IQRs) were 5.0 to 7.0, 3.9 to 6.0, and 5.0 to 9.0%, respectively (table 2). Manually entered values were higher than values in the automated recorded cohort (median, 6.0 [IQR, 5.0 to 7.0] mmHg vs. 4.7 [3.9 to 6.0] mmHg; P < 0.001; table 1). To adjust for the differences in patient characteristics between the manually and automated SPV, we used multivariable linear regression and found that manually entered values were slightly higher than automated values (0.76 mmHg ± 0.13 SE; P < 0.001).
Median SPV and PPV values for all cases included in the study formed distributions as shown in figures 2 and 3. The three sensitivity analyses (exclusion of cases with epidurals used to supplement general anesthesia, exclusion of SPV/PPV measurements after surgical incision, and a separate population of cases with anesthetic induction end time documented between 9:00 AM and 1:00 PM) yielded distributions with median/IQR values within 1.0 mmHg for SPV values and within 1.0% for PPV values (Supplemental Digital Content 1, https://links.lww.com/ALN/B341, Supplemental Digital Content 2, https://links.lww.com/ALN/B342, and Supplemental Digital Content 3, https://links.lww.com/ALN/B343).
By univariate analysis, factors associated with differences in SPV values included tidal volume, positioning (manually entered only), median PIP, preoperative β blocker, and preoperative diuretic (table 3). Statistically significant differences in PPV values were associated with age and preoperative angiotensin-converting enzyme inhibitor/angiotensin receptor blocker, β blocker, and diuretic (table 4).
No variables included in the univariate analysis were shown to have high collinearity. After adjustment for other factors using multivariable linear regression, nonsupine positioning and lack of preoperative β blockade were independently associated with increased SPV and PPV values, whereas elevated median ventilator tidal volume and elevated median PIP were each independently associated with increased SPV values only (Supplemental Digital Content 4, https://links.lww.com/ALN/B344, Supplemental Digital Content 5, https://links.lww.com/ALN/B345, and Supplemental Digital Content 6, https://links.lww.com/ALN/B346). MAP more than 20% below intraoperative baseline was independently associated only with PPV. Finally, a lack of preoperative diuretic was independently associated with increased SPV among manually entered values only. Although statistically significant, all independent associations identified were small (less than or equal to 1 mmHg or less than or equal to 2%). All other clinical factors demonstrated no significant independent associations.
Among 1,791 general anesthetics at our tertiary care academic institution, per-case median SPV and PPV values formed nonparametric distributions with medians of 4.7 mmHg (automated), 6.0 mmHg (manually entered), and 7.0%, respectively. Nonsupine positioning and lack of preoperative β blocker were independently associated with increased SPV and PPV, whereas elevated median ventilator tidal volume and elevated median PIP were each independently associated with increased SPV only. Based upon 25th to 75th percentiles, we conclude a “normal” range of manually entered SPV 5.0 to 7.0 mmHg, automated recorded SPV 3.9 to 6.0 mmHg, and PPV 5.0 to 9.0% among adult patients undergoing general anesthetics for elective noncardiac surgery. The SPV/PPV histograms derived from this study establish essential real-world reference distributions lacking in the current perioperative literature. These measurements allow for validation of previous studies and provide a prerequisite reference distribution for use in future studies assessing SPV/PPV-driven goal-directed therapies, as called for in the previous literature.21,39–41
In contrast to analyses of patients in the postcardiac surgery or intensive care unit (ICU) setting, our mean SPV and PPV values obtained were generally lower. In a meta-analysis of PPV among cardiac surgical and ICU patients, mean baseline PPV values of 7.1% for fluid challenge nonresponders and 16.6% for responders were observed.13 In studies analyzing SPV in cardiac surgical/ICU patients, mean baseline values were 8 to 14 mmHg.42–44 Compared to these patients, our elective general surgical population may have possessed factors (e.g., euvolemic state, absence of cardiopulmonary comorbidities, and absence of vasoactive infusions) predisposing to lower SPV/PPV values. Further studies are needed to investigate these differences; however, this difference illustrates the context sensitivity of SPV/PPV values and the value of our study investigating a broader, generally healthier elective surgical population.
Our study reports no independent association between obesity (BMI more than or equal to 30 kg/m2) and SPV/PPV values obtained intraoperatively; this finding sheds new light on a currently controversial concept. A theory that challenges measurement validity in the setting of obesity relates to increased intraabdominal pressure in the supine position for obese patients, for which a decrease in venous return results in an increase in arterial waveform pressure variation (rightward shift of SPV/PPV distributions). This relationship between intraabdominal pressure and arterial waveform pressure variation has been described in previous studies.20,45 Despite this relationship in controlled settings, our results question the impact of obesity on the SPV/PPV normal range. However, given that few patients had a BMI more than 40, we cannot draw conclusions about what impact extreme levels of obesity may have on SPV and PPV.
The statistically significant association between elevated tidal volumes and small increases in SPV values agrees with the previous literature.46–48 We were unable to demonstrate an association between tidal volumes more than 8 ml/kg ideal body weight and increased PPV values, in contrast to previous studies.47,48 While elevated PEEP failed to demonstrate any association, elevated median PIP demonstrated a stronger independent association compared to tidal volumes, with small increases in SPV values but not PPV values. Recent studies note similar findings, placing a stronger emphasis on the relationship between driving pressure and dynamic waveform indices compared to tidal volumes.16,49 In an era of lung-protective ventilation favoring lower tidal volumes and plateau pressures, our findings suggest a similar trend in SPV and PPV values that anesthesiologists may need to consider; indeed, the utility of an arterial pressure variation indexed to tidal volume has been demonstrated.48 Despite this association, a study also demonstrates that the ability of dynamic waveform indices to predict fluid responsiveness remains intact even at lung-protective ventilation levels and continues to greatly exceed predictive capabilities of static measures.18
Our finding of a positive association between nonsupine positioning and SPV values offers insight into the currently contrasting literature, which have examined alterations in dynamic arterial waveform indices associated with variations in surgical positioning, including prone,19,24,32 Trendelenburg,23,50 and reverse Trendelenburg.50 SVV and PPV demonstrated increases when prone,19 whereas SPV demonstrated no association.32 With regard to table positioning (e.g., Trendelenburg) as opposed to patient position, a lack of standardized documentation in our intraoperative record precluded this analysis and thus was a limitation to our study. Although the anesthesiologist must be aware of the possible effects of surgical positioning on SPV and PPV, the current literature is promising in that such indices remain effective measures of predicting fluid responsiveness in these settings.19,23,24,32
Our study demonstrated an independent association between hypotension and increased PPV values but not SPV values. These results elucidate previous conflicting studies, which have either shown a decrease in PPV associated with a MAP augmented via norepinephrine infusion51 or no association between MAP and PPV—and conversely, an association between MAP and SPV.28 Compared to previous studies, our analysis was performed retrospectively using clinical data, rather than in a controlled experimental environment, and thus, unmeasured confounders associated with hypotension may have influenced our results. Nonetheless, the observed associations represent valuable information for the anesthesiologist assessing SPV and PPV values in a clinical setting. As with other variables studied, it should also be noted that the independent effect size of hypotension on dynamic waveform indices was small (less than 1% for PPV).
Finally, our study noted lower SPV and PPV values for patients receiving preoperative β blockers and lower manually entered SPV values for patients receiving diuretics. Although the response of SPV/PPV measurements from medications altering vascular tone (e.g., phenylephrine and adenosine) has been investigated in the acute setting,28,52 no study to the authors’ knowledge has assessed the independent association between chronic cardiovascular medical therapies and SPV/PPV measurements. It is possible that such medical therapies are markers for underlying cardiovascular disease responsible for these associations; however, it should be noted that patients with heart failure diagnoses (among other comorbidities), as well as patients requiring vasopressor/inotrope infusions, were excluded from the study. Further prospective studies are needed to investigate this relationship. Nonetheless, chronic cardiovascular medical therapies serve as an important covariate in our multivariable analysis, enabling adjustment for other variables studied.
We also found a small difference (0.3 to 1.5 mmHg; table 2) in the distribution of SPV values between manually and automated entry (fig. 2). As the measurements were not made simultaneously in the same patients and persisted after multivariable analysis (0.76 ± 0.13 mmHg), these small differences may have reflected other unmeasured patient characteristics or differences between human calculations and computer algorithms. Further study is needed to clarify this issue.
Comparing SPV to PPV, PPV values in general were less affected by confounding factors. This finding in favor of PPV as a more robust measure aligns with the current literature, which supports PPV as a better predictor of fluid responsiveness.13 A limitation of our study was a lack of SVV measurement analysis. However, while studies have shown that SVV has a similar accuracy at predicting fluid responsiveness as SPV, it is less accurate than PPV.21–24
Our study possessed several additional limitations. Our study was performed retrospectively and was subject to limitations inherent to study design. Data were available as charted within the perioperative database; additional details beyond the scope of clinical care were unavailable. Similarly, given our study design, patient responses to therapeutic interventions were not assessed; this warrants future prospective study investigation. In addition, although standard clinical practices were employed, no specific study protocols—including standard methods for manually entered SPV measurement—were utilized. Despite anesthesiologist training on data entry, manual auditing of sample cases, and the use of data validation algorithms, no standard of data quality assurance could be applied at the point of care. All PPV values were obtained via the algorithm employed by the CARESCAPE B850 monitor; automated PPV calculation variations may exist across different manufacturers. Although our study identified a reference distribution of SPV/PPV values lacking in the current literature, potential confounders—including those used as exclusion criteria, those studied in univariate analyses, and those remaining unidentified—limit the generalizability of our results. Conversely, the reference intervals identified by our study were generated by a broad surgical population, and despite consideration of modifying factors, application to individual surgical cases with multiple modifying factors must be made with caution. Additionally, although measures to target a euvolemic patient population included restricting to a morning case time, elective procedures, and patients with minimal cardiopulmonary comorbidities, no rigorous measures were implemented to confirm a particular patient volume status. We performed our study at a single academic tertiary care center, and institutional variations, scope of practice, and geographic region may bias our results. Our study population drawing from convenience samples for SPV and PPV measurements represents the largest population currently studied for this purpose and was able to detect clinically meaningful associations between SPV/PPV values and clinical factors. However, a larger sample size may be required to uncover smaller, yet statistically significant associations for other factors.
Despite these limitations, our study successfully addresses a shortcoming in the current literature through establishing reference distributions of SPV and PPV values for patients undergoing elective surgical procedures. By providing a means of validation for previous studies, as well as a platform for assessing the impact of subsequent interventions, as well as goal-directed therapies, our study furthers the knowledge of two widely used clinical measures associated with fluid responsiveness and offers a standard reference range to the anesthesiologist for volume assessment and intraoperative fluid management therapies.
The authors acknowledge Hyeon Joo, M.S., University of Michigan Health System, Ann Arbor, Michigan, for his contributions in data acquisition and electronic search query programming for this project.
Supported in part by the National Center for Research Resources (grant no. UL1RR024986; Bethesda, Maryland), which is now at the National Center for Advancing Translational Sciences (grant no. 4UL1TR000433; Bethesda, Maryland). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health (Bethesda, Maryland).
The authors declare no competing interests.