Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Analysis of mutations in the E6 oncogene of human papillomavirus 16 in cervical cancer isolates from Moroccan women

  • Zineb Qmichou1, 2,
  • Meriem Khyatti3,
  • Mohamed Berraho4,
  • My Mustapha Ennaji2,
  • Laila Benbacer1,
  • Chakib Nejjari4,
  • Noureddine Benjaafar5,
  • Abdellatif Benider6,
  • Mohammed Attaleb1 and
  • Mohammed El Mzibri1Email author
Contributed equally
BMC Infectious Diseases201313:378

https://doi.org/10.1186/1471-2334-13-378

Received: 3 May 2013

Accepted: 15 August 2013

Published: 16 August 2013

Abstract

Background

Worldwide, cervical cancer is the second most common cancer in women. High-risk human papillomavirus (HPV) play a crucial role in the etiology of cervical cancer and the most prevalent genotype is HPV16. HPV 16 intratypic variants have been reported to differ in their prevalence, biological and biochemical properties. The present study was designed to analyze and identify HPV type 16 E6 variants among patients with cervical cancer in Morocco.

Methods

A total of 103 HPV16 positive samples were isolated from 129 cervical cancer cases, and variant status was subsequently determined by DNA sequencing of the E6 gene.

Results

Isolates from patients were grouped into the European (E), African (Af) and North-American (NA1) phylogenetic clusters with a high prevalence of E lineage (58.3%). The Af and NA1 variants were detected in 31.1% and 11.6% of the HPV16 positive specimens, respectively, whereas, only 3% of cases were prototype E350T. No European-Asian (EA), Asian (As) or Asian-American (AA) variants were observed in our HPV16-positive specimens. At the amino acid level, the most prevalent non-synonymous variants were L83V (T350G), H78Y (C335T), E113D (A442C), Q14D (C143G/G145T) and R10I (G132T), and were observed respectively in 65%, 41.8%, 38.8%, 30.1% and 23.3% of total samples.

Moreover, HPV16 European variants were mostly identified in younger women at early clinical diagnosis stages. Whereas, HPV16 Af variants were most likely associated with cervical cancer development in older women with pronounced aggressiveness.

Conclusion

This study suggests a predominance of E lineage strains among Moroccan HPV 16 isolates and raises the possibility that HPV16 variants have a preferential role in progression to malignancy and could be associated with the more aggressive nature of cervical cancer.

Keywords

Cervical cancer Human papillomavirus 16 Variants E6 Morocco

Background

Cervical cancer is the second most predominant cancer worldwide, and more than 85% of the global burden occurs in developing countries, where it accounts for 13% of all female cancer [1]. In Morocco, as it is the case in the other North African countries, cervical cancer is the second most common cancer among women and its incidence is the highest in this region with an age-standardized incidence rate (ASR) of 13.5 per 100 000 women [2]. Persistent infection with high risk Human papillomavirus (HPV) is the main aetiological factor in the development of cervical cancer. Globally, HPV16 is the most frequent HPV type found in cervical cancer and is identified in approximately 65% of cases [3]. In Morocco, our previous reports have shown that the HPV16 is also the most frequent type of HPV found in cervical cancer patients [46].

Given that the prevalence of cervical cancer varies in different regions and countries, a number of studies have addressed the possible association of HPV16 variant status with the risk for progression to malignancy. HPV variants are defined as isolates with primary DNA sequence differences that total no more than 2% of the L1 open reading frame (ORF) of the prototype sequence [7]. The sequence variations of the HPV16E6 ORF have been found to correctly classify the HPV16-variants [8]. Therefore, most studies on HPV16 variants have focused on E6, the major transforming protein that inhibit apoptosis and promote proliferation and their association to cancer aggressiveness [911]. Worldwide, variants of the HPV16 genome in cervical carcinoma, precancerous lesions and normal or inflammatory cervical cases, from different geographical areas have been classified into phylogenetic lineages as European (E), Asian (As), Asian-American (AA), African -1 and -2 (Af1 and Af2), and North American1 (NA1) [8, 12, 13]. Overall, multiple studies have shown that HPV16 variants differ in risk for progression to high grade intraepithelial lesions [1416], in their association with the development of cervical cancer [15, 1719], viral persistence [17, 2026] and the frequency of recurrence of cervical disease [15].

Many authors confirm that distribution of HPV variants is related to geographic or race distribution [27]. However, although numerous studies on HPV16 variants in cervical cancer have been carried out so far worldwide, only one HPV variant research in patients from Tunisia has been performed in North Africa [28]. To our knowledge, nothing is known concerning HPV variants in patients from Morocco, where there is a significantly high occurrence of this tumor. In the present study entirely done in Morocco, we have examined whether naturally occurring sequence variations in the HPV16 E6 gene exist in cervical cancer from Moroccan subjects with HPV 16 infection. The identification of more persistent or oncogenic HPV16 variants would be of great interest in the global strategy in the development and implementation of improved therapeutic vaccines against HPV.

Furthermore, HPV sequence variations studies can be used in epidemiological studies as a marker to track the spread of the virus in contact networks through populations.

Methods

Clinical specimens

One hundred twenty nine fresh frozen uterine cervix biopsies were recruited for the variant analysis in the present study. The tumor samples were from biopsies collected for the routine diagnosis of cervical cancer from Gynaecology department of the National Institute of Oncology in Rabat and Oncology Centre Mohammed VI pour le Traitement des Cancers, CHU Ibn Rochd in Casablanca, that receive patients from all regions of Morocco. The mean age of patients was 52 ± 12 years, with extreme ages at 28 and 80. Each sample was divided into two portions: one portion was put in neutral buffered formalin and processed for routine histopathological examination; the other portion was stored at −80°C immediately after surgical removal, for HPV genotyping. Protocols were approved by Ethical Review Committees of the Faculty of Medicine of Fez and the Hassan-II University Hospital, and written informed consent was obtained from each study subject.

Cell lines

CaSki and SiHa cervical cancer cell lines were obtained from the American Type Culture Collection (ATCC) and used as positive controls. These cell lines were maintained in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal calf serum (FCS).

Pathology

Samples were fixed in 10% buffered-formaldehyde. Pathology diagnosis was made using the routine hematoxylin-eosin stain on 5 μm paraffin sections. Histopathological results of cervical cancer lesions were classified according to the International Federation of Gynaecology and Obstetrics (FIGO).

DNA extraction

