Introduction
============

According to WHO statistics, invasive cervical cancer (ICC) still ranks as the fourth most prevalent cancer in women and more than 80% cases arise in developing countries [@B1],[@B2]. In recent years, the worldwide incidence and mortality of ICC have generally decreased, owing to the successful implement of cervical cancer screening program and great achievements in early diagnosis and treatment of premalignant conditions. However, the incidence of ICC in young women has been reported to be increasing steadily both in China and in other countries [@B3]-[@B7]. In most reports, young women seemed to have a more aggressive form of the disease and a poorer prognosis than older women [@B3]-[@B6]. Therefore, studies aimed to investigate the possible pathogenic factors responsible for young patients with ICC are of great concern for their early prevention.

Genital infection of human papillomaviruses (HPVs), particularly the persistent infection status of high-risk HPVs is recognized as a causal and necessary factor in the pathogenesis of ICC. It is well known that high-risk HPV infection is very popular among women at reproductive age. Human papillomavirus early proteins E6 and E7 subvert normal cell control mechanisms through interactions with two oncosuppressor proteins, p53 and Rb, respectively [@B8]. Human papillomavirus type 16 infection was found to account for about 50% of ICC worldwide [@B2],[@B9],[@B10]. However, majority of high-risk HPV infections are transient and disappear spontaneously, only a small percentage of HPV infections persist and progress to precancerous diseases or ICC, indicating that additional factors are needed during the pathogenic process. Both viral and host genetic factors of infected individuals were considered to be important factors in determining progression to neoplasia of high-risk HPV infection.

Recently, different HPV16 E6 variants were observed to be related to an increased risk for cervical ICC in a certain population [@B11]-[@B13]. Distinct geographic diversity was observed in the distribution of HPV16 E6 variants. Currently, the variants were classified into four lineages [@B14]: A, which includes A1-A3 (previously called European, EUR) and A4 (previously named Asian, As); B (previously named African 1, AF-1); C (previously named African 2, AF-2); D, which includes Asian-American (AA) and North-American (NA) lineages named previously. Based on nucleotide polymorphism at nt350, the A1-A3 lineage were further classified into 350T and 350G.The A4 lineage was seldom reported in the west, but in some Asian countries such as Japanese, Korean and the southern area of China, it has becoming the most common variant found in cervical cancer cases [@B15]-[@B18]. We investigated the distribution of HPV16 E6 variants in cervical squamous cell carcinoma (SCC) patients from a coastal city located in the southern part of China and found that the A4 lineage was the most common variant type [@B19]. In addition, we also found that the distribution of HPV16 E6 variants was significantly different between younger cervical SCCs (≤35ys) and older patients (\>35ys). In young patients group, we only detected the A4 lineage. However, other linages of HPV 16 E6 variants were found in the older cervical SCCs in addition to the A4 lineage. The unmatched distribution of the A4 lineage observed in young SCC patients indicated the possible oncogenicity of this variant type might be more aggressive in young cervical SCCs of this study population.

On the other hand, the host genetic diversity may also play an important part in the natural process of HPV infections. The distinct characteristics of HPV16 E6 variants distribution observed in cervical diseases among ethnically different populations suggest a possible relationship between E6 variants and human leukocyte antigens (HLA), which are population related. As we know, HLA plays a key role during the clearance of virally infected cells by presenting xenogenic antigens to the cytotoxic T cells. The highly polymorphic character of HLA genes may display different ability in the clearance of HPV-infected cells, and hence the progression to cervical cancer. HLA class II provides the main genetic susceptibility loci to cervical cancer. The strongest effects are from the HLA II DQB1 and DRB1 loci [@B20],[@B21]. Several researches have analyzed the possible correlation of HPV16 E6 variants with HLA class II polymorphisms in patients with ICC from different populations and found certain HLA alleles showed to be correlated to a changed risk for ICC among different ethnic population [@B22]-[@B24]. Although our previous study has implied that the uniform presence and distinct raised E6 A4 lineage in young age cervical SCCs may be closely related to the early occurred cervical SCC, so far little information is available about the HLA class II polymorphism in relation with the HPV16 E6 variants in young women aged 35 or younger with cervical SCC. In the present study, we examined the difference of DRB1, DQB1 polymorphisms between young cervical SCCs and older ones and the association with E6 variants, with the aim of discovering the possible interaction between E6 variants and HLA class II variability during the obviously accelerated progression to cervical cancer in young women.

Materials and Methods
=====================

Study Population
----------------

The specimens of cervical SCCs in this study were collected at the Gynecology ward of 1st Affiliated Hospital of Wenzhou Medical University between June 2009 and June 2015. Sixty-five women aged ≤35 while 713 patients aged ≥36 at the time of diagnosis. All members of this study had access to identifying information during data collection. Each patient enrolled in this study needed to be a native resident in Wenzhou city, which is a coastal city located in the south of China. Sixty-one women aged ≤35 were enrolled in group I, 85 patients aged ≥36 (selected by an 8-patient interval sampling design) in group II. The cut-off limit of 35 years for young cervical cancer has been applied by many researchers [@B3],[@B6]. The mean ages of the two groups were 32.4±2.9 years and 51.1±6.7 years respectively. At the same time, we invited 130 age-matched healthy volunteers who came for routing cervical disease screening to take part in our study as healthy controls to provide their cervical cellular specimen and blood samples. The ethnicity, marital status, and smoking habits of controls were similar to the patients. The average age of the control group was 39.8±10.3 years.

