Skip to main content

Identification of novel sublingual parameters to analyze and diagnose microvascular dysfunction in sepsis: the NOSTRADAMUS study

Abstract

Background

The availability of handheld, noninvasive sublingual video-microscopes allows for visualization of the microcirculation in critically ill patients. Recent studies demonstrate that reduced numbers of blood-perfused microvessels and increased penetration of erythrocytes into the endothelial glycocalyx are essential components of microvascular dysfunction. The aim of this study was to identify novel microvascular variables to determine the level of microvascular dysfunction in sepsis and its relationship with clinical variables.

Methods

This observational, prospective, cross-sectional study included 51 participants, of which 34 critically ill sepsis patients were recruited from intensive care units of a university hospital. Seventeen healthy volunteers served as controls. All participants underwent sublingual videomicroscopy by sidestream darkfield imaging. A new developed version of the Glycocheck™ software was used to quantify vascular density, perfused boundary region (PBR-an inverse variable of endothelial glycocalyx dimensions), red blood cell (RBC) velocity, RBC content, and blood flow in sublingual microvessels with diameters between 4 and 25 µm.

Results

A detailed analysis of adjacent diameter classes (1 µm each) of vessels between 4 and 25 µm revealed a severe reduction of vascular density in very small capillaries (5–7 µm), which correlated with markers of sepsis severity. Analysis of RBC velocity (VRBC) revealed a strong dependency between capillary and feed vessel VRBC in sepsis patients (R2 = 0.63, p < 0.0001) but not in healthy controls (R2 = 0.04, p = 0.43), indicating impaired capillary (de-)recruitment in sepsis. This finding enabled the calculation of capillary recruitment and dynamic capillary blood volume (CBVdynamic). Moreover, adjustment of PBR to feed vessel VRBC further improved discrimination between sepsis patients and controls by about 50%. By combining these dynamic microvascular and glycocalyx variables, we developed the microvascular health score (MVHSdynamic™), which decreased from 7.4 [4.6–8.7] in controls to 1.8 [1.4–2.7] in sepsis patients (p < 0.0001) and correlated with sepsis severity.

Conclusion

We introduce new important diameter-specific quantification and differentiated analysis of RBC kinetics, a key to understand microvascular dysfunction in sepsis. MVHSdynamic, which has a broad bandwidth to detect microvascular (dys-) function, might serve as a valuable tool to detect microvascular impairment in critically ill patients.

Background

Alterations of microvascular perfusion caused by endothelial cell dysfunction, glycocalyx degradation, increased leukocyte adhesion, microthrombus formation, and regional redistribution of blood flow contribute to organ failure in critically ill patients [1,2,3]. The use of handheld, noninvasive sublingual video microscopes, such as incident dark field (IDF) or sidestream dark field (SDF) microscopy, allows researchers to visualize and analyze microvascular perfusion alterations in vivo [4, 5]. Over the past decade, several measures of microvascular perfusion variables have been proposed to improve risk stratification, prognostication and eventually to individualize the therapy, especially in critically ill patients with sepsis and septic shock [6,7,8,9,10,11,12,13]. Despite tremendous efforts to improve video quality and standardize reporting, the quantitative analysis of functional measures remains challenging and is often performed manually [6]. Several visual scoring systems and specialized (semi-)automated analysis tools have been developed to shorten time-to-result and improve accuracy [14,15,16,17]. Common to all these approaches is the (automated) selection and pooled analysis of microvessels with a diameter < 20 (rarely < 10) µm. An improved spatiotemporal resolution within this < 20 µm range is desirable to further dissect and understand functional microvascular alterations in sepsis and critical illness.

We have previously used the Glycocheck™ System to analyze the perfused boundary region (PBR), an inverse variable of endothelial glycocalyx (eGC) dimensions, in sublingual microvessels from sepsis patients [18,19,20]. The software automatically detects, records, and analyzes the PBR according to diameter classes (each 1 µm) in microvessels with diameters between 4 and 25 μm. The methodology showed very good inter- and intra-observer reproducibility under real-life conditions [19] and excellent accuracy compared to other in vitro and ex vivo estimations of eGC thickness [18, 20, 21]. The aim of this study was to identify novel dynamic microvascular variables generated by a new diameter-class-wise approach to determine the level of microvascular dysfunction in sepsis and its relationship with established clinical variables and the number of dysfunctional organs.

Materials and methods

Study population

This prospective, observational, cross-sectional study took place in the medical and operative ICUs of the University Hospital Münster. The study was performed in accordance with the Declaration of Helsinki and was approved by the competent ethics committee (amendment to 2016-073-f-S). Some of the participants were already included in previous studies on eGC damage in sepsis [18, 20].

After written informed consent was obtained from the patients or their legal representatives, 34 adult ICU patients with sepsis were enrolled non-consecutively after initial resuscitation. Sepsis was defined by the sepsis-3 criteria published by the ESICM-SCCM Sepsis Redefinitions Task Force [22]. Exclusion criteria were age < 18 years, pregnancy, or oral mucosal inflammation or injury, which could locally influence the sublingual microvasculature. Seventeen apparently healthy volunteers served as controls.

Demographic variables, routine chemistry tests, and physiological variables, including the Sequential Organ Failure Assessment (SOFA) score [22] and a contemporary version of the Charlson Comorbidity Index (CCI) [23], were obtained for each subject at the time of sublingual videomicroscopy (Table 1). Each videomicroscopy set consisted of two complete measurements (see below) which were averaged to account for spatial heterogeneity of the sublingual microvasculature.

Table 1 Baseline characteristics

Video acquisition and analysis of the sublingual microvasculature

Figure 1 provides an overview of the process of video acquisition, data analysis and post-processing. Bedside intravital microscopy was performed with a sidestream dark field (SDF) camera (CapiScope HVCS, KK Technology, Honiton, UK) to visualize the sublingual microvasculature according to a standardized procedure as described in detail previously [20]. A physician experienced in the SDF technique, trained to recognize and avoid pressure artifacts performed the image acquisition based on the current round table recommendations [6]. The SDF camera uses green light-emitting stroboscopic diodes (540 nm) to detect the hemoglobin of passing red blood cells (RBCs). Using a 5 × objective with a 0.2 numerical aperture, images were captured, providing a 325-fold magnification in 720 × 576 pixels at 23 frames per second as described in detail previously [19, 24, 25].

Fig. 1
figure 1

Flowchart showing the process of video acquisition and data analysis. D = Vessel Diameter, in µm

