Emilie Hultin

Deep sequencing extends the diversity of human papillomaviruses in human skin

Most viruses in human skin are known to be human papillomaviruses (HPVs). Previous sequencing of skin samples has identified 273 different cutaneous HPV types, including 47 previously unknown types. In the present study, we wished to extend prior studies using deeper sequencing. This deeper sequencing without prior PCR of a pool of 142 whole genome amplified skin lesions identified 23 known HPV types, 3 novel putative HPV types and 4 non-HPV viruses. The complete sequence was obtained for one of the known putative types and almost the complete sequence was obtained for one of the novel putative types. In addition, sequencing of amplimers from HPV consensus PCR of 326 skin lesions detected 385 different HPV types, including 226 previously unknown putative types. In conclusion, metagenomic deep sequencing of human skin samples identified no less than 396 different HPV types in human skin, out of which 229 putative HPV types were previously unknown.

Human papillomavirus type 197 is commonly present in skin tumors.

Non‐melanoma skin cancers commonly contain Human Papillomavirus (HPV), but the types found have varied depending on the polymerase chain reaction (PCR) primer systems used. Whole genome amplified DNA (not amplified by any specific PCR primers) from 91 skin lesions [41 squamous cell skin carcinomas (SCCs), 8 keratoacanthomas, 22 actinic keratoses, 3 basal cell carcinomas and 17 SCCs in situ] were sequenced. All samples were sequenced both at 160 Mb and 1.8 Gb sequencing depth per sample. The sequences from 10 different HPVs in 47/91 specimens were found. Sequences represented four established HPV types (HPV types 16, 22, 120, 124), two previously known putative types (present in GenBank) and four previously unknown HPV sequences (new putative types). The most commonly detected virus was cloned, sequenced and designated as HPV197. Type‐specific real‐time PCR detected HPV197 in 34/91 specimens. For comparison, a pool of the same samples after general primer PCR amplification was also sequenced. This revealed 40 different HPVs, but only two HPV types were detected both with sequencing without prior PCR and with sequencing PCR amplicons, suggesting that sequencing without prior PCR gives a more unbiased representation of the HPVs present. In summary, it was found that HPV can be sequenced from most skin disease specimens and HPV197 appeared to be the most commonly present virus.

Next generation sequencing for human papillomavirus genotyping

BACKGROUND Human papillomavirus (HPV) genotyping using next generation sequencing (NGS) could be useful to study the HPV variant-specific epidemiology, including monitoring for possible emergence of new HPV variants after introduction of HPV vaccination programs. OBJECTIVES We wished to design and validate a method for rapid HPV detection, typing and sequencing in clinical samples. STUDY DESIGN Plasmids of 15 different HPV types were mixed and serially diluted in human DNA in concentrations ranging from 1 to 100 copies per sample, amplified using the HPV general PCR primer pair PGMY and sequenced using 454 technology. Sixty cervical samples were tested both with the NGS-based method and with a comparison method based on genotyping using type-specific probes bound to fluorescent beads (Luminex). Thirty-three clinical samples were repeat tested using NGS to evaluate reproducibility. RESULTS The NGS-based method correctly identified all 15 mixed HPV types when present in 100 copies/sample and 13/15 types when present in 10 copies/sample. For 36/60 cervical samples genotyping results using NGS and Luminex were identical. For 12/60 samples the NGS method was more sensitive than the Luminex test and most of the remaining discrepancies could be explained by the different type coverage of the assays. Reproducibility testing found complete or partial concordance in 30/33 samples. CONCLUSIONS NGS provides a sensitive and accurate method for genotyping of HPV. The fact that also the amplimer sequence is obtained could be important for studying the epidemiology of viral variants and monitoring of HPV vaccination.

Diversity of human papillomaviruses in skin lesions

Pools of frozen biopsies from patients with squamous cell carcinoma (SCC) (n=29) actinic keratosis (AK) (n=31), keratoacanthoma (n=91) and swab samples from 84 SCCs and 91 AKs were analysed with an extended HPV general primer PCR and high-throughput sequencing of amplimers. We found 273 different HPV isolates (87 known HPV types, 139 previously known HPV sequences (putative types) and 47 sequences from novel putative HPV types). Among the new sequences, five clustered in genus Betapapillomavirus and 42 in genus Gammapapillomavirus. Resequencing of the three pools between 21 to 70 times resulted in the detection of 283 different known or putative HPV types, with 156 different sequences found in only one of the pools. Type-specific PCRs for 37 putative types from an additional 296 patients found only two of these putative types. In conclusion, skin lesions contain a large diversity of HPV types, but most appeared to be rare infections.

Phylogenetically diverse TT virus viremia among pregnant women

Infections during pregnancy have been suggested to be involved in childhood leukemias. We used high-throughput sequencing to describe the viruses most readily detectable in serum samples of pregnant women. Serum DNA of 112 mothers to leukemic children was amplified using whole genome amplification. Sequencing identified one TT virus (TTV) isolate belonging to a known type and two putatively new TTVs. For 22 mothers, we also performed TTV amplification by general primer PCR before sequencing. This detected 39 TTVs, two of which were identical to the TTVs found after whole genome amplification. Altogether, we found 40 TTV isolates, 29 of which were putatively new types (similarities ranging from 89% to 69%). In conclusion, high throughput sequencing is useful to describe the known or unknown viruses that are present in serum samples of pregnant women.

