Global and regional prevalence and outcomes of COVID-19 in people living with HIV: A cutting-edge systematic review and meta-analysis

Background The relationship between HIV infection and COVID-19 clinical outcome is uncertain, with conflicting data and hypotheses. We aimed to assess the prevalence and risk of severe COVID-19 and death in people living with HIV (PLWH) on the global and continental level. Methods Electronic databases were systematically searched in July 2021. Studies were screened and then extracted following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. Narratives were synthesised and data pooled for global and continental prevalence and relative risk of severity and mortality in HIV-infected COVID-19 patients using random-effect model. Risk of bias was assessed using the Newcastle-Ottawa score, Egger’s test and presented as funnel plots. Results A total of 46 studies were included involving 18,034,947 COVID-19 cases of which 31,269 were PLWH. The global prevalence of PLWH with SARS-CoV-2 infection was 1% (95% CI = 0.9% - 1.1%) with the highest prevalence observed in sub-Saharan Africa. The relative risk (RR) of COVID-19 severity was significant only in Africa (RR, 95% CI = 1.14, 1.08 – 1.24) while risk of COVID-19 mortality was 1.53% (95% CI = 1.45 – 2.03) globally. The prevalence of PLWH in COVID-19 cases was significantly low, and the calculated global risk ratio show that HIV infection may be linked with increased COVID-19 death. The between-studies heterogeneity was significantly high while risk of publication bias was not significant.


Introduction
The 2019 coronavirus  pandemic caused by SARS-CoV-2 remains a global public health challenge that has affected over 186 million people and caused over 4 million deaths globally (1).
While most cases of COVID-19 are clinically mild or asymptomatic, older age and certain underlying illness such as cardiovascular, respiratory, and digestive diseases have been reported to increase risk of severe COVID-19 cases or death (2)(3)(4). Such comorbidities are associated with an increased fatality rate and present a challenge for intensive care management of COVID-19 patients (5,6).
Human Immunodeficiency Virus (HIV) belongs to a genus of zoonotic lentiviruses that cause acute immune-deficiency syndrome (AIDS) (7). Data from the Joint United Nations Programme on HIV/AIDS (UNAIDS) puts the number of people living with HIV (PLWH) at 38 million globally with 1.5 million new infections in 2020 and about 6 million not knowing their HIV infection status (8). Accordingly, the number of PLWH is projected to increase due to treatment availability and the associated reduction in AIDS-related deaths (9).
HIV is associated with a dysregulation of the immune system which predisposes patients to opportunistic infectious diseases (10). Indeed, most HIV-related death have been linked to secondary infections and abnormal inflammatory response resulting from AIDS (11) This is especially so in patients with uncontrolled HIV replication, high viral load and low CD4/CD8 count. Giving the immune-compromised state of most PLWH and increased possibility of secondary dysfunctions, an increased risk of infection, severity, and death due to COVID-19 may be expected. However, attenuated immune response may also protect against cytokine release storm and corresponding acute respiratory distress syndrome (ARDS) linked with SARS-CoV-2 infection (12). Hence, the association between HIV infection and COVID-19 remain unclear with sparse and conflicting reports (13,14). Aside from the heterogeneity from difference in epidemiology of HIV between countries and . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint continents, variability exist in treatment, management, and behaviour of PLWH. These amongst other factors, determines the rate of spread as well as the availability and uptake of HIV preventative and treatment measures (15). This is further exacerbated by the disruption to clinical care of various chronic diseases due to the diversion of medical resources to manage the burgeoning COVID-19 cases around the world (16). This review aimed to provide an updated insight into the global and continental prevalence, potential risk of severe illness and death due to COVID-19 in PLWH by conducting a meta-analysis grouped by continents.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint

Methods
The protocol of this systematic review was registered prospectively to PROSPERO (CRD42021264151). Following the Preferred Reporting in Systematic Reviews and Meta-Analyses (PRISMA) guidelines (17), Medline and Embase databases were searched on 2 nd of July 2021 using keywords and MeSH terms (Fig. S1). Further, a search of pre-printed papers was performed because of the rapidly developing nature of the topic. Studies retrieved from the search were imported into EndNote software and the duplicates removed. The resulting duplicated-free studies were then uploaded to Rayyan software and title, abstract and whole text screening was carried out (RPR, AK).
Manual reference search was also performed to cover papers that were not covered by the search strategy ( Fig. 1).

➢ Inclusion and exclusion criteria
The inclusion/exclusion of studies included followed the PECO (Population, Exposure, Comparison and Outcome) model and only involve studies that presents the clinical characteristics and/or composite endpoints of COVID-19 in people living with HIV (PLWH). Studies that combined HIV and other immunosuppressive diseases and conditions (cancer, congenital, or medically induced), non-English language, reviews, case reports, qualitative studies, editorials and were excluded.

