00421244284Int J CancerInt. J. CancerInternational journal of cancer0020-71361097-021531215037689872610.1002/ijc.32260NIHMS1017836ArticleCervical cancer risk in women living with HIV across four continents: a multicohort studyRohnerEliane1BütikoferLukas2SchmidlinKurt1SengayiMazvita3MaskewMhairi4GiddyJanet5TaghaviKatayoun1MooreRichard D.6GoedertJames J.7GillM. John8SilverbergMichael J.9D’SouzaGypsyamber10PatelPragna11CastilhoJessica L.12RossJeremy13SohnAnnette13Bani-SadrFirouze14TaylorNinon15PaparizosVassilios16BonnetFabrice1718VerbonAnnelies19VehreschildJörg Janne2021PostFrank A.22SabinCaroline23MocroftAmanda23DrondaFernando24ObelNiels25GrabarSophie262728SpagnuoloVincenzo29Quiros-RoldanEugenia30MussiniCristina31MiroJosé M.32MeyerLaurence3334HasseBarbara35KonopnickiDeborah36RocaBernardino37BargerDiana18CliffordGary M.38FranceschiSilvia39EggerMatthias140BohliusJulia1for the AIDS-defining Cancer Project Working Group of IeDEA and COHERE in EuroCoordInstitute of Social and Preventive Medicine, University of Bern, SwitzerlandCTU Bern, University of Bern, SwitzerlandNational Cancer Registry, National Health Laboratory Service, Johannesburg, South AfricaHealth Economics and Epidemiology Research Office, Department of Internal Medicine, School of Clinical Medicine, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South AfricaDepartment of Medicine, McCord Hospital, Durban, South AfricaJohns Hopkins University, School of Medicine, Baltimore, MarylandDivision of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MarylandUniversity of Calgary, Alberta, CanadaDivision of Research, Kaiser Permanente Northern California, Oakland, USAJohns Hopkins Bloomberg School of Public Health, Baltimore, MarylandDivision of Global HIV and TB, Centers for Disease Control and Prevention, Atlanta, GeorgiaDivision of Infectious Diseases, Vanderbilt University Medical Center, Nashville, USATREAT Asia/amfAR - The Foundation for AIDS Research, Bangkok, ThailandReims Champagne-Ardenne University, Faculté de médecine, CHU Reims, Hôpital Robert Debré, Tropical and Infectious Diseases, Reims, FranceIIIrd Medical Department with Haematology, Medical Oncology, Haemostaseology, Infectious Diseases and Rheumathology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria, Present address: Department of Dermatology, University Hospital Salzburg, Paracelsus Medical University, Salzburg, AustriaAIDS Unit, Clinic of Venereologic and Dermatologic Diseases, Athens Medical School, “Syngros” Hospital, Athens, GreeceCHU de Bordeaux, Service de Médecine Interne et Maladies Infectieuses, Hôpital Saint-André, Bordeaux, FranceUniv. Bordeaux, ISPED, Centre INSERM U1219-Bordeaux Population Health, F-33000 Bordeaux, FranceDepartment Medical Microbiology and Infectious Diseases, Erasmus Medical Center, Rotterdam, The NetherlandsDepartment I of Internal Medicine, University Hospital of Cologne, Cologne, GermanyGerman Centre for Infection Research, partner site Bonn-Cologne, Cologne, GermanyKing’s College Hospital NHS Foundation Trust, London, UKInstitute for Global Health, UCL, London, United KingdomDepartment of Infectious Diseases, Hospital Ramón y Cajal, Madrid, SpainDepartment of Infectious Diseases, Copenhagen University Hospital, Copenhagen, DenmarkSorbonne Universités, UPMC Univ Paris 06, UMR_S 1136, Institut Pierre Louis d’Epidémiologie et de Santé Publique, F-75013, Paris, FranceINSERM, UMR_S 1136, Institut Pierre Louis d’Epidémiologie et de Santé Publique, F-75013, Paris, FranceUniversité Paris Descartes et Assistance Publique-Hôpitaux de Paris, Groupe hospitalier Cochin Hôtel-Dieu, Paris, FranceVita-Salute San Raffaele University, Faculty of Medicine and Surgery, Milan, ItalyInfectious and Tropical Diseases Institute, University of Brescia, Brescia, ItalyInfectious Diseases Clinics, University Hospital, Modena, ItalyInfectious Diseases Service, Hospital Clinic – IDIBAPS, University of Barcelona, Barcelona, SpainINSERM, U1018, Epidemiology of HIV, Reproduction, Paediatrics, CESP, University Paris-Sud, Paris, FranceDepartment of Public Health and Epidemiology, Bicêtre Hospital, AP-HP, Le Kremlin Bicêtre, Paris, FranceDivision of Infectious Diseases and Hospital Epidemiology, University Hospital Zurich, University of Zurich, SwitzerlandDepartment of Infectious Diseases, St Pierre University Hospital, Université Libre de Bruxelles, Brussels, BelgiumHospital General Universitario, Castellón, SpainInternational Agency for Research on Cancer, Lyon, FranceCentro di Riferimento Oncologico di Aviano (CRO)-IRCCS, Aviano, ItalyCentre for Infectious Disease Epidemiology and Research, University of Cape Town, Cape Town, South Africa.Corresponding author: Eliane Rohner, MD, Institute of Social and Preventive Medicine, University of Bern, Mittelstrasse 43, 3012 Bern, Switzerland; eliane.rohner@ispm.unibe.ch14320191962019012202001220211463601609

