|
|
||||||||
ORIGINAL RESEARCH |
From the Department of Psychiatry, and the Department of Obstetrics and Gynecology, Division of Gynecologic Oncology, the Lauren V. Ackerman Laboratory of Surgical Pathology, and the Department of Genetics, Washington University School of Medicine, St. Louis, Missouri; the HLA Laboratory, Barnes-Jewish Hospital, St. Louis, Missouri; and the Department of Pathology, University of Alabama, Birmingham, Alabama.
Address reprint requests to: Rosalind J. Neuman, PhD, Washington University Medical Center, Box 8134, 4940 Childrens Place, St. Louis, MO 63110, E-mail: roz{at}gretta.wustl.edu
| Abstract |
|---|
|
|
|---|
Methods: Ninety-six patients with pathologically confirmed invasive cervical cancer were ascertained. Human papillomavirus types were determined in 80 patients, of whom 81.25% were HPV-positive, and 18.75% were HPV-negative. Deoxyribonucleic acid was extracted from samples, taken from patients and their parents, and sequenced to determine DQB1 genotypes. Nuclear family data were used to test whether the DQB1 locus is associated with invasive cervical cancer while controlling for high-risk HPV-positive patients. The transmission/disequilibrium test evaluates whether the frequency of transmission of parental marker alleles to their affected offspring deviates from the expected Mendelian frequency of 50%.
Results: The HLA DQB1 locus showed evidence for allelic association with invasive cervical cancer in high-risk HPV-positive patients (P = .006). The transmission/disequilibrium test showed that the DQB1*0303 allele was transmitted to high-risk HPV patients more often than expected by chance,
21 = 8.0, P = .005 (P = .035 when correcting for multiple tests). Tests of association were negative when applied to all 96 patients, irrespective of HPV status. No significant differences were found in the distribution of the DQB1 alleles among HPV-positive patients compared with those who were HPV-negative, indicating that HLA alleles are not associated with susceptibility to HPV infection.
Conclusion: These results suggest that the DQB1*0303 allele increases the risk for invasive cervical cancer in women who are HPV-positive.
Cervical cancer is the second leading cause of cancer deaths in women worldwide and the leading cause in many developing countries.1 The age-adjusted incidence of cervical cancer is estimated to be between ten and 40 per 100,000.2 In the United States, the age-adjusted incidence and mortality rates of cervical cancer between 1991 and 1995 were 8.0 and 2.8, respectively, per 100,000,3 with the reduction due to the extensive use of Papanicolaou-stained smears as a screening tool.
Infection with human papillomavirus (HPV) is the risk factor most often associated with cervical neoplasms with the highest risk attributed to HPV types 16 and 18, followed by HPV 31 and 33.4 Prevalence estimates of genital HPV infection vary according to the method of HPV diagnosis and the population studied. Estimates ranging from 10% to 40% have been reported for the United States.1,2,5 Although HPV infection is an important risk factor for development of cervical neoplasms, it is clear that it is neither necessary nor sufficient. Although up to 80% of cervical tumors harbor HPV DNA,1,6 most infected women do not develop such tumors. Thus, additional environmental or genetic factors must be involved in the predisposition to cervical carcinogenesis.
Compelling evidence has implicated genetic factors as another important risk factor for susceptibility to invasive cervical cancer. Most, but not all, such studies have pointed to specific HLA class II alleles as candidates for involvement in onset of cervical intraepithelial neoplasia (CIN) or invasive cervical cancer. Some have reported a strong positive association with alleles at the DQB1 locus, dependent upon HPV type,712 whereas others report no association.13,14 Haplotypes formed from the DR and DQ loci have produced positive15,16 and negative findings.12 Confounded with the genetic evidence is the increase in high-risk HPV status observed in CIN and invasive cervical cancer patients. The exact relationship between these two risk factors is not well understood. However, it is clear that cervical cancer does not follow a simple Mendelian pattern of inheritance but, rather, is most likely the result of a complex interplay of genetic and environmental risk factors.
Human leukocyte antigen (HLA) genetic studies of cervical cancer have been primarily traditional case-control studies in which the frequencies of alleles among cases and controls are compared. The pitfalls of case-control studies have been well documented. Problems arise from the use of inappropriate control groups giving rise to spurious evidence for association due to artifacts such as population stratification. Furthermore, causal inferences based on genetic differences between cases and controls are often difficult to replicate. In response, family-based controls or haplotype relative-risk (RR) methods1719 are recommended for genetic association studies.20
The present study used the transmission/disequilibrium test19 for a family-based genetic association study of invasive cervical cancer and alleles at the DQB1 locus. This test focuses on transmission of a parental marker allele to an affected offspring, comparing this transmission with those alleles that are not transmitted. Consequently, a perfectly matched control is created for each case. In particular, we applied a likelihood-based form of the transmission/disequilibrium test to these data.21 The likelihood form of this statistical test is particularly appropriate because of the large number of alleles at the DQB1 locus, making traditional haplotype RR methods less reliable because statistical problems arise when multiple tests are performed.
| Materials and Methods |
|---|
|
|
|---|
Deoxyribonucleic acid was extracted from whole blood by first lysing the red cells in a solution containing 0.32 M sucrose, 10 mM Tris-Cl pH 7.5, 5 mM MgCl2, and 1% Triton X-100 (Sigma, St. Louis, MO). White blood cells were pelleted and washed twice in the same solution to remove hemoglobin. The cell pellet was then resuspended and incubated in 100 µL 10 mM Tris-Cl, 1 mM ethylenediaminetetracetic acid (EDTA), pH 8.0, 20 µg/mL proteinase K, and 0.5% NP40 at 50C for 2 hours. The proteinase K was inactivated by boiling for 20 minutes. The buccal cells were scraped from the brush, pelleted, and resuspended in 2040 µL of 10 mM Tris-Cl, pH 8.0. Four µL of buccal cell solution was added to 10 mM Tris-HCl pH 9.0, 50 mM KCl, 2.0 mM MgCl2, and 14.8 µg of proteinase K for a total volume of 50 µL. The tube was incubated at 60C for 20 minutes and then the proteinase K was inactivated by boiling for 20 minutes. Five µL of this solution was used for polymerase chain reaction (PCR) amplification.
Deoxyribonucleic acid from blood or buccal cells was amplified by PCR using the DQB1 exon-2-specific primers DQBAMP-A and DQBAMP-4 or DQBAMP-6.22 The DQBAMP-4 and DQBAMP-6 primers selectively amplified the *05 and *06 or *02, *03, and *04 alleles, respectively. Therefore, two amplifications were done for each individual. The 2030 µL PCR reaction contained 2050 ng of blood DNA or 4 µL of buccal cell solution, 10 mM Tris-HCl pH 9.0, 50 mM KCl, 2.0 mM MgCl2, 0.2 mM deoxyribonucleotide triphosphates, 3.2 pmole of each primer, and 0.25 units Taq (Promega, Madison, WI) or AmpliTaq (Perkin-Elmer, Applied Biosystems Division, Foster City, CA) DNA polymerase and .1% Triton X-100 (when Taq polymerase was used). The reaction mixture was held at 94C for 10 minutes, followed by 35 cycles at 93C for 1 minute, 60C for 1 minute, and 72C for 1 minute in an Omnigene Thermal Cycler (Hybaid, Middlesex, UK). Five µL of each PCR was run on a 2% agarose gel and the presence of a DNA fragment was detected by ethidium bromide staining. Amplified samples were then used as sequencing templates.
Purification of PCR products was done by adding 2 to 8 µL of the PCR reaction to 10 units of exonuclease I and 2 units of shrimp alkaline phosphatase (Amersham Life Science/United States Biochemical, Cleveland, OH). The reaction was held at 37C for 25 minutes and 80C for 15 minutes. The sequencing reaction was done using the ABI PRISM Dye Terminator Cycle Sequencing Ready Reaction Kit (Perkin-Elmer) with AmpliTaq DNA polymerase, FS, according to instructions. After purifying the reaction product by spinning it through a Sephadex G50 column, the column effluent was precipitated, dried, and denatured before it was loaded on an ABI model 373 or 377 sequencing gel system (Perkin-Elmer). Genotypes for DQB1 were determined by reading one or two 180 bp sequence from exon 2 on each individual. By this method, 21 alleles could be identified. There were no discriminating sequences to distinguish between the 0601106012, 060510609, 02010202, 0303103032, or 04010402 alleles. Because the polymorphism is high at this locus, typing served as a check for paternity. In cases of nonpaternity, if the transmitted or nontransmitted alleles from the mother could be resolved unambiguously, they were used in the analysis. Homozygous samples and samples that could not be sequenced were also typed by dot-blot assay. Briefly, genomic DNA was amplified by PCR using DQB1-specific primers and hybridized with 16 oligonucleotide probes.23 Detection of the DQB1 alleles was performed by colorimetric reverse-dot-blot hybridization ( Duffy BF, Begovich AB, Novotny JF, Mohanakumar T. Allele level HLA DRB typing with Amplicor typing kit using an extended panel of primers and probes [abstract]. Hum Immunol 1997;55:137). This method also distinguished the genotypes 0301/0302 from 03032/0304 and 0603/0604 from 0607/0608 because heterozygous peaks in the sequencing trace could not be discerned as distinct alleles.
Tumor specimens were obtained for pathologic evaluation and HPV typing. One to three 25 µm sections of formalin-fixed paraffin-embedded tumors were added to 200 µL of digestion buffer (50 mM Tris-Cl pH 8.5, 1 mM EDTA, 0.5% Tween 20, 200 µg proteinase K/mL) and incubated overnight at 65C. After incubation, 100 µL of Chelex 100 (BioRad, Richmond, CA), as a deionized water slurry, was added and the samples were boiled for 10 minutes and spun at 14,000 rpm in a microcentrifuge for 5 minutes. The aqueous phase (DNA) above the Chelex beads and below the paraffin surface was removed and transferred to sterile micro-centrifuge tubes. Human papillomavirus typing was done by PCR amplification of the E6-E7 gene using the degenerate consensus primers pU-1M/pU-2R and restriction-fragment-length polymorphism analysis (RFLP).24,25
We used family-based data to test for association or linkage disequilibrium between the DQB1 locus and invasive cervical cancer. The term family-based indicates that the data consist of nuclear families rather than unrelated individuals designated as cases and controls. In the family-based design, the transmission of each parental allele is tracked to affected offspring, and the alleles that are not transmitted are considered to belong to the control. Figure 1
indicates that the A alleles were transmitted to the affected offspring and the B alleles were not transmitted. The strength of the family-based design is that it provides a perfectly matched control to the affected offspring. Traditional case-control designs may give rise to spurious disease associations because of ethnic admixture within a population or other problems relating to inappropriate controls. A family-based test neatly circumvents these problems by providing a pseudo or internal control perfectly matched to the affected offspring.
|
21) of the form (b - c)2/(b + c) is computed, where b is number of times the allele of interest is transmitted, and c is the number times it is not transmitted. To avoid multiple statistical tests, one for each allele at a highly polymorphic marker such as at the DQB1 locus (resulting in loss of power), we used a multi-allelic likelihood-ratio based approach to test for linkage disequilibrium.21 This multi-allelic approach jointly tests whether any of the alleles is excessively transmitted by formulating the likelihood of the data in terms of one parameter which, under the null hypothesis of no allelic excess transmission, would be equal to zero. The advantage of this statistic is that it has only one degree of freedom, irrespective of the number of alleles of the marker, and therefore, provides a more powerful test of linkage disequilibrium than traditional tests. This test was first applied to all 96 nuclear families and then applied only to those families in which the proband had a high-risk HPV type.
| Results |
|---|
|
|
|---|
|
|
|
21 = 8.0, P = .005 (P = .035, corrected for multiple tests). To eliminate the possibility that the HLA locus conferred susceptibility to HPV infection, we examined the distribution of the DQB1 alleles in HPV-positive compared with negative cases. We found no significant difference in the distribution of alleles between groups (P = .12).
| Discussion |
|---|
|
|
|---|
These findings are consistent with those of other studies in which a positive association was detected between DQB1*03 and either invasive cervical cancer or the precursor lesion, CIN. Wank and Thomssen7 first reported an association between DQw3 antigens and cervical squamous cell carcinoma in a German population. A predominance of DQB1*0301 and DQB1*0303 alleles were then found in the same set of patients when the initial study was extended using molecular typing.26 These results have also been replicated in populations from Norway,10 Japan,9 England,27 and France,28 as well as in a black sample.11 In white British patients, Odunsi et al12,29 found a significant association between DQB1*03 and CIN. Vandenwelde et al30 reported a significantly elevated frequency of the DQB1*03 alleles in high-risk HPV-positive Belgian women diagnosed with CIN grades I and II, but not in those with CIN-III. However, David et al31 found a strong association between DQB1*03 and CIN-III in a British sample. Moderate increases of DQB1*0303 among HPV16-positive cases have been reported in a New Mexico Hispanic population32 and in a French population.28
The next most frequent allele found to be associated with cervical cancer is DQB1*0602, which has been reported in women from Tanzania8 and Hispanics from the United States.32 Sanjeevi et al33 found DQB1*0602 to be strongly associated with CIN in HPV16-positive Swedish patients compared with HPV16-positive controls. Several haplotypes containing the DQB1*0602 allele also have been found to be over-represented in patient populations. In a Norwegian study, Helland et al16 found an increase of the DQA1*0102-DQB1*0602 haplotype among HPV-positive CIN grade II and III cases compared with controls. Apple et al32,15 found a significant increase in haplotype DRB1*1501-DQB1*0602 in a New Mexico Hispanic population diagnosed with invasive cervical cancer and CIN, particularly in HPV16-positive patients. Sanjeevi et al33 found the DQA1*0102-DQB1*0602 haplotype significantly increased among HPV16-positive CIN patients compared with HPV16-positive controls. Similarly, a study from Norway showed the same DQA1*0102-DQB1*0602 haplotype associated with an increased risk of developing CIN when infected with HPV16 without influencing progression to cancer.16 It should be noted that linkage disequilibrium exists among many of the HLA class II haplotypes34 and therefore, the increased haplotype frequency, in fact, may be due to only one of the loci.
Several negative findings have been reported. Helland et al16 found no increase in DQB1*03 alleles in a British population of CIN patients. No association between cervical carcinoma with alleles encoding the DQ3 antigen was found in a population of Swedish patients.35 Glew et al13,14 found no positive correlation between invasive cervical cancer and HLA DR-DQ haplotypes in their British population samples, irrespective of HPV status. The conflicting reports on HLA associations with invasive cervical cancer and CIN may be due to differences in the ethnic origin of patients and controls among the many various studies, selection of poorly matched control groups within each study, differences in methods of HLA typing, interacting environmental factors, genetic heterogeneity, and differences in statistical methodologies.
In general, overwhelming evidence indicates that, in addition to high-risk HPV types, the HLA complex contains genetic factors that increase the propensity for invasive cervical cancer among some populations. However, many women infected with high risk HPV types never progress to invasive cervical cancer. This phenomenon may be partially explained by the dependence of the immune systems response to HPV infections on specific genotypes at the DQB1 locus or on as yet undefined additional genetic or environmental factors. There are also some women with cervical cancer, albeit a small minority, who do not show any detectable signs of HPV infections. These cases may represent an etiologically separate group of patients with unknown genetic or environmental or both agents that cause invasive cervical malignancies. Consequently, it has become increasingly clear that invasive cervical cancer arises from a complex interplay between genetic and environmental factors, which may differ across disparate populations.
| Footnotes |
|---|
Received March 15, 1999. Received in revised form June 23, 1999. Accepted July 1, 1999.
| References |
|---|
|
|
|---|
2. Bosch FX, Manos MM, Munoz N, Sherman M, Jansen AM, Peto J, et al. Prevalence of human papillomavirus in cervical cancer: A worldwide perspective. International Biological Study on Cervical Cancer (IBSCC) Study Group. J Natl Cancer Inst 1995;87:796802.
3. Ries LAG, Kosary CL, Hankey BF, Miller BA, Edwards BK, eds. SEER cancer statistics review 19731995. Bethesda: National Cancer Institute, 1998.
4. Schiffman MH, Brinton LA. The epidemiology of cervical carcinogenesis. Cancer 1995;76:1888901.[Medline]
5. Koutsky L. Epidemiology of genital human papillomavirus infection. Am J Med 1997;102:38.[Medline]
6. Beutner KR, Tyring S. Human papillomavirus and human disease. Am J Med 1997;102:915.[Medline]
7. Wank R, Thomssen C. High risk of squamous cell carcinoma of the cervix for women with HLA-DQw3. Nature 1991;352:7235.[Medline]
8. Wank R, Meulen J, Luande J, Enerhardt HC, Pawlita M. Cervical intraepithelial neoplasia, cervical carcinoma, and risk for patients with HLA-DQB1*0602, *301, *0303 alleles. Lancet 1993;341:1215.
9. Nawa A, Nishiyama Y, Kobayashi T, Wakahara Y, Okamoto T, Kikkawa F, et al. Association of human leukocyte antigen-B1*03 with cervical cancer in Japanese women aged 35 years and younger. Cancer 1995;75:51821.[Medline]
10. Helland A, Borresen AL, Kristensen G, Ronningen KS. DQA1 and DQB1 genes in patients with squamous cell carcinoma of the cervix: Relationship to human papillomavirus infection and prognosis. Cancer Epidemiol Biomarkers Prev 1994;3:47986.[Abstract]
11. Gregoire L, Lawrence D, Kukuruga D, Eisenbrey AB, Lancaster WD. Association between HLA-DQB1 alleles and risk for cervical cancer in African-American women. Int J Cancer 1994;57:5047.[Medline]
12. Odunsi K, Terry G, Ho L, Bell J, Cuzick J, Ganesan TS. Association between HLA DQB1*03 and cervical intra-epithelial neoplasia. Mol Med 1995;1:16171.[Medline]
13. Glew SS, Duggan-Keen M, Cabrera T, Stern PL. HLA class II antigen expression in human papillomavirus-associated cervical cancer. Cancer Res 1992;52:400916.
14. Glew SS, Duggan-Keen M, Ghosh AK, Ivinson A, Sinnott P, Davidson J, et al. Lack of association of HLA polymorphisms with human papillomavirus-related cervical cancer. Hum Immunol 1993;37:15764.[Medline]
15. Apple RJ, Becker TM, Wheeler CM, Erlich HA. Comparison of human leukocyte antigen DR-DQ disease associations found with cervical dysplasia and invasive cervical carcinoma. J Natl Cancer Inst 1995;87:42736.
16. Helland A, Olsen AO, Gjoen K, Akselsen HE, Sauer T, Magnus P, et al. An increased risk of cervical intra-epithelial neoplasia grade II-III among human papillomavirus-positive patients with the HLA-DQA1*0102-DQB1*0602 haplotype: A population-based case-control study of Norwegian women. Int J Cancer 1998;76:1924.[Medline]
17. Falk CT, Rubinstein P. Haplotype relative risks: An easy reliable way to construct a proper control sample for risk calculations. Ann Hum Genet 1987;51:22733.[Medline]
18. Terwilliger J, Ott J. A haplotype-based haplotype relative risk approach to detecting allelic associations. Hum Hered 1992;42:33746.[Medline]
19. Spielman RS, McGinnis RE, Ewens WJ. Transmission test for linkage disequilibrium: The insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet 1993;52: 50616.[Medline]
20. Lander ES, Schork NJ. Genetic dissection of complex traits. Science 1994;265:203748.
21. Terwilliger J. A powerful likelihood method for the analysis of linkage disequilibrium between trait loci and one or more polymorphic marker loci. Am J Hum Genet 1995;56:77787.[Medline]
22. Spurkland A, Knutsen I, Markussen G, Vartdal F, Egeland T, Thorsby E. HLA matching of unrelated bone marrow transplant pairs: Direct sequencing of in vitro amplified HLA-DRB1 and -DQB1 genes using magnetic beads as solid support. Tissue Antigens 1993;41:15564.[Medline]
23. Fernandez-Vina MA, Falco M, Xiaojiang G, Cerra M, Sun Y, Raimondi E, et al. DqA1*03 subtypes have different associations with DRB1 and DQB1 alleles. Hum Immunol 1994;39:2908.[Medline]
24. Fujinaga Y, Shimada M, Okazawa K, Fukushima M, Kato I, Fujinaga K. Simultaneous detection and typing of genital human papillomavirus DNA using the polymerase chain reaction. J Gen Virol 1991;72:103944.
25. Rader JS, Kamarasova T, Huettner PC, Li L, Li Y, Gerhard DS. Allelotyping of all chromosomal arms in invasive cervical cancer. Oncogene 1996;13:273741.[Medline]
26. Wank R, Schendel DJ, Thomssen C. HLA antigens and cervical carcinoma. Nature 1992;356:223.
27. Mehal WZ, Lo YMD, Herrington CS, Evans MF, Papadopoulos MC, Odunsi K, et al. Role of human papillomavirus in determining the HLA associated risk of cervical carcinogenesis. J Clin Pathol 1994;47:107781.
28. Sastre-Garau X, Loste MN, Vincent-Salomon A, Favre M, Mouret E, de la Rochefordiere A, et al. Decreased frequency of HLA-DRB1 13 alleles in French women with HPV-positive carcinoma of the cervix. Int J Cancer 1996;69:15964.[Medline]
29. Odunsi K, Terry G, Ho L, Bell J, Cuzick J, Ganesan TS. Susceptibility to human papillomavirus-associated cervical intra-epithelial neoplasia is determined by specific HLA DR-DQ alleles. Int J Cancer 1996;67:595602.[Medline]
30. Vandenvelde C, De Foor M, Van Beers D. HLA-DQB1*03 and cervical intraepithelial neoplasia grades IIII. Lancet 1993;341: 4423.
31. David ALM, Taylor GM, Gokhale D, Aplin JD, Seif MW, Tindall VR. HLA-DQB1*03 and cervical intraepithelial neoplasia type III. Lancet 1992;340:52.[Medline]
32. Apple RJ, Erlich HA, Klitz W, Manos M, Becker TM, Wheeler CM. HLA DR-DQ associations with cervical carcinoma show papillomavirus-type specificity. Nat Genet 1994;6:15762.[Medline]
33. Sanjeevi CB, Hjelmstrom P, Hallmans G, Wiklund F, Lenner P, Angstrom T, et al. Different HLA-DR-DQ haplotypes are associated with cervical intraepithelial neoplasia among human papillomavirus type-16 seropositive and seronegative Swedish women. Int J Cancer 1996;68:40914.[Medline]
34. Begovich AB, McClure GR, Suraj VC, Helmuth RC, Fildes N, Bugawan TL, et al. Polymorphism, recombination, and linkage disequilibrium within the HLA class II region. J Immunol 1992; 148:24958.[Abstract]
35. Allen M, Kalantari M, Ylitalo N, Pettersson B, Hagmar B, Scheibenpflug L, et al. HLA DQ-DR haplotype and susceptibility to cervical carcinoma: Indications of increased risk for development of cervical carcinoma in individuals infected with HPV 18. Tissue Antigens 1996;48:327.[Medline]
This article has been cited by other articles:
![]() |
Z. Zhang, I. Borecki, L. Nguyen, D. Ma, K. Smith, P. C. Huettner, D. G. Mutch, T. J. Herzog, R. K. Gibb, M. A. Powell, et al. CD83 Gene Polymorphisms Increase Susceptibility to Human Invasive Cervical Cancer Cancer Res., December 1, 2007; 67(23): 11202 - 11208. [Abstract] [Full Text] [PDF] |
||||
![]() |
M Zoodsma, I M Nolte, M Schipper, E Oosterom, G van der Steege, E G E de Vries, G J te Meerman, and A G J van der Zee Analysis of the entire HLA region in susceptibility for cervical cancer: a comprehensive study J. Med. Genet., August 1, 2005; 42(8): e49 - e49. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Lin, L. A. Koutsky, C. W. Critchlow, R. J. Apple, S. E. Hawes, J. P. Hughes, P. Toure, A. Dembele, and N. B. Kiviat HLA Class II DR-DQ and Increased Risk of Cervical Cancer among Senegalese Women Cancer Epidemiol. Biomarkers Prev., October 1, 2001; 10(10): 1037 - 1045. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |