Multicenter benchmarking of short and long read wet lab protocols for clinical viral metagenomics.

Metagenomics is gradually being implemented for diagnosing infectious diseases. However, in-depth protocol comparisons for viral detection have been limited to individual sets of experimental workflows and laboratories. In this study, we present a benchmark of metagenomics protocols used in clinical diagnostic laboratories initiated by the European Society for Clinical Virology (ESCV) Network on NGS (ENNGS). A mock viral reference panel was designed to mimic low biomass clinical specimens. The panel was used to assess the performance of twelve metagenomic wet lab protocols currently in use in the diagnostic laboratories of participating ENNGS member institutions. Both Illumina and Nanopore, shotgun and targeted capture probe protocols were included. Performance metrics sensitivity, specificity, and quantitative potential were assessed using a central bioinformatics pipeline. Overall, viral pathogens with loads down to 104 copies/ml (corresponding to CT values of 31 in our PCR assays) were detected by all the evaluated metagenomic wet lab protocols. In contrast, lower abundant mixed viruses of CT values of 35 and higher were detected only by a minority of the protocols. Considering the reference panel as the gold standard, optimal thresholds to define a positive result were determined per protocol, based on the horizontal genome coverage. Implementing these thresholds, sensitivity and specificity of the protocols ranged from 67 to 100 % and 87 to 100 %, respectively. A variety of metagenomic protocols are currently in use in clinical diagnostic laboratories. Detection of low abundant viral pathogens and mixed infections remains a challenge, implying the need for standardization of metagenomic analysis for use in clinical settings.

Evaluation of human papillomavirus DNA in colorectal cancer and adjacent mucosal tissue samples.

BACKGROUND: Although the role of viral agents, such as human papillomavirus (e.g. HPV16, HPV18) in colorectal cancer (CRC) has been previously investigated, results remain inconclusive. METHODS: To further evaluate the involvement of oncogenic HPV types in CRC, 40 frozen neoplastic and 40 adjacent colonic tissues collected from Italian patients were analyzed by Luminex-based assays that detect a broad spectrum of HPV types, i.e. Alpha (n = 21), Beta (n = 46) and Gamma HPVs (n = 52). In addition, 125 frozen CRC samples and 70 surrounding mucosal tissues were collected from Czech patients and analyzed by broad spectrum PCR protocols: (i) FAP59/64, (ii) FAPM1 and (iii) CUT combined with Next Generation Sequencing (NGS). RESULTS: Using Luminex-basedassays, DNA from HPV16 was detected in 5% (2/40) CRC tissues from Italian patients. One HPV16 DNA-positive CRC case was subsequently confirmed positive for E6*I mRNA. Cutaneous beta HPV types were detected in 10% (4/40) adjacent tissues only, namely HPV111 (n = 3) and HPV120 (n = 1), while gamma HPV168 (n = 1) and HPV199 (n = 1) types were detected in adjacent and in tumor tissues, respectively. The NGS analysis of the CRC Czech samples identified HPV sequences from mucosal alpha-3 (HPV89), alpha-7 (HPV18, 39, 68 and 70) and alpha-10 species (HPV11), as well as cutaneous beta-1 (HPV20, 24, 93, 98, 105,124) beta-2 (HPV23), beta-3 (HPV49) and gamma-1 species (HPV205). CONCLUSIONS: Our findings indicate that HPV types belonging to the mucosal alpha, and the 'cutaneous' beta and gamma genera can be detected in the colonic mucosal samples with a low prevalence rate and a low number of HPV reads by Luminex and NGS, respectively. However, additional studies are required to corroborate these findings.

HPV38 impairs UV-induced transcriptional activation of the IL-18 pro-inflammatory cytokine.

IMPORTANCE: Here, we demonstrate that the direct binding of p53 on the IL-18 promoter region regulates its gene expression. However, the presence of E6 and E7 from human papillomavirus type 38 impairs this mechanism via a new inhibitory complex formed by DNA methyltransferase 1 (DNMT1)/PKR/ΔNp73α, which binds to the region formerly occupied by p53 in primary keratinocytes.

