ABSTRACT
Objectives To define the incidence of clinically-detected COVID-19 in people with HIV (PWH) in the US and evaluate how racial and ethnic disparities, comorbidities, and HIV-related factors contribute to risk of COVID-19.
Design Observational study within the CFAR Network of Integrated Clinical Systems cohort in 7 cities during 2020.
Methods We calculated cumulative incidence rates of COVID-19 diagnosis among PWH in routine care by key characteristics including race/ethnicity, current and lowest CD4 count, and geographic area. We evaluated risk factors for COVID-19 among PWH using relative risk regression models adjusted with disease risk scores.
Results Among 16,056 PWH in care, of whom 44.5% were Black, 12.5% were Hispanic, with a median age of 52 years (IQR 40-59), 18% had a current CD4 count < 350, including 7% < 200; 95.5% were on antiretroviral therapy, and 85.6% were virologically suppressed. Overall in 2020, 649 PWH were diagnosed with COVID-19 for a rate of 4.94 cases per 100 person-years. The cumulative incidence of COVID-19 was 2.4-fold and 1.7-fold higher in Hispanic and Black PWH respectively, than non-Hispanic White PWH. In adjusted analyses, factors associated with COVID-19 included female sex, Hispanic or Black identity, lowest historical CD4 count <350 (proxy for CD4 nadir), current low CD4/CD8 ratio, diabetes, and obesity.
Conclusions Our results suggest that the presence of structural racial inequities above and beyond medical comorbidities increased the risk of COVID-19 among PWH
PWH with immune exhaustion as evidenced by lowest historical CD4 or current low CD4:CD8 ratio had greater risk of COVID-19.
Introduction
The COVID-19 pandemic has had profound direct and structural effects on the health of people living with HIV (PWH) in the US[1-3]. Health disparities among Black and Hispanic populations in the US are closely related to structural inequities in social determinants of health[4]. The burden of the COVID-19 epidemic in Black and Hispanic communities is the result of deep socioeconomic disparities and inequities that affect health and access to healthcare, increased occupational and household exposure, and physiologic risks resulting from heightened allostatic loads and a disproportionate prevalence of medical comorbidities[4-10].
Black and Hispanic Americans are overrepresented among PWH due to shared structural inequities that also increase risk of HIV acquisition[11]. Such inequities may explain, in part, the observations from large population-based studies that showed PWH have a higher risk of severe disease and mortality due to COVID-19 than those without HIV[12-16]. Structural inequities in areas of residence, education, and wealth are more prevalent among PWH than in the general population[17], and unevenly distributed by race, ethnicity, and geographic region in the US[18-20].
Little is known about the effect of racial and ethnic disparities on risk of COVID-19 disease in PWH. Structural vulnerabilities among PWH likely contribute to increased risk of exposure to SARS-CoV-2 as well as differential access to testing and healthcare, making infections more or less likely to come to clinical attention. PWH also have a high prevalence of non-HIV-specific comorbidities associated with worse outcomes of SARS-CoV-2 infection. There are significant geographic, racial/ethnic, and income disparities in the distribution of these comorbidities among PWH compared to the general population[21]. It is essential to elucidate the role these factors play in COVID-19 incidence, separately from factors that may increase severity of COVID-19 disease in PWH[22]. We conducted this study in the longstanding Centers for AIDS Research Network of Integrated Clinical Systems (CNICS) multi-site US cohort of PWH.
Methods
CNICS is a prospective observational cohort study of adult PWH in routine clinical care at 10 academic institutions across the United States[23]. We studied all PWH engaged in care (defined as one or more in-person or virtual HIV primary care visits) between September 1, 2018, and December 31, 2020, and alive on March 1, 2020, at seven participating CNICS sites: Johns Hopkins University, Case Western Reserve University, Fenway Health, University of Alabama at Birmingham, University of California-San Diego, University of North Carolina at Chapel Hill, and University of Washington. CNICS research has been approved by the institutional review boards at each site. Sites were anonymized per regulatory requirements of some sites.
Methods of data collection for the CNICS cohort have been previously reported. Briefly, comprehensive clinical data collected through electronic medical records and other institutional data systems undergo rigorous data quality assessment and are harmonized in a central data repository that is updated quarterly[23]. Demographic data, including risk factors for HIV acquisition, are collected at cohort enrollment. Smoking status (current vs. former vs. never tobacco use) was determined through Patient-Reported Outcomes (PRO) collected through tablet-based surveys conducted every 4-6 months in conjunction with primary care visits[24] and supplemented by diagnosis codes.
SARS-CoV-2 infections
Candidate COVID-19 cases were identified through laboratory test results and provider documented diagnoses (ICD-10) recorded between March 1 and December 31, 2020, and verified through medical record review. SARS-CoV-2 test results and hospitalization records available to CNICS sites from external health systems were reviewed for data completeness. Medical records for all cases identified by diagnosis code without supporting laboratory tests performed within the CNICS system were reviewed by clinicians at each site using a standardized protocol. All appropriately-adjudicated cases were included, regardless of symptom status, which was not always documented for laboratory-diagnosed cases.
Covariates
We examined demographic covariates (age, sex at birth, self-reported race and Hispanic ethnicity, CNICS site) and the following chronic comorbid conditions: diabetes defined using a previously validated approach as hemoglobin A1c (HbA1c) ≥6.5%, a prescription of a diabetes-specific medication, or a diagnosis of diabetes with diabetes-related prescription[25]; treated hypertension as a diagnosis of hypertension with prescription of an anti-hypertensive medication; obesity as body mass index (BMI) ≥30 kg/m2; hepatitis C virus (HCV) coinfection by the presence of positive HCV antibody or detectable RNA or genotype; and chronic obstructive pulmonary disease (COPD) using a previously validated approach[26] as a diagnosis of COPD and ≥90-day continuous supply of long-acting controller medications. We used laboratory test data to compute clinical measures of chronic kidney disease (CKD) defined as last glomerular filtration rate (eGFR) <60 using CKD-EPI without race adjustment[27, 28], and risk scores for atherosclerotic cardiovascular disease (ASCVD[29]) and hepatic fibrosis (Fibrosis-4: FIB-4[30-32]). We examined CD4 counts (cells/μL) as both lowest historical value, as a proxy for CD4 nadir, and current CD4 count and CD4/CD8 ratio, as well as HIV viral load (VL; copies/mL). To avoid protopathic bias all laboratory test results were censored one week prior to COVID-19 diagnosis.
Statistical analysis
We calculated the cumulative incidence and incidence rate of first clinically-detected SARS-CoV-2 infection among PWH by key characteristics from March 1, 2020 through December 31, 2020. SARS-CoV-2 incidence was contrasted by demographic factors including race and ethnicity and geographic location, as well as clinical factors including lowest and current CD4 cell counts. Relative risks for COVID-19 were calculated using relative risk regression[33]. We accounted for potential confounding by adjusting models using disease risk scores (DRS)[34-36], the prognostic analogue of propensity scores useful when studying a limited number of exposed patients and outcomes, and a relatively large number of potential confounders (conceptual model, Supplemental Figure 1). DRS were constructed independently for each exposure of interest using logistic regression in the full cohort for all non-duplicative covariates[35] (e.g., FIB-4 not adjusted for age, a component of the FIB-4 calculation) including age, birth sex, race/ethnicity, smoking status (current vs. former vs. never use), diabetes, hypertension, and CNICS city. All analyses were conducted in Stata version 17 (StataCorp, College Station, TX).
Conceptual model of pathway to identification of a case of COVID-19 – simplified
* pathway to case identification (unmeasured components of what is ultimately measured)
† end of pathway: Case of COVID-19 (measured)
‡ modifiable COVID-related behaviors (unmeasured)
§ non-modifiable demographic characteristics
□ clinical/biological factors (measured)
¶ unmeasured/latent individual properties/experiences that substantially affect a person’s health and behaviors Shade of grey arrows is proportional to degree of effect
We employed disease risk scores (DRS) rather than standard multivariate adjustment to prevent overcorrection for colinear effects, especially when attempting to correct for the disproportionate burden of medical comorbidities on Black and Hispanic PWH to evaluate the residual effects of racism/marginalized identities on the association of racial/ethnic identity on COVID-19 diagnosis.
Abbreviations: NPI, non-pharmacologic interventions; BMI, body mass index; COPD, chronic obstructive pulmonary disorder; HTN, hypertension; CKD, chronic kidney disease
Results
Among 16,056 PWH in care during the study period, 44.5% identified as non-Hispanic Black, 37.9% as non-Hispanic White, 12.5% Hispanic, and 5.2% mixed, other, or unreported race/ethnicity. The median age was 52 years (IQR 40-59), 20.9% were female, and 53.6% were MSM. Most (95.5%) were on antiretroviral therapy (ART), and 85.6% had an undetectable HIV VL (Tables 1, 2). Overall in 2020, 649 PWH were diagnosed with COVID-19, for a cumulative incidence of clinically-detected infections of 4.0% (95% CI 3.7-4.3). of by the end of 2020. The majority (77%) of COVID-19 cases were confirmed by laboratory testing (SARS-CoV-2 RT-PCR) performed within the CNICS system. The remaining 23% of cases were verified from a variety of sources, including PCR and antigen testing performed in community or public health settings or at non-CNICS-affiliated medical sites. As previously reported, 108 cases (16.3%) were hospitalized, and 12 (1.9%) died[22].
COVID-19 incidence among PWH increased from 3.29/100 person-years (174 events, 5283 PY) during the second quarter of 2020 to 4.94/100 PYs (95% CI 4.57-5.33) by the end of 2020 By the end of 2020, the cumulative incidence in Black PWH was 5.64/100PY, 7.95/100PY in Hispanic PWH, and 3.30/100 PY for non-Hispanic White (p<0.0001 for both comparisons to non-Hispanic White (Table 3 and Figure 1)). Similarly, the IRR for COVID-19 was 1.71 (95% CI 1.42-2.07) for Black and 2.40 (95% CI 1.91-3.01) for Hispanic relative to non-Hispanic White PWH. There were geographic differences in the cumulative incidence of COVID-19 across sites, with cumulative incidence rates ranging from 3.03 to 8.32 per PY by site. COVID-19 cumulative incidence was highest among those with prior low CD4 cell count <350. We did not observe differences according to current CD4 cell count.
Sites were anonymized based on regulatory requirements of this cohort to not present individual site data. Y axes are cumulative incidences, presented as the proportion of the population under study. Results for current or lowest historical CD4 T cell count wre not different when stratified using 200 cells/mL as the cutoff.
In unadjusted analyses, we observed sociodemographic disparities in risk of COVID-19 also seen in the US general population, namely women and Black and Hispanic PWH were at increased risk of COVID-19 compared to men and non-Hispanic White PWH, respectively (p<0.001 for all). However, older PWH were not at higher risk of COVID-19 (Table 1). PWH who had diabetes or a BMI ≥30 had greater risk of COVID-19 than PWH without those conditions (p<0.001, Table 2). We observed a trend towards higher risk of COVID-19 in those with hypertension (p=0.10) and the “smoker’s paradox”[37], with ever-smokers having lower risk of COVID-19 than never smokers (50.2% of cases vs. 55.2% non-cases, p=0.013). CKD, ASCVD risk score, COPD, HCV, and elevated FIB-4 score were not associated with COVID-19 in unadjusted analyses. PWH with a lowest historical CD4 count <350 cell/mm3 had a higher risk of COVID-19 than those with a lowest CD4 ≥350 (26.7% and 44.8% of COVID-19 cases vs. 23.7% and 42.6% of non-COVID cases, respectively). The most recent CD4 count measured a median interval of 8.3 (IQR 3.2-19.0) months before COVID-19 diagnosis in cases, was not associated with COVID-19 (p=0.73), nor was ART status or detectable VL, although few in the cohort were not on ART (4.5%) or unsuppressed (15.4%).
In adjusted analyses, racial and ethnic identity was a strong predictor of COVID-19 in PWH (Figure 2). Hispanic (RR 2.03, 95% CI 1.62-2.55, p<0.001) and Black PWH (RR 1.33, 95% CI 1.10-1.61, p=0.003) had the highest risks of COVID-19. Other characteristics associated with COVID-19 after adjustment included female sex (RR 1.32, 95% CI 1.12-1.56, p=0.001), diabetes (RR 1.37, 95% CI 1.15-1.63, p<0.001), and having a BMI ≥30 (RR 1.58, 95% CI 1.36-1.83, p<0.001). The smoker’s paradox remained in adjusted analysis (RR 0.77, 95% CI 0.66-0.89, p=0.001). Neither CKD, HTN, HCV, COPD, nor ASCVD or FIB-4 risk scores predicted COVID-19 diagnosis. PWH whose lowest CD4 count was <350 had significantly higher risk of COVID-19 compared with those with lowest count ≥350 (RR 1.20, 95% CI 1.02-1.42), and we did not find an association with current CD4 count, at either a cut-point of <200 or <350. In addition, PWH with higher current CD4:CD8 ratio had decreased risk of COVID-19 (RR 0.92 for each 1 SD increase in ratio, 95% CI 0.86-0.99, p=0.042). Neither ART status nor detectable VL were predictive of COVID-19 in adjusted analysis (p=0.84 and 0.70, respectively).
Abbreviations: ASCVD, atherosclerotic cardiovascular disease; BMI, body mass index; CKD, chronic kidney disease; COPD, chronic obstructive pulmonary disease; eGFR, estimated glomerular filtration rate; FIB-4, Fibrosis-4 scoring system for liv r fibrosis; HCV, Hepatitis C virus; PWH, people with HIV; RR, relative risk VL, viral load.
a Relative risk regression models adjusted for demographic and clinical characteristics using disease risk scores, except for the ASCVD risk score analysis which is unadjusted. Disease risk scores were constructed independently for each exposure variable of interest using all non-duplicative covariates.
Discussion
In this large multi-site cohort of PWH receiving care in the US in 2020 prior to availability of SARS-CoV-2 vaccines, persons who were Hispanic, Black or female had significantly higher risk of COVID-19, as were those who had a lowest historical CD4 count <350 or lower current CD4/CD8 ratio. As seen in the general population and prior analyses of PWH, risk of COVID-19 diagnosis was higher in PWH with obesity and diabetes[12, 38-40]. The higher rates of COVID-19 among Black and Hispanic PWH persisted throughout 2020 after adjustment for other clinical and demographic risks. Results of our study in this diverse population of PWH suggest that structural disparities affect COVID-19 incidence above and beyond the disproportionate burden of medical comorbidities. Our findings in PWH, carefully adjusted for measurable demographic and clinical risk factors, expand on observations in the general population that have highlighted how structural inequities in housing, wealth, healthcare access, and occupations are associated with both increased risk of exposure to SARS-CoV-2 and of COVID-19 disease. The nearly three-fold risk difference of being diagnosed with COVID-19 between sites reflects differences in exposure risk and access to testing and care.
The finding that lowest CD4 count <350 cells/mL or low current CD4/CD8 ratio was more predictive of COVID-19 than current absolute CD4 count suggests that HIV-associated immune exhaustion is a driver of symptomatic SARS-CoV-2 infection[41]. These results align with other observations showing that SARS-Cov-2 specific T-cell responses correlate with CD4:CD8 ratio and percentage of naïve CD4-T Cells in PWH[42]. Neither ART status nor virologic suppression were associated with COVID-19, but few people in our cohort were not on ART or virologically unsuppressed. Notably, we found no association between age and COVID-19, which may reflect stricter adherence to COVID-19 precautions in older adults with HIV. Older PWH who have lived through the first decades of the AIDS pandemic, especially those with personal experience of AIDS, may have been more likely to limit social interactions or request work accommodations during 2020[43, 44]. However, once infected, older PWH were more likely to experience severe outcomes of COVID-19 infection[15, 22] as we previously reported [22].
Limitations of our study include lack of systematic surveillance testing, wherein case identification is limited by access to testing and availability of results, and ascertainment limited to PWH engaged in care. PWH with asymptomatic or minimally symptomatic COVID-19 were likely undertested and under-reported in medical records, a limitation also present in both cohort studies and public health data. Strengths of our study include comprehensive data on more than 16,000 PWH in care across the US with well-characterized comorbidities and HIV-specific clinical data. The latter allowed greater measurement precision to distinguish between the contribution of comorbidities, race/ethnic identity, and HIV factors to risk of COVID-19 diagnosis than studies using public health data. During the earliest months of the pandemic, people with mild symptoms or lacking risk factors for severe disease were discouraged from accessing clinical care, which may have introduced some bias in who was tested during those months[45]. Facility-based testing was often performed in emergency departments and hospitals, thus, people with advanced age or comorbidities may have been overrepresented amongst reported cases, although we do not expect that there was differential testing by CD4 count, ART status, sex, or race/ethnicity among people who presented to care. Conversely, persons with a history of AIDS or other comorbidities, who perceived themselves at greater risk of poor outcomes, may have been more likely to test after exposure or with mild symptoms, introducing potential confounding by indication.[45, 46] If present, these biases likely underestimate the effect of racial/ethnic identity on COVID-19 risk, as marginalized populations had less access to testing early in the pandemic, as well as differential health-seeking behaviors due to medical mistrust. While PWH who reported a history of IDU as a risk factor for HIV were not at increased risk of COVID-19, this historical risk category is a poor measure of current high-risk behaviors. We plan future analyses with a subsequent year of COVID-19 cases and PRO data in our cohort to study the effect of substance use including current IDU, and geocoding-based social determinants of health, on COVID-19 risk.
Race and ethnic identity are imperfect, albeit accessible proxies for the experience of structural racism, not directly assessable in this study. Other forms of stigmatization, including anti-LGBTQ bias, also increase allostatic loading[4, 47-50] and compound the structural barriers and health inequities faced by many PWH with marginalized identities. The increased incidence of COVID-19 among PWH with racialized identities and living in certain geographic regions further highlights how disparities in the US collide in the HIV and COVID-19 pandemics[3].
Conclusions
The results of our study in this diverse US population of PWH suggest that structural disparities differentially impact COVID-19 incidence, with higher and disproportionate rates of COVID-19 in women and Black and Hispanic PWH, and vast risk differences by geography. PWH with lowest CD4 count <350 cells/mL or lower current CD4/CD8 ratio had greater risk of COVID-19. Incidence estimates reflect clinically-detected COVID-19 diagnoses and provide conservative estimates for SARS-CoV-2 infection in the absence of surveillance testing.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Funding and Acknowledgments
We gratefully acknowledge the CNICS cohort participants and their providers for contributing to the study. Funding for this study and for CNICS came from the National Institute of Allergy and Infectious Diseases (NIAID) [CNICS R24 AI067039; UW CFAR NIAID Grant P30 AI027757; UAB CFAR grant P30 AI027767; UNC CFAR grant P30 AI50410; UCSD CFAR grant P30 AI036214 and by the Office of the President of the University of California COVID-19 seed grant; Case Western Reserve University CFAR grant P30 AI036219; Fenway Health/Harvard CFAR grant P30 AI060354, UCSF CFAR grant P30AI027763 and JHU CFAR grant P30 AI094189] and the National Institute on Drug Abuse (NIDA) [R01DA047045].
Footnotes
COI: RBI receives consulting fees from SeaGen and Abbvie, unrelated to this topic; EC holds Investigator-initiated grant support and medical advisory board membership, Gilead Sciences; HC receives grant support from ViiV, unrelated to this topic