The genetic variant rs55986091 HLA-DQB1 is associated with a protective effect against cervical cancer

Introduction: Cervical cancer (CC) is a prevalent malignancy affecting women globally. The primary causative factor of CC is the high-risk oncogenic human papillomavirus (HR-HPV). However, it is noteworthy that not all women infected with HR-HPV develop cancer, indicating the potential involvement of genetic predisposition in the development of CC. This study aims to identify genetic risks and their distribution in groups of women with different epidemiological features of HR-HPV.

Materials and methods: A comparison was conducted among four groups of women, comprising 218 HPV-negative women, 120 HPV-positive women, 191 women diagnosed with cervical intraepithelial neoplasia (CIN) grade 2 or 3, and 124 women diagnosed with CC. The analysis focused on four single nucleotide polymorphisms (SNPs): rs55986091 in HLA-DQB1, rs138446575 in TTC34, rs1048943 in CYP1A1, and rs2910164 in miRNA-146a.

Results: The rs55986091-A allele exhibited a protective effect within the “CC” group when compared to the “HPV-Negative” group (OR = 0.4, 95% CI= 0.25-0.65) using a log-additive model. Additionally, similar protective effects were observed in the “CIN 2/3” group compared to the “HPV-Negative” group (OR = 0.47, 95% CI = 0.28-0.79).

Conclusion: The data obtained emphasize the importance of developing PCR-based diagnostic kits for the identification of SNP alleles, particularly for rs55986091, among HR-HPV-positive women within the Russian population.

1 Background

Cervical cancer (CC) is globally recognized as the fourth most common and fatal oncological disease among women. In 2020, approximately 604,000 new cases of CC were reported, resulting in over 340,000 fatalities and contributing to 7.7% of all cancer-related deaths (1). The incidence of CC in Russia has witnessed a concerning rise over the past decade (2009 to 2019), with a nearly 12-13% increase (from 14,000 to 17,000 cases and from 112 to 127 cases per 100,000 population). This escalating trend highlights the societal significance of this disease. Notably, women aged 30 to 44 years who are in their socially active and reproductive years account for 32.4% of CC cases (2).

According to Russian cancer statistics from 2019, the existing preventive measures for cervical precancerous conditions have demonstrated limited effectiveness, particularly concerning vaccination programs. Furthermore, there is a lack of awareness among the population regarding the significance and accessibility of CC screening and a lack of motivation for undergoing preventive examinations. Consequently, the mortality rate for CC remains largely unchanged (2).

Infection with high-risk oncogenic human papillomavirus (HR-HPV) has been established as a confirmed carcinogenic factor for CC. Within this study, 14 HPV types were identified as highly carcinogenic, namely, types 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68 (3). The association between HR-HPV and CC has been found to be stronger than the association between smoking and lung cancer (4). Interestingly, among women who test positive for HPV, the incidence of CC is only 0.015%, suggesting that genetic susceptibility may play a role (5).

It is worth noting that HPV has the potential to induce the development of various tumors, not only in women but also in men. These include anogenital lesions, head and neck tumors, and other types of tumors. As a result, the social and economic burden associated with HPV-associated diseases may have been significantly underestimated (6).

The identification of single nucleotide polymorphisms (SNP) alleles associated with diseases enables the detection of potential inherited predispositions to the development of pathological conditions during the asymptomatic phase. This aids in the timely implementation of diagnostic and treatment measures (7).

Given the lengthy asymptomatic period, sexual transmission, and the impact on women of reproductive age, assessing the genetic risk of developing CC holds significant clinical importance. This is particularly relevant for high-risk groups, including women infected with HIV (8).

In our previous study, we conducted an analysis of SNPs associated with CC risk derived from genome-wide association studies (GWAS) and case-control studies, including a meta-analysis (9). For this particular publication, we focused on examining four specific SNPs: rs55986091 in HLA-DQB1, rs138446575 in TTC34, rs1048943 in CYP1A1, and rs2910164 in miRNA-146a.