The E2F4/p130 Repressor Complex Cooperates with Oncogenic ΔNp73α To Inhibit Gene Expression in Human Papillomavirus 38 E6/E7-Transformed Keratinocytes and in Cancer Cells.

Tumor suppressor p53 and its related proteins, p63 and p73, can be synthesized as multiple isoforms lacking part of the N- or C-terminal regions. Specifically, high expression of the ΔNp73α isoform is notoriously associated with various human malignancies characterized by poor prognosis. This isoform is also accumulated by oncogenic viruses, such as Epstein-Barr virus (EBV), as well as genus beta human papillomaviruses (HPV) that appear to be involved in carcinogenesis. To gain additional insight into ΔNp73α mechanisms, we have performed proteomics analyses using human keratinocytes transformed by the E6 and E7 proteins of the beta-HPV type 38 virus as an experimental model (38HK). We find that ΔNp73α associates with the E2F4/p130 repressor complex through a direct interaction with E2F4. This interaction is favored by the N-terminal truncation of p73 characteristic of ΔNp73 isoforms. Moreover, it is independent of the C-terminal splicing status, suggesting that it could represent a general feature of ΔNp73 isoforms (α, ß, γ, δ, ε, ζ, θ, η, and η1). We show that the ΔNp73α-E2F4/p130 complex inhibits the expression of specific genes, including genes encoding for negative regulators of proliferation, both in 38HK and in HPV-negative cancer-derived cell lines. Such genes are not inhibited by E2F4/p130 in primary keratinocytes lacking ΔNp73α, indicating that the interaction with ΔNp73α rewires the E2F4 transcriptional program. In conclusion, we have identified and characterized a novel transcriptional regulatory complex with potential implications in oncogenesis. IMPORTANCE The TP53 gene is mutated in about 50% of human cancers. In contrast, the TP63 and TP73 genes are rarely mutated but rather expressed as ΔNp63 and ΔNp73 isoforms in a wide range of malignancies, where they act as p53 antagonists. Accumulation of ΔNp63 and ΔNp73, which is associated with chemoresistance, can result from infection by oncogenic viruses such as EBV or HPV. Our study focuses on the highly carcinogenic ΔNp73α isoform and uses a viral model of cellular transformation. We unveil a physical interaction between ΔNp73α and the E2F4/p130 complex involved in cell cycle control, which rewires the E2F4/p130 transcriptional program. Our work shows that ΔNp73 isoforms can establish interactions with proteins that do not bind to the TAp73α tumor suppressor. This situation is analogous to the gain-of-function interactions of p53 mutants supporting cellular proliferation.

Integration of Sequencing and Epidemiologic Data for Surveillance of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infections in a Tertiary-Care Hospital.

BACKGROUND: The ongoing coronavirus disease 2019 pandemic significantly burdens hospitals and other healthcare facilities. Therefore, understanding the entry and transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is critical for effective prevention and preparedness measures. We performed surveillance and analysis of testing and transmission of SARS-CoV-2 infections in a tertiary-care hospital in Germany during the second and third pandemic waves in fall/winter 2020. METHODS: Between calendar week 41 in 2020 and calendar week 1 in 2021, 40%, of all positive patient and staff samples (284 total) were subjected to full-length viral genome sequencing. Clusters were defined based on similar genotypes indicating common sources of infection. We integrated phylogenetic, spatial, and temporal metadata to detect nosocomial infections and outbreaks, uncover transmission chains, and evaluate containment measures' effectiveness. RESULTS: Epidemiologic data and contact tracing readily recognize most healthcare-associated (HA) patient infections. However, sequencing data reveal that temporally preceding index cases and transmission routes can be missed using epidemiologic methods, resulting in delayed interventions and serially linked outbreaks being counted as independent events. While hospital-associated transmissions were significantly elevated at a moderate rate of community transmission during the second wave, systematic testing and high vaccination rates among staff have led to a substantial decrease in HA infections at the end of the second/beginning of the third wave despite high community transmissions. CONCLUSIONS: While epidemiologic analysis is critical for immediate containment of HA SARS-CoV-2 outbreaks, integration of genomic surveillance revealed weaknesses in identifying staff contacts. Our study underscores the importance of high testing frequency and genomic surveillance to detect, contain and prevent SARS-CoV-2-associated infections in healthcare settings.

