Variation sites at the HLA-G 3’ untranslated region confer differential susceptibility to HIV/HPV co-infection and aneuploidy in cervical cell

Abstract

Post-transcriptional regulatory elements associated with transcript degradation or transcript instability have been described at the 3’ untranslated region (3’UTR) of the HLA-G gene. Considering that HPV infection and aneuploidy, which causes gene instability, are associated with cervical cell malignancy, as well as the fact that HIV infection and HLA-G may modulate the immune response, the present study aimed to compare the frequencies of HLA-G 3′UTR polymorphic sites (14-base pair insertion/deletion, +3142C/G, and +3187A/G) between 226 HIV+ women co-infected (n = 82) or not with HPV (n = 144) and 138 healthy women. We also evaluated the relationship between those HLA-G 3’UTR variants and aneuploidy in cervical cells. HPV types and HLA-G polymorphisms were determined by PCR and sequencing of cervical samples DNA. Aneuploidy in cervical cell was measured by flow cytometry. The HLA-G 3’UTR 14-bp ins/del was not associated with either HIV nor HIV/HPV co-infection. The +3142G allele (p = 0.049) and +3142GG genotype (p = 0.047) were overrepresented in all HIV-infected women. On the other hand, the +3187G allele (p = 0.028) and the +3187GG genotype (p = 0.026) predominated among healthy women. The +3142G (p = 0.023) and +3187A (p = 0.003) alleles were associated with predisposition to HIV infection, irrespective of the presence or not of HIV/HPV co-infection. The diplotype formed by the combination of the +3142CX (CC or CG) and +3187AA genotype conferred the highest risk for aneuploidy in cervical cell induced by HPV. The HLA-G 3’UTR +3142 and +3187 variants conferred distinct susceptibility to HIV infection and aneuploidy.

About the Authors:

Fernanda Silva Medeiros

Roles Conceptualization, Data curation, Formal analysis, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Immunology, Laboratory of Immunogenetics, Aggeu Magalhães Institute, Oswaldo Cruz Foundation, Recife, Brazil

Albert Eduardo Silva Martins

Roles Conceptualization, Data curation, Investigation, Resources, Supervision, Visualization, Writing – original draft

Affiliation: Department of Immunology, Laboratory of Immunogenetics, Aggeu Magalhães Institute, Oswaldo Cruz Foundation, Recife, Brazil

Renan Garcia Gomes

Roles Conceptualization, Data curation, Investigation, Resources, Visualization

Affiliation: Department of Immunology, Laboratory of Immunogenetics, Aggeu Magalhães Institute, Oswaldo Cruz Foundation, Recife, Brazil

Sávio Augusto Vieira de Oliveira

Roles Resources, Visualization

Affiliation: Department of Immunology, Laboratory of Immunogenetics, Aggeu Magalhães Institute, Oswaldo Cruz Foundation, Recife, Brazil

Stefan Welkovic

Roles Investigation, Resources

Affiliation: Integrated Health Centre Amaury de Medeiros (CISAM), University of Pernambuco, Recife, Brazil

Magda Maruza

Roles Investigation, Resources

Affiliation: Hospital Correia Picanço, Health Secretariat of Pernambuco, Recife, Brazil

Maria Luiza Bezerra Menezes

Roles Investigation, Resources

Affiliation: Maternal and Child Department of the Faculty of Medical Sciences—University of Pernambuco (UPE), Recife, Brazil

Ricardo Arraes de Alencar Ximenes

Roles Investigation, Methodology

Affiliation: Department of Tropical Medicine, Federal University of Pernambuco (UFPE), Recife, Brazil

George Tadeu Nunes Diniz

Roles Formal analysis, Methodology, Visualization

Affiliation: Department of Public Health, Laboratory Computational Methods, Aggeu Magalhães Institute, Recife, Brazil

Eduardo Antônio Donadi

Roles Conceptualization, Funding acquisition, Methodology, Resources, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Medicine, School of Medicine of Ribeirão Preto, University of São Paulo (USP), Ribeirão Preto, Brazil

Norma Lucena-Silva

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliations Department of Immunology, Laboratory of Immunogenetics, Aggeu Magalhães Institute, Oswaldo Cruz Foundation, Recife, Brazil, Pediatric Oncology Service, IMIP Hospital, Recife, Brazil

ORCID logo http://orcid.org/0000-0002-6140-1794

Introduction

The human leukocyte antigen-G (HLA-G) is a non-classical major histocompatibility complex molecule, characterized by low protein variability and restricted tissue distribution. HLA-G was initially identified in cytotrophoblast cells and plays an important role in the modulation of the immune response by inhibiting maternal NK and CD8+ T cytotoxic cells, and inducing IL-10 producing T regulatory cells at the maternal-fetal interface, which are crucial events for fetus survival [1,2]. Low level of placental HLA-G has been associated with the development of eclampsia and repeated abortion episodes [3]. The immune checkpoint HLA-G molecule has a well-recognized role on the inhibition of several functions of the innate and adaptive immune responses [4].