We compared invasive cervical cancer (ICC) incidence rates in Europe, South Africa, Latin and North America among women living with HIV who initiated antiretroviral therapy (ART) between 1996 and 2014. We analyzed cohort data from the International epidemiology Databases to Evaluate AIDS (IeDEA) and the Collaboration of Observational HIV Epidemiological Research in Europe (COHERE) in EuroCoord. We used flexible parametric survival models to determine regional ICC rates and risk factors for incident ICC. We included 64,231 women from 45 countries. During 320,141 person-years (pys), 356 incident ICC cases were diagnosed (Europe 164, South Africa 156, North America 19, Latin America 17). Raw ICC incidence rates per 100,000 pys were 447 in South Africa (95% confidence interval [CI] 382–523), 136 in Latin America (95% CI 85–219), 76 in North America (95% CI 48–119), and 66 in Europe (95% CI 57–77). Compared with European women ICC rates at 5 years after ART initiation were more than double in Latin America (adjusted hazard ratio [aHR] 2.43, 95% CI 1.27–4.68) and 11-times higher in South Africa (aHR 10.66, 95% CI 6.73–16.88), but similar in North America (aHR 0.79, 95% CI 0.37–1.71). Overall, ICC rates increased with age (>50 years versus 16–30 years, aHR 1.57, 95% CI 1.03–2.40) and lower CD4 cell counts at ART initiation (per 100 cell/µl decrease, aHR 1.25, 95% CI 1.15–1.36). Improving access to early ART initiation and effective cervical cancer screening in women living with HIV should be key parts of global efforts to reduce cancer-related health inequities.

Cervical cancerHIVincidence ratecohort studyIntroduction

Vast global inequities in the burden of invasive cervical cancer (ICC) exist.1,2 While access to effective screening and treatment of pre-cancerous cervical lesions has substantially reduced the risk of developing ICC in high-income countries, ICC remains a common cause of premature mortality and morbidity in women in low- and middle-income countries.1,2 ICC disproportionally affects women living with human immunodeficiency virus (HIV), who are more likely to have persistent co-infection with high-risk human papillomavirus (HPV) types,3 to develop pre-cancerous cervical lesions,4 and to progress to ICC than HIV-negative women.4 The advent and scale-up of combination antiretroviral therapy (ART) has led to a dramatic decline in morbidity and mortality from many HIV-associated diseases,5 but these decreases have not occurred for ICC.6,7 Indeed, as life expectancy after starting ART increases, there is more time for pre-cancerous cervical lesions to develop into ICC, but early initiation of ART seems to lower HPV co-infection rates and improve control of pre-cancerous cervical lesions.8

Global inequities in ICC incidence rates among women living with HIV have not been assessed previously. Our aim was to assess such inequities by comparing ICC incidence rates across different geographic regions among women who had initiated ART. Additionally, we examined risk factors for developing ICC in these women.

MethodsDatabases

We analyzed routinely collected clinical, demographic, laboratory, and treatment data of women enrolled in observational HIV cohorts that participate in the International epidemiology Databases to Evaluate AIDS (IeDEA) or the Collaboration of Observational HIV Epidemiological Research in Europe (COHERE) in EuroCoord. IeDEA has regional data centers in the Asia-Pacific, Australia, North America, Latin America and four African regions. Cohorts from the following IeDEA regions initially contributed data to this study: the North American AIDS Cohort Collaboration on Research and Design (NA-ACCORD),9 the Caribbean, Central and South America network for HIV epidemiology (CCASAnet),10 IeDEA Southern Africa,11 and IeDEA Asia-Pacific.12 For the IeDEA Southern Africa region, we restricted the analysis to two cohorts from South Africa that reduced under-reporting of cancer cases in the HIV cohorts through record linkages with the National Cancer Registry.13 COHERE is a collaboration of observational HIV cohorts across Europe.14 It contributed data from 24 cohorts, covering 36 countries. All cohorts obtained ethical approval from local ethics committees or institutional review boards, and the Cantonal Ethics Committee of Bern (number 028/2015) also granted ethical approval for this study.

Inclusion criteria and definitions

We restricted the analysis to cohorts that systematically collected cancer data or had enhanced their data through record linkages with cancer registries. We included women living with HIV who started ART after 1995 at 16 years or older. We excluded women who started ART before enrolment into cohort, women without follow-up after ART initiation, and women without any CD4 cell count measurements at ART initiation or during follow-up. We also excluded cohorts with less than 100 eligible women and the Asia-Pacific region because of small sample size (post-hoc decision). We analyzed ICC cases diagnosed any time after ART initiation as incident cases and excluded women diagnosed with ICC before or at ART initiation (prevalent ICC cases) from the analysis. We defined ART as a combination of at least three antiretroviral drugs from any class, including protease inhibitors (PIs), nucleoside reverse transcriptase inhibitors, and non-nucleoside reverse transcriptase inhibitors (NNRTIs). We assumed that women remained on ART and did not consider treatment interruptions and terminations. CD4 cell count at ART initiation was defined as the cell count closest to ART initiation, during the period within 180 days before to seven days after initiation.

Statistical analyses

