HPV types 16/18 L1 E6 and E7 proteins seropositivity and cervical cancer risk in HIV-positive and HIV-negative black South African women

Background In populations with high rates of human immunodeficiency virus (HIV)-coinfection, the nature of the relationship between human papillomavirus (HPV)-16 and -18 (L1, E6 and E7) antibodies and cervical cancer is still uncertain. We measured the association between seropositivity to HPV (L1, E6 and E7) proteins and cervical cancer among black South African women with and without HIV co-infection. Methods We used questionnaire data and serum collected from consecutively recruited patients with a newly diagnosed cancer from the Johannesburg Cancer Study from 1346 cervical cancer cases and 2532 controls (diagnosed with other non-infection related cancers). Seropositivity to HPV proteins was measured using a multiplex serological assay based on recombinant glutathione S-transferase (GST) fusion proteins. We measured associations between their presence and cervical cancer using unconditional logistic regression models and evaluated the sensitivity and specificity of these HPV biomarkers. Results Among controls, HIV-negative women from rural areas compared to urban had significantly higher HPV seroprevalence, HPV16 E7 (8.6% vs 3.7%) and HPV18 E7 (7.9% vs 2.0%). HPV16 E6 and E7 antibodies were positively associated with cervical cancer in HIV-positive (Adjusted Odds Ratio (AOR) = 33; 95% CI 10–107) and HIV-negative women (AOR = 97; 95% CI 46–203). In HIV-positive women, HPV E6/E7 antibodies had low sensitivity (43.0%) and high specificity (90.6%) for cervical cancer detection. In HIV-negative women, HPV E6/E7 antibodies sensitivity was 70.6% and specificity was 89.7%. Conclusions Our data show that HPV (L1, especially E6 and E7) antibody positivity is associated with cervical cancer in both HIV-positive and HIV-negative women. Nonetheless, being HIV-positive plays an important role in the development of cervical cancer. Supplementary Information The online version contains supplementary material available at 10.1186/s13027-022-00418-2.


Introduction
Globally, in 2018, about 70% of cervical cancers were attributable to high-risk human papillomavirus (HPV) types 16 and 18 [1]. In South Africa, HPV16 and 18 are Open Access *Correspondence: elviras@nicd.ac.za 1 National Cancer Registry, National Health Laboratory Service, Sandringham, Johannesburg, South Africa Full list of author information is available at the end of the article important causes of cervical cancer and are included in the currently available HPV vaccine.
The HPV early (E6 and E7) oncoproteins and the late (L1) protein are encoded by the HPV genome. The HPV oncoproteins play an important role in the tumorigenesis of cervical cancer [2]. Humoral immune response against major capsid late protein (L1) is generated during infection with HPV; which is a marker of past or present infection [3]. Furthermore, integration of the HPV genome in the host cell results in overexpression of E6 and E7 oncoproteins, which are involved in the transformation and progression of cervical and other HPV-related cancers [4]. Different serological assays have been used to screen for HPV-16 and -18 (E6 and E7) [5][6][7][8] oncoprotein antibodies. Serological markers for HPV-antibodies provide useful data about past and current HPV infections and their relationship to cervical cancer [9,10].
Less is known about the role of high-risk HPV oncoproteins among cervical cancer patients in black African women who are HIV-positive and HIV-negative. In South Africa, despite high HIV prevalence (24.1%) among black South African women [15], studies on HPV16 and 18 (E6 and E7) and late (L1) seropositivity and the risk of cervical cancer are still lacking. We therefore assessed the seropositivity of HPV -16 and -18 early (E6 and E7), and late (L1) antibodies among cervical cancer cases and infection unrelated cancer controls using a multiplex HPV serology assay [16], and estimated their association with cervical cancer risk among HIV-positive and HIVnegative women.

