Ecological diversity profiles of non-vaccine-targeted HPVs after gender-based community vaccination efforts

The long-term effect of population-level human papillomavirus (HPV) vaccination on the viral ecology of the untargeted HPVs is poorly understood. We performed an 8-year follow-up of 33 communities randomized to gender-neutral HPV16/18 vaccination, girls-only HPV16/18 vaccination, and control communities without HPV vaccination. The 1992/93 and 1994 birth cohorts were invited in school years 2007/8 and 2008/9. Follow-up cervico-vaginal sampling at 18 and 22 years of age, 4 and 8 years post-vaccination, respectively, were attended by 11,396 and 5,602 participants. HPV types 6/11/16/18/31/33/35/39/45/51/52/56/58/59/66/68 were genotyped and used for the community-level ecological diversity estimations. Gender-neutral vaccination communities with a stronger herd immunity than girls-only vaccination communities show a significantly increased HPV α-diversity (p = 1.1 × 10-8) from 4 to 8 years post-vaccination, despite the clearance of the vaccine-targeted HPVs in these communities. This likely sign of niche occupation by the non-vaccine-targeted HPVs will potentially affect the future cervical cancer screening programs but should not interfere with the WHO mission to eliminate cervical cancer.

Assessing the risk of cervical neoplasia in the post-HPV vaccination era

This review is based on the recent EUROGIN scientific session: "Assessing risk of cervical cancer in the post-vaccination era," which addressed the demands of cervical intraepithelial neoplasia (CIN)/squamous intraepithelial lesion (SIL) triage now that the prevalence of vaccine-targeted oncogenic high-risk (hr) human papillomaviruses (HPVs) is decreasing. Change in the prevalence distribution of oncogenic HPV types that follows national HPV vaccination programs is setting the stage for loss of positive predictive value of conventional but possibly also new triage modalities. Understanding the contribution of the latter, most notably hypermethylation of cellular and viral genes in a new setting where most oncogenic HPV types are no longer present, requires studies on their performance in vaccinated women with CIN/SIL that are associated with nonvaccine HPV types. Lessons learned from this research may highlight the potential of cervical cells for risk prediction of all women's cancers.

Biophysical Properties of Bifunctional Phage-Biosensor

Biosensor research is a swiftly growing field for developing rapid and precise analytical devices for biomedical, pharmaceutical, and industrial use and beyond. Herein, we propose a phage-based biosensor method to develop a sensitive and specific system for biomedical detection. Our method is based on in vitro selected phages and their interaction with the targeted analytes as well as on optical properties that change according to the concentration of the model analyte. The green fluorescent protein (GFP) was chosen as our model analyte as it has its own well-known optical properties. Brilliant green was used as a reporter component for the sensor. Its presence enables a color intensity (absorbance) change when the analyte is present in the solution. Furthermore, the reporter dye functioned as a quencher for an additional lanthanide label in our assay. It mediated the specific phage-derived interference in the signal measured with the time-resolved luminescence. Most importantly, our results confirmed that the presented bifunctional phage with its liquid crystal properties enabled the measurement of GFP in a concentration-dependent, quantitative manner with a limit of detection of 0.24 µg/mL. In the future, our novel method to develop phage-based biosensors may provide highly sensitive and specific biosensors for biomedical or otherwise-relevant targets.

Drug-induced tumor-specific cytotoxicity in a whole tissue ex vivo model of human pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer. PDAC has a dismal prognosis and an inherent resistance to cytostatic drugs. The lack of reliable experimental models is a severe limitation for drug development targeting PDAC. We have employed a whole tissue ex vivo culture model to explore the effect of redox-modulation by sodium selenite on the viability and growth of PDAC. Drug-resistant tumors are more vulnerable to redox-active selenium compounds because of high metabolic activity and redox imbalance. Sodium selenite efficiently and specifically reduced PDAC cell viability (p <0.02) (n=8) and decreased viable de novo tumor cell outgrowth (p<0.05) while preserving non-neoplastic tissues. Major cellular responses (damaged tumor cells > 90%, tumor regression grades III-IV according to Evans) were observed for sodium selenite concentrations between 15-30 µM. Moreover, selenium levels used in this study were significantly below the previously reported maximum tolerated dose for humans. Transcriptome data analysis revealed decreased expression of genes known to drive PDAC growth and metastatic potential (CEMIP, DDR2, PLOD2, P4HA1) while the cell death-inducing genes (ATF3, ACHE) were significantly upregulated (p<0.0001). In conclusion, we report that sodium selenite has an extraordinary efficacy and specificity against drug-resistant pancreatic cancer in an organotypic slice culture model. Our ex vivo organotypic tissue slice culture model can be used to test a variety of drug candidates for swift and reliable drug responses to individual PDAC cases.

