Cholinergic signaling through nicotinic acetylcholine receptors stimulates the proliferation of cervical cancer cells: an explanation for the molecular role of tobacco smoking in cervical carcinogenesis?

We have analyzed the expression of mRNAs encoding nicotinic acetylcholine receptors (nAChRs) in CaSki, SiHa and HeLa cell lines, which are derived from two squamous and one adenocarcinoma of the cervix, respectively. We detected with reverse transcription-polymerase chain reaction mRNAs for ten of the 16 nAChR subunits, namely strong signals for alpha-5, alpha-7, alpha-9, beta-1 and epsilon, and weak signals for alpha-4, beta-2, beta-4, gamma and delta. We confirmed the translation of alpha-5 and beta-1, corresponding to the two strongest RNA signals, in SiHa and HeLa cells by Western blotting, and the localization of these proteins to the plasma membrane by immunofluorescence. The beta-1 subunit was detected membrane-associated in normal and neoplastic squamous epithelia of the cervix in situ, but appeared to be absent from the underlying mesenchyme and even from adjacent columnar epithelia. These observations suggest that normal and neoplastic cervical squamous epithelial cells express several combinations of the pentameric nAChRs. We also measured that the proliferation of SiHa and HeLa cells is stimulated by nicotine. This indicates that cholinergic signaling under normal physiological conditions and stimulated by nicotine in tobacco users affects epithelial homeostasis and neoplastic progression at the cervix in a way similar to the known effects on epithelia of the mouth, the airways and the lung. Since tobacco smoking is established as a risk factor in cervical carcinogenesis, and since nicotine and its derivatives become concentrated in cervical mucus, nAChR-dependent signaling is apparently an important molecular cofactor of human papillomavirus-dependent cervical carcinogenesis.

Human papillomavirus-16 and -18 in penile carcinomas: DNA methylation, chromosomal recombination and genomic variation.

Penile carcinomas are frequently associated with high risk human papillomavirus (HPV) types. Because little is known about the molecular biology of this association, we investigated three properties of HPV genomes in penile carcinomas from Brazilian patients: (i) HPV DNA methylation, (ii) junctions between HPV and cellular DNA and (iii) genomic variation. In cervical carcinogenesis, recombination between HPV and chromosomal DNA is frequent and likely necessary for progression, and DNA hypermethylation-specifically of the L1 gene-is a biomarker for cancerous progression. The same mechanisms apparently occur during penile carcinogenesis, because 95 HPV-16 molecules derived from 19 penile lesions had 58% of the CpGs in L1 and 22% in the 5' part of the long control region methylated, more than the percentages found in cervical carcinomas. In addition, 2 out of 3 HPV-18 infections, all present in double infections with HPV-16, showed L1 specific methylation typical of malignant cervical lesions. In 11 out of 15 HPV-16 lesions, we confirmed chromosomal integration by reverse ligation inverted PCR, while 4 samples had concatemeric integrations or episomes. Nine of 17 penile carcinomas contained HPV-16 AA variants, and 8 E variants. As AA variants are relatively rare in Brazilian cohorts of asymptomatic women, the high prevalence in penile carcinomas may indicate a higher risk of progression of AA lesions, as suspected for cervical infections. Our observations of frequent viral DNA methylation, chromosomal integration and the prevalence of high risk variants suggest that HPV-dependent carcinogenesis of the penis and cervix follows similar etiological and epidemiological parameters.

Prevalence and 3-year persistence of human papillomavirus serotypes in asymptomatic patients in Northern Mexico.

OBJECTIVE: To investigate clinical outcomes and 3-year persistence of human papillomavirus (HPV) infections among women in Mexico. METHODS: A prospective study enrolled sexually active women attending primary healthcare clinics in metropolitan Monterrey, Mexico, between June 3 and August 30, 2002. Baseline data were collected and participants underwent HPV screening. Patients with HPV infections were asked to attend a repeat screening appointment after 3 years, when the same screening data were gathered. Descriptive analyses were performed and the prevalence of cervical lesions and viral infections were examined. RESULTS: In total, 1188 patients who underwent initial HPV screening were included. Cervical lesions were detected in 5 (0.4%) patients and 239 (20.1%) patients had HPV infections; 129 (54.0%) of these patients attended 3-year follow-up. Among the 357 HPV serotypes identified, the most prevalent serotypes were HPV-59, HPV-52, HPV-16, and HPV-56, detected 62 (17.4%), 38 (10.6%), 27 (7.6%), and 18 (5.0%) times, respectively. Of the 129 patients attending 3-year follow-up, 104 (80.6%) were clear from HPV infections, 13 (10.1%) patients had persistent HPV infections, and 12 (9.3%) had HPV infections with different HPV types. CONCLUSIONS: The HPV prevalence was 20.1% in the present study; the most prevalent infections were HPV-59, HPV-52, HPV-16, and HPV-56. At 3-year follow-up, 25 (19.4%) patients had HPV infections.

