|Year : 2014 | Volume
| Issue : 1 | Page : 14-22
Polymerase chain reaction and deoxyribonucleic acid-sequencing based study on distribution of human papillomavirus 16/18 among histopathological types of cervical intra-epithelial neoplasia and primary invasive cervical carcinoma: A scenario in North Bengal, India
Prithwijit Ghosh1, Damayanti Das Ghosh2, Amita Majumdar (Giri)3, Sharmila Sengupta4, Chandana Das5, Indranil Mukhopadhyay6
1 Department of Pathology, North Bengal Medical College and Hospital, Siliguri; Department of Pathology, DESUN (NEON) Lab, Kolkata, West Bengal, India
2 Human Genetics Unit, Biological Sciences Division, Indian Statistical Institute; Cancer Biology and Inflammatory Disorder Division, Indian Institute of Chemical Biology, Kolkata, West Bengal, India
3 Department of Pathology, North Bengal Medical College and Hospital, Siliguri, West Bengal, India
4 Human Genetics Unit, Biological Sciences Division, Indian Statistical Institute, Kolkata; National Institute of Biomedical Genomics, Kalyani, West Bengal, India
5 Department of Gynecology and Obstetrics, North Bengal Medical College and Hospital, Siliguri; Nil Ratan Sircar Medical College, Kolkata, West Bengal, India
6 Human Genetics Unit, Biological Sciences Division, Indian Statistical Institute, Kolkata, West Bengal, India
|Date of Web Publication||24-Feb-2014|
25, 2nd Road, East End Park, Kolkata - 700 099, West Bengal
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Introduction: Human papillomavirus (HPV) types 16/18 are reportedly most common in cervical cancer (CaCx) with geographical variation of genotypes. HPV16 predominates both in squamous cell carcinoma (SCC) and adenocarcinoma in India, contrary to reported global predominance of HPV18 in the latter. Our study was aimed to determine the occurrence of HPV16/18 among histopathological types of cervical intra-epithelial neoplasia (CIN) and invasive CaCx from North Bengal, India and to identify any major deviation from the known Indian scenario of distribution of HPV16/18 genotypes in cases of SCC and adenocarcinoma.
Materials and Methods: This was a retrospective, cross-sectional, case-only type of study, in which 40 cases were histopathologically diagnosed as CIN/CaCx, on which polymerase chain reaction (PCR), deoxyribonucleic acid (DNA)-sequencing and bioinformatics by basic search local alignment tool were performed for HPV-genotyping.
Statistical Analysis: The distribution of HPV genotypes among cases of SCC and adenocarcinoma was compared by Fisher's exact-test.
Results: HPV was detected in 97.5% (39/40) cases. HPV16-infected cases (32/39; 82.05%) predominated over HPV18-infected ones (7/39; 17.95%). However, HPV18-only infection was significantly (P = 0.0045, one-sided Fisher's exact test) more among adenocarcinoma (3/4; 75%) than SCC (2/26; 7.69%) contrary to HPV16-only infection (SCC = 24/26, 92.31%; adenocarcinoma = 1/4; 25%) whereas both CIN3 cases were HPV16-positive.
Conclusion: Predominance of HPV18 over HPV16 in cases of adenocarcinoma in this region was contrasting to that of earlier Indian studies suggesting research on HPV18 related cervical carcinogenesis. PCR and DNA-sequencing could prove to be highly effective tools in HPV detection and genotyping. The study reported HPV16/18 infection in almost 98% of the cases, the knowledge about which might prove useful in future population based studies on HPV genotyping and designing of appropriate HPV-vaccines for this region.
Keywords: Cervical cancer, histopathology, human papillomavirus 16/18, polymerase chain reaction, North Bengal
|How to cite this article:|
Ghosh P, Ghosh DD, Majumdar (Giri) A, Sengupta S, Das C, Mukhopadhyay I. Polymerase chain reaction and deoxyribonucleic acid-sequencing based study on distribution of human papillomavirus 16/18 among histopathological types of cervical intra-epithelial neoplasia and primary invasive cervical carcinoma: A scenario in North Bengal, India. J Mid-life Health 2014;5:14-22
|How to cite this URL:|
Ghosh P, Ghosh DD, Majumdar (Giri) A, Sengupta S, Das C, Mukhopadhyay I. Polymerase chain reaction and deoxyribonucleic acid-sequencing based study on distribution of human papillomavirus 16/18 among histopathological types of cervical intra-epithelial neoplasia and primary invasive cervical carcinoma: A scenario in North Bengal, India. J Mid-life Health [serial online] 2014 [cited 2021 May 12];5:14-22. Available from: https://www.jmidlifehealth.org/text.asp?2014/5/1/14/127786
| Introduction|| |
Uterine cervical cancer (CaCx) is the third most common cancer among women world-wide with more than 85% of the global burden reported from the developing countries and India shares about one-fourth of the burden (http://globocan.iarc.fr/). The major etiological factor is persistent infection with oncogenic human papillomaviruses (HPVs), whereas early age of marriage, multiparty, smoking, etc., are probable cofactors.  HPVs are small, double-stranded deoxyribonucleic acid (DNA) viruses that belong to family papillomaviridae and they have been categorized into high-risk types HPV (HR-HPV), such as 16, 18, 31, 33, etc. and low-risk types HPV (LR-HPV), such as 6, 11, etc., with geographical variation in genotypic distribution.  The viral genome encodes for late (L1, L2) and early (oncogenic E6, E7, etc.) transcribing genes. 
Early precancerous cervical intra-epithelial neoplasias (CIN1/2) reportedly have higher chances of spontaneous regression compared with advanced CIN3, which progresses more frequently to invasive cancer.  Two major histopathological (H/P) types of CaCx are squamous cell carcinoma (SCC; 80-85%) and adenocarcinoma (10%). HPV16 is reportedly the most prevalent type both in cervical SCC and adenocarcinoma in India, while preferential occurrence of HPV18 has been reported in adenocarcinoma globally.  Classification systems such as the "World Health Organization classification of Epithelial Tumors of the Uterine Cervix", Modified Broder's, etc., have been used for H/P classification of cervical neoplasm. , Light microscopy can detect "koilocyte", the morphological marker of HPV infection.  More sensitive and specific molecular methods, like polymerase chain reaction (PCR), DNA-sequencing coupled with the application of bioinformatics like basic local alignment search tool (BLAST), have been used for HPV detection and genotyping. ,,
The currently available prophylactic HPV-vaccines reportedly provide about 80% protection, the remaining non-coverage being accounted by rarer HR-HPV types.  HPV genotyping is thus vital to procure knowledge about regional genotype distribution to help in development of a cost-effective, efficacious HPV-vaccine and can be used as an adjunct to Pap-smear screening, which has its own limitations. 
North Bengal covers a large area comprising of six districts namely Uttar Dinajpur, Dakshin Dinajpur, Coochbehar, Siliguri, Jalpaiguri and Darjeeling. To the best of our knowledge, previous published data on HPV genotyping and its distribution among H/P types of CIN/CaCx from North Bengal is unavailable. Thus, our study was an attempt to determine the same and to identify any major deviation in distribution of HPV genotypes 16/18 in cases of SCC and adenocarcinoma in this region.
| Materials and Methods|| |
This was a North-Bengal based retrospective, cross-sectional, case-only type of study on distribution of HPV genotypes among H/P types of CIN/invasive CaCx, conducted over a period of 1 year (February, 2010 to January, 2011). Cervical biopsy tissue collection and H/P examination were performed respectively at Departments of Gynecology and Obstetrics and Pathology of a tertiary hospital, in North Bengal, while molecular methods for HPV genotyping, such as PCR and DNA-sequencing, were performed at the Human Genetics Unit (HGU) of Indian Statistical Institute, Kolkata with informed patient-consent and approval from the concerned ethical committees. Initially, elective cervical punch biopsy tissue collection was made from those patients at their first visit, presenting with clinically suspicious unhealthy cervix and H/P examination was carried out on those cases. Cases with an already previous H/P diagnosis of CIN/CaCx on or off treatment were excluded from the study. Out of initial 49 tissues, 40 were histopathologically confirmed as CIN/CaCx, while nine had chronic cervicitis and were excluded from the study. Thus, the study population finally comprised of 40 cases on which molecular methods like PCR, DNA sequencing coupled with BLAST were performed to determine the HPV genotype.
Types of specimen
The cervical biopsy tissues were divided into two parts in the operation theatre. One part was submitted in 10% buffered formalin and received at Department of Pathology for H/P examination and, the other part was collected in "RNAlater" (Ambion, Life Technologies Corporation, USA) and carried in ice-pack to HGU for HPV genotyping. HPV infection was screened by PCR with L1-specific primers MY09/11, which if negative, was reconfirmed by nested PCR with GP5/6 primers. Cases positive for MY09/11 or GP5/6 were typed for HPV16 and HPV18 by type-specific PCR (TS-PCR) with respective E6-specific primers. In the absence of both HPV16 and 18, the HPV-infected samples were subjected to DNA sequencing followed by bioinformatics-based genotyping using BLAST.  An overview of the study design has been depicted in [Figure 1].
|Figure 1: Flow chart depicting study design. TS-PCR: Type-specific polymerase chain reaction, HPV: Human papillomavirus, H and E: Hematoxylin and eosin, PAS: Periodic Acid-Schiff|
Click here to view
Histopathology of biopsy tissue
The formalin-fixed biopsy tissues were routinely processed, paraffin-embedded and stained with hematoxylin and eosin (H and E) and the slides were examined under light microscope (Olympus CH20i, Olympus Opto Systems, India).  The cases diagnosed as cervical adenocarcinoma with H and E staining, were further stained with Periodic Acid-Schiff stain (PAS) to demonstrate mucin. 
Genomic DNA of both host and the virus was together isolated and purified from the samples by using QIAamp DNA mini kit (QIAGEN, Hilden, Germany) following manufacturer's protocol. The concentration and purity of DNA were estimated by spectrophotometry.
HPV genotyping using PCR, DNA-sequencing and BLAST
The reaction mixture (20 μl) for L1-specific conventional PCR contained ×10 PCR buffer (Roche Diagnostics, Switzerland), 1.5 mM MgCl 2 , 200 μM dNTP, 60 ng primer (MY09/11), 1 U thermostable Fast Start Taq polymerase (Roche Diagnostics, Switzerland) and 100 ng DNA.  The PCR program consisted of initial denaturation (95°C; 5 min), 30 cycles of denaturation (95°C; 1 min), annealing (57°C; 2 min) and elongation (72°C; 1 min 30 s) and, final elongation (72°C; 5 min). The MY09/11-negative PCR-products were confirmed by nested PCR in a reaction mixture (20 μl) containing ×10 PCR buffer, 1 mM MgCl 2 , 100 μM dNTP, 50 ng primer (GP5/6), 0.5 U FS Taq and 0.3 μl of MY09/11-product.  The PCR program consisted of initial denaturation (95°C; 5 min), 40 cycles of denaturation (95°C; 30 s), annealing (42°C; 1 min) and elongation (72°C; 30 s) and, final elongation (72°C; 5 min).
The HPV16-specific PCR mixture (20 μl) contained ×10 PCR buffer, 2 mM MgCl 2 , 200 μM dNTP, 4 ng primers, 0.5 U FS Taq and 100 ng DNA. The touch-down PCR-program consisted of initial 14 cycles of denaturation (95°C; 30 s), annealing (61°C; 30 s; reduction by 0.5°C/cycle) and elongation (72°C; 30 s), followed by 23 cycles of denaturation (95°C; 30 s), annealing (55°C; 30 s) and elongation (72°C; 30 s), along with initial denaturation and final elongation as above. The HPV18-specific PCR mixture (20 μl) contained ×10 PCR buffer, 2 mM MgCl 2 , 200 μM dNTP, 20 ng primers, 0.5 U FS Taq and 100 ng DNA. The PCR program consisted of initial denaturation (95°C; 5 min), 33 cycles of denaturation (95°C; 30 s), annealing (53°C; 30 s) and elongation (72°C; 30 s) and, final elongation (72°C; 3 min). HPV16-positive CaSki cell line DNA or HPV18-positive HeLa cell line DNA as a positive control, as applicable and HPV-negative genomic DNA and doubly distilled water as negative controls. The PCR-products were checked by gel-electrophoresis. The primer-sequences and amplicon sizes are provided in [Table 1].
HPV-DNA sequencing for HPV16/18-negative cases
Following Sanger's method, the sequencing-PCR mixture (10 μl) consisted of ×5 sequencing buffer and 1 μl BidDye® Terminator (Applied Biosystems, Life Technologies Corporation, USA), 25 ng primer and 2 μl diluted PCR-product (MY09/11 or GP5/6).  The PCR program consisted of initial denaturation (96°C; 6 min), 35 cycles of denaturation (96°C; 6 min), annealing (55°C; 30 s) and elongation (72°C; 30 s) and, final elongation (72°C; 7 min). The sequencing-PCR products were precipitated by applying 3 M sodium acetate (pH 5.2) and absolute alcohol and washed with 70% ethanol. The final precipitate was resuspended in Hi-Di TM formamide (Applied Biosystems) and subjected to capillary electrophoresis (3100 PRISM, Applied Biosystems). The electropherograms were studied using DNASTAR (www.dnastar.com).
The DNA-sequences were analyzed, using the bioinformatics application, called BLAST, to detect the matching HPV genotypes from the database (http://blast.ncbi.nlm.nih.gov/). , Chances of matching for each HPV genotype could be determined by a scoring system based on 'Expect-value' (E-value). Lower E-value signified higher probability of matching. The HPV genotypes matching with the query-sequence were displayed in increasing order of E-value by using BLAST.
Occurrence of HPV genotypes among different H/P types of intra-epithelial and invasive cervical neoplasia was compared by using Fisher's exact test. Other risk factors like smoking habit and patient's age were compared against different H/P types and HPV genotypes by t-test. The statistical analyzes were carried out by statistical software package SPSS v16.0 for windows (IBM, India). The level of significance (α) for every test was fixed at 0.05.
| Results|| |
Data on demographic variables, clinical history and clinical findings
In this study, the mean age of cases, mean age of marriage, mean age of the first child birth (FCB) and the mean parity were respectively 49.22 years, 15.95 years, 18.35 years and 4.1 [Table 2]. Four cases were grand multipara (parity >6). Majority of the cases was Hindu (31/40; 77.5%) by religion. Eighty-five percent (34/40) of the cases was illiterate. Tobacco-users accounted for 27.5% (11/40) cases. Majority of them developed infiltrating SCC (10/11; 90.9%; relative risk [RR] = 1.146; 95% confidence interval [CI] = 0.88-1.49), whereas, only a single patient had infiltrating adenocarcinoma (1/11; 9.1%; RR = 0.66; 95% CI = 0.08-5.27). Age of detection for HPV infected cases belonged mostly to the range of 40-49 years (21/39; 53.85%), followed by 50-59 years and ≥60 years (7/39; 17.95% in each) and lastly, <40 years (4/39; 10.25%). The lowest age detected (35 years) was a case of HPV16+18 co-infection. There was no significant difference in mean age of detection between infiltrating SCC (mean age = 48.97 ± 8.77 years) and either infiltrating adenocarcinoma (mean age = 48.8 ± 9.49 years; P = 0.99; t-test) or CIN3 (mean = 54.5 ± 0.71 years; P = 0.39; t-test). The most common clinical symptom and sign were respectively serosanguinous discharge (29/40; 72.5%) and clinically evident cervical growth with contact-bleeding (27/40; 67.5%).
H/P types and subtypes
Out of the 40 cases, 33 (82.5%) were infiltrating SCC, 5 (12.5%) were infiltrating adenocarcinoma (mucinous) and 2 (5%) were CIN3. Twenty-two SCC cases (66.67%) were infiltrating non-keratinizing SCC (INKSCC) and 11 (33.33%) were infiltrating keratinizing SCC (IKSCC). Ten IKSCC cases (10/11; 90.91%) were moderately differentiated (IMDKSCC) and one (1/11; 9.09%) was well-differentiated. Four adenocarcinoma cases (4/5; 80%) were infiltrating moderately differentiated (IMD adenocarcinoma) and one (1/5; 20%) was infiltrating poorly differentiated (IPD adenocarcinoma). [Figure 2]a represents IKSCC showing characteristic keratin pearl formations. [Figure 2]b and c represents INKSCC and CIN3, respectively. [Figure 3]a and b represents H and E and corresponding PAS stained sections of adenocarcinoma.
|Figure 2: (a) Infiltrating keratinizing squamous cell carcinoma of cervix with keratin pearl formation (×100). (b) Infiltrating non-keratinizing squamous cell carcinoma of cervix (×100). (c) Cervical intraepithelial neoplasm, grade 3 (×400). All sections stained with H and E|
Click here to view
|Figure 3: Infiltrating mucinous cervical adenocarcinoma (a) hematoxylin and eosin stain (×100). (b) Periodic acid-Schiff stain showing apical cytoplasmic positivity (×400)|
Click here to view
HPV was detected in 39 (97.5%) out of 40 cases by PCR. There was a single HPV (−)ve case (1/40; 2.5%) undetected by PCR and sequencing. HPV genotypes were categorized into six groups, namely, (i) HPV16 only (27/39, 69.23%), (ii) HPV18 only (5/39, 12.82%), (iii) HPV16 + HPV18 (1/39, 2.56%), (iv) HPV16 + others (4/39, 10.26%), (v) HPV18 + others (1/39, 2.56%) and (vi) unidentified HPV type (1/39, 2.56%) [Figure 4]. Altogether, 82.05% HPV (+)ve cases carried HPV16 and 17.95% HPV (+)ve cases carried HPV18. Cases belonging to groups (i), (ii) and (iii) were positive for HPV16/18/both by TS-PCR. Cases belonging to groups (iv) and (v) were negative for HPV16/18 by TS-PCR, but were positive for the same along with other LR-HPV types (co-infection) by more sensitive DNA-sequencing coupled with BLAST, possibly due to low copies of HPV16/18.  Group (vi) failed to match with known HPV genotypes by BLAST.  Representative pictures of gel-electrophoreses for conventional PCR (L1-My09/11), nested PCR (GP5/6), HPV16/18 TS-PCR products are shown in [Figure 5], [Figure 6], [Figure 7] and [Figure 8], respectively, whereas [Figure 9] and [Figure 10] respectively depict DNA sequence-based electropherograms and BLAST-output.
|Figure 4: Percentage distribution of six different categories of human papillomavirus infection|
Click here to view
|Figure 5: L1 (My09/11) polymerase chain reaction products. Lane 1: 50 bp deoxyribonucleic acid ladder; lanes 2-7: Samples showing specific band for L1 at 450 bp; lane 8: L1 (−)ve sample; lane 9: Positive control; lane 10: Negative control|
Click here to view
|Figure 6: GP5/6 polymerase chain reaction products. Lane 6: 50 bp deoxyribonucleic acid ladder; lanes 3, 4, 7, 8: Samples showing specific band for GP5/6 at 150 bp; lane 2, 5: GP5/6 (−)ve sample; lane 1: Positive control; lane 9: Negative control|
Click here to view
|Figure 7: Human papillomavirus (HPV)16 type-specific polymerase chain reaction products. Lane 4: 50 bp ladder; lanes 3, 5-8: HPV16 (+) ve samples showing specific band size at 116 bp; lane 9: HPV16 (−)ve sample; lane 2: Positive control; lane 1: Negative control|
Click here to view
|Figure 8: Human papillomavirus (HPV)18 type-specific polymerase chain reaction products. Lane 7: 50 bp ladder; lanes 1-2 : HPV18 (+) ve samples showing specific band size at 100 bp; lanes 3-5: HPV18 (−)ve sample; lane 6: Positive control; lane 8: Negative control|
Click here to view
|Figure 9: Representative electropherograms showing aligned sequences of L1 obtained by polymerase chain reaction with MY09/11 primers|
Click here to view
|Figure 10: Representative basic local alignment search tool analysis. Human papillomavirus (HPV)-genotypes matching with the query-sequence are displayed in decreasing order of identity determined by a specific scoring system (E-value). The "color key" denotes "query coverage" (percentage of the query sequence overlapping with the matched sequence). In this figure, among all the above hits, probability of HPV16 and HPV6b are the highest|
Click here to view
Distribution of different H/P types and subtypes with respect to categories of HPV genotypes
The SCC cases harbored HPV groups (i) (24/33; 72.73%), (iv) (4/33; 12.12%), (ii) (2/33; 6.06%), (iii) (1/33; 3.03%) and (v) (1/33; 3.03%) in decreasing order of frequency. The HPV (−)ve case also belonged to SCC category. The adenocarcinomas harbored HPV groups (ii) (3/5; 60%), (i) (1/5; 20%) and (vi) (1/5; 20%) in decreasing order of frequency. CIN3 cases harbored group (i) only. The distributions are depicted in [Figure 11].
|Figure 11: Distribution of different histopathological types of cervical neoplasm with respect to human papillomavirus genotypes. Inf: Infiltrating, adenoca: Adenocarcinoma, SCC: Squamous cell carcinoma, CIN: Cervical intraepithelial neoplasia|
Click here to view
Thus, HPV16 only group (24/26; 92.31%) was predominant over HPV18 only (2/26; 7.69%) in SCC cases. Contrarily, HPV18 only group (3/4; 75%) was predominant over HPV16 only (1/4; 25%) in adenocarcinoma. HPV18-only infection was significantly more among adenocarcinoma than SCC (P = 0.0045, one-sided Fisher's exact test). All the CIN3 cases were infected with HPV16 only. The distributions have been depicted in [Figure 12].
|Figure 12: Distribution of human papillomavirus 16/18-only categories among histopathological types of cervical neoplasia. Inf: Infiltrating, adenoca: Adenocarcinoma, SCC: Squamous cell carcinoma, CIN: Cervical intraepithelial neoplasia|
Click here to view
Out of the 22 INKSCCs and 11 IKSCCs, respectively 17 and 9 harbored HPV16/18-only infections. HPV16 only group was predominant in INKSCC cases (15/17; 88.23%) over HPV18 only (2/17; 11.77%). All the nine IKSCC cases harbored HPV16-only infection. There was no statistically significant difference in HPV16-only infection between INKSCC and IKSCC in comparison to the HPV18-only infection (P = 0.42, Fisher's exact test). HPV18 only group (2/3; 66.67%) was predominant over HPV16-only (1/3; 33.33%) in IMD adenocarcinoma, while the single case of IPD adenocarcinoma harbored HPV18-only. The distributions have been depicted in [Figure 13].
|Figure 13: Distribution of human papillomavirus 16/18-only categories among histopathological subtypes of cervical neoplasia. INKSCC: Infiltrating non-keratinizing squamous cell carcinoma, IKSCC: Infiltrating keratinizing squamous cell carcinoma, IMD: Infiltrating moderately differentiated, IPD: Infiltrating poorly differentiated, adenoca: Adenocarcinoma|
Click here to view
| Discussion|| |
Our study detected HPV infection in 97.5% of cases of cervical neoplasia, which reiterates the pivotal role of HPV in cervical carcinogenesis. The study supported previous reports on role of co-factors such as early age of marriage (mean = 15.95 years), FCB (mean 18.35 years) and high parity (mean = 4.1) in CaCx.  The RR due to smoking was found to be higher for infiltrating SCC than infiltrating adenocarcinoma supporting earlier reports.  There was no significant difference (P = 0.99) between mean ages of infiltrating SCC (48.97 years) and adenocarcinoma (48.8 years) similar to many earlier reports. , However, contrasting reports exist on the increasing incidence of infiltrating adenocarcinoma among much younger women.  Supporting earlier reports, the present study has found, unlike the western scenario, a much later age of detection within 40-49 years, pointing toward the existing limitations of Pap-smear screening in developing countries.  Another peak of detection found at ≥60 years could possibly be due to a post-menopausal loss of immunity.  The lowest age of detection was found to be 35 years in a case of HPV16+18 co-infection. Previous studies have reported a lower age of detection in cases of co-infections due to a possible synergistic effect of multiple HR-HPV genotypes in CaCx. 
The present study found infiltrating SCC cases (82.5%) to predominate over adenocarcinomas (12.5%) supporting earlier reports from India and abroad. , However, it differed from certain previous reports on increasing incidence of adenocarcinoma.  CIN3 was detected in only 5% of the cases, indicating the need to strengthen Pap-screening programs. Cases of CIN1/2 were not found, possibly due to their higher chances of spontaneous regression.  Nevertheless, detection of HR-HPV types even in lower grade lesions would necessitate proper follow-up to monitor progression to invasive cancer.  No case of adenocarcinoma in situ (AIS) was found, possibly due to inadequate endocervical sampling.  However, progression to invasive cancer is reportedly much faster in AIS than squamous intra-epithelial lesion.  Cases of INKSCC (66.67%) were found to predominate over IKSCC (33.33%) and the moderately differentiated subtype predominated both amongst cases of IKSCCs (90.01%) and adenocarcinomas (80%), supporting earlier findings. 
The present study employing PCR, DNA-sequencing and BLAST with L1-specific MY09/11 and GP5/6 and, HPV16/18-specific primers, detected HPV in all but one case, reflecting high sensitivity of the methodology.  DNA-sequencing is reported as the gold-standard method for HPV genotyping.  The HPV-negative case was histopathologically IMDKSCC. HPV-negative CaCx can arise from transforming host genetic alterations.  However, HPV false-negatives can result from low viral load, disruption in HPV-genome during viral integration, less sensitive PCR conditions, poor quality of specimen, sampling error, absence of episomal forms in glandular lesions.  There was a single case of unidentified HPV-type histopathologically detected as infiltrating adenocarcinoma. Unidentified HPV types based on results of DNA-sequencing coupled with bioinformatics (BLAST) can be due to poor quality of sequence-data or represent novel HPV type. 
The present study found 97.44% (38/39) of the HPV-positive cases to be infected with HPV16/HPV18/both with HPV16 (32/39; 82.05%) predominating over HPV18 (7/39; 17.95%) infection, supporting previous reports. , Other studies have found HR-HPV types, like 31, 33 besides HPV16 predominance.  However HPV52 predominance has been reported from several South East Asian studies, reflecting geographical variation of HPV genotypes.  There was a single case of HPV16+18 co-infection, while other co-infections harbored LR-HPV types 6b, 32, 26, etc., along with HR-HPV16/18 similar to previous reports.  However, many studies have failed to find a synergistic effect of LR - along with HR-HPV in increasing severity of lesions. 
Considering the HPV16/18-only cases, 92.31% HPV16-only and 7.69% HPV18-only infections were found among SCC cases, while 25% HPV16-only and 75% HPV18-only among adenocarcinoma cases. HPV18-only infection was significantly higher among adenocarcinomas than SCCs (P = 0.0045, one-sided Fisher's exact test) supporting some earlier studies.  In contrast, HPV16 has been earlier reported to be the most predominant type, irrespective of H/P types of CaCx in India.  HPV16-only infection (88.23%) predominated over HPV18-only (11.7%) in INKSCC cases and was the sole detected type in IKSCC cases supporting earlier reports. 
| Conclusion|| |
Our study has provided possibly the first report from North Bengal, India on HPV-genotyping in cases of CIN/CaCx. Predominance of HPV18 over HPV16 in cases of adenocarcinoma in this region was contrasting to that of earlier Indian studies. This finding was new and contrary to that of the Indian scenario and should encourage future studies on HPV-18 related cervical carcinogenesis to validate this finding. PCR and DNA-sequencing coupled with BLAST could prove to be highly effective tools in HPV detection and genotyping, which might encourage adoption of the same in future such studies. The study reported HPV16/18 infection in almost 98% of the cases of cervical neoplasia, knowledge about which might prove useful in future population based studies on HPV genotyping and designing of appropriate HPV-vaccines for this region.
A relatively short 1-year time period, lack of adequate funding, low turnover of operative cases in the hospital, in availability of colposcopic facilities causing a possible sampling error and consequent false-negative diagnosis, occasional per- or post-operative loss of samples due to incorrect handling of tissue required for molecular studies and last, but not the least, inadequate managerial and infrastructural facilities all were limitations of the study and contributed to reduction in sample size.
Well-funded, larger sample-sized, case-control population based studies, which was beyond the scope of our study, could further validate present findings and put more insight on HPV genotype distribution and to develop a cost-effective, efficacious HPV vaccine for this region.
| References|| |
|1.||Bharadwaj M, Hussain S, Nasare V, Das BC. HPV & HPV vaccination: Issues in developing countries. Indian J Med Res 2009;130:327-33. |
|2.||zur Hausen H. Papillomaviruses causing cancer: Evasion from host-cell control in early events in carcinogenesis. J Natl Cancer Inst 2000;92:690-8. |
|3.||Rebolj M, Lynge E, Bonde J. Human papillomavirus testing and genotyping in cervical screening. Expert Rev Anticancer Ther 2011;11:1023-31. |
|4.||Shukla S, Bharti AC, Mahata S, Hussain S, Kumar R, Hedau S, et al. Infection of human papillomaviruses in cancers of different human organ sites. Indian J Med Res 2009;130:222-33. |
|5.||Böcker W. WHO classification of breast tumors and tumors of the female genital organs: Pathology and genetics. Verh Dtsch Ges Pathol 2002;86:116-9. |
|6.||Kurman RJ. Blaustein's Pathology of the Female Genital Tract. 5 th ed. Springer-Verlag, New York Inc., USA; 2002. |
|7.||Wei YC, Chou YS, Chu TY. Detection and typing of minimal human papillomavirus DNA in plasma. Int J Gynaecol Obstet 2007;96:112-6. |
|8.||Jenson AB, Rosenthal JD, Olson C, Pass F, Lancaster WD, Shah K. Immunologic relatedness of papillomaviruses from different species. J Natl Cancer Inst 1980;64:495-500. |
|9.||Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol 1990;215:403-10. |
|10.||Bhatla N, Joseph E. Cervical cancer prevention & the role of human papillomavirus vaccines in India. Indian J Med Res 2009;130:334-40. |
|11.||Sowjanya AP, Jain M, Poli UR, Padma S, Das M, Shah KV, et al. Prevalence and distribution of high-risk human papilloma virus (HPV) types in invasive squamous cell carcinoma of the cervix and in normal women in Andhra Pradesh, India. BMC Infect Dis 2005;5:116. |
|12.||Bancroft JD. Observations on the effect on histochemical reactions of different processing methods. J Med Lab Technol 1966;23:105-8. |
|13.||Bauer HM, Ting Y, Greer CE, Chambers JC, Tashiro CJ, Chimera J, et al. Genital human papillomavirus infection in female university students as determined by a PCR-based method. JAMA 1991;265:472-7. |
|14.||Sanger F, Donelson JE, Coulson AR, Kössel H, Fischer D. Determination of a nucleotide sequence in bacteriophage f1 DNA by primed synthesis with DNA polymerase. J Mol Biol 1974;90:315-33. |
|15.||Feoli-Fonseca JC, Oligny LL, Brochu P, Simard P, Falconi S, Yotov WV. Human papillomavirus (HPV) study of 691 pathological specimens from Quebec by PCR-direct sequencing approach. J Med Virol 2001;63:284-92. |
|16.||Hussain SK, Madeleine MM, Johnson LG, Du Q, Malkki M, Wilkerson HW, et al. Cervical and vulvar cancer risk in relation to the joint effects of cigarette smoking and genetic variation in interleukin 2. Cancer Epidemiol Biomarkers Prev 2008;17:1790-9. |
|17.||Khair MM, Mzibri ME, Mhand RA, Benider A, Benchekroun N, Fahime EM, et al. 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-84. |
|18.||Alfsen GC, Thoresen SO, Kristensen GB, Skovlund E, Abeler VM. Histopathologic subtyping of cervical adenocarcinoma reveals increasing incidence rates of endometrioid tumors in all age groups: A population based study with review of all nonsquamous cervical carcinomas in Norway from 1966 to 1970, 1976 to 1980, and 1986 to 1990. Cancer 2000;89:1291-9. |
|19.||Agarossi A, Ferrazzi E, Parazzini F, Perno CF, Ghisoni L. Prevalence and type distribution of high-risk human papillomavirus infection in women undergoing voluntary cervical cancer screening in Italy. J Med Virol 2009;81:529-35. |
|20.||Yemelyanova A, Vang R, Seidman JD, Gravitt PE, Ronnett BM. Endocervical adenocarcinomas with prominent endometrial or endomyometrial involvement simulating primary endometrial carcinomas: Utility of HPV DNA detection and immunohistochemical expression of p16 and hormone receptors to confirm the cervical origin of the corpus tumor. Am J Surg Pathol 2009;33:914-24. |
|21.||Mariani L, Monfulleda N, Alemany L, Vizza E, Marandino F, Vocaturo A, et al. Human papillomavirus prevalence and type-specific relative contribution in invasive cervical cancer specimens from Italy. BMC Cancer 2010;10:259. |
|22.||Ursin G, Peters RK, Henderson BE, d'Ablaing G 3 rd , Monroe KR, Pike MC. Oral contraceptive use and adenocarcinoma of cervix. Lancet 1994;344:1390-4. |
|23.||Zuna RE, Allen RA, Moore WE, Mattu R, Dunn ST. Comparison of human papillomavirus genotypes in high-grade squamous intraepithelial lesions and invasive cervical carcinoma: Evidence for differences in biologic potential of precursor lesions. Mod Pathol 2004;17:1314-22. |
|24.||Denehy TR, Gregori CA, Breen JL. Endocervical curettage, cone margins, and residual adenocarcinoma in situ of the cervix. Obstet Gynecol 1997;90:1-6. |
|25.||Gharizadeh B, Ghaderi M, Donnelly D, Amini B, Wallin KL, Nyrén P. Multiple-primer DNA sequencing method. Electrophoresis 2003;24:1145-51. |
|26.||Hanahan D, Weinberg RA. Hallmarks of cancer: The next generation. Cell 2011;144:646-74. |
|27.||Evans MF, Adamson CS, Simmons-Arnold L, Cooper K. Touchdown General Primer (GP5+/GP6+) PCR and optimized sample DNA concentration support the sensitive detection of human papillomavirus. BMC Clin Pathol 2005;5:10. |
|28.||Bhatla N, Dar L, Rajkumar Patro A, Kumar P, Pati SK, Kriplani A, et al. Human papillomavirus-type distribution in women with and without cervical neoplasia in north India. Int J Gynecol Pathol 2008;27:426-30. |
|29.||Sandri MT, Riggio D, Salvatici M, Passerini R, Zorzino L, Boveri S, et al. Typing of human papillomavirus in women with cervical lesions: Prevalence and distribution of different genotypes. J Med Virol 2009;81:271-7. |
|30.||Cheung JL, Cheung TH, Tang JW, Chan PK. Increase of integration events and infection loads of human papillomavirus type 52 with lesion severity from low-grade cervical lesion to invasive cancer. J Clin Microbiol 2008;46:1356-62. |
|31.||Tachezy R, Smahelova J, Salakova M, Arbyn M, Rob L, Skapa P, et al. Human papillomavirus genotype distribution in Czech women and men with diseases etiologically linked to HPV. PLoS One 2011;6:e21913. |
|32.||Venturoli S, Ambretti S, Cricca M, Leo E, Costa S, Musiani M, et al. Correlation of high-risk human papillomavirus genotypes persistence and risk of residual or recurrent cervical disease after surgical treatment. J Med Virol 2008;80:1434-40. |
|33.||Hadzisejdæ I, Kraseviæ M, Haller H, Grahovac B. Distribution of human papillomavirus types in different histological subtypes of cervical adenocarcinoma. Coll Antropol 2007;31 Suppl 2:97-102. |
|34.||Sigurdsson K, Taddeo FJ, Benediktsdottir KR, Olafsdottir K, Sigvaldason H, Oddsson K, et al. HPV genotypes in CIN 2-3 lesions and cervical cancer: A population-based study. Int J Cancer 2007;121:2682-7. |
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9], [Figure 10], [Figure 11], [Figure 12], [Figure 13]
[Table 1], [Table 2]