OUP accepted manuscript

For head-to-head comparison of human papillomavirus (HPV) antibody levels induced by different vaccines, 25-year-old vaccine-naive women were given either the bivalent (n = 188) or the nonavalent HPV vaccine (n = 184). Six months after vaccination antibodies against pseudovirions from 17 different HPV types (HPV6/11/16/18/31/33/35/39/45/51/52/56/58/59/66/68/73) were measured. Antibodies against HPV16/18 were higher after bivalent HPV vaccination (mean international units [IU] 1140.1 and 170.5 for HPV16 and 18, respectively) than after nonavalent vaccination (265.1 and 22.3 IUs, respectively). The bivalent vaccine commonly induced antibodies against the nonvaccine HPV types 31/33/35/45 or 58. The nonavalent vaccine induced higher antibodies against HPV6/11/31/33/45/52/58 and 35.

Potential SARS-CoV-2 infectiousness among asymptomatic healthcare workers

A majority of SARS-CoV-2 infections are transmitted from a minority of infected subjects, some of which may be symptomatic or pre-symptomatic. We aimed to quantify potential infectiousness among asymptomatic healthcare workers (HCWs) in relation to prior or later symptomatic disease. We previously (at the onset of the SARS-CoV-2 epidemic) performed a cohort study of SARS-CoV-2 infections among 27,000 healthcare workers (HCWs) at work in the capital region of Sweden. We performed both SARS-CoV-2 RT-PCR and serology. Furthermore, the cohort was comprehensively followed for sick leave, both before and after sampling. In the present report, we used the cohort database to quantify potential infectiousness among HCWs at work. Those who had sick leave either before or after sampling were classified as post-symptomatic or pre-symptomatic, whereas the virus-positive subjects with no sick leave were considered asymptomatic. About 0.2% (19/9449) of HCW at work were potentially infectious and pre-symptomatic (later had disease) and 0.17% (16/9449) were potentially infectious and asymptomatic (never had sick leave either before nor after sampling). Thus, 33% and 28% of all the 57 potentially infectious subjects were pre-symptomatic or asymptomatic, respectively. When a questionnaire was administered to HCWs with past infection, only 10,5% of HCWs had had no indication at all of having had SARS-CoV-2 infection ("truly asymptomatic"). Our findings provide a unique quantification of the different groups of asymptomatic, potentially infectious HCWs.

Human exposome assessment platform

The Human Exposome Assessment Platform (HEAP) is a research resource for the integrated and efficient management and analysis of human exposome data. The project will provide the complete workflow for obtaining exposome actionable knowledge from population-based cohorts. HEAP is a state-of-the-science service composed of computational resources from partner institutions, accessed through a software framework that provides the world's fastest Hadoop platform for data warehousing and applied artificial intelligence (AI). The software, will provide a decision support system for researchers and policymakers. All the data managed and processed by HEAP, together with the analysis pipelines, will be available for future research. In addition, the platform enables adding new data and analysis pipelines. HEAP's final product can be deployed in multiple instances to create a network of shareable and reusable knowledge on the impact of exposures on public health.

Risk for SARS-CoV-2 infection in healthcare workers outside hospitals: A real-life immuno-virological study during the first wave of the COVID-19 epidemic

OBJECTIVES

Most COVID-19 related infections and deaths may occur in healthcare outside hospitals. Here we explored SARS-CoV-2 infections among healthcare workers (HCWs) in this setting.

DESIGN

All healthcare providers in Stockholm, Sweden were asked to recruit HCWs at work for a study of past or present SARS-CoV-2 infections among HCWs. Study participants This study reports the results from 839 HCWs, mostly employees of primary care centers, sampled in June 2020.

RESULTS

SARS-CoV-2 seropositivity was found among 12% (100/839) of HCWs, ranging from 0% to 29% between care units. Seropositivity decreased by age and was highest among HCWs <40 years of age. Within this age group there was 19% (23/120) seropositivity among women and 11% (15/138) among men (p<0.02). Current infection, as measured using PCR, was found in only 1% and the typical testing pattern of pre-symptomatic potential "superspreaders" found in only 2/839 subjects.

CONCLUSIONS

