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

The Lineage and Sublineage Investigation of Human Papillomavirus Type 16 in Tehran, Iran, During 2022-2023: A Cross-Sectional Study

Background and Aims

Diverse HPV 16 variants have been known to vary in geographic distribution and oncogenic potential. On this point, the nucleotide changes of the E6 gene were studied to find lineages of HPV 16 in normal, premalignant, and malignant stages of uterine cervical samples.

Methods

In this study, 120 HPV 16-infected samples were investigated using PCR and sequencing.

Results

In our samples, three lineages A, C, and D were found and lineage D was predominant (79.2%) followed by lineages A (20%) and C (0.8%). Concerning the association between histopathological stages and lineages, no statistically important differences were found.

Conclusion

Our finding revealed that two lineages A and D were circulating, with the dominancy of lineage, in Tehran, Iran. This finding emphasizes the geographic diversity of distinct HPV 16 (sub)lineages is different in the world.

Why HPV16? Why, now, HPV42? How the discovery of HPV42 in rare cancers provides an opportunity to challenge our understanding about the transition between health and disease for common members of the healthy microbiota

In 2022, a bioinformatic, agnostic approach identified HPV42 as causative agent of a rare cancer, later confirmed experimentally. This unexpected association offers an opportunity to reconsider our understanding about papillomavirus infections and cancers. We have expanded our knowledge about the diversity of papillomaviruses and the diseases they cause. Yet, we still lack answers to fundamental questions, such as what makes HPV16 different from the closely related HPV31 or HPV33; or why the very divergent HPV13 and HPV32 cause focal epithelial hyperplasia, while HPV6 or HPV42 do not, despite their evolutionary relatedness. Certain members of the healthy skin microbiota are associated to rare clinical conditions. We propose that a focus on cellular phenotypes, most often transient and influenced by intrinsic and extrinsic factors, may help understand the continuum between health and disease. A conceptual switch is required towards an interpretation of biology as a diversity of states connected by transition probabilities, rather than quasi-deterministic programs. Under this perspective, papillomaviruses may only trigger malignant transformation when specific viral genotypes interact with precise cellular states. Drawing on Canguilhem's concepts of normal and pathological, we suggest that understanding the transition between fluid cellular states can illuminate how commensal-like infections transition from benign to malignant.

Lineage and sublineage analysis of human papillomavirus type 58 in iranian women

BACKGROUND

Variant analysis of distinct HPV types is important from different aspects including epidemiology, pathogenicity, and evolution.

METHODS

For this reason, the full sequence of the E6 and E7 genes of HPV 58 was examined in 130 HPV 58-infected cervical samples using PCR and sequencing.

RESULTS

Our results revealed that three lineages A, B, and D were found in this study; among which the B lineage was more common (91.50%). About sublineages, all samples of the B lineage belonged to the B1 sublineage, and samples that were classified as the A and D lineages were found to belong to the A1 (0.77%), A2 (5.38%), A3 (1.50%), and D2 (0.77%) sublineages. No statistically significant differences were found between lineages and stages of disease or amino acid changes (P > 0.05).

CONCLUSION

Our results showed that lineage B, sublineage B1, was dominant in Iran. However, more studies with larger sample sizes from different parts of Iran are essential for assessing the pathogenicity risk of HPV 58 lineages in Iranian women with cervical cancer.

Variant analysis of human papillomavirus type 52 in Iranian women during 2018-2020: A case-control study

Background and Aims

Knowing the regional variants of distinct human papillomavirus (HPV) types is valuable as it can be beneficial for studying their epidemiology, pathogenicity, and evolution. For this reason, the sequence variations of the E6 gene of HPV 52 were investigated among women with normal cervical cytology and premalignant/malignant cervical samples.

Methods

Sixty-four HPV 52-positive samples were analyzed using semi-nested PCR and sequencing.

Results

Our findings showed that all samples belonged to lineage A (61%) or B (39%). Among samples that were infected with the A lineage, sublineages A1 and A2 were detected and sublineage A1 was dominant. No association was found between lineages and stage of disease (p > 0.05).

Conclusion

Our results revealed that the A lineage, sublineage A1, and B lineage were common in Iranian women. Nevertheless, more studies with larger sample sizes are required to estimate the pathogenicity risk of HPV 52 lineages in Iranian women with cervical cancer.

Global evaluation of lineage-specific human papillomavirus capsid antigenicity using antibodies elicited by natural infection

Human Papillomavirus (HPV) type variants have been classified into lineages and sublineages based upon their whole genome sequence. Here we have examined the specificity of antibodies generated following natural infection with lineage variants of oncogenic types (HPV16, 18, 31, 33, 45, 52 and 58) by testing serum samples assembled from existing archives from women residing in Africa, The Americas, Asia or Europe against representative lineage-specific pseudoviruses for each genotype. We have subjected the resulting neutralizing antibody data to antigenic clustering methods and created relational antigenic profiles for each genotype to inform the delineation of lineage-specific serotypes. For most genotypes, there was evidence of differential recognition of lineage-specific antigens and in some cases of a sufficient magnitude to suggest that some lineages should be considered antigenically distinct within their respective genotypes. These data provide compelling evidence for a degree of lineage specificity within the humoral immune response following natural infection with oncogenic HPV.

Whole-genome analysis and evolutionary characterization of cervical and oral human papillomavirus 16

High-throughput genome-wide sequencing has revealed high genomic variability of HPV16 in different geographic regions which is the most predominant genotype in human papillomavirus (HPV)-associated malignancies. Analysis of the HPV16 by whole-genome sequence (WGS) is an advanced method for the identification of mutations in the genome. There is limited information about HPV16 diversity in Pakistan, especially at the genomic level. Till now, WGS for HPV16 has not been previously reported in Pakistan. The current study has sequenced three HPV16 viral genomes, from two cervical and one oral cavity positive sample of women presented with general gynecological problems without any evidence of precancerous or cancerous lesions using an ion ampliseq customized panel. Sequencing analysis detected 38 variations, including single-nucleotide polymorphisms (SNPs) and two Indels, across three samples with the highest number of SNPs present in E1, E2, and L2, respectively. A total of 20 non-synonymous and 11 synonymous mutations with amino acid substitutions (T1421C, G1515A, T2223C, T1389C, G1483A, and T2191C) were identified. The phylogenetic analysis revealed the genomes of HPV16 are closely associated with those reported from Thailand and the United States. These are the first HPV16 WGS from Pakistan. However, more research is needed with a large sample size from diversified areas to assess the carcinogenic consequences and impact of HPV vaccinations.

Infectious diseases may have arrested the southward advance of microblades in Upper Palaeolithic East Asia