The analysis of all variables was carried out exclusively with the GlycoCheck™ Software (Microvascular Health Solutions Inc., Salt Lake City, UT, USA) software. Briefly, the GlycoCheck™ Software allows video acquisition after predefined image quality criteria (motion, intensity, and focus) are fulfilled. Specifically, light intensity range is automatically set to prevent under- or over-exposure of the video. The acceptable motion range is set so that every vessel can be tracked in all frames of each video. Each complete measurement consists of at least ten 2-s videos (40 frames/video), containing a total of about 3000 vascular segments of 10 μm length each. All videos are deliberately obtained from different positions to counterbalance spatial heterogeneity of the sublingual microcirculation. The software automatically subjects the obtained vascular segments to a strict quality check, as previously described in detail [25]. Briefly, vascular segments are considered valid only if: (1) red blood cell content is ≥ 50%, (2) RBC width is ≥ 2 µm, (3) segments are neither curved, (4) nor too close to each other. Invalid vascular segments are marked yellow and are automatically discarded, while all valid vascular segments (green lines) are further analyzed. Therefore, all further parameters are calculated based on these valid segments. Screenshots of the sublingual mucosa of five randomly chosen individuals are shown in Additional file 1: Fig. S1 (with and without automatically acquired vessel detection and quality check). Moreover, two additional movie files illustrate the sublingual mucosa of a healthy control and a sepsis individual with and without automatic vessel detection and subsequent quality check in more detail (Additional file 2: Video 1, Additional file 3: Video 2). Finally, the software obtains up to 840 radial intensity profiles for each valid vascular segment based on the RBC column width (RBCW), and automatically groups vessels from 4 to 25 µm diameter in 22 separate diameter classes (1 µm each). This concept has been successfully used and validated in the past [25]. Data from two complete measurements (hereafter referred to as “measurement set”) were extracted, analyzed, and averaged offline to avoid sampling error and to counterbalance spatiotemporal heterogeneity of the sublingual microcirculation [19, 20].

RBC velocity

Valid red blood cell (RBC) velocities (VRBC) are determined in individual vessel segments in an automatic fashion by cross correlation of longitudinal RBC intensity profiles between frames of recorded videos. A correlation coefficient of ≥ 0.80 was required to ensure accurate estimates of longitudinal RBC displacement. RBC velocity is determined by dividing RBC displacement by the time between video frames and expressed as µm/s. The automatic estimation of the VRBC correlated excellently with a manual analysis (R2 = 0.94, p < 0.0001; Additional file 4: Fig. S2).

Perfused vascular and capillary density

An absolute measure for valid perfused vascular density (mm/mm2) can be determined from the number of vascular segments containing RBCs multiplied by capillary segment length (each 10 µm). All detected RBC-containing vessel segments (RBC content ≥ 50%) with VRBC ≥ 0 µm/s were automatically counted in the video recordings of each subject. Vascular density was normalized to tissue surface area. As non-perfused vessels (i.e., without RBCs present), as well as vessels not meeting the quality criteria cannot be detected by our methodology, therefore vascular density in this manuscript refers to valid perfused vascular density (hereafter vascular density). Capillary density was defined as vascular density of vessels with a diameter equal or smaller than the diameter of a single red blood cell (RBC diameter ~ 7 to 8 µm [26]; capillary density D ≤ 7 µm).

Absolute and static capillary blood volume

An absolute measure for valid capillary blood volume (CBVabsolute) can be determined from the number of capillary segments multiplied by capillary segment length (i.e., capillary density (mm/mm2)) and segment-specific capillary cross-sectional area (\(\pi\) * radius2). In addition to counting the number of RBC-containing capillary segments, a functional estimate of CBV relative to larger vessel blood volume can be determined by measuring average VRBC in capillaries and larger blood vessels (hereafter feed vessels). The corresponding VRBC ratio in feed vessels (diameter 10 to 25 µm, hereafter D ≥ 10 µm) over VRBC in capillaries (diameter 4 to 7 µm, hereafter D ≤ 7 µm) denotes the CBV ratio (\(= V_{RBC} \left( {D \ge 10\,\upmu {\text{m}}} \right)/V_{RBC} \left( {D \le 7\,\upmu {\text{m}}} \right)\)). In short, an increase in capillary blood volume relative to feed vessel blood volume will reduce capillary VRBC and this will increase the CBV ratio. Multiplying the CBVabsolute with the CBV ratio gives the static capillary blood volume (CBVstatic), which is defined as: \(CBV_{static} = CBV_{absolute} *V_{RBC} \left( {D \ge 10\,\upmu {\text{m}}} \right)/V_{RBC} \left( {D \le 7\,\upmu {\text{m}}} \right)\).

Capillary recruitment and dynamic capillary blood volume

To take the ability to recruit additional capillaries into account, an estimate can be made of the capillary recruitment (CR) by measuring the slope of the relationship between VRBC (D ≤ 7 µm) and VRBC (D ≥ 10 µm). When the number of blood perfused capillaries increases upon and increase in VRBC (D ≥ 10 µm), the accompanying increase in VRBC (D ≤ 7 µm) will be less than proportional (i.e., the regression slope will be < 1) and capillary recruitment can be defined as 1 − slope(VRBC (D ≤ 7 µm), VRBC (D ≥ 10 µm)). Two extreme examples of this concept are shown in Additional file 5: Fig. S3: In case capillary blood volume doubles when large vessel RBC velocity increases twofold, the slope(VRBC (D ≤ 7 µm), VRBC (D ≥ 10 µm)) will be 0 and CR = 1 − slope 0 = 1 = 100%. In the absence of changes in capillary blood volume when VRBC (D ≥ 10 µm) increases twofold, capillary RBC velocities are expected to also change proportionally by twofold, and the slope of VRBC (D ≤ 7 µm) vs. VRBC (D ≥ 10 µm) will be 1 and CR = 1 − slope 1 = 0 = 0%. Multiplying the CBVstatic * (1 + CR) gives the dynamic capillary blood volume (CBVdynamic).

Static and dynamic perfused boundary region (PBR)

The software calculates the dynamic lateral movement of RBCs into the permeable part of the eGC layer, expressed as the PBR (in μm). An impaired eGC permits a greater number of RBCs to penetrate deep into the endothelium, which is translated as an increase in the PBR value. The radial distribution of RBCs in each valid segment defines the median RBC width (RBCW), as well as the outer edge of the RBC-perfused vessel diameter (Dperf). The PBRstatic is defined as the distance between the RBCW and Dperf and is calculated using the following formula: \(\left( {Dperf {-} RBCW} \right)/2\) [25].

In line with the above-described measures of velocity-dependent increases in capillary blood volume, it is possible that also penetration of RBCs into the luminal glycocalyx surface (as reflected by PBR) is velocity-dependent. To minimize possible flow-dependent variability in PBR estimation, the slope(PBRstatic, VRBC (D ≥ 10 µm)) can be used to PBRdynamic under equal conditions, i.e., the absence of RBC velocity (VRBC (D ≥ 10 µm) set to 0 µm/s) (Additional file 6: Fig. S4 for details).

Statistical analysis