Remdesivir-induced emergence of SARS-CoV2 variants in patients with prolonged infection.

We here investigate the impact of antiviral treatments such as remdesivir on intra-host genomic diversity and emergence of SARS-CoV2 variants in patients with a prolonged course of infection. Sequencing and variant analysis performed in 112 longitudinal respiratory samples from 14 SARS-CoV2-infected patients with severe disease progression show that major frequency variants do not generally arise during prolonged infection. However, remdesivir treatment can increase intra-host genomic diversity and result in the emergence of novel major variant species harboring fixed mutations. This is particularly evident in a patient with B cell depletion who rapidly developed mutations in the RNA-dependent RNA polymerase gene following remdesivir treatment. Remdesivir treatment-associated emergence of novel variants is of great interest in light of current treatment guidelines for hospitalized patients suffering from severe SARS-CoV2 disease, as well as the potential use of remdesivir to preventively treat non-hospitalized patients at high risk for severe disease progression.

Clinical Evaluation of a Fully-Automated High-Throughput Multiplex Screening-Assay to Detect and Differentiate the SARS-CoV-2 B.1.1.529 (Omicron) and B.1.617.2 (Delta) Lineage Variants.

BACKGROUND: The recently emerged SARS-CoV-2 B.1.1.529 lineage and its sublineages (Omicron variant) pose a new challenge to healthcare systems worldwide due to its ability to efficiently spread in immunized populations and its resistance to currently available monoclonal antibody therapies. RT-PCR-based variant tests can be used to screen large sample-sets rapidly and accurately for relevant variants of concern (VOC). The aim of this study was to establish and validate a multiplex assay on the cobas 6800/8800 systems to allow discrimination between the two currently circulating VOCs, Omicron and Delta, in clinical samples. METHODS: Primers and probes were evaluated for multiplex compatibility. Analytic performance was assessed using cell culture supernatant of an Omicron variant isolate and a clinical Delta variant sample, normalized to WHO-Standard. Clinical performance of the multiplex assay was benchmarked against NGS results. RESULTS: In silico testing of all oligos showed no interactions with a high risk of primer-dimer formation or amplification of human DNA/RNA. Over 99.9% of all currently available Omicron variant sequences are a perfect match for at least one of the three Omicron targets included in the multiplex. Analytic sensitivity was determined as 19.0 IU/mL (CI95%: 12.9-132.2 IU/mL) for the A67V + del-HV69-70 target, 193.9 IU/mL (CI95%: 144.7-334.7 IU/mL) for the E484A target, 35.5 IU/mL (CI95%: 23.3-158.0 IU/mL) for the N679K + P681H target and 105.0 IU/mL (CI95%: 80.7-129.3 IU/mL) for the P681R target. All sequence variances were correctly detected in the clinical sample set (225/225 Targets). CONCLUSION: RT-PCR-based variant screening compared to whole genome sequencing is both rapid and reliable in detecting relevant sequence variations in SARS-CoV-2 positive samples to exclude or verify relevant VOCs. This allows short-term decision-making, e.g., for patient treatment or public health measures.

Comparing susceptibility and contagiousness in concurrent outbreaks with a non-VOC and the VOC SARS-CoV-2 variant B.1.1.7 in daycare centers in Hamburg, Germany.