An unsolved archaeological puzzle of the East Asian Upper Palaeolithic is why the southward expansion of an innovative lithic technology represented by microblades stalled at the Qinling-Huaihe Line. It has been suggested that the southward migration of foragers with microblades stopped there, which is consistent with ancient DNA studies showing that populations to the north and south of this line had differentiated genetically by 19 000 years ago. Many infectious pathogens are believed to have been associated with hominins since the Palaeolithic, and zoonotic pathogens in particular are prevalent at lower latitudes, which may have produced a disease barrier. We propose a mathematical model to argue that mortality due to infectious diseases may have arrested the wave-of-advance of the technologically advantaged foragers from the north.

Retrospective Phylodynamic and Phylogeographic Analysis of the Human Papillomavirus 16 E6 Gene in the Mediterranean Region

Human papillomavirus 16 (HPV16) is considered to be strongly correlated with the development of cervical cancer. Among the 8 HPV16 genes, the E6 constitutes a remarkable marker to follow the evolutionary history and spatial phylodynamics of HPV16 in the Mediterranean basin. Thus, this work aims to decipher the major evolutionary events and crosstalks in the Mediterranean basin with a focus on Tunisian strains regarding the E6 oncogene. In this study, we first extracted the available and annotated Mediterranean strains of HPV16 E6 gene sequences (n = 155) from the NCBI nucleotide database. These sequences were aligned, edited, and used for the downstream phylogenetic analyses. Finally, a Bayesian Markov Chain Monte Carlo approach was applied to reconstruct the evolutionary history of HPV16 migration. Our results showed that the HPV circulating in Tunisia derived from a Croatian ancestor around the year 1987. This starting point spreads to most European countries to reach northern Africa through the Moroccan gateway in 2004.

Evolutionary analysis of JC polyomavirus in Misiones' population yields insight into the population dynamics of the early human dispersal in the Americas

BACKGROUND

JC polyomavirus (JCV) has an ethno-geographical distribution across human populations.

OBJECTIVE

Study the origins of the population of Misiones (Argentina) by using JCV as genetic marker.

METHODS

Viral detection and characterization was conducted by PCR amplification and evolutionary analysis of the intergenic region sequences.

RESULTS

22 out of 121 samples were positive for JCV, including 5 viral lineages: MY (n = 8), Eu-a (n = 7), B1-c (n = 4), B1-b (n = 2) and Af2 (n = 1). MY sequences clustered within a branch of Native American origin that diverged from its Asian counterpart about 21,914 years ago (HPD 95% interval 15,383-30,177), followed by a sustained demographic expansion around 5000 years ago.

CONCLUSIONS

JCV in Misiones reflects the multiethnic origin of the current population, with an important Amerindian contribution. Analysis of the MY viral lineage shows a pattern consistent with the arrival of early human migrations to the Americas and a population expansion by the pre-Columbian native societies.

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.

Human papillomavirus genomics: Understanding carcinogenicity

Human papillomavirus (HPV) causes virtually all cervical cancers and many cancers at other anatomical sites in both men and women. However, only 12 of 448 known HPV types are currently classified as carcinogens, and even the most carcinogenic type - HPV16 - only rarely leads to cancer. HPV is therefore necessary but insufficient for cervical cancer, with other contributing factors including host and viral genetics. Over the last decade, HPV whole genome sequencing has established that even fine-scale within-type HPV variation influences precancer/cancer risks, and that these risks vary by histology and host race/ethnicity. In this review, we place these findings in the context of the HPV life cycle and evolution at various levels of viral diversity: between-type, within-type, and within-host. We also discuss key concepts necessary for interpreting HPV genomic data, including features of the viral genome; events leading to carcinogenesis; the role of APOBEC3 in HPV infection and evolution; and methodologies that use deep (high-coverage) sequencing to characterize within-host variation, as opposed to relying on a single representative (consensus) sequence. Given the continued high burden of HPV-associated cancers, understanding HPV carcinogenicity remains important for better understanding, preventing, and treating cancers attributable to infection.

Synonymous nucleotide changes drive papillomavirus evolution

Papillomaviruses have been evolving alongside their hosts for at least 450 million years. This review will discuss some of the insights gained into the evolution of this diverse family of viruses. Papillomavirus evolution is constrained by pervasive purifying selection to maximize viral fitness. Yet these viruses need to adapt to changes in their environment, e.g., the host immune system. It has long been known that these viruses evolved a codon usage that doesn't match the infected host. Here we discuss how papillomavirus genomes evolve by acquiring synonymous changes that allow the virus to avoid detection by the host innate immune system without changing the encoded proteins and associated fitness loss. We discuss the implications of studying viral evolution, lifecycle, and cancer progression.

Disease-causing human viruses: novelty and legacy

About 270 viruses are known to infect humans. Some of these viruses have been known for centuries, whereas others have recently emerged. During their evolutionary history, humans have moved out of Africa to populate the world. In historical times, human migrations resulted in the displacement of large numbers of people. All these events determined the movement and dispersal of human-infecting viruses. Technological advances have resulted in the characterization of the genetic variability of human viruses, both in extant and in archaeological samples. Field studies investigated the diversity of viruses hosted by other animals. In turn, these advances provided insight into the evolutionary history of human viruses back in time and defined the key events through which they originated and spread.

Human Papillomavirus: Challenges and Opportunities for the Control of Cervical Cancer

Viruses are the most abundant and genetically diverse entities on the planet, infect all life forms and have evolved with their hosts. To date, 263 viral species have been identified that infect humans, of which only seven are considered type I oncogenic. Human papillomavirus (HPV) is the main virus associated with cancer and is responsible for practically all cases of cervical carcinoma. Screening tests for early detection have been available since the 1960s. Undoubtedly, the entailment between knowledge of HPV biology and the natural history of cervical cancer has contributed to the significant advances that have been made for its prevention since the 21^st century, with the development of prophylactic vaccines and improved screening strategies. Therefore, it is possible to eradicate invasive cervical cancer as a worldwide public health problem, as proposed by the WHO with the 90-70-90 initiative based on vaccination coverage, screening, and treatment, respectively. In addition, the emerging knowledge of viral biology generates opportunities that will contribute to strengthening prevention and treatment strategies in HPV-associated neoplasms.

Regional differences in human papillomavirus type 52 prevalence among Japanese women with cervical intraepithelial neoplasia†