Data are presented as absolute numbers, percentages, or medians with corresponding 25th and 75th percentiles (interquartile range; IQR), as appropriate. The nonparametric Mann–Whitney U test and the chi-square test were used to compare variables between patients and controls. To correct for multiple testing in comparisons of microcirculation variables per-diameter class we used the false discovery rate (FDR) approach of Benjamini, Krieger and Yekutieli, setting a qvalue < 0.05 as significant. Spearman rank correlation coefficient was used to assess correlations between clinical and microvascular variables. Associations between microvascular measures were evaluated using simple and adjusted linear regression models. All the tests used were two-sided, and statistical significance was set at p < 0.05. Our study was powered to detect a moderate correlation (Spearman correlation coefficient = 0.5) between microvascular health score and SOFA score in the septic cohort with 85% power given a two-sided alpha of 0.05 [27]. SPSS version 26 (IBM Corporation, Armonk, NY, USA) and GraphPad Prism version 8.4.3 (GraphPad Prism Software Inc., San Diego, CA, USA) were used for statistical analyses and preparation of figures.

Results

The clinical and demographic characteristics of the 51 study participants are shown in Table 1 and Additional file 7: Table A1. 50% of the patients had a respiratory focus of infection, while the remainder demonstrated other etiologies (Additional file 7: Table A1). From a total of 34 sepsis patients, 22 (65%) required organ replacement therapy (mechanical ventilation or/and renal replacement therapy) and 26 (76%) were vasopressor-dependent at study inclusion. Our sepsis cohort had a median (IQR) SOFA score of 10 (8–13), indicating moderate disease severity. The median time of sepsis begin was 2 days (1–4) and 6 (18%) patients had septic shock (Sepsis-3 definitions: lactate > 2 mmol/l and need for vasopressors to maintain a mean arterial pressure ≥ 65 mm Hg) at the time of study inclusion.

Analysis of sublingual microcirculation according to vessel diameter class

First, we compared vascular density, VRBC and PBRstatic between healthy controls and sepsis patients in a diameter-class-wise fashion. This approach revealed a statistically significant decrease in vascular density only in the diameter classes 5, 6 and 7 µm in sepsis patients (− 63%, − 42%, and − 28% compared to controls), whereas the remaining diameter classes from 8 to 25 µm were not different between the groups (Fig. 2a and Additional file 8: Fig. S5). However, vascular density was not different between the groups in a pooled analysis (4–25 µm range) (12.3 vs. 12.8 mm/mm2, p = 0.34). To check our results for plausibility, we correlated, in an explorative manner, vascular density with clinical variables, such as interleukin-6 (IL-6), procalcitonin (PCT), lactate and SOFA score. We found a robust inverse correlation in the 4 to 7 µm diameter range, while these associations were absent or even partially reversed in larger vessels (Fig. 2b and Additional file 7: Table A2).

Fig. 2
figure 2

Analysis of sublingual microcirculation according to vessel diameter class. a, c, e Median and IQR values of vascular density, PBRstatic values and RBC velocity of healthy controls and sepsis patients according to diameter class from 4 to 25 µm. b, d, f Bar charts showing the correlation coefficient (Spearman) between microvascular and clinical variables. IQR inter quartile range, IL6 interleukin 6, PBR perfused boundary region, PCT procalcitonin, RBC red blood cell, SOFA score sequential organ failure assessment score. Q value (adjusted p value): *q < 0.05, **q < 0.01, ***q < 0.001

VRBC showed a huge overlap between patients and controls. Median VRBC trended to be slightly lower in sepsis patients, especially in smaller vessels and was inversely correlated with variables of disease severity (Fig. 2c, d and Additional file 7: Table A2). A significant sepsis-induced increase in PBRstatic was particularly apparent in feed vessels, probably due to loss of the affected smaller vessels (lower capillary density) (Fig. 2e and Additional file 9: Fig. S6). PBRstatic values of individual diameter classes showed weak-to-moderate positive correlations with variables of sepsis severity (Fig. 2f and Additional file 7: Table A2).

Taken together, the diameter class-wise data analysis revealed individual sepsis-induced perturbations in vascular density, VRBC and PBRstatic.

Derivation of capillary recruitment and dynamic capillary blood volume.

Having observed a significant decrease in vascular density in very small capillaries (diameter (D) ≤ 7 µm) in sepsis patients, we determined and focused on the capillary blood volume (Fig. 3). CBVabsolute decreased from 16.5 [10.3–19.4] 103μm3 in healthy controls to 7.9 [5.9–14.5] 103μm3 in sepsis patients (p = 0.006). To further improve discrimination between the groups, we added some functional variables to the CBV calculation.

Fig. 3
figure 3

Derivation of capillary recruitment and dynamic capillary blood volume. a Box plots showing capillary blood volume (CBV) ratio in healthy controls (green) and sepsis patients (red). CBV ratio denotes the RBC velocity (VRBC) in feed vessels (D ≥ 10 µm) over VRBC in capillaries (D ≤ 7 µm). b Scatter dot plots and simple linear regression (slope) with 95% confidence intervals of VRBC in capillaries (D ≤ 7 µm) plotted against VRBC in feed vessels (D ≥ 10 µm). Different states at the ends of the slope lines (indicated by green/red bold letters A-D) are further explained in Fig. 5. c Bar charts showing the capillary recruitment (CR), defined as 1 − slope (VRBC (D ≤ 7 µm) vs. VRBC (D ≥ 10 µm)) per group. d Box plots showing the development of different measures and estimates of CBV. Left: CBVabsolute is determined from the number of capillary segments multiplied by capillary segment length (i.e., capillary density (mm/mm2)) and segment-specific capillary cross-sectional area (π * radius2). Briefly, as each vascular segment can be considered a cylinder, the segment-specific capillary cross-sectional area can be calculated with the mathematical type π * radius2 (circle’s area). The radius is estimated every 10 µm (segment’s length) and recorded accordingly. Middle: CBVstatic is calculated as CBVabsolute * VRBC (D ≥ 10 µm)/VRBC (D ≤ 7 µm). Right: CBVdynamic is defined as CBVstatic * (1 + CR). RBC red blood cell, CBV capillary blood volume, CR capillary recruitment, D diameter, V velocity

First, we calculated CBVstatic, by multiplying CBVabsolute with a ratio derived from VRBC in feed vessels relative to capillaries (VRBC (D ≥ 10 µm)/VRBC (D ≤ 7 µm)). However, the suspected decrease in the CBV ratio (due to decrease in capillary density) in sepsis patients did not occur (Fig. 3a). The reason for this only became apparent after plotting capillary VRBC (D ≤ 7 µm) as a function of large vessel VRBC (D ≥ 10 µm) (Fig. 3b). This type of analysis revealed a strong dependency between VRBC (D ≤ 7 µm) and VRBC (D ≥ 10 µm) in sepsis patients (R2 = 0.63, p < 0.0001), indicating impaired capillary (de-)recruitment in this group. In contrast, capillary VRBC is relatively constant in healthy controls as reported before [28], indicating functioning (de-)recruitment of CBV associated with changes of feed vessel blood flow in healthy subjects. CBVstatic decreased from 18.5 [11.1–22.3] 103μm3 in healthy controls to 10.2 [7.1–17.1] 103μm3 in sepsis patients (p = 0.005).