Previously, the rs55986091 SNP (with the protective allele A) was examined in a single GWAS (10). The study revealed a significant association, with an OR of 0.66, a 95% CI of 0.60-0.72, and a p-value of less than 1×10^-11, specifically for the group with cervical intraepithelial neoplasia grade 3 (CIN3) or invasive CC (10). Many reviews acknowledge the significance of rs55986091 as the SNP most strongly associated with CC, referring to GWAS findings (11, 12).

The rs138446575 SNP (with the risk allele T) was also examined in the same GWAS study (10). The findings revealed a significant association, with an OR of 2.39 and a 95% CI of 1.75-3.27, at a p-value of 4.97×10⁻⁸, specifically within the groups diagnosed with invasive CC. The rs138446575 SNP is located between the TTC34 and PRDM16 genes in a region whose functions have not yet been fully identified.

The rs1048943 SNP (with the risk allele C) has been investigated in meta-analyses, revealing significant associations with CC. Among women from India, the meta-analysis showed an OR of 2.34, with a 95% CI of 1.37-3.99 (13). For other populations, the analysis demonstrated an OR of 2, with a 95% CI of 1.33-3 (14). In the analysis conducted on the European population (white population of Europe and the United States), a positive trend in CC development was observed, with an OR of 2.22 and a 95% CI of 1.48-3.33 (9). The rs1048943 SNP causes an amino acid substitution in the CYP1A1 gene, resulting in reduced enzymatic activity and being associated with the development of other types of cancers (14, 15).

The rs2910164 (with the risk allele C) has been examined in previous articles (16, 17), and statistically significant associations have been identified. In one study, the analysis revealed an OR of 1.35, with a 95% CI of 1.06-1.74 (16). Another study reported an OR of 1.41, with a 95% CI of 1.17-1.70 (17). These significant associations were observed specifically within the Asian population. Additionally, the rs2910164 has been included in the assessment of risks for the European population.

The primary objective of this study was to determine the genetic risks and their distribution among different groups of women with varying HR-HPV clinical presentations.

2 Materials and methods

2.1 Oligonucleotide design

In this study, a real-time PCR-based approach was performed for the discrimination of SNP alleles (18). The design of the probes used in this approach was based on the methodology employed in the study (19). Primers were selected for pyrosequencing of all SNPs, and an additional Sanger sequencing was chosen specifically for rs55986091. The specific sequences of the oligonucleotides used in the study can be found in Table 1.

TABLE 1

Table 1 Oligonucleotide sequences.

2.2 DNA extraction, PCR, and sequencing

DNA extraction from whole blood samples and cervical scrapings (exocervix and endocervix) collected in the BD SurePath transport medium (BD Diagnostics, USA) was performed using the RIBO-prep (Registration Certificate FSR 2008/03147 of 15/9/2008) and AmpliSens^® DNA-sorb-D (Registration Certificate RZN 2015/3503 of 27/3/2019) kits. The investigation of HPV types was conducted using the AmpliSens^® HR HPV screen-titre-14-FL (Registration Certificate RZN 2017/5387 of 22/2/2019) and AmpliSens^® HR HPV geno-type-titre-FL (Registration Certificate RZN 2017/6533 of 27/3/2019) kits. These kits enabled HR-HPV genotyping for 14 specific types, namely, 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68.

The PCR reactions were performed using a reaction mix of 25 μl, which included primers, probes, and dNTP (0.44 mM) in a volume of 10 μl. Additionally, 0.5 μl of TaqF polymerase and 4.5 μl of RT-PCR-mix-2 FEP/FRT reagents were added. The extracted DNA (10 μl) was also included in the reaction mix. Real-time PCR was carried out using the Rotor-Gene Q cycler (Qiagen, Germany) and DTprime (DNA-Technology, Russia) instruments. The following cycling conditions were employed: an initial denaturation step at 95°C for 15 minutes (1 cycle), followed by 45 cycles of denaturation at 95°C for 5 seconds, annealing at 60°C for 20 seconds, and extension at 72°C for 10 seconds. During the amplification process, a fluorescent signal was detected at 60°C. The amplification results were analyzed by setting the threshold at 10% of the highest value of the fluorescent signal for each channel. All the reagents and kits used for PCR and DNA extraction were manufactured by AmpliSens, Russia, and provided by the Central Research Institute of Epidemiology of Rospotrebnadzor.