➢ Data collections
Two authors (RPR, AK) independently screened titles and abstracts of potentially eligible studies and any disagreements were resolved by discussing with a third reviewer (TO). Full texts of the studies that were included based on the abstract were fully read and assess against the inclusion/exclusion criteria.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint

➢ Data extraction and analysis
The authors and year of publication, study design and period, country of study, sample size, number of PLWH and non-HIV COVID-19 patients as well as their clinical outcomes and biometrics data were extracted into a table. Clinical outcomes identified were severity of the disease and death. Severe COVID-19 was described as prolonged hospital stay, ICU admission and requirement of mechanical ventilation (MV). Analyses were performed using Stata/MP 17; prevalence was calculated by "metaprop" procedure. Forest plots were used to present the pooled prevalence of PLWH in COVID-19 cases grouped by continents of studies. Continent-grouped effect sizes (95% Confidence Intervals, CIs) and the test results of between-studies heterogeneity (I 2 Statistic, p value) were also computed.
The 'metan' procedure was used to assess risk of severity and mortality in PLWH-COVID-19 patients compared with general population in the included studies and the risk ratios were grouped by continents to further assess the intercontinental variation in these risks.

➢ Quality assessment
A modified version of the Newcastle-Ottawa Score (NOS) was used to assess the risk of bias in the included studies (18). This includes 3 domains and 9 questions scored accordingly with a star if satisfied. The selection domain assessed the randomness, multicentre involvement in selection of study population as well as sample size. Sample size of ≥ 100 was decided based on previous studies' estimate of ~1% prevalence of HIV infection in COVID-19 cases (19,20). The standard ascertainment of COVID-19 and HIV were also assessed against the WHO guideline (21). Finally, the follow up time (≥2 weeks), mode of outcome confirmation and whether all patients were accounted for were assessed (Table S1). Studies with ≥5 stars (> 50%) were considered unbiased. To further assess publication bias in the studies pooled for prevalence and the risk of severity and mortality in PLWH-COVID-19, funnel plots and Egger test were computed using 'metafunnel' and 'metabias' procedures respectively in STATA. Statistical significance was set at 95% (p < 0.05).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint

Identification
Screening . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Results
The systematic search of databases including preprints and reference search generated an initial total of 966 studies including 245 duplicates to give a total of 721 studies. Initial title and abstract screening led in the exclusion of 652 studies followed by full-text review and further exclusion of 11 studies to give a total of 46 studies which satisfy the inclusion/exclusion criteria (Fig. 1).  (Table 1). The risk of bias assessment showed low bias in the included studies with 87% (40/46) of the studies below the bias threshold (Table S1).

➢ General description of studies included
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