We describe two outbreaks of SARS-CoV-2 in daycare centers in the metropolitan area of Hamburg, Germany. The outbreaks occurred in rapid chronological succession, in neighborhoods with a very similar sociodemographic structure, thus allowing for cross-comparison of these events. We combined classical and molecular epidemiologic investigation methods to study infection entry, spread within the facilities, and subsequent transmission of infections to households. Epidemiologic and molecular evidence suggests a superspreading event with a non-variant of concern (non-VOC) SARS CoV-2 strain at the root of the first outbreak. The second outbreak involved two childcare facilities experiencing infection activity with the variant of concern (VOC) B.1.1.7 (Alpha). We show that the index cases in all outbreaks had been childcare workers, and that children contributed substantially to secondary transmission of SARS-CoV-2 infection from childcare facilities to households. The frequency of secondary transmissions in households originating from B.1.1.7-infected children was increased compared to children with non-VOC infections. Self-reported symptoms, particularly cough and rhinitis, occurred more frequently in B.1.1.7-infected children. Especially in light of the rapidly spreading VOC B.1.617.2 (Delta), our data underline the notion that rigorous SARS-CoV-2 testing in combination with screening of contacts regardless of symptoms is an important measure to prevent SARS-CoV-2 infection of unvaccinated individuals in daycare centers and associated households.

Rapid Automated Screening for SARS-CoV-2 B.1.617 Lineage Variants (Delta/Kappa) through a Versatile Toolset of qPCR-Based SNP Detection.

BACKGROUND: The recent emergence of distinct and highly successful SARS-CoV-2 lineages has substantial implications for individual patients and public health measures. While next-generation-sequencing is routinely performed for surveillance purposes, RT-qPCR can be used to rapidly rule-in or rule-out relevant variants, e.g., in outbreak scenarios. The objective of this study was to create an adaptable and comprehensive toolset for multiplexed Spike-gene SNP detection, which was applied to screen for SARS-CoV-2 B.1.617 lineage variants. METHODS: We created a broad set of single nucleotide polymorphism (SNP)-assays including del-Y144/145, E484K, E484Q, P681H, P681R, L452R, and V1176F based on a highly specific multi-LNA (locked nucleic acid)-probe design to maximize mismatch discrimination. As proof-of-concept, a multiplex-test was compiled and validated (SCOV2-617VOC-UCT) including SNP-detection for L452R, P681R, E484K, and E484Q to provide rapid screening capabilities for the novel B.1.617 lineages. RESULTS: For the multiplex-test (SCOV2-617VOC-UCT), the analytic lower limit of detection was determined as 182 IU/mL for L452R, 144 IU/mL for P681R, and 79 IU/mL for E484Q. A total of 233 clinical samples were tested with the assay, including various on-target and off-target sequences. All SNPs (179/179 positive) were correctly identified as determined by SARS-CoV-2 whole genome sequencing. CONCLUSION: The recurrence of SNP locations and flexibility of methodology presented in this study allows for rapid adaptation to current and future variants. Furthermore, the ability to multiplex various SNP-assays into screening panels improves speed and efficiency for variant testing. We show 100% concordance with whole genome sequencing for a B.1.617.2 screening assay on the cobas6800 high-throughput system.

MinION nanopore sequencing and assembly of a complete human papillomavirus genome.

BACKGROUND: The MinION sequencer belongs to the third generation of sequencing technology that allows for the generation of ultra-long reads, representing a potentially more effective approach to characterize entire viral genome sequences than other time-consuming and low-throughput methodologies. METHODS: We report the use of the MinION nanopore sequencer to sequence the full-length genome of human papillomavirus (HPV)-ICB2 (7441 bp), which was previously characterized in our laboratory. Three independent MinION libraries were prepared and sequenced using either three consecutive 12 -h runs (Protocol A) or a single run of 48 h starting from a pool of three barcoded DNA libraries (Protocol B). A fully automated bioinformatics pipeline was developed for the reconstruction of the viral genome. RESULTS: Protocols A and B generated 9,354,933 and 3,255,879 reads, respectively. Read length N50 values ranged between 6976 and 7360 nucleotides over the four sequencing runs. Bioinformatics analysis showed that both protocols allowed for the reconstruction of the whole viral genome, with pairwise percentages of identity to HPV-ICB2 of 100 % for protocol A and 99.98 % for protocol B. CONCLUSION: Our results show that the use of the MinION nanopore sequencer represents an effective strategy for whole-genome sequencing of HPVs with a minimal error rate.