The PCR products were sequenced using the PyroMark Q24 (Qiagen, Germany) and 3500xL (Applied Biosystems, United States) instruments, and reagent kits were recommended by the equipment manufacturers.

2.3 Description of the sample groups

A retrospective study was conducted on biological samples collected between 2017 and 2021 at the Central Institute of Epidemiology. Informed consent was obtained from all patients regarding the use of their DNA samples, and samples without informed consent were excluded from the study. The collected samples were categorized into five distinct groups:

The “Population Control” group consisted of randomly selected whole-blood samples obtained from donors whose sex, age, and HPV infection status were unknown.

The “HPV-Negative” group included women who tested negative for high-risk human papillomavirus (HR-HPV) and had negative results for intraepithelial lesion of malignancy (NILM) during cervical screening over a 5-year monitoring period. The cervical scrapings were used for this assessment.

The “HPV-Positive” group consisted of women who tested positive for high-risk human papillomavirus (HR-HPV) and had negative results for intraepithelial lesion of malignancy (NILM) during cervical screening over a 5-year monitoring period. Within this group, 29% of women showed elimination of the virus, while 71% demonstrated persistent HR-HPV infection.

The “CIN 2/3” group consisted of women who tested positive for high-risk human papillomavirus (HR-HPV) and had cytologically and/or histologically confirmed CIN grade 2 or 3.

The “CC” group comprised women who tested positive for high-risk human papillomavirus (HR-HPV) and had histologically confirmed carcinoma in situ or squamous cell carcinoma of the cervix.

Cytological and histological evaluations, as well as the diagnosis of carcinoma, were conducted based on investigations of cervical scrapings. All the samples included in the study were collected from residents of Moscow, Russia. Samples from the “HPV-Negative,” “HPV-Positive,” “CIN 2/3,” and “CC” groups were recruited between 2017 and 2021. Samples from the “Population Control” group were collected between 2018 and 2020.

Detailed information about the samples is given in Table 2.

TABLE 2

Table 2 Description of the groups.

2.4 Statistical data

The obtained allele data were compared with the allele frequencies from the Ensembl database for Europeans and validated using the Hardy-Weinberg equilibrium and Pearson’s χ2 test. The associations of the SNPs were assessed using logistic regression and expressed as OR with 95% CI for standard genetic models in case-control studies (20). The Akaike information criterion was calculated to determine the most preferred model. The statistical analysis of the results was performed using the R program (21). The linkage disequilibrium analysis, which examines the non-random association of alleles at different loci, was conducted using LDMatrix from the National Cancer Institute (22). To account for multiple testing, the p-values were adjusted using the Bonferroni correction by multiplying the obtained p-values by the number of SNPs and comparisons.

3 Results

The results of genotyping of the samples are given in Table 3.

TABLE 3

Table 3 The results of genotyping of the groups.

In the initial stage of the study, we conducted genotyping of the Population Control group for four markers and compared the allele frequencies with the data obtained from the Ensembl database (23). To determine if there were any significant differences in allele frequencies between the European and Russian populations, a χ2 test was performed. The analysis revealed that the p-value obtained from the χ2 test was greater than 0.2, indicating that no statistically significant differences were observed in the allele frequencies between the European and Russian populations for these markers.

Based on the analysis, no significant differences were found for the markers rs2910164, rs138446575, and rs1048943 between the groups. However, significant differences were observed for the marker rs55986091 in HLA-DQB1 across all four inheritance models. The codominant and recessive models specifically showed the lowest OR of 0.16. The calculation of OR for the five genetic models can be found in Table 4.

TABLE 4

Table 4 Comparison of the «HPV-Negative» and «CC» groups by rs55986091 HLA-DQB1.