We calculated raw ICC incidence rates by dividing the number of incident ICC cases by person-years (pys) at risk. Time at risk was measured from ART initiation to ICC diagnosis, last follow-up visit, death, or database closure, whichever happened first. We used proportional hazard flexible parametric survival models15 to estimate regional ICC incidence rates and to identify risk factors for developing ICC. We compared ICC rates at 2 years and 5 years after ART initiation across geographic regions. We used restricted cubic splines with 4 degrees of freedom and allowed for time-dependent region-effects with 2 degrees of freedom to model the baseline hazard. We performed likelihood ratio tests to test interactions between risk factors and regions. We assessed the following potential risk factors in the analysis: age at ART initiation (16–30, 31–50, >50 years); first-line ART regimen (NNRTI-based, PI-based, other); calendar period of ART initiation (1996–1998, 1999–2003, 2004–2007, 2008–2014); CD4 cell count at ART initiation; and current (time-updated) CD4 cell count. We treated CD4 cell count at ART initiation and current CD4 cell count as continuous variables. Analyses including CD4 cell count at ART initiation were restricted to women with available data on this variable. HIV RNA load at ART initiation was assessed in descriptive analyses.

We fit a crude model that included only the time-dependent region-effects, resulting in region-specific baseline hazards, and no other risk factors. The main adjusted model included region, CD4 cell count at ART initiation, age at ART initiation, first-line ART regimen, and calendar period of ART initiation. From the main adjusted model, we predicted ICC incidence rates for women with a specific set of risk factors, i.e. for women who initiated an NNRTI-based regimen between 2008–2014 at age 31–50 years with a CD4 cell count of 200 cells/µl. In a sensitivity analysis, the adjusted model included current (time-updated) CD4 cell count instead of CD4 cell count at ART initiation. In a second sensitivity analysis, we excluded ICC cases diagnosed within the first three months after ART initiation as prevalent cases and women with less than three months of follow-up. Results are presented as medians with interquartile ranges (IQR), number and percentages of women, incidence rates per 100,000 pys and hazard ratios (HRs) with 95% confidence intervals (CIs). We used Stata 14 (Stata Corporation, College Station, Texas, USA) and R (R Foundation, Vienna, Austria) for our analyses.

ResultsDescriptive analyses

The merged dataset included information on 126,063 women living with HIV. We excluded 44,419 women because they did not receive ART and another 14,413 women for reasons detailed in Supplementary Figures S1-S5. We made a post-hoc decision to exclude the Asia-Pacific region because too few eligible women remained after applying our exclusion criteria.

We included data on 64,231 women living with HIV, drawn from 36 cohorts and 45 countries across Europe, North America, Latin America, and South Africa (Figure 1). Overall, median age at ART initiation was 34.9 years (IQR 29.3–41.9), and was higher in North America (38.6 years) than in other regions (Table 1). Median CD4 cell count at ART initiation was 115 cells/µl (IQR 50–182) in South Africa, 178 cells/µl (IQR 74–281) in Latin America, and 241 cells/µl in both North America and Europe (Table 1). In South Africa, less than 1% of women started ART before 2004, but 26% of women in Latin America, 40% in Europe and 70% in North America initiated ART between 1996 and 2003. Most women in South Africa (93%) and Latin America (70%) received an NNRTI-based first-line regimen, but the majority of women in the European (55%) and North American (60%) cohorts received a PI-based first-line regimen. Median follow-up after ART initiation was around 5 years in Europe, North, and Latin America, but shorter in South Africa (2.1 years).

Over 320,141 pys of follow-up, 356 incident ICC cases were diagnosed (164 in Europe, 156 in South Africa, 19 in North America, and 17 in Latin America). In women who developed ICC, median time from ART initiation to ICC diagnosis was 1.9 years (IQR 0.7–4.2), and it ranged from 1.7 years in South Africa and North America to 2.6 years in Latin America (Supplementary Table S1). Median age at ICC diagnosis was 33 years in Latin America, and 38–40 years in South Africa, North America, and Europe. Median CD4 cell count at ICC diagnosis ranged from 275 cells/µl in Latin America to 370 cells/µl in North America.

Comparing ICC risk across regions

The raw ICC incidence rate was highest in South Africa, with 447/100,000 pys (95% CI 382–523), followed by Latin America (136/100,000 pys; 95% CI 85–219), North America (76/100,000 pys; 95% CI 48–119), and Europe (66/100,000 pys; 95% CI 57–77). In Europe, North America, and Latin America, there was some evidence for a decrease in crude and adjusted ICC incidence rates after more than one year on ART, except in South Africa (Figure 2). In crude analyses, ICC rates at 5 years after ART initiation were 11-times higher in women living with HIV in South Africa than in their European counterparts (HR 11.06, 95% CI 7.80–15.68). The much higher ICC rate in South African women was not explained by differences in CD4 cell count at ART initiation, age at ART initiation, first-line ART regimen, or calendar period of ART initiation (adjusted HR [aHR] 10.66, 95% CI 6.73–16.88; see Table 2). In crude (HR 2.32, 95% CI 1.24–4.31) and adjusted analyses (aHR 2.43, 95% CI 1.27–4.68), ICC rates at 5 years after ART initiation were more than twice as high in Latin American as in European women. In North American and European women, ICC rates after ART initiation were comparable in crude (HR 0.98, 95% CI 0.48–1.99) and adjusted analyses (aHR 0.79, 95% CI 0.37–1.71). The regional comparisons of ICC rates were similar at 2 years after ART initiation (Table 2). Also at 2 years after ART initiation, ICC rates were much higher in South Africa than in Europe (aHR 6.23, 95% CI 4.29–9.05). When we excluded ICC cases diagnosed within the first three months after ART initiation in a sensitivity analysis, results did not meaningfully change (Supplementary Table S2).