SARS-CoV-2 Reinfection in a Healthcare Worker Despite the Presence of Detectable Neutralizing Antibodies.

So far, only a few reports about reinfections with SARS-CoV-2 have been published, and they often lack detailed immunological and virological data. We report about a SARS-CoV-2 reinfection with a genetically distinct SARS-CoV-2 variant in an immunocompetent female healthcare worker that has led to a mild disease course. No obvious viral escape mutations were observed in the second virus variant. The infectious virus was shed from the patient during the second infection episode despite the presence of neutralizing antibodies in her blood. Our data indicate that a moderate immune response after the first infection, but not a viral escape, did allow for reinfection and live virus shedding.

Diversity of human papillomavirus in the anal canal of HIV-positive and HIV-negative men.

BACKGROUND: To characterize the HPV diversity in the anal mucosa of men with different sexual behavior and HIV status by next-generation sequencing (NGS). METHODS: Anal swabs from HIV-positive (n = 94; mean age, 38 years) and HIV-negative (n = 100; mean age, 37.5 years) men who have sex with men (MSM) and HIV-negative men (predominantly men who have sex with women, MSW) (n = 99; mean age, 38.2 years) were analyzed by broad-spectrum PCR protocols combined with NGS. FINDINGS: Alpha HPV types (n = 74) were detected mainly in the MSM groups (HPV6, 11, and 43 were the most abundant types) compared with MSW (n = 16) (HPV11, 32, and 87 were among the most abundant). In contrast, beta HPVs were more abundantly detected among MSW (n = 45) than in the HIV-positive (n = 16) and HIV-negative (n = 26) MSM groups. Gamma HPVs were detected almost equally in HIV-positive MSM (n = 62), HIV-negative MSM (n = 58), and MSW (n = 57). In addition, 31 putative novel PV types were identified. CONCLUSIONS: Our data show that beta and gamma HPV types are present in the anal mucosa, thus reinforcing the existing evidence that they can be detected at anatomical sites other than skin. Alpha and beta HPV distribution among these three groups appears to vary according to sexual behavior.

SARS-CoV-2 outbreak investigation in a German meat processing plant.

We describe a multifactorial investigation of a SARS-CoV-2 outbreak in a large meat processing complex in Germany. Infection event timing, spatial, climate and ventilation conditions in the processing plant, sharing of living quarters and transport, and viral genome sequences were analyzed. Our results suggest that a single index case transmitted SARS-CoV-2 to co-workers over distances of more than 8 m, within a confined work area in which air is constantly recirculated and cooled. Viral genome sequencing shows that all cases share a set of mutations representing a novel sub-branch in the SARS-CoV-2 C20 clade. We identified the same set of mutations in samples collected in the time period between this initial infection cluster and a subsequent outbreak within the same factory, with the largest number of confirmed SARS-CoV-2 cases in a German meat processing facility reported so far. Our results indicate climate conditions, fresh air exchange rates, and airflow as factors that can promote efficient spread of SARS-CoV-2 via long distances and provide insights into possible requirements for pandemic mitigation strategies in industrial workplace settings.

Human papillomavirus type 38 alters wild-type p53 activity to promote cell proliferation via the downregulation of integrin alpha 1 expression.