The analysis of the AICs and BICs indicated that the log-additive model was the most likely inheritance model for the rs55986091 genetic variant. This suggests that carrying the A allele may provide a protective effect against CC development. On the other hand, the G allele was found to increase the risk of CC (OR=2.49, 95% CI=1.55-4.01).

Specifically, in the Dominant model, individuals carrying the GG and GA genotypes had an increased risk of CC development with an OR of 6.13 (95% CI=2.13-17.65). In the “CC” group, there were 120 individuals with the GG and GA genotypes compared to 4 individuals with the AA genotype. In the “HPV-Negative” group, there were 181 individuals with the GG and GA genotypes compared to 37 individuals with the AA genotype. These findings suggest a significant association between the GG and GA genotypes and an increased risk of CC.

The frequencies of observed genotypes deviated from the Hardy-Weinberg equilibrium (p-value<0.0001).

In order to compare the HPV-Negative and CIN 2/3 groups, we examined four genetic markers. Among these markers, rs55986091 showed a consistent association, with the dominant model having the lowest OR of 0.47. After adjusting for multiple comparisons using the Bonferroni correction, the statistically significant differences remained significant for the dominant and log-additive genetic models.

For further details and specific OR values, refer to Table 5.

TABLE 5

Table 5 «HPV-Negative» and «CIN 2/3» by rs5598609 HLA-DQB1.

The analysis using the Akaike information criterion (AIC) and Bayesian information criterion (BIC) reaffirmed that the log-additive inheritance model is the most likely model for this particular locus.

Furthermore, we compared the HPV-Negative and HPV-Positive groups using four genetic markers. While some statistical differences were observed in the recessive model for rs55986091, these differences were considered insignificant after applying the Bonferroni correction to account for multiple comparisons. For more comprehensive information, please refer to Table 6.

TABLE 6

Table 6 «HPV-Negative» and «HPV-Positive» by rs55986091 HLA-DQB1.

The analysis of allele frequencies revealed that the “Population Control” group had a slightly higher frequency of the A allele compared to women without HPV, but this difference was not statistically significant (p-value>0.05). However, the frequency of the A allele in women with CIN 2/3 was similar to that in the “HPV-positive” group and significantly different from the “HPV-Negative” group (p-value<0.05). For a more detailed overview of the allele frequencies, please refer to Figure 1.

FIGURE 1

Figure 1 Allele frequencies for rs55986091 in all studied groups.

The dominant HPV types in the HPV-infected group were 68, 16, and 31 (in 14%, 13%, and 13% of people, respectively). In the CIN 2/3 group, 16, 31, and 33 HPV types were most often detected (in 53, 19, and 17% of people, respectively). In the group of women with CC, the HPV types often identified were 16, 18, and 33 (in 53, 13, and 8% of people, respectively). Detailed information about the percentage of people infected with different types of HPV in three groups is provided in Supplementary Table S1. The association of HPV type and risk genotype in the three groups was also considered. However, no statistically significant differences were found in the distribution of HPV types depending on the genotype in each of the groups.

4 Discussion

The rs55986091-A allele exhibited a protective effect in the “HPV-Negative” group compared to the “CC” group. The rs55986091 is located in the regulatory region between the HLA-DQA1 and HLA-DQB1 genes (23). In particular, this sequence functions as a super-enhancer for the HLA-DQB1 gene, which is responsible for the transcription and subsequent differentiation of cells (24). The HLA-DQB1 gene is part of the human leukocyte antigen (HLA) complex class II, which plays a crucial role in presenting antigens from external sources to the immune system (25). The upregulation of HLA class II expression and the high concentration of antigen presented by HLA class II molecules on antigen-presenting cells trigger the activation of CD4+ T cells, which is vital for mounting an immune response against pathogens (26).

We propose two hypotheses regarding the influence of this SNP on the development of CC in women:

1. There is a genetic linkage between the observed SNP and a specific allele of the human leukocyte antigen (HLA) gene. This HLA allele is believed to enhance the effectiveness and/or prolong the duration of the immune response in individuals who have received the HPV vaccine.