Genomic DNA was extracted from fresh frozen tissue and cervical cancer cell lines using a standard technique of digestion with proteinase K in the presence of sodium dodecyl sulfate (SDS) at 37°C overnight, followed by phenol/chloroform purification [29]. The tissue DNA was precipitated with 2/5 volumes of 7.4 M ammonium acetate and 2 volumes of 100% ethanol, followed by incubation at −20°C and centrifugation at top speed (13000 relative centrifugal force). DNA was then resuspended in 25–30 μL of sterile distilled water and stored at −20°C until use. In order to evaluate the efficiency of DNA extraction, all samples were amplified by PCR using PC04 and GH20 primers specific for human β-globin gene (Table 1).
Table 1

List of primers used for PCR amplification and DNA sequencing

 

Primer

Fragment generated size

Sequence

Tm (°C)

β-globin

PC04

268 bp

5′-CAACTTCATCCACGTTCACC-3′

54.5

GH20

5′-GAAGAGCCAAGGACAGGTAC-3′

HPV detection & genotyping

MY09

450 bp

5′-CGT CCM ARR GGA WAC TGATC- 3′

55

MY11

5′-GCM CAG GGW CAT AAY AAT GG-3′

GP5+

150 pb

5′-TTTGTTACTGTGGTAGATACTAC-3′

55

GP6+

5′-CTTATACTAAATGTCAAATAAAAA-3′

HPV16 E6 variants

F

524 bp

5′-CGAAACCGGTTAGTATAA-3′

55

R

5′-GTATCTCCATGCATGATT-3′

F: Forward primer; R: Reverse primer

HPV genotyping

HPV detection and genotyping of cervical cancer cases were performed using a nested PCR, with the consensus MY09/11 and GP5+/6+ primers (Table 1), and direct DNA sequencing [30]. All specimens were first subjected to PCR amplification with HPV consensus primers MY09/11 [31]. For nested PCR, 2 μL of the MY09/11 PCR products was used as a template for PCR amplification using GP5+/6+ general primers [30, 32]. For every reaction, a negative control, in which DNA template was omitted from the amplification mixture, and a positive control, with DNA extracted from SiHa and Caski cell lines, were included. PCR products were analyzed by electrophoresis on 2% agarose gels followed by staining with ethidium bromide (10 mg/mL).