Next, we calculated capillary recruitment (CR) as 1- the slope(VRBC (D ≤ 7 µm), VRBC (D ≥ 10 µm)). CR per group decreased from 78% in healthy controls to 8% in sepsis patients (Fig. 3c). Finally, we calculated CBVdynamic by multiplying CBVstatic with (1 + CR). CBVdynamic showed a significant higher and much wider normal range in healthy controls compared to CBVabsolute and CBVstatic, respectively. CBVdynamic decreased from 32.8 [19.7–39.5] 103μm3 in healthy controls to 11.1 [7.7–18.5] 103μm3 in sepsis patients (p < 0.0001). Overall, we were able to show that by adding dynamic variables to CBV calculation, the discrimination between the groups could be significantly improved (Fig. 3d).

Derivation of an RBC velocity-adjusted perfused boundary region

We have previously shown that PBR values (averaged across all diameter classes from 5 to 25 μm) increase in sepsis patients [18,19,20], indicating sepsis-induced damage of the eGC. In this study, PBRstatic (D 4 to 25 µm) increased from 2.24 [2.13–2.35] µm in healthy controls to 2.48 [2.33–2.62] µm in sepsis patients (p < 0.0001) (Fig. 4a). In sepsis patients, PBRstatic showed a tendency to increase at low levels of feed vessel velocity, whereas in healthy controls, PBRstatic showed a tendency to stay constant or even slightly decrease at low feed vessel velocity (Fig. 4b). To eliminate the influence of different RBC velocities on PBRstatic, we estimated PBRdynamic at a feed vessels VRBC of 0 µm/s based on the different slopes(PBRstatic, VRBC (D ≥ 10 µm)) in healthy subjects and sepsis patients (Fig. 4b and Additional file 6: Fig. S4). PBRdynamic increased from 1.95 [1.84–2.08] µm in healthy controls to 2.58 [2.43–2.76] µm in sepsis patients (p < 0.0001) (Fig. 4c). In summary, adjustment of the PBR to VRBC improved discrimination between the groups by about 50%.

Fig. 4
figure 4

Derivation of an RBC velocity-adjusted perfused boundary region. a Box plots showing PBRstatic (D 4 to 25 µm) in healthy controls (green) and sepsis patients (red). b Scatter dot plots and simple linear regression (slope) with 95% confidence intervals of PBRstatic (D 4 to 25 µm) plotted against VRBC in feed vessels (D ≥ 10 µm). Different states at the ends of the slope lines (indicated by green/red bold letters A-D) are further explained in Fig. 6. c Box plots of PBRdynamic, a velocity-adjusted estimate of the PBR where the VRBC (D ≥ 10 µm) was set to 0 µm/s. RBC red blood cell, CBV capillary blood volume, D diameter, V velocity

Development of the Microvascular Health Score

We have previously shown, that glycocalyx damage and microcirculatory impairment do neither coincide, nor do they occur in proportion in every sepsis patient [20]. This un-coupling of changes in PBR and CBV can also be reproduced in the group of sepsis patients (Fig. 5a–c). To account for this finding, we combined the different variables into one Microvascular Health Score (MVHS™) (Additional file 10: Fig. S7). Because CBV decreases and PBR increases in sepsis, we formed the quotient of CBV/PBR to calculate the MVHS in a static and a dynamic version (MVHSstatic = CBVstatic (D 4 to 6 μm)/PBRstatic (D 4 to 25 μm) and MVHSdynamic = CBVdynamic (D 4 to 6 μm)/PBRdynamic (D 4 to 25 μm)).

Fig. 5
figure 5

Association of the CBVdynamic, PBRdynamic and dynamic Microvascular Health Score (MVHSdynamic) with disease severity. Association of a CBVdynamic and b PBRdynamic with sequential organ failure assessment score (SOFA) score after dichotomizing (median) the sepsis group. c Scatter dot plots and simple linear regression (slope) with 95% confidence intervals of CBVdynamic plotted against PBRdynamic (D 4 to 25 µm) in the septic population. The red dotted lines represent the median values of CBV and PBR, respectively. d Association of MVHSdynamic with sequential organ failure assessment score (SOFA) score after dichotomizing (median) the group

The median MVHSstatic was significantly lower in sepsis patients 1.70 [1.32–2.53] points compared to healthy controls 4.03 [2.11–4.33] points (p < 0.001). The median MVHSdynamic stressed the difference between sepsis patients and healthy controls (1.78 [1.38–2.67] vs. 7.43 [4.65–8.73] points, p < 0.0001) even more than the MVHSstatic.

Association of the Microvascular Health Score (MVHS) with disease severity

Both, MVHSstatic and MVHSdynamic correlated moderate-to-strong with SOFA score, number of dysfunctional organs, lactate, CRP, IL-6 and PCT (all p < 0.001) in a pooled analysis that included healthy controls and sepsis patients (Table 2). These observations remained unchanged in additional sensitivity analyses (Additional file 7: Table A3).

Table 2 Correlations of MVHS with clinical variables

In the subgroup of sepsis patients, the correlations of MVHSstatic and MVHSdynamic with SOFA score and number of dysfunctional organs remained significant (all p ≤ 0.01) (Table 2 and Additional file 7: Table A4). MVHSdynamic was significantly different in sepsis subgroups stratified by either median SOFA score or median number of dysfunctional organs, respectively (Fig. 5d and Additional file 11: Fig. S8).

Discussion

In this study, we employed a differentiated, diameter class-wise analysis of RBC kinetics to identify novel variables of microvascular dysfunction in sepsis. This approach revealed several size- and group-specific characteristics of the measured variables within the diameter-range from 4 to 25 µm. Accordingly, we propose the MVHSdynamic, which has a broad bandwidth to detect microvascular (dys-)function in healthy subjects and critically ill patients. To our knowledge, the < 20 µm range has neither been divided into single µm-steps nor individually examined and compared with clinical information from healthy and septic subjects before.

Capillary recruitment

