Precancerous Lesion Progression and Regression
Using empirical data from HIPVIRG,14 we calculated precancer progression and regression probabilities for each nadir CD4+ category and adjusted to preserve nadir CD4+ relationships (see Web Appendix, http://links.lww.com/QAI/B174). We allowed baseline progression and regression probabilities to be the same across HPV genotypes and relied on the calibration process to adjust the baseline probabilities for each HPV genotype.
Among HIPVIRG participants, 2 individuals with nadir CD4+ <200 (1345 person-months) who did not receive treatment for AIN2/3 progressed to invasive cancer, resulting in a monthly progression probability of 0.15% (CI: 0.02% to 0.54%). We set the baseline probability of progressing to cancer for individuals with a nadir CD4+ ≥200 to half (reflecting slower progression) of the baseline probability of progressing to cancer for individuals with a nadir CD4+ <200.16
Anal Cancer Progression, Detection, and Mortality
Among the men who develop invasive anal cancer, their cancer may progress to more severe stages (classified using AJCC sixth edition staging grouped by stage I; stages II/III; and stage IV), facing an increasing probability of being detected through symptoms in each advanced stage (Table 2).21 In addition to background mortality due to HIV and other causes, individuals who develop anal cancer face excess mortality due to anal cancer death. In the ART-era, anal cancer prognosis for an HIV-infected individual is similar to an HIV-negative individual22–24; therefore, we used population-based data from the National Cancer Institute's (NCI) Surveillance, Epidemiology, and End Results (SEER) program for cancer registries (ie, not HIV-specific). We extracted the 5-year cancer cause-specific25 survival probabilities stratified by stage and time since diagnosis (year 1, years 2–3, and years 4–20) to derive the monthly probability of dying from anal cancer (Web Appendix Table 2, http://links.lww.com/QAI/B174).
HIV-infected individuals face a monthly probability of initiating ART that increases as nadir CD4+ count declines (Web Appendix Table 3, http://links.lww.com/QAI/B174). In 2002 at the baseline of the HIPVIRG cohort, 93% of individuals had initiated ARTs. In the absence of ART initiation, individuals may progress to lower CD4+ counts over time26 and face elevated CD4+-specific risks of dying from HIV-related causes (Web Appendix Table 3, http://links.lww.com/QAI/B174). On initiating ART (to ensure the model-simulated cohort corresponds to that of the HIPVIRG cohort), individuals stay within their CD4+ nadir category for their remaining lifetime but face a lower HIV-related mortality (remains slightly elevated compared with the general population27).
After identification of baseline inputs of the model, we calibrated (or fit) the model for specific model transitions that were unobservable or varied from setting to setting. We identified and allowed key uncertain parameter inputs to vary to characterize joint uncertainty in the natural history of anal cancer and to ensure good fit to empirical data from HIV-infected MSM study populations in North America. In particular, we relied on calibration to infer HPV genotype–specific variations in age- and CD4+ nadir–specific input parameters.
After a similar calibration approach used for a natural history model of HPV-induced cervical carcinogenesis,28,29 we identified plausible ranges around baseline input parameter values and repeated model simulations by drawing uniformly across the predefined parameter search space for each model run, resulting in unique combinations of natural history parameters. For specific parameters, we constrained the search space to preserve well-understood relationships between HPV genotypes (Table 2; Web Appendix, http://links.lww.com/QAI/B174). To identify the parameter sets that maximized correspondence between model-projected outcomes and empirical data (ie, calibration targets), we used a likelihood-based approach that involved calculating the likelihood score for each unique parameter set using a likelihood ratio test to identify “good-fitting” parameter sets.
The calibration target data (Web Appendix Table 4, http://links.lww.com/QAI/B174) included data formats not directly used as model inputs (eg, age-, HPV genotype–, and CD4+–specific data). We used genotype-specific prevalence of HPV infection by age, prevalence of AIN1 and AIN2/3 precursors by age, and the frequency of HPV genotypes in AIN1 and AIN2/3 from the HIPVIRG cohort study.14 To mimic the 6-month detection rate schedule in the HIPVIRG cohort study, we simulated a population similar to the study population; specifically, model detection rates reflected biannual screening visits. Screening, however, did not involve the treatment of precursors or screen detection of anal cancer. In addition, we used data from the most recent meta-analysis evaluating HPV-type distribution in anal cancer tissues among HIV-infected men, which showed a lower contribution of HPV-16 to anal cancers among HIV-infected men compared with HIV-uninfected men (ie, 57% vs. 77%, respectively).30 Age- and CD4+ nadir–specific anal cancer incidence rates were derived from the multicohort collaboration North American AIDS Cohort Collaboration on Research and Design (NA-ACCORD) study (see Web Appendix, http://links.lww.com/QAI/B174). For calibration, we restricted the analytic cohort to cancers diagnosed among men aged 30 years and older between 2000 and 2007 (n = 57 cases and 47,366 person-years) to correspond with the HIPVIRG longitudinal trial period.
After calibration, we selected a sample (ie, 50) of the “good-fitting” parameter sets with the highest aggregate likelihood score to capture probabilistic uncertainty in the natural history process. We used the sample of “good-fitting” sets to assess the face validity and projective validity of the model fit to external data not used in the calibration process, such as age of cancer diagnosis, and life expectancy by CD4+ nadir category. Model outcomes reflect projections for men diagnosed with HIV before 2002, which capture the distribution of nadir CD4+ counts and ART status initiations reported in the HIPVIRG cohort study.
Model outcomes were reported as the median, minimum, and maximum values across the 50 “good-fitting” sets. Finally, we used the good-fitting parameter sets to infer the posterior (ie, postcalibration process) progression probabilities for the progression from AIN2/3 to invasive anal cancer by CD4+ nadir and HPV genotype.
Among the 200,000 randomly sampled parameter sets, the 50 “good-fitting” sets that produced the highest agreement to our 55 calibration targets, generally yielded reasonable fits to empirical data on age-specific HPV prevalence (Fig. 2A); age-specific prevalence of AIN1 and AIN2/3 (Web Appendix Figure 1, http://links.lww.com/QAI/B174); genotype distribution of HPV infections in AIN1, AIN2/3, and cancer (Fig. 2B); age- and CD4+ nadir–specific anal cancer incidence (Fig. 2C); and stage of cancer detection (Web Appendix Figure 2, http://links.lww.com/QAI/B174). Cancer incidence rates were three- to six-fold higher among men with CD4+ nadir <200 compared to men with CD4+ nadir >500. For example, the median cancer incidence peaked at 176 per 100,000 men among individuals aged 45–49 years with CD4+ nadir <200 and 28 per 100,000 men among individuals aged 45–49 years with CD4+ nadir >500.
Among these “good-fitting” sets, the model estimated a lifetime risk of developing detected anal cancer of 2.5% (range: 0.4%–5.7%) with a mean age of cancer diagnosis of 46.1 years (range: 44.9–47.3 years) and a median diagnosis of 44.5 years (range: 43.0–46.0 years), which are consistent with empirical data from other HIV-infected MSM cohorts that reported cumulative incidence of 1.7% (2005–2009),31 and mean ages of 45.422 and 49 years,21 or a median age of 38 years1 (Table 3). Using anal cancer survival probabilities derived from SEER registries, we projected anal cancer mortality rates that varied by age; median age-specific rates ranged from 9 to 42 per 100,000 HIV-infected MSM (Web Appendix Figure 4, http://links.lww.com/QAI/B174).
Among all men and irrespective of CD4+ nadir, the total (undiscounted) life expectancy was 59.7 years (range: 59.5–59.8) but also varied by CD4+ nadir category. For example, the life expectancy among men with CD4+ nadir <200 was 53.8 years (range: 53.5–54.0), whereas the life expectancy among men with CD4+ nadir >500 was 72.8 years (range: 72.6–72.9); these projections were generally consistent with life expectancies reported from HIV-infected MSM in NA-ACCORD32 and Kaiser Permanente33 before 2006.
Progression to Anal Cancer
After calibration, the resulting median monthly progression probability across the top-fitting parameter sets differed substantially by HPV genotype and nadir CD4+ status. For example, the median monthly progression probability from an HPV-16 AIN2/3 to invasive cancer was 0.027% (ie, 3.19% 10-year probability) among individuals with CD4+ nadir >500 and 0.054% (ie, 6.28% 10-year probability) among individuals with CD4+ nadir <200 (See Web Appendix Figure 3, http://links.lww.com/QAI/B174). By contrast, the median monthly progression probability from an AIN2/3 to invasive cancer without detectable HPV was 0.002% (ie, 0.24% 10-year probability) among individuals with CD4+ nadir >500 and 0.004% (ie, 0.48% 10-year probability) among individuals with CD4+ nadir <200.
We developed an individual-based natural history simulation model by synthesizing primary and secondary data from the most recent North American studies of anal HPV-related neoplasia in HIV-infected MSM to address alternative anal cancer screening approaches in a timely manner. This calibrated model can be used to superimpose various existing and proposed scenarios of anal cancer screening to determine the long-term population-level health and economic consequence, often unavailable from empirical trials, as well as to determine which screening strategies provide the most value for money. Our natural history model, which was consistent with empirical data from other HIV-infected MSM anal cancer studies, also provides important insights on the progression potential of high-grade precancers to invasive disease by HPV genotype and CD4+ nadir. Such endpoints are difficult to observe in empirical studies and have not been directly observed to date.
Several groups17,18,34 have developed disease simulation models to evaluate the value of screening for anal precancers among HIV-infected MSM. Among these studies, the range of recommendations included annual anal Pap screening, which would trigger anoscopy and biopsy if abnormal in the United States17 to no screening in England.18 However, these studies did not account for independent HPV genotypes, and the potential sustained role of nadir CD4+ after ART initiation on disease progression. In addition, Lam et al34 did not evaluate long-term outcomes such as cancer or life expectancy. A recent study by Ong et al19 developed a Markov model evaluation of the cost effectiveness of digital anal rectal examination (DARE) every 1–5 years; however, the model did not include precancerous health states and made simplifying assumptions around the age-specific probability of developing anal cancer, ignoring important age differences in cancer risk between the general population and HIV-infected MSM.
To our knowledge, the model we describe here is the first with capabilities to evaluate nuanced screening and vaccination policies that captures complex interactions between anal carcinogenesis and HIV natural history, which are important to capture when evaluating potential anal cancer screening strategies. To date, no other anal cancer models have incorporated HPV genotype–specific data from large cohort studies, nor used the North American cohorts of HIV and anal cancer diagnosis to inform model transition probabilities.
There are limitations to our model structure, data availability, and calibration approach. Our model is inherently tied to the empirical data that informed it. Similar to many natural history studies, our model used the composite outcome using cytology/histology results as true underlying health state. Although cytology alone may not be highly sensitive, high-resolution anoscopy was performed on all men. In addition, HIV-infected MSM often have multiple AIN lesions, and tracking the status of individual lesions is difficult in clinical practice. For example, a new AIN2/3 diagnosis could either be a previously detected lesion that progressed, or a new lesion that appeared, which may have been truly absent or simply missed at the previous visit. To minimize missing lesions, the HIPVIRG study required chart review where the worst lesion seen was confirmed. In addition, anoscopists completed at each session a diagram to indicate changes in the lesions they previously saw. Similarly, many HIV-infected MSM have multiple HPV infections,7,14 and microdissection to determine the HPV genotype attributable to individual AIN1s and AIN2/3s was not performed. In our model, we assumed each HPV type was acquired and progressed independently, driven by the HPV incidence probabilities observed in the HIPVIRG cohort study.
There may be other mechanisms by which ART affects HPV natural history other than through preventing a lower absolute CD4+ nadir. For example, ART initiation affects HPV natural history, most notably through immune reconstitution. It is possible that absolute CD4+ nadir, length of time at CD4+ nadir, timing of this CD4+ nadir in relation to the present, and extent of immune reconstitution (ie, current CD4+) all impact probability of progression from AIN2/3 to invasive cancer. The Swiss Cohort study9 suggested that CD4+ nadir, CD4+ count 6–7 years before present, and current CD4+ count were each associated with cancer incidence; however, the authors did not present a multivariable analysis, so it is unclear which factors are independently associated with cancer incidence, based on their data. Kesselring et al35 published an adjusted hazard ratio for anal cancer of 1.52 per year of CD4+ count below 200, suggesting that it is not just absolute CD4+ nadir, but rather time at CD4+ nadir, that has an impact on anal cancer incidence. Although our model attempts to capture the role nadir CD4+ count plays in the risk of developing anal cancer and does so using what the literature suggests is likely the most robust clinical predictor of progression to invasive cancer, ie, absolute CD4+ nadir,8,9,11–14,16,35–38 the model is not able to comprehensively integrate the complex—and not fully understood—interplay between CD4+ count through time and risk of anal cancer. In addition, we reduced the baseline probability of progressing from AIN2/3 to invasive cancer for individuals with a CD4+ nadir ≥200 by half compared to individuals with a CD4+ nadir <200 but acknowledge that the reference group from the multivariable analysis from D'Souza and colleagues16 was compared with total cancer incidence irrespective of lesion status.
No data are available to inform the progression between stages of cancer and symptom stage detection probabilities. We used a panel of experts from Harvard-affiliated teaching hospitals to inform baseline estimates of the duration between stage I, stages II/III, and stage IV cancers, and calibrated these baseline inputs to fit stage distribution observed in SEER among the general population (Web Appendix, http://links.lww.com/QAI/B174). Extensive sensitivity analysis will be performed on all policy analyses to evaluate the robustness of screening policy and cancer progression and detection assumptions. In addition, we assumed that, for undetected cancers, progression to later stages was independent of CD4+ nadir. We plan to explore the role of differential cancer progression by CD4+ nadir in future analyses.
Our model does not explicitly capture the dynamic, sexual patterns of HPV acquisition but relies on age-specific incidence calibrated to observed prevalence to serve as a proxy for sexual behavior. HPV prevalence from the HIPVIRG cohort suggests that the HPV prevalences simulated in our model are representative of other HIV-infected MSM populations.7 In addition, we assume constant transition probabilities that do not capture the role of persistent HPV infections and lesions. Study cohort size restricted our ability to stratify HPV clearance, and lesion progression and regression by duration; however, ongoing longitudinal studies such as SPANC39,40 and ANCHOR (Clinical trials identifier: NCT02135419) may provide insight once trial data are available.
Using primary and secondary data from North American HIV-infected MSM cohorts and registries, we synthesized existing evidence into a single natural history model that reflects the most recent understanding of anal carcinogenesis, accounting for interactions between HPV type-specific infections, CD4+ lymphocyte strata, and ART status. Such models can be used to evaluate the population-level health and economic consequences of complex national screening policies currently on the table in the United States.
The authors thank Catherine Regan at the Harvard Center for Health Decision Science for her help in processing model simulations and output. They also thank Jennifer Y. Wo, MD, and Theodore S. Hong, MD from Massachusetts General Hospital Cancer Center and Harvard Medical School (HMS), Harvey J. Mamon, MD, PhD, and Neil E. Martin, MD, MPH, from Dana-Farber/Brigham and Women's Cancer Center and HMS, and Anand Mahadevan, MD, from Geisinger Cancer Center (formerly at Beth Israel Deaconess Medical Center and HMS) for their clinical input and expertise regarding the behavior of invasive anal cancer.
1. Silverberg MJ, Lau B, Justice AC, et al. Risk of anal cancer in HIV-infected and HIV-uninfected individuals in North America. Clin Infect Dis. 2012;54:1026–1034.
2. Grulich AE, van Leeuwen MT, Falster MO, et al. Incidence of cancers in people with HIV/AIDS compared with immunosuppressed transplant recipients: a meta-analysis. Lancet. 2007;370:59–67.
3. Shiels M, Cole S, Kirk G, et al. A meta-analysis of the incidence of non-AIDS cancers in HIV-infected individuals. J Acquir Defic Syndr. 2009;52:11.
4. Bedimo R, McGinnis K, Dunlap M, et al. Incidence of non-AIDS-defining malignancies in HIV-infected versus noninfected patients in the HAART era: impact of immunosuppression. J Acquir Immune Defic Syndr. 2009;52:5.
5. Wells JS, Holstad MM, Thomas T, et al. An Integrative Review of Guidelines for anal cancer screening in HIV-infected persons. AIDS Patient Care STDs. 2014;28:350–357.
6. Berry JM, Jay N, Cranston RD, et al. Progression of anal high-grade squamous intraepithelial lesions to invasive anal cancer among HIV-infected men who have sex with men. Int J Cancer J Int du Cancer. 2014;134:1147–1155.
7. Machalek DA, Poynten M, Jin F, et al. Anal human papillomavirus infection and associated neoplastic lesions in men who have sex with men: a systematic review and meta-analysis. Lancet Oncol. 2012;13:487–500.
8. Palefsky J, Holly E, Ralston M, et al. Anal squamous intraepithelial lesions in HIV-positive and HIV-negative homosexual and bisexual men: prevalence and risk factors. Acquir Immune Defic Syndr Hum Retrovirol. 1998;17:6.
9. Bertisch B, Franceschi S, Lise M, et al. Risk factors for anal cancer in persons infected with HIV: a nested case-control study in the Swiss HIV Cohort Study. Am J Epidemiol. 2013;178:877–884.
10. Wilkin TJ, Palmer S, Brudney KF, et al. Anal intraepithelial neoplasia in heterosexual and homosexual HIV-positive men with access to antiretroviral therapy. J Infect Dis. 2004;190:1685–1691.
11. Powles T, Robinson D, Stebbing J, et al. Highly active antiretroviral therapy and the incidence of non–AIDS-defining cancers in people with HIV infection. J Clin Oncol. 2009;27:884–890.
12. Egger M, May M, Chêne G, et al. Prognosis of HIV-1-infected patients starting highly active antiretroviral therapy: a collaborative analysis of prospective studies. Infect Dis Clin Pract. 2002;11:10.
13. Crum-Cianflone NF, Hullsiek KH, Marconi VC, et al. Anal cancers among HIV-infected persons: HAART is not slowing rising incidence. AIDS. 2010;24:8.
14. de Pokomandy A, Rouleau D, Ghattas G, et al. HAART and progression to high-grade anal intraepithelial neoplasia in men who have sex with men and are infected with HIV. Clin Infect Dis. 2011;52:1174–1181.
15. Palefsky JM, Holly EA, Efirdc JT, et al. Anal intraepithelial neoplasia in the highly active antiretroviral therapy era among HIV-positive men who have sex with men. AIDS. 2005;19:1407–1414.
16. D'Souza G, Wiley DJ, Li X, et al. Incidence and epidemiology of anal cancer in the multicenter AIDS Cohort Study (MACS). J Acquir Immune Defic Syndr. 2008;48:491–499.
17. Goldie SJ, Kuntz KM, Weinstein MC, et al. The clinical effectiveness and cost-effectiveness of screening for anal squamous intraepithelial lesions in homosexual and bisexual hiv-positive men. JAMA. 1999;281:1822–1829.
18. Czoski-Murray C, Karnon J, Jones R, et al. Cost-effectiveness of screening high-risk HIV-positive men who have sex with men (MSM) and HIV-positive women for anal cancer. Health Technology Assess. 2010;14:131.
19. Ong JJ, Fairley CK, Carroll S, et al. Cost-effectiveness of screening for anal cancer using regular digital ano-rectal examinations in men who have sex with men living with HIV. J Int AIDS Soc. 2016;19:20514.
20. de Pokomandy A, Rouleau D, Ghattas G, et al. Prevalence, clearance, and incidence of anal human papillomavirus infection in hiv-infected men: the HIPVIRG cohort study. J Infect Dis. 2009;199:965–973.
21. Read T, Huson K, Millar J, et al. Size of anal squamous cell carcinomas at diagnosis: a retrospective case series. Int J STD AIDS. 2013;24:879–882.
22. Abramowitz L, Mathieu N, Roudot-Thoraval F, et al. Epidermoid anal cancer prognosis comparison among HIV+ and HIV− patients. Aliment Pharmacol Ther. 2009;30:414–421.
23. Chiao EY, Giordano TP, Richardson P, et al. Human immunodeficiency virus–associated squamous cell cancer of the anus: epidemiology and outcomes in the highly active antiretroviral therapy era. J Clin Oncol. 2008;26:474–479.
24. Marcus JL, Chao C, Leyden WA, et al. Survival among HIV-infected and HIV-uninfected individuals with common non–AIDS-defining cancers. Cancer Epidemiol Biomarkers Prev. 2015;24:1167–1173.
25. Howlader N, Ries LA, Mariotto AB, et al. Improved estimates of cancer-specific survival rates from population-based data. JNCI J Natl Cancer Inst. 2010;102:1584–1598.
26. Lodi S, Phillips A, Touloumi G, et al. Time from human immunodeficiency virus seroconversion to reaching CD4+ cell count thresholds <200, <350, and <500 cells/mm3
: assessment of need following changes in treatment guidelines. Clin Infect Dis. 2011;53:817–825.
27. McManus H, O'Connor CC, Boyd M, et al. Long-term survival in HIV positive patients with up to 15 years of antiretroviral therapy. PLoS One. 2012;7:e48839.
28. Campos NG, Burger EA, Sy S, et al. An updated natural history model of cervical cancer: derivation of model parameters. Am J Epidemiol. 2014;180:545–555.
29. Kim JJ, Kuntz KM, Stout NK, et al. Multiparameter calibration of a natural history model of cervical cancer. Am J Epidemiol. 2007;166:137–150.
30. Clifford GM, Lin C, Franceschi S. A meta-analysis of HPV
type distribution across the full spectrum of Anal disease, by HIV status. Paper presented at: international papillomavirus; 2017; Cape Town, South Africa.
31. Silverberg MJ, Lau B, Achenbach CJ, et al. Cumulative incidence of cancer among persons with HIV in North America: a cohort study. Ann Intern Med. 2015;163:507–518.
32. Samji H, Cescon A, Hogg RS, et al. Closing the gap: increases in life expectancy among treated HIV-positive individuals in the United States and Canada. PLoS One. 2013;8:e81355.
33. Marcus JL, Chao CR, Leyden WA, et al. Narrowing the gap in life expectancy between hiv-infected and hiv-uninfected individuals with access to care. J Acquir Immune Defic Syndr. 2016;73:39–46.
34. Lam JM, Hoch JS, Tinmouth J, et al. Cost-effectiveness of screening for anal precancers in HIV-positive men. AIDS. 2011;25:635–642.
35. Kesselring A, Gras L, Smit C, et al. Immunodeficiency as a risk factor for non-AIDS-defining malignancies in HIV-1-infected patients receiving combination antiretroviral therapy. Clin Infect Dis. 2011;52:1458–1465.
36. Richel O, Van Der Zee RP, Smit C, et al. Brief report: anal cancer in the hiv-positive population: slowly declining incidence after a decade of cART. J Acquir Immune Defic Syndr. 2015;69:602–605.
37. Chaturvedi AK, Madeleine MM, Biggar RJ, et al. Risk of human papillomavirus-associated cancers among persons with AIDS. J Natl Cancer Inst. 2009;101:1120–1130.
38. Legarth R, Helleberg M, Kronborg G, et al. Anal carcinoma in HIV-infected patients in the period 1995–2009: a Danish nationwide cohort study. Scand J Infect Dis. 2013;45:453–459.
39. Machalek DA, Grulich AE, Hillman RJ, et al. The Study of the Prevention of Anal Cancer (SPANC): design and methods of a 3-year prospective cohort study. BMC public health. 2013;13:946.
40. Duncan KC, Chan KJ, Chiu CG, et al. HAART slows progression to anal cancer in HIV-infected MSM. AIDS. 2015;29:305–311.
HPV; anal neoplasms; decision analysis; mathematical model
Supplemental Digital Content
Copyright © 2018 Wolters Kluwer Health, Inc. All rights reserved.