Clinical characteristics of recruited patients
A study outline is shown in Fig. 2. A total of 258 samples, consisting of BCs and ddPCR performed simultaneously, were collected from 228 patients from October 2022 to June 2023. In these samples, one contributed to two samples in 7 patients, 3 samples in 4 patients, and even more samples were collected from certain patients. As presented in Table 1, the median age of the patients was 78 years (IQR, 70–85 years), and 55.04% (146) were male. In terms of inflammatory indicators, the average plasma levels of CRP, Interleukine-6 (IL-6) and PCT were 61.56mg/L (IQR, 23.17–104.47 mg/L), 34.12 pg/L (IQR, 12.28–103.945 pg/L) and 0.47 ng/L (IQR, 0.106–2.28 ng/L), respectively. In view of coagulation function, the levels of fibrinogen and D-dimer were 4.47g/L (IQR, 3.24–5.75 g/L) and 2.1mg/L (IQR, 1.045–4.83 mg/L), respectively. Moreover, the platelet (PLT) was recorded as 187.05 ± 94.85 10^9/L. The severity of the disease was also assessed on Day 1, with the mean SOFA and APACHE II scores were 3.36 ± 3.27 and 13.01 ± 6.3, respectively. Among these patients with a cumulative 28-day mortality rate of 24.42%, 27.52% experienced acute kidney injury (AKI), 1.6% required renal replacement therapy (RRT), and 17.1% needed mechanical ventilation. Furthermore 15.1% received vasopressors and 66.7% were treated with combination antibiotic therapy. In addition, analysis of the 28-day survivors and non-survivors revealed no significant difference in the prevalence of hypertension and diabetes (P > 0.05). However, the survivors exhibited a younger age (P < 0.01), lower prevalence of AKI and coronary heart disease, and a reduced need of mechanical ventilation (P < 0.001), RRT (P = 0.018), vasoactive drug usage (P < 0.001), and immunosuppression (P < 0.001) when compared with non-survivors (Table 1).
Table 1
Clinical baseline characteristics of patients with BSIs
Characteristics
|
All
|
28-day survivors
|
28-day non-survivors
|
p value
|
number
|
258
|
195
|
63
|
-
|
General characteristics
|
|
|
|
|
Age,years
|
76.12 ± 12.32
|
74.43 ± 13.02
|
81.37 ± 7.90
|
< 0.001
|
Male,n(%)
|
142 ( 55.04 )
|
109 ( 55.9 )
|
33 ( 52.4 )
|
0.626
|
Comorbidities
|
|
|
|
|
Hypertension, n (%)
|
141 ( 54.65 )
|
108 ( 55.4 )
|
33 ( 52.4 )
|
0.677
|
Diabetes, n (%)
|
106 ( 41.09 )
|
83 ( 42.6 )
|
23 ( 36.5 )
|
0.396
|
Coronary heart disease, n (%)
|
149 ( 57.75 )
|
104 ( 53.3 )
|
45 ( 71.4 )
|
0.011
|
AKI n (%)
|
71 ( 27.52 )
|
39 ( 20.0 )
|
32 ( 50.8 )
|
< 0.001
|
Suspected infection site
|
|
|
|
|
Lower respiratory tract, n (%)
|
210 (81.4)
|
154 ( 79.0)
|
56 ( 88.9)
|
0.079
|
Urinary tract, n (%)
|
63 (24.42)
|
49 ( 25.1)
|
14 ( 22.2)
|
0.641
|
Intra-abdominal infection, n (%)
|
30 (11.63)
|
21 ( 10.8)
|
9 (14.3)
|
0.449
|
Skin and soft tissue, n (%)
|
15 ( 5.81 )
|
10 ( 5.1 )
|
5 ( 7.9 )
|
0.408
|
Abscess, n (%)
|
13 ( 5.04 )
|
12 ( 6.2 )
|
1 ( 1.6)
|
0.155
|
Clinical scores
|
|
|
|
|
MSS score,mean (SD)
|
2.81 ± 1.03
|
2.53 ± 0.77
|
3.68 ± 1.25
|
< 0.001
|
APACHEII score,mean (SD)
|
13.01 ± 6.3
|
11.22 ± 4.86
|
18.66 ± 6.97
|
< 0.001
|
SOFA score,mean (SD)
|
3.36 ± 3.27
|
2.31 ± 1.96
|
6.66 ± 4.23
|
< 0.001
|
MEWS score,mean (SD)
|
2.09 ± 1.48
|
1.66 ± 0.98
|
3.58 ± 1.89
|
< 0.001
|
Blood laboratory examination
|
|
|
|
|
WBC(10^9/L),mean (SD)
|
10.11 ± 6.32
|
9.34 ± 6.07
|
12.57 ± 6.54
|
< 0.001
|
Neutrophil percentage(%),mean (SD)
|
76.86 ± 12.38
|
74.92 ± 12.06
|
83.17 ± 11.25
|
< 0.001
|
Neutrophil cell count(10^9/L)
|
6.72 ( 3.9,11.15 )
|
5.98 ( 3.69,9.21 )
|
9.1 ( 6.23,14.25 )
|
0.045
|
Lymphocyte percentage(%)
|
12.6 ( 7.2,18.8 )
|
13.9 ( 7.9,20.1 )
|
8.6 ( 4.25,14.6 )
|
< 0.001
|
Lymphocyte cell count(10^9/L)
|
1.005 ( 0.66,1.435 )
|
1.04 ( 0.71,1.51 )
|
0.87 ( 0.465,1.275 )
|
0.178
|
PLT(10^9/L),mean (SD)
|
187.05 ± 94.85
|
195.75 ± 90.63
|
159.25 ± 103.16
|
0.008
|
RDW(%),mean (SD)
|
14.37 ± 2.79
|
14.10 ± 2.74
|
15.26 ± 2.80
|
0.004
|
Neutrophil Lymphocyte count ratio (NLCR)
|
6.02 ( 3.68,11.73 )
|
5.52 ( 3.22,10.38 )
|
9.7 ( 5.5,21.5 )
|
< 0.001
|
PLR
|
171.57 ( 113.43,248.69 )
|
172.28 ( 116.59,256.7 )
|
162.26 ( 90.96,243.1 )
|
0.678
|
RPR
|
0.08 ( 0.06,0.12 )
|
0.08 ( 0.06,0.1 )
|
0.11 ( 0.07,0.19 )
|
< 0.001
|
CRP(mg/L)
|
61.56 ( 23.17,104.47 )
|
55.27 ( 20.44,99.62 )
|
69.19 ( 34.83,144.97 )
|
0.003
|
IL-6(pg/mL)
|
34.12 ( 12.28,103.95 )
|
26.41 ( 10.86,76.92 )
|
75.32 ( 22.02,221.2 )
|
< 0.001
|
SAA(mg/L)
|
149.119 ( 56.82,288 )
|
151.67 ( 52.23,288 )
|
143.05 ( 62.41,219.12 )
|
0.670
|
PCT(ng/mL )
|
0.47 ( 0.11,2.28 )
|
0.28 ( 0.09,1.66 )
|
1.47 ( 0.49,5.97 )
|
0.694
|
HBP(ng/mL)
|
59.99 ( 27.33,150.64 )
|
59.84 ( 26.72,144.05 )
|
61.64 ( 29.95,156.68 )
|
0.860
|
Lac(mmol/L)
|
2 ( 1.6,2.6 )
|
1.9 ( 1.5,2.4 )
|
2.4 ( 1.8,3.4 )
|
0.616
|
SCr(umol/L),mean (SD)
|
101.56 ± 74.91
|
85.34 ± 43.39
|
152.39 ± 118.57
|
< 0.001
|
ESR(mm/h),mean (SD)
|
41.53 ± 24.35
|
43.45 ± 24.94
|
32.20 ± 19.34
|
0.103
|
Ferritin(ug/L)
|
492 ( 285,767 )
|
477 ( 184.5,778 )
|
541.5 ( 305.25,750 )
|
0.935
|
IL-1β(pg/mL)
|
2.5 ( 2.5,2.66 )
|
2.5 ( 2.5,2.68 )
|
2.5 ( 2.5,2.59 )
|
0.897
|
IL-2(pg/mL)
|
2.5 ( 2.5,2.5 )
|
2.5 ( 2.5,2.5 )
|
2.5 ( 2.5,2.5 )
|
0.370
|
IL-4(pg/mL)
|
2.5 ( 2.5,2.55 )
|
2.5 ( 2.5,2.5 )
|
2.5 ( 2.5,2.98 )
|
0.371
|
IL-5(pg/mL)
|
2.5 ( 2.5,2.5 )
|
2.5 ( 2.5,2.5 )
|
2.5 ( 2.5,2.5 )
|
0.193
|
IL-8(pg/mL)
|
17.17 ( 8.94,41.68 )
|
14 ( 7.49,30.5 )
|
37.67 ( 15.65,87.47 )
|
< 0.001
|
IL-10(pg/mL)
|
4.77 ( 2.97,8.82 )
|
4.44 ( 2.73,7.44 )
|
7.16 ( 4.18,14.68 )
|
0.253
|
IL-12P70(pg/mL)
|
2.5 ( 2.5,2.68 )
|
2.5 ( 2.5,2.54 )
|
2.5 ( 2.5,2.83 )
|
0.360
|
IL-17(pg/mL)
|
10 ( 10,14.74 )
|
10 ( 10,14.79 )
|
10 ( 10,14.09 )
|
0.270
|
IFN-α(pg/mL)
|
2.5 ( 2.5,2.5 )
|
2.5 ( 2.5,2.5 )
|
2.5 ( 2.5,2.5 )
|
0.909
|
IFN-γ(pg/mL)
|
2.57 ( 2.5,3.95 )
|
2.64 ( 2.5,3.99 )
|
2.5 ( 2.5,3.82 )
|
0.414
|
TNF-α(pg/mL)
|
2.5 ( 2.5,3.09 )
|
2.5 ( 2.5,3.17 )
|
2.5 ( 2.5,2.81 )
|
0.643
|
CD3(uL),mean (SD)
|
619.65 ± 381.86
|
641.89 ± 377.35
|
551.88 ± 391.95
|
0.181
|
CD4(uL),mean (SD)
|
398.83 ± 278.92
|
399.15 ± 261.83
|
397.86 ± 329.02
|
0.979
|
CD8(uL),mean (SD)
|
193.01 ± 136.07
|
210.23 ± 143.91
|
140.56 ± 91.79
|
0.003
|
CD4/CD8(%)
|
1.98 ( 1.31,3.63 )
|
1.9 ( 1.3,3.17 )
|
2.2 ( 1.6,4.4 )
|
0.231
|
CD19(uL)
|
106.5 ( 59.25,204 )
|
122 ( 61,217 )
|
92 ( 49,178 )
|
0.347
|
CD16 + CD56+(%)
|
157.5 ( 87.75,263 )
|
170 ( 92,260 )
|
107 ( 66,274 )
|
0.930
|
HLA-DR + CD3+/CD3+(%),mean (SD)
|
40.56 ± 21.05
|
41.84 ± 21.67
|
35.41 ± 17.93
|
0.212
|
Regulatory T cell (%),mean (SD)
|
3.16 ± 1.54
|
3.16 ± 1.58
|
3.16 ± 1.40
|
0.993
|
nCD64 index
|
5.84 ( 1.21,28.67 )
|
9.15 ( 1.14,32.88 )
|
1.77 ( 1.25,8.55 )
|
0.037
|
C1q(mg/L ),mean (SD)
|
157.09 ± 42.58
|
161.69 ± 41.82
|
141.10 ± 41.82
|
0.005
|
C3(mg/L )
|
0.98 ( 0.83,1.25 )
|
1.065 ( 0.85,1.28 )
|
0.88 ( 0.67,1.03 )
|
0.484
|
C4(mg/L ),mean (SD)
|
0.28 ± 0.13
|
0.29 ± 0.14
|
0.27 ± 0.13
|
0.489
|
IgG (g/L),mean (SD)
|
11.45 ± 3.52
|
11.54 ± 3.49
|
11.20 ± 3.64
|
0.624
|
IgA (g/L),mean (SD)
|
2.79 ± 1.29
|
2.69 ± 1.24
|
3.05 ± 1.39
|
0.157
|
IgM (g/L),mean (SD)
|
0.8 ± 0.45
|
0.86 ± 0.49
|
0.64 ± 0.28
|
0.012
|
IgE (g/L)
|
91.49 ( 23.78,210.38 )
|
94.81 ( 24.89,220.95 )
|
71.29 ( 21.22,167.83 )
|
0.774
|
NGAL(ng/mL)
|
87 ( 70.5,98.5 )
|
82 ( 62,91.5 )
|
842 ( 99,1114 )
|
0.636
|
TAT(ng/mL)
|
4.06 ( 2.56,6.84 )
|
3.75 ( 2.31,6.25 )
|
6.24 ( 3.51,7.88 )
|
0.079
|
tPAIC (ng/mL)
|
5.16 ( 3.09,8.24 )
|
4.88 ( 2.98,7.55 )
|
10.42 ( 4.06,25.46 )
|
< 0.001
|
TM(TU/mL),mean (SD)
|
15.04 ± 7.8
|
14.16 ± 7.01
|
20.43 ± 10.17
|
0.001
|
PIC (ug/mL )
|
1.33 ( 1,1.92 )
|
1.32 ( 1.00,1.93 )
|
1.33 ( 0.93,1.78 )
|
0.637
|
D-dimer (mg/L)
|
2.1 ( 1.05,4.83 )
|
1.68 ( 0.94,4.00 )
|
4.45 ( 2.17,8.39 )
|
0.004
|
Fibrinogen (g/L),mean (SD)
|
4.58 ± 1.73
|
4.77 ± 1.70
|
3.97 ± 1.69
|
0.002
|
APTT (g/L),mean (SD)
|
31.41 ± 8.42
|
30.32 ± 6.99
|
34.94 ± 11.31
|
< 0.001
|
PT (S),mean (SD)
|
14.34 ± 4.86
|
13.94 ± 4.56
|
15.62 ± 5.59
|
0.019
|
Clinical characteristics during hospitalization
|
|
|
|
|
Renal replacement therapy,n(%),mean (SD)
|
4 (1.6)
|
1 (0.5)
|
3 (4.8)
|
0.018
|
Use of vasoactive drugs,n(%),mean (SD)
|
39 (15.1)
|
8 (4.1)
|
31 (49.2)
|
< 0.001
|
Mechanical ventilation,n(%),mean (SD)
|
44 (17.1)
|
19 (9.7)
|
25 ( 39.7)
|
< 0.001
|
Combination antibiotic therapy,n(%),mean (SD)
|
172 (66.7)
|
115(59.0)
|
57 ( 91.9)
|
< 0.001
|
Immunosuppression,n(%),mean (SD)
|
76 (29.5)
|
38 ( 19.6)
|
38 (61.3)
|
< 0.001
|
Outcomes
|
|
|
|
|
Hospitalization expenses(RMB)
|
42275.04 ( 26497.07,86109.64 )
|
38839.34 ( 25408.35,58302.21 )
|
88737.45 ( 33113.53,193645 )
|
< 0.001
|
Antibiotic costs(RMB)
|
5127.63 ( 2643.1,10311.25 )
|
4678.245 ( 2575.923,7935.403 )
|
8952.88 ( 2857.26,45982 )
|
< 0.001
|
hospital stay, n (%)
|
14 ( 10,21.5 )
|
14 ( 10,19 )
|
13.5 ( 10,28.75 )
|
0.004
|
ICU days, n (%),
|
0 ( 0,4.5 )
|
0 ( 0,0 )
|
2.5 ( 0,13 )
|
< 0.001
|
Data were expressed as a mean ± standard deviation for normally distributed continuous variables or median (interquartile range) for non-normally distributed continuous variables. Categorical variables were expressed as n (%) |
AKI acute kidney injury; WBC white blood cell; PLT platelet; RDW red blood cell volume distribution width ;CRP C-reactive protein ;IL Interleukin; SAA Serumamyloid A ;PCT procalcitonin; HBP Heparin-Binding Protein ;ESR erythrocyte sedimentation rate; tPAIC tissue Plasminogen Tctivator-inhibitor Complex; TM thrombomodulin; PT prothrombinTime; APTT activated partial thromboplastin time; APACHE II Acute Physiology and Chronic Health Evaluation II; SOFA Sequential Organ Failure Assessment; MSS Modified Shapiro Score; ddPCR droplet digital PCR |
Performance of the ddPCR testing and the concordance between ddPCR and BC
In general, as illustrated in Table 2 and Fig. 3, the etiological diagnosis revealed that the ddPCR yielded 147 positive results from a total of 258 blood samples, with a positive rate of 56.98% (Fig. 3). Among them, bacteria accounted for 51.3%, viruses accounted for 43.2%, and 5.5% for fungi (Table 2). Of all the bacteria detected, the proportion of Gram-positive (G+) bacteria and Gram-negative(G−) bacteria were 50.4% and 49.6%, respectively. In contrast, BC only detected 18 positives. Among all pathogens detected by ddPCR, Streptococcus (n = 38) and EBV (n = 70) were the most frequently identified. Moreover, 60 G- bacteria were detected, with the top three strains being E. coli (n = 21), K. pneumoniae (n = 19), and A. baumannii (n = 14). Furthermore, the ddPCR assay revealed the presence of 61 G+ pathogens, with Streptococcus (n = 38), Enterococcus (n = 17), and Staphylococcus aureus (n = 6) being the predominant species. Additionally, among the remaining 13 strains, Candida (n = 5), Aspergillus (n = 6), and Pneumocystis yerbii (n = 2) were the most frequently detected fungi. As shown in Table 2, Table S2 and Fig. 3, results of BC and ddPCR were concordantly positive in 16 episodes with 13 identical pathogens and 3 different pathogens, concordantly negative in 106 episodes, and discordant in 220 episodes. In comparison with BCs, with the most common causative agents of culture-proven BSI being Escherichia coli (22.2%), Klebsiella pneumoniae (16.7%), Enterococcus (16.7%), and Staphylococcus aureus (5.56%), pathogens included in the ddPCR panel were identified in 88.8% (16 out of 18) of positive BCs (Table S2). Notably, the target detection listed by ddPCR did not encompass Clostridium fusiforme, Bacteroides fragilis, and Proteus mirabilis, which were identified through BC.
Table 2
Performance of ddPCR results for targeted organisms
|
BC+/ddPCR+,n
|
BC+/ddPCR-,n
|
BC-/ddPCR+,n
|
BC-/ddPCR-,n
|
Pathogens (all)
|
13
|
5
|
215
|
106
|
Klebsiella pneumonia
|
3
|
0
|
16
|
-
|
Acinetobacter baumannii
|
1
|
0
|
13
|
-
|
Escherichia coli
|
3
|
1
|
18
|
-
|
Pseudomonas aeruginosa
|
0
|
0
|
6
|
-
|
Staphylococcus aureus
|
1
|
0
|
5
|
-
|
Enterococcus
|
3
|
0
|
14
|
-
|
Streptococcus
|
1
|
0
|
37
|
-
|
EBV
|
0
|
0
|
70
|
-
|
CMV
|
0
|
0
|
25
|
-
|
VZV
|
0
|
0
|
2
|
-
|
HSV-1
|
0
|
0
|
5
|
-
|
Candida
|
1
|
0
|
4
|
-
|
Aspergillus
|
0
|
0
|
6
|
-
|
Pneumocystis yerbii
|
0
|
0
|
2
|
-
|
Clostridium fusiforme
|
0
|
1
|
0
|
-
|
Bacteroides fragilis
|
0
|
2
|
0
|
-
|
Proteus mirabilis
|
0
|
1
|
0
|
-
|
EBV epstein-barr virus; CMV cytomegalovirus; VZV varicella-zoster virus; HSV-1 Herpes simplex virus 1
When considering BSIs with comprehensive microbiological testing, the ddPCR testing demonstrated a sensitivity of 91.73% and a specificity of 81.6% (Table 3). Furthermore, the positive predictive value (PPV) and negative predictive value (NPV) were found to be 84.14% and 90.27%, respectively. Additionally, the sensitivity of ddPCR in detecting G+ bacteria, G− bacteria, and fungi was determined to be 96%, 90.7%, and 92.31%, respectively. The optimal diagnostic power for quantifying BSI through ddPCR is achieved with a copy cutoff of 166, which strikes a balance between sensitivity in detecting positive BSI patients and specificity in identifying case controls. The area under the receiver operating characteristic (AUROC) curves was determined to be 0.853 [95% confidence interval 0.756–0.951] (Table S3, Fig. 4). These preliminary data suggested that ddPCR had potential to rapidly identify targeted pathogens with high specificity and specificity.
Table 3
Sensitivity and Specificity of ddPCR in detecting different types of pathogens
|
sample(n = 258)
|
ddPCR+
|
ddPCR-
|
Sensitivity(%)
|
Specificity(%)
|
PPV(%)
|
NPV(%)
|
Total(MSS ≥ 2)
|
Positive by all microbiological testing
|
122
|
3
|
91.73
|
81.60
|
84.14
|
90.27
|
|
Negative by all microbiological testing
|
23
|
110
|
|
|
|
|
G-
|
Positive by all microbiological testing
|
48
|
2
|
96.00
|
53.37
|
33.10
|
98.23
|
|
Negative by all microbiological testing
|
97
|
111
|
|
|
|
|
G+
|
Positive by all microbiological testing
|
39
|
4
|
90.70
|
50.70
|
26.90
|
96.46
|
|
Negative by all microbiological testing
|
106
|
109
|
|
|
|
|
Fungi
|
Positive by all microbiological testing
|
12
|
1
|
92.31
|
45.71
|
8.28
|
99.12
|
|
Negative by all microbiological testing
|
133
|
112
|
|
|
|
|
Virus
|
Positive by Virus antibody
|
37
|
18
|
67.27
|
46.80
|
25.52
|
84.07
|
|
Negative by Virus antibody
|
108
|
95
|
|
|
|
|
ddPCR droplet digital PCR; PPV positive predictive value; NPV negative predictive value; G- gram-negative bacteria; G + gram-positive bacteria
In addition to pathogen identification, the AMR genes panel was utilized to identify seven AMR genes, namely blaKPC, mecA, blaOXA-48, blaNDM, blaIMP, vanA, and vanM. However, only the blaKPC, mecA, and blaNDM genes were found to be positive through ddPCR testing, as shown in Table 4. The ddPCR analysis revealed that there were 5 episodes with a positive result for blaKPC and 1 for blaNMP. Among these episodes, the simultaneous detection of K. pneumoniae and the AMR gene occurred in 31.6% of cases, which held significant clinical implications. In comparison to the results obtained from the BC, it was observed that two instances of K. pneumoniae reported in the BC exhibited resistance towards carbapenems. Furthermore, the ddPCR results indicated the presence of the bla KPC gene in these strains. Notably, the mecA positive sample and blaNDM positive sample were not subjected to pathogen testing. From a therapeutic respective, the identification of drug resistance genes within a span of three hours facilitated the selection of sensitive antibiotics for the target pathogen as determined by the initial day ddPCR assay. Consequently, the patient’s condition exhibited gradual improvement, accompanied by a reduction in both pathogen load and AMR gene load (Table S4).
Table 4
AMR genes detected by ddPCR and the related pathogens
AMR genes
|
Pathogens
|
counts
|
blaKPC(n = 5)
|
Klebsiella pneumoniae
|
5
|
blaNMP(n = 2)
|
Klebsiella pneumoniae
|
1
|
|
none
|
1
|
mecA (n = 10)
|
Staphylococcus aureus
|
3
|
|
none
|
7
|
AMR genes antimicrobial resistance genes.
Clinical potential value of ddPCR for Epstein-Barr virus (EBV)infection
One hundred and two positive viral pathogens including EBV, CMV, VZV and HSV-1 were detected through the ddPCR, with a positive rate of 43.2% (Table 2). Interestingly, EBV was the most frequently identified virus in our study. As shown in Table 5, a total of 258 episodes from 228 patients with BSIs underwent testing for EBV antibody and mcfDNA using ddPCR. Of these, 70 (27.13%) tested positive for EBV reactivation. Among these 70 cases, 20 (28.57%) were found to have concurrent COVID-19 infection and immunosuppression (Table S5). When it comes to EBV antibody, we found that the EBV antibodies (VCA-IgM, VCA-IgG, and EBNA-IgG) in blood were related to copy number of ddPCR in BSI patients. Our results indicated that the copy number of ddPCR with VCA-IgM negative, VCA-IgG negative, and EBNA-IgG negative in blood was significantly higher than that of VCA-IgM negative, VCA-IgG positive, and EBNA-IgG positive. Consequently, the group characterized by VCA-IgM negative, VCA-IgG negative, and EBNA-IgG negative was considered to exhibit immunologic unresponsiveness associated with immunosuppression (Table 5). The correlation analysis conducted between the number of EBV copies as determined by ddPCR and immune indicators revealed a statistically significant correlation between the EBV copy number and the CD4+/CD8 + ratio (r = -0.312, p = 0.029) (Fig. 5).
Table 5
EBV detected by ddPCR assay and related comorbidities VCA viral capsid antigen; EBNA Epstein Barr Virus Nuclear Antigen; COVID-19 Corona Virus Disease 2019
Anti-EBV Antibodies
|
number
|
Explaination
|
EBV(copy number)
|
COVID-19
|
Immunosuppression
|
VCA-IgM
|
VCA-IgG
|
EBNA-IgG
|
negative
|
negative
|
negative
|
44
|
No immune response
|
4671.48 ± 13195.01*
|
13 (29.55)
|
18 (40.91)*
|
positive
|
negative
|
negative
|
1
|
Acute infection
|
1472
|
0
|
0
|
negative
|
positive
|
positive
|
23
|
Previous infection
|
611.87 ± 876.21*
|
7 (30.43)
|
2 (8.70)*
|
negative
|
positive
|
negative
|
2
|
Acute infection/Previous infection
|
19515 ± 27431.50
|
1
|
1
|
VCA viral capsid antigen; EBNA Epstein Barr Virus Nuclear Antigen; COVID-19 Corona Virus Disease 2019
Clinical potential value of ddPCR for liver abscess
In our research, it was observed that ddPCR exhibited a high level of sensitivity and specificity in detecting liver abscess in patients (Supplementary Table S6). Furthermore, the ratio of pathogens detected by ddPCR to those detected by pus culture was found to be 100% (7 out of 7 cases). Among these cases, Klebsiella pneumoniae was detected in 5 episodes, while Escherichia coli was detected in two. In stark contrast BC did not yield any relevant pathogenic bacteria. Moreover, liver abscess patients underwent antimicrobial de-escalation therapy following negative results obtained from the ddPCR assay conducted on the third day. Fig. S1 displays scatter plots of liver abscess representative chip analysis results from a clinical case that was dynamically examined and clinically improved after antibiotic treatment.
Correlative analysis between ddPCR and the biomarkers
Sepsis develops as a consequence of a complicated, dysregulated host response to infection, which is characterized not only by increased inflammation but also mainly by abnormal coagulation function as well as immune suppression. To explore whether the correlation exists between the copies of pathogens identified through ddPCR in BSIs and various markers of inflammation, coagulation, and immunity, spearman correlation coefficient was utilized to describe the relationship. The results are presented in Fig. 5. From the perspective of inflammatory markers, we observed the correlation between the levels of following inflammatory indicators and the pathogen load identified by ddPCR (PCT: Spearman’s rho = 0.309, P<0.001; CRP: Spearman’s rho = 0.242, P = 0.004; Neutrophil Lymphocyte count ratio (NLCR): Spearman’s rho = 0.221, P = 0.009; White blood cell (WBC): Spearman’s rho = 0.254, P = 0.002; Neutrophil percentage: Spearman’s rho = 0.294 P<0.001; Neutrophil cell count: Spearman’s rho = 0.242 P = 0.004; Lymphocyte percentage: Spearman’s rho = -0.196 P = 0.021), while no relationship was found with lymphocyte count. Additionally, a correlation was discovered between the PCT level and the pathogen load of G− bacteria detected by ddPCR (Spearman’s rho = 0.589, P < 0.001). Furthermore, Klebsiella pneumoniae exhibited a correlation with PCT (Spearman’s rho = 0.757, P < 0.001) as well as blaKPC (Spearman’s rho = 0.928, P < 0.01). However, no correlation was found between the pathogen load with heparin-binding protein (HBP) or Serumamyloid A (SAA). Cytokines are important indicators of inflammation, indeed, when considering copies of pathogens, a positive correlation was found with cytokines such as IL-6, IL-8 and IL-10, with correlation coefficients of 0.192, 0.241 and 0.240, respectively. For other cytokines like Interleukin-1β (IL-1β), IL-17, Tumor Necrosis Factor-a (TNF-a), IL-4 and IL-5 were not found to be related to the pathogen load in our research. From the view point of coagulation, the ddPCR assay revealed a significant positive correlation between the copies of pathogens detected and coagulation parameters, including tissue Plasminogen Activator-inhibitor Complex (tPAIC) (Spearman’s rho = 0.421, P < 0.001), thrombomodulin (TM) (Spearman’s rho = 0.364, P < 0.01), D-dimer (Spearman’s rho = 0.271, P < 0.001), ProthrombinTime (PT) (Spearman’s rho = 0.248, P < 0.01), and Activated Partial Thromboplastin Time (APTT) (Spearman’s rho = 0.291, P < 0.01). In the respect of the correlation between the copies of pathogens detected and immunity indicators, it was observed that copies of pathogens were not related to cellular immunity-related markers as CD3, CD4, CD8, Treg, Human Leukocyte Antigen-DR (HLA-DR)/CD3 and humoral immunity-associated proteins as CD19 and IgG, IgA and IgM. We also found no correlations between copies of pathogens and complement system biomarkers as C1q, C3 and C4. As a widely used biomarker reflecting the severity of sepsis, lactate was found to be related to the pathogen load detected by ddPCR (Spearman’s rho = 0.19, P < 0.05).
ddPCR performance for 28-day prognosis
The aforementioned results revealed that pathogens identified through ddPCR was closely related to indicators reflecting the severity of infection (Fig. 5). Subsequently, correlation analysis between severity of illness score and pathogen load detected by ddPCR was conducted, and they exhibited a positive correlation (SOFA: Spearman's rho = 0.322, P<0.001; APACHE II: Spearman's rho = 0.217, P<0.01; MEWS: Spearman's rho = 0.244, P<0.01) (Fig. 5). To assess the performance of ddPCR as a continuous metric in comparison to other established prognostic biomarkers, we conducted calculations of AUROCs for 28-day mortality. As presented in Table S7 and Fig. 6, our findings indicated that ddPCR exhibited an AUROC of 0.714 (95% CI, 0.621–0.807) for the identification and prediction of 28-day prognosis. The result demonstrated a sensitivity of 63.4% and a specificity of 72.3%, with a derived cut-off of 1046 copies/mL. The corresponding Kaplan-Meier curves and the outcomes of log-rank tests for ddPCR copy numbers above or below 1046 copies/mL were also presented. The results of univariate analysis for patients with BSIs who survived for 28 days (n = 195) and those who did not survive for 28 days (n = 63) in the development cohort are presented in Table S8. Multivariate analysis revealed that NLCR, C1q, PLT, and D-dimer were identified as independent risk factors for the 28-day mortality in patients with BSI in the development cohort (Table 6). A nomogram was constructed based on the aforementioned equation. The calibration plot of the nomogram demonstrated a satisfactory fit within the development cohort, when considering the predicted probability or the actual probability. Additionally, it exhibited strong statistical consistency in predicting the 28-day mortality caused by BSI, as evidenced by a C value of 0.796. In comparison with the outcomes observed, the nomogram displayed a sensitivity of 74.4% and a specificity of 76.2% in predicting 28-day mortality (Fig.S2). As we had found that the ddPCR was related to the severity of disease, we further assessed the ICU time and hospital stays between positive and negative groups. To our surprise, we found no statistical difference between two groups (Fig.S3).
Table 6
Logistic regression analyses of factors associated with 28-day prognosis focus on ddPCR
|
p
|
HR
|
95.0% CI
|
|
Lower
|
Upper
|
copy number
|
0.009
|
4.94
|
1.487
|
16.415
|
Use of vasoactive drugs
|
0.022
|
13.487
|
1.467
|
123.996
|
MEWS score
|
0.001
|
2.666
|
1.503
|
4.727
|
NLCR
|
0.012
|
1.067
|
1.014
|
1.122
|
C1q(mg/L )
|
0.001
|
1.036
|
1.014
|
1.058
|
PLT(10^9/L)
|
0.022
|
1.006
|
1.001
|
1.011
|
D-dimer (mg/L)
|
0.048
|
1.085
|
1.001
|
1.176
|
RDW(%)
|
0.001
|
1.718
|
1.234
|
2.393
|
MEWS Modified Early Warning Score;NLCR neutrophil lymphocyte count ratio ;PLT platelet; RDW red blood cell volume distribution width
Health economic evaluation of the ddPCR
The ddPCR assay had several advantages in health economic evaluation which is exhibited in Fig. 7. Based on the results of the microbiological test and clinical assessment, patients diagnosed with BSI were categorized into negative and positive groups using ddPCR assay. Notably, the negative groups identified through ddPCR exhibited comparatively lower hospitalization expenses when compared to the positive groups identified through ddPCR (P < 0.05). When it comes to the antibiotics cost, negative groups were significantly lower than that of positive groups (P < 0.001). In addition, we evaluated the percentage of antibiotics costs within the total hospitalization expenses, with the result that negative groups was comparatively lower than the positive groups (P < 0.001).