HIV and risk of COVID-19 death: a population cohort study from the Western Cape Province, South Africa
======================================================================================================

* Mary-Ann Davies

## Abstract

**Background** The effect of HIV co-infection on COVID-19 outcomes in sub-Saharan Africa is unknown.

**Methods** We conducted a population cohort study using linked data from adults attending public sector health facilities in the Western Cape, South Africa. We used Cox-proportional hazards models adjusted for age, sex, location and comorbidities to examine the association between HIV and COVID-19 death among (i) public sector “active patients” (≥1 health visit in the 3 years before March 2020), (ii) laboratory-diagnosed COVID-19 cases and (iii) hospitalized COVID-19 cases. COVID-19 was diagnosed with SARS-CoV-2 PCR tests. We calculated the standardized mortality ratio (SMR) for COVID-19 comparing HIV positive vs. negative adults using modelled population estimates.

**Results** Among 3,460,932 public sector patients (16% HIV positive), 22,308 were diagnosed with COVID-19, of whom 625 died. In adjusted analysis, HIV increased risk of COVID-19 mortality (adjusted hazard ratio [aHR]:2.14; 95% confidence interval [CI]:1.70; 2.70), with similar risks across strata of viral load and immunosuppression. increased HIV-associated risk of COVID-19 death remained when restricting to COVID-19 cases (aHR:1.70; 95%CI:132; 2.18) or hospitalized cases (aHR:1.45; 95%CI:1.14; 1.84). Current and previous tuberculosis also increased COVID-19 mortality risk (aHR [95%CI]:2.70 [1.81; 4.04] and 1.51 [1.18; 1.93] respectively). The SMR for COVID-19 death associated with HIV was 2.39 (95% CI:1.96; 2.86); population attributable fraction 8.5% (95%CI:6.1; 11.1).

**Conclusion** HIV was associated with a doubling of COVID-19 mortality risk. While our findings may over-estimate the HIV-associated risk COVID-19 death due to residual confounding, people with HIV should be considered a high-risk group for COVID-19 management.

## Introduction

The effects of the intersecting pandemics of HIV, tuberculosis and coronavirus disease-19 (COVID-19) in sub-Saharan Africa are unknown. Studies to date suggest no increased risk of adverse outcomes in HIV co-infected patients, but these are small case series from Europe and North America, often limited to hospitalized patients, and may not be relevant to sub-Saharan Africa where people living with HIV (PLWH) are younger with different comorbidities, frequently including tuberculosis.1-7 PLWH may experience more severe COVID-19 disease due to HIV-related immune suppression, which may be exacerbated by transient immune deficiency from coronaviruses.8,9 In support of this hypothesis a large UK cohort study reported increased risk of COVID-19 death with immunosuppressive comorbidity, including PLWH.10 Two factors may reduce risk of severe COVID-19 in PLWH: dysfunctional cellular immunity could protect against the cytokine storm,11,12 and some antiretroviral drugs (tenofovir and lopinavir-ritonavir) have *in vitro* activity against coronaviruses, with better outcomes reported for PLWH receiving tenofovir disoproxil fumarate (TDF) vs. other antiretrovirals.12,13

It is important to establish if HIV and tuberculosis increase risk of COVID-19 death so that patients with these conditions can be provided with augmented prevention and potential therapeutic interventions. We used linked data from adults attending public sector health facilities in the Western Cape Province, South Africa, to identify factors associated with COVID-19 death, including HIV and tuberculosis.

## Methods

### Study design

We conducted a cohort study using de-identified data from the Western Cape Provincial Health Data Centre (WCPHDC) of public sector patients aged ≥20 years with documented sex and not known to have died before March 1, 2020 (before the first diagnosed COVID-19 case in South Africa and several weeks before the first COVID-19 death) and included all follow up through June 9, 2020. The outcome was COVID-19 death. Our main analysis examined risk of COVID-19 death in the general population, so all patients were included regardless of SARS-CoV-2 testing. The study was approved by the University of Cape Town and Stellenbosch University Health Research Ethics Committees and the Western Cape Province Department of Health. Individual informed consent requirement was waived for this secondary analysis of de-identified data.