Tumor suppressors can exert pro-proliferation functions in specific contexts. In the beta human papillomavirus type 38 (HPV38) experimental model, the viral proteins E6 and E7 promote accumulation of a wild-type (WT) p53 form in human keratinocytes (HKs), promoting cellular proliferation. Inactivation of p53 by different means strongly decreases the proliferation of HPV38 E6/E7 HKs. This p53 form is phosphorylated at S392 by the double-stranded RNA-dependent protein kinase PKR, which is highly activated by HPV38. PKR-mediated S392 p53 phosphorylation promotes the formation of a p53/DNMT1 complex, which inhibits expression of integrin alpha 1 (ITGA1), a repressor of epidermal growth factor receptor (EGFR) signaling. Ectopic expression of ITGA1 in HPV38 E6/E7 HKs promotes EGFR degradation, inhibition of cellular proliferation, and cellular death. Itga1 expression was also inhibited in the skin of HPV38 transgenic mice that have an elevated susceptibility to UV-induced skin carcinogenesis. In summary, these findings reveal the existence of a specific WT p53 form that displays pro-proliferation properties.

Oral Infection by Mucosal and Cutaneous Human Papillomaviruses in the Men Who Have Sex with Men from the OHMAR Study.

Both mucosal and cutaneous Human Papillomaviruses (HPVs) can be detected in the oral cavity, but investigations regarding the epidemiology of cutaneous HPVs at this site are scarce. We assessed mucosal (alpha) and cutaneous (beta and gamma) HPV infection in oral samples of HIV-infected and uninfected men who have sex with men (MSM). Oral rinse-and-gargles were collected from 310 MSM. Alpha HPVs were detected using the Linear Array, whereas beta and gamma HPVs were detected using multiplex PCR and Luminex technology. An amplicon-based next-generation sequencing (NGS) protocol was applied to a subset of samples collected from 30 HIV-uninfected and 30 HIV-infected MSM. Beta HPVs were significantly more common than alpha types (53.8% vs. 23.9% for HIV-infected subjects, p < 0.0001; 50.3% vs. 17.1% for HIV-uninfected subjects, p < 0.0001). Gamma HPVs were also frequently detected (30.8% and 25.9% in HIV-infected and uninfected MSM, respectively). NGS produced 2,620,725 reads representative of 146 known HPVs (16 alpha-PVs, 53 beta-PVs, 76 gamma-PVs, one unclassified) and eight putative new HPVs, taxonomically assigned to the beta genus. The oral cavity contains a wide spectrum of HPVs, with beta types representing the predominant genus. The prevalence of beta and gamma HPVs is high even in immunorestored HIV-infected individuals. NGS confirmed the abundance of cutaneous HPVs and identified some putative novel beta HPVs. This study confirms that cutaneous HPVs are frequently present at mucosal sites and highlights that their pathological role deserves further investigation since it may not be limited to skin lesions.

Transforming Properties of Beta-3 Human Papillomavirus E6 and E7 Proteins.

