Our search algorithm in the MEDLINE, OVID, PubMed and MedRXIV databases retrieved 1441 records. We screened based on the title and abstract and excluded 1377 studies. Sixty-four articles were assessed for eligibility. Finally, we identified 29 observational cohort studies that qualified according to our predefined inclusion and exclusion criteria (Figure 1), which included a total of 2095 patients with COVID-19 admitted to an ICU.
Characteristics of included studies
Of the 29 studies, 8 (27.59%) were from France [2, 13, 21-26], 4 (13.79%) were from the United Kingdom (UK) [27-32], 4 (13.79%) were from the Netherlands [33, 34], 2 (6.90%) were from Germany [35, 36], 2 (6.90%) were from Spain [15, 37], 2 (6.90%) were from Switzerland [14, 38], 1 (3.45%) was from the USA [39], 1 (3.45%) was from Belgium [40], 1 (3.45%) was from China[6], 1(3.45%) was from Pakistan [41], 1(3.45%) was from Denmark [42], 1 (3.45%) was from Italy [43] and 1(3.45%) from Mexico [44]. Twelve studies dealt with patients who received mechanical ventilation, eight studies were of patients with ARDS, one study was on extracorporeal membrane oxygenation (ECMO), 1 study looked at acute respiratory failure (ARF), and seven studies included other patients with severe COVID-19 (See additional file 1, supplementary 1).
Quality control of included studies
The risk of bias included in the cohort studies was assessed using the Newcastle‑Ottawa-Scale (NOS) as presented in Table 1. The overall score was 141 out of 232 (60.78%), which is considered to be indicative of moderate quality. All the eligible studies included patients with severe COVID-19 admitted to ICU. Among the studies, 34.48% (10/29) did not describe the determination of exposure. In addition, 72.41% (21/29) studies did not control the comparability of the cohorts on the basis of the design or analysis. Although all studies (29/29) followed up the patients for long enough for defined outcomes to occur, only 51.72% (15/29) demonstrated clearly the adequacy of follow up of the cohorts. All studies assessed the outcomes of interest based on medical record linkage (See additional file 1, supplementary 2).
Incidence of CAPA
In the studies, 2095 patients with severe COVID-19 admitted to an ICU were investigated and 264 cases of CAPA were reported. The incidence of CAPA was calculated as 0.14 (95% confidence interval [CI] = 0.11–0.17, I2 = 81.2%). In the sub‑group analysis by study design, the incidence of CAPA was calculated at 0.14 (95% CI = 0.09–0.19, I 2= 87.4%) in the prospective studies and 0.14 (95% CI = 0.10–0.18, I2 = 74.1%) in the retrospective studies (Table 1). When stratified analysis was undertaken by research centre (multi-centre or single centre), the CAPA incidence was 0.12 (95% CI = 0.10–0.15, I2 = 77.3%) for the single centre studies and 0.20 (95% CI = 0.12–0.28, I2 = 81.2%) for the multi-centre studies. We found a high heterogeneity for calculation of the incidence of CAPA and significant publication bias according to Egger’s regression test (See additional file 1, supplementary 3 Figure A). Remarkably, we found that patients with COVID-19 admitted to an ICU would develop CAPA after approximately 7.28 days of mechanical ventilation (nine studies, 95% CI = 5.48-–9.08, I2 = 46%) (Figure 2).
Colonisation by Aspergillus
In five cohort studies including 365 patients with COVID-19 admitted to an ICU, 21 patients were assessed for Aspergillus colonisation. The calculated colonization rate of was 5.75% (21/365) (See additional file 1, supplementary 1), suggesting that clinicians should comprehensively analyse Aspergillus colonisation combined with other clinical evidence (such as imaging and inflammatory factors, etc.) to judge whether the patient might develop CAPA and require further anti-fungal treatment.
All-cause mortality of CAPA
Twenty-three studies were included in the final analysis, with 119 deceased patients and 2786 surviving patients. The overall mortality of CAPA was 0.07 (95% CI = 0.05–0.09, I2 = 64.2%) among patients with severe COVID-19. In the subgroup analysis for study design, the calculated all-cause mortality of CAPA was 0.07 (95% CI = 0.04–0.10, I2 = 67.1%) for the prospective studies, 0.08 (95% CI = 0.05-0.11, I2 = 62.9) for the retrospective studies, 0.06 (95% CI = 0.09–0.09, I2 = 0%) for the single centre studies, and 0.11 (95% CI = 0.06–0.16, I2 = 0%) for the multi-centre studies. There was significant publication bias on Egger’s regression test (Table 1).
CFR of CAPA
Twenty-three studies were eligible for further analysis of all-cause CFR with 119 deceased patients and 145 surviving patients. We found that the pooled CFR was 0.51 (95% CI, 0.44–0.58, I2 = 86.5%). In the subgroup analysis by study design, the calculated CFRs were 0.49 (95% CI, 0.39–0.59, I2 = 88.4) for the prospective studies, 0.53 (95% CI = 0.43–0.63, I2 = 85.9%) for the retrospective studies, 0.50 (95% CI = 0.42–0.58, I2 = 86.5%) for the single centre studies, and 0.58 (95% CI = 0.42–0.73, I2 = 89.5%) for the multi-centre studies (Table 1).
Comparing patients with CAPA with those without CAPA (control)
Risk factors for CAPA
Eleven studies were pooled to investigate and compared patients with and without CAPA (127 vs. 788 patients) (Table 2). No significant differences were observed regarding general population characteristics (age, sex, hypertension, diabetes, obesity, cardiovascular disease, history of cancer, chronic pulmonary disease, chronic kidney disease, an immunosuppression), blood examination (white blood cell [WBC] count, neutrophil count, lymphocyte count, creatinine, C-reactive protein [CRP], and lactate dehydrogenase (LDH)), and respiratory sample characteristics at baseline (Table 2). Compared to patients without CAPA (control), patients with CAPA had a significantly lower median BMI (four studies, 27.32 vs. 28.97 kg/m2, P = 0.034) and higher median creatinine (four studies, 127.94 vs. 88.23, P = 014). No differences in the use of Lopinavir/ritonavir, Hydroxychloroquine (HCQ), Azithromycin, Tocilizumab, ventilated prone, and ECMO were observed in the patients with CAPA compared with those without (Table 2). However, compared with the controls, patients with CAPA were more likely to have received corticosteroids during admission (eight studies, 41.0% [41/100] vs. 38.0% [273/719], risk ratio [RR] = 1.98, 95% CI = 1.08–3.63, I2 = 74.10%) and renal replacement therapy (two studies, 42.0% [29/69] vs. 28.2% [90/319], RR = 1.61, 95% CI = 1.04–2.50, I2 = 32.5) (Table 2).
Outcome with CAPA or without CAPA
Eleven studies were included. The analysis showed that patients with CAPA were associated significantly with a 1.66-fold higher risk for mortality (risk ratio [RR] = 1.66, 95% CI = 1.31–2.12, I2 = 22.5%) (Figure 3) without significant heterogeneity and publication bias (See supplementary 3). In the subgroup analysis by study design, patients with CAPA were associated significantly with a 1.53-fold elevated risk of mortality (six studies, RR = 1.53, 95% CI = 1.02–2.21, I2 = 47.9) for the retrospective cohort studies, a 1.93-fold elevated risk of mortality (five studies, RR = 1.93, 95% CI = 1.37–2.71, I2 = 0%) for the prospective cohort studies (Figure 3 A), a 1.52-fold elevated risk of mortality (nine studies, RR = 1.52, 95% CI = 1.12–2.07, I2 = 24.7%) the for single centre studies, and a 2.03-fold elevated risk of mortality (two studies, RR = 2.03, 95% CI = 1.46–2.81, I2 = 0%) for the multi-centre studies (Figure 3 B). No statistically significant difference was found in the length of ICU stay for CAPA and non-CAPA patients (MD = 5.58, 95% CI = -1.93 to 13.08, P = 0.145) (Table 2).