Observing that RBC dynamics differ between capillaries (D ≤ 7 µm) and feeding vessels (D ≥ 10 µm) led us to relate capillaries and feeding vessels to each other—which resulted in the calculation of capillary recruitment. The concept of capillary recruitment in skeletal muscle was first proposed by August Krogh in 1919 [29]. Krogh hypothesized that the opening of previously closed muscle capillaries, i.e., capillary recruitment, would allow capillary RBC velocities to remain low, and capillary oxygen extraction thereby efficient, despite large increases in blood supply. While the binary distinction between "closed" and "open" capillaries, was key to Krogh's model argument, modern observations suggest that capillary recruitment should be regarded rather in the context of continuous than binary changes in RBC distributions and velocities among capillaries [30, 31]. In this respect, our quantitative data show that healthy individuals (at least under resting conditions) have relatively constant median capillary RBC velocities, which seems to be independent of the associated VRBC in the feed vessels. Even though the underlying (auto-) regulation mechanisms remain incompletely understood, compliance with a narrow RBC velocity range seems plausible to guarantee an optimal supply–demand ratio in capillaries. Conceptually, constant capillary RBC velocities in the face of varying VRBC in larger feeding vessels is consistent with a) increasing numbers of RBC-perfused capillaries at tissue sites with higher metabolic demand and increased blood supply in feeding vessels, as well as b) reduced numbers of RBC-perfused capillaries at sites with lower metabolic demand and reduced feeding vessel blood supply. The failure to maintain constant RBC velocities in capillaries of patients with sepsis and the fact that capillary RBC velocities in our septic cohort change proportionally with RBC velocity changes in feedings vessels, reflects that the number of perfused capillaries in sepsis patients is fixed and insensitive to local variations in tissue metabolic demand. A simplified version of the pathophysiologic concept of capillary (de-)recruitment observed in this study is visualized in Fig. 6.

Fig. 6
figure 6

Cartoon showing the pathophysiologic concept of capillary (de-)recruitment in healthy conditions (a, b) and during sepsis (c, d). The respective ends of the spectrum of actual measured RBC velocities are indicated in the red circles. Conditions shown in A-D refer to corresponding green/red bold letters in Fig. 3B). In theory, higher RBC velocity states (low supply demand ratio, b, d) should go along with slightly thicker glycocalyx than lower RBC velocity states (high supply demand ration a, b). Furthermore, RBC flux is actually controlled by arterioles rather than precapillary sphincters. These details have been omitted for the sake of clarity

Estimates of capillary recruitment in our cohort are based on per group analysis of relations between RBC velocities in capillaries and feeding vessels. New recording and analysis strategies that allow for calculation of capillary recruitment on a per-patient basis are currently under investigation. This modification will allow for capturing enough low- and high flow situations per subject to generate intra-individual regression slopes. Thus, it will be possible to determine the capillary recruitment per individual subject, which might further improve discrimination of dynamic variables, such as CBVdynamic and MVHSdynamic.

Endothelial glycocalyx damage

The sepsis-associated PBR increase measured in this study is consistent with changes observed in previous studies [18,19,20, 32,33,34]. The current study shows that adjustment of PBR estimates to VRBC improved discrimination between the groups by about 50%. This gain in discriminatory power is important, as PBR estimates have so far shown a considerable overlap between healthy controls and sepsis patients. This adjustment accounts for RBC velocity dependency of the glycocalyx [35, 36]. Especially in sepsis patients, PBR tends to increase when RBC velocity decreases, suggesting that glycocalyx extends towards a more porous and permeable compartment in the absence of exposure to fluid/RBC shear forces. The fact that PBR increases are more extreme in the sepsis patients, is consistent with previous studies showing that glycocalyx is damaged and more permeable in sepsis, which likely contributes to its increased RBC velocity dependency. However, it is difficult to estimate if the higher PBR observed in sepsis patients is caused solely by real glycocalyx damage (i.e., enzymatic shedding [18, 37]) or if temporary and functional changes that enable nutrient and solute supply to the interstitial space contribute to that PBR increase. We have previously shown that PBR estimates correlate excellent with glycocalyx thickness measured by atomic-force microscopy in cultured endothelial cells exposed to concordant sera from healthy subjects and sepsis patients [18, 20]. Glycocalyx damage in vitro was completely abolished when enzymatic activity of heparanase, a heparan sulphate specific endo-beta-D-glucuronidase, was blocked by addition of heparin [18]. Furthermore, circulating levels of Syndecan-1, a core protein of the intact glycocalyx, increased tenfold in sepsis patients and correlated with PBR values [20], indicating that enzymatic shedding is a dominant mechanism of PBR increase in sepsis. Elegant work in the cremasteric microcirculation showed that, inflammation-induced glycocalyx shedding increases effective capillary diameter. This effect is, however, defeated by enhanced white blood cell (WBC)-endothelial interactions and subsequent venular WBC obstruction, resulting in a net reduction of capillary VRBC [38]. Glycocalyx damage may therefore partly explain the trend towards lower VRBC in capillaries.

Microvascular health score

In the current cross-sectional cohort, the proposed MVHSstatic and MVHSdynamic did correlate with SOFA score and number of dysfunctional organs in the subgroup of sepsis patients. Correlations between classical variables of the microcirculation and the SOFA score have been rarely reported. The SOFA correlated negatively with microvascular flow indices (MFI) as well as total and perfused vascular densities (TVD and PVD) in sepsis patients receiving activated protein C [33]. However, we couldn´t detect any significant correlations of MFI, TVD or PVD with the SOFA score in a previous sepsis study [20]. The finding that the MVHS performed better than its individual components is a strong argument for combining different measures of microvascular dysfunction such as vascular blood volume (CBV) and glycocalyx properties (PBR) into one score. In our view, in particular the combination of CBV and PBR makes sense, because CBV and PBR appear to behave independently of each other [20], and may thus indicate different aspects of microvascular dysfunction. To clarify whether the MVHS can indeed predict relevant outcome in sepsis patients, we initiated two prospective, observational, longitudinal studies to evaluate the MVHS in the emergency room (Early Detection of Glycocalyx Damage in Emergency Room Patients—the EDGE Study, Clinicaltrials.gov Identifier: NCT03126032) and in the ICU (Analysis of Sublingual Glycocalyx Damage at ICU Admission to Predict Risk of Death—the ASGARD Study, Clinicaltrials.gov Identifier: NCT03847493). In these studies, serial measurements can possibly reveal differential changes in individual components of the MVHS in response to the initial therapy.

Limitations

We acknowledge some limitations of our study. First, it is a single-center study with a limited sample size. Therefore, the findings cannot be directly generalized for all sepsis patients. However, the interdependencies between our novel variables and their associations with clinical variables remained significant after multiple adjustments and additional sensitivity analyses. Second, this study was neither designed nor powered to test the performance of novel variables for outcome prediction. However, our findings add interesting new aspects to the lively field of intravital microscopy research. Further longitudinal studies are needed to evaluate the MVHS for outcome prediction. Third, since all calculations are ultimately based on the flow properties of RBCs, we can only analyze microvessels in which a minimal number of RBCs are present and the predefined quality criteria are met. Vessels without RBCs or invalid vascular segments are therefore not detected by the software. That could have an impact not only on the density calculations, but also on the estimation of the PBR, as severely affected capillaries might not be accessible to RBCs anymore. However, past studies have shown, that the average PBR (4–25 µm) of healthy and septic subjects measured in vivo correlates excellently with the decrease in glycocalyx thickness of endothelial cells after incubation with the concordant serum samples, further supporting the accuracy of PBR values. Fourth, technical limitations, as, e.g., pressure artifacts or low spatial resolution might have affected the results. However, all measurements were performed by only one very experienced investigator specially trained to avoid pressure artifacts. Moreover, the camera has a sufficiently high resolution.

