Human Leukocyte Antigen (HLA) Class II -DRB1 and -DQB1 Alleles and the Association with Cervical Cancer in HIV/HPV Co-Infected Women in South Africa

Human Leukocyte Antigen-Allelic Variations May Influence the Age at Cancer Diagnosis in Lynch Syndrome

Lynch syndrome (LS) is an inherited cancer predisposition disorder associated with an elevated risk of developing various solid cancers, but mostly colorectal cancer (CRC). Despite having the same germline pathogenic variant (PV) in one of the mis-match repair genes or the EPCAM gene, Lynch syndrome variant heterozygotes (LSVH) exhibit a remarkable phenotypic variability in the risk of developing cancer. The role of human leukocyte antigen (HLA) in modifying cancer development risk prompted our hypothesis into whether HLA variations act as potential genetic modifiers influencing the age at cancer diagnosis in LSVH. To investigate this, we studied a unique cohort of 426 LSVH carrying the same germline PV in the hMLH1 gene (MLH1:c.1528C > T) in South Africa. We intuitively selected 100 LSVH with the greatest diversity in age at cancer diagnosis (N = 80) and the oldest cancer unaffected LSVH (N = 20) for a high-throughput HLA genotyping of 11 HLA class I and class II loci using the shotgun next-generation sequencing (NGS) technique on the Illumina MiSeq platform. Statistical analyses employed Kaplan-Meier survival analyses with log-rank tests, and Cox proportional hazards using binned HLA data to minimize type I error. Significant associations were observed between young age at cancer diagnosis and HLA-DPB1*04:02 (mean age: 37 y (25-50); hazard ratio (HR) = 3.37; corrected p-value (q) = 0.043) as well as HLA-DPB1 binned alleles (including HLA-DPB1*09:01, HLA-DPB1*10:01, HLA-DPB1*106:01, HLA-DPB1*18:01, HLA-DPB1*20:01, HLA-DPB1*26:01, HLA-DPB1*28:01, HLA-DPB1*296:01, and HLA-DPB1*55:01) (mean age: 37 y (17-63); HR = 2.30, q = 0.045). The involvement of HLA-DPB1 alleles in the age at cancer diagnosis may highlight the potential role of HLA class II in the immune response against cancer development in LSVH. When validated in a larger cohort, these high-risk HLA-DPB1 alleles could be factored into cancer risk prediction models for personalized cancer screening in LSVH.

Network modeling suggests HIV infection phenocopies PI3K-AKT pathway mutations to enhance HPV-associated cervical cancer

Women coinfected with human immunodeficiency virus type 1 (HIV-1) and human papillomavirus (HPV) are six times as likely to develop invasive cervical carcinoma compared to those without HIV. Unlike other HIV-associated cancers, the risk of cervical cancer development does not change when HPV/HIV coinfected women begin antiretroviral therapy, suggesting HIV-associated immune suppression is not a key driver of cervical cancer development in coinfected women. Here, we investigated whether the persistent secretion of inflammatory factors in HIV-positive patients on antiretroviral therapy could enhance cancer signaling in HPV-infected cervical cells via endocrine mechanisms. We integrated previously reported HIV-induced secreted inflammatory factors (Hi-SIFs), HIV and HPV virus-human protein interactions, and cervical cancer patient genomic data using network propagation to understand the pathways underlying disease development in HPV/HIV coinfection. Our results pinpointed the PI3K-AKT signaling pathway to be enriched at the interface between Hi-SIFs and HPV-host molecular networks, in alignment with PI3K pathway mutations being prominent drivers of HPV-associated, but HIV independent, cervical cancer development. Furthermore, we experimentally stimulated cervical cells with 14 Hi-SIFs to assess their ability to activate PI3K-AKT signaling. Strikingly, we found 8 factors (CD14, CXCL11, CXCL9, CXCL13, CXCL17, AHSG, CCL18, and MMP-1) to significantly upregulate AKT phosphorylation (pAKT-S473) relative to a phosphate buffered saline control. Our findings suggest that Hi-SIFs cooperate with HPV infection in cervical cells to over-activate PI3K-AKT signaling, effectively phenocopying PI3K-AKT pathway mutations, resulting in enhanced cervical cancer development in coinfected women. Our insights could support the design of therapeutic interventions targeting the PI3K-AKT pathway or neutralizing Hi-SIFs in HPV/HIV coinfected cervical cancer patients.