For DNA sequencing, the PCR products were purified by the ExoSaP-IT clean up system (USB, USA) and sequenced directly on an ABI 3130XL DNA analyzer (Applied Biosystems, Foster city, CA, USA), using GP6+ primer as the sequencing primer and BigDye® Terminator v3.1 Cycle Sequencing Kit (AppliedBiosystems, Foster city, CA, USA), according to manufacturer’s protocol. For each sample, PCR amplification and DNA sequencing were performed twice. Nucleotides sequences were aligned and compared with those of known HPV types available through GenBank by using the online BLAST 2.0 software server (http://www.ncbi.nih.gov/blast). In accordance with established guidelines, a nucleotide sequence was assigned to an HPV type if it corresponded with a known HPV genotype by >90% [33].

Determination of HPV16 E6 variants

Samples with a single HPV16 genotype were selected for E6 variants investigation. Determination of HPV16 E6 variants was performed by E6 specific PCR, using primers flanking the encoding region of HPV16 E6 gene (nt 52–575) (Table 1), and DNA sequencing.

PCR amplifications were performed according to Lizano et al. [34]. For every reaction, positive controls using DNA extracted from SiHa and Caski cell lines, as well as a negative control without template DNA, were included. PCR products were tested on an ethidium bromide stained 2% agarose gel. PCR products were purified and sequenced with the same primers as those used for the amplification [34]. The HPV16 E6 sequences obtained were aligned with those of European prototypes of HPV16 sequence (HPV16-R) (GenBank accession numbers: K02718, NC_001526), available through the GenBank database (NCBI, National Institute of health, Bethesda, MD, USA), using the BLAST 2.0 software server (http://www.ncbi.nih.gov/BLAST/), MEGA v4.0 and SeqScape v2.6 (Applied Biosystems, Foster city, USA) programs.

Statistical analysis

The HPV16 variants according to other clinical features were compared using a Fisher’s exact test using the MedCalc statistical software (http://www.medcalc.org). The statistical relationship was considered as significant if the derived p value was < 0.05.

Results

Histopathological data

The pathological analysis revealed the predominance of Squamous Cell Carcinoma (SCC) that represents 96.89% of cases (125/129). The remaining 4 cases (3.1%) were diagnosed as adenocarcinoma. According to stage, to the most important percentages were found in stage II and III representing 48.8% (61 cases) and 38.4% (48 cases), respectively. Moreover, 14 cases were classified as stage I (11.2%), whereas only 2 cases (1.6%) were classified as stage IV.

Subtype distribution of HPV DNA

The presence of amplifiable DNA was confirmed for all 129 cases by PCR based-technique using primers for a fragment of β-globin gene and therefore all DNA samples were adequate for further analysis. Molecular detection of HPV DNA, using nested PCR amplification of a conserved region of the HPV L1 gene DNA with the consensus MY09/11 and GP5+/6+primers, revealed the presence of viral DNA in 91.47% (115 of 129 SCC and 3 of 4 adenocarcinoma).

In this study, we have identified 5 carcinogenic HPV genotypes: 16, 18, 31, 33 and 35. The most frequently found was HPV16, detected in 87.29% of all positive samples (103 cases), followed by HPV 18, present in 5.08% of samples (6 cases), HPV 31 in 2.54% of samples (3 cases), HPV 33 in 1.69% of samples (2 cases) and HPV 35 in 1 specimen (0.85% of cases). Undetermined HPV genotypes were identified in 3SCC cases.

Sequence variations of HPV16 E6 gene

The 103 HPV16 single infected specimens were selected for E6 intratypic variation analysis. Overall, DNA sequences were compared to the HPV16 reference sequence (Accession Number: K02718, NC_001526). Accordingly, HPV16-positive specimens were classified into E, Af1, Af2 and NA branches (Table 2).
Table 2

Characterization of HPV16 E6 variants and their lineage classification

 

HPV 16 variants

N

T109C

A131G

G132T

C143G

G145T

C285G

T286A

A289G

T295G

T325C

C335T

T350G

A403G

A442C

Predicted substitution

 

HPV16 Reference

 

T

A

G

C

G

C

T

A

T

T

C

T

A

A

 

Lineage

Sub-lineage

Class/sub-class

                

EUR

Ep

E-T350

3

-

-

-

-

-

-

-

-

-

-

-

-

-

-

 

E-350G

E-G350

16

-

-

-

-

-

-

-

-

-

-

-

G

-

-

L83V

E-C442/G350

40

-

-

-

-

-

-

-

-

-

-

-

G

-

C

L83V/E113D

E-G350/T145

1

-

-

-

-

T

-

-

-

-

-

-

G

-

-

Q14H/L83V

AFR

Af1

Af1-d/G295

5

-

G

-

G

T

-

a

g

G

-

T

G

 

-

R10G/Q14D/D64E/H78Y/L83V

Af1-d/G295/C325

2

-

G

-

G

T

-

a

g

G

c

T

G

-

-

R10G/Q14D/D64E/H78Y/L83V

Af2

Af2-a/C109/G403

13

c

-

T

G

T

-

a

g

-

-

T

-

g

-

R10I/Q14D/H78Y

Af2-a/G285

1

-

-

T

G

T

G

a

g

-

-

T

-

g

-

R10I/Q14D/A61G/H78Y

Af2a/r

10

-

-

T

G

T

-

a

g

-

-

T

-

g

-

R10I/Q14D/H78Y

NA

NA1

NA1

9

-

-

-

-

-

-

a

g

-

-

T

-

g

-

H78Y

NA1

2

-

-

-

-

-

-

a

g

-

-

T

G

-

-

H78Y/L83V

NA1-b/r

1

-

-

-

-

T

-

a

g

-

-

T

G

-

-

Q14D/H78Y/L83V

N*

 

13

7

24

31

33

1

43

43

7

2

43

67

33

40

 

Sequence alterations relative to the E6 open reading frame (ORF) of the reference HPV16 sequence [35]. For subclasses, the number and the letter represent the nucleotide change at that position. Lower case represents a nucleotide change at that position without an amino acid change. Capital letters indicate variants with an amino acid change. In the ‘Predicted substitution’ column, the letter preceding the amino acid position refers to the reference HPV16 sequence and the letter after it refer to the substitution. Dashes indicate no variation. * Total cases harboring the point mutation. EUR: European lineage, AFR: African lineage, NA: Nord American lineage.

By comparison with the prototype sequence of HPV16 E6, it was found that only 2.9% of samples (3 out 103) revealed the prototype sequence, while 97.1% of samples (100 out 103) showed at least one specific nucleotide variation in the HPV16 E6 region (Table 2). Moreover, in all analyzed specimens, there was no evidence of premature stop codons, bases insertions or deletions.

Overall, 9 variants with predicted E6 amino acid changes were identified (Table 2). The most prevalent non-synonymous variants were L83V (T350G), H78Y (C335T), E113D (A442C) and Q14D (C143G/G145T) and were observed respectively in 65%, 41.8%, 38.8% and 30.1% of total samples. The other non-synonymous variants were R10I (G132T), R10G (A131G), D64E (T295G), Q14H (G145T) and A61G (C285G) and were seen in 23.3%, 6.8%, 6.8%, 1.9% and 1% of cases respectively.

A total of 12 different intratypic variants were identified among the 103 HPV16 specimens isolated from Moroccan women. These intratypic variants are grouped into three different variant classes, including European (E), African (Af) and the North-American (NA1) variant classes. The European variants HPV16 were the most frequently detected isolates in our population and represent 58.3% of total HPV16 isolates (n=60), followed by Af characterized in 30.1% of cases (n=31). Whereas the NA variants were found only in 11.6% of cases (n=12) (Table 2). In the European class, two sub-classes prevail, E-350G/442C identified in 40 specimens (38.8%) and E-350G identified in 16 specimens (15.5%). The common no-E variants are represented by Af2, NA1 and Af1 sub-classes, which are identified in 24 (23.3%), 12 (11.6%) and 7 (6.8%) cases, respectively. Asian (As) and Asian-American (AA) HPV16 variants were not detected in the present study.

The distribution of the different E6 variations in the HPV16 positive samples with respect to FIGO staging of cervical neoplasia is reported in Table 3. E class variants, the most common intratypes, were mainly associated with cancers detected in early clinical stages. The majority of cervical neoplasia with HPV16 E variant lesions presented at stage I and II with 75% (9/12) and 60% (30/50) respectively. Conversely, Af variants seem to be most likely detected in advanced cancer stages, there were detected in 37.5% (15/40) of cases at stage III and only in 16.67% (2/12) and 28% (14/50) of cases at stages I and II respectively. NA1 variants were detected with similar frequencies in early as well as advanced stages. However, no significant statistical difference between HPV variant and stage (early I and II, advanced III and IV) of cancer was detected (p > 0.05).
Table 3

Distribution of HPV16 E6 variants according to clinical stage

Clinical stage

N

E

Af

NA

E350G

E350G/442C

E-G350/G145

Ep

Af1

Af2

NA1

I

12

2

7

-

-

1

1

1

9

2

II

50

10

18

1

1

2

12

6

30

14

III

40

4

14

-

2

4

11

5

20

15

IV

1

-

1

-

-

-

-

-

1

 

Subtotal

16

40

1

3

7

24

12

Total

103

60

31

12

The distribution of HPV16 E6 variants according to patients’ age was also assessed and is reported in Table 4. Results showed an overrepresentation of Af class in the age group over 45 years old whereas E and NA variants were most identified in the age group under 45 years. Moreover, among women under 45 years, E350G/442C is the most predominant variant with 48.15% of cases (13/27). However, in women over 45 years, two variants prevail: E350G/442C and Af2 which are present in 35.53% (27/76) and 25% (19/76) respectively, statistical analyses were significant for only HPV16 E and Af variants related to patients’ age (p < 0.0001).
Table 4

Distribution of HPV16 E6 variants according to patients’ age

Patients’ age

N

E

Af

NA

E350G

E350G/442C

E-G350/G145

Ep

Af1

Af2

NA1

≤ 45 years

27

4

13

-

-

-

5

5

17

5

5

> 45 years

76

12

27

1

3

7

19

7

43

26

7

Subtotal

16

40

1

3

7

24

12

Total

103

60

31

12

Correlation between HPV16 E6 variants and malignant phenotype showed that E variants are mostly associated with high degree of differentiation, moderately and well differentiated carcinoma, whereas Af variants are more frequently found in poorly differentiated carcinoma (Table 5).
Table 5

Correlation between malignant phenotype and HPV16 E6 variants

Differentiation

N

E

Af

NA

E350G

E350G/442C

E-G350/G145

Ep

Af1

Af2

NA1

Poorly

17

3

4

-

-

1

7

2

7

8

Moderately

39

5

18

1

2

4

7

2

26

11

Well

31

6

13

-

-

1

6

5

19

7

Unknown

16

2

5

-

1

1

4

3

8

5

    

Subtotal

16

40

1

3

7

24

12

Total

103

60

31

12

Discussion

The present study, entirely realized in Morocco, highlighted the distribution patterns of intratypic variants of HPV16 among women with cervical cancer. A total of 129 samples of cervical cancer cases were collected. Molecular detection of HPV, using nested PCR amplification of a conserved region of the HPV L1 gene with the consensus MY09/11 and GP5+/6+ primers, revealed the presence of viral DNA in 91.47% of cases (118/129). We identified 5 carcinogenic HPV genotypes (16, 18, 31, 33 and 35). Globally, HPV prevalence and distribution is in concordance with previously reported data in Morocco [4, 6], and tally with the overall world distribution of HPV types in cervical cancer [36, 37]. Moreover, HPV16 was the HPV genotype most frequently found in our study population, with 87.29% of all positive samples, which is in agreement with the worldwide reported data [3, 18, 37, 38]. However, previous studies conducted in Morocco have reported lower prevalence of HPV16.

Lalaoui et al. [4] have found HPV 16 in 49% of HPV positive cases, whereas in the study of Meftah El Khair et al.[31], HPV16 was found in 71% of HPV positive cases, including mono-infected (37%) and co-infected (34%) cases. This difference may be related to the sampling bias but also could be due to the molecular technique used for the HPV genotyping. In the previous studies, HPV genotyping was realized by combining consensus PCR and dot blot hybridization with specific probes. However, in our study, the HPV genotyping was performed by DNA sequencing of the hyper-variable region in L1 fragment which is reported to be more accurate and give a high sensitivity [30, 33].

To examine the prevalence of naturally occurring sequence variations in the HPV16 E6 gene in cervical cancer from Moroccan subjects, all HPV16 positive samples (n=103) were subjected to intratypic variant analysis by nucleotide sequencing of the E6 gene and comparison with the prototype sequence of HPV16 E6.

Of all HPV16 positive samples, only 3 cases exhibited a prototype E6 gene sequence and 97.1% of cases (100/103) showed at least one specific nucleotide variation in the HPV16 E6 region. This finding is in agreement with published results reporting that approximately 90% of the E6 genes in invasive cervical cancer contained variations [39, 40].

It is widely accepted that the distribution of HPV variants is related to geographic or race distribution [21, 27] and therefore, owing to the geographic location of Morocco in the North Africa, we expected the predominance of African variants. Unexpectedly, however, the E6 gene sequencing analyses highlighted that the most frequently HPV16 variants detected in our population belonged to the European lineage (58.3%), whereas the African lineage was detected only in 30.1% of cases. The NA1 lineage was less represented and limited to 11.6% of cases. The frequent circulation of the European lineage in Morocco most likely reflects the immigration flow that characterized this region for a long time in the past. Morocco, due to its geographic position, is recognized as a border region between Europe and Africa, and therefore was a corner stone in all exchanges between the two continents. Authors have already reported that HPV strains could be introduced as a result of migration leading to the spread of the molecular variants of HPV from one continent to another [41]. Krennhrubec et al.[28] have reported the same results in Tunisia, a neighbor country in North Africa, with 64% of HPV16 infections being of European lineage.

HPV16 European variant is a highly spread variant around the world, a meta-analysis conducted on HPV16 positive tumors samples collected globally from different ethnical group reported that European lineage was the most predominant lineage and was well distributed among the different geographical regions. This widespread presence suggests biological relevance of this variant [42].

A large number of HPV16 variants circulating in our population belong to African branches (30.1% of cases). The association of Af variants with an increased risk of cervical cancer development is controversial. Some authors have reported a greater association between infection with HPV16 Af variants and an increased risk of cervical carcinogenesis [15, 43, 44]. However, Tu et al.[45] have reported no evidence that Af variants could be associated with an increased risk of neoplasia. Within the Af lineage, the two main subgroups, Af1 and Af2 were detected in 7 (6.8%) and 24 (23.3%) cases respectively. These results are in contrast to previous studies that reported the predominance of Af1 variants, among the Af lineage, in North Africa [12, 42].

Interestingly, we have detected NA1 lineage in 11.6% of our patients’ population. The occurrence of NA1 classes, which is a derivative of the AA and Af2 variants [12], is difficult to explain. The NA1 lineage, identified in 1990s [12], was recently found to be particularly frequent in samples from North Africa [42]. We believe that NA1 lineage could be introduced recently in Morocco and spread in the population. That could explain the lower frequency of this lineage in the study population.

According to clinico-pathological data, HPV16 European variants (350G alone or along with other mutations) were mostly identified in younger women at early clinical stages (I and II) at diagnosis. However, most of the cancer patients, harboring the HPV16 Af variants, were older (over 45 years) and exhibited a pronounced aggressiveness of the disease, as the majority of cases were at advanced stage (III) with poorly differentiated carcinoma. These results raise the possibility that HPV16 variants have a preferential role in progression to malignancy and could be associated with the more aggressive nature of cervical cancer. Similar results were obtained in Tunisia reporting a higher proportion of non-European HPV16 variants in advanced stage (III). In contrast, the non-E variants were more frequently detected at earlier average age at diagnosis than European variants [28]. This difference could not be due to genetic diversity of the host, as the two populations are very closely related, but could be due to a sampling bias and biological behavior of the host.

Most nucleotide changes reported in our study have been previously described and some of them are of particular interest and may have considerable biological impact [10, 46]. The most prevalent non-synonymous variants were L83V (T350G), H78Y (C335T), E113D (A442C), Q14D (C143G/G145T) and R10I (G132T), and were observed respectively in 65%, 41.8%, 37.9%, 30.1% and 23.3%of total samples. Against the above background of higher pathogenicity associated with the L83V variant, the finding of our study reiterates the same. Indeed, the most frequently observed mutation in our study was L83V which is detected in 65% of cases (67/103), either alone (15.5% of cases) or along with other mutations (49.5% of cases). This mutation is the most commonly found in HPV16 E6 gene from cervical invasive cancers [44], and has been linked to an increased risk for cervical disease progression [20]. Many authors addressed a geographic variation of the oncogenicity of this variant [22, 39, 47]. In some European populations, HPV16 E350G variant was strongly associated with the oncogenicity and persistency of HPV16 infection [20, 39]. Moreover, HPV16 E350G variant was found to display more efficient degradation of Bax (protein regulated by p53) and binding to E6-BP (E6 binding protein) and decreased binding to human discs large protein (hDlg) [48]. This variant appears to have surpassed the E6 prototype in enhancing the mitogen activated protein kinase (MAPK) signaling evolved in oncogenic ras-mediated transformation and in the cooperative transformation with the deregulated Notch 1 pathway [49]. Also, this L83V polymorphism appears to interact with natural human variations in the p53 genes (in particular codon 72 polymorphism) to confer differences in cervical cancer risk [50].

Additionally, E6 variants R10G/L83V and Q14H/H78Y/L83V which are frequently detected in our study, were more prone to undergo cell-detachment-induced apoptosis than E6 prototype in a model of human normal immortalized keratinocytes (NIKS) [51] and has a greater in vitro ability to suppress keratinocyte differentiation responses and to induce p53 degradation than the prototype [52].

Interestingly, the E113D (A442C) was detected in 37.9% of samples and was exclusively found in the European variant (E350G/442C). This point mutation is largely reported in Asian populations and is implied in invasive cervical carcinoma [5355].

The majorities of the E6 amino acid switches in our samples were located at amino acid residues10, 14, 78 and 83 and were shown to be under positive selective pressure [56, 57]. This selection pressure, driven by protein–protein polymorphic interactions, could facilitate the life cycle and pathogenicity of HPV associated disease and might be under immune selection [57]. Most recently, Bayesian analysis identified several important amino acid positions that may be driving adaptive selection in the HPV 16 population, including R10G, L83V, and E113D in the E6 gene [48]. Therefore, the positive selection at these codons might be a contributing factor responsible for the phenotypic differences in carcinogenesis and immunogenicity among cervical cancers in Morocco. Indeed, it has been suggested that women with three distinct HLA class I alleles, namely HLA-B*44, HLA-B*51, or HLA-B*57 who were infected with the HPV16 E6 variant L83V had an approximately four to five fold increased risk for cancer compared with controls [49, 58]. Additionally, in another study, the HPV 16 E6 variant G131 (R10G) was demonstrated to alter a B*07 binding epitope such that it may influence immune recognition by cytotoxic T lymphocytes [59]. Also, the immunologic relevance of the HPV16 E6 N-terminal region and variant positions E6 amino acids 10 and 14 is supported by the demonstration of an endogenously processed HLA A*0201-restricted E6 peptide (E6 amino acids 11 to 19) as well as of an overlapping HLA B-7-restricted E6 peptide (E6 amino acids 8 to 15) in this region [59, 60].

Therefore, the observed amino acid variations in HPV16 E6 detected in our study may impact on the viral life cycle, which is tightly linked to the differentiation program of the cell [57].

Conclusion

The introduction of molecular sub-typing of HPV allowed better genotyping and therefore a better characterization of strains circulating in Morocco. Thereby, the identification of HPV16 E, Af and NA1 variants suggests a different sexual mixing behavior. Their circulation in our population may involve a change in the biological properties and a change in the transforming potential, increasing the risk of cervical cancer development. Currently, Morocco is a destination of many sub-Saharan migrants looking for better life, and is also a gateway to Europe, that warrants the need for monitoring a wider range of HPV variants involved in cervical cancer in Morocco. Thus, in order to better understand the dynamics of the disease and to optimize current control programs and policies for a better survey of cervical cancer, it is useful to develop routine epidemiological surveillance, including molecular sub-typing.

Authors’ information

Senior authors: Mohammed Attaleb and Mohammed El Mzibri.

Notes

Declarations

Acknowledgments

Authors would like to thank the Plateforme de Génomique Fonctionnelle, Unité d’Appui Technique à la Recherche Scientifique, Centre National de la Recherche Scientifique et Technique (Morocco) for their excellent technical cooperation. This study was partially supported by the EU FP7/2007-2013 grant 260715.

Authors’ Affiliations

(1)
Unité de Biologie et Recherche Médicale, Centre National de l'Energie, des Sciences et des Techniques Nucléaires, (CNESTEN)
(2)
Laboratoiore de Microbiologie, Hygiène et Virologie, FST Mohammedia
(3)
Laboratoire d’Oncovirologie. Institut Pasteur du Maroc
(4)
Laboratoire d'Epidémiologie, Recherche clinique et Santé Communautaire, Faculté de Médecine et de Pharmacie
(5)
Service de radiothérapie, Institut National d'Oncologie de Rabat
(6)
Centre Mohammed VI pour le Traitement des Cancers

References

  1. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM: Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010, 127: 2893-2917. 10.1002/ijc.25516.View ArticlePubMedGoogle Scholar
  2. Benider A, Bennani Othmani M, Harif M, Karkouri M, Quessar A, Sahraoui S, Sqali S: Registre des cancers de la région du grand Casablanca, année 2004. Edition. 2007, 73-Google Scholar
  3. Smith JS, Lindsay L, Hoots B, Keys J, Franceschi S, Winer R, Clifford GM: Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: a meta-analysis update. Int J Cancer. 2007, 121: 621-632. 10.1002/ijc.22527.View ArticlePubMedGoogle Scholar
  4. Lalaoui K, El Mzibri M, Amrani M, Belabbas MA, Lazo PA: Human papillomavirus DNA in cervical lesions from Morocco and its implications for cancer control. Clin Microbiol Infect. 2003, 9: 144-148. 10.1046/j.1469-0691.2003.00494.x.View ArticlePubMedGoogle Scholar
  5. El Mzibri M, Amrani M, Khair MM, Attaleb M, Jaddi H, Ennaji MM, Khan B, Belabbas MA: Prevalence of human papilloma virus infection in malignant lesions of the uterine cervix in Morocco. World J Nucl Med. 2009, 8: 192-196.Google Scholar
  6. Meftah El khair M, Ennaji MM, El kebbaj R, Ait Mhand R, Attaleb M, El Mzibri M: p53 Codon 72 polymorphism and risk of cervical carcinoma in Moroccan women. Med Oncol. 2010, 27: 861-866. 10.1007/s12032-009-9297-6.View ArticlePubMedGoogle Scholar
  7. de Villiers E, Fauquet C, Broker TR, Bernard HU, zur Hausen H: Classification of papillomaviruses. Virology. 2004, 324: 17-27. 10.1016/j.virol.2004.03.033.View ArticlePubMedGoogle Scholar
  8. Yamada T, Wheeler CM, Halpern AL, Stewart AC, Hildesheim A, Jenison SA: Human papillomavirus type 16 variant lineages in United States populations characterized by nucleotide sequence analysis of the E6, L2, and L1 coding segments. J Virol. 1995, 69: 7743-7753.PubMedPubMed CentralGoogle Scholar
  9. Zur Hausen H: Papillomaviruses and cancer: from basic studies to clinical application. Nat Rev Cancer. 2002, 2: 342-350. 10.1038/nrc798.View ArticlePubMedGoogle Scholar
  10. Lizano M, Berumen J, Garcia-Carranca A: HPV-related carcinogenesis: basic concepts, viral types and variants. Arch Med Res. 2009, 40: 428-434. 10.1016/j.arcmed.2009.06.001.View ArticlePubMedGoogle Scholar
  11. Kim MK, Kim HS, Kim SH, Oh JM, Han JY, Lim JM, Juhnn YS, Song YS: Human papillomavirus type 16 E5 oncoprotein as a new target for cervical cancer treatment. Biochem Pharmacol. 2010, 80: 1930-1935. 10.1016/j.bcp.2010.07.013.View ArticlePubMedGoogle Scholar
  12. Yamada T, Manos MM, Peto J, Greer CE, Munoz N, Bosch FX, Wheeler CM: Human papillomavirus type 16 sequence variation in cervical cancers: a worldwide perspective. J Virol Methods. 1997, 71: 2463-2472.Google Scholar
  13. Huertas-Salgado A, Martin-Gamez DC, Moreno P, Murillo R, Bravo MM, Villa L, Molano M: E6 Molecular variants of human papillomavirus (HPV) type 16: an updated and unified criterion for clustering and nomenclature. Virology. 2011, 410: 201-215. 10.1016/j.virol.2010.10.039.View ArticlePubMedGoogle Scholar
  14. Hildesheim A, Schiffman M, Bromley C, Wacholder S, Herrero R, Rodriguez A, Bratti MC, Sherman ME, Scarpidis U, Lin QQ, et al: Human papillomavirus type 16 variants and risk of cervical cancer. J Natl Cancer Inst. 2001, 93: 315-318. 10.1093/jnci/93.4.315.View ArticlePubMedGoogle Scholar
  15. Xi LF, Koutsky LA, Hildesheim A, Galloway DA, Wheeler CM, Winer RL, Ho J, Kiviat NB: Risk for high-grade cervical intraepithelial neoplasia associated with variants of human papillomavirus types 16 and 18. Cancer Epidemiol Biomarkers Prev. 2007, 16: 4-10. 10.1158/1055-9965.EPI-06-0670.View ArticlePubMedGoogle Scholar
  16. Villa LL, Sichero L, Rahal P, Caballero O, Ferenczy A, Rohan T, Franco EL: Molecular variants of human papillomavirus types 16 and 18 preferentially associated with cervical neoplasia. J Gen Virol. 2000, 81: 2959-2968.View ArticlePubMedGoogle Scholar
  17. Berumen J, Ordonez RM, Lazcano E, Salmeron J, Galvan SC, Estrada RA, Yunes E, Garcia-Carranca A, Gonzalez-Lira G, Madrigal-de la Campa A: Asian-American variants of human papillomavirus 16 and risk for cervical cancer: a case–control study. J Natl Cancer Inst. 2001, 93: 1325-1330. 10.1093/jnci/93.17.1325.View ArticlePubMedGoogle Scholar
  18. Li N, Franceschi S, Howell-Jones R, Snijders PJ, Clifford GM: Human papillomavirus type distribution in 30,848 invasive cervical cancers worldwide: variation by geographical region, histological type and year of publication. Int J Cancer. 2010, 128: 927-935.View ArticleGoogle Scholar
  19. Kammer C, Tommasino M, Syrjanen S, Delius H, Hebling U, Warthorst U, Pfister H, Zehbe I: Variants of the long control region and the E6 oncogene in European human papillomavirus type 16 isolates: implications for cervical disease. Br J Cancer. 2002, 86: 269-273. 10.1038/sj.bjc.6600024.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Grodzki M, Besson G, Clavel C, Arslan A, Franceschi S, Birembaut P, Tommasino M, Zehbe I: Increased risk for cervical disease progression of French women infected with the human papillomavirus type 16 E6-350G variant. Cancer Epidemiol Biomarkers Prev. 2006, 15: 820-822. 10.1158/1055-9965.EPI-05-0864.View ArticlePubMedGoogle Scholar
  21. Xi LF, Kiviat NB, Hildesheim A, Galloway DA, Wheeler CM, Ho J, Koutsky LA: Human papillomavirus type 16 and 18 variants: racerelated distribution and persistence. J Natl Cancer Inst. 2006, 98: 1045-1052. 10.1093/jnci/djj297.View ArticlePubMedGoogle Scholar
  22. Lee K, Magalhaes I, Clavel C, Briolat J, Birembaut P, Tommasino M, Zehbe I: Human papillomavirus 16 E6, L1, L2 and E2 gene variants in cervical lesion progression. Virus Res. 2008, 131: 106-110. 10.1016/j.virusres.2007.08.003.View ArticlePubMedGoogle Scholar
  23. Burk RD, Terai M, Gravitt PE, Brinton LA, Kurman RJ, Barnes WA, Greenberg MD, Hadjimichael OC, Fu L, McGowan L, et al: Distribution of human papillomavirus types 16 and 18 variants in squamous cell carcinomas and adenocarcinomas of the cervix. Cancer Res. 2003, 63: 7215-7220.PubMedGoogle Scholar
  24. Zuna RE, Moore WE, Shanesmith RP, Dunn ST, Wang SS, Schiffman M, Blakey GL, Teel T: Association of HPV16 E6 variants with diagnostic severity in cervical cytology samples of 354 women in a US population. Int J Cancer. 2009, 125: 2609-2613. 10.1002/ijc.24706.View ArticlePubMedPubMed CentralGoogle Scholar
  25. Quint KD, de Koning MN, van Doorn LJ, Quint WG, Pirog EC: HPV genotyping and HPV16 variant analysis in glandular and squamous neoplastic lesions of the uterine cervix. Gynecol Oncol. 2010, 117: 297-301. 10.1016/j.ygyno.2010.02.003.View ArticlePubMedGoogle Scholar
  26. Schiffman M, Rodriguez AC, Chen Z, Wacholder S, Herrero R, Hildesheim A, Desalle R, Befano B, Yu K, Safaeian M, et al: A population-based prospective study of carcinogenic human papillomavirus variant lineages, viral persistence, and cervical neoplasia. Cancer Res. 2010, 70: 3159-3169. 10.1158/0008-5472.CAN-09-4179.View ArticlePubMedPubMed CentralGoogle Scholar
  27. de Araujo Souza PS, Sichero L, Macia PC: HPV variants and HLA polymorphisms: the role of variability on the risk of cervical cancer. Future Oncol. 2009, 5: 359-370. 10.2217/fon.09.8.View ArticlePubMedGoogle Scholar
  28. KrennHrubec K, Mrad K, Sriha B, Ben Ayed F, Bottalico DM, Ostolaza J, Smith B, Tchaikovska T, Soliman AS, Burk RD: HPV types and variants among cervical cancer tumors in three regions of Tunisia. J Med Virol. 2011, 83: 651-657. 10.1002/jmv.22011.View ArticlePubMedPubMed CentralGoogle Scholar
  29. Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory manual 2nd ed. NY’ Cold Spring Harbor Laboratory. 1989, 1659-PressGoogle Scholar
  30. Lee SH, Vigliotti VS, Vigliotti JS, Pappu S: Routine human papillomavirus genotyping by DNA sequencing in community hospital laboratories. Infect Agent Cancer. 2007, 2: 11-10.1186/1750-9378-2-11.View ArticlePubMedPubMed CentralGoogle Scholar
  31. Meftah El khair M, El Mzibri M, Mhand RA, Benider A, Benchekroun N, Fahime EM, Benchekroun MN, Ennaji MM: Molecular detection and genotyping of human papillomavirus in cervical carcinoma biopsies in an area of high incidence of cancer from Moroccan women. J Med Virol. 2009, 81: 678-684. 10.1002/jmv.21279.View ArticlePubMedGoogle Scholar
  32. Tornesello ML, Losito S, Benincasa G, Fulciniti F, Botti G, Greggi S, Buonaguro L, Buonaguro FM: Human papillomavirus (HPV) genotypes and HPV16 variants and risk of adenocarcinoma and squamous cell carcinoma of the cervix. Gynecol Oncol. 2011, 121: 32-42. 10.1016/j.ygyno.2010.12.005.View ArticlePubMedGoogle Scholar
  33. Lee SH, Vigliotti VS, Vigliotti JS, Pappu S: Validation of human papillomavirus genotyping by signature DNA sequence analysis. BMC Clin Pathol. 2009, 9: 3-10.1186/1472-6890-9-3.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Lizano M, De la Cruz-Hernandez E, Carrillo-Garcia A, Garcia-Carranca A, Ponce de Leon-Rosales S, Duenas-Gonzalez A, Hernandez-Hernandez DM, Mohar A: Distribution of HPV16 and 18 intratypic variants in normal cytology, intraepithelial lesions, and cervical cancer in a Mexican population. Gynecol Oncol. 2006, 102: 230-235. 10.1016/j.ygyno.2005.12.002.View ArticlePubMedGoogle Scholar
  35. Seedorf K, Krammer G, Durst M, Suhai S, Rowekamp WG: Human papillomavirus type 16 DNA sequence. Virology. 1985, 145: 181-185. 10.1016/0042-6822(85)90214-4.View ArticlePubMedGoogle Scholar
  36. Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV, Snijders PJ, Meijer CJ: Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003, 348: 518-527. 10.1056/NEJMoa021641.View ArticlePubMedGoogle Scholar
  37. Clifford GM, Smith JS, Aguado T, Franceschi S: Comparison of HPV type distribution in high-grade cervical lesions and cervical cancer: a meta-analysis. Br J Cancer. 2003, 89: 101-105. 10.1038/sj.bjc.6601024.View ArticlePubMedPubMed CentralGoogle Scholar
  38. de Sanjose S, Quint WG, Alemany L, Geraets DT, Klaustermeier JE, Lloveras B, Tous S, Felix A, Bravo LE, Shin HR, et al: Human papillomavirus genotype attribution in invasive cervical cancer: a retrospective cross-sectional worldwide study. Lancet Oncol. 2010, 11: 1048-1056. 10.1016/S1470-2045(10)70230-8.View ArticlePubMedGoogle Scholar
  39. Zehbe I, Voglino G, Delius H, Wilander E, Tommasino M: Risk of cervical cancer and geographical variations of human papillomavirus 16 E6 polymorphisms. Lancet. 1998, 352: 1441-1442. 10.1016/S0140-6736(05)61263-9.View ArticlePubMedGoogle Scholar
  40. Pande S, Jain N, Prusty BK, Bhambhani S, Gupta S, Sharma R, Batra S, Das BC: Human papillomavirus type 16 variant analysis of E6, E7, and L1 genes and long control region in biopsy samples from cervical cancer patients in north India. J Clin Microbiol. 2008, 46: 1060-1066. 10.1128/JCM.02202-07.View ArticlePubMedPubMed CentralGoogle Scholar
  41. Chan SY, Ho L, Ong CK, Chow V, Drescher B, Durst M, ter Meulen J, Villa L, Luande J, Mgaya HN, Bernard HU: Molecular variants of human papillomavirus type16 from four continents suggest ancient pandemic spread of the virus and its coevolution with humankind. J Virol. 1992, 66: 2057-2066.PubMedPubMed CentralGoogle Scholar
  42. Cornet I, Gheit T, Franceschi S, Vignat J, Burk RD, Sylla BS, Tommasino M, Clifford GM, the IARC HPV Variant Study Group: Human papillomavirus type 16 genetic variants: phylogeny and classification based on E6 and LCR. J Virol. 2013, 86: 6855-6861.View ArticleGoogle Scholar
  43. Pista A, Oliviera A, Barateiro A, Costa H, Verdasca N, Paixão : Molecular variants of human papillomavirus type 16 and 18 and risk for cervical neoplasia in Portugal. J Med Virol. 2007, 79: 1889-1897. 10.1002/jmv.21002.View ArticlePubMedGoogle Scholar
  44. Tornesello ML, Duraturo ML, Salatiello I, Buonaguro L, Losito S, Botti G, Stellato G, Greggi S, Piccoli R, Pilotti S, et al: Analysis of human papillomavirus type-16 variants in Italian women with cervical intraepithelial neoplasia and cervical cancer. J Med Virol. 2004, 74: 117-126. 10.1002/jmv.20154.View ArticlePubMedGoogle Scholar
  45. Tu JJ, Kuhn L, Denny L, Beattie KJ, Lorincz A, Wright TC: Molecular variants of human papillomavirus type 16 and risk for cervical neoplasia in south Africa. Int J Gynecol Cancer. 2006, 16: 736-742. 10.1111/j.1525-1438.2006.00401.x.View ArticlePubMedGoogle Scholar
  46. Chopjitt P, Ekalaksananan T, Pientong C, Kongyingyoes B, Kleebkaow P, Charoensri N: Prevalence of human papillomavirus type 16 and its variants in abnormal squamous cervical cells in northeast Thailand. Int J Infect Dis. 2009, 13: 212-219. 10.1016/j.ijid.2008.06.017.View ArticlePubMedGoogle Scholar
  47. Hu X, Guo Z, Tianyun P, Pont en F, Wilander E, Andersson S, Pont en J: HPV typing and HPV16 E6-sequence variations in synchronous lesions of cervical squamous-cell carcinoma from Swedish patients. Int J Cancer. 1999, 83: 34-37. 10.1002/(SICI)1097-0215(19990924)83:1<34::AID-IJC7>3.0.CO;2-Q.View ArticlePubMedGoogle Scholar
  48. Sun M, Gao L, Liu Y, Zhao Y, Wang X, Pan Y, Ning T, Cai H, Yang H, Zhai W, Ke Y: Whole genome sequencing and evolutionary analysis of human papillomavirus type 16 in central China. PloS one. 2012, 7: e36577-10.1371/journal.pone.0036577.View ArticlePubMedPubMed CentralGoogle Scholar
  49. Chakrabarti O, Veeraraghavalu K, Tergaonkar V, Liu Y, Androphy EJ, Stanley MA, Krishna S: Human papillomavirus type 16 E6 amino acid 83 variants enhance E6-mediated MAPK signaling and differentially regulate tumorigenesis by notch signaling and oncogenic Ras. J Virol. 2004, 78: 5934-5945. 10.1128/JVI.78.11.5934-5945.2004.View ArticlePubMedPubMed CentralGoogle Scholar
  50. van Duin M, Snijders PJ, Vossen MT, Klaassen E, Voorhorst F, Verheijen RH, Helmerhorst TJ, Meijer CJ, Walboomers JM: Analysis of human papillomavirus type 16 E6 variants in relation to p53 codon 72 polymorphism genotypes in cervical carcinogenesis. J Gen Virol. 2000, 81: 317-325.View ArticlePubMedGoogle Scholar
  51. Zehbe I, Richard C, DeCarlo CA, Shai A, Lambert PF, Lichtig H, Tommasino M, Sherman L: Human papillomavirus 16 E6 variants differ in their dysregulation of human keratinocyte differentiation and apoptosis. Virology. 2009, 383: 69-77. 10.1016/j.virol.2008.09.036.View ArticlePubMedGoogle Scholar
  52. Stöppler H, Koval D, Schlegel R: The serine protease inhibitors TLCK and TPCK inhibit the in vitro immortalization of primary human keratinocytes by HPV-18 DNA. Oncogene. 1996, 13: 1545-1548.PubMedGoogle Scholar
  53. Wu Y, Chen Y, Li L, Yu G, He Y, Zhang Y: Analysis of mutations in the E6/E7 oncogenes and L1 gene of human papillomavirus 16 cervical cancer isolates from China. J Gen Virol. 2006, 87: 1181-1188. 10.1099/vir.0.81649-0.View ArticlePubMedGoogle Scholar
  54. Tonon SA, Basiletti J, Badano I, Alonio LV, Villa LL, Teyssie AR, Picconi MA: Human papillomavirus type 16 molecular variants in Guarani Indian women from Misiones, Argentina. Int J Infect Dis. 2007, 11: 76-81. 10.1016/j.ijid.2006.01.002.View ArticlePubMedGoogle Scholar
  55. Picconi MA, Alonio LV, Sichero L, Mbayed V, Villa LL, Gronda J, Campos R, Teyssié A: Human papillomavirus type-16 variants in Quechua aboriginals from Argentina. J Med Virol. 2003, 69: 546-555. 10.1002/jmv.10343.View ArticlePubMedGoogle Scholar
  56. DeFilippis RA, Goodwin EC, Wu L, DiMaio D: Endogenous human papillomavirus E6 and E7 proteins differentially regulate proliferation, senescence, and apoptosis in HeLa cervical carcinoma cells. J Virol. 2003, 77: 1551-1563. 10.1128/JVI.77.2.1551-1563.2003.View ArticlePubMedPubMed CentralGoogle Scholar
  57. Chen Z, Terai M, Fu L, Herrero R, DeSalle R, Burk RD: Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses. J Virol. 2005, 79: 7014-7023. 10.1128/JVI.79.11.7014-7023.2005.View ArticlePubMedPubMed CentralGoogle Scholar
  58. Zehbe I, Mytilineos J, Wikström I, Henriksen R, Edler L, Tommasino M: Association between human papillomavirus 16 E6 variants and human leukocyte antigen class I polymorphism in cervical cancer of Swedish women. Hum Immunol. 2003, 64: 538-542. 10.1016/S0198-8859(03)00033-8.View ArticlePubMedGoogle Scholar
  59. Ellis JR JKP, Baird J, Hounsell EF VRD, Rowe M, Hopkins D, Duggan-Keen MF, Bartholomew JS, Young LS: The association of an HPV16 oncogene variant with HLA-B7 has implications for vaccine design in cervical cancer. Nat Med. 1995, 1: 464-470. 10.1038/nm0595-464.View ArticlePubMedGoogle Scholar
  60. Bartholomew JS, Stacey SN, Coles B, Burt DJ, Arrand JR, Stern PL: Identification of a naturally processed HLA A0201-restricted viral peptide from cells expressing human papillomavirus type 16 E6 oncoprotein. Eur J Immunol. 1994, 24: 3175-3179. 10.1002/eji.1830241239.View ArticlePubMedGoogle Scholar
  61. Pre-publication history

    1. The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/13/378/prepub

Copyright

© Qmichou et al.; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.