Two Human papillomavirus 11 complete genomes recovered from inverted sinonasal papillomas in humans

We communicate here two complete Human papillomavirus 11 (HPV11) genomes recovered from one transitional and from one squamous inverted sinonasal papilloma, a rare proliferative disease in humans. Both genomes belong to the HPV11_A2 sublineage.

Evolution and diversity of nucleotide and dinucleotide composition in poxviruses

Poxviruses (family Poxviridae) have long dsDNA genomes and infect a wide range of hosts, including insects, birds, reptiles and mammals. These viruses have substantial incidence, prevalence and disease burden in humans and in other animals. Nucleotide and dinucleotide composition, mostly CpG and TpA, have been largely studied in viral genomes because of their evolutionary and functional implications. We analysed here the nucleotide and dinucleotide composition, as well as codon usage bias, of a set of representative poxvirus genomes, with a very diverse host spectrum. After correcting for overall nucleotide composition, entomopoxviruses displayed low overall GC content, no enrichment in TpA and large variation in CpG enrichment, while chordopoxviruses showed large variation in nucleotide composition, no obvious depletion in CpG and a weak trend for TpA depletion in GC-rich genomes. Overall, intergenome variation in dinucleotide composition in poxviruses is largely accounted for by variation in overall genomic GC levels. Nonetheless, using vaccinia virus as a model, we found that genes expressed at the earliest times in infection are more CpG-depleted than genes expressed at later stages. This observation has parallels in betahepesviruses (also large dsDNA viruses) and suggests an antiviral role for the innate immune system (e.g. via the zinc-finger antiviral protein ZAP) in the early phases of poxvirus infection. We also analysed codon usage bias in poxviruses and we observed that it is mostly determined by genomic GC content, and that stratification after host taxonomy does not contribute to explaining codon usage bias diversity. By analysis of within-species diversity, we show that genomic GC content is the result of mutational biases. Poxvirus genomes that encode a DNA ligase are significantly AT-richer than those that do not, suggesting that DNA repair systems shape mutation biases. Our data shed light on the evolution of poxviruses and inform strategies for their genetic manipulation for therapeutic purposes.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.

Sweat Gland Tumors Arising on Acral Sites: A Molecular Survey

Recurrent oncogenic drivers have been identified in a variety of sweat gland tumors. Recently, integration of human papillomavirus type 42 (HPV42) has been reported in digital papillary adenocarcinoma (DPA). The main objectives of the present study were (i) to provide an overview of the prevalence of previously identified oncogenic drivers in acral sweat gland tumors and (ii) to genetically characterize tumors in which no recurrent genetic alteration has been identified yet. Cases of acral sweat gland tumors were identified from the database of the French network CARADERM. After histologic review, the presence of previously identified genetic alterations was investigated in the entire cohort (n=79) using a combination of immunohistochemistry and targeted DNA and RNA sequencing. Tumor entities with no recurrent genetic alterations were submitted to whole-transcriptome sequencing. CRTC1::MAML2 fusion was identified in cases of hidradenoma and hidradenocarcinoma (n=9/12 and n=9/12). A p.V600E mutation of BRAF was observed in all cases of tubular adenoma (n=4). YAP1:MAML2 and YAP1::NUTM1 fusions were observed in poroid tumors (n=15/25). ETV6::NTRK3 and TRPS1::PLAG1 fusion transcripts were identified in secretory carcinoma (n=1/1) and cutaneous mixed tumors (n=3/4), respectively. The HPV42 genome was detected in most cases of DPA (n=10/11) and in 1 adnexal adenocarcinoma not otherwise specified. Finally, whole-transcriptome analysis revealed BRD3::NUTM1 or NSD3::NUTM1 fusions in 2 cases of NUT adnexal carcinoma and NCOA4::RET and CCDC6::RET fusion transcripts in 2 cystadenoma/hidrocystoma-like tumors. Our study confirms distinctive cytogenetic abnormalities in a wide number of acral adnexal neoplasms and supports the use of molecular analysis as a valuable aid in the diagnosis of these rare and often difficult to diagnose group of neoplasms.

Digital Papillary Adenocarcinoma in Nonacral Skin: Clinicopathologic and Genetic Characterization of 5 Cases

Digital papillary adenocarcinoma (DPA) is a rare sweat gland neoplasm that has exceptionally been reported outside acral locations. Recently, human papillomavirus 42 was identified as the main oncogenic driver of DPA. Herein, we report 5 tumors arising in extra-acral locations predominantly in the female anogenital skin. Four patients were female and 1 patient was male. The mean age at the diagnosis time was 65 years (range: 55 to 82 y). Tumors were located on the vulva (n=3), perianal area (n=1), and forearm (n=1). Histologically, all tumors were lobular and mainly solid and composed of sheets of cells with rare focal papillae and frequent glandular structures in a "back-to-back" pattern and lined by atypical basophilic cells. Immunohistochemistry showed diffuse positivity for SOX10. Epithelial membrane antigen and carcinoembryonic antigen highlighted the luminal cells and staining for p63 and p40 revealed a consistent and continuous myoepithelial component around glandular structures. Follow-up was available in 3 cases (mean duration: 12 mo [range: 8 to 16 mo]). One patient developed local recurrence and 1 experienced regional lymph node metastases. HPV Capture Next-generation sequencing revealed the presence of the HPV42 genome in all samples. Viral reads distributions were compatible in the 5 cases with an episomal nature of the viral genome, with a recurrent deletion in the E1 and/or E2 open reading frames. In conclusion, this study demonstrates that digital DPA may rarely present in nonacral locations mainly in the female anogenital area, usually with a more solid pattern as compared with those cases presenting on the digits and it is also associated with HPV42.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

Does exposure to different types of menstrual protections affect the vaginal environment?

The vaginal ecosystem is a key component of women's health. It also represents an ideal system for ecologists to investigate the consequence of perturbations on species diversity and emerging properties between organisational levels. Here, we study how exposure to different types of menstrual protection products is linked to microbial, immunological, demographic, and behavioural measurements in a cohort of young adult women who reported using more often tampons (n = 107) or menstrual cups (n = 31). We first found that cup users were older and smoked less than tampon users. When analysing health indicators, we detected potential associations between cups use reporting and fungal genital infection. A multivariate analysis confirmed that, in our cohort, reporting using cups over tampons was associated with the higher odds ratio to report a fungal genital infection diagnosis by a medical doctor within the last 3 months. We did not detect significant differences between groups in terms of their bacterial vaginal microbiota composition and found marginal differences in the level of expression of 20 cytokines. However, a multivariate analysis of these biological data identified some level of clustering based on the type of menstrual production preferred (cups or tampons). These results suggest that exposure to different types of menstrual products could influence menstrual health. Larger studies and studies with a more powered setting are needed to assess the robustness of these associations and identify causal mechanisms.

Multiresistant Enterobacteriaceae in yellow-legged gull chicks in their first weeks of life