Conclusion and outlook

We used a highly differentiated analysis to quantify RBC dynamics in sublingual microvessels and found that taking the interdependence between capillary blood volume, capillary recruitment and PBR into account, is key to better understand microvascular dysfunction in sepsis. The MVHSdynamic reflects an impairment of very small capillaries and has a broad bandwidth to detect microvascular dysfunction in critically ill patients. Future clinical studies should evaluate the prognostic value of MVHS in sepsis.

Availability of data and materials

The datasets used and/or analyses during the current study are available from the corresponding author on reasonable request.

Abbreviations

ASGARD:

Analysis of Sublingual Glycocalyx Damage at ICU Admission to Predict Risk of Death

BMI:

Body mass index

CBV:

Capillary blood volume

CCI score:

Charlson Comorbidity Index

CR:

Capillary recruitment

CRP:

C-reactive protein

D:

Diameter

EDGE:

Early detection of glycocalyx damage in emergency room patients

FDR:

False discovery rate

Hb:

Hemoglobin

ICU:

Intensive care unit

IL-6:

Interleukin-6

IQR:

Interquartile range

MAP:

Mean arterial pressure

MFI:

Microvascular flow index

MVHS:

Microvascular health score

PBR:

Perfused boundary region

PCT:

Procalcitonin

PVD:

Perfused vascular density

RBC:

Red blood cell

RBCW:

Red blood cell width

SDF:

Sidestream dark field

SOFA score:

Sequential Organ Failure Assessment score

TVD:

Total vascular density

UKM:

University Hospital Münster

VRBC :

Red blood cell velocity

WBC:

White blood cell