Although geographical differences in the distribution of human papillomavirus genotypes have been observed worldwide, no studies have reported on national differences in the prevalence of human papillomavirus types in Japan. Here, we report a cross-sectional study to explore regional differences in the prevalence of human papillomavirus types among Japanese women with cervical intraepithelial neoplasia or invasive cervical cancer. Using human papillomavirus genotyping data from the nationwide prospective study on human papillomavirus vaccine effectiveness, we compared the frequency of detection of 15 high-risk and two low-risk human papillomavirus types in each disease category between the women who visited hospitals located in eastern Japan and those who visited hospitals located in western Japan. The risk of cervical intraepithelial neoplasia progression was assessed by calculating a prevalence ratio of each human papillomavirus type for cervical intraepithelial neoplasia grade 2/3 versus grade 1. Among the human papillomavirus types studied, human papillomavirus 52 was detected significantly more frequently in western hospitals than in eastern hospitals in cervical intraepithelial neoplasia grade 1 patients, but was less frequent in cervical intraepithelial neoplasia grade 2/3. The prevalence of particular human papillomavirus types was not significantly different between patients in hospitals in eastern Japan and those in hospitals in western Japan for invasive cervical cancer. In both eastern and western hospitals, a higher risk of cervical intraepithelial neoplasia progression was observed in patients infected with human papillomavirus 16, 31 or 58. In contrast, there was a significantly higher prevalence of human papillomavirus 52 infection in women with cervical intraepithelial neoplasia grade 2/3 than in those with cervical intraepithelial neoplasia grade 1 in eastern hospitals (prevalence ratio, 1.93; 95% confidence interval, 1.48-2.58), but not in western hospitals (prevalence ratio, 1.03; 95% confidence interval, 0.83-1.30). Regional differences of human papillomavirus 52 prevalence in cervical intraepithelial neoplasia lesions may exist and emphasize the importance of continuous monitoring of human papillomavirus type prevalence throughout the country in order to accurately assess the efficacy of human papillomavirus vaccines.

Ancient Evolutionary History of Human Papillomavirus Type 16, 18 and 58 Variants Prevalent Exclusively in Japan

Human papillomavirus (HPV) is a sexually transmitted virus with an approximately 8-kilo base DNA genome, which establishes long-term persistent infection in anogenital tissues. High levels of genetic variations, including viral genotypes and intra-type variants, have been described for HPV genomes, together with geographical differences in the distribution of genotypes and variants. Here, by employing a maximum likelihood method, we performed phylogenetic analyses of the complete genome sequences of HPV16, HPV18 and HPV58 available from GenBank (n = 627, 146 and 157, respectively). We found several characteristic clusters that exclusively contain HPV genomes from Japan: two for HPV16 (sublineages A4 and A5), one for HPV18 (sublineage A1) and two for HPV58 (sublineages A1 and A2). Bayesian phylogenetic analyses of concatenated viral gene sequences showed that divergence of the most recent common ancestor of these Japan-specific clades was estimated to have occurred ~98,000 years before present (YBP) for HPV16 A4, ~39,000 YBP for HPV16 A5, ~38,000 YBP for HPV18 A1, ~26,000 for HPV58 A1 and ~25,000 YBP for HPV58 A2. This estimated timeframe for the divergence of the Japan-specific clades suggests that the introduction of these HPV variants into the Japanese archipelago dates back to at least ~25,000 YBP and provides a scenario of virus co-migration with ancestral Japanese populations from continental Asia during the Upper Paleolithic period.

Molecular archeology of human viruses

The evolution of human-virus associations is usually reconstructed from contemporary patterns of genomic diversity. An intriguing, though still rarely implemented, alternative is to search for the genetic material of viruses in archeological and medical archive specimens to document evolution as it happened. In this chapter, we present lessons from ancient DNA research and incorporate insights from virology to explore the potential range of applications and likely limitations of archeovirological approaches. We also highlight the numerous questions archeovirology will hopefully allow us to tackle in the near future, and the main expected roadblocks to these avenues of research.

Phylogenomic Analysis of Human Papillomavirus Type 31 and Cervical Carcinogenesis: A Study of 2093 Viral Genomes

Human papillomavirus (HPV) type 31 (HPV31) is closely related to the most carcinogenic type, HPV16, but only accounts for 4% of cervical cancer cases worldwide. Viral genetic and epigenetic variations have been associated with carcinogenesis for other high-risk HPV types, but little is known about HPV31. We sequenced 2093 HPV31 viral whole genomes from two large studies, one from the U.S. and one international. In addition, we investigated CpG methylation in a subset of 175 samples. We evaluated the association of HPV31 lineages/sublineages, single nucleotide polymorphisms (SNPs) and viral methylation with cervical carcinogenesis. HPV31 A/B clade was >1.8-fold more associated with cervical intraepithelial neoplasia grade 3 and cancer (CIN3+) compared to the most common C lineage. Lineage/sublineage distribution varied by race/ethnicity and geographic region. A viral genome-wide association analysis identified SNPs within the A/B clade associated with CIN3+, including H23Y (C626T) (odds ratio = 1.60, confidence intervals = 1.17-2.19) located in the pRb CR2 binding-site within the E7 oncogene. Viral CpG methylation was higher in lineage B, compared to the other lineages, and was most elevated in CIN3+. In conclusion, these data support the increased oncogenicity of the A/B lineages and suggest variation of E7 as a contributing risk factor.

Human papillomaviruses: diversity, infection and host interactions

Human papillomaviruses (HPVs) are an ancient and highly successful group of viruses that have co-evolved with their host to replicate in specific anatomical niches of the stratified epithelia. They replicate persistently in dividing cells, hijack key host cellular processes to manipulate the cellular environment and escape immune detection, and produce virions in terminally differentiated cells that are shed from the host. Some HPVs cause benign, proliferative lesions on the skin and mucosa, and others are associated with the development of cancer. However, most HPVs cause infections that are asymptomatic and inapparent unless the immune system becomes compromised. To date, the genomes of almost 450 distinct HPV types have been isolated and sequenced. In this Review, I explore the diversity, evolution, infectious cycle, host interactions and disease association of HPVs.

Lineage analysis of human papillomavirus type 39 in cervical samples of Iranian women

Background

The data with regards to the regional variants of distinct HPV types is of great value. Accordance with this, this study aimed to investigate the sequence variations of E6 gene and long control region of HPV 39 among normal, premalignant and malignant cervical samples in order to characterize the frequent HPV 39 variants circulating in Tehran, Iran.

Methods

In total, 70 cervical samples (45 normal, 16 premalignant, and 9 malignant samples) infected with HPV 39 were analyzed by nested-PCR and sequencing.

Results

Our results revealed that all samples belonged to A lineage. Almost all sequences (98.6%) were classified in A1 sublineage and only one sample (1.4%) was A2 sub lineage.

Conclusions

Our findings showed that lineages A, sublineage A1, is dominant in Tehran, Iran. However, the small sample size was the most important limitations of this study. Further studies with larger sample size from different geographical regions of Iran are necessary to estimate the pathogenicity risk of HPV 39 variants in this population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01619-8.

Association of cervical carcinogenesis risk with HPV16 E6 and E7 variants in the Taizhou area, China

Background

Human papillomavirus (HPV) type 16 accounts for a larger share of cervical cancer and has been a major health problem worldwide for decades. The progression of initial infection to cervical cancer has been linked to viral sequence properties; however, the role of HPV16 variants in the risk of cervical carcinogenesis, especially with longitudinal follow-up, is not fully understood in China.