Methylation of human papillomavirus 16, 18, 31, and 45 L2 and L1 genes and the cellular DAPK gene: Considerations for use as biomarkers of the progression of cervical neoplasia.

During progression of cervical cancer, human papillomavirus genomes and cellular tumor suppressor genes can become methylated. Toward a better understanding of these biomarkers, we studied 104 samples with HPV16, 18, 31, and 45 representing five pathological categories from asymptomatic infection to cancer. We grouped all samples by HPV type and pathology and measured the overall methylation of informative amplicons of HPV late genes and the cellular DAPK gene. Methylation of all four HPV types as well as of the DAPK gene is lowest in asymptomatic infection and increases successively in all four pathological categories during progression to cancer. 27 out of 28 cancer samples showed methylation both in the L2/L1 genes as well as in DAPK, but a much lower fraction in all other pathological categories. We discuss the problem to develop diagnostic tests based on complex methylation patterns that make it difficult to classify amplicons as "methylated" or "unmethylated".

Muscarinic cholinergic signaling in cervical cancer cells affects cell motility via ERK1/2 signaling.

AIMS: The etiology of cervical cancer depends primarily on infection with human papillomaviruses, but tobacco smoking is the most important behavioral risk factor for this cancer. Therefore, we have previously confirmed involvement of nicotinic acetylcholine receptors (nAChRs) in cervical cancer biology. In order to comprehensively evaluate the role of cholinergic signaling in cervical cells, we have addressed additional participation of muscarinic acetylcholine receptors (mAChRs). MAIN METHODS: We have studied the expression of mAChRs and cholinergic system components by reverse transcription PCR and Western blots, the motility of cervical cancer cells in cell culture, and the signaling from mAChRs via the ERK1/2 signaling pathway. KEY FINDINGS: The cervical cancer cells HeLa, SiHa and CaSki express four of the five mAChRs, M1, M3, M4, and M5, and the acetylcholine (ACh) synthesizing and degrading enzymes choline acetyltransferase (ChAT), acetylcholinesterase (AChE), and butyrylcholinesterase (BChE), and vesicular ACh transporter (VAChT). mAChR-dependent signaling induces cervical cell motility, which requires ERK1/2 activation, and could be abrogated by mAChR antagonists. SIGNIFICANCE: The epidemiological finding that tobacco smoke raises the prevalence of cervical cancer has led to analysis of the cholinergic signaling in cervical biology and carcinogenesis. Cervical cancer cells express several nAChRs and mAChRs, whose activation leads to changes of cellular properties such as increased motility and proliferation that favor a carcinogenic phenotype. The signaling involves intracellular phosphorylation cascades including ERK1/2.

Laser capture microdissection of cervical human papillomavirus infections: copy number of the virus in cancerous and normal tissue and heterogeneous DNA methylation.

Research on the pathogenicity of human papillomaviruses (HPVs) during cervical carcinogenesis often relies on the study of homogenized tissue or cultured cells. This approach does not detect molecular heterogeneities within the infected tissue. It is desirable to understand molecular properties in specific histological contexts. We asked whether laser capture microdissection (LCM) of archival cervical tumors in combination with real-time polymerase chain reaction and bisulfite sequencing permits (i) sensitive DNA diagnosis of small clusters of formalin-fixed cells, (ii) quantification of HPV DNA in neoplastic and normal cells, and (iii) analysis of HPV DNA methylation, a marker of tumor progression. We analyzed 26 tumors containing HPV-16 or 18. We prepared DNA from LCM dissected thin sections of 100 to 2000 cells, and analyzed aliquots corresponding to between nine and 70 cells. We detected nine to 630 HPV-16 genome copies and one to 111 HPV-18 genome copies per tumor cell, respectively. In 17 of the 26 samples, HPV DNA existed in histologically normal cells distant from the margins of the tumors, but at much lower concentrations than in the tumor, suggesting that HPVs can infect at low levels without pathogenic changes. Methylation of HPV DNA, a biomarker of integration of the virus into cellular DNA, could be measured only in few samples due to limited sensitivity, and indicated heterogeneous methylation patterns in small clusters of cancerous and normal cells. LCM is powerful to study molecular parameters of cervical HPV infections like copy number, latency and epigenetics.