References

  1. Ince C, De Backer D, Mayeux PR. Microvascular Dysfunction in the Critically Ill. Crit Care Clin. 2020;36(2):323–31. https://doi.org/10.1016/j.ccc.2019.11.003.

    Article  PubMed  Google Scholar 

  2. Pool R, Gomez H, Kellum JA. Mechanisms of Organ Dysfunction in Sepsis. Crit Care Clin. 2018;34(1):63–80. https://doi.org/10.1016/j.ccc.2017.08.003.

    Article  PubMed  Google Scholar 

  3. Rossaint J, Zarbock A. Pathogenesis of Multiple Organ Failure in Sepsis. Crit Rev Immunol. 2015;35(4):277–91. https://doi.org/10.1615/critrevimmunol.2015015461.

    Article  PubMed  Google Scholar 

  4. Bezemer R, Bartels SA, Bakker J, Ince C. Clinical review: Clinical imaging of the sublingual microcirculation in the critically ill–where do we stand? Crit Care. 2012;16(3):224. https://doi.org/10.1186/cc11236.

    Article  PubMed  PubMed Central  Google Scholar 

  5. De Backer D. Is microcirculatory assessment ready for regular use in clinical practice? Curr Opin Crit Care. 2019;25(3):280–4. https://doi.org/10.1097/MCC.0000000000000605.

    Article  PubMed  Google Scholar 

  6. Ince C, Boerma EC, Cecconi M, De Backer D, Shapiro NI, Duranteau J, Pinsky MR, Artigas A, Teboul JL, Reiss IKM, et al. Second consensus on the assessment of sublingual microcirculation in critically ill patients: results from a task force of the European Society of Intensive Care Medicine. Intensive Care Med. 2018;44(3):281–99. https://doi.org/10.1007/s00134-018-5070-7.

    Article  PubMed  Google Scholar 

  7. Pranskunas A, Koopmans M, Koetsier PM, Pilvinis V, Boerma EC. Microcirculatory blood flow as a tool to select ICU patients eligible for fluid therapy. Intensive Care Med. 2013;39(4):612–9. https://doi.org/10.1007/s00134-012-2793-8.

    Article  CAS  PubMed  Google Scholar 

  8. Massey MJ, Hou PC, Filbin M, Wang H, Ngo L, Huang DT, Aird WC, Novack V, Trzeciak S, Yealy DM, et al. Microcirculatory perfusion disturbances in septic shock: results from the ProCESS trial. Crit Care. 2018;22(1):308. https://doi.org/10.1186/s13054-018-2240-5.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Hutchings SD, Naumann DN, Hopkins P, Mellis C, Riozzi P, Sartini S, Mamuza J, Harris T, Midwinter MJ, Wendon J. Microcirculatory impairment is associated with multiple organ dysfunction following traumatic hemorrhagic shock: the MICROSHOCK study. Crit Care Med. 2018;46(9):e889–96. https://doi.org/10.1097/CCM.0000000000003275.

    Article  CAS  PubMed  Google Scholar 

  10. Scorcella C, Damiani E, Domizi R, Pierantozzi S, Tondi S, Carsetti A, Ciucani S, Monaldi V, Rogani M, Marini B, et al. MicroDAIMON study: Microcirculatory DAIly monitoring in critically ill patients: a prospective observational study. Ann Intensive Care. 2018;8(1):64. https://doi.org/10.1186/s13613-018-0411-9.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Vellinga NA, Boerma EC, Koopmans M, Donati A, Dubin A, Shapiro NI, Pearse RM, Machado FR, Fries M, Akarsu-Ayazoglu T, et al. International study on microcirculatory shock occurrence in acutely ill patients. Crit Care Med. 2015;43(1):48–56. https://doi.org/10.1097/CCM.0000000000000553.

    Article  PubMed  Google Scholar 

  12. Vellinga NAR, Boerma EC, Koopmans M, Donati A, Dubin A, Shapiro NI, Pearse RM, van der Voort PHJ, Dondorp AM, Bafi T, et al. Mildly elevated lactate levels are associated with microcirculatory flow abnormalities and increased mortality: a microSOAP post hoc analysis. Crit Care. 2017;21(1):255. https://doi.org/10.1186/s13054-017-1842-7.

    Article  PubMed  PubMed Central  Google Scholar 

  13. De Backer D, Donadello K, Sakr Y, Ospina-Tascon G, Salgado D, Scolletta S, Vincent JL. Microcirculatory alterations in patients with severe sepsis: impact of time of assessment and relationship with outcome. Crit Care Med. 2013;41(3):791–9. https://doi.org/10.1097/CCM.0b013e3182742e8b.

    Article  CAS  PubMed  Google Scholar 

  14. Watchorn JC, Fargaly H, Gilani M, Assadi J, Deitchman AR, Naumann DN, Wollborn J, Goebel U, McCurdy MT, Hutchings SD. The reproducibility of the point of care microcirculation (poem) score when used to assess critically ill patients: a multicenter prospective observational study. Shock. 2019. https://doi.org/10.1097/SHK.0000000000001474.

    Article  Google Scholar 

  15. Sardinha J, MacKinnon S, Lehmann C. Rapid clinical assessment of the sublingual microcirculation—visual scoring using microVAS in comparison to standard semi-automated analysis. Clin Hemorheol Microcirc. 2019;72(3):229–38. https://doi.org/10.3233/CH-180427.

    Article  PubMed  Google Scholar 

  16. Carsetti A, Aya HD, Pierantozzi S, Bazurro S, Donati A, Rhodes A, Cecconi M. Ability and efficiency of an automatic analysis software to measure microvascular parameters. J Clin Monit Comput. 2017;31(4):669–76. https://doi.org/10.1007/s10877-016-9928-3.

    Article  PubMed  Google Scholar 

  17. Hilty MP, Guerci P, Ince Y, Toraman F, Ince C. MicroTools enables automated quantification of capillary density and red blood cell velocity in handheld vital microscopy. Commun Biol. 2019;2:217. https://doi.org/10.1038/s42003-019-0473-8.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Drost CC, Rovas A, Kusche-Vihrog K, Van Slyke P, Kim H, Hoang VC, Maynes JT, Wennmann DO, Pavenstadt H, Linke W, et al. Tie2 activation promotes protection and reconstitution of the endothelial glycocalyx in human sepsis. Thromb Haemost. 2019;119(11):1827–38. https://doi.org/10.1055/s-0039-1695768.

    Article  PubMed  Google Scholar 

  19. Rovas A, Lukasz AH, Vink H, Urban M, Sackarnd J, Pavenstadt H, Kumpers P. Bedside analysis of the sublingual microvascular glycocalyx in the emergency room and intensive care unit - the GlycoNurse study. Scand J Trauma Resusc Emerg Med. 2018;26(1):16. https://doi.org/10.1186/s13049-018-0483-4.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Rovas A, Seidel LM, Vink H, Pohlkotter T, Pavenstadt H, Ertmer C, Hessler M, Kumpers P. Association of sublingual microcirculation parameters and endothelial glycocalyx dimensions in resuscitated sepsis. Crit Care. 2019;23(1):260. https://doi.org/10.1186/s13054-019-2542-2.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Hesse B, Rovas A, Buscher K, Kusche-Vihrog K, Brand M, Di Marco GS, Kielstein JT, Pavenstadt H, Linke WA, Nofer JR, et al. Symmetric dimethylarginine in dysfunctional high-density lipoprotein mediates endothelial glycocalyx breakdown in chronic kidney disease. Kidney Int. 2020;97(3):502–15. https://doi.org/10.1016/j.kint.2019.10.017.

    Article  CAS  PubMed  Google Scholar 

  22. Singer M, Deutschman CS, Seymour CW, Shankar-Hari M, Annane D, Bauer M, Bellomo R, Bernard GR, Chiche JD, Coopersmith CM, et al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA. 2016;315(8):801–10. https://doi.org/10.1001/jama.2016.0287.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Quan H, Li B, Couris CM, Fushimi K, Graham P, Hider P, Januel JM, Sundararajan V. Updating and validating the Charlson comorbidity index and score for risk adjustment in hospital discharge abstracts using data from 6 countries. Am J Epidemiol. 2011;173(6):676–82. https://doi.org/10.1093/aje/kwq433.

    Article  PubMed  Google Scholar 

  24. Dane MJ, Khairoun M, Lee DH, van den Berg BM, Eskens BJ, Boels MG, van Teeffelen JW, Rops AL, van der Vlag J, van Zonneveld AJ, et al. Association of kidney function with changes in the endothelial surface layer. Clin J Am Soc Nephrol. 2014;9(4):698–704. https://doi.org/10.2215/CJN.08160813.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Lee DH, Dane MJ, van den Berg BM, Boels MG, van Teeffelen JW, de Mutsert R, den Heijer M, Rosendaal FR, van der Vlag J, van Zonneveld AJ, et al. Deeper penetration of erythrocytes into the endothelial glycocalyx is associated with impaired microvascular perfusion. PLoS ONE. 2014;9(5):e96477. https://doi.org/10.1371/journal.pone.0096477.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Diez-Silva M, Dao M, Han J, Lim CT, Suresh S. Shape and biomechanical characteristics of human red blood cells in health and disease. MRS Bull. 2010;35(5):382–8. https://doi.org/10.1557/mrs2010.571.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Hulley SB CS, Browner WS, Grady D, Newman TB. Designing clinical research : an epidemiologic approach. In., 4th edn. Philadelphia, PA: Lippincott Williams & Wilkins; 2013: Appendix 6C, p. 79.

  28. Zhang CE, Staals J, van Oostenbrugge RJ, Vink H. Uncoupling of microvascular blood flow and capillary density in vascular cognitive impairment. Front Neurol. 2019;10:1268. https://doi.org/10.3389/fneur.2019.01268.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Krogh A. The supply of oxygen to the tissues and the regulation of the capillary circulation. J Physiol. 1919;52(6):457–74. https://doi.org/10.1113/jphysiol.1919.sp001844.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Angleys H, Ostergaard L. Krogh’s capillary recruitment hypothesis, 100 years on: Is the opening of previously closed capillaries necessary to ensure muscle oxygenation during exercise? Am J Physiol Heart Circ Physiol. 2020;318(2):H425–47. https://doi.org/10.1152/ajpheart.00384.2019.

    Article  CAS  PubMed  Google Scholar 

  31. Poole DC, Copp SW, Ferguson SK, Musch TI. Skeletal muscle capillary function: contemporary observations and novel hypotheses. Exp Physiol. 2013;98(12):1645–58. https://doi.org/10.1113/expphysiol.2013.073874.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Donati A, Damiani E, Domizi R, Romano R, Adrario E, Pelaia P, Ince C, Singer M. Alteration of the sublingual microvascular glycocalyx in critically ill patients. Microvasc Res. 2013;90:86–9. https://doi.org/10.1016/j.mvr.2013.08.007.

    Article  PubMed  Google Scholar 

  33. Donati A, Damiani E, Botticelli L, Adrario E, Lombrano MR, Domizi R, Marini B, Van Teeffelen JW, Carletti P, Girardis M, et al. The aPC treatment improves microcirculation in severe sepsis/septic shock syndrome. BMC Anesthesiol. 2013;13(1):25. https://doi.org/10.1186/1471-2253-13-25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Pranskunas A, Tamosuitis T, Balciuniene N, Damanskyte D, Sneider E, Vitkauskiene A, Sirvinskas E, Pilvinis V, Boerma EC. Alterations of conjunctival glycocalyx and microcirculation in non-septic critically ill patients. Microvasc Res. 2018;118:44–8. https://doi.org/10.1016/j.mvr.2018.02.004.

    Article  PubMed  Google Scholar 

  35. Wang G, Kostidis S, Tiemeier GL, Sol W, de Vries MR, Giera M, Carmeliet P, van den Berg BM, Rabelink TJ. Shear stress regulation of endothelial glycocalyx structure is determined by glucobiosynthesis. Arterioscler Thromb Vasc Biol. 2020;40(2):350–64. https://doi.org/10.1161/ATVBAHA.119.313399.

    Article  CAS  PubMed  Google Scholar 

  36. Harding IC, Mitra R, Mensah SA, Nersesyan A, Bal NN, Ebong EE. Endothelial barrier reinforcement relies on flow-regulated glycocalyx, a potential therapeutic target. Biorheology. 2019;56(2–3):131–49. https://doi.org/10.3233/BIR-180205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Wiesinger A, Peters W, Chappell D, Kentrup D, Reuter S, Pavenstadt H, Oberleithner H, Kumpers P. Nanomechanics of the endothelial glycocalyx in experimental sepsis. PLoS ONE. 2013;8(11):e80905. https://doi.org/10.1371/journal.pone.0080905.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Lipowsky HH, Gao L, Lescanic A. Shedding of the endothelial glycocalyx in arterioles, capillaries, and venules and its effect on capillary hemodynamics during inflammation. Am J Physiol Heart Circ Physiol. 2011;301(6):H2235-2245. https://doi.org/10.1152/ajpheart.00803.2011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to acknowledge the support of the Open Access Publication Fund of the University of Münster.