Wild animal species living in anthropogenic areas are commonly carriers of antimicrobial‐resistant bacteria (AMRB), but their role in the epidemiology of these bacteria is unclear. Several studies on AMRB in wildlife have been cross‐sectional in design and sampled individual animals at only one point in time. To further understand the role of wildlife in maintaining and potentially transmitting these bacteria to humans and livestock, longitudinal studies are needed in which samples are collected from individual animals over multiple time periods. In Europe, free‐ranging yellow‐legged gulls (Larus michahellis) commonly live in industrialized areas, forage in landfills, and have been found to carry AMRB in their feces. Using bacterial metagenomics and antimicrobial resistance characterization, we investigated the spatial and temporal patterns of AMRB in a nesting colony of yellow‐legged gulls from an industrialized area in southern France. We collected 54 cloacal swabs from 31 yellow‐legged gull chicks in 20 nests on three dates in 2016. We found that AMRB in chicks increased over time and was not spatially structured within the gull colony. This study highlights the complex occurrence of AMRB in a free‐ranging wildlife species and contributes to our understanding of the public health risks and implications associated with ARMB‐carrying gulls living in anthropogenic areas.

Detecting anal human papillomavirus infection in men who have sex with men living with HIV: implications of assay variability

BACKGROUND

Incidence of anal cancer (AC) caused by persistent human papillomavirus (HPV) infection has risen in the last years in men who have sex with men (MSM) living with HIV. There is consensus that this population should be screened for anal precancerous lesions, but the role of HPV DNA testing in AC screening programmes is still under debate.

OBJECTIVES

This study employed two molecular test to detect anal HPV DNA and compared assay performance and prognostic value for the diagnosis of histology proven high-grade intraepithelial anal lesions.

METHODS

MSM living with HIV attended their regular check-up visits consisting of detection of anal HPV infection, anal cytology, digital anorectal examination and high resolution anoscopy. HPV DNA was detected using Hybrid Capture 2 High-Risk test (HC2, total assay) and LINEAR ARRAY HPV Genotyping Test (LA, type-specific assay) RESULTS: Among 274 participant, prevalence of HPV DNA was 48.5% by HC2 and 89.4% by LA. HPV16 (30.6%) and HPV6 (19.6%) were the most common genotypes identified. Prevalence of multiple HPV infections was 56.2%. Agreement between HPV DNA assays was 75.2% (κ=0.51; 95% CI 0.42 to 0.60). Total HPV detection demonstrated high sensitivity (90%; 95% CI 68.3 to 98.8) and moderate specificity (58.4%; 95% CI 50.2 to 66.3), while type-specific HPV16/18 genotyping provided an increase in specificity and showed the highest area under the curve (0.81; 95% CI 0.74 to 0.89) and Youden's index (0.63).

CONCLUSIONS

Both methodologies identified a high prevalence of anal HPV infection and multiple HPV infections in MSM living with HIV, showing a moderate overall agreement between them. Either total HPV detection or type-specific HPV16/18 detection together with a threshold ≥atypical squamous cells of undetermined significance for abnormal cytology showed an acceptable diagnostic accuracy.

Concomitant and productive genital infections by HSV-2 and HPV in two young women: A case report

Human papillomaviruses (HPVs), the most oncogenic virus known to humans, are often associated with Herpes Simplex Virus-2 (HSV-2) infections. The involvement of the latter in cervical cancer is controversial but its long-term infections might modulate the mucosal microenvironment in a way that favors carcinogenesis. We know little about coinfections between HSV-2 and HPVs, and studying the immunological and microbiological dynamics in the early stages of these infections may help identify or rule out potential interactions. We report two cases of concomitant productive, although asymptomatic, HSV-2 and HPV infections in young women (aged 20 and 25). The women were followed up for approximately a year, with clinical visits every two months and weekly self-samples. We performed quantitative analyses of their HSV-2 and HPV viral loads, immunological responses (IgG and IgM antibodies and local cytokines expression profiles), vaginal microbiota composition, as well as demographic and behavior data. We detect interactions between virus loads, immune response, and the vaginal microbiota, which improve our understanding of HSV-2 and HPVs’ coinfections and calls for further investigation with larger cohorts.

SARS-CoV-2 viral RNA levels are not 'viral load'

Ct values are commonly used as proxies of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) 'viral load'. Since coronaviruses are positive single-stranded RNA [(+)ssRNA] viruses, current reverse transcription (RT)-qPCR target amplification does not distinguish replicative from transcriptional RNA. Although analyses of Ct values remain informative, equating them with viral load may lead to flawed conclusions as it is presently unknown whether (and to what extent) variation in Ct reflects variation in viral load or in gene expression.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

Human DNA decays faster with time than viral dsDNA: an analysis on HPV16 using pathology archive samples spanning 85 years

Background

Quality of the nucleic acids extracted from Formalin Fixed Paraffin Embedded (FFPE) samples largely depends on pre-analytic, fixation and storage conditions. We assessed the differential sensitivity of viral and human double stranded DNA (dsDNA) to degradation with storage time.

Methods

We randomly selected forty-four HPV16-positive invasive cervical cancer (ICC) FFPE samples collected between 1930 and 1935 and between 2000 and 2004. We evaluated through qPCR the amplification within the same sample of two targets of the HPV16 L1 gene (69 bp, 134 bp) compared with two targets of the human tubulin-β gene (65 bp, 149 bp).

Results

Both viral and human, short and long targets were amplified from all samples stored for 15 years. In samples archived for 85 years, we observed a significant decrease in the ability to amplify longer targets and this difference was larger in human than in viral DNA: longer fragments were nine times (CI 95% 2.6–35.2) less likely to be recovered from human DNA compared with 1.6 times (CI 95% 1.1–2.2) for viral DNA.

Conclusions

We conclude that human and viral DNA show a differential decay kinetics in FFPE samples. The faster degradation of human DNA should be considered when assessing viral DNA prevalence in long stored samples, as HPV DNA detection remains a key biomarker of viral-associated transformation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01529-9.

Genomic and phylogenetic characterization of ChPV2, a novel goat PV closely related to the Xi-PV1 species infecting bovines

BACKGROUND

Papillomaviruses (PVs) infecting artiodactyls are very diverse, and only second in number to PVs infecting primates. PVs associated to lesions in economically important ruminant species have been isolated from cattle and sheep.

METHODS

Potential PV DNA from teat lesions of a Damascus goat was isolated, cloned and sequenced. The PV genome was analyzed using bioinformatics approaches to detect open reading frames and to predict potential features of encoded proteins as well as putative regulatory elements. Sequence comparison and phylogenetic analyses using the concatenated E1E2L2L1 nucleotide and amino acid alignments was used to reveal the relationship of the new PV to the known PV diversity and its closest relevants.

RESULTS

We isolated and characterized the full-genome of novel Capra hircus papillomavirus. We identified the E6, E7, E1, E2, L2, L1 open reading frames with protein coding potential and putative active elements in the ChPV2 proteins and putative regulatory genome elements. Sequence similarities of L1 and phylogenetic analyses using concatenated E1E2L2L1 nucleotide and amino acid alignments suggest the classification as a new PV type designated ChPV2 with a phylogenetic position within the XiPV genus, basal to the XiPV1 species. ChPV2 is not closely related to ChPV1, the other known goat PV isolated from healthy skin, although both of them belong confidently into a clade composed of PVs infecting cervids and bovids. Interestingly, ChPV2 contains an E6 open reading frame whereas all closely related PVs do not CONCLUSION: ChPV2 is a novel goat PV closely related to the Xi-PV1 species infecting bovines. Phylogenetic relationships and genome architecture of ChPV2 and closely related PV types suggest at least two independent E6 losses within the XiPV clade.

Loop 1 of APOBEC3C Regulates its Antiviral Activity against HIV-1