Risk factors for incident ICC

We did not find evidence of regional variation in the effect of CD4 cell count at ART initiation, age at ART initiation, first-line ART regimen, or calendar period of ART initiation on the risk of developing ICC (all p-values for interaction ≥ 0.13, see Table 3). Across all regions combined, the risk of developing ICC increased among women who initiated ART at lower CD4 cell counts (per 100 cell/µl decrease, aHR 1.25, 95% CI 1.15–1.36), and with older age at ART initiation (>50 years versus 16–30 years, aHR 1.57, 95% CI 1.03–2.40). There was no association between type of first-line ART regimen and the risk of developing incident ICC (PI-based versus NNRTI-based, aHR 1.05, 95% CI 0.79–1.41), and we did not observe a relevant decline in ICC rates by calendar period of ART initiation. The effects of the risk factors assessed in the main adjusted model remained similar when we excluded ICC cases diagnosed within the first three months after ART initiation from the analysis (Supplementary Table S3).

In a sensitivity analysis, we assessed the effect of current CD4 cell count on the risk of developing ICC and found that it varied across regions (p-value for interaction = 0.017). In analyses adjusted for age, first-line ART regimen, and calendar period of ART initiation, we did not find an association between current CD4 cell count and risk of developing ICC in South Africa (per 100 cells/µl decrease, aHR 1.00, 95% CI 0.92–1.10) or North America (aHR 1.08, 95% CI 0.90–1.30). However, a decrease of 100 cells/µl in current CD4 cell count increased the risk of developing ICC by 18% in European women (aHR 1.18, 95% CI 1.10–1.27) and 41% in Latin American women (aHR 1.41, 95% CI 1.07–1.86; see Supplementary Table S4 and Supplementary Figure S6).

Discussion

Across geographic regions, we found large inequities for cervical cancer incidence in women living with HIV. ICC incidence rates were high in women living with HIV in all regions studied, but the risk of developing ICC was much higher in women who had initiated ART in South Africa or Latin America than in women who had initiated ART in Europe or North America. Across all regions combined, the risk of developing ICC increased with older age and lower CD4 cell counts at ART initiation.

We believe this is the first study to provide a comparison of ICC incidence rates among women living with HIV across several geographic regions. To improve comparability of results across regions, we applied the same inclusion criteria and statistical methods across the whole dataset. With more than 60,000 women and 356 ICC cases included, this is also the largest study of ICC incidence in women living with HIV. However, several limitations of our study need to be acknowledged. Less than 20 ICC cases each were recorded in Latin America and North America. Thus, our comparison of ICC rates between those regions and Europe are of limited precision. ICC case identification and validation are likely to vary across regions and may have affected observed regional differences in ICC rates. Our results for South Africa may not be generalizable to Southern Africa as a region, given that we restricted our analyses to two urban cohorts in South Africa, which had been linked with the National Cancer Registry to reduce under-reporting of ICC cases. Because we included all women who started ART, irrespective of whether they remained in treatment, our results may not be representative of women who stayed continuously on ART. HIV RNA measurements at ART initiation were missing for one-third of women included in Latin America, and almost 80% of women from South Africa. Therefore, we could not use HIV RNA load to evaluate treatment response over time. Information on duration of HIV infection, HPV co-infection status, cervical cancer screening history, and smoking status was generally not available. Thus, we could not explore their effects on the risk of being diagnosed with ICC or adjust the regional comparisons for these potential confounders. Furthermore, as data on history of hysterectomy were not available, we could not exclude women who were no longer at risk of developing ICC. It would also have been interesting to assess ICC-related inequities in more depth, but we did not have data on ICC stage at diagnosis, for example.

We found that across all regions women living with HIV were at high risk of developing ICC after ART initiation. Most previous studies did not restrict their analyses to women who had initiated ART, but rather report ICC incidence estimates for women living with HIV irrespective of ART use.6,16–18 The raw ICC incidence rates in women living with HIV who had initiated ART, ranging from 66/100,000 pys in Europe to 447/100,000 pys in South Africa, were substantially higher than the ICC incidence rates reported for women from the general population in the included regions (≤30/100,000 pys).2 In our study, ICC rates after ART initiation were by far highest in South Africa, followed by Latin America, and they were lower in women who had started ART in North America or Europe. These findings corroborate the regional ICC incidence rate pattern in the general population,2 but the difference between South Africa and other regions is even more pronounced among women living with HIV.