High-risk HPV genotypes in Zimbabwean women with cervical cancer: Comparative analyses between HIV-negative and HIV-positive women

BACKGROUND

High-risk human papillomavirus HPV (HR-HPV) modifies cervical cancer risk in people living with HIV, yet African populations are under-represented. We aimed to compare the frequency, multiplicity and consanguinity of HR-HPVs in HIV-negative and HIV-positive Zimbabwean women.

METHODS

This was a cross-sectional study consisting of women with histologically confirmed cervical cancer attending Parirenyatwa Group of Hospitals in Harare, Zimbabwe. Information on HIV status was also collected for comparative analysis. Genomic DNA was extracted from 258 formalin fixed paraffin embedded tumour tissue samples, and analysed for 14 HR-HPV genotypes. Data was analysed using Graphpad Prism and STATA.

RESULTS

Forty-five percent of the cohort was HIV-positive, with a median age of 51 (IQR = 42-62) years. HR-HPV positivity was detected in 96% of biospecimens analysed. HPV16 (48%), was the most prevalent genotype, followed by HPV35 (26%), HPV18 (25%), HPV58 (11%) and HPV33 (10%), irrespective of HIV status. One third of the cohort harboured a single HPV infection, and HPV16 (41%), HPV18 (21%) and HPV35 (21%) were the most prevalent. HIV status did not influence the prevalence and rate of multiple HPV infections (p>0.05). We reported significant (p<0.05) consanguinity of HPV16/18 (OR = 0.3; 95% CI = 0.1-0.9), HPV16/33 (OR = 0.3; 95% CI = 0.1-1.0), HPV16/35 (OR = 3.3; 95% CI = 2.0-6.0), HPV35/51 (OR = 6.0; 95%CI = 1.8-15.0); HPV39/51 (OR = 6.4; 95% CI = 1.8-15), HPV31/52 (OR = 6.2; 95% CI = 1.8-15), HPV39/56 (OR = 11 95% CI = 8-12), HPV59/68 (OR = 8.2; 95% CI = 5.3-12.4), HPV66/68 (OR = 7; 95% CI = 2.4-13.5), independent of age and HIV status.

CONCLUSION

We found that HIV does not influence the frequency, multiplicity and consanguinity of HR-HPV in cervical cancer. For the first time, we report high prevalence of HPV35 among women with confirmed cervical cancer in Zimbabwe, providing additional evidence of HPV diversity in sub-Saharan Africa. The data obtained here probes the need for larger prospective studies to further elucidate HPV diversity and possibility of selective pressure on genotypes.

Association of HLA mismatch and MTOR inhibitor regimens with malignancy and mortality after kidney transplantation