APOBEC3 deaminases (A3s) provide mammals with an anti-retroviral barrier by catalyzing dC-to-dU deamination on viral ssDNA. Within primates, A3s have undergone a complex evolution via gene duplications, fusions, arms race, and selection. Human APOBEC3C (hA3C) efficiently restricts the replication of viral infectivity factor (vif)-deficient Simian immunodeficiency virus (SIVΔvif), but for unknown reasons, it inhibits HIV-1Δvif only weakly. In catarrhines (Old World monkeys and apes), the A3C loop 1 displays the conserved amino acid pair WE, while the corresponding consensus sequence in A3F and A3D is the largely divergent pair RK, which is also the inferred ancestral sequence for the last common ancestor of A3C and of the C-terminal domains of A3D and A3F in primates. Here, we report that modifying the WE residues in hA3C loop 1 to RK leads to stronger interactions with substrate ssDNA, facilitating catalytic function, which results in a drastic increase in both deamination activity and in the ability to restrict HIV-1 and LINE-1 replication. Conversely, the modification hA3F_WE resulted only in a marginal decrease in HIV-1Δvif inhibition. We propose that the two series of ancestral gene duplications that generated A3C, A3D-CTD and A3F-CTD allowed neo/subfunctionalization: A3F-CTD maintained the ancestral RK residues in loop 1, while diversifying selection resulted in the RK → WE modification in Old World anthropoids' A3C, possibly allowing for novel substrate specificity and function.

Papillomaviruses infecting cetaceans exhibit signs of genome adaptation following a recombination event

Papillomaviruses (PVs) have evolved through a complex evolutionary scenario where virus-host co-evolution alone is not enough to explain the phenotypic and genotypic PV diversity observed today. Other evolutionary processes, such as host switch and recombination, also appear to play an important role in PV evolution. In this study, we have examined the genomic impact of a recombination event between distantly related PVs infecting Cetartiodactyla (even-toed ungulates and cetaceans). Our phylogenetic analyses suggest that one single recombination was responsible for the generation of extant 'chimeric' PV genomes infecting cetaceans. By correlating the phylogenetic relationships to the genomic content, we observed important differences between the recombinant and non-recombinant cetartiodactyle PV genomes. Notably, recombinant PVs contain a unique set of conserved motifs in the upstream regulatory region (URR). We interpret these regulatory changes as an adaptive response to drastic changes in the PV genome. In terms of codon usage preferences (CUPrefs), we did not detect any particular differences between orthologous open reading frames in recombinant and non-recombinant PVs. Instead, our results are in line with previous observations suggesting that CUPrefs in PVs are rather linked to gene expression patterns as well as to gene function. We show that the non-coding URR of PVs infecting cetaceans, the central regulatory element in these viruses, exhibits signs of adaptation following a recombination event. Our results suggest that also in PVs, the evolution of gene regulation can play an important role in speciation and adaptation to novel environments.

Detection of Alpha, Beta, Gamma, and Unclassified Human Papillomaviruses in Cervical Cancer Samples From Mexican Women

Background: Cervical cancer (CC) is associated to high-risk human papillomavirus (HPV) infections, for this reason it is crucial to have sensitive and accurate HPV diagnostic tests. To date, most research is focused on HPVs within the Alphapapillomavirus (α-PVs) genus and little attention has been paid to cervical infections with other HPV genotypes, like those of the Betapapillomavirus (β-PVs) and Gammapapillomavirus (γ-PVs) genera. The aim of this study was to determine the HPV genotypes from different genera in women with CC using Next-Generation Sequencing (NGS). Methods: The study comprised 48 HPV positive CC samples evaluated with the Linear Array HPV Genotyping test and individually sequenced by 454 NGS using PGMY09/11 and FAP primers. To determine the HPV genotypes present in each sample, the obtained sequences were compared with all HPV L1 gene reference sequences from the Papillomavirus Episteme database (PaVE). Moreover, 50 HPV positive low-grade cervical lesion samples individually genotyped with NGS were also included to determine the genotypes present preferentially in CC patients. Results: Among the 48 CC samples, 68.75% consisted of multiple HPV infections, 51 different genotypes were detected, of which 7 are still unclassified, 28 belong to α-PVs (6, 11, 16, 18, 26, 30, 33, 35, 39, 42, 43, 44, 45, 51, 52, 53, 54, 59, 62, 66, 68, 69, 70, 71, 74, 81, 102, 114), 10 to β-PVs (5, 12, 21, 37, 38b, 47, 80, 107, 118, 122), and 6 to γ-PVs (101, 103, 123, 135, 147, 214). Among them, HPV16 was the most prevalent genotype (54.2%), followed by HPV18 (16.7%), HPV38b (14.6%), and HPVs 52/62/80 (8.3%). Some genotypes were exclusively found in CC when compared with Cervical Intraepithelial Neoplasia grade 1 (CIN1) samples, such as HPVs 5, 18, 38b, 107, 122, FA39, FA116, mSK_120, and mSK_136. Conclusions: This work demonstrates the great diversity of HPV genotypes detected by combining PGMY and FAP primers with NGS in cervical swabs. The relatively high attribution of β- and γ- PVs in CC samples suggest their possible role as carcinogenic cofactors, but deeper studies need to be performed to determine if they have transforming properties and the significance of HPV-coinfections.

Towards a multi-level and a multi-disciplinary approach to DNA oncovirus virulence

One out of 10 cancers is estimated to arise from infections by a handful of oncogenic viruses. These infectious cancers constitute an opportunity for primary prevention through immunization against the viral infection, for early screening through molecular detection of the infectious agent, and potentially for specific treatments, by targeting the virus as a marker of cancer cells. Accomplishing these objectives will require a detailed understanding of the natural history of infections, the mechanisms by which the viruses contribute to disease, the mutual adaptation of viruses and hosts, and the possible viral evolution in the absence and in the presence of the public health interventions conceived to target them. This issue showcases the current developments in experimental tissue-like and animal systems, mathematical models and evolutionary approaches to understand DNA oncoviruses. Our global aim is to provide proximate explanations to the present-day interface and interactions between virus and host, as well as ultimate explanations about the adaptive value of these interactions and about the evolutionary pathways that have led to the current malignant phenotype of oncoviral infections. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

COUSIN (COdon Usage Similarity INdex): A Normalized Measure of Codon Usage Preferences