HLA-G expression may be differentially modulated by variation sites especially those at the regulatory regions of the HLA-G gene (promoter and 3’ untranslated sequence) [5–7]. Three variation sites at the HLA-G 3’ untranslated region (3’UTR) are important for the post-transcriptional regulation of HLA-G: the 14-base pair (14-bp) insertion/deletion variants influence HLA-G mRNA stability [8,9], the +3142C/G variants regulate HLA-G expression mediated by specific microRNAs [10,11], and the +3187A/G variants affect HLA-G mRNA degradation [12]. Polymorphic sites observed at HLA-G 3’UTR have also been associated with risk to develop cancer [13], autoimmunity disorders [14,15], and infectious diseases [11,16], among other conditions.

HLA-G 3’UTR alleles have been associated with susceptibility to human immunodeficiency virus (HIV) infection in adults [6,17] and in perinatal HIV transmission [18]. HIV+ women have a higher risk of developing human papillomavirus (HPV) co-infection, which is the major cause of human cervical cancer [19,20]. HPV infection [9,11,21–23] and HPV-associated cervical cancer [24–26] have also been associated with the presence of polymorphic sites at the HLA-G gene, irrespective of HIV infection.

In a previous study, we reported a high frequency (47.5%) of HPV infection among HIV+ women from the state of Pernambuco, northeast Brazil, and 59% of those women harbored high-risk HPV types. Despite the high frequency of HPV infection, we observed a low occurrence of cervical lesions in the HIV+ women, which was attributed to the high adherence to antiretroviral therapy [27]. In addition, we observed an association between cellular intraepithelial neoplasia and the presence of aneuploid cells, which causes genomic instability in cervical cells of HIV+/HPV+ women [28].

Considering that, i) the HPV infection is associated with several stages of cervical cell differentiation, ii) the HIV infection and HLA-G may modulate the immune response, iii) the HLA-G 3’UTR segment has several post-transcriptional control elements for HLA-G expression, and iv) aneuploidy is one of the first events in cervical cancer development, the present study aimed to compare the frequencies of HLA-G 3′UTR polymorphic sites (14-bp ins/del, +3142C/G and +3187A/G) between HIV+ women co-infected or not with HPV and healthy women, and to evaluate the association of these HLA-G 3’UTR variants with aneuploidy in cervical cells.

Methodology

Study design and ethical consideration

We performed a case-control study comprising 226 HIV+ women followed-up at the Correia Picanço and the Integrated Health Centre Amaury de Medeiros (CISAM) Hospitals of Recife, capital of the state of Pernambuco, northeast Brazil. Of the 226 HIV+ patients, 82 (36.3%) were co-infected with HPV and 144 (63.7%) were not co-infected with HPV. In parallel, we studied 138 HIV- healthy women followed-up at CISAM for routine examination. The study was approved by the Research Ethics Committee of CISAM under the protocol CAEE: 0011.0.250.000–05. All women agreed to participate in the study and signed a free and informed consent before any assessment.

Data collection

All women answered a questionnaire for clinical and epidemiological assessment, before undergoing colposcopy and cytological examination. Colposcopy was performed according to the criteria of the International Nomenclature of Colposcopy Aspects of Rome [27–28]. Endocervical samples collected using cytobrush were used for cytological analysis, detection of HPV co-infection and immunogenetic studies. Cervical biopsy was performed whenever colposcopic alterations were present and the cervical tissue fragments were preserved in 10% formalin for histopathological analysis. The clinical procedures and the analysis and interpretation of the cervical cancer screening followed the recommendations of the Brazilian guidelines for cervical cancer screening [29]. Details on clinical characterization of HIV+ infected patients have been previously described [27]. Healthy women were not tested for HPV.

Detection of HIV/HPV co-infection