### Study population and data sources

The Western Cape has nearly 7 million inhabitants, of whom ∼520,000 are PLWH with more than 90% of them dependent on public sector health services. The WCPHDC has been described in detail.14 Briefly, WCPHDC consolidates administrative, laboratory, and pharmacy data from routine electronic clinical information systems used in all public sector health facilities with linkage through a unique identifier. Multiple data sources are triangulated to enumerate health conditions such as diabetes mellitus (“diabetes”), hypertension, tuberculosis and HIV, with high or moderate certainty evidence assigned for each inferred condition (Supplementary Table 1). Patients are deemed to be PLWH (high certainty) if they have a positive HIV diagnostic test and/or HIV-RNA test and/or have received triple antiretroviral therapy (ART) and/or are registered in the HIV chronic disease management system; moderate certainty is assigned for those with only a CD4 count measure and/or two antiretroviral drugs prescribed and/or ICD-10 diagnosis code of HIV. HIV testing coverage is high; >90% of PLWH in the province know their HIV diagnosis.15 Comorbidities were based on high or moderate certainty evidence for the main analysis and restricted to high certainty evidence in sensitivity analyses. The virologic, immunologic and ART status of PLWH on March 1, 2020 was categorized, based on most recent measures, as “confirmed virologically suppressed on ART” (HIV-RNA<1000 copies/ml in last 15 months and ART dispensed in last 6 months), “likely virologically suppressed on ART” (HIV-RNA <1000 copies/ml 15-24 months previously or HIV-RNA <1000 copies/ml >24 months previously if ART dispensed in last 6 months), “viraemic or immunosuppressed” (HIV-RNA >1000 copies/ml in last 15 months or CD4 count <200 cells/µl within 18 months before March 2020) or “unknown”. Among hospitalized COVID-19 cases we examined the association between COVID-19 death and CD4 count at COVID-19 diagnosis/admission. Until January 2020, adult first-line ART was TDF plus emtricitabine/lamivudine plus efavirenz, with abacavir replacing TDF for patients with kidney disease; zidovudine plus emtricitabine/lamivudine plus a protease inhibitor was used for second-line ART. Dolutegravir was introduced in first and second-line since January 2020. Diabetic control was categorized according to glycosylated haemoglobin (HbA1c) measurement within the last 2 years as <7% (controlled); 7-8.9% (poorly controlled), ≥9% (uncontrolled).

### COVID-19 diagnosis

All COVID-19 diagnoses were based on positive SARS-CoV-2 PCR tests. Testing was available for all patients with COVID-19 symptoms until June 1, 2020; thereafter public sector laboratory testing was restricted to patients requiring admission or older than 55 years or with comorbidities, due to a temporary limit in testing capacity. Asymptomatic contacts of cases were not routinely tested. Hospital admissions and all deaths in SARS-CoV-2 positive cases are recorded and reviewed daily in the WCPHDC.

### Statistical analysis

We used Cox-proportional hazards models adjusted for age, sex and other comorbidities to examine the association between HIV and COVID-19 death among (i) all public sector patients with ≥1 health visit in the 3 years before March 1, 2020 (considered “active patients”), (ii) laboratory-diagnosed COVID-19 cases and (iii) hospitalized COVID-19 cases. We adjusted for location within Cape Town vs. rest of the province and subdistrict of residence for Cape Town patients to account for geographical variation in infection rates and as a proxy for socio-economic status. Patients were censored on the date of death if deceased without a COVID-19 diagnosis, or on June 9, 2020, whichever was earliest. Database closure was 7 days later to allow death reporting delays. For the analysis of COVID-19 death in laboratory-diagnosed cases we included cases diagnosed before June 1, 2020 when testing was available for all patients with COVID-19 symptoms, but included all patients diagnosed by June 9, 2020 in a sensitivity analysis. The proportional-hazard assumption was assessed with Schoenfeld residuals. All analyses were conducted using Stata 15.1.