Setting and participants
The study population were participants recruited into the Johannesburg Cancer Study (JCS) between 1995 and 2016. The aims of JCS included measuring the relative importance of known and emerging risk factors for cancer in a black African population in Johannesburg, South Africa. The details of the JCS have been described elsewhere [17]. Briefly, the JCS collected serum samples from 26,000 consecutive consenting black South African patients with newly diagnosed cancers (> 90% histopathology confirmed), referred to the medical oncology and radiation therapy wards of (the main tertiary public) Charlotte Maxeke Johannesburg Academic Hospital and associated referral hospitals and clinics. Self-reported data on demographics and key lifestyle risk factors were collected using a structured questionnaire. The JCS also collected venous blood in one serum separation tube (SST) and one Ethylenediaminetetraacetic acid (EDTA) tube. SST serum was separated from whole blood within two days of sample collection and the sample was divided into a maximum of 4 aliquots. HIV testing was done using the Vironostika HIV Ag/Ab kit [18]. The study was approved by the University of the Witwatersrand Human Research Ethics Committee (Medical) (certificate number: M200252).

Cases and controls
The JCS is amenable to a case-control design and analysis by selecting controls that are unrelated to the exposures of interest. The current study was restricted to women aged 25 to 54 years because a large proportion of the older women did not have serology samples available. Cases were 1346 women with newly diagnosed cervical cancer (C53). Controls were 2532 women with newly diagnosed infection unrelated cancers (breast (C50) (n = 1953), colon (C18-20) (n = 197), oesophagus (C15) (n = 48), endometrium (C54-55) (n = 75), lung (C33-34) (n = 76), pancreas (C25) (n = 21) and minor cancer types (n = 162) (Additional file 1: Table S1)). Figure 1 outlines the criteria used to select cases and controls.