Previous SARS-CoV-2 infections were common among younger HCWs in this setting. Pre-symptomatic infection was uncommon, in line with the strong variability in SARS-CoV-2 exposure between units. Prioritizing infection prevention and control including sufficient and adequate personal protective equipment, and vaccination for all HCWs are important to prevent nosocomial infections and infections as occupational injuries during an ongoing pandemic.

How infectious diseases arrived in the colonial Americas

Analysis of viral DNA from human remains suggests that the transatlantic slave trade may have introduced new pathogens that contributed to the devastating disease outbreaks in colonial Mexico.

Differences in risk for SARS-CoV-2 infection among healthcare workers

•

Clustering of SARS-CoV-2 exposure among HCWs in specific hospital wards was found.

•

SARS-CoV-2 infection was most common among nurse assistants in COVID-19 wards.

•

HCWs in wards with infectious diseases experience did not have increased risk.

Healthcare workers (HCWs) are a risk group for SARS-CoV-2 infection, but which healthcare work that conveys risk and to what extent such risk can be prevented is not clear. Starting on April 24th, 2020, all employees at work (n = 15,300) at the Karolinska University Hospital, Stockholm, Sweden were invited and 92% consented to participate in a SARS-CoV-2 cohort study. Complete SARS-CoV-2 serology was available for n = 12,928 employees and seroprevalences were analyzed by age, sex, profession, patient contact, and hospital department. Relative risks were estimated to examine the association between type of hospital department as a proxy for different working environment exposure and risk for seropositivity, adjusting for age, sex, sampling week, and profession. Wards that were primarily responsible for COVID-19 patients were at increased risk (adjusted OR 1.95 (95% CI 1.65–2.32) with the notable exception of the infectious diseases and intensive care units (adjusted OR 0.86 (95% CI 0.66–1.13)), that were not at increased risk despite being highly exposed. Several units with similar types of work varied greatly in seroprevalences. Among the professions examined, nurse assistants had the highest risk (adjusted OR 1.62 (95% CI 1.38–1.90)). Although healthcare workers, in particular nurse assistants, who attend to COVID-19 patients are a risk group for SARS-CoV-2 infection, several units caring for COVID-19 patients had no excess risk. Large variations in seroprevalences among similar units suggest that healthcare work-related risk of SARS-CoV-2 infection may be preventable.

High Amounts of SARS-CoV-2 Precede Sickness Among Asymptomatic Health Care Workers

Background

Whether SARS-CoV-2 positivity among symptomatic subjects reflects past or future disease may be difficult to ascertain.

Methods

We tested a cohort of 9449 employees at work at the Karolinska University Hospital, Stockholm, Sweden for SARS-CoV-2 RNA and antibodies, linked the testing results to sick leave records and determined associations with past or future sick leave using multinomial logistic regression. ClinicalTrials.gov NCT04411576.

Results

Subjects with high amounts of SARS-CoV-2 virus, as indicated by the Cycle threshold (Ct) value in the PCR, had the highest risk for sick leave in the two weeks after testing (OR 11.97 (CI 95% 6.29-22.80)) whereas subjects with low amounts of virus had the highest risk for sick leave in the past three weeks before testing (OR 6.31 (4.38-9.08)). Only 2.5% of employees were SARS-CoV-2 positive while 10.5% were positive by serology and 1.2% were positive in both tests. Serology-positive subjects were not at excess risk for future sick leave (OR 1.06 (95% CI, 0.71-1.57)).

Conclusions

High amounts of SARS-CoV-2 virus, as determined using PCR Ct values, associates with development of sickness in the next few weeks. The results support the concept that PCR Ct may be informative when testing for SARS-CoV-2 is performed.

Estimating Total Excess Mortality During a Coronavirus Disease 2019 Outbreak in Stockholm, Sweden

Total excess mortality peaked during a coronavirus disease 2019 (COVID-19) outbreak in Stockholm, but 25% of these deaths were not recognized as COVID-19 related nor occurred in hospitals. Estimate of total excess mortality may give a more comprehensive picture of the total disease burden during a COVID-19 outbreak, and may facilitate managing future outbreaks.

Vaccination With Moderate Coverage Eradicates Oncogenic Human Papillomaviruses If a Gender-Neutral Strategy Is Applied

BACKGROUND

Human papillomavirus (HPV) vaccination of girls with very high (>90%) coverage has the potential to eradicate oncogenic HPVs, but such high coverage is hard to achieve. However, the herd effect (HE) depends both on the HPV type and the vaccination strategy.

METHODS

We randomized 33 Finnish communities into gender-neutral HPV16/18 vaccination, girls-only HPV16/18 vaccination, and hepatitis B virus vaccination arms. In 2007-2010, 11 662 of 20 513 of 40 852 of 39 420 resident boys/girls from 1992 to 1995 birth cohorts consented. In 2010-2014, cervicovaginal samples from vaccinated and unvaccinated girls at age 18.5 years were typed for HPV6/11/16/18/31/33/35/39/45/51/52/56/58/59/66/68. Vaccine efficacy for vaccinated girls, HE for unvaccinated girls, and the protective effectiveness (PE) for all girls were estimated. We extended the community-randomized trial results about vaccination strategy with mathematical modeling to assess HPV eradication.

RESULTS

The HE and PE estimates in the 1995 birth cohort for HPV18/31/33 were significant in the gender-neutral arm and 150% and 40% stronger than in the girls-only arm. Concordantly, HPV18/31/33 eradication was already predicted in adolescents/young adults in 20 years with 75% coverage of gender-neutral vaccination. With the 75% coverage, eventual HPV16 eradication was also predicted, but only with the gender-neutral strategy.

CONCLUSIONS

Gender-neutral vaccination is superior for eradication of oncogenic HPVs.

Long-term follow-up of human papillomavirus type replacement among young pregnant Finnish females before and after a community-randomised HPV vaccination trial with moderate coverage

Large scale human papillomavirus (HPV) vaccination against the most oncogenic high-risk human papillomavirus (HPV) types 16/18 is rapidly reducing their incidence. However, attempts at assessing if this leads to an increase of nonvaccine targeted HPV types have been hampered by several limitations, such as the inability to differentiate secular trends. We performed a population-based serological survey of unvaccinated young women over 12 years. The women were under 23-years-old, residents from 33 communities which participated in a community-randomised trial (CRT) with approximately 50% vaccination coverage. Serum samples were retrieved pre-CRT and post-CRT implementation. Seropositivity to 17 HPV types was assessed. HPV seroprevalence ratios (PR) comparing the postvaccination to prevaccination era were estimated by trial arm. This was also assessed among the sexual risk-taking core group, where type replacement may occur more rapidly. In total, 8022 serum samples from the population-based Finnish Maternity Cohort were retrieved. HPV types 16/18 showed decreased seroprevalence among the unvaccinated in communities only after gender-neutral vaccination (PR_16/18A = 0.8, 95% CI 0.7-0.9). HPV6/11 and HPV73 were decreased after gender-neutral vaccination (PR_6/11A = 0.8, 95% CI 0.7-0.9, PR_73A = 0.7, 95% CI 0.6-0.9, respectively) and girls-only vaccination (PR_6/11B = 0.8, 95% CI 0.7-0.9, PR_73B = 0.9, 95% CI 0.8-1.0). HPV68 alone was increased but only after girls-only vaccination (PR_68B = 1.3, 95% CI 1.0-1.7, PR_core68B = 2.8, 95% CI 1.2-6.3). A large-scale, long-term follow-up found no type replacement in the communities with the strongest reduction of vaccine HPV types. Limited evidence for an increase in HPV68 was restricted to girls-only vaccinated communities and may have been due to secular trends (ClinicalTrials.gov number: NCT00534638).

Whole-genome sequence analysis of a Pan African set of samples reveals archaic gene flow from an extinct basal population of modern humans into sub-Saharan populations

Background

Population demography and gene flow among African groups, as well as the putative archaic introgression of ancient hominins, have been poorly explored at the genome level.

Results

Here, we examine 15 African populations covering all major continental linguistic groups, ecosystems, and lifestyles within Africa through analysis of whole-genome sequence data of 21 individuals sequenced at deep coverage. We observe a remarkable correlation among genetic diversity and geographic distance, with the hunter-gatherer groups being more genetically differentiated and having larger effective population sizes throughout most modern-human history. Admixture signals are found between neighbor populations from both hunter-gatherer and agriculturalists groups, whereas North African individuals are closely related to Eurasian populations. Regarding archaic gene flow, we test six complex demographic models that consider recent admixture as well as archaic introgression. We identify the fingerprint of an archaic introgression event in the sub-Saharan populations included in the models (~ 4.0% in Khoisan, ~ 4.3% in Mbuti Pygmies, and ~ 5.8% in Mandenka) from an early divergent and currently extinct ghost modern human lineage.

Conclusion

The present study represents an in-depth genomic analysis of a Pan African set of individuals, which emphasizes their complex relationships and demographic history at population level.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1684-5) contains supplementary material, which is available to authorized users.