The high ICC incidence rates we found in women from South Africa are similar to ICC incidence rates in women living with HIV in the United States in the early 1990s.7 In the United States, ICC incidence rates in women living with HIV had already dropped in the mid-1990s, before ART became available, and this drop has partly been attributed to better screening and more effective treatment of pre-cancerous cervical lesions.7 The extent to which ART protects women living with HIV from developing ICC is still being explored. Although ART reduces the prevalence of high-risk HPV in women living with HIV and promotes regression of cervical lesions,8 many women in our analyses, notably in South Africa, may have started ART too late, when potentially irreversible pre-cancerous cervical lesions were already present. Furthermore, not all women in our study would have achieved sustained suppression of HIV RNA, and high HIV RNA loads have been associated with an increased risk of HPV infection and cervical pre-cancerous lesions.19 Low CD4 cell counts have also been associated with a higher risk of HPV infection3 and development of severe cervical lesions.20 Accordingly, several studies showed an increased ICC risk in women with low nadir,21 baseline,18 or current CD4 cell counts.17 It remains a matter of debate at what stage of cervical carcinogenesis the effect of HIV-related immunodeficiency is largest. Across all regions combined, we found that the risk of developing ICC increased in women who initiated ART at low CD4 cell counts. High current CD4 cell counts had a protective effect in Latin America and Europe, but not in North America and South Africa.

Our analyses revealed massive regional differences in ICC rates in women living with HIV. Several factors could account for this finding. HPV prevalence in women living with HIV in sub-Saharan Africa or Latin America is higher than in North America or Europe,22 and this may contribute to the increased ICC burden in South African and Latin American women living with HIV. Women in South Africa and Latin America also tended to initiate ART at lower CD4 cell counts than women in Europe, and low CD4 cell counts at ART initiation increased the risk of developing ICC. Nevertheless, in our analyses large regional differences in ICC rates persisted after adjusting for CD4 cell counts. Therefore, inequities in access to effective cervical cancer screening and treatment of pre-cancerous cervical lesions are likely to be the main driver of regional variation in ICC rates in women living with HIV. Substantial global efforts are needed to improve cervical cancer screening and treatment for women living with HIV, and to promote national HPV vaccination programmes. Unfortunately, most Southern African countries and some regions of Latin America lack the resources to treat ICC.23,24

The availability of HPV vaccination and the long natural history from HPV infection through cervical intraepithelial neoplasia to invasive cancer make ICC particularly amenable to primary and secondary prevention.25 However, it has been estimated that in 2014 less than one-third of female adolescents aged 10–20 years in high-income countries and only 1% in low-income countries had received the full course of HPV vaccine.26 At present, data on HPV vaccination coverage among women living with HIV are lacking.27 Access to screening services with early detection and treatment of pre-cancerous cervical lesions remains key for ICC prevention in women living with HIV. However, there are extensive regional differences in access to effective cervical cancer screening. Less than 10% of women living in low-income countries have access to effective cervical cancer screening as compared to more than 60% in high-income countries.28 Integrating cervical cancer screening services into established HIV care programmes may facilitate screening access for women living with HIV and improve sustainability of screening programmes.29 Yet, it remains unclear how many women living with HIV actually receive regular screening for pre-cancerous cervical lesions. HIV cohorts and integrated cervical cancer screening services often do not systematically collect patient-level data on screening and treatment of pre-cancerous cervical lesions.30 Rigorous patient-level monitoring of cervical cancer screening and treatment programmes is essential to identify coverage gaps and target interventions.30

Conclusion

Our finding that women living with HIV who initiated ART in South Africa or Latin America were at much higher risk of developing ICC than women in North America or Europe reveals drastic global health inequities. ICC prevention through early ART initiation and scale up of effective cervical cancer screening services for women living with HIV, alongside the promotion of global access to HPV vaccination should be key parts of international efforts to reduce cancer-related health inequities.

Supplementary MaterialAcknowledgements

We thank all patients, care providers and data managers in the different IeDEA regions and COHERE in EuroCoord. We would also like to acknowledge Kali Tal for her editorial suggestions. More detailed Acknowledgement concerning the participating consortia can be found in the supplementary material.

Funding

Research reported in this publication was supported by the National Institute of Allergy and Infectious Diseases (NIAID), the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), the National Cancer Institute (NCI), the National Institute of Mental Health, and the National Institute on Drug Abuse of the U.S. National Institutes of Health (NIH) under Award Number U01AI069924 (Southern Africa), U01AI069907 (Asia-Pacific), U01AI069923 (Caribbean, Central, and South America), and U01-AI069918 (North America). National Institutes of Health (NIH) under Award Number U01AI069924 (Southern Africa), U01AI069907 (Asia-Pacific), U01AI069923 (Caribbean, Central, and South America), U01-AI069918 (North America), and U01A1096186 (the IeDEA Network Coordinating Center at Vanderbilt). The North American AIDS Cohort Collaboration on Research and Design (NA-ACCORD) was also supported by NIH grants F31DA037788, G12MD007583, K01AI093197, K23EY013707, K24AI065298, K24AI118591, K24DA000432, KL2TR000421, M01RR000052, N01CP01004, N02CP055504, N02CP91027, P30AI027757, P30AI027763, P30AI027767, P30AI036219, P30AI050410, P30AI094189, P30AI110527, P30MH62246, R01AA016893, R01CA165937, R01DA011602, R01DA012568, R01 AG053100, R24AI067039, U01AA013566, U01AA020790, U01AI031834, U01AI034989, U01AI034993, U01AI034994, U01AI035004, U01AI035039, U01AI035040, U01AI035041, U01AI035042, U01AI037613, U01AI037984, U01AI038855, U01AI038858, U01AI042590, U01AI068634, U01AI068636, U01AI069432, U01AI069434, U01AI103390, U01AI103397, U01AI103401, U01AI103408, U01DA03629, U01DA036935, U01HD032632, U10EY008057, U10EY008052, U10EY008067, U24AA020794, U54MD007587, UL1RR024131, UL1TR000004, UL1TR000083, UL1TR000454, UM1AI035043, Z01CP010214 and Z01CP010176; contracts CDC-200–2006-18797 and CDC-200–2015-63931 from the Centers for Disease Control and Prevention, USA; contract 90047713 from the Agency for Healthcare Research and Quality, USA; contract 90051652 from the Health Resources and Services Administration, USA; grants CBR-86906, CBR-94036, HCP-97105 and TGF-96118 from the Canadian Institutes of Health Research, Canada; Ontario Ministry of Health and Long Term Care; and the Government of Alberta, Canada. Additional support was provided by the National Cancer Institute, National Institute for Mental Health and National Institute on Drug Abuse. The COHERE study group has received unrestricted funding from: Agence Nationale de Recherches sur le SIDA et les Hépatites Virales (ANRS), France; HIV Monitoring Foundation, The Netherlands; and the Augustinus Foundation, Denmark. The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7/2007–2013) under EuroCoord grant agreement no. 260694. A list of the funders of the participating cohorts can be found at www.COHERE.org. JMM received a personal 80:20 research grant from the Institut d’Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Barcelona, Spain during 2017–19. JLC was supported by NIH grant K23AI120875. This study was also made possible by the generous support of the American people through the United States Agency for International Development (INROADS USAID-674-A-12–00029), and by a grant from the Swiss National Science Foundation (Ambizione-PROSPER PZ00P3_160407 to JB special project funding grant 174281 to ME). The content is solely the responsibility of the authors and does not necessarily represent the official views of the funders.