Codon Usage Preferences (CUPrefs) describe the unequal usage of synonymous codons at the gene, chromosome, or genome levels. Numerous indices have been developed to evaluate CUPrefs, either in absolute terms or with respect to a reference. We introduce the normalized index COUSIN (for COdon Usage Similarity INdex), that compares the CUPrefs of a query against those of a reference and normalizes the output over a Null Hypothesis of random codon usage. The added value of COUSIN is to be easily interpreted, both quantitatively and qualitatively. An eponymous software written in Python3 is available for local or online use (http://cousin.ird.fr). This software allows for an easy and complete analysis of CUPrefs via COUSIN, includes seven other indices, and provides additional features such as statistical analyses, clustering, and CUPrefs optimization for gene expression. We illustrate the flexibility of COUSIN and highlight its advantages by analyzing the complete coding sequences of eight divergent genomes. Strikingly, COUSIN captures a bimodal distribution in the CUPrefs of human and chicken genes hitherto unreported with such precision. COUSIN opens new perspectives to uncover CUPrefs specificities in genomes in a practical, informative, and user-friendly way.

Genome Plasticity in Papillomaviruses and De Novo Emergence of E5 Oncogenes

The clinical presentations of papillomavirus (PV) infections come in many different flavors. While most PVs are part of a healthy skin microbiota and are not associated to physical lesions, other PVs cause benign lesions, and only a handful of PVs are associated to malignant transformations linked to the specific activities of the E5, E6, and E7 oncogenes. The functions and origin of E5 remain to be elucidated. These E5 open reading frames (ORFs) are present in the genomes of a few polyphyletic PV lineages, located between the early and the late viral gene cassettes. We have computationally assessed whether these E5 ORFs have a common origin and whether they display the properties of a genuine gene. Our results suggest that during the evolution of Papillomaviridae, at least four events lead to the presence of a long noncoding DNA stretch between the E2 and the L2 genes. In three of these events, the novel regions evolved coding capacity, becoming the extant E5 ORFs. We then focused on the evolution of the E5 genes in AlphaPVs infecting primates. The sharp match between the type of E5 protein encoded in AlphaPVs and the infection phenotype (cutaneous warts, genital warts, or anogenital cancers) supports the role of E5 in the differential oncogenic potential of these PVs. In our analyses, the best-supported scenario is that the five types of extant E5 proteins within the AlphaPV genomes may not have a common ancestor. However, the chemical similarities between E5s regarding amino acid composition prevent us from confidently rejecting the model of a common origin. Our evolutionary interpretation is that an originally noncoding region entered the genome of the ancestral AlphaPVs. This genetic novelty allowed to explore novel transcription potential, triggering an adaptive radiation that yielded three main viral lineages encoding for different E5 proteins, displaying distinct infection phenotypes. Overall, our results provide an evolutionary scenario for the de novo emergence of viral genes and illustrate the impact of such genotypic novelty in the phenotypic diversity of the viral infections.

Evolutionary Changes after Translational Challenges Imposed by Horizontal Gene Transfer

Genes acquired by horizontal gene transfer (HGT) may provide the recipient organism with potentially new functions, but proper expression level and integration of the transferred genes in the novel environment are not granted. Notably, transferred genes can differ from the receiving genome in codon usage preferences, leading to impaired translation and reduced functionality. Here, we characterize the genomic and proteomic changes undergone during experimental evolution of Escherichia coli after HGT of three synonymous versions, presenting very different codon usage preference, of an antibiotic resistance gene. The experimental evolution was conducted with and without the corresponding antibiotic and the mutational patterns and proteomic profiles after 1,000 generations largely depend on the experimental growth conditions (e.g., mutations in antibiotic off-target genes), and on the synonymous gene version transferred (e.g., mutations in genes responsive to translational stress). The transfer of an exogenous gene extensively modifies the whole proteome, and these proteomic changes are different for the different version of the transferred gene. Additionally, we identified conspicuous changes in global regulators and in intermediate metabolism, confirmed the evolutionary ratchet generated by mutations in DNA repair genes and highlighted the plasticity of bacterial genomes accumulating large and occasionally transient duplications. Our results support a central role of HGT in fuelling evolution as a powerful mechanism promoting rapid, often dramatic genotypic and phenotypic changes. The profound reshaping of the pre-existing geno/phenotype allows the recipient bacteria to explore new ways of functioning, far beyond the mere acquisition of a novel function.

Natural history, dynamics, and ecology of human papillomaviruses in genital infections of young women: protocol of the PAPCLEAR cohort study

INTRODUCTION

Human papillomaviruses (HPVs) are responsible for one-third of all cancers caused by infections. Most HPV studies focus on chronic infections and cancers, and we know little about the early stages of the infection. Our main objective is to better understand the course and natural history of cervical HPV infections in healthy, unvaccinated and vaccinated, young women, by characterising the dynamics of various infection-related populations (virus, epithelial cells, vaginal microbiota and immune effectors). Another objective is to analyse HPV diversity within hosts, and in the study population, in relation to co-factors (lifestyle characteristics, vaccination status, vaginal microbiota, human genetics).

METHODS AND ANALYSIS

The PAPCLEAR study is a single center longitudinal study following 150 women, aged 18-25 years, for up to 2 years. Visits occur every 2 or 4 months (depending on HPV status) during which several variables are measured, such as behaviours (via questionnaires), vaginal pH, HPV presence and viral load (via qPCR), local concentrations of cytokines (via MesoScale Discovery technology) and immune cells (via flow cytometry). Additional analyses are outsourced, such as titration of circulating anti-HPV antibodies, vaginal microbiota sequencing (16S and ITS1 loci) and human genotyping. To increase the statistical power of the epidemiological arm of the study, an additional 150 women are screened cross-sectionally. Finally, to maximise the resolution of the time series, participants are asked to perform weekly self-samples at home. Statistical analyses will involve classical tools in epidemiology, genomics and virus kinetics, and will be performed or coordinated by the Centre National de la Recherche Scientifique (CNRS) in Montpellier.

ETHICS AND DISSEMINATION

This study has been approved by the Comité de Protection des Personnes Sud Méditerranée I (reference number 2016-A00712-49); by the Comité Consultatif sur le Traitement de l'Information en matière de Recherche dans le domaine de la Santé (reference number 16.504); by the Commission Nationale Informatique et Libertés (reference number MMS/ABD/AR1612278, decision number DR-2016-488) and by the Agence Nationale de Sécurité du Médicament et des Produits de Santé (reference 20160072000007). Results will be published in preprint servers, peer-reviewed journals and disseminated through conferences.

TRIAL REGISTRATION NUMBER

NCT02946346; Pre-results.

Origin and evolution of papillomavirus (onco)genes and genomes

Papillomaviruses (PVs) are ancient viruses infecting vertebrates, from fishes to mammals. Although the genomes of PVs are small and show conserved synteny, PVs display large genotypic diversity and ample variation in the phenotypic presentation of the infection. Most PV genomes contain two small early genes E6 and E7. In a bunch of closely related human papillomaviruses (HPVs), the E6 and E7 proteins provide the viruses with oncogenic potential. The recent discoveries of PVs without E6 and E7 in different fish species place a new root on the PV tree, and suggest that ancestral PVs consisted of the minimal PV backbone E1-E2-L2-L1. Bayesian phylogenetic analyses date the most recent common ancestor of the PV backbone to 424 million years ago (Ma). Common ancestry tests on extant E6 and E7 genes indicate that they share a common ancestor dating back to at least 184 Ma. In AlphaPVs infecting Old World monkeys and apes, the appearance of the E5 oncogene 53-58 Ma concurred with (i) a significant increase in substitution rate, (ii) a basal radiation and (iii) key gain of functions in E6 and E7. This series of events was instrumental to construct the extant phenotype of oncogenic HPVs. Our results assemble the current knowledge on PV diversity and present an ancient evolutionary timeline punctuated by evolutionary innovations in the history of this successful viral family. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Human papillomavirus DNA detected in fingertip, oral and bathroom samples from unvaccinated adolescent girls in Tanzania

OBJECTIVE

Human papillomavirus (HPV) DNA has been detected in vaginal samples from adolescent girls who report no previous sex and, in high-income settings, from fingertips, raising the possibility of non-sexual transmission. No such studies originate from East Africa which bears among the highest cervical cancer incidence and HPV prevalence worldwide. HPV-related oral cancer incidence is increasing, but oral HPV prevalence data from East Africa are limited. We aimed to describe the HPV DNA prevalence in genital and non-genital sites and in the bathroom of unvaccinated adolescent girls, and examine genotype concordance between sites.

METHODS

We nested a cross-sectional study of HPV in genital and extragenital sites within a cohort study of vaginal HPV acquisition. Unvaccinated girls age 16-18 years in Tanzania, who reported ever having had sex, were consented, enrolled and tested for the presence of HPV DNA in vaginal samples collected using self-administered swabs, oral samples collected using an oral rinse, and on fingertips and bathroom surfaces collected using a cytobrush.

RESULTS

Overall, 65 girls were enrolled and 23 (35%, 95% CI 23% to 47%) had detectable vaginal HPV. Adequate (β-globin positive) samples were collected from 36 girls' fingertips and HPV was detected in 7 (19%, 95% CI 6% to 33%). 63 girls provided adequate oral samples, 4 (6%, 95% CI 0% to 13%) of which had HPV DNA detected. In bathroom samples from 58 girls, 4 (7%, 95% CI 0% to 14%) had detectable HPV DNA. Of the 11 girls with extragenital HPV, six had the same genotype in >1 site.

CONCLUSION

We found a high prevalence of HPV in non-genital sites in adolescent girls and in their bathrooms, in this region with a high cervical cancer incidence. Concordance of genotypes between sites supports the possibility of autoinoculation.

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.

The Use of HPV16-E5, EGFR, and pEGFR as Prognostic Biomarkers for Oropharyngeal Cancer Patients

Background: Anti-epidermal-growth-factor-receptor (EGFR) therapies in combination with radiotherapy are being studied on de-escalation clinical trials for HPV-related oropharyngeal cancer (OPC) patients. The HPV16-E5 oncoprotein increases recycling of activated EGFR to the cell surface, enhancing factor signal transduction. Our aim was to evaluate viral HPV16-E5 oncogene expression as well as EGFR and phosphorylated-EGFR (pEGFR), protein levels as biomarkers for clinical outcome in a retrospective cohort of OPC patients.

Methods: Formalin-fixed-paraffin-embedded OPCs were collected from 1990 to 2013. OPC samples containing HPV-DNA were subject to viral E6^*I mRNA detection and p16^INK4a immunohistochemistry (IHC). HPV16-positive cases were evaluated for HPV16-E5 (RT-PCR) and EGFR/pEGFR (IHC). A stratified and matched random sample of HPV-negative samples was used as control and evaluated for EGFR/pEGFR. Overall survival (OS) and disease free survival (DFS) estimates were assessed for locally advanced OPC patients (stage III, IVa,b 7th edition).

Results: Among 788 OPC patient samples, 53 were double positive for HPV16-DNA/p16^INK4a. HPV16-E5 expression was found in 41 of 53 samples (77.4%). EGFR expression was observed in 37.7 vs 70.8% of HPV16-positive vs HPV-negative samples, respectively; (adjusted OR = 0.15) 5% CI = 0.04–0.56]). Expression of pEGFR followed an inverse pattern with 39.6 and 24.9% detection in HPV16-positive and HPV-negative samples; (adjusted OR = 1.58 [95% CI = 0.48–5.17]). Within HPV16-positive cases, no association between HPV16-E5/EGFR nor pEGFR was observed. With a median follow-up of 39.36 months (min = 0.03 – max = 272.07), the combination of HPV status and EGFR or pEGFR expression were predictors of better OS (p < 0.001, for both) and DFS (p < 0.001 for EGFR and p = 0.003 for pEGFR).

