Genomic Diversity of HPV6 and HPV11 in Recurrent Respiratory Papillomatosis: Association with Malignant Transformation in the Lungs and Clinical Outcomes.

Recurrent respiratory papillomatosis (RRP) is a rare, proliferative disease caused by human papillomavirus 6 (HPV6) and HPV11. RRP can occasionally spread and undergo malignant transformation. We analysed samples across time for five RRP patients with malignant transformation and four with highly recurrent, non-malignant RRP by applying high-throughput sequencing. Patients with malignant transformation were infected by HPV11_A1/A2, while most non-malignant cases were associated with HPV6. Transient multiple infections with HPV6 and HPV11 were found in two patients, and resolved later to single infections. Viral genome loads were homogeneous across groups (median=78 viral genomes per human genome). Within-patient, we did not observe differences between the viral sequences in the papillomatous lesions and in the malignant tissue. Genetic analysis of the NLRP1 gene revealed no known mutations linked to idiopathic RRP, though some novel variants merit to be explored in larger cohorts. HPV11 infections appear associated with RRP malignant transformation in young patients. Multiple infections can occur in RRP, but within-patient viral diversity is minimal for a given genotype. Our results confirm the importance of viral genotype in disease prognosis and are consistent with growing evidence of HPV11 infections to be differentially associated with RRP malignant transformation in young patients.

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.

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.

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.

High Plasma Levels of Fibroblast Growth Factor 23 Are Associated with Increased Risk of COVID-19 in End-Stage Renal Disease Patients on Hemodialysis: Results of a Prospective Cohort.

End-stage renal disease (ESRD) patients are a population with high rates of COVID-19 and mortality. These patients present a low response to anti-SARS-CoV-2 immunization, which is associated with immune dysfunction. ESRD patients also present high plasma titers of Fibroblast Growth Factor 23 (FGF23), a protein hormone that reduces immune response in vivo and in vitro. Increased FGF23 levels associate with higher infection-related hospitalizations and adverse infectious outcomes. Thus, we evaluated whether ESRD patients with high FGF23 titers have an increased rate of SARS-CoV-2 infection. METHODS: We performed a prospective cohort of ESRD patients in hemodialysis who had measurements of plasma intact FGF23 in 2019. We determined COVID-19 infections, hospitalizations, and mortality between January 2020 and December 2021. RESULTS: We evaluated 243 patients. Age: 60.4 ± 10.8 years. Female: 120 (49.3%), diabetes: 110 (45.2%). During follow-up, 45 patients developed COVID-19 (18.5%), 35 patients were hospitalized, and 12 patients died (mortality rate: 26.6%). We found that patients with higher FGF23 levels (defined as equal or above median) had a higher rate of SARS-CoV-2 infection versus those with lower levels (18.8% versus 9.9%; Hazard ratio: 1.92 [1.03-3.56], p = 0.039). Multivariate analysis showed that increased plasma FGF23 was independently associated with SARS-CoV-2 infection and severe COVID-19. DISCUSSION: Our results suggest that high plasma FGF23 levels are a risk factor for developing COVID-19 in ESRD patients. These data support the potential immunosuppressive effects of high circulating FGF23 as a factor implicated in the association with worse clinical outcomes. Further data are needed to confirm this hypothesis.

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 menstrual products 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 organizational levels. Here, we study how exposure to different types of menstrual 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 menstrual product type 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.

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.

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.

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.

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.

Dynamic Insulin Basal Needs Estimation and Parameters Adjustment in Type 1 Diabetes.

Technology advances have made possible improvements such as Continuous Glucose Monitors, giving the patient a glucose reading every few minutes, or insulin pumps, allowing more personalized therapies. With the increasing number of available closed-loop systems, new challenges appear regarding algorithms and functionalities. Several of the analysed systems in this paper try to adapt to changes in some patients' conditions and, in several of these systems, other variables such as basal needs are considered fixed from day to day to simplify the control problem. Therefore, these systems require a correct adjustment of the basal needs profile which becomes crucial to obtain good results. In this paper a novel approach tries to dynamically determine the insulin basal needs of the patient and use this information within a closed-loop algorithm, allowing the system to dynamically adjust in situations of illness, exercise, high-fat-content meals or even partially blocked infusion sites and avoiding the need for setting a basal profile that approximately matches the basal needs of the patient. The insulin sensitivity factor and the glycemic target are also dynamically modified according to the situation of the patient. Basal insulin needs are dynamically determined through linear regression via the decomposition of previously dosed insulin and its effect on the patient's glycemia. Using the obtained value as basal insulin needs and other mechanisms such as basal needs modification through its trend, ISF and glycemic targets modification and low-glucose-suspend threshold, the safety of the algorithm is improved. The dynamic basal insulin needs determination was successfully included in a closed-loop control algorithm and was simulated on 30 virtual patients (10 adults, 10 adolescent and 10 children) using an open-source python implementation of the FDA-approved (Food and Drug Administration) UVa (University of Virginia)/Padova Simulator. Simulations showed that the proposed system dynamically determines the basal needs and can adapt to a partial blockage of the insulin infusion, obtaining similar results in terms of time in range to the case in which no blockage was simulated. The proposed algorithm can be incorporated to other current closed-loop control algorithms to directly estimate the patient's basal insulin needs or as a monitoring channel to detect situations in which basal needs may differ from the expected ones.

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.

HPV cervical infections and serological status in vaccinated and unvaccinated women.

Understanding genital infections by Human papillomaviruses (HPVs) remains a major public health issue, especially in countries where vaccine uptake is low. We investigate HPV prevalence and antibody status in 150 women (ages 18 to 25) in Montpellier, France. At inclusion and one month later, cervical swabs, blood samples and questionnaires (for demographics and behavioural variables) were collected. Oncogenic, non-vaccine genotypes HPV51, HPV66, HPV53, and HPV52 were the most frequently detected viral genotypes overall. Vaccination status, which was well-balanced in the cohort, showed the strongest (protective) effect against HPV infections, with an associated odds ratio for alphapapillomavirus detection of 0.45 (95% confidence interval: [0.22;0.58]). We also identified significant effects of age, number of partners, body mass index, and contraception status on HPV detection and on coinfections. Type-specific IgG serological status was also largely explained by the vaccination status. IgM seropositivity was best explained by HPV detection at inclusion only. Finally, we identify a strong significant effect of vaccination on genotype prevalence, with a striking under-representation of HPV51 in vaccinated women. Variations in HPV prevalence correlate with key demographic and behavioural variables. The cross-protective effect of the vaccine against HPV51 merits further investigation.

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.