Research Process and Sample Collection
--------------------------------------

The research plan was permitted by the research ethics committee of the 1st Affiliated Hospital of Wenzhou Medical University (20040912). All patients and controls gave written informed consent. The medical, gynecological, sexual and childbearing history were provided by all patients and controls. Physical and gynecological examinations were performed afterwards. For all cervical SCC patients, peripheral serum samples and cervical cytology samples were collected before surgical treatment. The cervical cytology samples were sent for HPV genotyping. Their cervical tissue samples were collected during the surgery and were placed immediately in liquid nitrogen until required for analysis. Of the healthy controls, peripheral serum samples and cervical cytology samples were collected. The cervical cytology samples were sent for Liquid-based cytological test and HPV genotyping. Blood samples drawn from all cervical SCCs and healthy controls were analyzed for DRB1 and DQB1 loci.

HPV Genotyping
--------------

The cervical cytology samples of all cervical SCCs and healthy controls were examined for HPV genotypes as described earlier [@B19].

HPV 16 E6 Variants Analysis
---------------------------

The cervical DNA samples were extracted in HPV16 positive cervical SCC patients and sequenced as described earlier [@B19]. The BLAST 2.0 software server (<http://www.ncbi.nih.gov/blast>) was used to analyze DNA sequence. The variants were classified into four lineages as A (A1-3, A4), B, C and D. Based on nucleotide polymorphism at nt350, the A1-3 lineage were further classified into 350T and 350G.

Human Leukocyte Antigen Typing
------------------------------

Genomic DNA of patients and healthy controls were isolated from peripheral venous blood samples by using the DNeasy Blood & Tissue Kit (QIAGEN, Germany). Genotyping of HLA-DRB1 and HLA-DQB1 as low-resolution allele distribution was determined by PCR sequence specific primers (SSP) method. DRB1and DQB1 PCR-SSP Morgan^TM^ Kits were supplied by Texas BioGene Inc. (Texas, USA). Amplification reactions were done according to manufacturer\'s protocols. Results were photographed with a gel-imaging system and analyzed by Morgan™ SSPal HLA Typing Analysis Software.

Statistical Analyses
--------------------

Statistical analyses were conducted using Statistical Package for the Social Sciences (SPSS) 22.0 software. Allele frequencies of DRB1and DQB1 were calculated by direct counting. Chi-squared (*x*^2^) comparisons were performed to identify any differences on categorical variables. Odds ratios (ORs), 95% confidence intervals (95% CIs) and *p* values were calculated. According to Bonferroni\'s correction, corrected *p* value (Pc) was calculated for multiple comparisons, which multiplied the *p* value by the number of alleles tested for each locus. All *p* values were assumed significant as *p*\<0.05.

Results
=======

Characteristics of the Cervical SCC Patients
--------------------------------------------

Demographic characteristics between two cervical SCC groups were compared (Table [1](#T1){ref-type="table"}). No significant difference was found except for the number of pregnancies.

Details of HPV genotype in cervical cytology samples of 146 cervical SCCs and 130 healthy controls were summed up (Table [2](#T2){ref-type="table"}). The infection rate of HPV was 96.6% (141/146) and 11.5% (15 /130) in cervical SCCs and healthy controls respectively. The most common genotype found in this population was HPV 16. 45.21% (66/146) of cervical SCC patients and 4.6% (6/130) of controls were HPV 16 positive. The percentage of HPV16-positive patients in group I (57.4%, 35/61) was significantly higher than that of group II (36.5%, 31/85) (*x*^2^=6.267, *p*\<0.05, OR=2.345, 95% CI: 1.197-4.595).

Characteristics of HPV16 E6 variants
------------------------------------

66 cervical SCCs and 6 healthy controls with HPV16 infection were further determined for E6 sequence variation (Table [3](#T3){ref-type="table"}).

The A4 lineage was the most predominant variant found in this study. 63.6% (42/66) cervical SCCs and 16.7% (1/6) healthy controls contained this variant. The other non-synonymous variations detected were the B and A1-3 lineages. The E6 variants distribution was further analyzed in cervical SCCs. In group I, the A4 lineage was the unique variant detected. In addition to the A4 lineage, the B and A1-3 lineages were also observed in group II. Although the A4 lineage was the most common variant observed in both groups, the distribution was obviously different between the two groups. The occurrence rate was 82.9% (29/35) in group I, whilst 41.9% (13/31) in group II (*x*^2^=11.897, *p*\< 0.05).

Analysis of HLA-DRB1 and HLA-DQB1 Polymorphisms
-----------------------------------------------

Data of the HLA-DRB1 and HLA-DQB1 alleles from 146 cervical SCCs and 130 healthy controls were summarized (Table [4](#T4){ref-type="table"}). Hardy-Weinberg equilibrium was found in the allele frequencies of the 2 loci among all healthy controls. The most frequent allele for the HLA-DRB1 loci found in cervical SCCs was DRB1\*09 (15.1%, 44/292), whilst DRB1\*15 (13.5%, 35/260) in healthy controls. The allele frequencies among cervical SCCs were obviously higher than that among healthy controls for DRB1\*11 and DQB1\*0601 alleles, with a Pc value of 0.013 in the former and 0.016 in the latter. However, the allele frequencies for DRB1\*13 and DQB1\*02 were significantly decreased in cervical SCCs as compared with healthy controls, with a Pc value of 0.000 in the former and 0.04 in the latter. Comparison of the other DRB1 and DQB1 allele frequencies between cervical SCCs and healthy controls revealed no significant statistical difference (*p*/Pc\>0.05).

Meanwhile, differences of DRB1 and DQB1 polymorphisms distribution between HPV16-positive cervical SCCs and healthy controls were compared simultaneously (Table [4](#T4){ref-type="table"}). Among HPV16-positive cervical SCCs, the most frequent alleles for DRB1 and DQB1 loci observed were DRB1\*08 (17.4%, 23/132) and DQB1\*0601 (22.7%, 30/132). Negative associations were found for DRB1\*13 and DQB1\*02 alleles with respect to HPV16-positive cervical SCCs, with a Pc value of 0.026 in the former and 0.04 in the latter. DQB1\*0601 allele was observed to be positively related to HPV16 infection, which occurred significantly more often in HPV16-positive cervical SCCs (Pc=0.000). In addition, a positive relationship of DRB1\*11 with HPV16-positive cervical SCCs was suggestive (*p*=0.017), but no significant difference was kept after Bonferroni\'s correction (Pc=0.221). No significance was observed for the other DRB1/ DQB1 alleles in HPV16-positive cervical SCCs (*p*/Pc\>0.05).

Differences of HLA-DRB1 and HLA-DQB1 Polymorphisms between Group I and Group II
-------------------------------------------------------------------------------

In order to explore a putative link between specific HLA-DRB1 or HLA-DQB1 alleles and the early onset cervical SCCs, differences between younger cervical SCCs and older patients in allele distribution and frequency were further compared (Table [4](#T4){ref-type="table"}). The most frequent allele detected for the DRB1 loci in group I was DRB1\*09 (20.5%, 25/122), whereas DRB1\*08 (14.1%, 24/170) and DRB1\*12 (14.1%, 24/170) were the two most frequent alleles found in group II. It is noteworthy that the allele frequencies of DRB1\*04 and DRB1\*09 were significantly higher in group I compared with those in group II, whereas DRB1\*07 frequency in group I was significantly lower than that in group II. However, statistical significances were all disappeared after the application of Bonferroni\'s correction(Pc\>0.05).

For the HLA-DQB1 loci, the allelic distribution found that DQB1\*0301(30.3%, 37/122) was the most frequent allele in group I, whereas DQB1\*0501 (23.5%, 40/170) was the most frequent allele in group II. In addition to this, the frequency of DQB1\*0501 was obviously lower in group I than that in group II, with a Pc value of 0.000. A higher frequency of DQB1\*0301 was found in group I than that in group II (*p*=0.016), but no significant difference was kept after Bonferroni\'s correction (Pc=0.128).

The Relation of HLA-DRB1 and HLA-DQB1 Polymorphisms with HPV16 E6 Variants
--------------------------------------------------------------------------

With the aim to identify a possible relationship of HLA-DRB1, DQB1 polymorphisms with specific E6 variants, we further compared the allele distribution in HPV16-positive cervical SCCs and healthy controls (Table [5](#T5){ref-type="table"}).

Special associations between three alleles and certain E6 variants were observed. DRB1\*09 and DQB1\*0601 frequencies were remarkably increased among cervical SCCs with the A4 lineage compared with those among healthy controls, with a Pc value of 0.039 in the former and 0.000 in the latter. The allele frequency of DRB1\*11 was obviously higher among cervical SCCs with the B variant compared with that among healthy controls, with a Pc value of 0.026. Additionally, we found that DRB1\*08 and DQB1\*0301 frequencies were more common in cervical SCCs with the A1-3(350T) lineage compared with healthy controls (*p*\<0.05), but no significant difference was kept after Bonferroni\'s correction (Pc\>0.05). Similar tendencies were also observed in DQB1\*0303 among cervical SCCs with the A4 lineage and in DRB1\*12 in cervical SCCs with the B lineage compared with healthy controls (*p*\<0.05). But likewise, no significant difference was kept after Bonferroni\'s correction (Pc\>0.05).

The Distribution of HLA-DRB1 and HLA-DQB1 Alleles Associated with the HPV16 A4 Lineage between Young Cervical SCCs and Older Patients
-------------------------------------------------------------------------------------------------------------------------------------

As the A4 lineage was the most prevalent type observed among patients with HPV16-positive cervical SCC in both groups, and was the exclusive genotype detected in group I, we further compared the distribution of DRB1, DQB1 alleles associated with the A4 lineage among cervical SCC patients (Table [6](#T6){ref-type="table"}).

A significant association was observed for DRB1\*03, with a Pc value of 0.039. This allele was not found among cervical SCCs positive for the A4 lineage in group I, whereas it was detected among those in group II with the percentage of 19.2% (5/26). It is worthwhile to note that DRB1\*08 frequency was significantly increased among cervical SCCs positive for the A4 lineage in group I compared with those in group II (*p*\<0.05). By contrast, DQB1\*02 frequency was decreased among patients with the A4 lineage in group I (*p*\<0.05). However, statistical significances were lost after adjustment made for the other alleles (Pc\>0.05). Additionally, allele frequencies of DQB1\*04 and DQB1\*0601were increased among cervical SCCs with the A4 lineage in group I, but no statistical significance was observed (*p*\>0.05). Although HLA-DRB1\*09 frequency was significantly higher in cervical SCCs with A4 lineage than controls, no significant difference was found between group I and group II (*p*=0.946).

Discussion
==========

It is generally acknowledged that the time course from high-risk HPV infection to the development of invasive cervical cancer usually takes about 15\~20 years [@B25]. The obvious earlier onset age of cervical SCC among young patients suggests the existence of some other important factors associated with the accelerated carcinogenesis. It is possible that genetic factors influence viral persistence, and specific more virulent variants of HPV may play important roles in early onset cervical SCCs. In the present study, we analyzed the distribution of HPV16 E6 sequence variants and HLA-DRB1, -DQB1 allele frequencies among a cervical SCC population in China. The dissimilarities of DRB1, DQB1 polymorphisms between young cervical SCCs and older cases and the correlation with HPV16 E6 variants were also analyzed.

In agreement with our previous study, the A4 lineage was also the most prevalent genotype observed (63.6%). We observed that the E6 variants distributed significantly different between young cervical SCC patients and older cases. The A4 lineage was the exclusive type observed in young cervical SCCs, and its prevalence was significantly higher than that of older cases (82.9% vs. 41.9%, *p*\<0.05). The diverse distribution of the A4 lineage between younger cervical SCCs and older ones occurred in the same region indicated that this lineage was probably more carcinogenic in young women. Results of a few studies [@B19],[@B26] also indicated that young women appeared to be more vulnerable or susceptible to certain E6 variant genotypes in a specific population. The host genetic diversity, for example the HLA alleles polymorphisms may be used to explain this difference. Some HLA molecules might be more efficient in presenting epitopes from certain E6 variant genotypes and thus determine its potential carcinogenicity [@B23].

Previous studies from diverse ethnic populations have analyzed the possible correlation of HLA II molecules with cervical cancer [@B22],[@B27],[@B28]. Although the same HLA class II alleles may be detected frequently to be associated with the occurrence of cervical cancer in different populations, there\'s still no consensus on which HLA II alleles may increase the risk of cervical carcinoma in any population. In the present study, allele frequencies of DRB1\*11 and DQB1\*0601 increased significantly in cervical SCCs when compared with healthy controls, whereas DRB1\*13 and DQB1\*02 were significant decreased in cervical SCCs. In addition, a similar association was observed for DRB1\*13 and DQB1\*0601 in HPV16-positive cervical SCCs when compared with healthy controls. These results suggested the protective effect of DRB1\*13 and the predisposing effect of DQB1\*0601 may be specific for HPV16 in this study population.

To further explore possible accelerating factors associated with early onset cervical SCC in young population, the investigation of both HPV intra-typic variations and HLA-DRB1, HLA-DQB1 polymorphisms between young cervical SCCs and older patients, as demonstrated in this study, is necessary. To the best of our current knowledge, this study is the first comprehensive investigation regarding the difference of DRB1, DQB1 polymorphisms between young cervical SCCs and older patients. The results indicated different distribution of certain DRB1, DQB1 alleles between young cervical SCCs and older ones. The allele frequencies for DRB1\*04, DRB1\*09 and DQB1\*0301 among young cervical SCCs were significantly increased when compared with older ones (*p*\<0.05). On the other hand, DRB1\*07 and DQB1\*0501 in young cervical SCCs were significantly decreased when compared with older ones (*p*\<0.05). After adjustment for the other alleles was made, statistical significance was kept only for DQB1\*0501(Pc\<0.05), thereby indicating its protective effect. This result suggests that young women who are carriers of HLA-DQB1\*0501 allele may have a relatively lower risk for the development of an early onset cervical SCC.

Based on the above findings between young cervical SCCs and older patients, we hypothesize that certain DRB1, DQB1 alleles may be associated to special HPV16 E6 variants, and this putative relationship may contribute to the obviously accelerated carcinogenesis in this young Chinese population with cervical SCC.

Concerning the possible correlation between the E6 variant genotype and HLA polymorphisms in the development of cervical cancer, a number of studies from different ethnic populations have found that certain HLA class II alleles related with specific HPV variants showed a predisposing effect or a protective effect in cervical carcinogenesis [@B23],[@B24],[@B29],[@B30]. However, we know little about a possible relation between specific HPV16 E6 variant and HLA-DRB1, HLA-DQB1 alleles in the obviously accelerated carcinogenesis among young women with cervical SCC.

In the present study, we found that the allele frequencies of DRB1\*09 and DQB1\*0601 were significantly higher among cervical SCC patients with the A4 lineage when compared with healthy controls. This result suggested that carriers of DRB1\*09 or DQB1\*0601 alleles may be at increased risk for developing cervical SCC once infected with the A4 lineage. Additionally, DRB1\*11 was found to be positively associated with the B variant.

As we have noticed before, the A4 lineage was not only the most common type of E6 variant observed among this study population, but it was also the unique variant found in young cervical SCCs. Then we further extended the type-specific analysis concerning the distribution of DRB1, DQB1 alleles associated with the A4 lineage between young cervical SCCs and older patients. The DRB1\*03 was not detected among young patients positive for the A4 lineage. A significant higher frequency of DRB1\*03 was found among older patients positive for the A4 lineage compared with that among young patients (Pc\<0.05). This result indicates that young women who are carriers of HLA-DRB1\*03 allele may at lower risk of developing an early onset cervical SCC associated with the HPV16 A4 lineage in this Chinese population. In addition to this, an increased frequency of DRB1\*08 and a decreased frequency of DQB1\*02 were also observed among young patients positive for the A4 lineage compared with those among older patients, with *p* values of \<0.05. However, statistical significance was not kept upon Bonferroni\'s correction. These results suggested that these alleles may display a possible predisposing or protective effect in the carcinogenesis of an early onset cervical SCC associated with the HPV16 A4 lineage in this Chinese population. Failure to get a significant correlation may due to the modest sample size.

Conclusion
==========

The present data not only provided evidence for the difference of DRB1, DQB1 polymorphisms between young patients with cervical SCC and older ones, but also gave clues on the possible relationship between certain HPV16 E6 variants and HLA class II alleles. In addition, our data further implied that certain HLA-DRB1 and HLA-DQB1 alleles may be related to the A4 lineage among young cervical SCCs, which was the unique HPV16 E6 variant found in young patients. These results suggested that the early onset cervical cancer in this study population was, at least in part, affected by specific association between certain HPV variant and genetic polymorphisms of HLA. In order to further verify the possible link between certain HPV16 E6 variant and HLA-DRB1, HLA-DQB1 alleles and the relationship with the early occurred cervical SCCs in young population, studies on larger sample size will be required. Such information may help to elucidate the pathogenic factors that associated with young women cervical SCCs, and thus to provide more accurate approaches for better prophylaxis and treatment.

This work was supported by grants from the Natural Science Foundation of Zhejiang Province, China (Grant No. Y204080).

HPV

:   human papillomavirus

HLA

:   human leukocyte antigen

ICC

:   invasive cervical cancer

SCC

:   squamous cell carcinoma

SSP

:   sequence specific primers.

###### 

Comparison of demographic characteristics between young cervical SCC patients and older ones.

  Characteristics                     Group I (n=61)        Group II (n=85)
  ----------------------------------- --------------------- ---------------------
  **Ethnicity**                       100% Han population   100% Han population
  **Mean age**                        32.4±2.9years         51.1±6.7years
  **Mean age at first intercourse**   20.6±1.7years         21.1±2.2 years
  **Sexual partner numbers**                                
  **1**                               86.7%                 85.1%
  **≥ 1**                             13.3%                 14.9%
  **Number of pregnancies**           1∼3                   1∼8
  **Smoking habits**                                        
  **Yes**                             2.3%                  3.1%
  **No**                              97.7%                 96.9%

###### 

Comparison of HPV Genotyping between Cervical SCCs and Healthy Controls.

  HPV Genotype                             Cervical SCCs   Healthy Controls       
  ---------------------------------------- --------------- ------------------ --- ------
  **Single Infection:**                                                           
  **high-risk type**                                                              
  **16**                                   57              40.4               6   40.0
  **18**                                   9               6.4                1   6.7
  **52**                                   13              9.2                3   20.0
  **58**                                   12              8.5                1   6.7
  **31**                                   10              7.1                0   0.0
  **33**                                   8               5.7                1   6.7
  **45**                                   4               2.8                0   0.0
  **51**                                   3               2.1                0   0.0
  **56**                                   1               0.7                0   0.0
  **35**                                   2               1.4                0   0.0
  **59**                                   1               0.7                0   0.0
  **39**                                   3               2.1                1   6.7
  **68**                                   2               1.4                0   0.0
  **probable high-risk type**                                                     
  **53**                                   0               0.0                0   0.0
  **66**                                   1               0.7                0   0.0
  **low-risk type**                                                               
  **6/11/42/43/44/81**                     3               2.1                1   6.7
  **Multiple Infection with type 16**                                             
  **16+52**                                3               2.1                0   0.0
  **16+18**                                2               1.4                0   0.0
  **16+33**                                2               1.4                0   0.0
  **16+58**                                1               0.7                0   0.0
  **16+31**                                1               0.7                0   0.0
  **Multiple Infection without type 16**   3               2.1                1   6.7

###### 

Comparison of E6 Variants between Cervical SCCs and Healthy controls.

  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Nucleotide position   Predicted amino acid changes   Variants   Cervical SCCs   Healthy controls\   Total\                                                                      
                                                                                  N (%)               \                                                                           
                                                                                                      N (%)                                                                       
  --------------------- ------------------------------ ---------- --------------- ------------------- -------- ------- ------------ ------------ ---------- ---------- ---------- ----------
  **C**                 **T**                          **T**      **C**           **T**               **G**    **A**   \-           A1-3(350T)   6(17.1)    5(16.1)    3(50.0)    14(19.4)

  **-**                 **-**                          **G**      **-**           **-**               **-**    **-**   \-           A1-3(350T)   0          2(6.5)     1(16. 7)   3(4.2)

  **-**                 **G**                          **-**      **-**           **-**               **-**    **-**   D25E         A4           23(65.7)   12(38.7)   1(16. 7)   36(50.0)

  **-**                 **G**                          **-**      **-**           **-**               **T**    **-**   D25E/E97D    A4           6(17.1)    1(3.2)     0          7(9.7)

  **-**                 **-**                          **-**      **-**           **G**               **-**    **-**   L83V         A1-3(350G)   0          5(16.1)    1(16. 7)   6(8.3)

  **-**                 **G**                          **-**      **-**           **G**               **-**    **-**   T22S/L83V    A1-3(350G)   0          1(3.2)     0          1(1.4)

  **-**                 **-**                          **-**      **T**           **-**               **-**    **C**   H78Y/E113D   B            0          5(16.2)    0          5(6.9)
  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Abbreviations: SCC, squamous cervical cancer

###### 

HLA-DRB1 and HLA-DQB1 Allele Frequencies among Cervical SCCs and Healthy Controls.

  ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Alleles      Healthy Controls\   Cervical SCCs\   HPV16+\    Group I\   Group II\   Cervical SCCs\        HPV16+ SCCs\            Group I vs Group II                                                 
               2n=260\             2n=292\          SCCs\      2n=122\    2n=170\     vs\                   vs\                                                                                         
               (%)                 (%)              2n=132\    (%)        (%)         Healthy Controls      Healthy Controls                                                                            
                                                    (%)                                                                                                                                                 
  ------------ ------------------- ---------------- ---------- ---------- ----------- --------------------- ----------------------- --------------------- ---------------------- ---------------------- ----------------------
  **DRB1**                                                                                                                                                                                              

  **\*01**     8(3.1)              2(0.7)           1(0.8)     0(0.0)     2(1.2)      0.074                 4.603(0.969\~21.876)    0.275                 4.159(0.515\~33.608)   ---                    ---

  **\*03**     23(8.8)             25(8.6)          9(6.8)     11(9.0)    14(8.2)     0.906                 1.036(0.573\~1.875)     0.488                 1.326(0.595\~2.954)    0.814                  1.104(0.483\~2.523)

  **\*04**     16(6.2)             24(8.2)          10(7.6)    15(12.3)   9(5.3)      0.350                 0.732(0.380\~1.411)     0.593                 0.800(0.353\~1.815)    [0.032]{.ul}^a^        2.508(1.059\~5.937)

  **\*07**     21(8.1)             17(5.8)          8(6.1)     3(2.5)     14(8.2)     0.296                 1.421(0.733\~2.757)     0.471                 1.362(0.586\~3.163)    [0.038]{.ul}^b^        0.281(0.079\~1.000)

  **\*08**     29(11.2)            38(13.0)         23(17.4)   14(11.5)   24(14.1)    0.504                 0.839(0.501\~1.405)     0.084                 0.595(0.329\~1.076)    0.508                  0.789(0.390\~1.595)

  **\*09**\    26(10.0)            44(15.1)         19(14.4)   25(20.5)   19(11.2)    0.074                 0.626(0.374\~1.050)     0.197                 0.661(0.351\~1.244)    [0.028]{.ul}^c^        2.048(1.071\~3.919)
  **\*10**                                                                                                                                                                                              

  8(3.1)       4(1.4)              2(1.5)           1(0.8)     3(1.8)     0.170       2.286(0.680\~7.681)   0.354                   2.063(0.432\~9.858)   0.079                  3.220(0.816\~12.706)   

  **\*11**     10(3.8)             34(10.7)         13(9.8)    13(11.5)   21(12.4)    [0.001]{.ul}^d^       0.304(0.147\~0.627)     [0.017]{.ul}^e^       0.366(0.156\~0.859)    0.656                  0.846(0.406\~1.764)

  **\*12**     29(11.2)            42(14.4)         18(13.6)   18(14.8)   24(14.1)    0.258                 0.747(0.451\~1.239)     0.475                 0.795(0.424\~1.4920    0.879                  1.053(0.544\~2.039)

  **\*13**     18(6.9)             2(0.7)           0(0.0)     2(1.6)     0(0.0)      [0.000]{.ul}^f^       10.785(2.478\~46.945)   [0.002]{.ul}^g^       0.931(0.900\~0.962)    ---                    ---

  **\*14**     21(8.1)             16 (5.5)         7(5.3)     3(2.5)     13(7.6)     0.268                 1.461(0.745\~2.865)     0.314                 1.569(0.649\~3.792)    0.055                  0.304(0.085\~1.093)

  **\*15**     35(13.5)            31(10.6)         16(12.1)   14(11.5)   17(10.0)    0.304                 1.310(0.782\~2.192)     0.709                 1.128(0.599\~2.123)    0.686                  1.167(0.552\~2.467)

  **\*16**     16(6.2)             13(4.45)         6(4.5)     3(2.5)     10(5.9)     0.371                 1.407(0.664\~2.984)     0.513                 1.377(0.526\~3.606)    0.162                  0.403(0.109\~1.498)

  **DQB1**                                                                                                                                                                                              

  **\*02**     55(21.2)            36(12.3)         13(9.8)    14(11.5)   22(12.9)    [0.005]{.ul}^h^       1.908(1.206\~3.018)     [0.005]{.ul}^i^       2.456(1.288\~4.682)    0.707                  0.872(0.427\~1.782)

  **\*0301**   65(25.0)            68(23.3)         27(20.5)   37(30.3)   31(18.2)    0.639                 1.098(0.743\~1.622)     0.316                 1.296(0.780\~2.154)    [0.016]{.ul}^j^        1.952(1.128\~3.377)

  **\*0302**   21(8.1)             18(6.2)          7(5.3)     5(4.1)     13(7.6)     0.381                 1.338(0.696\~2.570)     0.314                 1.569(0.649\~3.792)    0.214                  0.516(0.179\~1.488)

  **\*0303**   31(11.9)            44(15.1)         23(17.4)   18(14.8)   26(15.3)    0.282                 0.763(0.466\~1.250)     0.135                 0.642(0.357\~1.152)    0.899                  0.959(0.500\~1.839)

  **\*0306**   8(3.1)              4(1.4)           0(0.0)     2(1.6)     2(1.2)      0.170                 2.286(0.080\~7.681)     0.097                 0.969(0.948\~0.990)    0.737                  1.400(0.194\~10.078)

  **\*04**     21(8.1)             19(6.5)          9(6.8)     12(9.8)    7(4.1)      0.478                 1.262(0.663\~2.405)     0.658                 1.201(0.534\~2.701)    0.051                  2.540(0.970\~6.654)

  **\*0501**   34(13.1)            48(16.4)         23(17.4)   8(6.6)     40(23.5)    0.268                 0.765(0.476\~1.230)     0.249                 0.713(0.401\~1.269)    [0.000]{.ul}^k^        0.228(0.103\~0.507)

  **\*0601**   25(9.6)             55(18.8)         30(22.7)   26(21.3)   29(17.1)    [0.002]{.ul}^l^       0.458(0.276\~0.760)     [0.000]{.ul}^m^       0.362(0.203\~0.646)    0.359                  1.317(0.730\~2.374)
  ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Abbreviations: HPV, human papillomavirus; SCC, squamous cervical cancer

\*The valid sample size equaled to 2n for each case provided 2 separate alleles.

"---": values were not counted due to limited samples.

\_^a^: Pc=0.416; \_ ^b^: Pc=0.494; \_ ^c^: Pc=0.364; \_ ^d^: Pc=0.013; \_ ^e^: Pc=0.221; \_ ^f^: Pc=0.000; \_ ^g^: Pc=0.026; \_ ^h^: Pc=0.04; \_ ^i^: Pc=0.04; \_ ^j^: Pc=0.128; \_ ^k^: Pc=0.000; \_ ^l^: Pc=0.016; \_ ^m^: Pc=0.000.

###### 

The Relationship between HLA-DRB1, DQB1 Polymorphisms and HPV16 E6 Variants.

  -------------------------------------------------------------------------------------------------------------------------------------
  Alleles   Healthy Controls\   HPV16-positive Cervical SCCs\   E6 Variants                                         
            2n=260\             2n=132\                                                                             
            (%)                 (%)                                                                                 
  --------- ------------------- ------------------------------- -------------------- --------- -------------------- -------------------
  DRB1                                                                                                              

  \*01      8(3.1)              1(0.8)                          0(0.0)               0(0.0)    0(0.0)               1(10.0)

  \*03      23(8.8)             9(6.8)                          4(15.4)              0(0.0)    5(6.0)               0(0.0)

  \*04      16(6.2)             10(7.6)                         2(7.7)               0(0.0)    8(9.5)               0(0.0)

  \*07      21(8.1)             8(6.1)                          1(3.8)               1(8.3)    6(7.1)               0(0.0)

  \*08      29(11.2)            23(17.4)                        [7(27.0)]{.ul}^a^    4(33.3)   11(13.1)             1(10.0)

  \*09      26(10.0)            19(14.4)                        0(0.0)               0(0.0)    [19(22.6)]{.ul}^b^   0(0.0)

  \*10      8(3.1)              2(1.5)                          0(0.0)               0(0.0)    2(2.4)               0(0.0)

  \*11      10(3.8)             13(9.8)                         3(11.5)              2(16.7)   5(6.0)               [3(30.0)]{.ul}^c^

  \*12      29(11.2)            18(13.6)                        6(23.1)              2(16.7)   6(7.1)               [4(40.0)]{.ul}^d^

  \*13      18(6.9)             0(0.0)                          0(0.0)               0(0.0)    0(0.0)               0(0.0)

  \*14      21(8.1)             7(5.3)                          2(7.7)               3(25.0)   2(2.4)               0(0.0)

  \*15      35(13.5)            16(12.1)                        0(0.0)               0(0.0)    16(19.0)             0(0.0)

  \*16      16(6.2)             6(4.5)                          1(3.8)               0(0.0)    4(4.8)               1(10.0)

  DQB1                                                                                                              

  \*02      55(21.2)            13(9.8)                         2(7.7)               2(16.7)   7(8.3)               2(20.0)

  \*0301    65(25.0)            27(20.5)                        [12(46.2)]{.ul}^e^   5(41.6)   7(8.3)               3(30.0)

  \*0302    21(8.1)             7(5.3)                          0(0.0)               0(0.0)    7(8.3)               0(0.0)

  \*0303    31(11.9)            23(17.4)                        3(11.5)              1(8.3)    [18(21.4)]{.ul}^f^   1(10.0)

  \*0306    8(3.1)              0(0.0)                          0(0.0)               0(0.0)    0(0.0)               0(0.0)

  \*04      21(8.1)             9(6.8)                          0(0.0)               0(0.0)    9(10.7)              0(0.0)

  \*0501    34(13.1)            23(17.4)                        7(27.0)              2(16.7)   12(14.3)             2(20.0)

  \*0601    25(9.6)             30(22.7)                        2(7.7)               2(16.7)   [24(28.6)]{.ul}^g^   2(20.0)
  -------------------------------------------------------------------------------------------------------------------------------------

Abbreviations: HPV, human papillomavirus; SCC, squamous cervical cancer

\*The valid sample size was calculated as 2n for each case provided 2 separate alleles.

\_ a: p=0.045, Pc=0.585; \_ b: p =0.003, Pc=0.039; \_ c: p =0.002, Pc=0.026; \_ d: p =0.025, Pc=0.325; \_ e: p =0.020, Pc=0.160; \_ f: p =0.030, Pc=0.240; \_ g: p =0.000, Pc=0.000.

###### 

The Distribution of DRB1, DQB1 Alleles Associated with the A4 lineage between Group I and Group II.

  ------------------------------------------------------------------------------------------
  Alleles   Group I 2n=58(%)   Group II\   *X*^2^   *p*               OR(95%CI)
                               2n=26(%)                               
  --------- ------------------ ----------- -------- ----------------- ----------------------
  DRB1                                                                

  \*03      0(0.0)             5(19.2%)    8.673    [0.003]{.ul}^a^   1.238(1.026\~1.494)

  \*04      7(12.1%)           1(3.8%)     0.616    0.433             3.431(0.400\~29.435)

  \*07      3(5.2%)            3(11.5%)    0.347    0.556             0.418(0.079\~2.227)

  \*08      11(19.0%)          0(0.0)      4.130    [0.042]{.ul}^b^   0.810(0.715\~0.918)

  \*09      13(22.4%)          6(23.1%)    0.005    0.946             0.963(0.320\~2.897)

  \*10      2(3.4%)            0(0.0)      \-       \-                \-

  \*11      2(3.4%)            3(11.5%)    0.903    0.342             0.274(0.043\~1.748)

  \*12      5(8.6%)            1(3.8%)     0.107    0.743             2.358(0.262\~21.266)

  \*14      2(3.4%)            0(0.0)      \-       \-                \-

  \*15      9(15.5%)           7(26.9%)    0.865    0.352             0.499(0.163\~1.529)

  \*16      4(6.9%)            0(0.0)      0.669    0.413             0.931(0.868\~0.999)

  DQB1                                                                

  \*02      2(3.4%)            5(19.2%)    3.970    [0.046]{.ul}^c^   0.150(0.027\~0.833)

  \*0301    4(6.9%)            3(11.5%)    0.081    0.776             0.568(0.418\~2.742)

  \*0302    4(6.9%)            3(11.5%)    0.081    0.776             0.568(0.418\~2.742)

  \*0303    11(19.0%)          7(26.9%)    0.675    0.411             0.635(0.214\~1.884)

  \*04      9(15.5%)           0(0.0)      3.042    0.081             \-

  \*0501    8(13.8%)           4(15.4%)    0.000    1.000             0.880(0.240\~3.231)

  \*0601    20(34.5%)          4(15.4%)    3.208    0.073             2.895(0.876\~9.563)
  ------------------------------------------------------------------------------------------

[^1]: \*Equal contributors and co-first authors.

[^2]: Competing Interests: The authors have declared that no competing interests exist.