Genomic DNA was extracted from cervical samples using proteinase K (Invitrogen, CA, USA) following the protocol described by Martins et al. 2014 [27]. Subsequently, the DNA samples were amplified by polymerase chain reaction (PCR), using the human glyceraldehyde phosphate dehydrogenase (GAPDH) gene to evaluate sample quality. GAPDH-positive samples were then subjected to another PCR using the degenerate primers, MY09 and MY11, complementary to the L1 region of the HPV genome [30]. PCR products were visualized on a 2% agarose gel stained with 0.5 μg/mL ethidium bromide. The presence of a 450 bp DNA fragment confirmed HIV/HPV co-infection. HPV typing was determined by sequencing the HPV-amplified DNA using the Big Dye Terminator kit and the ABI3100 Genetic Analyzer (Applied Biosystems, Foster City, CA). Sequencing chromatograms were visualized using the MEGA 5.0 software [31] for evaluation of the sequence quality. Analysis of sequence data and sequence similarity searches were performed using the Blast-N tool of the National Center for Biotechnology Information (NCBI) (http://blast.ncbi.nlm.nih.gov/Blast.cgi).

Aneuploidy in cervical cell

Cervical cells were lysed using Pharm Lyse lysing buffer (Becton Dickinson, New Jersey, USA) for 30 minutes prior incubation with RNA solution (100 μg/mL) and propidium iodide. DNA fluorescence was measured by flow cytometry (laser excitation at 488 nm and emission above 600 nm), and DNA index was estimated through by comparison the ratio of the DNA content of cells analyzed with labeled blood diploid cells using the ModFitLT V3.0 software (Verity Software House Inc., Topsham, USA). The presence of two peaks on a histogram with a DNA index greater than 1.16 (hyperploidy) or less than 1.00 (hypoploidy), each with more than 10% of the cell population analyzed in the area corresponding to G0–G1 of the cell cycle in the sample, was considered to be aneuploidy [28].

Polymorphism of the HLA-G 3’ untranslated region

Polymorphic sites between the +2945 and +3259 positions of the HLA-G gene were identified as described by Castelli et al. (2010) [32]. Briefly, amplification was performed in a final volume of 25 μL containing 25 pmol of each primer (HG08F and HG08R), 50 ng of genomic DNA, 1x PCR buffer (20 mM Tris-HCl, 50 mM KCl, 1.5 mM MgCl2), 0.3 mM dNTP, and 1 U Taq DNA polymerase (BioTools, Madrid, Spain). The cycling conditions included one cycle at 94°C for 3 min, followed by 40 cycles at 94°C for 30 s, 60°C for 60 s, 72°C for 30 s, and a final extension step at 72°C for 7 min. A 350 bp amplified fragment was sequenced with the reverse primer, HG08R. All polymorphic sites were identified using SeqMan software version 7.0.0 (DNAStar Inc., Madison, USA) and individually annotated.

Data analyses

The clinical-epidemiological and laboratory data were analyzed using the Excel 2007 software. HPV types were classified according to their phylogenetic characteristics, as previously reported by de Villiers et al. (2004) [33]. Thirty (43.0%) viral isolates were classified as low-risk (HPV 6, 11, 61, 54, 62, 71, 72, 81, 84, 85, 86) and forty (57.0%) were classified as high-risk (HPV 16, 18, 31, 33, 45, 52, 53, 56, 58, 59, 66, 69, 70, 82). We were able to present the frequency of virus types in only 70 patients. The lack of HPV genotype information for 14.6% of the cases was due to the low quality of the sequencing chromatogram, including cases of mixed-HPV subtype infection [27]. However, for the study of the frequencies of HLA-G 3'UTR variants, we were able to evaluate 82 HIV+/HPV+ patients (not all HPV-typed). Data were analyzed using the statistical software Stata 10.0 (Stata Corp LP, USA). Odds ratios (OR) and their respective 95% confidence intervals (CI95%) and p-values (χ2 test and the likelihood ratio) were estimated in the univariate analysis. Associations with p < 0.05 were considered significant. Variables with p < 0.25 were included in the multivariate analysis. The non-parametric Mann-Whitney U test was used when appropriate. The allele and genotype frequencies of the three polymorphic sites 14-bp ins-del, +3142C/G, and +3187A/G were estimated by direct counting using the GENEPOP software version 3.4 [34]. The two-tailed Fisher exact test was used to compare the allele and genotype frequencies between HIV+ and healthy women, and between HIV+ women co-infected with HPV or not. The association between diplotype frequencies and aneuploidy were estimated by logistic regression. Results were considered significant at p < 0.05. Bonferroni correction was performed when necessary. Contingency tables were analyzed using the software GraphPad Prism version 5.01 for Windows (GraphPad Software, San Diego, CA, www.graphpad.com).

Results

Clinical features of the studied patients

In the present study, we evaluated 226 women from a HIV-infected cohort that had been previously studied by our research group and had shown low frequency of high-grade cervical lesion and high adherence to antiretroviral therapy [27], which is provided at no cost to the patient by the Brazilian government. Unfortunately, the research group did not have access to the histological reports of all women that participated in the study because the program for cervical cancer screening subsidized by the Brazilian Ministry of Health determines that histological evaluation must only be performed in women who had cytological abnormalities and atypical colposcopy findings what was not the case for several of the studied patients. Details of demographic and clinical features of HIV+ women co-infected or not with HPV are shown in Table 1. HIV+/HPV+ women presented a 5.5-fold higher risk of altered cervical cytology (p < 0.000) and a 2.0-fold higher risk of chromosomal aberration in cervical cells (p = 0.036) than HIV+/HPV- women.

[Figure omitted. See PDF.]

Table 1. Clinical and epidemiological characteristics of HIV-infected women according the HPV co-infection status attending three reference centers for HIV/AIDS in Recife, northeast of Brazil.

https://doi.org/10.1371/journal.pone.0204679.t001

HLA-G 3’ UTR polymorphisms

First, we compared the frequency of the HLA-G 3’UTR variants between all the HIV+ women (n = 226) and the healthy women (n = 138). The +3142G allele (p = 0.049) and the +3142GG genotype (p = 0.047) were overrepresented in the HIV+ women, while the +3187G allele (p = 0.028) and the +3187GG genotype (p = 0.026) predominated among healthy women. No significant differences were observed for the 14-bp ins/del polymorphism between HIV+ and healthy women (Table 2).

[Figure omitted. See PDF.]

Table 2. Allelict and genotypic frequencies of HLA-G 3’ untranslated region (3’UTR) polymorphic sites of the HLA-G gene in the controls group and in HIV-infected women stratified according to the presence of HIV/HPV co-infection.

https://doi.org/10.1371/journal.pone.0204679.t002

Secondly, we compared the frequency of the same HLA-G 3’UTR variants between HIV+/HPV+ women (n = 82) and healthy controls (n = 138), and no significant results were observed. In contrast, when HIV+/HPV- women (n = 144) were compared to healthy controls, several differences were observed: the +3142G allele (p = 0.023) and the +3142GG genotype (p = 0.038) were associated with susceptibility to HIV infection, and the +3187A alleles (p = 0.003) and +3187AA genotype (p = 0.036) were overrepresented in HIV+/HPV- patients. The power of association between the +3142C/G and +3187A/G polymorphisms and susceptibility to HIV infection increased when the HIV+/HPV+ women were excluded from the HIV+ group. The +3187G allele was overrepresented in HIV+/HPV+ women (p = 0.036). The +3187A/G allelic frequencies of HIV+/HPV- women and healthy women were significantly different (p = 0.018) even after correction for multiple comparison using the Bonferroni test, considering three polymorphic sites and two alleles. No significant differences were observed regarding the 14-bp ins/del polymorphism between HIV+/HPV- and control women. The HLA-G +3142GG/+3187AA diplotype showed a trend increase (p = 0.071) in the group of HIV+ women compared to healthy women (Table 3).

[Figure omitted. See PDF.]

Table 3. Diplotypes frequencies of HLA-G 3’ untranslated region (3’UTR) polymorphic sites of the HLA-G gene in the controls group and in HIV-infected women.

https://doi.org/10.1371/journal.pone.0204679.t003

Association of HLA-G 3′ UTR polymorphisms with aneuploidy in cervical cell

We recently reported a lack of association between aneuploidy in cervical cells and HIV+/HPV+ co-infection (p = 0.22), low-risk (p = 0.43) or high-risk HPV (p = 0.17), and cytological alteration (p = 0.35) in the HIV-infected women cohort [28]. However, aneuploidy was associated with histological lesion (p = 0.03) [28]. Considering that aneuploidy is a phenomenon associated with cervical cell transformation, and HLA-G is an immune checkpoint protein with inhibitory effect on several immune cells, we hypothesized that women harboring HLA-G 3’UTR alleles associated with high HLA-G expression may also have higher risk to develop aneuploidy.

Thirty four samples were excluded from HIV+ group (n = 226) because they either had no data on cellular ploidy (n = 24) or lacked polymorphism information (n = 10). Comparing the diplotype distribution of women with cervical cells presenting or not aneuploidy, we observed that the +3142GG/+3187AA diplotype was the most prevalent in both aneuploid (45.45%) and diploid (43.92%) cervical cells. Considering the double heterozygous +3142CG/+3187AG diplotype as reference and using the logistic regression model, we observed that: (i) diplotypes defined by +3142CC/+3187GX, where X is either the G or A allele, showed no differences related to ploidy (Table 4), (ii) the double heterozygous diplotype (+3142CG/+3187AG) was the second most prevalent in the diploid group (29.05%) and was underrepresented in the aneuploidy group (15.91%), and (iii) the frequency of the +3142CX/+3187AA diplotype was significantly higher than the frequency of the double heterozygous haplotypes in the aneuploid group, and this association was strengthened when HPV infection was considered in the model. Unfortunately, we were not able to investigate the influence of HPV-risk classification on the relationship between +3142CX/+3187AA diplotype and aneuploidy using the logistic regression model due to a small sample size after data stratification by HPV risk.

[Figure omitted. See PDF.]

Table 4. Linear regression between the HLA-G +3142/+3187 diplotypes and presence of aneuploidy in cervical cells from HIV-infected women.

https://doi.org/10.1371/journal.pone.0204679.t004

We also evaluated whether oncotic cytology was associated with HLA-G polymorphisms in HIV+ women (Table 5). No significant differences were observed between the +3142/+3187 diplotypes and the presence of cytological alterations (p = 0.407), even after data stratification by presence (p = 0.324) or absence (p = 0.639) of HPV infection. No association (p = 0.603) was observed between the +3142/+3187 diplotypes and CD4+ T lymphocyte count (Table 6).

[Figure omitted. See PDF.]

Table 5. HLA-G +3142/+3187 diplotypes, cytological abnormalities and stratified case analysis according to the risk of HPV infection in HIV-infected women.

https://doi.org/10.1371/journal.pone.0204679.t005

[Figure omitted. See PDF.]

Table 6. Frequency distribution of the HLA-G +3142/+3187 diplotypes in relation to the CD4+ T lymphocyte count values in HIV-infected women.

https://doi.org/10.1371/journal.pone.0204679.t006

Discussion

The role of HLA-G in HIV infection is controversial due to several reasons, including disease diversity, presence of underlying disorders, patient ethnicity and treatment regimens. Although there are eight well-known polymorphic sites in the HLA-G 3’UTR region, in our study, the frequencies of only the +3142C/G and +3187A/G sites were statistically different between healthy controls and HIV+ patients. We believe that by including a screening step for selecting only HLA-G polymorphic sites that contribute to susceptibility to viral infection, we minimize the inclusion of genetic factors not related to it. These two polymorphic sites are target for microRNA, important post-transcriptional gene expression regulators, which expression profiles have been demonstrated to diverge according to cell type and stimulus [35]. Nevertheless, we also discussed the 14-bp ins/del allele and genotype frequencies, because of its influence on HLA-G levels in HIV+ patients, as has been demonstrated in previous studies [6, 8, 17].

We also understand that environmental pressures may affect the evolution of the immune system differentially and therefore the polymorphic sites we studied may not be representative for other populations from different locations. In a previous study, we evaluated HLA-G 3’UTR polymorphic sites allele, genotype and haplotype frequencies in healthy blood donors from a northeastern and southeastern Brazilian regions, and we found that the distribution of the HLA-G 3'UTR 14bp-ins/del, +3142C/G, and +3187A/G were similar [36]. However, comparing these published data with the distribution frequency of these polymorphic sites observed in healthy women group reported in the present study, no differences were found for the alleles and genotypes regarding to the northeastern population; in contrast, we found a higher frequency of +3187GG genotype (14.4% vs. 4.5%, respectively) in the healthy women group of the present study than in healthy blood donors from southeast Brazil [32]. This findings could be related to the contribution of different ethnic groups to the constitution of the Brazilian populations [36]. It is noteworthy that a study on the Brazilian ancestry based on mitochondrial DNA showed a major contribution of African ancestry in northeast Brazil [37].

There are few reports in the literature regarding the association between presence of HLA-G 3’UTR polymorphisms to susceptibility to HIV and HPV infection, and the 14-bp ins/del site is by far the most studied [8, 17, 18]. The frequency of +3187A allele in homozygosis was found in only one study on HIV vertical transmission, and it was not appropriated for comparison [38]. Future studies may confirm whether the presence of the +3142G and +3187A alleles are associated with risk of HIV infection in different populations.

Nevertheless, HIV infection increases the production of HLA-G by naïve T CD8+ cells and increases effector and memory cell lines [39]. The virus itself induces high levels of sHLA-G [40] and the antiretroviral therapy, during the immunological reconstitution, down-regulates HLA-G [41]. The HLA-G 14-bp ins allele, which is associated with low levels of sHLA-G in healthy Caucasian individuals [7], has been associated with low levels of sHLA-G in African HIV-infected individuals [6]. In contrast, the presence of the 14-bp del/del genotype, which has been associated with high sHLA-G levels in healthy Caucasians [7], was associated with high expression of HLA-G mRNA, high viral load, lower CD4+ T lymphocytes counts and low survival rate in HIV-infected Zimbabwean women [8]. The HLA-G 3’UTR 14-bp ins allele, the ins/ins genotype, the 14-bp insertion/+3142G (insG) haplotype and the insG/insG diplotype were overrepresented in African-derived patients HIV+ from Southern Brazil, irrespective of the underlying disorders; but, the same was not observed in European-derived Brazilians [17]. In the present study, the associations with the 14-bp polymorphism was not observed in the HIV+ women co-infected or not with HPV from the state of Pernambuco, northeast Brazil, which have a higher African-ancestry than Southern Brazilians.

The role of HLA-G on HPV mono-infection has not been completely elucidated yet. Gimenes et al. (2014) [42] recently reviewed studies correlating the association between HLA-G polymorphisms and cervical carcinogenesis. The presence of high levels of HLA-G mRNA in cervical lesion was reported to contribute to their progression in HPV 16/18 infected patients [21,43,44]. The HLA-G 3’UTR diversity has only been evaluated in patients with HPV mono-infection, including the 14-bpdel/+3142C haplotype in Italian women [9], and the +3142C allele in Taiwanese women [45], which were associated with increased risk for developing HPV cervical cancer. For the HIV/HPV co-infection, we observed that the +3142G and +3187A alleles were associated with predisposition to HIV infection independent of the presence of HIV/HPV co-infection. In addition, the double doses, +3142GG and +3187AA, also exhibited a trend increased risk for HIV infection, independent of the HIV/HPV co-infection. Considering that the +3142G allele reduces HLA-G mRNA expression by the action of specific microRNAs [10], and that the +3187A allele decreases the HLA-G mRNA stability due to the proximity to an AU motif [12], our results indicate that susceptibility to HIV, irrespective of HPV infection, may be favored by low HLA-G expression.

Since HIV and HPV present distinct biological behavior, we further analyzed the relationship between the HLA-G 3’UTR polymorphic sites and aneuploidy, which may represent the initial step for progression towards cervical malignancy [46]. In the present study, the +3142CX/+3187AA diplotype was associated with aneuploidy of cervical cells; i.e., the presence of +3142CC and +3142CG increased by 7.40-fold and 3.21-fold the risk for cervical lesion, respectively. Noteworthy, the association with +3142C/G polymorphic site was only observed when the +3187AA genotype was present as well.

Although the +3142C allele was underrepresented in the HIV+ patients compared to healthy women, this allele was associated with aneuploidy in cervical cell in the HIV+ women. Studies have shown that the +3142C allele and +3142CC genotype increase sHLA-G production, which has been associated with HPV infection and risk to develop cervical cancer [9,45]. In this context, we showed that the presence of HPV infection increased the contribution of the +3142C allele to the susceptibility of the +3142CX/+3187AA diplotype to aneuploidy by 20.5% (in homozygosis) and 4.6% (in heterozygosis) (Table 4); however, it was not associated with cytological abnormalities. We have previously shown the relationship between the presence of cervical tissue lesions and the DNA index, but cytological alteration was not predictive of progression of cervical lesion [28].

In resume, we showed that the +3142G and +3187A alleles, which are related to low HLA-G expression [10,12], were associated with predisposition to HIV infection independent of the presence of HIV/HPV co-infection, but HIV induces the production of HLA-G in untreated patients through unclear mechanisms [41]. Perhaps, the increased HLA-G level dependent of HIV infection is one of the mechanisms behind the association between HIV infection and high-risk of HPV co-infection [47,48], and increased risk of cytological abnormalities and cervical cancer [49,50]. Indeed, we showed that in HIV+ women the +3142C allele and the +3142CC genotype, which are related to increased sHLA-G production, were associated to the presence of aneuploidy in cervical cells in an allele dose-dependent effect. Our results corroborate previous studies that showed that the +3142C allele is associated to HPV infection and to a higher risk to develop cervical cancer in HIV- women [9,45]. Taken together, we conclude that aneuploidy, one of the first steps towards the development of cervical malignancies [46], is favored by HLA-G alleles related to high HLA-G production, because the inhibitory effect of the HLA-G on immune cells modulates the immune cellular response, which then contribute to replication and persistence of HPV infections, and increases the chance of cervical lesion progression.

Our study dealt with limitations: we did not quantify the HLA-G levels, as well as, the fact that histological reports were not available for all women that participated in the study because the Brazilian Ministry of Health’s protocol for screening cervical cancer, defines that histological evaluation must only be performed in women who had cytological abnormalities and atypical colposcopy findings. In addition, unfortunately, some women did not return to the ambulatory clinic for the histological exams or reassessment, and, consequently, their exam results could not be provided to the research group, including seven patients with altered cytology for whom the histology results were not recorded.

Acknowledgments

We thank Viviane Carvalho for invaluable technical assistance and the Program for Technological Development in Tools for Health (PDTIS-FIOCRUZ).

Citation: Medeiros FS, Martins AES, Gomes RG, de Oliveira SAV, Welkovic S, Maruza M, et al. (2018) Variation sites at the HLA-G 3’ untranslated region confer differential susceptibility to HIV/HPV co-infection and aneuploidy in cervical cell. PLoS ONE 13(10): e0204679. https://doi.org/10.1371/journal.pone.0204679

References

1. Carosella ED, Favier B, Rouas-Freiss N, Moreau P, Lemaoult J. Beyond the increasing complexity of the immunomodulatory HLA-G molecule. Blood. 2008;111: 4862–4870. pmid:18334671

2. Donadi EA, Castelli EC, Arnaiz-Villena A, Roger M, Rey D, Moreau P. Implications of the polymorphism of HLA-G on its function, regulation, evolution and disease association. Cell Mol Life Sci. 2011;68: 369–395. pmid:21107637

3. Hviid T, Rizzo R, Christiansen O, Melchiorri L, Lindhard A, Baricordi O. HLA-G and IL-10 in serum in relation to HLA-G genotype and polymorphisms. Immunogenetics. 2004;56: 212–217. pmid:15133645

4. Carosella ED. The tolerogenic molecule HLA-G. Immunol Lett. 2011;138: 22–24. pmid:21333680

5. Baricordi O, Stignani M, Melchiorri L, Rizzo R. HLA-G and Inflammatory Diseases. Inflamm Allergy-Drug Targets. 2008;7: 67–74. pmid:18691135

6. Lajoie J, Massinga Loembe M, Poudrier J, Guédou F, Pépin J, Labbé A-C, et al. Blood soluble human leukocyte antigen G levels are associated with human immunodeficiency virus type 1 infection in Beninese commercial sex workers. Hum Immunol. 2010;71: 182–185. pmid:19913587

7. Martelli-Palomino G, Pancotto JA, Muniz YC, Mendes-Junior CT, Castelli EC, Massaro JD, et al. Polymorphic Sites at the 3’ Untranslated Region of the HLA-G Gene Are Associated with Differential hla-g Soluble Levels in the Brazilian and French Population. PLoS One. 2013;8: e71742. pmid:24204558

8. Larsen MH, Zinyama R, Kallestrup P, Gerstoft J, Gomo E, Thørner LW, et al. HLA-G 3′ Untranslated Region 14–Base Pair Deletion: Association With Poor Survival in an HIV-1–Infected Zimbabwean Population. J Infect Dis. 2013;207: 903–906. pmid:23264673

9. Bortolotti D, Gentili V, Rotola A, Di Luca D, Rizzo R. Implication of HLA-G 3’ untranslated region polymorphisms in human papillomavirus infection. Tissue Antigens. 2014;83: 113–118. pmid:24397524

10. Tan Z, Randall G, Fan J, Camoretti-Mercado B, Brockman-Schneider R, Pan L, et al. Allele-Specific Targeting of microRNAs to HLA-G and Risk of Asthma. Am J Hum Genet. 2007;81: 829–834. pmid:17847008

11. Xu H-H, Shi W-W, Lin A, Yan W-H. HLA-G 3′ untranslated region polymorphisms influence the susceptibility for human papillomavirus infection. Tissue Antigens. 2014;84: 216–222. pmid:24758208

12. Yie S -m., Li L -h., Xiao R, Librach CL. A single base-pair mutation in the 3’-untranslated region of HLA-G mRNA is associated with pre-eclampsia. Mol Hum Reprod. 2008;14: 649–653. pmid:18952696

13. Li T, Huang H, Liao D, Ling H, Su B, Cai M. Genetic polymorphism in HLA-G 3′UTR 14-bp ins/del and risk of cancer: a meta-analysis of case–control study. Mol Genet Genomics. 2015;290: 1235–1245. pmid:25573780

14. Lucena-Silva N, de Souza VSB, Gomes RG, Fantinatti A, Muniz YCN, de Albuquerque RS, et al. HLA-G 3’ Untranslated Region Polymorphisms Are Associated with Systemic Lupus Erythematosus in 2 Brazilian Populations. J Rheumatol. 2013;40: 1104–1113. pmid:23678155

15. de Albuquerque RS, Mendes-Junior CT, Lucena-Silva N, da Silva CLL, Rassi DM, Veiga-Castelli LC, et al. Association of HLA-G 3’ untranslated region variants with type 1 diabetes mellitus. Hum Immunol. 2016;77: 358–364. pmid:26883941

16. Lucena-Silva N, Teixeira MAG, De A, Ramos L, De Albuquerque RS, Diniz GTN, et al. The +3187A/G HLA-G polymorphic site is associated with polar forms and reactive reaction in leprosy. Mol Genet Genomic Med. 2013;1: 123–130. pmid:24498610

17. Da Silva GK, Vianna P, Veit TD, Crovella S, Catamo E, Aparicio Cordero EA, et al. Influence of HLA-G polymorphisms in human immunodeficiency virus infection and hepatitis C virus co-infection in Brazilian and Italian individuals. Infect Genet Evol. 2014;21: 418–423. pmid:24389119

18. Segat L, Zupin L, Kim H-Y, Catamo E, Thea DM, Kankasa C, et al. HLA-G 14 bp deletion/insertion polymorphism and mother-to-child transmission of HIV. Tissue Antigens. 2014;83: 161–167. pmid:24571474

19. Naucler P, Mabota da Costa F, Leopoldo da Costa J, Ljungberg O, Bugalho A, Dillner J, et al. Human papillomavirus type-specific risk of cervical cancer in a population with high human immunodeficiency virus prevalence: case–control study. J Gen Virol. 2011;92: 2784–2791. pmid:21900420

20. Crosbie EJ, Einstein MH, Franceschi S, Kitchener HC. Human papillomavirus and cervical cancer. Lancet. 2013;382: 889–899. pmid:23618600

21. Dong D, Hong Yang, Ke Li, Gang Xu, Song L, Fan X, et al. Human Leukocyte Antigen-G (HLA-G) Expression in Cervical Lesions: Association With Cancer Progression, HPV 16/18 Infection, and Host Immune Response. Reprod Sci. 2010;17: 718–723. pmid:20445010

22. Ferguson R, Ramanakumar A V., Richardson H, Tellier PP, Coutlée F, Franco EL, et al. Human leukocyte antigen (HLA)-E and HLA-G polymorphisms in human papillomavirus infection susceptibility and persistence. Hum Immunol. 2011;72: 337–341. pmid:21256910

23. Metcalfe S, Roger M, Faucher M-C, Coutlée F, Franco EL, Brassard P. The association between human leukocyte antigen (HLA)-G polymorphisms and human papillomavirus (HPV) infection in Inuit women of northern Quebec. Hum Immunol. 2013;74: 1610–1615. pmid:23994586

24. Simões RT, Gonçalves MAG, Castelli EC, Júnior CM, Bettini JSR, L Discorde M, et al. HLA-G polymorphisms in women with squamous intraepithelial lesions harboring human papillomavirus. Mod Pathol. 2009;22: 1075–1082. pmid:19407850

25. Ferguson R, Ramanakumar A V, Koushik A, Coutl Ee F, Franco E, Roger M. Human leukocyte antigen G polymorphism is associated with an increased risk of invasive cancer of the uterine cervix. Int J cancer. 2012;131: 312–319. pmid:22095460

26. Silva ID, Muniz YCN, Sousa MCPS, Silva KR, Castelli EC, Filho JCG, et al. HLA-G 3′UTR polymorphisms in high grade and invasive cervico-vaginal cancer. Hum Immunol. 2013;74: 452–458. pmid:23228396

27. Martins AES, Lucena-Silva N, Garcia RG, Welkovic S, Barboza A, Menezes MLB, et al. Prevalence of human papillomavirus infection, distribution of viral types and risk factors in cervical samples from human immunodeficiency virus-positive women attending three human immunodeficiency virus-acquired immune deficiency syndrome reference cent. Mem Inst Oswaldo Cruz. 2014;109: 738–747. pmid:25317701

28. Martins AES, Lucena-Silva N, Garcia RG, Welkovic S, Barbosa A, Menezes MLB, et al. Prognostic Evaluation of DNA Index in HIV-HPV Co-Infected Women Cervical Samples Attending in Reference Centers for HIV-AIDS in Recife. Tang J, editor. PLoS One. 2014;9: e104801. pmid:25144309

29. Schneider ML; Schneider V. Atlas de diagnóstico diferencial em citologia ginecológica. Germany: Revinter; 1998; pp 165.

30. Manos M, Ting Y, Wright D, Lewis A, Broker T, Wolinsky S. Use of polymerase chain reaction amplification for the detection of genital human papillomaviruses. Cancer Cells 7: Molecular Diagnostics of Human Cancer. 1989. pp. 209–214.

31. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011;28: 2731–2739. pmid:21546353

32. Castelli E, Mendes-Junior C, Deghaide N, de Albuquerque R, Muniz Y, Simões R, et al. The genetic structure of 3'untranslated region of the HLA-G gene: polymorphisms and haplotypes. Genes Immun. 2010;11: 134–141. pmid:19798077

33. de Villiers E-M, Fauquet C, Broker TR, Bernard H-U, zur Hausen H. Classification of papillomaviruses. Virology. 2004;324: 17–27. pmid:15183049

34. Rousset F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8: 103–106. pmid:21585727

35. Porto IO, Mendes-Junior CT, Felício LP, Georg RC, Moreau P, Donadi EA, et al. MicroRNAs targeting the immunomodulatory HLA-G gene: a new survey searching for microRNAs with potential to regulate HLA-G. Mol Immunol. 2015;65: 230–241. https://doi.org/10.1016/j.molimm.2015.01.030 pmid:25700346

36. Lucena‐Silva N, Monteiro AR, de Albuquerque RS, Gomes RG, Mendes‐Junior CT, Castelli EC, et al. Haplotype frequencies based on eight polymorphic sites at the 3′ untranslated region of the HLA‐G gene in individuals from two different geographical regions of Brazil. Tissue antigens. 2012;79: 272–278. pmid:22283419

37. Alves-Silva J, da Silva Santos M, Guimarães PE, Ferreira AC, Bandelt HJ, Pena SD, et al. The ancestry of Brazilian mtDNA lineages. Am J Hum Genet. 2000;67:444–61. pmid:10873790

38. Hong HA, Paximadis M, Gray GE, Kuhn L, Tiemessen CT. (2015). Maternal human leukocyte antigen-G (HLA-G) genetic variants associate with in utero mother-to-child transmission of HIV-1 in Black South Africans. Infect Genet Evol. 2015;30: 147–158. pmid:25541520

39. Lozano JM, González R, Luque J, Frias M, Rivero A, Peña J. CD8 + HLA-G + Regulatory T Cells Are Expanded in HIV-1-Infected Patients. Viral Immunol. 2009;22: 463–465. pmid:19951184

40. Huang J, Burke P, Yang Y, Seiss K, Beamon J, Cung T, et al. Soluble HLA-G Inhibits Myeloid Dendritic Cell Function in HIV-1 Infection by Interacting with Leukocyte Immunoglobulin-Like Receptor B2. J Virol. 2010;84: 10784–10791. pmid:20702625

41. Murdaca G, Contini P, Setti M, Cagnati P, Lantieri F, Indiveri F, et al. Behavior of non-classical soluble HLA class G antigens in human immunodeficiency virus 1-infected patients before and after HAART: Comparison with classical soluble HLA-A, -B, -C antigens and potential role in immune-reconstitution. Clin Immunol. Elsevier Inc.; 2009;133: 238–244. pmid:19762282

42. Gimenes F, Teixeira JJV, de Abreu ALP, Souza RP, Pereira MW, da Silva VRS, et al. Human leukocyte antigen (HLA)-G and cervical cancer immunoediting: A candidate molecule for therapeutic intervention and prognostic biomarker? Biochim Biophys Acta. 2014; 1846:576–89. pmid:25453366

43. Yoon BS, Kim YT, Kim JW, Kim SH, Kim JH, Kim SW. Expression of human leukocyte antigen-G and its correlation with interleukin-10 expression in cervical carcinoma. Int J Gynecol Obstet. 2007;98: 48–53. pmid:17490670

44. Rodríguez JA, Galeano L, Palacios DM, Gómez C, Serrano ML, Bravo MM, et al. Altered HLA Class I and HLA-G Expression Is Associated with IL-10 Expression in Patients with Cervical Cancer. Pathobiology. 2012;79: 72–83. pmid:22213066

45. Yang Y-C, Chang T-Y, Chen T-C, Lin W-S, Chang S-C, Lee Y-J. Human leucocyte antigen-G polymorphisms are associated with cervical squamous cell carcinoma risk in Taiwanese women. Eur J Cancer. 2014;50: 469–474. pmid:24262586

46. Melsheimer P. DNA Aneuploidy and Integration of Human Papillomavirus Type 16 E6/E7 Oncogenes in Intraepithelial Neoplasia and Invasive Squamous Cell Carcinoma of the Cervix Uteri. Clin Cancer Res. 2004;10: 3059–3063. pmid:15131043

47. Garbuglia AR, Piselli P, Lapa D, Sias C, Del Nonno F, Baiocchini A, et al. Frequency and multiplicity of human papillomavirus infection in HIV-1 positive women in Italy. J Clin Virol. 2012;54: 141–146. pmid:22437054

48. Rocha-Brischiliari SC, Gimenes F, De Abreu ALP, Irie MMT, Souza RP, Santana RG, et al. Risk factors for cervical HPV infection and genotypes distribution in HIV-infected South Brazilian women. Infect Agent Cancer. 2014;9: 1–6.

49. Ananworanich J, Prasitsuebsai W, Kerr SJ, Hansudewechakul R, Teeratakulpisarn N, Saisawat K, et al. Cervical cytological abnormalities and HPV infection in perinatally HIV-infected adolescents. J Virus Erad. 2015;1: 30–37. pmid:26005716

50. Abraham AG, Strickler HD, Jing Y, Gange SJ, Sterling TR, Silverberg M, et al. Invasive cervical cancer risk among HIV-infected women: A North American multi-cohort collaboration prospective study. J Acquir Immune Defic Syndr. 2013;62: 405–413. pmid:23254153