Occurrence of human papillomavirus (HPV) type replacement by sexual risk-taking behaviour group: Post-hoc analysis of a community randomized clinical trial up to 9 years after vaccination (IV)

Oncogenic non-vaccine human papillomavirus (HPV) types may conceivably fill the vacated ecological niche of the vaccine types. The likelihood of this may differ by the risk of acquiring HPV infections. We examined occurrence of HPV types among vaccinated and unvaccinated subgroups of 1992-1994 birth cohorts with differing acquisition risks up to 9 years post-implementation of HPV vaccination in 33 Finnish communities randomized to: Arm A (gender-neutral HPV16/18 vaccination), Arm B (girls-only HPV16/18 vaccination and hepatitis B-virus (HBV) vaccination of boys), and Arm C (gender-neutral HBV vaccination). Out of 1992-1994 born resident boys (31,117) and girls (30,139), 8,618 boys and 15,615 girls were vaccinated, respectively, with 20-30% and 50% coverage in 2007-2009. In 2010-2013, 8,868 HPV16/18 and non-HPV vaccinated females, and in 2014-2016, 5,574 originally or later (2010-2013) HPV16/18 vaccinated females attended two cervical sampling visits, aged 18.5 and 22-years. The samples were typed for HPV6/11/16/18/31/33/35/39/45/51/52/56/58/59/66/68 using PCR followed by MALDI-TOF MS. HPV prevalence ratios (PR) between Arms A/B vs. C were calculated for Chlamydia trachomatis positives (core-group), and negatives (general population minus core group). At both visits the vaccine-protected HPV type PRs did not significantly differ between the core-group and non-core group. Among the vaccinated 18-year-olds, HPV51 occurrence was overall somewhat increased (PR_core = 1.4, PR_non-core. = 1.4) whereas the HPV52 occurrence was increased in the core-group only (PR_core = 2.5, PR_non-core = 0.8). Among the non-HPV vaccinated 18-year-olds, the HPV51/52 PRs were higher in the core-group (PR_core = 3.8/1.8, PR_non-core = 1.2/1.1). The 22-year-olds yielded no corresponding observations. Monitoring of the sexual risk-taking core-group may detect early tendencies for HPV type replacement.