Conflicts of interest

FB received fees from ViiV Healthcare, Janssen, BMS, Gilead, and MSD for educational presentations and research grants from Gilead and Janssen. MJG has served as ad hoc member on Advisory HIV Boards to Merck, ViiV and Gilead in the past three years. AM has received honoraria, lecture fees, consultancy or travel support from Gilead and ViiV. MS has received research grants to his institution from Merck and Gilead. AS has received grants and travel support to her institution from ViiV Healthcare. JJV has personal fees from Merck/MSD, Gilead, Pfizer, Astellas Pharma, Basilea, German Centre for Infection Research (DZIF), University Hospital Freiburg/ Congress and Communication, Academy for Infectious Medicine, University Manchester, German Society for Infectious Diseases (DGI), Ärztekammer Nordrhein, University Hospital Aachen, Back Bay Strategies, German Society for Internal Medicine (DGIM) and grants from Merck/MSD, Gilead, Pfizer, Astellas Pharma, Basilea, German Centre for Infection Research (DZIF), German Federal Ministry of Education and Research (BMBF).

Abbreviations:aHR

adjusted hazard ratio

ART

antiretroviral therapy

CI

confidence interval

COHERE

Collaboration of Observational HIV Epidemiological Research in Europe

HIV

human immunodeficiency virus

HPV

human papillomavirus

ICC

invasive cervical cancer

IeDEA

International epidemiology Databases to Evaluate AIDS

IQR

interquartile range

pys

person-years

NNRTI

non-nucleoside reverse transcriptase inhibitor

PI

protease inhibitor