Effects of cellular differentiation, chromosomal integration and 5-aza-2'-deoxycytidine treatment on human papillomavirus-16 DNA methylation in cultured cell lines.

Human papillomavirus-16 (HPV-16) genomes in cell culture and in situ are affected by polymorphic methylation patterns, which can repress the viral transcription. In order to understand some of the underlying mechanisms, we investigated changes of the methylation of HPV-16 DNA in cell cultures in response to cellular differentiation, to recombination with cellular DNA, and to an inhibitor of methylation. Undifferentiated W12E cells, derived from a precancerous lesion, contained extrachromosomal HPV-16 DNA with a sporadically methylated enhancer-promoter segment. Upon W12E cell differentiation, the viral DNA was demethylated, suggesting a link between differentiation and the epigenetic state of HPV-16 DNA. The viral genomes present in two W12I clones, in which individual copies of the HPV-16 genome have integrated into cellular DNA (type 1 integrants), were unmethylated, akin to that seen in the cervical carcinoma cell line SiHa (also a type 1 integrant). This finding is consistent with hypomethylation being necessary for continued viral gene expression. In contrast, two of three type 2 integrant W12I clones, containing concatemers of HPV-16 genomes integrated into the cellular DNA contained hypermethylated viral DNA, as observed in the cervical carcinoma cell line CaSki (also a type 2 integrant). A third, type 2, W12I clone, interestingly with fewer copies of the viral genome, contained unmethylated HPV-16 genomes. Epithelial differentiation of W12I clones did not lead to demethylation of chromosomally integrated viral genomes as was seen for extrachromosomal HPV-16 DNA in W12E clones. Hypomethylation of CaSki cells in the presence of the DNA methylation inhibitor 5-aza-2'-deoxycytidine reduced the cellular viability, possibly as a consequence of toxic effects of an excess of HPV-16 gene products. Our data support a model wherein (i) the DNA methylation state of extrachromosomal HPV16 replicons and epithelial differentiation are inversely coupled during the viral life cycle, (ii) integration of the viral genome into the host chromosome events leads to an alteration in methylation patterns on the viral genome that is dependent upon the type of integration event and possibly copy number, and (iii) integration universally results in the viral DNA becoming refractory to changes in methylation state upon cellular differentiation that are observed with extrachromosomal HPV-16 genomes.

Human papillomavirus-16 DNA methylation patterns support a causal association of the virus with oral squamous cell carcinomas.

Infection with human papillomavirus-16 (HPV-16) is the cause of most anogenital carcinomas. This virus is also detected in about 20% of all head and neck squamous cell carcinomas. While there is strong evidence for a causal etiological role in the case of tonsillar carcinomas, causal association with malignant lesions of the oral cavity is not yet conclusive. Our previous investigations of HPV-16 DNA methylation in anogenital sites have identified hypermethylation of the L1 gene and part of the long control region in many malignant lesions, but rarely in asymptomatic infections and low-grade precancerous lesions. Here, we report hypermethylation of this diagnostically important segment of the viral DNA in 10 out of 12 HPV-16 positive oral carcinomas from Mexican patients. These data indicate epigenetic changes of HPV-16 in oral carcinomas similar to those in anogenital carcinomas, suggesting carcinogenic processes under the influence of HPV-16 in most if not all of these oral malignant lesions.

Genome variation of human papillomavirus types: phylogenetic and medical implications.

Human papillomaviruses (HPVs) are described as "types" based on their genome sequences and identified by a number. For example, HPV-6 is associated with genital warts, and HPV-16 with anogenital cancers. The genomes of many HPV types have been reisolated, sequenced and compared to reference "prototypes" countless times by laboratories throughout the world. It was found that each HPV type occurs in the form of "variants", identified by about 2% nucleotide differences in most genes and 5% in less conserved regions. Less than 100 variants of any HPV type have been detected, a scenario that is very different from the quasi-species formed by many RNA viruses. The variants of each HPV type form phylogenetic trees, and variants from specific branches are often unique to specific ethnic groups. Immigrant populations contain, depending on their respective ethnic origins, mixtures of variants. The absence of HPV genomes intermediate to specific types show that all HPV types existed already when humans became a species. Consequently, humans had always suffered from lesions like anogenital cancer, genital warts and common warts. A growing number of epidemiological, etiological and molecular data suggest that variants of the same HPV type are biologically distinct and may confer differential pathogenic risks. Since the distribution of some variants of HPV-16 and 18 correlates with the distribution of human populations that have an increased risk to develop anogenital cancer, the study of HPV type variation may point to one of the reasons for the higher incidence rates of these lesions in specific cohorts.