Distinct geographic clustering of oncogenic human papillomaviruses multiple infections in cervical cancers: Results from a worldwide cross-sectional study

Coinfections by multiple Human Papillomaviruses (HPVs) are observed in approximately 6-8% of invasive cervical cancer (ICC) cases worldwide. But neither the presence of persistent HPVs coinfections nor their etiological role in the development of ICC is well understood. Cervical HPVs coinfections have been observed randomly, mostly in women with preneoplastic lesions, and only few studies have globally analyzed ICC cases. Here we explored the HPVs multiple infection patterns in a large worldwide sample of cross-sectional ICC cases. Paraffin-embedded ICC biopsy samples were tested using stringent HPV genotyping. Logistic regression models were used to identify the most likely pairwise HPV types in multiple infections. Multivariate analysis was applied to detect significant HPV coinfection patterns beyond pairwise HPVs comparison. Among 8780 HPV DNA-positive ICC cases worldwide, 6.7% (N = 587) contained multiple HPVs. Pairwise analysis revealed that HPV16|74, HPV31|33, HPV31|44, HPV33|44 and HPV45|70 pairs were significantly more frequently found together in multiple infections compared to any other HPV type combination, which supports the occasional role of Alpha-10 LR-HPVs in cervical cancers. In contrast, HPV16|31, HPV16|45, HPV16|51 and HPV18|HPV45 pairs were significantly less frequently found together than with any other HPV pair combination. Multivariate analysis sustained the results and revealed for the first time a distinct coinfection pattern in African ICCs stemming from the clustering of oncogenic HPV51/35/18/52 coinfections in African women. We suggest that the differential geographic HPVs coinfections clustering observed might be compatible with a specific modulation of the natural history/oncogenic potential of particular HPVs multiple infections and warrant monitoring for post-vaccinated.

Transmission between Archaic and Modern Human Ancestors during the Evolution of the Oncogenic Human Papillomavirus 16