Methods

We aimed to investigate the genetic variability of HPV16 E6 and E7 oncogenes in isolates from cervical exfoliated cells. Between December 2012 and December 2014, a total of 310 single HPV16-positive samples were selected from women living in the Taizhou area, China. Sequences of all E6 and E7 oncogenes were analysed by PCR-sequencing assay. Detailed sequence comparison, genetic heterogeneity analyses and maximum-likelihood phylogenetic tree construction were performed with BioEdit Sequence Alignment Editor and MEGA X software. Data for cytology tests and histological diagnoses were obtained from our Taizhou Area Study with longitudinal follow-up for at least 5 years. The relationship between HPV16 variants and cervical carcinogenesis risk was analysed by the chi-square test or Fisher’s exact test.

Results

In this study, we obtained 64 distinct variation patterns with the accession GenBank numbers MT681266-MT681329. Phylogenetic analysis revealed that 98.3% of HPV16 variants belong to lineage A, in which the A4 (Asian) sublineage was dominant (64.8%), followed by A2 (12.1%), A1 (11.4%), and A3 (10.0%). The A4 (Asian) sublineage had a higher risk of CIN2+ than the A1–3 (European) sublineages (OR = 2.69, 95% CI = 1.04–6.97, P < 0.05). Furthermore, nucleotide variation in HPV16 E6 T178G is associated with the development of cervical cancer.

Conclusion

These data could provide novel insights into the role of HPV16 variants in cervical carcinogenesis risk in China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08531-y.

Lineage analysis of human papillomavirus type 18 based on E6 region in cervical samples of Iranian women

Distinct human papillomavirus (HPV) 18 variants are thought to differ in oncogenic potential and geographic distribution. As such, understanding the regional variants of HPV 18 would be of great importance for evolutionary, epidemiological, and biological analysis. In this regard, the sequence variations of E6 gene were investigated to characterize more common variants of HPV 18 in normal cells, premalignant, and malignant samples collected from the cervix. In total, 99 samples of HPV 18 were analyzed by polymerase chain reaction and sequencing. In overall, lineages A was identified in all study subjects, among which sublineage A4 was dominant although the difference observed was not statistically significant with regard to different stages of disease. Sublineage A4 comprised 90.9% of samples and the remaining were belonged to sublineages A1, A2, A3, and A5 at the frequency of 6.1%, 1%, 1%, and 1%, respectively. In conclusion, our findings clearly highlight the sublineage A4 of HPV 18 as the most dominant variant in Iran.

Association of HPV35 with cervical carcinogenesis among women of African ancestry: Evidence of viral-host interaction with implications for disease intervention

HPV35 has been found in only ∼2% of invasive cervical cancers (ICC) worldwide but up to 10% in Sub-Saharan Africa, warranting further investigation and consideration of impact on preventive strategies. We studied HPV35 and ethnicity, in relation to the known steps in cervical carcinogenesis, using multiple large epidemiologic studies in the U.S. and internationally. Combining five U.S. studies, we measured HPV35 positivity and, in Northern California, observed HPV35 type-specific population prevalence and estimated 5-year risk of developing precancer when HPV35-positive. HPV35 genetic variation was examined for differences in carcinogenicity in 1053 HPV35+ cervical specimens from a U.S. cohort and an international collection. African-American women had more HPV35 (12.1% vs 5.1%, P < .001) and more HPV35-associated precancers (7.4% vs 2.1%, P < .001) compared to other ethnicities. Precancer risks after HPV35 infection did not vary by ethnicity (global P = .52). The HPV35 A2 sublineage showed an increased association with precancer/cancer in African-Americans (OR = 5.6 vs A1, 95% CI = 1.3-24.8) and A2 was more prevalent among ICC in Africa than other world regions (41.9% vs 10.4%, P < .01). Our analyses support a strong link between HPV35 and cervical carcinogenesis in women of African ancestry. Current HPV vaccines cover the majority of cervical precancer/cancer across all ethnic groups; additional analyses are required to determine whether the addition of HPV35 to the already highly effective nine-valent HPV vaccine would provide better protection for women in Africa or of African ancestry.

You Will Never Walk Alone: Codispersal of JC Polyomavirus with Human Populations

JC polyomavirus (JCPyV) is one of the most prevalent human viruses. Findings based on the geographic distribution of viral subtypes suggested that JCPyV codiverged with human populations. This view was however challenged by data reporting a much more recent origin and expansion of JCPyV. We collected information on ∼1,100 worldwide strains and we show that their geographic distribution roughly corresponds to major human migratory routes. Bayesian phylogeographic analysis inferred a Subsaharan origin for JCPyV, although with low posterior probability. High confidence inference at internal nodes provided strong support for a long-standing association between the virus and human populations. In line with these data, pairwise FST values for JCPyV and human mtDNA sampled from the same areas showed a positive and significant correlation. Likewise, very strong relationships were found when node ages in the JCPyV phylogeny were correlated with human population genetic distances (nuclear-marker based FST). Reconciliation analysis detected a significant cophylogenetic signal for the human population and JCPyV trees. Notably, JCPyV also traced some relatively recent migration events such as the expansion of people from the Philippines/Taiwan area into Remote Oceania, the gene flow between North-Eastern Siberian and Ainus, and the Koryak contribution to Circum-Arctic Americans. Finally, different molecular dating approaches dated the origin of JCPyV in a time frame that precedes human out-of-Africa migration. Thus, JCPyV infected early human populations and accompanied our species during worldwide dispersal. JCPyV typing can provide reliable geographic information and the virus most likely adapted to the genetic background of human populations.

Molecular epidemiology of human papillomavirus 18 infections in Japanese Women

Among the oncogenic genotypes of human papillomavirus (HPV), HPV18 is the second most common type detected in cervical cancer worldwide and is primarily involved in the generation of cervical adenocarcinoma. Although HPV intra-type variants confer different risks of cervical carcinogenesis, there is little information on the genetic diversity of HPV18 compared to the most prevalent type, HPV16. In this study, we investigated the genetic variation of HPV18 in cervical specimens obtained from Japanese women with normal cervices or cervical cancers and precancers. Of the 101 HPV18-positive samples analyzed, viral whole genome amplification followed by next-generation sequencing led to the determination of viral complete genome sequences of 18 samples. Phylogenetic analysis of these HPV18 whole genome sequences identified a distinct variant cluster consisting of only Japanese samples (n = 7) belonging to sublineage A1. Viral genome sequences were also analyzed for the E6/E7 (n = 66) and E2 (n = 27) genes by Sanger sequencing. Phylogenetic analyses of these regions showed that the variant distribution among Japanese women was strongly biased toward sublineage A1 (72 of 87; 82.8%). No significant differences were observed in the prevalence of specific sublineages between cervical cancer/precancer cases and controls, and between squamous cell carcinoma and adenocarcinoma cases. These data contribute to our understanding of the genetic diversity of HPV18 in Japanese women.