We also calculated the standardized mortality ratio (SMR) of the actual number of COVID-19 deaths in people with HIV vs. the number that would be expected if PLWH had the same risk of COVID-19 death as HIV-negative people of the same age and sex. We used data on the age, sex and HIV status of all COVID-19 deaths (public and private sector) and the Thembisa Western Cape HIV model to estimate the Western Cape population size and HIV prevalence, by age and sex, in 2020.16 We used 1000 bootstrap replications to calculate 95% confidence intervals (CI) for the SMR (See Supplementary Appendix 2).

Since individual socio-economic status and some comorbidities are not recorded in WCPHDC, including body mass index (BMI), we calculated E-values to determine the minimum strength of association that an unmeasured confounder (such as raised BMI or differences in socio-economic status) would need to have with both HIV positive status and COVID-19 death to fully account for any association between HIV and COVID-19 death.17 We conducted quantitative bias analysis to assess the impact of potential confounding by obesity.

## Results

### Patient characteristics

Among 3,460,932 “active patients” aged ≥20 years or older on March 1, 2020, 22,308 were diagnosed with COVID-19, of whom 625 (2.8%) died (Table 1). The proportion of men was lower among COVID-19 cases vs. non-cases (31% vs. 42%), likely due to initial cases being among essential workers in sectors employing predominantly women. Diabetes and hypertension were common in all patients, with higher prevalence among COVID-19 cases than non-cases, (diabetes:14% vs. 8%; hypertension: 23% vs 16%), while proportions of other comorbidities were similar. COVID-19 deceased cases were older than surviving cases (median age [interquartile range] 63 years [54-71] vs. 37 [30-48]); several comorbidities were more common among deceased patients.