Every human suffers through life a number of papillomaviruses (PVs) infections, most of them asymptomatic. A notable exception are persistent infections by Human papillomavirus 16 (HPV16), the most oncogenic infectious agent for humans and responsible for most infection-driven anogenital cancers. Oncogenic potential is not homogeneous among HPV16 lineages, and genetic variation within HPV16 exhibits some geographic structure. However, an in-depth analysis of the HPV16 evolutionary history was still wanting. We have analyzed extant HPV16 diversity and compared the evolutionary and phylogeographical patterns of humans and of HPV16. We show that codivergence with modern humans explains at most 30% of the present viral geographical distribution. The most explanatory scenario suggests that ancestral HPV16 already infected ancestral human populations and that viral lineages co-diverged with the hosts in parallel with the split between archaic Neanderthal-Denisovans and ancestral modern human populations, generating the ancestral HPV16A and HPV16BCD viral lineages, respectively. We propose that after out-of-Africa migration of modern human ancestors, sexual transmission between human populations introduced HPV16A into modern human ancestor populations. We hypothesize that differential coevolution of HPV16 lineages with different but closely related ancestral human populations and subsequent host-switch events in parallel with introgression of archaic alleles into the genomes of modern human ancestors may be largely responsible for the present-day differential prevalence and association with cancers for HPV16 variants.

Squamous intraepithelial lesions of the anal squamocolumnar junction: Histopathological classification and HPV genotyping

Background

Human papillomavirus (HPV)-related anal cancer lesions are often found adjacent to the squamocolumnar junction (SCJ). We have assessed the histopathology and associated HPV genotypes in anal SCJ lesions in surgically excised anal warts in HIV-negative and –positive patients.

Methods

Histopathology identified 47 squamous intraepithelial lesions (SILs) adjacent to the SCJ amongst a total of 145 cases of clinically diagnosed anal condylomata. The anal SCJ lesions were further analyzed with p16, CK7 and p63 immunohistochemistry and HPV genotyping.

Results

Sixteen (16/47) of the excised anal wart lesions contained HSIL; Three were HSIL and exclusively associated with oncogenic HPVs. A further thirteen (13/47) were mixed lesions. Of these eight were HSILs with LSIL and six were HSILs with papillary immature metaplasia (PIM); Ten of the mixed lesions were associated with one or more oncogenic HPVs, while three cases were exclusively associated with HPV6.

Conclusions

Clinically diagnosed anal warts cannot be assumed to be limited to low-grade lesions as anal warts of the SCJ often show heterogeneous lesions, with coexistence of LSIL, PIM, and HSIL. Lesions showing PIM, however, may mimic HSIL, because they are hypercellular, but lack the nuclear atypia and conspicuous mitotic activity of HSIL; and are p16 negative.

•

Anal warts of the SCJ are often heterogeneous; with LSIL, including PIM and HSIL.

•

Anal HSILs only with HPV6 may indicate that LR-HPVs have some role in oncogenesis.

•

Recognition of PIM is important given its potential for being confused with HSIL.

•

Clinical diagnosis of anal warts cannot be assumed to be limited to LSIL.

•

Microscopic examination of anal warts is encouraged with appropriate immunostains.

HPV16 variants distribution in invasive cancers of the cervix, vulva, vagina, penis, and anus

Human papillomavirus (HPV)16 is the most oncogenic human papillomavirus, responsible for most papillomavirus‐induced anogenital cancers. We have explored by sequencing and phylogenetic analysis the viral variant lineages present in 692 HPV16‐monoinfected invasive anogenital cancers from Europe, Asia, and Central/South America. We have assessed the contribution of geography and anatomy to the differential prevalence of HPV16 variants and to the nonsynonymous E6 T350G polymorphism. Most (68%) of the variance in the distribution of HPV16 variants was accounted for by the differential abundance of the different viral lineages. The most prevalent variant (above 70% prevalence) in all regions and in all locations was HPV16_A1‐3, except in Asia, where HPV16_A4 predominated in anal cancers. The differential prevalence of variants as a function of geographical origin explained 9% of the variance, and the differential prevalence of variants as a function of anatomical location accounted for less than 3% of the variance. Despite containing similar repertoires of HPV16 variants, we confirm the worldwide trend of cervical cancers being diagnosed significantly earlier than other anogenital cancers (early fifties vs. early sixties). Frequencies for alleles in the HPV16 E6 T350G polymorphism were similar across anogenital cancers from the same geographical origin. Interestingly, anogenital cancers from Central/South America displayed higher 350G allele frequencies also within HPV16_A1‐3 lineage compared with Europe. Our results demonstrate ample variation in HPV16 variants prevalence in anogenital cancers, which is partly explained by the geographical origin of the sample and only marginally explained by the anatomical location of the lesion, suggesting that tissue specialization is not essential evolutionary forces shaping HPV16 diversity in anogenital cancers.