Serology data and laboratory methods
Serum aliquots were stored at − 25 °C and one aliquot was shipped on dry ice to the German Cancer Research Center (Deutsches Krebsforschungszentrum (DKFZ)) in Heidelberg, for antibody testing for HPV16 and 18 (L1, E6 and E7) using a multiplex serological assay. This was based on a glutathione S-transferase (GST) capture immunoassay in combination with fluorescent beads on a Luminex platform [16]. The final net (bead and GST background subtracted) Median Fluorescence Intensity (MFI) values of all samples were analysed at a dilution of 1:1000. Previously established standard cut-offs were generated at 1:100 serum dilution [19,20] and were thus not applicable. Using the Visual Inflection Point (VIP) method [

Statistical analysis
Data on demographic characteristics were summarized using frequencies and percentages for the categorical variables, medians and the interquartile range (IQR) for the continuous variable that was not normally distributed. A Pearson's Chi-squared test (for categorical variables) and t-test (for age) were used for the comparison of demographics between cases and controls. To assess cross-reactivity between the proteins, a tetrachoric correlation coefficient rank test [22] with Bonferroni adjustments (to take account of multiple comparisons) was used to analyze the correlation between each HPV16 and 18 (L1, E6 and E7) antibodies and HIV antibodies (See Additional file 1: Fig. S2). In the control group, we calculated the seroprevalence of antibodies to HPV16 and 18 (E6 and E7) oncoproteins and L1 protein across different demographic and lifestyle factors, such as age (25-34, 35-44, 45-54 years), place of residence (rural, urban), period of interview (1995-1999, 2000-2004, 2005-2009, 2010-2016), marital status (never married, married, previously married), educational level (none, primary, secondary and above), number of sexual partners (0-1, 2-5, 6+) HIV-status (negative and positive) and parity (0, 1, 2, 3, 4+). We have previously shown that cancer types unrelated to the exposure of interest, resemble background population prevalence for that exposure [18,23]. We performed a test for heterogeneity on categorical variables and a score test for trends in proportions for ordinal categorical variables. On the assumption that the selected controls should have similar seroprevalences to each other, we compared heterogeneity in seroprevalence of HPV16 and 18 (L1, E6 and E7) antibodies across different control cancer types using the Cochran-Mantel-Haenszel test (See Additional file 1: Fig. S1).
We used an unconditional logistic regression model to assess the association by calculating adjusted Odds Ratios (AOR) and 95% confidence intervals between HPV16 and 18 (E6 and E7) and L1 proteins seropositivity and the risk of cervical cancer. We stratified by HIV status to assess if the seroprevalence of the HPV proteins differs in HIV-positive and HIV-negative patients. We adjusted for age, education level, number of sexual partners, marital status, period of interview and place of residence. To assess whether the clinical performance of HPV-16 and -18 antibodies as a diagnostic marker for cervical cancer we performed the same in HIV-positive and HIVnegative cancer patients. We calculated the sensitivity, specificity of the various HPV antibodies for detecting cervical cancer using the "diagti" command in STATA. In addition, the area under the ROC curve (AUC) of the receiver operator characteristics was computed to compare the performance of the best combination of HPV E6 and E7 antibodies to discriminate cervical cancer from controls (See Additional file 1: Table S4). All the statistical tests were computed using STATA software version 16.0 (Stata Corp, college station, Tx) and SAS v 9.4 (SAS Institute, Cary, NC). All tests were considered significant at a two-tailed alpha of 5%.
In general, we did not observe differences in seroprevalence to HPV types 16 and 18 L1 and (E6 and E7) antibodies by demographic and sexual/reproductive factors (p value > 0.05). Among HIV-negative controls, the prevalence of HPV16 E7 antibody was higher among women that lived in rural areas (8.6%) as compared to women who lived in urban areas (3.7%) (p value for heterogeneity = 0.005) ( Table 2). A similar pattern was observed for HPV18 E7 (rural (7.9%) vs urban (2.0%); p value = 0.001). Among HIV-positive controls, the prevalence of HPV16 E6 &E7 antibodies decreased with an increase in number births up to 2 births (p value = 0.006) ( Table 3).

Discussion
In this study, we assessed antibody seropositivity to each of HPV-16 and -18 oncoproteins (E6 and E7), and L1 protein and calculated their association with cervical cancer. While the seroprevalence of HPV16 L1 in the previous JCS study among HIV-seronegative women was higher [14] than in this study, ORs between HPV16 L1 and cervical cancer were about the same (overall crude OR 2.2; 95% CI 1.8-2.6, vs 2.1; 1.5-2.7, Fig. 3). Previous data on HPV16 L1 seroprevalence in controls reflect different assay methods and cutoffs used, and perhaps some cross-reactivity with other HPV types. Our key findings, using larger sample size, more recent serological methods measuring a range of high-risk HPV types showed the important role of HPV-16 and -18 (E6 and E7) oncoproteins in cervical cancer risk in both HIV-positive and HIV-negative cancer patients. HPV-16 and -18 (E6 and E7) oncoproteins were associated in cervical cancer risk in both HIV-positive and HIV-negative cancer patients. HIV-positive patients compared to HIV-negative patients had lower sensitivity but with both high specificity for HPV16 and HPV 18(E6 and E7). Findings from our study confirmed that, in HIV-negative controls, HPV16 E7 and HPV18 E7 seroprevalence differed by place of residence. Among HIV-positive women, HPV16 E6&E7 seroprevalence increased with an increase in number of births.
The higher seroprevalence of HPV16 E6, HPV16 E7 and HPV18 E7 in women from rural areas might indicate a greater number of women with early undetected cervical or related cancer lesions. In rural areas of South Africa, women have poorer access to health care services, are less likely to present themselves to cervical cancer screening services and have less access to cancer diagnostic facilities, which are mainly found in tertiary hospitals in urban areas [24].
The high seroprevalence of HPV16 E6 (35.9%) antibodies in cervical cancer cases found in our study are similar to the seroprevalence of HPV16 E6 (32%) in cervical cancer cases reported from the United Kingdom (35%) and India and Algeria (32%) [5,10]. We found a relatively low      [25]. Our study demonstrates that the high seroprevalence of HPV16 E6 and E7 antibodies in cervical cancer development is in accordance with other studies from Russia and Italy [8,26]. Antibodies to the E6 and E7 oncoproteins are late markers of invasive cervical cancer [27]. In our study, seropositivity to each of the HPV-16 and -18 (E6 and E7) antibodies was associated with high risks of cervical cancer. Notably, seropositivity to HPV16 (E6 & E7) and HPV18 (E6 & E7) exhibited the highest ORs for cervical cancer (Fig. 3). Similar findings were reported in a casecontrol study of patients recruited from India and Algeria [10]. Our findings support those from other cross-sectional, case-control and prospective studies, albeit each using different assays [27][28][29][30].
Individuals with suppressed immune systems are susceptible to chronic infection, including HPV [31]. Thus, HIV infection increases the risk of HPV persistence [32]. Existing evidence of the association between HPV (E6 and E7) and L1 proteins antibodies and HIV has come from studies on men who have sex with men [10,33,34]. In contrast, most of the existing studies on the association between antibodies to HPV-16 and -18 L1, E6 and E7 proteins and cervical cancer risk have been conducted in low HIV prevalence settings [6,8,13,27]. In our study, access to HIV data made it possible to stratify our analysis by HIV status. Our finding showed that antibody seropositivity to each of HPV-16 and -18 (E6 and E7) oncoproteins was strongly associated with the risk of cervical cancer among HIV-positive women. Our data support the hypothesis that HIV plays an important role in the persistence of HPV infection.
Unexpected findings in our study were that HPV-16 and -18 (E6 and E7) oncoproteins seropositivity and cervical cancer risk were higher in HIV-negative cancer patients compared to HIV-positive cancer patients. Similar findings have been observed in previous studies of HPV antibodies among HIV-positive and HIV-negative patients [35][36][37][38]. In our study, the possible explanation could be that cervical cancer risk in HIV positive women may be related to more than the two main types HPV16 and HPV18 that we tested for so their relative contribution appears attenuated (Additional file 1: Table S3). Other studies on HPV genotype distribution by HIVstatus that were conducted in South Africa reported that HPV types 35, 58 and 33 were more commonly identified [39,40] in addition to types 16 and 18. However, if this is indeed true, then the current HPV vaccine (16/18) may be less effective for HIV positive women. Another possible explanation is that in HIV-positive patients, antibodies could be an indication of reactivation of a latent HPV while in HIV-negative patients antibodies could be reinfection with a new type [35]. Another unexpected finding in our study was that the number of sexual partners was not associated with HPV L1 and (E6 and E7) antibodies. The plausible explanation could be that peak acquisition of HPV occurs at an age earlier than the age of participants in our study [41].
In the new guidelines for screening cervical cancer, the World Health Organization (WHO) has included HPV antibodies and oncoproteins as a future potential screening test [42]. We, therefore, assessed the clinical performance of HPV E6 /E7 antibodies as possible screening tests. We found that the sensitivity of the HPV antibodies for detection of cervical cancer was low but the specificity was high among HIV-positive women. In contrast, sensitivity and specificity for the detection of cervical cancer were high among HIV-negative women. Predictive values would thus vary considerably depending on the background prior likelihood of disease. Our finding is in line with the previous findings on the sensitivity and specificity of HPV16/18 E6/E7 antibodies as a screening test [28,43].
There are strengths to our study. The seroprevalence of HPV antibodies was conducted on a large sample size, and we used a high throughput multiplex serology so all tests were done under the same laboratory conditions. The laboratory tests were conducted 'blind' without prior knowledge of the case/control status of the participants. Nonetheless, our study has some limitations. Serology data on exposure and cervical cancer diagnoses (of all stages) were collected at the same time. We did not have tumor HPV Deoxyribonucleic  Acid (DNA) data to determine the causative type(s). Since only 50-70% of women with detectable HPV DNA in the cervix seroconvert [44], our study underestimates the true prevalence of HPV infection. The JCS did not recruit women with cervical pre-cancerous lesions. Study results are based on black women aged 25-54 years who were recruited from the catchment area of the largest public tertiary hospital in Johannesburg. Therefore, our findings might not be generalizable to other South African regions and women aged 55 years and above.

Conclusions
Our data contribute to the evidence on the importance of HPV-16 and -18 E6 and E7 oncoproteins antibodies in discriminating cervical cancer from controls in a black South African population. HPV L1 antibodies show to be exposure markers of past HPV infection. HPV E6 and E7 show to be important markers of invasive cervical cancer. Furthermore, antibodies to HPV-16 and-18 E6 and E7 seropositivity are strongly associated with cervical cancer in both HIV-positive and HIV-negative patients.