2. This SNP affects the expression of the HLA-DQB1 gene, whose product may be associated with a protective effect against cancer.

The first hypothesis is supported by several pieces of evidence. Previous studies have demonstrated the protective effects of certain HLA alleles, such as HLA-DQB1*02, HLA-DQB1*0603, DQB1*03, DQB1*0604, DQB1*0501, and DQB1*0603, in CC and other precancerous lesions (27). Carriers of these HLA alleles may exhibit an enhanced immune response after receiving the HPV vaccine, which could contribute to a lower risk of CC development. Similar findings have been observed in the context of hepatitis B, where DQB1*06 variants were associated with a more durable immune response in individuals who received the hepatitis vaccine (28).

One limitation of the present study is the lack of information regarding HPV vaccination in the study group. Adjusting for vaccination status could provide a better understanding of the observed effects. In future research, this limitation should be addressed by incorporating data on HPV vaccination and employing bioinformatics tools to investigate the linkage disequilibrium of this SNP with HLA alleles DQB1 and DQA1.

The second hypothesis also has a basis in the available evidence and should be further investigated in upcoming studies. Given that this SNP is located in the regulatory region of HLA-DQB1, it is possible that it is associated with alterations in gene expression. Reduced or impaired expression of HLA molecules has been linked to the development of CIN and CC (27). Previous studies have also observed associations between polymorphic variants and HLA alleles, such as the linkage disequilibrium between rs2395029 and HLA-B*57:01, which are 100 kb apart (29). Although the connection of CC has been established specifically for rs55986091, it is noteworthy that there are more than 10 other SNPs (R2>0.5) (22) associated with rs55986091 within the super-enhancer region (NC_000006.12:c32651217-32660516). This suggests that variations in the nucleotide sequence spanning over 9,000 bases, particularly those involving the A allele in rs55986091, may have a higher affinity for binding RNA polymerase. Consequently, the immune response triggered in the presence of this allele might be more potent, not only in individuals infected with HPV but also in those affected by other infections (29, 30). These findings are further supported by our data indicating statistical differences in allele frequencies between HPV-negative women and women who were infected and subsequently cleared the virus.

In conclusion, the data indicate an association between this SNP and the immune response to HPV entry into the cell. However, further research is required on the deeper mechanisms underlying this.

Significant deviation from the Hardy-Weinberg equilibrium was observed for rs55986091 (p-value<0.001), indicating a departure from the expected genotype frequencies. Specifically, the frequency of observed heterozygotes was lower than expected, while the frequency of observed homozygotes in both genotypes was slightly higher than expected. Further analysis of the super-enhancer region (NC_000006.12:c32651217-32660516) revealed the presence of at least 100 SNPs with similar genotype distributions.

Based on the findings, it is plausible to suggest that homozygous SNP variants in this region exhibit a higher affinity for recognition by RNA polymerase, resulting in enhanced gene expression and a more robust immune response. Additionally, there may be a genetic linkage between rs55986091 and one of the SNPs situated within the exon, influencing the encoding of the peptide allele (31).

Graffelman et al. previously observed a deficiency of heterozygotes between genes as well as an excess within the genes for HLA class II, which they attributed to the presence of null alleles and sequence errors (32). However, in the case of our SNP, this explanation is not applicable as we did not find any differences between our sample and the Ensembl database. We also conducted pyrosequencing and Sanger sequencing to validate our data. Moreover, the disequilibrium observed persisted across all samples and populations in the Ensembl database. Therefore, we can rule out sequence errors and instead consider explanations such as natural selection or non-random mating as potential factors contributing to the observed disequilibrium. In future studies, we plan to further investigate this aspect by examining genotypes based on the aforementioned SNP in individuals with other infections, as well as measuring the gene expression levels associated with specific genotypes of this SNP.

Additionally, the association of the HPV type with one of the host genotype variants was considered. We did not find a statistically significant association for any of the types in each group. We assume that with an increase in samples, it will become more realistic to establish a relationship between this SNP and the type of HPV. The realization of the carcinogenic potential of the virus is influenced by many factors (smoking, lifestyle, age, and others), including the presence of a genetic predisposition, which is realized by a combination of unfavorable combinations of risk alleles in SNP.