View this table:
[Table 1:](http://medrxiv.org/content/early/2020/07/03/2020.07.02.20145185.1/T1)

Table 1: 
Characteristics of (i) Western Cape “active patients” aged ≥20 years in public sector (public sector health care visit in last 3 years before March 1, 2020) according to COVID-19 outcome (ii) COVID-19 cases in “active patients” and (iii) hospitalized COVID-19 cases in “active patients”.

### Patients with and without HIV

Although the proportion of PLWH was similar among surviving and deceased COVID-19 cases, a greater proportion of COVID-19 deaths were in patients aged <50 years in those with vs. without HIV (39% vs 13%) (Table 2). A substantial proportion of COVID-19 deceased PLWH had diabetes (50%) and hypertension (42%), however these conditions were more common in deceased people without HIV. Current and previous tuberculosis were more frequent in PLWH irrespective of COVID-19; a high proportion of COVID-19 deceased cases with HIV had current (14%) and/or previous tuberculosis (37%).

View this table:
[Table 2:](http://medrxiv.org/content/early/2020/07/03/2020.07.02.20145185.1/T2)

Table 2: 
Characteristics of (i) Western Cape “active patients” ≥20 years of age in public sector (public sector health care visit in last 3 years) with and without HIV according to COVID-19 outcome.

### COVID-19 death in all public sector patients

Among all public sector patients, the probability of COVID-19 death by 100 days since March 1, 2020 was 180/million (95% confidence interval [CI]: 167; 196). COVID-19 death was associated with male sex, increasing age, diabetes (with higher risk with elevated HbA1c), hypertension and chronic kidney disease (Table 3). Current tuberculosis increased the hazard of COVID-19 death (adjusted hazard ratio [aHR]:2.70; 95%CI:1.81; 4.04) with a smaller increased risk for previous tuberculosis (aHR:1.51; 95%CI:1.18; 1.93). The increased risk of COVID-19 death with current tuberculosis was present for both rifampicin-sensitive and resistant disease during intensive phase treatment (Supplementary Table 3).

View this table:
[Table 3:](http://medrxiv.org/content/early/2020/07/03/2020.07.02.20145185.1/T3)

Table 3: 
Univariate and multivariate hazard ratios (HRs) and 95% confidence intervals (CI) for associations with COVID-19 death from March 1 to June 9, 2020, among all public sector patients ≥20 years with a public sector health visit in the previous 3 years (n=3,460,932), using Cox-proportional hazards models

The hazard of COVID-19 death adjusted for age and sex was increased in PLWH compared to HIV-negative patients, and the association increased further when adjusting for all recorded comorbidities (aHR:2.14; 95%CI:1.70; 2.70). Associations with most comorbidities increased when restricting to those with high certainty comorbidity evidence (Supplementary Table 4). Increased mortality risk was similar in PLWH irrespective of viraemia/immunosuppression, but there were few viraemic or immunosuppressed patients, as reflected in the wide CIs for the hazard ratio in different groups.

### Death in COVID-19 cases and hospitalized patients

Among the 15,203 COVID-19 cases diagnosed before June 1, 2020, mortality was increased in men and older patients (Table 4; Figure 1). Although the increased risk of death associated with all comorbidities was attenuated compared to the population analysis, PLWH still had an increased hazard of death compared to HIV negative COVID-19 cases (aHR:1.70; 95%CI:1.32; 2.18), with similar results when including patients diagnosed after the change in testing criteria (Supplementary Table 5). Among COVID-19 cases in PLWH on ART, mortality was lower in patients on TDF (vs. abacavir/zidovudine) (aHR: 0.42; 95%CI: 0.22;0.78) with no difference for other antiretrovirals. Among hospitalized COVID-19 cases, the effects of all mortality risk factors including HIV were attenuated, but HIV remained associated with death (aHR:1.45; 95%CI:1.14; 1.84). Among hospitalized patients, CD4 <200 cells/µl at COVID-19 diagnosis or admission was associated with COVID-19 death (aHR vs HIV negative: 2.36; 95%CI:1.47; 3.78).

View this table:
[Table 4:](http://medrxiv.org/content/early/2020/07/03/2020.07.02.20145185.1/T4)

Table 4: 
Multivariate hazard ratios (HRs) and 95% confidence intervals (CI) for associations with COVID-19 death from Cox-proportional hazards models among (i) all adult COVID-19 cases diagnosed before June 1, 2020 (n=15,203) and (ii) all hospitalized adult COVID-19 cases (n=2,978).

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/03/2020.07.02.20145185.1/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2020/07/03/2020.07.02.20145185.1/F1)

Figure 1: 
Comparison of adjusted hazard ratios (HR) and 95% confidence intervals (CI) for associations with COVID-19 death from Cox-proportional hazards models among (i) all public sector patients ≥20 years with a public sector health visit in the previous 3 years (n=3,460,932) (ii) all adult COVID-19 cases diagnosed before June 1, 2020 (n=15,203) and (iii) all hospitalized COVID-19 cases (n=2, 978).

### Potential bias from unmeasured confounding

To assess whether the association between HIV and COVID-19 mortality could be due to residual unmeasured confounding e.g. by socio-economic status, or unrecorded comorbidities, we calculated the E-value for an unmeasured confounder. The E-value for the analysis among all public sector patients was 3.70 (and 2.79 for the lower bound of the CI), suggesting that only a strong association between both HIV and low socio-economic status (or another confounder such as raised BMI), and low socio-economic status and COVID-19 death would account for all of the observed association between HIV and COVID-19 death. Most public sector patients have low socio-economic status with limited variability, and the effect of HIV on COVID-19 death was similar when restricting to the poorest subdistrict with the highest HIV prevalence in Cape Town.15 Similarly, quantitative bias analysis showed that the HIV-associated increased risk of COVID-19 death was unlikely due to confounding by raised BMI (Supplementary Appendix 4).

### Standardized mortality ratio

Among all laboratory diagnosed COVID-19 cases, there were 135 deaths among an estimated ∼520,000 PLWH in the province (260 deaths/million) and 786 deaths among 6.36 million people without HIV (124 deaths/million). The SMR for COVID-19 mortality in PLWH, relative to HIV-negative people was 2.39 (95% CI: 1.96; 2.86) and the attributable fraction of COVID-19 deaths due to HIV was 8.5% (95% CI: 6.1; 11.1).

## Discussion

In this analysis of nearly 3.5 million adults (16% PLWH) in South Africa we found an approximately two-fold increased risk of COVID-19 death in PLWH, irrespective of viral suppression, and a similar increased risk for patients with current tuberculosis. PLWH receiving TDF had a lower risk of COVID-19 death compared to those on other antiretrovirals, however these results may be confounded in the context of a public health approach with TDF in recommended first-line ART unless contra-indicated. While the HIV- and tuberculosis-associated increased risk of COVID-19 death may be over-estimated if there is residual confounding due to socio-economic status or unrecorded comorbidities, our results support considering people with HIV and tuberculosis as being at elevated risk of severe COVID-19. Nonetheless, despite a high burden of advanced HIV in the province, the attributable fraction of all deaths which could be ascribed to HIV was less than 10%.

While most case series of HIV and SARS-CoV-2 co-infection have not shown high risk of poor outcomes in PLWH, 1-3,5-7 cohorts of hospitalized PLWH with COVID-19 in London and New York have reported substantial morbidity and mortality including among patients with a suppressed viral load on ART.18,19 A comparison of hospitalized COVID-19 PLWH and people without HIV in New York suggested that COVID-19 outcomes may be worse in PLWH, but the study was under-powered to detect differences by HIV status.20 Absence of increased mortality risk in hospitalized patients with comorbidities may be explained by selection bias; risk factors for COVID-19 death may be attenuated by restricting to the subset of hospitalized patients who are already at high risk of mortality.21 It is therefore expected that for all comorbidities, the increased risk of death was progressively attenuated when restricting to cases (people with sufficiently severe symptoms to be tested) and hospitalized patients.

Similar to our findings, several studies have reported a high prevalence of comorbidities among PLWH with severe COVID-19.3,6,7 In our study, the overall high prevalence of diabetes in people with and without HIV, high proportion with poor diabetic control and elevated risks for COVID-19 compared to those reported from other countries are particularly concerning.10 The high prevalence of comorbidities in deceased PLWH suggests that the effect of HIV may at least partly be due to an increased risk of comorbidities at younger ages, possibly due to the “aging effect” of HIV.2,7 including those not recorded in WCPHDC such as cardiovascular disease. Persistent immune dysfunction may also play a role in severe COVID-19 despite viral suppression, and the hazard ratio point estimates for association with COVID-19 death were greater in immunosuppressed or viraemic PLWH, although the numbers of these patients with COVID-19 were small with wide CIs. Our finding that among COVID-19 cases in PLWH on ART those on TDF-based ART (vs. other therapies) were protected against COVID-death concurs with a recently published cohort of PLWH on ART from Spain.13 The protective effect of TDF in our setting may be over-estimated as TDF is only used in first-line therapy and contra-indicated in patients with poor renal function, however our analysis was adjusted for kidney disease and viral suppression. CD4 <200 cells/µl during admission was associated with COVID-19 death; this is likely due to the well-described lymphopenia in severe COVID-19 which is known to be prognostic of poor outcomes.9

We found both current and previous tuberculosis to be associated with COVID-19 death. While tuberculosis may exacerbate COVID-19 with impaired immune responses and increased angiotensin converting enzyme 2 receptor expression in respiratory epithelial cells, COVID-19 pneumonia may enhance tuberculosis progression.22-24 Nonetheless, since current tuberculosis itself causes death, it is difficult to disentangle the effects of COVID-19 vs tuberculosis disease on mortality *per se*.

To our knowledge this is the largest report on SARS-CoV-2 and HIV co-infected patients. Strengths include the study size using population-level data, laboratory confirmed SARS-CoV-2 diagnoses in all COVID-19 patients, inclusion of hospitalized and non-hospitalized cases and deaths, and coherence of associations from population (cohort study and SMR) to diagnosed to hospitalised patients. Being an observational study, limitations include potential under-ascertainment of comorbidities in routine administrative data, lack of data on other potential risk factors including socio-economic status and BMI, and relatively large numbers of PLWH without recent viral load or CD4 count results. In particular, patients with no recent information on disease control (e.g. HIV-RNA or HBA1c) may have less contact with health services and not reside permanently in the province, with under-ascertainment of outcomes.

## Conclusion

While our findings of a doubling of mortality risk in COVID-19 PLWH may over-estimate the effect of HIV on COVID-19 death due to the presence of residual confounding, PLWH should nonetheless be considered for inclusion as a high-risk group for COVID-19 management irrespective of virologic suppression, especially if they have other comorbidities.

## Data Availability

Data from the study is not available as was conducted on routine patient data and the requirement for informed consent was waived. Analysis code is available on request.

## Acknowledgements

We would like to acknowledge all patients in the Western Cape and to thank the Western Cape Department of Health Provincial Health Data Centre, the Western Cape Department of Health COVID-19 Outbreak Response Team, the Western Cape Communicable Disease Control sub-directorate and Western Cape health care workers involved in the COVID-19 response for their contributions to this report. We acknowledge funding for the Western Cape Provincial Health Data Centre from the Western Cape Department of Health, the US National Institutes for Health (R01 HD080465, U01 AI069924), the Bill and Melinda Gates Foundation (1164272), and the Wellcome Trust (203135/Z/16/Z).

## Footnotes

*   Author: Western Cape Department of Health in collaboration with the National Institute for Communicable Diseases, South Africa.

*   The list of all contributing authors is in the supplementary appendix.

*   Received July 2, 2020.
*   Revision received July 3, 2020.
*   Accepted July 3, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Blanco JL, Ambrosioni J, Garcia F, et al. COVID-19 in patients with HIV: clinical case series. Lancet HIV 2020;7:e314–e6.
    
    

2.  2.Gervasoni C, Meraviglia P, Riva A, et al. Clinical features and outcomes of HIV patients with coronavirus disease 2019. Clinical Infectious Diseases 2020. doi: 10.1093/cid/ciaa579.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa579&link_type=DOI) 

3.  3.Härter G, Spinner CD, Roider J, et al. COVID-19 in people living with human immunodeficiency virus: a case series of 33 patients. Infection 2020:1–6.
    
    

4.  4.Jones R, Nelson M, Bracchi M, Asboe D, Boffito M. COVID-19 in patients with HIV. Lancet HIV 2020;7:e383.
    
    

5.  5.Ridgway JP, Farley B, Benoit JL, et al. A Case Series of Five People Living with HIV Hospitalized with COVID-19 in Chicago, Illinois. AIDS Patient Care and STDs 2020. doi:10.1089/apc.2020.0103.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/apc.2020.0103&link_type=DOI) 

6.  6.Shalev N, Scherer M, LaSota ED, et al. Clinical characteristics and outcomes in people living with HIV hospitalized for COVID-19. Clinical Infectious Diseases 2020. doi:10.1093/cid/ciaa635.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa635&link_type=DOI) 