The beta human papillomaviruses (HPVs) are subdivided into 5 species (beta-1 to beta-5), and they were first identified in the skin. However, the beta-3 species appears to be more highly represented in the mucosal epithelia than in the skin. Functional studies have also highlighted that beta-3 HPV49 shares some functional similarities with mucosal high-risk (HR) HPV16. Here, we describe the characterization of the in vitro transforming properties of the entire beta-3 species, which includes three additional HPV types: HPV75, HPV76, and HPV115. HPV49, HPV75, and HPV76 E6 and E7 (E6/E7), but not HPV115 E6 and E7, efficiently inactivate the p53 and pRb pathways and immortalize or extend the life span of human foreskin keratinocytes (HFKs). As observed for HR HPV16, cell cycle deregulation mediated by beta-3 HPV E6/E7 expression leads to p16INK4a accumulation, whereas no p16INK4a was detected in beta-2 HPV38 E6/E7 HFKs. As shown for HPV49 E6, HPV75 and HPV76 E6s degrade p53 by an E6AP/proteasome-mediated mechanism. Comparative analysis of cellular gene expression patterns of HFKs containing E6 and E7 from HR HPV16, beta-3 HPV types, and beta-2 HPV38 further highlights the functional similarities of HR HPV16 and beta-3 HPV49, HPV75, and HPV76. The expression profiles of these four HPV HFKs show some similarities and diverge substantially from those of beta-3 HPV115 E6/E7 and beta-2 HPV38 E6/E7 HFKs. In summary, our data show that beta-3 HPV types share some mechanisms with HR HPV types and pave the way for additional studies aiming to evaluate their potential role in human pathologies.IMPORTANCE Human papillomaviruses are currently classified in different genera. Mucosal HPVs belonging to the alpha genus have been clearly associated with carcinogenesis of the mucosal epithelium at different sites. Beta HPV types have been classified as cutaneous. Although findings indicate that some beta HPVs from species 1 and 2 play a role, together with UV irradiation, in skin cancer, very little is known about the transforming properties of most of the beta HPVs. This report shows the transforming activity of E6 and E7 from beta-3 HPV types. Moreover, it highlights that beta-3 HPVs share some biological properties more extensively with mucosal high-risk HPV16 than with beta-2 HPV38. This report provides new paradigms for a better understanding of the biology of the different HPV types and their possible association with lesions at mucosal and/or cutaneous epithelia.

PVAmpliconFinder: a workflow for the identification of human papillomaviruses from high-throughput amplicon sequencing.

BACKGROUND: The detection of known human papillomaviruses (PVs) from targeted wet-lab approaches has traditionally used PCR-based methods coupled with Sanger sequencing. With the introduction of next-generation sequencing (NGS), these approaches can be revisited to integrate the sequencing power of NGS. Although computational tools have been developed for metagenomic approaches to search for known or novel viruses in NGS data, no appropriate tool is available for the classification and identification of novel viral sequences from data produced by amplicon-based methods. RESULTS: We have developed PVAmpliconFinder, a data analysis workflow designed to rapidly identify and classify known and potentially new Papillomaviridae sequences from NGS amplicon sequencing with degenerate PV primers. Here, we describe the features of PVAmpliconFinder and its implementation using biological data obtained from amplicon sequencing of human skin swab specimens and oral rinses from healthy individuals. CONCLUSIONS: PVAmpliconFinder identified putative new HPV sequences, including one that was validated by wet-lab experiments. PVAmpliconFinder can be easily modified and applied to other viral families. PVAmpliconFinder addresses a gap by providing a solution for the analysis of NGS amplicon sequencing, increasingly used in clinical research. The PVAmpliconFinder workflow, along with its source code, is freely available on the GitHub platform: https://github.com/IARCbioinfo/PVAmpliconFinder.

Detection of a large spectrum of viral infections in conjunctival premalignant and malignant lesions.

To study the interaction between HIV and other carcinogenic infections in conjunctival squamous cell carcinoma (SCC), we evaluated the presence of a broad spectrum of human viruses in conjunctiva specimens. Beta Human papillomavirus (HPV; n = 46), gamma HPV (n = 52), polyomaviruses (n = 12) and herpes viruses (n = 3) was determined in DNA extracted from 67 neoplastic and 55 non-neoplastic conjunctival tissues of HIV-positive and HIV negative subjects by Luminex-based assays. Next-generation sequencing (NGS) was also used to further characterize the presence of cutaneous HPVs. Detection of beta-2 HPV infections was associated with the risk of neoplasia (adjusted odds ratio [aOR] 3.0; 95% confidence interval [CI] 1.3-6.8), regardless of HIV status (HIV positive, aOR 2.6, 95% CI 0.9-7.7; HIV negative, aOR 3.5, 95% CI 0.9-14.4). EBV was strongly associated with the risk of neoplasia (aOR 12.0, 95% CI 4.3-33.5; P < .01) mainly in HIV individuals (HIV positive, aOR 57.5; 95% CI: 10.1-327.1; HIV negative aOR 2.6; 95% CI: 0.2-34.7). NGS allowed to identify 13 putative novel HPVs in cases and controls. Our findings suggest a role of beta HPV types and EBV, in conjunctival SCC. However, additional studies of viral expression in tumor tissue are required to confirm the causal association.