One limitation of this study is the absence of adjustment for demographic factors such as age, smoking status, occupation, vaccination history, and specific HPV types. These factors can potentially influence the development of CC and could have confounding effects on the observed associations. Additionally, the lack of clinical data on patients further limits the comprehensive understanding of the relationships between genetic variants and disease outcomes. Future studies should consider incorporating these demographic factors and clinical data to obtain a more comprehensive analysis.

This study focused exclusively on CC and did not investigate other HPV-associated cancers. Therefore, it is important to replicate the findings in samples of patients with head and neck cancer, including male patients, to determine the association of this SNP with HPV persistence in different cancer types.

The identified associations between specific SNPs and the risk of developing CC have significant implications for the development of screening programs aimed at identifying individuals at higher risk. Additionally, these findings play a crucial role in personalized management approaches for patients infected with high-risk HPV, particularly for individuals with a positive HIV status (8). The personalized approach takes into account individual genetic profiles and can contribute to more effective management and treatment strategies for HPV-associated cancers.

5 Conclusion

In the Russian population, the presence of the rs55986091 variant was observed in AA genotype carriers with an OR of 0.16, and in GG and GA genotype carriers with an OR of 6.13. Identifying both individual and population risks of CC development will improve the monitoring of individuals infected with high-risk HPV (HR-HPV). This will increase their motivation for vaccination and cervical examinations, thereby making a significant contribution to reducing CC-related mortality rates.

The findings highlight the importance of developing and conducting clinical trials for an allele discrimination assay targeting rs55986091 in HR-HPV-positive women. Future research plans involve investigating the underlying mechanisms by measuring the expression level of HLA-DQB1 in individuals with different genotype variants. Additionally, testing other SNPs that have shown statistically significant results in previous studies will be pursued (9). The data obtained from these studies will help in the creation of an analytical algorithm based on machine learning for calculating individual relative risks, further improving risk assessment and personalized management strategies for CC.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics statement

The studies involving human participants were reviewed and approved by The Council on Bioethics of the Central Research Institute of Epidemiology. The patients/participants provided their written informed consent to participate in this study.

Author contributions

MV: Writing - Original Draft, Formal analysis, Investigation, Visualization. KM: Methodology, Data Curation, Writing - Review and Editing, Investigation. ED: Conceptualization, Methodology, Writing - Review and Editing, Data Curation. TR: Conceptualization, Investigation, Writing - Review and Editing. AP: Conceptualization, Data Curation, Writing - Review and Editing. VA: Supervision, Resources, Project administration, Writing - Review and Editing. All authors contributed to the article and approved the submitted version.

Acknowledgments

The team of authors thanks Vitaly I. Korchagin for making editorial edits and Irina A. Shvadlenko for help in translating the article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1207935/full#supplementary-material

Abbreviations

CC, Cervical cancer; HPV-HCR, High-risk oncogenic human papillomavirus; SNP, Single nucleotide polymorphisms; AIC, Akaike information criterion; BIC, Bayesian information criterion.

References

1. Global Cancer Observatory: Cancer Today. Lyon, France: International Agency for Research on Cancer. Available at: https://gco.iarc.fr/today (Accessed 10.07.2022).

Google Scholar

2. Kaprin AD, Starinskiy VV, Shakhzadova AO. The State of Oncological Care to the Population of Russia in 2019 (Sostoya-nie onkologicheskoy pomoshchi naseleniyu Rossii v 2019 go-du). – M.: MNIOI im. Russian: P.A. Herzen, Moscow (2021).

Google Scholar