Modelling the evolution of viral oncogenesis

Most human oncogenic viruses share several characteristics, such as being DNA viruses, having long (co)evolutionary histories with their hosts and causing either latent or chronic infections. They can reach high prevalences while causing relatively low case mortality, which makes them quite fit according to virulence evolution theory. After analysing the life histories of DNA oncoviruses, we use a mathematical modelling approach to investigate how the virus life cycle may generate selective pressures favouring or acting against oncogenesis at the within-host or at the between-host level. In particular, we focus on two oncoprotein activities, namely extending cell life expectancy and increasing cell proliferation rate. These have immediate benefits (increasing viral population size) but can be associated with fitness costs at the epidemiological level (increasing recovery rate or risk of cancer) thus creating evolutionary trade-offs. We interpret the results of our nested model in light of the biological features and identify future perspectives for modelling oncovirus dynamics and evolution. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

From Ancient to Emerging Infections: The Odyssey of Viruses in the Male Genital Tract

The male genital tract (MGT) is the target of a number of viral infections that can have deleterious consequences at the individual, offspring, and population levels. These consequences include infertility, cancers of male organs, transmission to the embryo/fetal development abnormalities, and sexual dissemination of major viral pathogens such as human immunodeficiency virus (HIV) and hepatitis B virus. Lately, two emerging viruses, Zika and Ebola, have additionally revealed that the human MGT can constitute a reservoir for viruses cleared from peripheral circulation by the immune system, leading to their sexual transmission by cured men. This represents a concern for future epidemics and further underlines the need for a better understanding of the interplay between viruses and the MGT. We review here how viruses, from ancient viruses that integrated the germline during evolution through old viruses (e.g., papillomaviruses originating from Neanderthals) and more modern sexually transmitted infections (e.g., simian zoonotic HIV) to emerging viruses (e.g., Ebola and Zika) take advantage of genital tract colonization for horizontal dissemination, viral persistence, vertical transmission, and endogenization. The MGT immune responses to viruses and the impact of these infections are discussed. We summarize the latest data regarding the sources of viruses in semen and the complex role of this body fluid in sexual transmission. Finally, we introduce key animal findings that are relevant for our understanding of viral infection and persistence in the human MGT and suggest future research directions.

Analysis of Human Papillomavirus (HPV) 16 Variants Associated with Cervical Infection in Italian Women

This study aims to evaluate HPV16 variants distribution in a population of Italian women living in two different regions (Lombardy and Sardinia) by sequence analyses of HPV16-positive cervical samples, in order to reconstruct the phylogenetic relationship among variants to identify the currently circulating lineages. Analyses were conducted starting from DNA isolated from 67 HPV16-positive cervical samples collected from two different Italian centres (31 from Lombardy and 36 from Sardinia) of women with normal and abnormal cervical cytology. The entire long control region (LCR) and 300 nt of the E6 gene was sequenced to identify intra-type variants. Sequence comparison and phylogenetic analysis were made using a distance-based neighbour joining method (NJ) and Kimura two-parameter model. Data obtained reported that Italian sequences mainly belonged to the European lineage, in particular sublineage A2. Only five sequences clustered in non-European branches: two in North American lineage (sublineage D1), two in African-1 (sublineage B1) and one in African-2. A new 27 nucleotide duplication in the central segment of the LCR region was found in a sequence obtained from a sample isolated in Sardinia. A predominance of European variants was detected, with some degree of variability among the studied HPV16 strains. This study contributes to the implementation of data regarding the molecular epidemiology of HPV16 variants.

Association of HLA-A, -B, DRB, and DQB Alleles with Persistent HPV-16 Infection in Women from Tamil Nadu, India

Women with persistent human papillomavirus (HPV) infections have a high risk of developing cervical cancer (CaCx). HPV-16 alone accounts for more than 60% of CaCx worldwide. Most of the HPV infections are transient and only a subset of women develop persistent HPV-16 infection. Many studies have shown associations of different human leukocyte antigen (HLA) alleles with HPV-mediated CaCx, but there are only a few studies globally that relate to persistent HPV-16 infection. Furthermore, such studies from India are sparse. Hence, we investigated the association of HLA-A, B, DRB, and DQB alleles with persistent HPV-16 infection and HPV-16-positive CaCx in south India (Tamil Nadu). HPV-16 persistent infection was observed in 7% of normal women. A total of 50 women with HPV-16-positive CaCx, 21 women with HPV-16 persistent infection, and 74 HPV-16-negative normal women were recruited for this study. Low-resolution typing of HLA-A, B, DRB, and DQB alleles was performed. HLA-B*44 and DRB1*07 showed a significant association with persistent HPV-16 infection (odds ratio, p-value = 26.3, 0.03 and 4.7, 0.01, respectively). HLA-B*27 and DRB1*12 were significantly associated with both HPV-16⁺ CaCx and persistent HPV-16 infection (23.8, 0.03; 52.9, 0.01; 9.8, 0.0009; and 13.8, 0.009; respectively). HLA-B*15 showed a negative association with HPV-16-positive CaCx (0.1, 0.01), whereas DRB1*04 exhibited protection to both HPV-16-positive CaCx and persistent HPV-16 infection (0.3, 0.0001 and 0.1, 0.0002, respectively). Thus, we show HLA allelic association with HPV-16 infection in Tamil Nadu. Larger studies on high-resolution HLA typing coupled with HPV-16 genome diversity will offer further insights into host/pathogen genome coevolution.

Disease transmission and introgression can explain the long-lasting contact zone of modern humans and Neanderthals

Neanderthals and modern humans both occupied the Levant for tens of thousands of years prior to the spread of modern humans into the rest of Eurasia and their replacement of the Neanderthals. That the inter-species boundary remained geographically localized for so long is a puzzle, particularly in light of the rapidity of its subsequent movement. Here, we propose that infectious-disease dynamics can explain the localization and persistence of the inter-species boundary. We further propose, and support with dynamical-systems models, that introgression-based transmission of alleles related to the immune system would have gradually diminished this barrier to pervasive inter-species interaction, leading to the eventual release of the inter-species boundary from its geographic localization. Asymmetries between the species in the characteristics of their associated ‘pathogen packages’ could have generated feedback that allowed modern humans to overcome disease burden earlier than Neanderthals, giving them an advantage in their subsequent spread into Eurasia.

Modern humans and Neanderthals coexisted in the Levant for tens of thousands of years before modern humans spread and replaced Neanderthals. Here, Greenbaum et al. develop a model showing that transmission of disease and genes can explain the maintenance and then collapse of this contact zone.