Methylation of the human papillomavirus-18 L1 gene: a biomarker of neoplastic progression?

Epigenetic transcriptional regulation plays an important role in the life cycle of human papillomaviruses (HPVs) and the carcinogenic progression of anogenital HPV associated lesions. We performed a study designed to assess the methylation status of the HPV-18 genome, specifically of the late L1 gene, the adjacent long control region (LCR), and part of the E6 oncogene in cervical specimens with a range of pathological diagnoses. In asymptomatic infections and infections with precancerous (precursor) lesions, HPV-18 DNA was mostly unmethylated, with the exception of four samples where hypermethylation of L1 was detected. In contrast, L1 sequences were strongly methylated in all cervical carcinomas, while the LCR and E6 remained unmethylated. HeLa cells, derived from a cervical adenocarcinoma, contain chromosomally integrated HPV-18 genomes. We found that L1 is hypermethylated in these cells, while the LCR and E6 are unmethylated. Treatment of HeLa cells with the methylation inhibitor 5-Aza-2'-deoxycytidine (5-Aza-CdR) led to the expected reduction of L1 methylation. After removal of 5-Aza-CdR, L1 methylation resumed and exceeded pretreatment levels. Unexpectedly, the LCR and E6 also became methylated under these conditions, albeit at lower levels than L1. We hypothesize that L1 is preferentially methylated after integration of the HPV genome into the cellular DNA, possibly since linearization prohibits its normal transcription, while the enhancer and promoter may be protected from methylation by transcription factors. Since our data suggest that HPV-18 L1 methylation can only be detected in carcinomas, except in some few precancerous lesions and asymptomatic infections, L1 methylation may constitute a powerful molecular marker for detecting this important step of neoplastic progression.

Worldwide genomic diversity of the human papillomaviruses-53, 56, and 66, a group of high-risk HPVs unrelated to HPV-16 and HPV-18.

Among more than 200 human papillomavirus (HPV) types presumed to exist, 18 "high-risk" HPV types are frequently found in anogenital cancer. The best studied types are HPV-16 and 18, which are only distantly related to one another and form two separate phylogenetic branches, each including six closely related types. HPV-30, 53, 56, and 66 form a third phylogenetic branch unrelated to HPV-16 and 18. Worldwide comparison of HPV-16 and 18 isolates revealed a distribution of variant genomes that correlated with the geographic origin and the ethnicity of the infected cohort and led to the concept of unique African, European, Asian, and Native American HPV-16 and 18 variants. Here, we address the question whether similar phylogenies are found for HPV-53, 56, and 66 by determining the sequence of the long control regions (LCR) of these HPVs in samples from Europe, Asia, and Africa, and from immigrant societies in North and South America. Phylogenetic trees calculated from point mutations and a few insertions/deletions affecting 2-4.2% of the nucleotide sequences were distinct for each of the three HPVs and divergent from HPV-16 and 18. In contrast to the "star-phylogenies" formed by HPV-16 and 18 variants, 44 HPV-53 isolates represented nine variants, which formed two deep dichotomic branches reminiscent of the beginning split into two new taxa, as recently observed for subtypes of HPV-44 and 68. A total of 66 HPV-56 isolates represented 17 variants, which formed three branches preferentially containing European, Asian, and African variants. Variants of a fourth branch, deeply separated from the other three, were characterized by a 25 bp insertion and created a dichotomy rather than star-like phylogeny. As it contained isolates from cohorts in all continents, it may have evolved before the spread of humans into all continents. 18 of 31 HPV-66 isolates represented the prototype clone, which was found in all parts of the world, while the remaining 13 clones formed 11 branches without any geographic association. Our findings confirm the notion of a quantitatively limited genomic diversity of each HPV type with some correlation to the geographic origin of the sample. In addition, we observed in some variants of these three HPV types mutations that affect the amino acid sequence of the E6 oncoproteins and the L1 capsid protein, supporting the possibility of immunogenic and oncogenic diversity between variants of any HPV type.