3. Carozzi F, Visioli CB, Confortini M, Iossa A, Mantellini P, Burroni E, et al. hr-HPV testing in the follow-up of women with cytological abnormalities and negative colposcopy. Br J Cancer (2013) 109(7):1766–74. doi: 10.1038/bjc.2013.519

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Okunade KS. Human papillomavirus and cervical cancer. J Obstet Gynaecol (2020) 40(5):602–8. doi: 10.1080/01443615.2019.1634030

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Duenas-Gonzalez A, Serrano-Olvera A, Cetina L, Coronel J. New molecular targets against cervical cancer. Int J Womens Health (2014) 6:1023–31. doi: 10.2147/IJWH.S49471

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Popova AA, Domonova EA, Vinogradova NA, Yu SO. Anogenital human papillomavirus infection in HIV-infected women (according to the results of a pilot study in the Moscow region). Epidemiologiya i infektsionnye bolezni. Aktu-al’nye voprosy (2021) 11(3):40–5. doi: 10.18565/epidem.2021.11.3.40-5

CrossRef Full Text | Google Scholar

7. Glotov AS, Kazakov SV, Vashukova ES, Pakin VS, Danilova MM, Nasykhova YA, et al. Targeted-sequencing analysis of ACVR2A gene identifies novel risk variants associated with preeclampsia. J Matern Fetal Neonatal Med (2019) 32(17):2790–6. doi: 10.1080/14767058.2018.1449204

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Popova AA, Domonova EA, Pokrovskaya AV, Shipulina OY, Pokrovsk VV. Prevalence of human papillomavirus (Papillomaviridae; Human papillomavirus) of high carcinogenic risk based on the results of screening of three anatomical loci in men stratified by sexual behavior and HIV status. Problems Virol (2021) 66(3):217–26. doi: 10.36233/0507-4088-53

CrossRef Full Text | Google Scholar

9. Vinokurov MA, Mironov KO, Korchagin VI, Popova AA. Genetic polymorphism associated with cervical cancer: a systematic review. J Microbiology Epidemiol immunobiology = Zhurnal mikrobiologii èpidemiologii i immunobiologii (2022) 99(3):353–61. doi: 10.36233/0372-9311-251

CrossRef Full Text | Google Scholar

10. Bowden SJ, Bodinier B, Kalliala I, Zuber V, Vuckovic D, Doulgeraki T, et al. Genetic variation in cervical preinvasiveand invasive disease: a genome-wide association study. LancetOncol (2021) 22(4):548–57. doi: 10.1016/S1470-2045(21)00028-0

CrossRef Full Text | Google Scholar

11. Ramachandran D, Dörk T. Genomic risk factors for cervical cancer. Cancers (Basel) (2021) 13(20):5137. doi: 10.3390/cancers13205137

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Espinoza H, Ha KT, Pham TT, Espinoza JL. Genetic predisposition to persistent human papillomavirus-infection and virus-induced cancers. Microorganisms (2021) 9(10):2092. doi: 10.3390/microorganisms9102092

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sengupta D, Guha U, Mitra S, Ghosh S, Bhattacharjee S, Sengupta M. Meta-analysis of polymorphic variants conferring genetic risk to cervical cancer in Indian women supports CYP1A1 as an important associated locus. Asian Pac J Cancer Prev (2018) 19(8):2071–81. doi: 10.22034/APJCP.2018.19.8.2071

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Ding B, Sun W, Han S, Cai Y, Ren M, Shen Y. Cytochrome P450 1A1 gene polymorphisms and cervical cancer risk: A systematic review and meta-analysis. Med (Baltimore) (2018) 97(13):e0210. doi: 10.1097/MD.0000000000010210

CrossRef Full Text | Google Scholar

15. Hoidy WH, Jaber FA, Al-Askari MA. Association of CYP1A1 rs1048943 polymorphism with prostate cancer in Iraqi men patients. Asian Pac J Cancer Prev (2019) 20(12):3839–42. doi: 10.31557/APJCP.2019.20.12.3839

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhou B, Wang K, Wang Y, Xi M, Zhang Z, Song Y, et al. Common genetic polymorphisms in pre-microRNAs and risk of cervical squamous cell carcinoma. Mol Carcinog (2011) 50(7):499–505. doi: 10.1002/mc.20740

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Yue C, Wang M, Ding B, Wang W, Fu S, Zhou D, et al. Polymorphism of the pre-miR-146a is associated with risk of cervical cancer in a Chinese population. Gynecol Oncol (2011) 122(1):33–7. doi: 10.1016/j.ygyno.2011.03.032