Background The association of mammalian target of rapamycin inhibitors (MTORI) with malignancies and mortality in kidney transplant recipients (KTR) with different degrees of human leukocyte antigen mismatch (HLA-mm) at transplant has not been previously studied. Methods Our observational cohort study included 166, 256 adult KTRs in 2000-2018. Immunosuppression in the first post-transplant year were MTORIs in 13,056 (7.85%) and non-MTORIs in 153,200 (92.15%). We used Cox multivariable regression models to determine the cause-specific hazard ratio (HRcs) of non-melanoma skin cancer (NMSC),solid organ malignancies (SOM)] and all-cause death (death_ac); and the HR of the composite outcomes of NMSC or death_ac and SOM or death_ac associated with MTORI versus non-MTORI regimens in the overall study sample and the 0, 1-3, and 4-6 HLA-A, B and DR mm subgroups. Results NMSC risk was lower with MTORI than non-MTORI in all HLA-mm subgroups [(0 mm, HRcs = 0.67; 95% CI = 0.46-0.97, 1-3 mm, HRcs = 0.73; 95% CI = 0.61-0.87, 4-6 mm, HRcs = 0.69; 95% CI = 0.62-0.76)]. SOM risks were similar between regimens in the 0 HLA mm subgroup (HRcs = 1.10 (95% CI = 0.78-1.57) and lower with MTORI than non-MTORI in the 1-3, and 4-6 HLA-mm subgroups, [(HR = 0.84; (95% CI = 0.71-0.99), and (HR = 0.86; 95% CI = 0.78-0.94); respectively]. Risks of death_ac and composite outcomes (NMSC or death_ac and SOM or death_ac) were higher with MTORI than non-MTORI in almost all HLA-mm subgroups. Conclusion MTORIs are associated with protection from NMSC and SOM in almost all HLA-mm subgroups ca; however, their association with increased all-cause mortality in adult kidney transplant recipients needs further investigation.

The association of human leukocyte antigen alleles with clinical disease progression in HIV-positive cohorts with varied treatment strategies

OBJECTIVES

The Strategic Timing of AntiRetroviral Treatment (START) and Strategies for Management of Antiretroviral Therapy (SMART) trials demonstrated that ART can partly reverse clinically defined immune dysfunction induced by HIV replication. As control of HIV replication is influenced by the HLA region, we explored whether HLA alleles independently influence the risk of clinical events in HIV+ individuals.

DESIGN

Cohort study.

METHODS

In START and SMART participants, associations between imputed HLA alleles and AIDS, infection-related cancer, herpes virus-related AIDS events, chronic inflammation-related conditions, and bacterial pneumonia were assessed. Cox regression was used to estimate hazard ratios for the risk of events among allele carriers versus noncarriers. Models were adjusted for sex, age, geography, race, time-updated CD4+ T-cell counts and HIV viral load and stratified by treatment group within trials. HLA class I and II alleles were analyzed separately. The Benjamini--Hochberg procedure was used to limit the false discovery rate to less than 5% (i.e. q value <0.05).

RESULTS

Among 4829 participants, there were 132 AIDS events, 136 chronic inflammation-related conditions, 167 bacterial pneumonias, 45 infection-related cancers, and 49 herpes virus-related AIDS events. Several associations with q value less than 0.05 were found: HLA-DQB1∗06:04 and HLA-DRB1∗13:02 with AIDS (adjusted HR [95% CI] 2.63 [1.5-4.6] and 2.25 [1.4-3.7], respectively), HLA-B∗15:17 and HLA-DPB1∗15:01 with bacterial pneumonia (4.93 [2.3-10.7] and 4.33 [2.0-9.3], respectively), and HLA-A∗69:01 with infection-related cancer (15.26 [3.5-66.7]). The carriage frequencies of these alleles were 10% or less.

CONCLUSION

This hypothesis-generating study suggests that certain HLA alleles may influence the risk of immune dysfunction-related events irrespective of viral load and CD4+ T-cell count.

New novel non-MHC genes were identified for cervical cancer with an integrative analysis approach of transcriptome-wide association study

Although genome-wide association studies (GWAS) have successfully identified multiple genetic variants associated with cervical cancer, the functional role of those variants is not well understood. To bridge such gap, we integrated the largest cervical cancer GWAS (N = 9,347) with gene expression measured in six human tissues to perform a multi-tissue transcriptome-wide association study (TWAS). We identified a total of 20 associated genes in the European population, especially four novel non-MHC genes (i.e. WDR19, RP11-384K6.2, RP11-384K6.6 and ITSN1). Further, we attempted to validate our results in another independent cervical cancer GWAS from the East Asian population (N = 3,314) and re-discovered four genes including WDR19, HLA-DOB, MICB and OR2B8P. In our subsequent co-expression analysis, we discovered SLAMF7 and LTA were co-expressed in TCGA tumor samples and showed both WDR19 and ITSN1 were enriched in "plasma membrane". Using the protein-protein interaction analysis we observed strong interactions between the proteins produced by genes that are associated with cervical cancer. Overall, our study identified multiple candidate genes, especially four non-MHC genes, which may be causally associated with the risk of cervical cancer. However, further investigations with larger sample size are warranted to validate our findings in diverse populations.

Fine Mapping of the Major Histocompatibility Complex Region and Association of the HLA-B*52:01 Allele With Cervical Cancer in Japanese Women

IMPORTANCE

Understanding the genetic contribution of the major histocompatibility complex (MHC) region to the risk of cervical cancer (CC) will help understand how immune responses to infection with human papillomaviruses are associated with CC.

OBJECTIVE

To determine whether the HLA-B*52:01 allele is associated with CC in Japanese women.

DESIGN, SETTING, AND PARTICIPANTS

This was a multicenter genetic association study. Genotype and phenotype data were obtained from BioBank Japan Project. Additional patients with CC were enrolled from the Aichi Cancer Center Research Institute. An MHC fine-mapping study was conducted on CC risk in the Japanese population by applying a human leukocyte antigen (HLA) imputation method to the large-scale CC genome-wide association study data of using the Japanese population-specific HLA reference panel. Participants included 540 women in BioBank Japan Project with CC or 39 829 women without gynecologic diseases, malignant neoplasms, and MHC-related diseases as controls. An additional 168 patients with CC were recruited from Aichi Cancer Center Research Institute. Histopathological subtypes and clinical stages were not considered. Participants with low genotype call rate, closely related participants, and outliers in the principal component analysis were excluded. Data analysis was performed from August 2018 to January 2020.

MAIN OUTCOMES AND MEASURES

Loci within the MHC region associated with CC risk, and the direction and size of association.

RESULTS

A total of 704 CC cases and 39 556 controls were analyzed. All participants were Japanese women with a median (range) age of 67 (18 to 100) years. One of the class I HLA alleles of HLA-B*52:01 was significantly associated with CC risk (odds ratio, 1.60; 95% CI, 1.38-1.86; P = 7.4 × 10-10). Allele frequency spectra of HLA-B*52:01 are heterogeneous among worldwide populations with high frequency in Japanese populations (0.109 in controls), suggesting its population-specific risk associated with CC. The conditional analysis suggested that HLA-B*52:01 could explain most of the MHC risk associated with CC because no other HLA alleles remained significant after conditioning on the HLA-B*52:01. The HLA amino acid residue-based analysis suggested that HLA-B p.Tyr171His located in the peptide-binding groove was associated with the most significant CC risk (odds ratio, 1.47; 95% CI, 1.30-1.66; P = 1.2 × 10-9).

CONCLUSIONS AND RELEVANCE

The results of this study contribute to understanding of the genetic background of CC. The results suggest that immune responses mediated by class I HLA molecules are associated with susceptibility to CC.

Impacts of CD152 polymorphisms on cervical cancer susceptibility

AIM

The objective of this study was to discuss the effect of CD152 polymorphisms rs231775, rs3087243 and rs5742909 on the susceptibility to cervical cancer.

METHODS

The databases of PubMed, EMBASE, Cochrane Library, ISI Web of Science, Google Scholar Web, CNKI and Wanfang were searched for eligible studies. Chi-square-based Q test examined heterogeneity between included studies, and when P_{heterogeneity} was less than 0.05, random-effect model was used to calculate odds ratios (ORs) with their 95 % confidence intervals (95 % CIs); or else, fixed-effect model was selected. Sensitivity analysis was implemented to determine the stability of final results through removing enrolled studies one at a time and then re-obtaining overall estimates. Publication bias among included studies was checked employing Begg's funnel plot and Egger's test.

RESULTS

CD152 polymorphism rs231775 decreased cervical cancer risk in total analysis under the genetic models of GG vs. AA, GG vs. AA + AG and G vs. A (OR = 0.73, 95 % CI = 0.59-0.91; OR = 0.78, 95 % CI = 0.65-0.94; OR = 0.92, 95 % CI = 0.87-0.98), and so did the polymorphism rs3087243 in total analysis under the comparisons of AA vs. GG, AA + GA vs. GG, AA vs. GG + GA, A vs. G and GA vs. GG (OR = 0.51, 95 % CI = 0.42-0.60; OR = 0.71, 95 % CI = 0.62-0.82; OR = 0.57, 95 % CI = 0.50-0.66; OR = 0.70, 95 % CI = 0.64-0.77; OR = 0.83, 95 % CI = 0.72-0.97). Besides, the polymorphism rs5742909 elevated the disease onset in total analysis under the contrasts of TT vs. CC, TT + CT vs. CC, TT vs. CC + CT, T vs. C and CT vs. CC (OR = 2.66, 95 % CI = 1.75-4.04; OR = 1.54, 95 % CI = 1.24-1.91; OR = 2.13, 95 % CI = 1.12-4.03; OR = 1.44, 95 % CI = 1.17-1.78; OR = 1.49, 95 % CI = 1.22-1.83).

CONCLUSION

CD152 polymorphisms rs231775 and rs3087243 significantly decrease the risk of cervical cancer, while rs5742909 may increase the disease risk.

Personalized Human Papillomavirus Vaccination for Persistence of Immunity for Cervical Cancer Prevention: A Critical Review With Experts' Opinions

The development of cervical cancer has been shown to involve both viral and host factors. The host factors are those that determine the specific response to human papillomavirus (HPV) infection by the patient's immune system. The immune responses to vaccines have been shown to be influenced by polymorphisms in genes involved in innate and adaptive immunity. The specific genetic variants that may influence the immune responses to HPV vaccine which may contribute to persistence of immunity (POI) have not been widely studied yet. In order to address the question as to “is it right to vaccinate all children, and all with equal dose?” we have critically examined the knowledge of common immunogenetic and immunogenomic variations that may influence the HPV vaccine POI across various populations. We have also identified a number of specific research questions that need to be addressed in future research into host molecular genetic variations and HPV vaccine POI in order to afford life-long protection against the development of cervical cancer. This work informs future insights for improved HPV vaccine designs based on common host molecular genetic variations.

Correlation of interleukin-18 gene polymorphism with the susceptibility of condyloma acuminatum in Chinese population

Host immunogenetic setting is involved in the regulation of human papillomavirus (HPV) infection and development of condyloma acuminatum (CA). We investigated the correlation of two common single nucleotide polymorphisms (SNPs) (−607C/A and −137G/C) of IL-18 with the susceptibility of CA in a large Chinese cohort. Out of 408 CA patients analyzed, 300 had HPV infection transmitted through sexual contact (SC) and 108 through non-sexual contact (NSC). In addition, 360 healthy volunteers were enrolled as controls. SNPs at positions −607C/A and −137G/C in IL-18 promoter were analyzed. Comparing CA patients to healthy controls, no dominant relevance was found between the IL-18 promoter −607 C/A or −137G/C polymorphisms and the CA disease either identified genotypically (p > 0.05) or by allelically (p > 0.05). However, the IL-18 promoter −137G/C polymorphism genotype and allele frequencies in the NSC CA group, but not between in the SC group, were significantly higher than in the controls. There was no dominant relevance between IL-18-607C/A polymorphism genotype and allele frequencies among SC, NSC CA patients, and controls. Our study demonstrates that polymorphism −137G/C in IL-18 promoter is significantly correlated with risk of CA in NSC patients.

Investigation of Cervical Tumor Biopsies for Chromosomal Loss of Heterozygosity (LOH) and Microsatellite Instability (MSI) at the HLA II Locus in HIV-1/HPV Co-infected Women

Background: A subgroup of women who are co-infected with human immunodeficiency virus type 1 (HIV-1) and human papillomavirus (HPV) progress rapidly to cervical disease regardless of high CD4 counts. Chromosomal loss of heterozygosity (LOH) and microsatellite instability (MSI) are early frequent genetic alterations occurring in solid tumors. Loss of an allele or part of a chromosome can have multiple functional effects on immune response genes, oncogenes, DNA damage-repair genes, and tumor-suppressor genes. To characterize the genetic alterations that may influence rapid tumor progression in some HIV-1-positive women, the extent of LOH and MSI at the HLA II locus on chromosome 6p in cervical tumor biopsy DNA samples with regard to HIV-1/HPV co-infection in South African women was investigated.

Methods: A total of 164 women with cervical disease were recruited for this study, of which 74 were HIV-1-positive and 90 were HIV-1-seronegative. DNA from cervical tumors and matched buccal swabs were used for analyses. Six fluorescently-labeled oligonucleotide primer pairs in a multiplex PCR amplification were used to study LOH and MSI. Pearson chi-squared test for homogeneity of proportions using an exact p value, a two-proportion Z-score test, ROC curves and a logistic regression model were used for statistical analyses. All p-values were corrected for false discovery rate (FDR) using the Benjamini-Hochberg test and the adjusted p-values (q-values) were reported. All tests were significant when both p and q < 0.05.

Results: Tumor DNA from HIV-1/HPV co-infected women demonstrated a higher frequency of LOH/MSI at the HLA II locus on chromosome 6p21.21 than tumor DNA from HIV-1-seronegative women (D6S2447, 74.2 vs. 42.6%; p = 0.001, q = 0.003), D6S2881 at 6p21.31 (78.3 vs. 42.9%; p = 0.002, q = 0.004), D6S2666 at 6p21.32 (79 vs. 57.1%; p = 0.035, q = 0.052), and D6S2746, at 6p21.33 (64.3 vs. 29.4%; p < 0.001, q < 0.001), respectively.

Conclusions: HPV infection alone can induce LOH/MSI at the HLA II locus in cervical tumor DNA, whereas HIV-1 co-infection exacerbates it, suggesting that this may accelerate cervical disease progression in a subgroup of HIV-1-positive women.