Conclusions: HPV16-E5 is highly expressed on HPV16-positive OPCs. Interestingly, HPV16-positive cases expressed significantly more pEGFR while HPV-negative cases expressed more EGFR. The combinations of HPV status and EGFR or pEGFR may be useful biomarkers for evaluating prognosis outcome in OPC patients.

Two Novel, Distantly Related Papillomaviruses Isolated from Healthy Skin of the Timor Deer (Rusa timorensis)

We report the complete genome sequences of Rusa timorensis papillomavirus 1 (RtimPV1) and Rusa timorensis papillomavirus 2 (RtimPV2), isolated from hair follicles of asymptomatic skin from the same Timor deer specimen. RtimPV1 and RtimPV2 are evolutionarily only distantly related. RtimPV1 lacks a canonical E2-binding site, and RtimPV2 does not carry an E6 gene.

Genome Sequences of Two Novel Papillomaviruses Isolated from Healthy Skin of Pudu puda and Cervus elaphus Deer

We report the complete genome sequences of Pudu puda papillomavirus1 (PpudPV1) and Cervus elaphus papillomavirus2 (CelaPV2), isolated from healthy skin hair follicles of a Southern pudu and a red deer, respectively. PpudPV1 is basal to the DyokappaPVs, whereas CelaPV2 is basal to the XiPVs (Beta-XiPV crown group).

Stably expressed APOBEC3H forms a barrier for cross-species transmission of simian immunodeficiency virus of chimpanzee to humans

APOBEC3s (A3s) are potent restriction factors of human immunodeficiency virus type 1/simian immunodeficiency viruses (HIV-1/SIV), and can repress cross-species transmissions of lentiviruses. HIV-1 originated from a zoonotic infection of SIV of chimpanzee (SIVcpz) to humans. However, the impact of human A3s on the replication of SIVcpz remains unclear. By using novel SIVcpz reporter viruses, we identified that human APOBEC3B (A3B) and APOBEC3H (A3H) haplotype II strongly reduced the infectivity of SIVcpz, because both of them are resistant to SIVcpz Vifs. We further demonstrated that human A3H inhibited SIVcpz by deaminase dependent as well independent mechanisms. In addition, other stably expressed human A3H haplotypes and splice variants showed strong antiviral activity against SIVcpz. Moreover, most SIV and HIV lineage Vif proteins could degrade chimpanzee A3H, but no Vifs from SIVcpz and SIV of gorilla (SIVgor) lineages antagonized human A3H haplotype II. Expression of human A3H hapII in human T cells efficiently blocked the spreading replication of SIVcpz. The spreading replication of SIVcpz was also restricted by stable A3H in human PBMCs. Thus, we speculate that stably expressed human A3H protects humans against the cross-species transmission of SIVcpz and that SIVcpz spillover to humans may have started in individuals that harbor haplotypes of unstable A3H proteins.

Cellular cytidine deaminases of the APOBEC3 (A3) family are potent restriction factors that are able to inhibit retroviruses, this A3 restriction is counteracted by lentivirus Vif proteins. Human APOBEC3H (A3H) represents the most evolutionarily divergent A3 gene; it includes seven haplotypes and several splice variants. The polymorphism of human A3H has relevance for HIV-1 infection and AIDS progression. HIV-1 originated from cross-transmission of SIVcpz to humans. However, little is known about how human A3s affect the replication or transmission of SIVcpz. In this study, we comprehensively analyzed the anti-SIVcpz activity of chimpanzee and human A3s. Human A3H haplotype II was identified as strong inhibitor against SIVcpz regardless of Vif. In addition, other stably expressed human A3H haplotypes and splice variants showed strong antiviral activity against SIVcpz. Moreover, based on the recent Great Ape Genome Project, we found that the polymorphism of chimpanzee A3H is lower compared with the diversity of human A3H. And chimpanzee A3H haplotypes identified in this study showed similar anti-SIVcpz activity and Vif sensitivity. Our results provide a model that stably expressed human A3H protects humans against the cross-species transmission of SIVcpz and that SIVcpz spillover to humans may have started in individuals that harbor haplotypes of unstable A3H proteins.

Why Human Papillomavirus Acute Infections Matter