Detection of DNA viruses in prostate cancer

We tested prostatic secretions from men with and without prostate cancer (13 cases and 13 matched controls) or prostatitis (18 cases and 18 matched controls) with metagenomic sequencing. A large number (>200) of viral reads was only detected among four prostate cancer cases (1 patient each positive for Merkel cell polyomavirus, JC polyomavirus and Human Papillomavirus types 89 or 40, respectively). Lower numbers of reads from a large variety of viruses were detected in all patient groups. Our knowledge of the biology of the prostate may be furthered by the fact that DNA viruses are commonly shed from the prostate and can be readily detected by metagenomic sequencing of expressed prostate secretions.

Viremia during pregnancy and risk of childhood leukemia and lymphomas in the offspring: Nested case-control study

A possible role for infections of the pregnant mother in the development of childhood acute leukemias and lymphomas has been suggested. However, no specific infectious agent has been identified. Offspring of 74,000 mothers who had serum samples taken during pregnancy and stored in a large‐scale biobank were followed up to the age of 15 years (750,000 person years) through over‐generation linkages between the biobank files, the Swedish national population and cancer registers to identify incident leukemia/lymphoma cases in the offspring. First‐trimester sera from mothers of 47 cases and 47 matched controls were retrieved and analyzed using next generation sequencing. Anelloviruses were the most common viruses detected, found in 37/47 cases and in 40/47 controls, respectively (OR: 0.6, 95% CI: 0.2–1.9). None of the detected viruses was associated with leukemia/lymphoma in the offspring. Viremia during pregnancy was common, but no association with leukemia/lymphoma risk in the offspring was found.

Are human papillomavirus DNA prevalences providing high-flying estimates of infection? An international survey of HPV detection on environmental surfaces

Most epidemiological studies of human papillomavirus (HPV) infections rely on HPV DNA detection. Recent studies have reported very high prevalence and acquisition rates in men.1 However, presence of HPV DNA is not proof of infection, as it might represent environmental contamination. Studies of HPV DNA on environmental surfaces could provide insights of the possible magnitude of this problem. We studied surfaces that frequently contact anogenital skin: toilet seats in airport restrooms. Apparently clean seats in 23 airports located in 13 countries (Belgium, Denmark, Finland, France, Germany, Italy, Jordan, the Netherlands, Russia, …

Extension of the viral ecology in humans using viral profile hidden Markov models

When human samples are sequenced, many assembled contigs are “unknown”, as conventional alignments find no similarity to known sequences. Hidden Markov models (HMM) exploit the positions of specific nucleotides in protein-encoding codons in various microbes. The algorithm HMMER3 implements HMM using a reference set of sequences encoding viral proteins, “vFam”. We used HMMER3 analysis of “unknown” human sample-derived sequences and identified 510 contigs distantly related to viruses (Anelloviridae (n = 1), Baculoviridae (n = 34), Circoviridae (n = 35), Caulimoviridae (n = 3), Closteroviridae (n = 5), Geminiviridae (n = 21), Herpesviridae (n = 10), Iridoviridae (n = 12), Marseillevirus (n = 26), Mimiviridae (n = 80), Phycodnaviridae (n = 165), Poxviridae (n = 23), Retroviridae (n = 6) and 89 contigs related to described viruses not yet assigned to any taxonomic family). In summary, we find that analysis using the HMMER3 algorithm and the “vFam” database greatly extended the detection of viruses in biospecimens from humans.

Viruses in cancers among the immunosuppressed

Most cancer forms known to be caused by viruses are increased among the immunosuppressed, but several cancer forms without established viral etiology are also increased, notably nonmelanoma skin carcinoma (NMSC). We followed all 13,429 solid organ transplantation patients in Sweden for cancer occurrence after transplantation. We requested these tumor specimens and sequenced the first 89 specimens received (62 NMSCs, 27 other cancers). The sequences were analyzed for viruses based on two bioinformatics algorithms (paracel‐blast (sensitive for detection of known viruses) and hidden Markov model (HMM; sensitive for distantly related viruses)). Among the 62 NMSCs, the virus family detected in the largest proportion of specimens was Mimiviridae (9/62 NMSCs). The majority of the virus‐related reads belonged to Papillomaviridae. The HMM analysis identified 86 additional previously not described viral contigs related to 11 virus families, with reads related to Mimiviridae being the most common (detected in 28/62 NMSCs) with the most prevalent contig (Mimivirus SE906, 1937 bp) detected in 17/62 NMSCs. Among the 27 other cancers, viral sequences were detected in only 5 specimens by blast analysis, compared to in all 27 specimens by HMM (Mimiviridae, Poxviridae, Phycodnaviridae and virus‐related sequences yet unclassified to any family). 99% of the virus reads belonged to a single previously not described sequence (Mimivirus SE996, 911 bp). A multitude of viruses is readily detectable in specimens with cancers occurring among the immunosuppressed, with sequences related to Mimiviridae being the most prevalent. Further research would be needed to elucidate the biological significance of the viruses.