Mutational landscape and intra-host diversity of human papillomavirus type 16 long control region and E6 variants in cervical samples

Human papillomavirus genotype 16 (HPV16) is the most frequent high-risk HPV (HR-HPV) identified in cervical precursor lesions and cervical cancer (CC) worldwide. The oncogenic potential of HPV16 is partly dependent on the lineage involved in the infection and the presence of clinically relevant mutations. In this report, we present the distribution of HR-HPV and the mutational profile and intra-host variability of HPV16 lineages, based on analysis of the long control region (LCR) and the E6 gene in samples with normal cytology (n = 39), squamous intraepithelial lesions (n = 25), and CC (n = 39). HR-HPV genotyping was performed using multiplex real-time PCR. HPV16 lineage assignments and mutation frequencies were determined by conventional PCR and Sanger DNA sequencing, and intra-patient viral populations were analyzed using next-generation sequencing (NGS). The most frequent HR-HPV type was HPV16, followed by HPV31 and HPV18. The frequency of HPV16 sublineages was A1/A2 > D2 > D3 and B1. Moreover, the most frequent mutations, both in samples from this study and in the available sequences from Mexican isolates in the GenBank database were LCR-G7518A, which is involved in carcinogenesis, and E6-T350G (producing L83V), associated with persistence of infection. Otherwise, deep sequencing revealed high conservation of viral lineages and mutations, independently of the stages studied. In conclusion, the high frequency and stability of these molecular markers, as well as the circulating viral lineages, could be related to the incidence of CC associated with HPV16. Hence, they deserve a broader analysis to determine the risk of specific populations for progression of the disease.

Cytomegalovirus distribution and evolution in hominines

Herpesviruses are thought to have evolved in very close association with their hosts. This is notably the case for cytomegaloviruses (CMVs; genus Cytomegalovirus) infecting primates, which exhibit a strong signal of co-divergence with their hosts. Some herpesviruses are however known to have crossed species barriers. Based on a limited sampling of CMV diversity in the hominine (African great ape and human) lineage, we hypothesized that chimpanzees and gorillas might have mutually exchanged CMVs in the past. Here, we performed a comprehensive molecular screening of all 9 African great ape species/subspecies, using 675 fecal samples collected from wild animals. We identified CMVs in eight species/subspecies, notably generating the first CMV sequences from bonobos. We used this extended dataset to test competing hypotheses with various degrees of co-divergence/number of host switches while simultaneously estimating the dates of these events in a Bayesian framework. The model best supported by the data involved the transmission of a gorilla CMV to the panine (chimpanzee and bonobo) lineage and the transmission of a panine CMV to the gorilla lineage prior to the divergence of chimpanzees and bonobos, more than 800,000 years ago. Panine CMVs then co-diverged with their hosts. These results add to a growing body of evidence suggesting that viruses with a double-stranded DNA genome (including other herpesviruses, adenoviruses, and papillomaviruses) often jumped between hominine lineages over the last few million years.

The HPV16 Genome Is Stable in Women Who Progress to In Situ or Invasive Cervical Cancer: A Prospective Population-Based Study

The human papillomavirus (HPV) rate of evolution is essential for cancer-preventive strategies targeting HPV. We analyzed variability over time in a prospective, population-based nested case-control study of in situ (CIS) and invasive squamous cervical cancer (SCC). Among 757,690 women who participated in cervical screening in Sweden during 1969 to 2002, we identified 94 women who had HPV16 persistence in two serial cervical screening samples (median 24 months apart, range 0.5-178 months) and later were diagnosed with CIS (n = 59), SCC (n = 32), or remained healthy (n = 3). Whole-HPV16-genome sequencing and comparison of sequences in the serial samples revealed that all women had the same HPV16 lineage, particularly lineage A, in both serial smears. Fifty-six percent of women had an identical 7,906 base pair HPV16 sequence in both samples, and no woman had more than 15 nucleotide substitutions. The median substitution rate was 0 substitutions/site/year (95% confidence interval, 0-0.00008), with no variation between quartiles of follow-up. We concluded that in most women with HPV16 persistence preceding disease, the nucleotide substitution rate was not measurable within up to 15-years follow-up. This slow rate of evolution has important implications for both HPV-based screening and HPV vaccination. SIGNIFICANCE: These findings show there is no genomic variation over time in HPV16 infections progressing to cervical cancer, which could influence risk stratification of women when screening for cervical cancer and inform HPV vaccination strategies.

Human papillomavirus 16 sub-lineage dispersal and cervical cancer risk worldwide: Whole viral genome sequences from 7116 HPV16-positive women

BACKGROUND

Human papillomavirus (HPV)16 can be separated into genetic sub-lineages (A1-4, B1-4, C1-4, D1-4) which may have differential cervical cancer risk.

METHODS

A next-generation sequencing assay was used to whole-genome sequence 7116 HPV16-positive cervical samples from well-characterised international epidemiological studies, including 2076 controls, 1878 squamous cell carcinoma (SCC) and 186 adenocarcinoma/adenosquamous cell carcinoma (ADC), and to assign HPV16 sub-lineage. Logistic regression was used to estimate region-stratified country-adjusted odds ratios (OR) and 95%CI.

RESULTS

A1 was the most globally widespread sub-lineage, with others showing stronger regional specificity (A3 and A4 for East Asia, B1-4 and C1-4 for Africa, D2 for the Americas, B4, C4 and D4 for North Africa). Increased cancer risks versus A1 were seen for A3, A4 and D (sub)lineages in regions where they were common: A3 in East Asia (OR=2.2, 95%CI:1.0-4.7); A4 in East Asia (6.6, 3.1-14.1) and North America (3.8, 1.7-8.3); and D in North (6.2, 4.1-9.3) and South/Central America (2.2, 0.8-5.7), where D lineages were also more frequent in ADC than SCC (3.2, 1.5-6.5; 12.1, 5.7-25.6, respectively).

CONCLUSIONS

HPV16 genetic variation can strongly influence cervical cancer risk. However, burden of cervical cancer attributable to different sub-lineages worldwide is largely driven by historical HPV16 sub-lineage dispersal.

Impact of naturally occurring variation in the human papillomavirus 52 capsid proteins on recognition by type-specific neutralising antibodies