Most infections by human papillomaviruses (HPVs) are `acute', that is non-persistent. Yet, for HPVs, as for many other oncoviruses, there is a striking gap between our detailed understanding of chronic infections and our limited data on the early stages of infection. Here we argue that studying HPV acute infections is necessary and timely. Focusing on early interactions will help explain why certain infections are cleared while others become chronic or latent. From a molecular perspective, descriptions of immune effectors and pro-inflammatory pathways during the initial stages of infections have the potential to lead to novel treatments or to improved handling algorithms. From a dynamical perspective, adopting concepts from spatial ecology, such as meta-populations or meta-communities, can help explain why HPV acute infections sometimes last for years. Furthermore, cervical cancer screening and vaccines impose novel iatrogenic pressures on HPVs, implying that anticipating any viral evolutionary response remains essential. Finally, hints at the associations between HPV acute infections and fertility deserve further investigation given their high, worldwide prevalence. Overall, understanding asymptomatic and benign infections may be instrumental in reducing HPV virulence.

Evolutionary analysis of Old World arenaviruses reveals a major adaptive contribution of the viral polymerase

The Old World (OW) arenavirus complex includes several species of rodent-borne viruses, some of which (i.e., Lassa virus, LASV and Lymphocytic choriomeningitis virus, LCMV) cause human diseases. Most LCMV and LASV infections are caused by rodent-to-human transmissions. Thus, viral evolution is largely determined by events that occur in the wildlife reservoirs. We used a set of human- and rodent-derived viral sequences to investigate the evolutionary history underlying OW arenavirus speciation, as well as the more recent selective events that accompanied LASV spread in West Africa. We show that the viral RNA polymerase (L protein) was a major positive selection target in OW arenaviruses and during LASV out-of-Nigeria migration. No evidence of selection was observed for the glycoprotein, whereas positive selection acted on the nucleoprotein (NP) during LCMV speciation. Positively selected sites in L and NP are surrounded by highly conserved residues, and the bulk of the viral genome evolves under purifying selection. Several positively selected sites are likely to modulate viral replication/transcription. In both L and NP, structural features (solvent exposed surface area) are important determinants of site-wise evolutionary rate variation. By incorporating several rodent-derived sequences, we also performed an analysis of OW arenavirus codon adaptation to the human host. Results do not support a previously hypothesized role of codon adaptation in disease severity for non-Nigerian strains. In conclusion, L and NP represent the major selection targets and possible determinants of disease presentation; these results suggest that field surveys and experimental studies should primarily focus on these proteins.

Plasmid and clonal interference during post horizontal gene transfer evolution

Plasmids are nucleic acid molecules that can drive their own replication in a living cell. They can be transmitted horizontally and can thrive in the host cell to high-copy numbers. Plasmid replication and gene expression consume cellular resources and cells carrying plasmids incur fitness costs. But many plasmids carry genes that can be beneficial under certain conditions, allowing the cell to endure in the presence of antibiotics, toxins, competitors or parasites. Horizontal transfer of plasmid-encoded genes can thus instantaneously confer differential adaptation to local or transient selection conditions. This conflict between cellular fitness and plasmid spread sets the scene for multilevel selection processes. We have engineered a system to study the short-term evolutionary impact of different synonymous versions of a plasmid-encoded antibiotic resistance gene. Applying experimental evolution under different selection conditions and deep sequencing allowed us to show rapid local adaptation to the presence of antibiotic and to the specific version of the resistance gene transferred. We describe the presence of clonal interference at two different levels: at the within-cell level, because a single cell can carry several plasmids, and at the between-cell level, because a bacterial population may contain several clones carrying different plasmids and displaying different fitness in the presence/absence of antibiotic. Understanding the within-cell and between-cell dynamics of plasmids after horizontal gene transfer is essential to unravel the dense network of mobile elements underlying the worldwide threat to public health of antibiotic resistance.

Differential HPV16 variant distribution in squamous cell carcinoma, adenocarcinoma and adenosquamous cell carcinoma

Human Papillomavirus 16 (HPV16) causes 70% of invasive cervical cancers (ICC) worldwide. Interaction between HPV16 genetic diversity, host genetics and target tissue largely determine the chances to trigger carcinogenesis. We have analyzed the differential prevalence of viral variants in 233 HPV16-monoinfected squamous (SCC), glandular (ADC) and mixed (ADSC) ICCs from four continents, assessing the contribution of geographical origin and cancer histology. We have further quantified the contribution of viral variants and cancer histology to differences in age at tumor diagnosis. The model fitted to the data explained 97% of the total variance: the largest explanatory factors were differential abundance among HPV16 variants (78%) and their interaction with cancer histology (9.2%) and geography (10.1%). HPV16_A1-3 variants were more prevalent in SCC while HPV16_D variants were increased in glandular ICCs. We confirm further a non-random geographical structure of the viral variants distribution. ADCs were diagnosed at younger ages than SCCs, independently of the viral variant triggering carcinogenesis. HPV16 variants are differentially associated with histological ICCs types, and ADCs are systematically diagnosed in younger women. Our results have implications for the implementation of cervical cancer screening algorithms, to ensure proper early detection of elusive ADCs.

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.

Searching beyond the usual papillomavirus suspects in squamous carcinomas of the vulva, penis and head and neck

Human Papillomaviruses (HPVs) are involved in the etiology of anogenital and head and neck cancers. The HPV DNA prevalence greatly differs by anatomical site. Indeed, the high rates of viral DNA prevalence in anal and cervical carcinomas contrast with the lower fraction of cancer cases attributable to HPVs in other anatomical sites, chiefly the vulva, the penis and head and neck. Here we analyzed 2635 Formalin Fixed Paraffin Embedded surgical samples that had previously tested negative for the presence of HPVs DNA using the SPF10/DEIA procedure, in order to identify the presence of other PVs not explicitly targeted by standard molecular epidemiologic approaches. All samples were reanalyzed using five broad-PV PCR primer sets (CP1/2, FAP6064/FAP64, SKF/SKR, MY9/MY11, MFI/MFII) targeting the main PV main clades. In head and neck carcinoma samples (n=1141), we recovered DNA from two BetaHPVs, namely HPV20 and HPV21, and from three cutaneous AlphaPVs, namely HPV2, HPV57 and HPV61. In vulvar squamous cell carcinoma samples (n=902), we found one of the samples containing DNA of one cutaneous HPV, namely HPV2, and 29 samples contained DNA from essentially mucosal HPVs. In penile squamous cell carcinoma samples (n=592), we retrieved the DNA of HPV16 in 16 samples. Our results show first that the SPF10/DEIA is very sensitive, as we recovered only 2.1% (55/2635) false negative results; second, that although the DNA of cutaneous HPVs can be detected in cancer samples, their relative contribution remains anyway minor (0.23%; 6/2635) and may be neglected for screening and vaccination purposes; and third, their contribution to malignancy is not necessarily warranted and needs to be elucidated.

Determinants of FIV and HIV Vif sensitivity of feline APOBEC3 restriction factors

Background

Feline immunodeficiency virus (FIV) is a global pathogen of Felidae species and a model system for Human immunodeficiency virus (HIV)-induced AIDS. In felids such as the domestic cat (Felis catus), APOBEC3 (A3) genes encode for single-domain A3Z2s, A3Z3 and double-domain A3Z2Z3 anti-viral cytidine deaminases. The feline A3Z2Z3 is expressed following read-through transcription and alternative splicing, introducing a previously untranslated exon in frame, encoding a domain insertion called linker. Only A3Z3 and A3Z2Z3 inhibit Vif-deficient FIV. Feline A3s also are restriction factors for HIV and Simian immunodeficiency viruses (SIV). Surprisingly, HIV-2/SIV Vifs can counteract feline A3Z2Z3.

Results

To identify residues in feline A3s that Vifs need for interaction and degradation, chimeric human–feline A3s were tested. Here we describe the molecular direct interaction of feline A3s with Vif proteins from cat FIV and present the first structural A3 model locating these interaction regions. In the Z3 domain we have identified residues involved in binding of FIV Vif, and their mutation blocked Vif-induced A3Z3 degradation. We further identified additional essential residues for FIV Vif interaction in the A3Z2 domain, allowing the generation of FIV Vif resistant A3Z2Z3. Mutated feline A3s also showed resistance to the Vif of a lion-specific FIV, indicating an evolutionary conserved Vif–A3 binding. Comparative modelling of feline A3Z2Z3 suggests that the residues interacting with FIV Vif have, unlike Vif-interacting residues in human A3s, a unique location at the domain interface of Z2 and Z3 and that the linker forms a homeobox-like domain protruding of the Z2Z3 core. HIV-2/SIV Vifs efficiently degrade feline A3Z2Z3 by possible targeting the linker stretch connecting both Z-domains.

Conclusions

Here we identified in feline A3s residues important for binding of FIV Vif and a unique protein domain insertion (linker). To understand Vif evolution, a structural model of the feline A3 was developed. Our results show that HIV Vif binds human A3s differently than FIV Vif feline A3s. The linker insertion is suggested to form a homeo-box domain, which is unique to A3s of cats and related species, and not found in human and mouse A3s. Together, these findings indicate a specific and different A3 evolution in cats and human.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-016-0274-9) contains supplementary material, which is available to authorized users.

Rapid acquisition of HPV around the time of sexual debut in adolescent girls in Tanzania

BACKGROUND

No reports exist on genotype-specific human papillomavirus (HPV) acquisition in girls after first sex in sub-Saharan Africa, despite high HPV prevalence and cervical cancer incidence.

METHODS

We followed 503 HP-unvaccinated girls aged 15-16 years in Mwanza, Tanzania, 3-monthly for 18 months with interviews and self-administered vaginal swabs. Swabs were tested for 13 higHRisk and 24 low-risk HPV genotypes. Incidence, clearance and duration of overall HPV and genotype-specific infections were calculated and associated factors evaluated.

RESULTS

A total of 106 participants reported first sex prior to enrolment (N = 29) or during follow-up (N = 77). One was HIV-positive at the final visit. The remaining 105 girls contributed 323 adequate specimens. Incidence of any new HPV genotype was 225/100 person-years (pys), and incidence of vaccine types HPV-6, -11, -16 and -18 were 12, 2, 2 and 7/100 pys, respectively. Reporting sex in the past 3 months and knowing the most recent sexual partner for a longer period before sex were associated with HPV acquisition. Median time from reported sexual debut to first HPVinfection was 5 months, and infection duration was 6 months.

CONCLUSION

This is the first description of HPV acquisition after first sex in sub-Saharan Africa where the incidence of cervical cancer is amongst the highest in the world. HPV incidence was very high after first sex, including some vaccine genotypes, and infection duration was short. This very high HPV incidence may help explain high cervical cancer rates, and supports recommendations that the HPV vaccine should be given to girls before first sex.

Assessing parallel gene histories in viral genomes

BACKGROUND

The increasing abundance of sequence data has exacerbated a long known problem: gene trees and species trees for the same terminal taxa are often incongruent. Indeed, genes within a genome have not all followed the same evolutionary path due to events such as incomplete lineage sorting, horizontal gene transfer, gene duplication and deletion, or recombination. Considering conflicts between gene trees as an obstacle, numerous methods have been developed to deal with these incongruences and to reconstruct consensus evolutionary histories of species despite the heterogeneity in the history of their genes. However, inconsistencies can also be seen as a source of information about the specific evolutionary processes that have shaped genomes.

RESULTS

The goal of the approach here proposed is to exploit this conflicting information: we have compiled eleven variables describing phylogenetic relationships and evolutionary pressures and submitted them to dimensionality reduction techniques to identify genes with similar evolutionary histories. To illustrate the applicability of the method, we have chosen two viral datasets, namely papillomaviruses and Turnip mosaic virus (TuMV) isolates, largely dissimilar in genome, evolutionary distance and biology. Our method pinpoints viral genes with common evolutionary patterns. In the case of papillomaviruses, gene clusters match well our knowledge on viral biology and life cycle, illustrating the potential of our approach. For the less known TuMV, our results trigger new hypotheses about viral evolution and gene interaction.

CONCLUSIONS

The approach here presented allows turning phylogenetic inconsistencies into evolutionary information, detecting gene assemblies with similar histories, and could be a powerful tool for comparative pathogenomics.

HPV Involvement in Head and Neck Cancers: Comprehensive Assessment of Biomarkers in 3680 Patients

BACKGROUND

We conducted a large international study to estimate fractions of head and neck cancers (HNCs) attributable to human papillomavirus (HPV-AFs) using six HPV-related biomarkers of viral detection, transcription, and cellular transformation.

METHODS

Formalin-fixed, paraffin-embedded cancer tissues of the oral cavity (OC), pharynx, and larynx were collected from pathology archives in 29 countries. All samples were subject to histopathological evaluation, DNA quality control, and HPV-DNA detection. Samples containing HPV-DNA were further subject to HPV E6*I mRNA detection and to p16(INK4a), pRb, p53, and Cyclin D1 immunohistochemistry. Final estimates of HPV-AFs were based on HPV-DNA, HPV E6*I mRNA, and/or p16(INK4a) results.

RESULTS

A total of 3680 samples yielded valid results: 1374 pharyngeal, 1264 OC, and 1042 laryngeal cancers. HPV-AF estimates based on positivity for HPV-DNA, and for either HPV E6*I mRNA or p16(INK4a), were 22.4%, 4.4%, and 3.5% for cancers of the oropharynx, OC, and larynx, respectively, and 18.5%, 3.0%, and 1.5% when requiring simultaneous positivity for all three markers. HPV16 was largely the most common type. Estimates of HPV-AF in the oropharynx were highest in South America, Central and Eastern Europe, and Northern Europe, and lowest in Southern Europe. Women showed higher HPV-AFs than men for cancers of the oropharynx in Europe and for the larynx in Central-South America.

CONCLUSIONS

HPV contribution to HNCs is substantial but highly heterogeneous by cancer site, region, and sex. This study, the largest exploring HPV attribution in HNCs, confirms the important role of HPVs in oropharyngeal cancer and drastically downplays the previously reported involvement of HPVs in the other HNCs.

High-level of viral genomic diversity in cervical cancers: A Brazilian study on human papillomavirus type 16

Invasive cervical cancer (ICC) is the third most frequent cancer among women worldwide and is associated with persistent infection by carcinogenic human papillomaviruses (HPVs). The combination of large populations of viral progeny and decades of sustained infection may allow for the generation of intra-patient diversity, in spite of the assumedly low mutation rates of PVs. While the natural history of chronic HPVs infections has been comprehensively described, within-host viral diversity remains largely unexplored. In this study we have applied next generation sequencing to the analysis of intra-host genetic diversity in ten ICC and one condyloma cases associated to single HPV16 infection. We retrieved from all cases near full-length genomic sequences. All samples analyzed contained polymorphic sites, ranging from 3 to 125 polymorphic positions per genome, and the median probability of a viral genome picked at random to be identical to the consensus sequence in the lesion was only 40%. We have also identified two independent putative duplication events in two samples, spanning the L2 and the L1 gene, respectively. Finally, we have identified with good support a chimera of human and viral DNA. We propose that viral diversity generated during HPVs chronic infection may be fueled by innate and adaptive immune pressures. Further research will be needed to understand the dynamics of viral DNA variability, differentially in benign and malignant lesions, as well as in tissues with differential intensity of immune surveillance. Finally, the impact of intralesion viral diversity on the long-term oncogenic potential may deserve closer attention.

Cancer, Warts, or Asymptomatic Infections: Clinical Presentation Matches Codon Usage Preferences in Human Papillomaviruses

Viruses rely completely on the hosts’ machinery for translation of viral transcripts. However, for most viruses infecting humans, codon usage preferences (CUPrefs) do not match those of the host. Human papillomaviruses (HPVs) are a showcase to tackle this paradox: they present a large genotypic diversity and a broad range of phenotypic presentations, from asymptomatic infections to productive lesions and cancer. By applying phylogenetic inference and dimensionality reduction methods, we demonstrate first that genes in HPVs are poorly adapted to the average human CUPrefs, the only exception being capsid genes in viruses causing productive lesions. Phylogenetic relationships between HPVs explained only a small proportion of CUPrefs variation. Instead, the most important explanatory factor for viral CUPrefs was infection phenotype, as orthologous genes in viruses with similar clinical presentation displayed similar CUPrefs. Moreover, viral genes with similar spatiotemporal expression patterns also showed similar CUPrefs. Our results suggest that CUPrefs in HPVs reflect either variations in the mutation bias or differential selection pressures depending on the clinical presentation and expression timing. We propose that poor viral CUPrefs may be central to a trade-off between strong viral gene expression and the potential for eliciting protective immune response.

Papillomaviruses: Viral evolution, cancer and evolutionary medicine

Papillomaviruses (PVs) are a numerous family of small dsDNA viruses infecting virtually all mammals. PVs cause infections without triggering a strong immune response, and natural infection provides only limited protection against reinfection. Most PVs are part and parcel of the skin microbiota. In some cases, infections by certain PVs take diverse clinical presentations from highly productive self-limited warts to invasive cancers. We propose PVs as an excellent model system to study the evolutionary interactions between the immune system and pathogens causing chronic infections: genotypically, PVs are very diverse, with hundreds of different genotypes infecting skin and mucosa; phenotypically, they display extremely broad gradients and trade-offs between key phenotypic traits, namely productivity, immunogenicity, prevalence, oncogenicity and clinical presentation. Public health interventions have been launched to decrease the burden of PV-associated cancers, including massive vaccination against the most oncogenic human PVs, as well as systematic screening for PV chronic anogenital infections. Anti-PVs vaccines elicit protection against infection, induce cross-protection against closely related viruses and result in herd immunity. However, our knowledge on the ecological and intrapatient dynamics of PV infections remains fragmentary. We still need to understand how the novel anthropogenic selection pressures posed by vaccination and screening will affect viral circulation and epidemiology. We present here an overview of PV evolution and the connection between PV genotypes and the phenotypic, clinical manifestations of the diseases they cause. This differential link between viral evolution and the gradient cancer-warts-asymptomatic infections makes PVs a privileged playground for evolutionary medicine research.

Human papillomavirus DNA prevalence and type distribution in anal carcinomas worldwide

Knowledge about human papillomaviruses (HPV) types involved in anal cancers in some world regions is scanty. Here, we describe the HPV DNA prevalence and type distribution in a series of invasive anal cancers and anal intraepithelial neoplasias (AIN) grades 2/3 from 24 countries. We analyzed 43 AIN 2/3 cases and 496 anal cancers diagnosed from 1986 to 2011. After histopathological evaluation of formalin-fixed paraffin-embedded samples, HPV DNA detection and genotyping was performed using SPF-10/DEIA/LiPA25 system (version 1). A subset of 116 cancers was further tested for p16(INK4a) expression, a cellular surrogate marker for HPV-associated transformation. Prevalence ratios were estimated using multivariate Poisson regression with robust variance in the anal cancer data set. HPV DNA was detected in 88.3% of anal cancers (95% confidence interval [CI]: 85.1-91.0%) and in 95.3% of AIN 2/3 (95% CI: 84.2-99.4%). Among cancers, the highest prevalence was observed in warty-basaloid subtype of squamous cell carcinomas, in younger patients and in North American geographical region. There were no statistically significant differences in prevalence by gender. HPV16 was the most frequent HPV type detected in both cancers (80.7%) and AIN 2/3 lesions (75.4%). HPV18 was the second most common type in invasive cancers (3.6%). p16(INK4a) overexpression was found in 95% of HPV DNA-positive anal cancers. In view of the results of HPV DNA and high proportion of p16(INK4a) overexpression, infection by HPV is most likely to be a necessary cause for anal cancers in both men and women. The large contribution of HPV16 reinforces the potential impact of HPV vaccines in the prevention of these lesions.

Novel papillomaviruses in free-ranging Iberian bats: no virus-host co-evolution, no strict host specificity, and hints for recombination

Papillomaviruses (PVs) are widespread pathogens. However, the extent of PV infections in bats remains largely unknown. This work represents the first comprehensive study of PVs in Iberian bats. We identified four novel PVs in the mucosa of free-ranging Eptesicus serotinus (EserPV1, EserPV2, and EserPV3) and Rhinolophus ferrumequinum (RferPV1) individuals and analyzed their phylogenetic relationships within the viral family. We further assessed their prevalence in different populations of E. serotinus and its close relative E. isabellinus. Although it is frequent to read that PVs co-evolve with their host, that PVs are highly species-specific, and that PVs do not usually recombine, our results suggest otherwise. First, strict virus–host co-evolution is rejected by the existence of five, distantly related bat PV lineages and by the lack of congruence between bats and bat PVs phylogenies. Second, the ability of EserPV2 and EserPV3 to infect two different bat species (E. serotinus and E. isabellinus) argues against strict host specificity. Finally, the description of a second noncoding region in the RferPV1 genome reinforces the view of an increased susceptibility to recombination in the E2-L2 genomic region. These findings prompt the question of whether the prevailing paradigms regarding PVs evolution should be reconsidered.

Ribosomal mutations affecting the translation of genes that use non-optimal codons

Genes that are laterally acquired by a new host species often contain codons that are non-optimal to the tRNA repertoire of the new host, which may lead to insufficient translational levels. Inefficient translation can be overcome by different mechanisms, such as incremental amelioration of the coding sequence, compensatory mutations in the regulatory sequences leading to increased transcription or increase in gene copy number. However, there is also a possibility that ribosomal mutations can improve the expression of such genes. To test this hypothesis, we examined the effects of point mutations in the endogenous ribosomal proteins S12 and S5 in Escherichia coli, which are known to be involved in the decoding of the mRNA, on the efficiency of translation of exogenous genes that use non-optimal codons, in vivo. We show that an S12 mutant in E. coli is able to express exogenous genes, with non-optimal codons, to higher levels than the wild-type, and explore the mechanisms underlying this phenomenon in this mutant. Our results suggest that the transient emergence of mutants that allow efficient expression of exogenous genes with non-optimal codons could also increase the chances of fixation of laterally transferred genes.