ReferencesGinsburgO, BrayF, ColemanMP, The global burden of women’s cancers: a grand challenge in global health. Lancet 2017; 389: 847–60.27814965FerlayJ, SoerjomataramI, ErvikM, GLOBOCAN 2012 v1.0, Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11 [Internet] Lyon, France: International Agency for Research on Cancer; 2013 Available from: http://globocan.iarc.fr, accessed on 18/04/2018.PalefskyJM, MinkoffH, KalishLA, Cervicovaginal Human Papillomavirus Infection in Human Immunodeficiency Virus-1 (HIV)-Positive and High-Risk HIV-Negative Women. J Natl Cancer Inst 1999; 91: 226–36.10037100DenslowSA, RositchAF, FirnhaberC, TingJ, SmithJS. Incidence and progression of cervical lesions in women with HIV: a systematic global review. Int J STD AIDS 2014; 25: 163–77.24216030PalellaFJ, DelaneyKM, MoormanAC, Declining Morbidity and Mortality among Patients with Advanced Human Immunodeficiency Virus Infection. N Engl J Med 1998; 338: 853–60.9516219BiggarRJ, ChaturvediAK, GoedertJJ, EngelsEA, HIV/AIDS Cancer Match Study. AIDS-related cancer and severity of immunosuppression in persons with AIDS. J Natl Cancer Inst 2007; 99: 962–72.17565153BlattnerWA, NowakRG. Epidemiology of AIDS-Defining Malignancies. In: Cancers in People with HIV and AIDS New York, NY: Springer New York, 2014: 17–30.KellyH, WeissHA, BenaventeY, Association of antiretroviral therapy with high-risk human papillomavirus, cervical intraepithelial neoplasia, and invasive cervical cancer in women living with HIV: a systematic review and meta-analysis. Lancet HIV 2018; 5: e45–58.29107561GangeSJ, KitahataMM, SaagMS, Cohort profile: the North American AIDS Cohort Collaboration on Research and Design (NA-ACCORD). Int J Epidemiol 2007; 36: 294–301.17213214McGowanCC, CahnP, GotuzzoE, Cohort Profile: Caribbean, Central and South America Network for HIV research (CCASAnet) collaboration within the International Epidemiologic Databases to Evaluate AIDS (IeDEA) programme. Int J Epidemiol 2007; 36: 969–76.17846055EggerM, EkoueviDK, WilliamsC, Cohort Profile: the international epidemiological databases to evaluate AIDS (IeDEA) in sub-Saharan Africa. Int J Epidemiol 2012; 41: 1256–64.21593078Asia-Pacific | IeDEA Available at: http://www.iedea.org/regions/asia-pacific; accessed on 18/04/2018.SengayiM, SpoerriA, EggerM, Record linkage to correct under-ascertainment of cancers in HIV cohorts: The Sinikithemba HIV clinic linkage project. Int J Cancer 2016; 139: 1209–16.27098265ChêneG, PhillipsA, CostagliolaD, Cohort Profile: Collaboration of Observational HIV Epidemiological Research Europe (COHERE) in EuroCoord. Int J Epidemiol 2017; 46: 797–797n.27864413LambertPC, RoystonP. Further development of flexible parametric models for survival analysis. Stata J 2009; 9: 265–90.PoleselJ, FranceschiS, SuligoiB, Cancer incidence in people with AIDS in Italy. Int J Cancer 2010; 127: 1437–1445.20049835GuiguetM, BouéF, CadranelJ, LangJ-M, RosenthalE, CostagliolaD. Effect of immunodeficiency, HIV viral load, and antiretroviral therapy on the risk of individual malignancies (FHDH-ANRS CO4): a prospective cohort study. Lancet Oncol 2009; 10: 1152–9.19818686AbrahamAG, D’SouzaG, JingY, Invasive cervical cancer risk among HIV-infected women: a North American multicohort collaboration prospective study. J Acquir Immune Defic Syndr 2013; 62: 405–13.23254153MinkoffH, ZhongY, BurkRD, Influence of adherent and effective antiretroviral therapy use on human papillomavirus infection and squamous intraepithelial lesions in human immunodeficiency virus-positive women. J Infect Dis 2010; 201: 681–90.20105077SilverbergMJ, LeydenWA, ChiA, Human Immunodeficiency Virus (HIV)- and Non-HIV-Associated Immunosuppression and Risk of Cervical Neoplasia. Obstet Gynecol 2018; 131: 47–55.29215531PatelP, HansonDL, SullivanPS, Incidence of types of cancer among HIV-infected persons compared with the general population in the United States, 1992–2003. Ann Intern Med 2008; 148: 728–36.18490686CliffordGM, TullyS, FranceschiS. Carcinogenicity of Human Papillomavirus (HPV) Types in HIV-Positive Women: A Meta-Analysis From HPV Infection to Cervical Cancer. Clin Infect Dis 2017; 64: 1228–35.28199532GossPE, LeeBL, Badovinac-CrnjevicT, The Lancet Oncology Commission Planning cancer control in Latin America and the Caribbean. Lancet Oncol 2013; 14: 391–436.23628188KinghamTP, AlatiseOI, VanderpuyeV, Treatment of cancer in sub-Saharan Africa. Lancet Oncol 2013; 14: e158–67.23561747World Health Organization. Guidelines for screening and treatment of precancerous lesions for cervical cancer prevention Geneva: World Health Organization, 2013.BruniL, DiazM, Barrionuevo-RosasL, Global estimates of human papillomavirus vaccination coverage by region and income level: a pooled analysis. Lancet Glob Health 2016; 4: e453–63.27340003KojicEM, RanaAI, Cu-UvinS. Human Papillomavirus Vaccination in HIV-infected Women: Need for Increased Coverage. Expert Rev Vaccines 2016; 15: 105–17.26599305GakidouE, NordhagenS, ObermeyerZ. Coverage of Cervical Cancer Screening in 57 Countries: Low Average Levels and Large Inequalities. PLoS Med 2008; 5: e132.18563963SigfridL, MurphyG, HaldaneV, Integrating cervical cancer with HIV healthcare services: A systematic review. PLoS One 2017; 12: e0181156.28732037DrummondJL, WereMC, ArrossiS, Wools-KaloustianK. Cervical cancer data and data systems in limited-resource settings: Challenges and opportunities. Int J Gynaecol Obstet 2017; 138 Suppl 1: 33–40.28691330

Map of countries that contributed cohort data; number of included women, by country.

Regional ICC incidence rates with 95% confidence intervals by time since ART initiation, predicted from the crude model (red) and the main adjusted model (blue) for women who initiated an NNRTI-based regimen between 2008–2014, at age 31–50 years, with a CD4 cell count of 200 cells/μl.

Characteristics of included women at ART initiation.

South AfricaLatin AmericaNorth AmericaEurope
N (%)N (%)N (%)N (%)
All women13,569 (100%)2,261 (100%)4,120 (100%)44,281 (100%)
Median follow-up time (IQR) [years]2.1 (0.8–4.1)5.0 (2.1–8.5)5.4 (2.3–9.9)4.8 (2.0–8.5)
Median age at ART initiation (IQR) [years]35.1 (29.7–41.7)35.4 (29.3–43.6)38.6 (32.4–44.8)34.6 (29.0–41.5)
Age at ART initiation [years]
 16–303,584 (26%)627 (28%)693 (17%)12,934 (29%)
 31–508,969 (66%)1,363 (60%)2,947 (72%)26,944 (61%)
 >501,016 (7%)271 (12%)480 (12%)4,403 (10%)