PubMed Abstract | CrossRef Full Text | Google Scholar

18. You Y, Moreira BG, Behlke MA, Owczarzy R. Design of LNA probes that improve mismatch discrimination. Nucleic Acids Res (2006) 34(8):e60. doi: 10.1093/nar/gkl175

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Axelrod EV, Mironov KO, Dunaeva EA, Shipulin GA. (The comparison of three molecular genetic techniques for identifying major mutations in gene HFE related to development of inherent hemochromatosis.). Klin Lab Diagn (2016) 61(5):316–20. doi: 10.18821/0869-2084-2016-61-5-316-320

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Horita N, Kaneko T. Genetic model selection for a case-control study and a meta-analysis. Meta Gene (2015) 5:1–8. doi: 10.1016/j.mgene.2015.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Stiglic G, Watson R, Cilar L. R you ready? Using the R programme for statistical analysis and graphics. Res Nurs Health (2019) 42(6):494–9. doi: 10.1002/nur.21990

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Machiela MJ, Chanock SJ. LDassoc: an online tool for interactively exploring genome-wide association study results and prioritizing variants for functional investigation. Bioinf (2018) 34(5):887–9. doi: 10.1093/bioinformatics/btx561

CrossRef Full Text | Google Scholar

23. Ensemble genome browser. Available at: https://www.ensembl.org/index.html/ (Accessed 02.04.2022).

Google Scholar

24. Khan A, Zhang X. dbSUPER: a database of super-enhancers in mouse and human genome. Nucleic Acids Res (2016) 44(D1):D164–71. doi: 10.1093/nar/gkv1002

PubMed Abstract | CrossRef Full Text | Google Scholar

25. UniProt database (2022). Available at: http://www.uniprot.org/ (Accessed 02.05. 2022).

Google Scholar

26. Shiina T, Hosomichi K, Inoko H. The HLA genomic loci map: expression, interaction, diversity and disease. J Hum Genet (2009) 54(1):15–39. doi: 10.1038/jhg.2008.5

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Paaso A, Jaakola A, Syrjänen S, Louvanto K. From HPV infection to lesion progression: the role of HLA alleles and host immunity. Acta Cytol (2019) 63(2):148–58. doi: 10.1159/000494985

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Wang WC, Lin YS, Chang YF, Yeh CC, Su CT, Wu JS, et al. Association of HLA-DPA1, HLA-DPB1, and HLA-DQB1 alleles with the long-term and booster immune responses of young adults vaccinated against the Hepatitis B virus as neonates. Front Immunol (2021) 12:710414. doi: 10.3389/fimmu.2021.710414

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Vince N, Douillard V, Geffard E, Meyer D, Castelli EC, Mack SJ, et al. SNP-HLA Reference Consortium (SHLARC): HLA and SNP data sharing for promoting MHC-centric analyzes in genomics. Genet Epidemiol (2020) 44(7):733–40. doi: 10.1002/gepi.22334

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Enz LS, Zeis T, Schmid D, Geier F, van der Meer F, Steiner G, et al. Increased HLA-DR expression and cortical demyelination in MS links with HLA-DR15. Neurol Neuroimmunol Neuroinflamm (2019) 7(2):e656. doi: 10.1212/NXI.0000000000000656

PubMed Abstract | CrossRef Full Text | Google Scholar

31. HLA Nomenclature (2020). Available at: https://hla.alleles.org/ (Accessed 02.06.2021).

Google Scholar

32. Graffelman J, Jain D, Weir B. A genome-wide study of Hardy-Weinberg equilibrium with next generation sequence data. Hum Genet (2017) 136(6):727–41. doi: 10.1007/s00439-017-1786-7

PubMed Abstract | CrossRef Full Text | Google Scholar