Worldwide genomic diversity of the high-risk human papillomavirus types 31, 35, 52, and 58, four close relatives of human papillomavirus type 16.

Among the more than one hundred formally described human papillomavirus (HPV) types, 18 are referred to as high-risk HPV types due to their association with anogenital cancer. Despite pathogenic similarities, these types form three remotely related taxonomic groups. One of these groups is called HPV species 9 and is formed by HPV-16, the most common and best-studied type, together with HPV-31, -33, -35, -52, -58, and -67. Previous worldwide comparisons of HPV-16 samples showed about 2% nucleotide diversity between isolates, which were subsequently termed variants. The distribution of divergent variants has been found to correlate frequently with the geographic origin and the ethnicity of the infected patients and led to the concept of unique African, European, Asian, and Native American HPV-16 variants. In the current study, we address the question of whether geography and ethnicity also correlate with sequence variations found for HPV-31, -35, -52, and -58. This was done by sequencing the long control region in samples derived from Europe, Asia, and Africa, and from immigrant populations in North and South America. We observed maximal divergence between any two variants within each of these four HPV types ranging from 1.8 to 3.6% based on nucleotide exchanges and, occasionally, on insertions and deletions. Similar to the case with HPV-16, these mutations are not random but indicate a relationship between the variants in form of phylogenetic trees. An interesting example is presented by a 16-bp insert in select variants of HPV-35, which appears to have given rise to additional variants by nucleotide exchanges within the insert. All trees showed distinct phylogenetic topologies, ranging from dichotomic branching in the case of HPV-31 to star phylogenies of the other three types. No clear similarities between these types or between these types and HPV-16 exist. While variant branches in some types were specific for Europe, Africa, or East Asia, none of the four trees reflected human evolution and spread to the extent illustrated by HPV-16. One possible explanation is that the rare HPV types that we studied spread and thereby diversified more slowly than the more abundant HPV-16 and may have established much of today's variant diversity already before the worldwide spread of humans 100,000 years ago. Most variants had prototypic amino acid sequences within the E6 oncoprotein and a segment of the L1 capsid protein. Some had one, two, or three amino acid substitutions in these regions, which might indicate biological and pathogenic diversity between the variants of each HPV type.

Papillomavirus subtypes are natural and old taxa: phylogeny of human papillomavirus types 44 and 55 and 68a and -b.

A human papillomavirus (HPV) type is defined as an HPV isolate whose L1 gene sequence is at least 10% different from that of any other type, while a subtype is 2 to 10% different from any HPV type. In order to analyze the phylogeny behind the subtype definition, we compared 49 isolates of HPV type 44 (HPV-44) and its subtype HPV-55, previously misclassified as a separate type, and 41 isolates of the subtype pair HPV-68a and -b, sampled from cohorts in four continents. The subtypes of each pair are separated by deep dichotomic branching, and three of the four subtypes have evolved large phylogenetic clusters of genomic variants forming a "star" phylogeny, with some branches specific for ethnically defined cohorts. We conclude that subtypes of HPV types are natural and old taxa, equivalent to types, which either diverged more recently than types or evolved more slowly.

Conserved methylation patterns of human papillomavirus type 16 DNA in asymptomatic infection and cervical neoplasia.

DNA methylation contributes to the chromatin conformation that represses transcription of human papillomavirus type16 (HPV-16), which is prevalent in the etiology of cervical carcinoma. In an effort to clarify the role of this phenomenon in the regulation and carcinogenicity of HPV-16, 115 clinical samples were studied to establish the methylation patterns of the 19 CpG dinucleotides within the long control region and part of the L1 gene by bisulfite modification, PCR amplification, DNA cloning, and sequencing. We observed major heterogeneities between clones from different samples as well as between clones from individual samples. The methylation frequency of CpGs was measured at 14.5%. In addition, 0.21 and 0.23%, respectively, of the CpA and CpT sites, indicators of de novo methylation, were methylated. Methylation frequencies exceeded 30% in the CpGs overlapping with the L1 gene and were about 10% for most other positions. A CpG site located in the linker between two nucleosomes positioned over the enhancer and promoter of HPV-16 had minimal methylation. This region forms part of the HPV replication origin and is close to binding sites of master-regulators of transcription during epithelial differentiation. Methylation of most sites was highest in carcinomas, possibly due to tandem repetition and chromosomal integration of HPV-16 DNA. Methylation was lowest in dysplasia, likely reflecting the transcriptional activity in these infections. Our data document the efficient targeting of HPV genomes by the epithelial methylation machinery, possibly as a cellular defense mechanism, and suggest involvement of methylation in HPV oncogene expression and the early-late switch.