7.  7.Vizcarra P, Pérez-Elías MJ, Quereda C, et al. Description of COVID-19 in HIV-infected individuals: a single-centre, prospective cohort. Lancet HIV 2020. S2352-3018(20)30164-8.
    
    

8.  8.Qin C, Zhou L, Hu Z, et al. Dysregulation of immune response in patients with COVID-19 in Wuhan, China. Lancet Infectious Disease 2020. (Accessed 1 June 2020; at: [http://dx.doi.org/10.2139/ssrn.3541136](http://dx.doi.org/10.2139/ssrn.3541136))
    
    

9.  9.Tan L, Wang Q, Zhang D, et al. Lymphopenia predicts disease severity of COVID-19: a descriptive and predictive study. Signal Transduction and Targeted Therapy 2020;5.
    
    

10. 10.Williamson E, Walker AJ, Bhaskaran K, et al. OpenSAFELY: factors associated with COVID-19-related hospital death in the linked electronic health records of 17 million adult NHS patients. 2020. (Accessed 8 May 2020, at [https://www.medrxiv.org/content/10.1101/2020.05.06.20092999v1](https://www.medrxiv.org/content/10.1101/2020.05.06.20092999v1).)
    
    

11. 11.Mascolo S, Romanelli A, Carleo MA, Esposito V. Could HIV infection alter the clinical course of SARS-CoV-2 infection? When less is better. Journal of Medical Virology 2020. doi: 10.1002/jmv.25881.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.25881&link_type=DOI) 

12. 12.Laurence J. Why Aren’t People Living with HIV at Higher Risk for Developing Severe Coronavirus Disease 2019 (COVID-19)? AIDS patient care and STDs 2020;34:247–8.
    
    

13. 13.Del Amo J, Polo R, Moreno S, et al. Incidence and Severity of COVID-19 in HIV-Positive Persons Receiving Antiretroviral Therapy: A Cohort Study. Annals of Internal Medicine 2020. doi: 10.7326/m20-3689.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/m20-3689&link_type=DOI) 

14. 14.Boulle A, Heekes A, Tiffin N, et al. Data Centre Profile: The Provincial Health Data Centre of the Western Cape Province, South Africa. International Journal of Population Data Science 2019;4.
    
    

15. 15.Western Cape Burden of Disease Rapid Review Update 2019. 2020. (Accessed 1 March 2020, at [https://www.westerncape.gov.za/assets/departments/health/burden\_of\_disease\_report\_2020.pdf](https://www.westerncape.gov.za/assets/departments/health/burden\_of_disease_report_2020.pdf).)
    
    

16. 16.Thembisa AIDS and Demographic Model for South Africa v4.3. 2020. (Accessed 28 June 2020 at [www.thembisa.org](http://www.thembisa.org).)
    
    

17. 17.VanderWeele TJ, Ding P. Sensitivity Analysis in Observational Research: Introducing the E-Value. Annals of Internal Medicine 2017;167:268–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M16-2607&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28693043&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F03%2F2020.07.02.20145185.1.atom) 

18. 18.Childs K, Post FA, Norcross C, et al. Hospitalized patients with COVID-19 and HIV: a case series. Clinical Infectious Diseases 2020. doi: 10.1093/cid/ciaa657.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa657&link_type=DOI) 

19. 19.Suwanwongse K, Shabarek N. Clinical features and outcome of HIV/SARS-CoV-2 co-infected patients in the Bronx, New York City. Journal of Medical Virology 2020. doi:10.1002/jmv.26077.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.26077&link_type=DOI) 

20. 20.Karmen-Tuohy S, Carlucci PM, Zacharioudakis IM, et al. Outcomes among HIV-positive patients hospitalized with COVID-19. Journal of Acquired Immune Deficiency Syndromes 2020. doi: 10.1097/QAI.0000000000002423.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/QAI.0000000000002423&link_type=DOI) 

21. 21.Griffiith G, Morris T, Tudball M, et al. Collider bias undermines our understanding of COVID-19 disease risk and severity. 2020. (Accessed 6 June 2020, 2020, at [https://www.medrxiv.org/content/10.1101/2020.05.04.20090506v3.full.pdf](https://www.medrxiv.org/content/10.1101/2020.05.04.20090506v3.full.pdf).)
    
    

22. 22.Wilkinson RJ. Tuberculosis and type 2 Diabetes Mellitus: an inflammatory danger signal in the time of COVID-19. Clinical Infectious Diseases. 2020. doi: 10.1093/cid/ciaa747.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa747&link_type=DOI) 

23. 23.Motta I, Centis R, D’Ambrosio L, et al. Tuberculosis, COVID-19 and migrants: preliminary analysis of deaths occurring in 69 patients from two cohorts. Pulmonology 2020. doi: 10.1016/j.pulmoe.2020.05.002.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pulmoe.2020.05.002&link_type=DOI) 

24. 24.Tadolini M, Codecasa LR, García-García JM, et al. Active tuberculosis, sequelae and COVID-19 co-infection: first cohort of 49 cases. European Respiratory Journal 2020. doi: 10.1183/13993003.01398-2020
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZXJqIjtzOjU6InJlc2lkIjtzOjEyOiI1Ni8xLzIwMDEzOTgiO3M6NDoiYXRvbSI7czo1MjoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8wMy8yMDIwLjA3LjAyLjIwMTQ1MTg1LjEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)