We investigated the impact of naturally occurring variation within the major (L1) and minor (L2) capsid proteins on the antigenicity of human papillomavirus (HPV) type 52 (HPV52). L1L2 pseudoviruses (PsVs) representing HPV52 lineage and sublineage variants A1, A2, B1, B2, C and D were created and tested against serum from naturally infected individuals, preclinical antisera raised against HPV52 A1 and D virus-like particles (VLPs) and neutralising monoclonal antibodies (MAbs) raised against HPV52 A1 VLP. HPV52 lineage D PsV displayed a median 3.1 (inter-quartile range 2.0-5.6) fold lower sensitivity to antibodies elicited following natural infection with, where data were available, HPV52 lineage A. HPV52 lineage variation had a greater impact on neutralisation sensitivity to pre-clinical antisera and MAbs. Chimeric HPV52 A1 and D PsV were created which identified variant residues in the FG (Q281K) and HI (K354T, S357D) loops as being primarily responsible for the reported differential sensitivities. Homology models of the HPV52 L1 pentamer were generated which permitted mapping these residues to a small cluster on the outer rim of the surface exposed pentameric L1 protein. These data contribute to our understanding of HPV L1 variant antigenicity and may have implications for seroprevalence or vaccine immunity studies based upon HPV52 antigens.

Phylogeny and Classification of Human Papillomavirus (HPV)16 and HPV18 Variants Based on E6 and L1 genes in Tunisian Women with Cervical Lesions

Background:

High-risk human papillomavirus (HPV) types are the main etiological factors for cervical cancer. HPV16 and HPV18 are generally the most common forms associated with development of high-grade cervical lesions. This study was undertaken to identify intratypic variants of HPV16 and HPV18 among women with cervical lesions in Tunisia.

Materials and Methods:

DNA was extracted from cervical samples collected from 49 women. using a PureLinkTM Genomic DNA mini Kit (Invitrogen). E6 and L1 open reading frames (ORF) were amplified by PCR and viral DNA amplicons were subjected to automated sequencing using Big Dye Terminators technology (Applied Biosystems). The obtained sequences were analyzed using an appropriate software program to allow phylogenetic trees to be generated.

Results:

HPV16 and HPV18 were detected in 15 and 5 cases, respectively. HPV16 E6 sequences clustered with the European German lineage (A2) whereas one isolate diverged differently in the L1 region and clustered with the African sub-lineage (B1). HPV 18 E6 sequences clustered with the European sub-lineage (A1) but L1 sequences clustered as a new clade which diverged from A1-A5.

Conclusions:

Our results suggest that the distribution of HPV16 and HPV18 sequences in women with cervical lesions in Tunisia is mainly related to European epidemiological conditions and point to the presence of recombinant HPV forms.

A meta-analysis of association of Human Leukocyte Antigens A, B, C, DR and DQ with Human Papillomavirus 16 infection

Human Papillomavirus (HPV) induced cervical cancer (CaCx) is a major health problem in women from both developing and developed regions of the world. This virus accounts for >95% of the CaCx cases with a preponderance of HPV type -16 (65%). Paradoxically HPV-16 is prevalent even in the cervix of healthier women and anti HPV-16 T-cell response is considered critical for the viral clearance. Studies on HLA association with HPV-16 infection and cervical cancer have yielded varied HLA associations in different epidemiological settings. To validate these associations, we performed a meta-analysis of HLA-A, B, C, DR and DQ association with HPV-16 infection. Of the 1409 studies retrieved, 26 qualified for meta-analysis based on stringent inclusion and exclusion criteria. HLA-B*47, B*57, DRB1*10, DRB1*15 and DQB1*0303 were significantly associated with HPV-16 infection (OR = 3.4, 1.8, 1.5, 1.1 and 1.5 respectively). HLA-B*49, B*39, A28 (serotype), C*04 and DRB1*13 were negatively associated with HPV-16 (OR = 0.5, 0.6, 0.7, 0.7, and 0.7 respectively). Certain HLA alleles such as B*07, DRB1*15, DRB1*11 and DRB1*07 showed weakly positive associations. A comprehensive analysis coupling HPV-16 antigenic diversity and the HLA variation in various global populations shall provide further insights into the immunogenetic predisposition to HPV-16 and shall help identify host-parasite co-evolution.

Niche adaptation and viral transmission of human papillomaviruses from archaic hominins to modern humans

Recent discoveries on the origins of modern humans from multiple archaic hominin populations and the diversity of human papillomaviruses (HPVs) suggest a complex scenario of virus-host evolution. To evaluate the origin of HPV pathogenesis, we estimated the phylogeny, timing, and dispersal of HPV16 variants using a Bayesian Markov Chain Monte Carlo framework. To increase precision, we identified and characterized non-human primate papillomaviruses from New and Old World monkeys to set molecular clock models. We demonstrate specific host niche adaptation of primate papillomaviruses with subsequent coevolution with their primate hosts for at least 40 million years. Analyses of 212 HPV16 complete genomes and 3582 partial sequences estimated ancient divergence of HPV16 variants (between A and BCD lineages) from their most recent common ancestors around half a million years ago, roughly coinciding with the timing of the split between archaic Neanderthals and modern Homo sapiens, and nearly three times longer than divergence times of modern Homo sapiens. HPV16 A lineage variants were significantly underrepresented in present African populations, whereas the A sublineages were highly prevalent in European (A1-3) and Asian (A4) populations, indicative of viral sexual transmission from Neanderthals to modern non-African humans through multiple interbreeding events in the past 80 thousand years. Remarkably, the human leukocyte antigen B*07:02 and C*07:02 alleles associated with increased risk in cervix cancer represent introgressed regions from Neanderthals in present-day Eurasians. The archaic hominin-host-switch model was also supported by other HPV variants. Niche adaptation and virus-host codivergence appear to influence the pathogenesis of papillomaviruses.

Intragenus (Homo) variation in a chemokine receptor gene (CCR5)