The human papillomavirus-18 genome is efficiently targeted by cellular DNA methylation.

Human papillomaviruses (HPVs) infect epithelia, including the simple and the squamous epithelia of the cervix, where they can cause cancer and precursor lesions. The molecular events leading from asymptomatic HPV infections to neoplasia are poorly understood. There is evidence that progression is modulated by transcriptional mechanisms that control HPV gene expression. Here, we report the frequent methylation of HPV-18 genomes in cell culture and in situ. DNA methylation is generally known to lead to transcriptional repression due to chromatin changes. We investigated two cell lines derived from cervical cancers, namely, C4-1, which contains one HPV-18 genome, and different clones of HeLa, with 50 HPV-18 genomes. By restriction cleavage, we detected strong methylation of the L1 gene and absence of methylation of parts of the long control region (LCR). A 3-kb segment of the HPV-18 genomes downstream of the oncogenes was deleted in both cell lines. Bisulfite sequencing showed that in C4-1 cells and two HeLa clones, 18 of the 19 CpG residues in the 1.2-kb terminal part of the L1 gene were methylated, whereas a third HeLa clone had only eight methylated CpG groups, indicating changes of the methylation pattern after the establishment of the HeLa cell line. In the same four clones, none of the 12 CpG residues that overlapped with the enhancer and promoter was methylated. In six HPV-18 containing cancers and five smears from asymptomatic patients, most of the CpG residues in the L1 gene were methylated. There was complete or partial methylation, respectively, of the HPV enhancer in three of the cancers, and lack of methylation in the remaining eight samples. The promoter sequences were methylated in three of the six cancers and four of the six smears, and unmethylated elsewhere. Our data show that epithelial cells efficiently target HPV-18 genomes for DNA methylation, which may affect late and early gene transcription.

Genomic diversity of human papillomavirus-16, 18, 31, and 35 isolates in a Mexican population and relationship to European, African, and Native American variants.

Cervical cancer, mainly caused by infection with human papillomaviruses (HPVs), is a major public health problem in Mexico. During a study of the prevalence of HPV types in northeastern Mexico, we identified, as expected from worldwide comparisons, HPV-16, 18, 31, and 35 as highly prevalent. It is well known that the genomes of HPV types differ geographically because of evolution linked to ethnic groups separated in prehistoric times. As HPV intra-type variation results in pathogenic differences, we analyzed genomic sequences of Mexican variants of these four HPV types. Among 112 HPV-16 samples, 14 contained European and 98 American Indian (AA) variants. This ratio is unexpected as people of European ethnicity predominate in this part of Mexico. Among 15 HPV-18 samples, 13 contained European and 2 African variants, the latter possibly due to migration of Africans to the Caribbean coast of Mexico. We constructed phylogenetic trees of HPV-31 and 35 variants, which have never been studied. Forty-six HPV-31 isolates from Mexico, Europe, Africa, and the United States (US) contained a total of 35 nucleotide exchanges in a 428-bp segment, with maximal distances between any two variants of 16 bp (3.7%), similar to those between HPV-16 variants. The HPV-31 variants formed two branches, one apparently the European, the other one an African branch. The European branch contained 13 of 29 Mexican isolates, the African branch 16 Mexican isolates. These may represent the HPV-31 variants of American Indians, as a 55% prevalence of African variants in Mexico seems incomprehensible. Twenty-seven HPV-35 samples from Mexico, Europe, Africa, and the US contained 11 mutations in a 893-bp segment with maximal distances between any two variants of only 5 mutations (0.6%), including a characteristic 16-bp insertion/deletion. These HPV-35 variants formed several phylogenetic clusters rather than two- or three-branched trees as HPV-16, 18, and 31. An HPV-35 variant typical for American Indians was not identifiable. Our research suggests type specific patterns of evolution and spread of HPV-16, 18, 31, and 35 both before and after the worldwide migrations of the last four centuries. The high prevalence of highly carcinogenic HPV-16 AA variants, and the extensive diversity of HPV-18, 31, and 35 variants with unknown pathogenic properties raise the possibility that HPV intra-type variation contributes to the high cervical cancer burden in Mexico.