Funding

Open Access funding enabled and organized by Projekt DEAL. This work was supported by the German Research Foundation (KFO342) to AR, PK (KU 2873/3-1) and JR (RO 4537/5-1). The funding body had no role in the design of the study, collection, analysis and interpretation of data or in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

AR designed the study, performed the measurements, recorded clinical data, analyzed the data, prepared the figures, and contributed to the manuscript; JS recruited patients and reviewed the manuscript; SK retrieved clinical and microbiologic data and reviewed the manuscript; HP and JR contributed to the design of the study, discussed the findings and reviewed the manuscript; HV developed and refined the Glycocheck™ methodology, invented the novel variables, analyzed the data and contributed to the manuscript; PK had the initial idea, supervised the study, prepared the figures and drafted the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Philipp Kümpers.

Ethics declarations

Ethics approval and consent to participate

The study was performed in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the General Medical Council Westfalen-Lippe and the WWU Münster, Germany (file number: amendment to 2016–073-f-S). Informed consent was obtained.

Consent for publication

Consent for publication has been obtained. Written informed consent has been obtained from all identifiable individuals.

Competing interests

AR, JS, JR, SK, HP, and PK declare that they have no competing interests. HV is Chief Science Officer of GlycoCheck™ BV, The Netherlands. GlycoCheck™ and MVHS™ are trademarks registered by Microvascular Health Solutions LLC (Alpine, UT, United States).

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1

. Screenshots of randomly selected videos of one healthy (A-B) and four sepsis (C-J) individuals. Left: Screenshots without showing automatic vessel detection. Right: Vessels with diameter between 4 and 25 µm of the images on the left column are automatically highlighted and subjected to an automatic quality check. Invalid vascular segments are marked yellow and are automatically discarded, while all valid vascular segments (green lines) are further analyzed (see Figure 1).

Additional file 2: Video 1. Measurement of a healthy individual with and without automatic vessel detection and segments’ analysis.

Additional file 3: Video 2. Measurement of a sepsis patient with and without automatic vessel detection and segments’ analysis.

Additional file 4: Fig. S2

. Correlation between manually and automatically measurement of VRBC. Average manually derived longitudinal RBC movement per vessel was plotted against automatically measured VRBC of the same specific vessel (Spearman correlation). For manual validation of VRBC, 15 vessels (diameter 5 to 17 µm) from 3 randomly selected movies were analyzed manually. Therefore, longitudinal movement of individual RBCs in that vessels were tracked and measured as pixels per frame through subsequent vascular segments. Measured pixels per frame were transformed to µm/sec using camera frame-rate and pixel size.

Additional file 5: Fig. S3

. Derivation of capillary recruitment. The capillary recruitment is calculated by the slope of the relationship between VRBC (D ≤ 7 µm) and VRBC (D ≥ 10 µm). Two examples: In case capillary blood volume doubles when large vessel RBC velocity increases 2-fold, the slope (VRBC (D ≤ 7 µm), VRBC (D ≥ 10 µm)) will be 0 and CR = 1 – slope 0 = 1 = 100%. In the absence of changes in capillary blood volume when VRBC (D ≥ 10 µm) increases 2-fold, capillary RBC velocities are expected to also change proportionally by 2-fold, and the slope of VRBC (D ≤ 7 µm) vs. VRBC (D ≥ 10 µm) will be 1 and CR = 1 – slope 1 = 0 = 0%.

Additional file 6: Fig. S4

. Derivation of PBRdynamic from PBRstatic. All dots are shifted to V = 0 µm/sec parallelly to the slopes of each group.

Additional file 7: Table A1

. Comorbidities, focus of infection and isolated pathogens in the septic cohort. Table A2. Correlation coefficient between microvascular and clinical variables shown in Figure 1B, D, F. Table A3. Sensitivity analysis in sepsis patients and healthy controls. Table A4. Sensitivity analysis in sepsis patients.

Additional file 8: Fig. S5

. Capillary dropout in sepsis patients. Bar charts showing the percentage of loss of vascular density in sepsis patients compared to healthy controls. Q value (adjusted P value): *q < 0.05, **q < 0.01, ***q < 0.001

Additional file 9: Fig. S6

. Pooled PBRstatic values in predefined diameter ranges Boxplots of PBRstatic values of healthy controls and septic patients based on the different microvascular diameter ranges as output by the previous version of the GlyoCheck™ software used in [20].

Additional file 10: Fig. S7

. Conceptional development of the dynamic version of Microvascular Health score (MVHSdynamic)

Additional file 11: Fig. S8

. Association of MVHSdynamic with numbers of dysfunctional organs in sepsis patients (red) after dichotomizing (median) the group.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rovas, A., Sackarnd, J., Rossaint, J. et al. Identification of novel sublingual parameters to analyze and diagnose microvascular dysfunction in sepsis: the NOSTRADAMUS study. Crit Care 25, 112 (2021). https://doi.org/10.1186/s13054-021-03520-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13054-021-03520-w

Keywords