➢ Prevalence of PLWH co-infected with Sars-CoV-2
Of the 46 studies included, 10 studies reported only PLWH-COVID-19 and were not included in the analysis (30,33,40,48,54,63). Another study in which patients were not randomly sampled (61) was also excluded from the analysis to reduce bias. Twelve of these studies were conducted in Europe  Fig 2b).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint  Fig 3a). This result holds true for both Europe and USA both of which are . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint associated with better prevention and management of HIV. However, the risk for severe  in PLWH from Africa was found to increase by 14% (RR (95%CI) = 1.14 (1.05 -1.24). Also, while the overall heterogeneity was significantly high (85%, p < 0.001), there was no between study variation in the studies from Africa (I 2 = 0%, p = 0.43). Indeed, funnel plot and Egger's test showed no publication bias in all included studies (T(95%CI) = -1.32(-3.02 -0.75), p = 0.21; Fig 3b). . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint  1.52(1.14 -2.03); Fig 4a). On the regional level, there was no significant increased risk of COVID-19 mortality in PLWH in Africa or Europe. But . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint a 2-fold increase in severity was observed in the USA for the studies included. Despite this difference in regional risk ratios, the computed funnel plot and Egger's test showed no publication bias in the included studies (T(95%CI) = -0.44 (-1.58 -2.41) ; Fig 4b). . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Discussion
This study provides an updated systematic assessment of the prevalence and clinical outcomes of COVID-19 in PLWH on both the global and regional level. To the best of our knowledge, this is a first systematic review and meta-analysis dedicated singly to the understanding of the clinical outcome of PLWH co-infected with SARS-CoV-2 in both global level and within continents with consideration of the variation in the prevalence, prevention, and management of HIV infections. We found a significantly low global prevalence of PLWH in COVID-19 cases (1%) compared to other comorbidities such as COPD, cardiovascular, chronic liver, and kidney disease (2,3,70,71). This is consistent with previous reports which estimate the prevalence of HIV coinfection at 1-2% of COVID-19 patients either admitted to the hospital or in the general population (72). Further, while the proportion was below 2% in Europe, North America and Asia, the prevalence of HIV infection in COVID-19 cases was high in Africa (11%). This is reflective of the global epidemiology of PLWH whereby more than half of the global cases are located within continental Africa. Interestingly, 67% of the studies from Africa are also within the East and Southern Africa region which accounts for over 54% of the global HIV cases (73). Our finding is consistent with earlier systematic reviews which showed similar global prevalence of PLWH in COVID-19 cases (19,74).
Our result also shows that PLWH may not be at higher risk of severe COVID-19 defined by admission to intensive care units or need for mechanical ventilation on the global level. Interestingly, this lack of association between HIV infection and COVID-19 severity hold true in Europe and the United States, but not in Africa. We found a 14% increase in risk of severity of COVID-19 for PLWH in Africa.
Moreover, 60% (3/5) of the studies analysed for the risk of severe COVID-19 in Africa were conducted in South Africa and all studies originated from sub-Saharan Africa: a region associated with high HIV infection rate and poorer antiretroviral treatment (ART) availability (73). Further, we found a 52% increase in the risk of death from COVID-19 in PLWH on the global level. However, only the North . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint American (United States, USA) continent showed a significant risk of mortality (2-fold) among the regions computed. The 2-fold increased risk found in the United States was attributable to the study by Shalev et al included in the meta-analysis. Interestingly, this significant risk ratio disappeared on the continental and the global risk of mortality on the global level fell to 30% when the Shalev et al study was removed (fig. S2a). However, the Egger's test showed no significant small study effect when the study was included ( fig. 4 vs fig. S2). It should be noted also that the average age of PLWH in the Shalev study is 60.7 years, 71% of whom have more than one COVID-19 outcome-defining comorbidities. Moreover, most studies (6/9, 67%) within the USA were conducted in Georgia and New York, both of which are on the top ten states with highest HIV infections and hardest hit by The observed risk of severe illness (Africa) and death (globally) from COVID-19 in these studies may be attributed to the interplay between several factors. Firstly, the availability of effective HIV management tools in the developed countries means that PLWH now live longer in these regions (83). Increased age is associated with senescence of the natural immune system which may combine with other immune-dampening features of chronic, untreated HIV infection to increase the risk of . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint severity and death from COVID-19. Also, PLWH especially those with undiagnosed or uncontrolled infections, low CD4 count, opportunistic infections and high viral load may present with severe COVID-19 and are at higher risk of death (84). Aside from CD4 and CD8 T-cells activation, effective and early immunoglobin G (IgG) generation results in effective SARS-CoV-2 clearance and improves clinical outcomes (16,85). However, uncontrolled HIV replication, may trigger increased CD8 T cells activation, inflammation, T cells exhaustion and dysfunction in B cells activities (86,87). The Moreover, the prevention (sensitisation and pre-exposure prophylaxis), diagnosis and management (ART) of HIV and other chronic diseases have been affected by the global shift in medical resources to contain the COVD-19 pandemic (16). This shift has been suggested to be a contributary factor to the susceptibility of affected groups to severe COVD-19 and death (91)(92)(93). Expectedly, the disruption to healthcare systems especially HIV clinics and the downstream effect have been relatively worse in developing countries possibly resulting in worse outcome for PLWH co-infected with COVID-19 (16,94). However, more data will be needed to establish the extent of these disruptions in regions already behind in the fight against HIV.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint Put together, our result show that while the risk of severity illness and death due to COVID-19 increases respectively in Africa and globally, the mechanistic link between HIV infection and the clinical course of COVID-19 may be more complex than previously thought. Firstly, the regional aggregation performed in this study showed that the prevalence of PLWH in COVID-19 cases can be best translated in the background of current global epidemiology of HIV infections. Indeed, the variation introduced by the differences in regional HIV infection rates makes estimation of global prevalence of HIV-SARS-COV-2 coinfection less reliable if not controlled for regional prevalence of This study has several limitations. Firstly, some of the included studies are case series reporting only the PLWH co-infected with SARS-CoV-2. Such studies were excluded from the prevalence analysis, and this may explain the high heterogeneity in the results. However, the random effect model was used to account for the variations in experimental design. Secondly, most studies did not report the distributions of comorbidities, race, age, CD4 and CD8 counts, duration of HIV infection, ART use among other confounders in the studied groups. Thus, we could not control for these parameters in this study. Also, most studies did not report the clinical outcomes (death and severity) of COVID-19 in both PLWH and patients who are HIV-negative and could not be included in the relative risk . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint computation. However, the Egger's test and funnel plots show that there is no publication bias in the analysed data.
Our findings have several clinical and research implications. First, it further widens the body of evidence by including more recent studies to report that PLWH may be at an increased risk of severe COVID-19 and death and which regions may be mostly at risk. Secondly, we show that the risk of . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 16, 2021. ; https://doi.org/10.1101/2021.07.12.21260361 doi: medRxiv preprint