First-line ART regimen
 NNRTI-based12,594 (93%)1,579 (70%)1,261 (31%)16,724 (38%)
 PI-based939 (7%)623 (28%)2,476 (60%)24,325 (55%)
 Other ART36 (<1%)59 (3%)383 (9%)3,232 (7%)
Year of ART initiation
 1996–19980 (0%)26 (1%)1,600 (39%)4,260 (10%)
 1999–200351 (<1%)560 (25%)1,276 (31%)13,236 (30%)
 2004–20076,791 (50%)801 (35%)864 (21%)12,821 (29%)
 2008–20146,727 (50%)874 (39%)380 (9%)13,964 (32%)
Median CD4 cell count at ART initiation (IQR) [cells/µl]115 (50–182)178 (74–281)241 (107–385)241 (129–363)
CD4 cell count at ART initiation [cells/µl]
 < 503,113 (23%)340 (15%)577 (14%)4,676 (11%)
 50–992,402 (18%)274 (12%)303 (7%)3,271 (7%)
 100–1994,733 (35%)465 (21%)663 (16%)7,890 (18%)
 200–3491,821 (13%)601 (27%)1,047 (25%)13,122 (30%)
 350–499276 (2%)154 (7%)587 (14%)6,146 (14%)
 500–699139 (1%)74 (3%)349 (8%)3,069 (7%)
 ≥ 70066 (<1%)30 (1%)170 (4%)1,638 (4%)
 Missing1,019 (8%)323 (14%)424 (10%)4,469 (10%)
Median HIV RNA at ART initiation (IQR) [log10 copies/ml]4.4 (2.6–5.2)4.8 (4.0–5.3)4.3 (3.4–5.0)4.5 (3.7–5.1)
HIV RNA at ART initiation [log10 copies/ml]
 < 2.7788 (6%)91 (4%)608 (15%)4,704 (11%)
 2.7–3.9404 (3%)260 (11%)806 (20%)7,319 (17%)
 4.0–4.9870 (6%)585 (26%)1,218 (30%)14,051 (32%)
 ≥ 5.0984 (7%)568 (25%)946 (23%)11,764 (27%)
 Missing10,523 (78%)757 (33%)542 (13%)6,443 (15%)

ART, antiretroviral therapy; HIV, human immunodeficiency virus; IQR, interquartile range; NNRTI, non-nucleoside reverse-transcriptase inhibitor; PI, protease-inhibitor; RNA, ribonucleic acid.

Comparison of ICC rates between different regions and Europe:

Crude and adjusted HRs for ICC at 2 years and 5 years after ART initiation in women living with HIV.

At 2 yearsAt 5 years
Crude HR (95% CI)Adjusted HR* (95% CI)Crude HR (95% CI)Adjusted HR* (95% CI)
Region
 Europe1.001.001.001.00
 North America1.81 (0.94 – 3.48)1.51 (0.73 – 3.12)0.98 (0.48 – 1.99)0.79 (0.37 – 1.71)
 Latin America1.93 (1.09 – 3.42)1.83 (0.99 – 3.37)2.32 (1.24 – 4.31)2.43 (1.27 – 4.68)
 South Africa6.84 (5.20 – 9.00)6.23 (4.29 – 9.05)11.06 (7.80 – 15.68)10.66 (6.73 – 16.88)

Adjusted for CD4 cell count at ART initiation, age at ART initiation, first-line ART regimen, and calendar period of ART initiation.

ART, antiretroviral therapy; CI, confidence interval; HR, hazard ratio; ICC, invasive cervical cancer.

Adjusted hazard ratios for the effect of different factors on the risk of developing incident ICC in women who have initiated ART.

Hazard ratio* (95% CI)p-value for interaction**
CD4 cell count at ART initiation0.76
 per 100 cells/µl decrease1.25 (1.15 – 1.36)
Age at ART initiation [years]0.34
 16–301.00
 31–501.38 (1.05 – 1.81)
 >501.57 (1.03 – 2.40)
First-line ART regimen0.21
 NNRTI-based1.00
 PI-based1.05 (0.79 – 1.41)
 Other ART0.57 (0.27 – 1.18)
Year of ART initiation0.13
 1996–19981.49 (0.92 – 2.42)
 1999–20031.19 (0.80 – 1.77)
 2004–20070.83 (0.61 – 1.14)
 2008–20141.00

Adjusted for region, CD4 cell count at ART initiation, age at ART initiation, calendar year of ART initiation, and first-line ART regimen.

Derived from likelihood ratio test comparing the adjusted model with and without the interaction of a specific variable with region.

ART, antiretroviral therapy; CI, confidence interval; ICC, invasive cervical cancer; NNRTI, non-nucleoside reverse-transcriptase inhibitor; PI, protease-inhibitor

Novelty and impact:

Our study compares invasive cervical cancer (ICC) rates among women living with HIV who initiated antiretroviral therapy (ART) across different geographic regions. It shows that ICC incidence rates are particularly high in women living with HIV in South Africa or Latin America. ICC prevention through early ART initiation, equitable access to effective cervical cancer screening, and promotion of human papillomavirus vaccination should be key elements of global efforts to reduce cancer-related health inequities.