Humans have a comparatively higher rate of more polymorphisms in regulatory regions of the primate CCR5 gene, an immune system gene with both general and specific functions. This has been interpreted as allowing flexibility and diversity of gene expression in response to varying disease loads. A broad expression repertoire is useful to humans-the only globally distributed primate-due to our unique adaptive pattern that increased pathogen exposure and disease loads (e.g., sedentism, subsistence practices). The main objective of the study was to determine if the previously observed human pattern of increased variation extended to other members of our genus, Homo. The data for this study are mined from the published genomes of extinct hominins (four Neandertals and two Denisovans), an ancient human (Ust'-Ishim), and modern humans (1000 Genomes). An average of 15 polymorphisms per individual were found in human populations (with a total of 262 polymorphisms). There were 94 polymorphisms identified across extinct Homo (an average of 13 per individual) with 41 previously observed in modern humans and 53 novel polymorphisms (32 in Denisova and 21 in Neandertal). Neither the frequency nor distribution of polymorphisms across gene regions exhibit significant differences within the genus Homo. Thus, humans are not unique with regards to the increased frequency of regulatory polymorphisms and the evolution of variation patterns across CCR5 gene appears to have originated within the genus. A broader evolutionary perspective on regulatory flexibility may be that it provided an advantage during the transition to confrontational foraging (and later hunting) that altered human-environment interaction as well as during migration to Eurasia and encounters with novel pathogens.

[A biological anthropology of the disappearance of the Neandertal Man: recent data]

What could have been the causes of the disappearance of Neanderthals? We will try here to make a synthesis between one of the fundamental questions of biological anthropology relating to human evolution (hypotheses on the causes of the extinction of Neanderthals) and evolutionary bio-medical concepts, some of which have recently been reformulated thanks to the progress of paleogenomics (ancestral inheritance of the current human immune system, paleo-microbiology, host-pathogen relationship…).

The Role of aDNA in Understanding the Coevolutionary Patterns of Human Sexually Transmitted Infections

Analysis of pathogen genome data sequenced from clinical and historical samples has made it possible to perform phylogenetic analyses of sexually transmitted infections on a global scale, and to estimate the diversity, distribution, and coevolutionary host relationships of these pathogens, providing insights into pathogen emergence and disease prevention. Deep-sequenced pathogen genomes from clinical studies and ancient samples yield estimates of within-host and between-host evolutionary rates and provide data on changes in pathogen genomic stability and evolutionary responses. Here we examine three groups of pathogens transmitted mainly through sexual contact between modern humans to provide insight into ancient human behavior and history with their pathogens. Exploring ancient pathogen genomic divergence and the ancient viral-host parallel evolutionary histories will help us to reconstruct the origin of present-day geographical distribution and diversity of clinical pathogen infections, and will hopefully allow us to foresee possible environmentally induced pathogen evolutionary responses. Lastly, we emphasize that ancient pathogen DNA research should be combined with modern clinical pathogen data, and be equitable and provide advantages for all researchers worldwide, e.g., through shared data.

Paleo-oncology in the Dakhleh Oasis, Egypt: Case studies and a paleoepidemiological perspective

This article describes six cases of cancer from the Dakhleh Oasis, Egypt. A mummy had a confirmed 'primary' diagnosis of adenocarcinoma of the rectum. The remaining diagnoses were based on the distribution and types of skeletal lesions in conjunction with age, sex, and/or the molecular phylogeny of the Human Papillomavirus (HPV). HPV is a confirmed cause of cancer of the uterine cervix (UC) and testes (TC) and it evolved in Africa long before Homo sapiens emerged. Today these cancers are common in young adult females and males, a fact which was pivotal in respectively including them in the differential diagnosis of UC and TC. The remaining diagnoses were acute lymphocytic leukemia in a 3-5 year old child and an older female with metastatic carcinoma. Due to problems of determining specific diagnoses and their prevalence in 'paleo' populations, we opted for a lifetime cancer risk statistic (LTCR). The LTCR in ancient Dakhleh was ∼5/1000 (6/1087). In modern Western societies the LTCR cancer approaches 50% (500/1000). Thus the LTCR in today's western societies is 100 times greater than in ancient Dakhleh. These cases demonstrate that oncogenes and their environmental cofactors were present in antiquity, but were significantly less pervasive than today.

Ancient cancers and infection-induced oncogenesis

Cancers have been reported in bone and soft tissue of ancient agricultural populations. Fossilized bones from prehistoric periods provide evidence of tumors but only one example of cancer. Difficulties in diagnosing the causes of lesions in mummified tissue and fossilized bone, and in interpreting the prevalence of cancers from remains, draw attention to the need for complementary approaches to assess the occurrence of cancer in ancient populations. This paper integrates current knowledge about pathogen induction of cancer with phylogenetic analyses of oncogenic pathogens, and concludes that pathogen-induced cancers were probably generally present in ancient historic and prehistoric human populations. Consideration of cancers in extant human populations and wildlife lends credence to this conclusion, with the caveat that the presence of cancers may depend on population-specific exposures to oncogenic parasites and carcinogens.

Intra-host sequence variability in human papillomavirus

Human papillomaviruses (HPVs) co-evolve slowly with the human host and each HPV genotype displays epithelial tropisms. We assessed the evolution of intra HPV genotype variants within samples, and their association to anogenital site, cervical cytology and HIV status. Variability in the L1 gene of 35 HPV genotypes was characterized phylogenetically using maximum likelihood, and portrayed by phenotype. Up to a thousand unique variants were identified within individual samples. In-depth analyses of the most prevalent genotypes, HPV16, HPV18 and HPV52, revealed that the high diversity was dominated by a few abundant variants. This suggests high intra-host mutation rates. Clades of HPV16, HPV18 and HPV52 were associated to anatomical site and HIV co-infection. Particularly, we observed that one HPV16 clade was specific to vaginal cells and one HPV52 clade was specific to anal cells. One major HPV52 clade, present in several samples, was strongly associated with cervical neoplasia. Overall, our data suggest that tissue tropism and HIV immunosuppression are strong shapers of HPV evolution.

Classification and evolution of human papillomavirus genome variants: Alpha-5 (HPV26, 51, 69, 82), Alpha-6 (HPV30, 53, 56, 66), Alpha-11 (HPV34, 73), Alpha-13 (HPV54) and Alpha-3 (HPV61)

HPV variants from the same type can be classified into lineages and sublineages based on the complete genome differences and the phylogenetic topologies. We examined nucleotide variations of twelve HPV types within the species Alpha-5 (HPV26, 51, 69, 82), Alpha-6 (HPV30, 53, 56, 66), Alpha-11 (HPV34, 73), Alpha-13 (HPV54) and Alpha-3 (HPV61) by analyzing 1432 partial sequences and 181 complete genomes from multiple geographic populations. The inter-lineage and inter-sublineage mean differences of HPV variants ranged between 0.9-7.3% and 0.3-0.9%, respectively. The heterogeneity and phylogenies of HPV isolates indicate an independent evolutionary history for each type. The noncoding regions were the most variable regions whereas the capsid proteins were relatively conserved. Certain variant lineages and/or sublineages were geographically-associated. These data provide the basis to further classify HPV variants and should foster future studies on the evolution of HPV genomes and the associations of HPV variants with cancer risk.

The Intersection of HPV Epidemiology, Genomics and Mechanistic Studies of HPV-Mediated Carcinogenesis

Of the ~60 human papillomavirus (HPV) genotypes that infect the cervicovaginal epithelium, only 12–13 “high-risk” types are well-established as causing cervical cancer, with HPV16 accounting for over half of all cases worldwide. While HPV16 is the most important carcinogenic type, variants of HPV16 can differ in their carcinogenicity by 10-fold or more in epidemiologic studies. Strong genotype-phenotype associations embedded in the small 8-kb HPV16 genome motivate molecular studies to understand the underlying molecular mechanisms. Understanding the mechanisms of HPV genomic findings is complicated by the linkage of HPV genome variants. A panel of experts in various disciplines gathered on 21 November 2016 to discuss the interdisciplinary science of HPV oncogenesis. Here, we summarize the discussion of the complexity of the viral–host interaction and highlight important next steps for selected applied basic laboratory studies guided by epidemiological genomic findings.