Detection of human papillomaviruses in paired healthy skin and actinic keratosis by next generation sequencing.

Actinic keratosis (AK) arises on photo-damaged skin and is considered to be the precursor lesion of cutaneous squamous cell carcinoma (cSCC). Many findings support the involvement of ß human papillomaviruses (HPVs) in cSCC, while very little is known on γ HPV types. The objective of this study was to characterize the spectrum of PV types in healthy skin (HS) and AK samples of the same immunocompetent individuals using next generation sequencing (NGS). Viral DNA of 244 AK and 242 HS specimens were amplified by PCR using two different sets of primers (FAP59/64 and FAPM1). Purified amplicons were pooled and sequenced using NGS. The study resulted in the identification of a large number of known ß and γ PV types. In addition, 27 putative novel ß and 16 γ and 4 unclassified PVs were isolated. HPV types of species γ-1 (e.g. HPV4) appeared to be strongly enriched in AK versus HS. The NGS analysis revealed that a large spectrum of known and novel PVs is present in HS and AK. The evidence that species γ-1 HPV types appears to be enriched in AK in comparison to HS warrants further studies to evaluate their role in development of skin (pre)cancerous lesions.

Merkel Cell Polyomavirus Downregulates N-myc Downstream-Regulated Gene 1, Leading to Cellular Proliferation and Migration.

Merkel cell polyomavirus (MCPyV) is the first human polyomavirus etiologically associated with Merkel cell carcinoma (MCC), a rare and aggressive form of skin cancer. Similar to other polyomaviruses, MCPyV encodes early T antigen genes, viral oncogenes required for MCC tumor growth. To identify the unique oncogenic properties of MCPyV, we analyzed the gene expression profiles in human spontaneously immortalized keratinocytes (NIKs) expressing the early genes from six distinct human polyomaviruses (PyVs), including MCPyV. A comparison of the gene expression profiles revealed 28 genes specifically deregulated by MCPyV. In particular, the MCPyV early gene downregulated the expression of the tumor suppressor gene N-myc downstream-regulated gene 1 (NDRG1) in MCPyV gene-expressing NIKs and hTERT-MCPyV gene-expressing human keratinocytes (HK) compared to their expression in the controls. In MCPyV-positive MCC cells, the expression of NDRG1 was downregulated by the MCPyV early gene, as T antigen knockdown rescued the level of NDRG1. In addition, NDRG1 overexpression in hTERT-MCPyV gene-expressing HK or MCC cells resulted in a decrease in the number of cells in S phase and cell proliferation inhibition. Moreover, a decrease in wound healing capacity in hTERT-MCPyV gene-expressing HK was observed. Further analysis revealed that NDRG1 exerts its biological effect in Merkel cell lines by regulating the expression of the cyclin-dependent kinase 2 (CDK2) and cyclin D1 proteins. Overall, NDRG1 plays an important role in MCPyV-induced cellular proliferation.IMPORTANCE Merkel cell carcinoma was first described in 1972 as a neuroendocrine tumor of skin, most cases of which were reported in 2008 to be caused by a PyV named Merkel cell polyomavirus (MCPyV), the first PyV linked to human cancer. Thereafter, numerous studies have been conducted to understand the etiology of this virus-induced carcinogenesis. However, it is still a new field, and much work is needed to understand the molecular pathogenesis of MCC. In the current work, we sought to identify the host genes specifically deregulated by MCPyV, as opposed to other PyVs, in order to better understand the relevance of the genes analyzed on the biological impact and progression of the disease. These findings open newer avenues for targeted drug therapies, thereby providing hope for the management of patients suffering from this highly aggressive cancer.