VirusTAP: Viral Genome-Targeted Assembly Pipeline

Intra-Patient Genomic Variations of Human Papillomavirus Type 31 in Cervical Cancer and Precancer

Human papillomavirus type 31 (HPV31) is detected less frequently in cervical cancer than two major causative types, HPV16 and HPV18. Here, we report a comprehensive analysis of HPV31 genome sequences in cervical lesions collected from Japanese women. Of 52 HPV31-positive cervical specimens analyzed by deep sequencing, 43 samples yielded complete genome sequences of around 7900 base pairs and 9 samples yielded partially deleted genome sequences. Phylogenetic analysis showed that HPV31 variant distribution was lineage A in 19 samples (36.5%), lineage B in 28 samples (53.8%), and lineage C in 5 samples (9.6%), indicating that lineage B variants are dominant among HPV31 infections in Japan. Deletions in the viral genome were found in the region from the E1 to L1 genes, but all the deleted genomes retained the E6/E7 genes. Among intra-patient nucleotide variations relative to a consensus genome sequence in each sample, C-to-T substitutions were most frequently detected, followed by T-to-C and C-to-A substitutions. High-frequency, intra-patient mutations (>10%) in cervical cancer samples were found in the E1, E2, and E7 genes, and all of them were nonsynonymous substitutions. The enrichment of high-frequency nonsynonymous substitutions strongly suggests that these intra-patient mutations are positively selected during the development of cervical cancer/precancer.

Pathogen Profiles in Outpatients with Non-COVID-19 during the 7th Prevalent Period of COVID-19 in Gunma, Japan

The identification of pathogens associated with respiratory symptoms other than the novel coronavirus disease 2019 (COVID-19) can be challenging. However, the diagnosis of pathogens is crucial for assessing the clinical outcome of patients. We comprehensively profiled pathogens causing non-COVID-19 respiratory symptoms during the 7th prevalent period in Gunma, Japan, using deep sequencing combined with a next-generation sequencer (NGS) and advanced bioinformatics technologies. The study included nasopharyngeal swabs from 40 patients who tested negative for severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) using immuno-chromatography and/or quantitative reverse transcription polymerase chain reaction (qRT-PCR) methods. Comprehensive pathogen sequencing was conducted through deep sequencing using NGS. Additionally, short reads obtained from NGS were analyzed for comprehensive pathogen estimation using MePIC (Metagenomic Pathogen Identification Pipeline for Clinical Specimens) and/or VirusTap. The results revealed the presence of various pathogens, including respiratory viruses and bacteria, in the present subjects. Notably, human adenovirus (HAdV) was the most frequently detected virus in 16 of the 40 cases (40.0%), followed by coryneforms, which were the most frequently detected bacteria in 21 of the 40 cases (52.5%). Seasonal human coronaviruses (NL63 type, 229E type, HKU1 type, and OC43 type), human bocaviruses, and human herpesviruses (human herpesvirus types 1-7) were not detected. Moreover, multiple pathogens were detected in 50% of the subjects. These results suggest that various respiratory pathogens may be associated with non-COVID-19 patients during the 7th prevalent period in Gunma Prefecture, Japan. Consequently, for an accurate diagnosis of pathogens causing respiratory infections, detailed pathogen analyses may be necessary. Furthermore, it is possible that various pathogens, excluding SARS-CoV-2, may be linked to fever and/or respiratory infections even during the COVID-19 pandemic.

Unveiling Mycoviromes Using Fungal Transcriptomes

Viruses infecting fungi are referred to as mycoviruses. Here, we carried out in silico mycovirome studies using public fungal transcriptomes mostly derived from mRNA libraries. We identified 468 virus-associated contigs assigned to 5 orders, 21 families, 26 genera, and 88 species. We assembled 120 viral genomes with diverse RNA and DNA genomes. The phylogenetic tree and genome organization unveiled the possible host origin of mycovirus species and diversity of their genome structures. Most identified mycoviruses originated from fungi; however, some mycoviruses had strong phylogenetic relationships with those from insects and plants. The viral abundance and mutation frequency of mycoviruses were very low; however, the compositions and populations of mycoviruses were very complex. Although coinfection of diverse mycoviruses in the fungi was common in our study, most mycoviromes had a dominant virus species. The compositions and populations of mycoviruses were more complex than we expected. Viromes of Monilinia species revealed that there were strong deviations in the composition of viruses and viral abundance among samples. Viromes of Gigaspora species showed that the chemical strigolactone might promote virus replication and mutations, while symbiosis with endobacteria might suppress virus replication and mutations. This study revealed the diversity and host distribution of mycoviruses.

Rapid genotyping of targeted viral samples using Illumina short-read sequencing data

The most important information about microorganisms might be their accurate genome sequence. Using current Next Generation Sequencing methods, sequencing data can be generated at an unprecedented pace. However, we still lack tools for the automated and accurate reference-based genotyping of viral sequencing reads. This paper presents our pipeline designed to reconstruct the dominant consensus genome of viral samples and analyze their within-host variability. We benchmarked our approach on numerous datasets and showed that the consensus genome of samples could be obtained reliably without further manual data curation. Our pipeline can be a valuable tool for fast identifying viral samples. The pipeline is publicly available on the project’s GitHub page (https://github.com/laczkol/QVG).

VGEA: an RNA viral assembly toolkit

Next generation sequencing (NGS)-based studies have vastly increased our understanding of viral diversity. Viral sequence data obtained from NGS experiments are a rich source of information, these data can be used to study their epidemiology, evolution, transmission patterns, and can also inform drug and vaccine design. Viral genomes, however, represent a great challenge to bioinformatics due to their high mutation rate and forming quasispecies in the same infected host, bringing about the need to implement advanced bioinformatics tools to assemble consensus genomes well-representative of the viral population circulating in individual patients. Many tools have been developed to preprocess sequencing reads, carry-out de novo or reference-assisted assembly of viral genomes and assess the quality of the genomes obtained. Most of these tools however exist as standalone workflows and usually require huge computational resources. Here we present (Viral Genomes Easily Analyzed), a Snakemake workflow for analyzing RNA viral genomes. VGEA enables users to map sequencing reads to the human genome to remove human contaminants, split bam files into forward and reverse reads, carry out de novo assembly of forward and reverse reads to generate contigs, pre-process reads for quality and contamination, map reads to a reference tailored to the sample using corrected contigs supplemented by the user's choice of reference sequences and evaluate/compare genome assemblies. We designed a project with the aim of creating a flexible, easy-to-use and all-in-one pipeline from existing/stand-alone bioinformatics tools for viral genome analysis that can be deployed on a personal computer. VGEA was built on the Snakemake workflow management system and utilizes existing tools for each step: fastp (Chen et al., 2018) for read trimming and read-level quality control, BWA (Li & Durbin, 2009) for mapping sequencing reads to the human reference genome, SAMtools (Li et al., 2009) for extracting unmapped reads and also for splitting bam files into fastq files, IVA (Hunt et al., 2015) for de novo assembly to generate contigs, shiver (Wymant et al., 2018) to pre-process reads for quality and contamination, then map to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences, SeqKit (Shen et al., 2016) for cleaning shiver assembly for QUAST, QUAST (Gurevich et al., 2013) to evaluate/assess the quality of genome assemblies and MultiQC (Ewels et al., 2016) for aggregation of the results from fastp, BWA and QUAST. Our pipeline was successfully tested and validated with SARS-CoV-2 (n = 20), HIV-1 (n = 20) and Lassa Virus (n = 20) datasets all of which have been made publicly available. VGEA is freely available on GitHub at: https://github.com/pauloluniyi/VGEA under the GNU General Public License.

Intertypic reassortment of mammalian orthoreovirus identified in wastewater in Japan

Mammalian orthoreovirus (MRV), a non-enveloped virus with a ten-segmented double-stranded RNA genome, infects virtually all mammals, including humans. Human infection with MRV seems to be common in early childhood, but is rarely symptomatic. Despite the ubiquitous presence of MRV in mammals as well as in environmental waters, the molecular characterisation of the MRV genome remains to be fully elucidated. In this study, two novel strains, MRV-2 THK0325 and MRV-1 THK0617, were unintentionally isolated from wastewater in Japan via an environmental surveillance of enteric viruses. Homology and phylogenetic analysis demonstrated that all the segments of THK0325 were closely related to the MRV-2 Osaka strains, which were recently proposed to have existed for at least two decades in Japan. Most of the segments in THK0617 also showed a close relationship with the MRV-2 Osaka strains, but the M2, S1, and S3 segments belong to another MRV cluster. According to the S1 sequence, the determinant of serotype THK0617 was classified as MRV-1, and both the M2 and S3 segments were closely related to MRV-1 and -3 from the tree shrew in China. These results suggest that the MRV-2 Osaka-like strain spread widely throughout Japan, accompanied by intertypic reassortment occurring in East Asia.

High Levels of Within-Host Variations of Human Papillomavirus 16 E1/E2 Genes in Invasive Cervical Cancer

Human papillomavirus type 16 (HPV16) is the most common HPV genotype found in invasive cervical cancer (ICC). Recent comprehensive genomics studies of HPV16 have revealed that a large number of minor nucleotide variations in the viral genome are present in each infected woman; however, it remains unclear whether such within-host variations of HPV16 are linked to cervical carcinogenesis. Here, by employing next-generation sequencing approaches, we explored the mutational profiles of the HPV16 genome within individual clinical specimens from ICC (n = 31) and normal cervix (n = 21) in greater detail. A total of 367 minor nucleotide variations (167 from ICC and 200 from the normal cervix) were detected throughout the viral genome in both groups, while nucleotide variations at high frequencies (>10% abundance in relative read counts in a single sample) were more prevalent in ICC (10 in ICC versus 1 in normal). Among the high-level variations found in ICC, six were located in the E1/E2 genes, and all of them were non-synonymous substitutions (Q142K, M207I, and L262V for E1; D153Y, R302T, and T357A for E2). In vitro functional analyses of these E1/E2 variants revealed that E1/M207I, E2/D153Y, and E2/R302T had reduced abilities to support viral replication, and that E2/D153Y and E2/R302T failed to suppress the viral early promoter. These results imply that some within-host variations of E1/E2 present at high levels in ICC may be positively selected for and contribute to cervical cancer development through dysfunction or de-stabilization of viral replication/transcription proteins.

Detection of severe fever with thrombocytopenia syndrome virus and other viruses in cats via unbiased next-generation sequencing

We used unbiased next-generation sequencing (NGS) to detect unknown viruses in cats. Serum or plasma samples were obtained from clinically ill cats with suspected acute viral infections. Nucleic acid was extracted from serum or plasma samples to construct a complementary DNA library for NGS. Comprehensive nucleotide sequencing analyses enabled detection of the genomes of various viruses, including the severe fever with thrombocytopenia syndrome virus, feline immunodeficiency virus, feline morbillivirus, parvovirus, and Torque teno felis virus. Our findings indicate that comprehensive nucleotide analyses of serum or plasma samples can be used to detect infections with unknown viruses in cats.

Comprehensive pathogen detection in sera of Kawasaki disease patients by high-throughput sequencing: a retrospective exploratory study

Background

Kawasaki disease (KD) is an idiopathic systemic vasculitis that predominantly damages coronary arteries in children. Various pathogens have been investigated as triggers for KD, but no definitive causative pathogen has been determined. As KD is diagnosed by symptoms, several days are needed for diagnosis. Therefore, at the time of diagnosis of KD, the pathogen of the trigger may already be diminished. The aim of this study was to explore comprehensive pathogens in the sera at the acute stage of KD using high-throughput sequencing (HTS).

Methods

Sera of 12 patients at an extremely early stage of KD and 12 controls were investigated. DNA and RNA sequences were read separately using HTS. Sequence data were imported into the home-brew meta-genomic analysis pipeline, PATHDET, to identify the pathogen sequences.

Results

No RNA virus reads were detected in any KD case except for that of equine infectious anemia, which is known as a contaminant of commercial reverse transcriptase. Concerning DNA viruses, human herpesvirus 6B (HHV-6B, two cases) and Anelloviridae (eight cases) were detected among KD cases as well as controls. Multiple bacterial reads were obtained from KD and controls. Bacteria of the genera Acinetobacter, Pseudomonas, Delfita, Roseomonas, and Rhodocyclaceae appeared to be more common in KD sera than in the controls.

Conclusion

No single pathogen was identified in serum samples of patients at the acute phase of KD. With multiple bacteria detected in the serum samples, it is difficult to exclude the possibility of contamination; however, it is possible that these bacteria might stimulate the immune system and induce KD.

High-Throughput Sequencing for Deciphering the Virome of Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.), also known as lucerne, is a major forage crop worldwide. In the United States, it has recently become the third most valuable field crop, with an estimated value of over $9.3 billion. Alfalfa is naturally infected by many different pathogens, including viruses, obligate parasites that reproduce only inside living host cells. Traditionally, viral infections of alfalfa have been considered by breeders, growers, producers and researchers to be diseases of limited importance, although they are widespread in all major cultivation areas. However, over the past few years, due to the rapid development of high-throughput sequencing (HTS), viral metagenomics, bioinformatics tools for interpreting massive amounts of HTS data and the increasing accessibility of public data repositories for transcriptomic discoveries, several emerging viruses of alfalfa with the potential to cause serious yield losses have been described. They include alfalfa leaf curl virus (family Geminiviridae), alfalfa dwarf virus (family Rhabdoviridae), alfalfa enamovirus 1 (family Luteoviridae), alfalfa virus S (family Alphaflexiviridae) and others. These discoveries have called into question the assumed low economic impact of viral diseases in alfalfa and further suggested their possible contribution to the severity of complex infections involving multiple pathogens. In this review, we will focus on viruses of alfalfa recently described in different laboratories on the basis of the above research methodologies.

HPV-EM: an accurate HPV detection and genotyping EM algorithm

Accurate HPV genotyping is crucial in facilitating epidemiology studies, vaccine trials, and HPV-related cancer research. Contemporary HPV genotyping assays only detect < 25% of all known HPV genotypes and are not accurate for low-risk or mixed HPV genotypes. Current genomic HPV genotyping algorithms use a simple read-alignment and filtering strategy that has difficulty handling repeats and homology sequences. Therefore, we have developed an optimized expectation–maximization algorithm, designated HPV-EM, to address the ambiguities caused by repetitive sequencing reads. HPV-EM achieved 97–100% accuracy when benchmarked using cell line data and TCGA cervical cancer data. We also validated HPV-EM using DNA tiling data on an institutional cervical cancer cohort (96.5% accuracy). Using HPV-EM, we demonstrated HPV genotypic differences in recurrence and patient outcomes in cervical and head and neck cancers.

Molecular epidemiology of human papillomavirus 18 infections in Japanese Women

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

Transcriptome analysis of HPV-induced warts and healthy skin in humans

Background

The human papillomaviruses (HPV) are a group of viruses that, depending on the strain, can cause cancer or the formation of benign growths known as warts. Scarce information exists with regard to the genetic nature of non-genital cutaneous warts induced by the human papillomavirus (HPV).

Methods

The main purpose of this study is to investigate the differences between the gene expression profiles of common warts and healthy skin in HPV-positive individuals by RNA sequencing on the Illumina HiSeq 2500. After obtaining shave biopsies of common warts and healthy skin from twelve Arab males, we were able to analyze the transcriptomes of 24 paired cases and controls.

Results

Common warts were found to possess a highly significant and unique molecular signature. Many of the most up-regulated (KRT16, EPGN, and ABCG4) and down-regulated genes (C15orf59, CYB561A3, and FCGRT) in warts were the subject of little investigation in the published literature. Moreover, the top 500 differentially expressed genes were found to be associated with immune and autoimmune pathways, such as the neutrophil degranulation, toll-like receptor 7/8 (TLR 7/8) cascade, toll-like receptor 9 (TLR9) cascade, and toll-like receptor 10 (TLR10) pathways, among others.

Conclusions

Our findings are particularly important because they serve as the most comprehensive to date with regard to the modulation of human skin gene expression by HPV infection.

DEN-IM: dengue virus genotyping from amplicon and shotgun metagenomic sequencing

Dengue virus (DENV) represents a public health threat and economic burden in affected countries. The availability of genomic data is key to understanding viral evolution and dynamics, supporting improved control strategies. Currently, the use of high-throughput sequencing (HTS) technologies, which can be applied both directly to patient samples (shotgun metagenomics) and to PCR-amplified viral sequences (amplicon sequencing), is potentially the most informative approach to monitor viral dissemination and genetic diversity by providing, in a single methodological step, identification and characterization of the whole viral genome at the nucleotide level. Despite many advantages, these technologies require bioinformatics expertise and appropriate infrastructure for the analysis and interpretation of the resulting data. In addition, the many software solutions available can hamper the reproducibility and comparison of results. Here we present DEN-IM, a one-stop, user-friendly, containerized and reproducible workflow for the analysis of DENV short-read sequencing data from both amplicon and shotgun metagenomics approaches. It is able to infer the DENV coding sequence (CDS), identify the serotype and genotype, and generate a phylogenetic tree. It can easily be run on any UNIX-like system, from local machines to high-performance computing clusters, performing a comprehensive analysis without the requirement for extensive bioinformatics expertise. Using DEN-IM, we successfully analysed two types of DENV datasets. The first comprised 25 shotgun metagenomic sequencing samples from patients with variable serotypes and genotypes, including an in vitro spiked sample containing the four known serotypes. The second consisted of 106 paired-end and 76 single-end amplicon sequences of DENV 3 genotype III and DENV 1 genotype I, respectively, where DEN-IM allowed detection of the intra-genotype diversity. The DEN-IM workflow, parameters and execution configuration files, and documentation are freely available at https://github.com/B-UMMI/DEN-IM).

Complete Genome Sequencing of Influenza A Viruses Using Next-Generation Sequencing

Recently, chain termination sequencing methods have been replaced by more efficient next-generation sequencing (NGS) methods. For influenza A, NGS allows for deep sequencing to characterize virus populations, efficient complete genome sequencing, and a non-sequence-dependent method to identify viral variants. There are numerous approaches to preparing samples for NGS and subsequent data processing methods that can be applied to influenza A sequencing. This chapter provides a brief overview of the process of NGS for influenza A and some useful bioinformatics tools for developing an NGS workflow for influenza A viruses.

Next generation sequencing and proteomics in plant virology: how is Colombia doing?

Crop production and trade are two of the most economically important activities in Colombia, and viral diseases cause a high negative impact to agricultural sector. Therefore, the detection, diagnosis, control, and management of viral diseases are crucial. Currently, Next-Generation Sequencing (NGS) and ‘Omic’ technologies constitute a right-hand tool for the discovery of novel viruses and for studying virus-plant interactions. This knowledge allows the development of new viral diagnostic methods and the discovery of key components of infectious processes, which could be used to generate plants resistant to viral infections. Globally, crop sciences are advancing in this direction. In this review, advancements in ‘omic’ technologies and their different applications in plant virology in Colombia are discussed. In addition, bioinformatics pipelines and resources for omics data analyses are presented. Due to their decreasing prices, NGS technologies are becoming an affordable and promising means to explore many phytopathologies affecting a wide variety of Colombian crops so as to improve their trade potential.

HPViewer: sensitive and specific genotyping of human papillomavirus in metagenomic DNA

Motivation

Shotgun DNA sequencing provides sensitive detection of all 182 HPV types in tissue and body fluid. However, existing computational methods either produce false positives misidentifying HPV types due to shared sequences among HPV, human and prokaryotes, or produce false negative since they identify HPV by assembled contigs requiring large abundant of HPV reads.

Results

We designed HPViewer with two custom HPV reference databases masking simple repeats and homology sequences respectively and one homology distance matrix to hybridize these two databases. It directly identified HPV from short DNA reads rather than assembled contigs. Using 100 100 simulated samples, we revealed that HPViewer was robust for samples containing either high or low number of HPV reads. Using 12 shotgun sequencing samples from respiratory papillomatosis, HPViewer was equal to VirusTAP, and Vipie and better than HPVDetector with the respect to specificity and was the most sensitive method in the detection of HPV types 6 and 11. We demonstrated that contigs-based approaches had disadvantages of detection of HPV. In 1573 sets of metagenomic data from 18 human body sites, HPViewer identified 104 types of HPV in a body-site associated pattern and 89 types of HPV co-occurring in one sample with other types of HPV. We demonstrated HPViewer was sensitive and specific for HPV detection in metagenomic data.

Availability and implementation

HPViewer can be accessed at https://github.com/yuhanH/HPViewer/.

Supplementary information

Supplementary data are available at Bioinformatics online.

TAR-VIR: a pipeline for TARgeted VIRal strain reconstruction from metagenomic data

Background

Strain-level RNA virus characterization is essential for developing prevention and treatment strategies. Viral metagenomic data, which can contain sequences of both known and novel viruses, provide new opportunities for characterizing RNA viruses. Although there are a number of pipelines for analyzing viruses in metagenomic data, they have different limitations. First, viruses that lack closely related reference genomes cannot be detected with high sensitivity. Second, strain-level analysis is usually missing.

Results

In this study, we developed a hybrid pipeline named TAR-VIR that reconstructs viral strains without relying on complete or high-quality reference genomes. It is optimized for identifying RNA viruses from metagenomic data by combining an effective read classification method and our in-house strain-level de novo assembly tool. TAR-VIR was tested on both simulated and real viral metagenomic data sets. The results demonstrated that TAR-VIR competes favorably with other tested tools.

Conclusion

TAR-VIR can be used standalone for viral strain reconstruction from metagenomic data. Or, its read recruiting stage can be used with other de novo assembly tools for superior viral functional and taxonomic analyses. The source code and the documentation of TAR-VIR are available at https://github.com/chjiao/TAR-VIR.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2878-2) contains supplementary material, which is available to authorized users.

High Throughput Sequencing For Plant Virus Detection and Discovery

Over the last decade, virologists have discovered an unprecedented number of viruses using high throughput sequencing (HTS), which led to the advancement of our knowledge on the diversity of viruses in nature, particularly unraveling the virome of many agricultural crops. However, these new virus discoveries have often widened the gaps in our understanding of virus biology; the forefront of which is the actual role of a new virus in disease, if any. Yet, when used critically in etiological studies, HTS is a powerful tool to establish disease causality between the virus and its host. Conversely, with globalization, movement of plant material is increasingly more common and often a point of dispute between countries. HTS could potentially resolve these issues given its capacity to detect and discover. Although many pipelines are available for plant virus discovery, all share a common backbone. A description of the process of plant virus detection and discovery from HTS data are presented, providing a summary of the different pipelines available for scientists' utility in their research.

Whole-Genome Analysis of Human Papillomavirus Type 16 Prevalent in Japanese Women with or without Cervical Lesions

Recent large-scale genomics studies of human papillomaviruses (HPVs) have shown a high level of genomic variability of HPV16, the most prevalent genotype in HPV-associated malignancies, and provided new insights into the biological and clinical relevance of its genetic variations in cervical cancer development. Here, we performed deep sequencing analyses of the viral genome to explore genetic variations of HPV16 that are prevalent in Japan. A total of 100 complete genome sequences of HPV16 were determined from cervical specimens collected from Japanese women with cervical intraepithelial neoplasia and invasive cervical cancer, or without cervical malignancies. Phylogenetic analyses revealed the variant distribution in the Japanese HPV16 isolates; overall, lineage A was the most prevalent (94.0%), in which sublineage A4 was dominant (52.0%), followed by sublineage A1 (21.0%). The relative risk of sublineage A4 for cervical cancer development was significantly higher compared to sublineages A1/A2/A3 (odds ratio = 6.72, 95% confidence interval = 1.78–28.9). Interestingly, a novel cluster of variants that branched from A1/A2/A3 was observed for the Japanese HPV16 isolates, indicating that unique HPV16 variants are prevalent among Japanese women.

virMine: automated detection of viral sequences from complex metagenomic samples

Metagenomics has enabled sequencing of viral communities from a myriad of different environments. Viral metagenomic studies routinely uncover sequences with no recognizable homology to known coding regions or genomes. Nevertheless, complete viral genomes have been constructed directly from complex community metagenomes, often through tedious manual curation. To address this, we developed the software tool virMine to identify viral genomes from raw reads representative of viral or mixed (viral and bacterial) communities. virMine automates sequence read quality control, assembly, and annotation. Researchers can easily refine their search for a specific study system and/or feature(s) of interest. In contrast to other viral genome detection tools that often rely on the recognition of viral signature sequences, virMine is not restricted by the insufficient representation of viral diversity in public data repositories. Rather, viral genomes are identified through an iterative approach, first omitting non-viral sequences. Thus, both relatives of previously characterized viruses and novel species can be detected, including both eukaryotic viruses and bacteriophages. Here we present virMine and its analysis of synthetic communities as well as metagenomic data sets from three distinctly different environments: the gut microbiota, the urinary microbiota, and freshwater viromes. Several new viral genomes were identified and annotated, thus contributing to our understanding of viral genetic diversity in these three environments.

Genome Detective: an automated system for virus identification from high-throughput sequencing data

SUMMARY

Genome Detective is an easy to use web-based software application that assembles the genomes of viruses quickly and accurately. The application uses a novel alignment method that constructs genomes by reference-based linking of de novo contigs by combining amino-acids and nucleotide scores. The software was optimized using synthetic datasets to represent the great diversity of virus genomes. The application was then validated with next generation sequencing data of hundreds of viruses. User time is minimal and it is limited to the time required to upload the data.

AVAILABILITY AND IMPLEMENTATION

Available online: http://www.genomedetective.com/app/typingtool/virus/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Near-Complete Genome Sequence of Infectious Bronchitis Virus Strain VFAR-047 (GI-16 Lineage), Isolated in Peru

Here, we report the near-complete genome sequence of the infectious bronchitis virus (IBV) strain VFAR-047, isolated in Peru in 2014. This strain was classified into GI lineage 16 (GI-16) based on both the genome and Spike 1 (S1) sequence analysis.

Here, we report the near-complete genome sequence of the infectious bronchitis virus (IBV) strain VFAR-047, isolated in Peru in 2014. This strain was classified into GI lineage 16 (GI-16) based on both the genome and Spike 1 (S1) sequence analysis. Furthermore, four potential recombination events with other GI-16 and GI-11 strains were identified.

Detection by Direct Next Generation Sequencing Analysis of Emerging Enterovirus D68 and C109 Strains in an Environmental Sample From Scotland

Background: Human enteroviruses (EVs) have been linked with severe disease and syndromes as varied as acute respiratory illness, myocarditis, and flaccid paralysis. With global polio eradication on sight the focus of clinical investigations has expanded to the identification of other EV serotypes associated with severe neurological conditions such as EV-D68, responsible for large outbreaks in 2014 and 2016 that spread worldwide and were related with severe respiratory disease leading to acute myelitis in some cases. New EV serotypes with epidemic potential continue to emerge such as EV-C104, EV-C105, EV-C109, and EV-C117 identified in respiratory samples in recent years.

Methods: We used a next generation sequencing (NGS) approach to detect multiple EV serotypes directly in a sewage concentrate from Glasgow (Scotland, United Kingdom) generating whole-capsid nucleotide sequences that were compared to sequences of cell culture isolates from this sewage sample and clinical EV isolates from GenBank.

Results: Thirteen different serotypes belonging to all four A, B, C, and D EV species were identified in the sewage concentrate. EV strains closely related to EV-D68 epidemic isolates of B3 lineage reported in the United States and Europe in 2016 and to EV-C109 respiratory isolates found in Denmark and Netherlands in 2015 were identified.

Conclusion: Environmental surveillance (ES) can effectively detect EV circulation in human populations. The use of NGS for ES can help overcoming the limitations of traditional cell culture and sequencing methods, which are selective and biased, and can contribute to the early detection and assessment of spread of emerging EV pathogens.

Identification of potential pathogenic viruses in patients with acute myocarditis using next-generation sequencing

Myocarditis is an inflammatory disease of the myocardium and leads to cardiac dysfunction and heart failure. Although viral infections are considered to be the most common etiology of myocarditis, the identification of the causative virus is still challenging. Recently, next-generation sequencing (NGS) has been applied in the diagnosis of infectious diseases. The aim of the current study was to comprehensively analyze potential pathogenic microorganisms using NGS in the sera of patients with myocarditis. Twelve pediatric and five adult patients hospitalized for acute myocarditis were included. Serum samples in the acute phase were obtained and analyzed using NGS to detect pathogen-derived DNA and RNA. Viral sequence reads were detected in 7 (41%) of the 17 myocarditis patients by NGS. Among these patients, detection of Epstein-Barr virus, human parvovirus B19, torque teno virus, and respiratory syncytial virus reads by NGS was consistent with polymerase chain reaction or antigen test results in one patient each. A large number of human pegivirus reads were detected from one patient by RNA sequencing; however, its pathogenicity to human is unknown. Conversely, the number of detected virus-derived reads was small in most cases, and the pathophysiological role of these viruses remains to be clarified. No significant bacterial or fungal reads other than normal bacterial flora was detected. These data indicate that comprehensive detection of virus-derived DNA and RNA using NGS can be useful for the identification of potential pathogenic viruses in myocarditis.

Maximal viral information recovery from sequence data using VirMAP

Accurate classification of the human virome is critical to a full understanding of the role viruses play in health and disease. This implies the need for sensitive, specific, and practical pipelines that return precise outputs while still enabling case-specific post hoc analysis. Viral taxonomic characterization from metagenomic data suffers from high background noise and signal crosstalk that confounds current methods. Here we develop VirMAP that overcomes these limitations using techniques that merge nucleotide and protein information to taxonomically classify viral reconstructions independent of genome coverage or read overlap. We validate VirMAP using published data sets and viral mock communities containing RNA and DNA viruses and bacteriophages. VirMAP offers opportunities to enhance metagenomic studies seeking to define virome-host interactions, improve biosurveillance capabilities, and strengthen molecular epidemiology reporting.

Within-Host Variations of Human Papillomavirus Reveal APOBEC Signature Mutagenesis in the Viral Genome

Persistent infection with oncogenic human papillomaviruses (HPVs) causes cervical cancer, accompanied by the accumulation of somatic mutations into the host genome. There are concomitant genetic changes in the HPV genome during viral infection; however, their relevance to cervical carcinogenesis is poorly understood. Here, we explored within-host genetic diversity of HPV by performing deep-sequencing analyses of viral whole-genome sequences in clinical specimens. The whole genomes of HPV types 16, 52, and 58 were amplified by type-specific PCR from total cellular DNA of cervical exfoliated cells collected from patients with cervical intraepithelial neoplasia (CIN) and invasive cervical cancer (ICC) and were deep sequenced. After constructing a reference viral genome sequence for each specimen, nucleotide positions showing changes with >0.5% frequencies compared to the reference sequence were determined for individual samples. In total, 1,052 positions of nucleotide variations were detected in HPV genomes from 151 samples (CIN1, n = 56; CIN2/3, n = 68; ICC, n = 27), with various numbers per sample. Overall, C-to-T and C-to-A substitutions were the dominant changes observed across all histological grades. While C-to-T transitions were predominantly detected in CIN1, their prevalence was decreased in CIN2/3 and fell below that of C-to-A transversions in ICC. Analysis of the trinucleotide context encompassing substituted bases revealed that TpCpN, a preferred target sequence for cellular APOBEC cytosine deaminases, was a primary site for C-to-T substitutions in the HPV genome. These results strongly imply that the APOBEC proteins are drivers of HPV genome mutation, particularly in CIN1 lesions.IMPORTANCE HPVs exhibit surprisingly high levels of genetic diversity, including a large repertoire of minor genomic variants in each viral genotype. Here, by conducting deep-sequencing analyses, we show for the first time a comprehensive snapshot of the within-host genetic diversity of high-risk HPVs during cervical carcinogenesis. Quasispecies harboring minor nucleotide variations in viral whole-genome sequences were extensively observed across different grades of CIN and cervical cancer. Among the within-host variations, C-to-T transitions, a characteristic change mediated by cellular APOBEC cytosine deaminases, were predominantly detected throughout the whole viral genome, most strikingly in low-grade CIN lesions. The results strongly suggest that within-host variations of the HPV genome are primarily generated through the interaction with host cell DNA-editing enzymes and that such within-host variability is an evolutionary source of the genetic diversity of HPVs.

Isolation and Complete Genome Sequencing of Zika Virus Imported from the Dominican Republic to Japan

Zika virus (ZIKV) infection has been documented within Central and South America, Asia, and Africa. Here we report the isolation of virus from a patient infected with ZIKV returning to Japan from the Dominican Republic. The ZIKV strain was imaged by electron microscopy and its complete genome sequence was analyzed. Phylogenetic analysis and molecular characterization revealed that the strain was of Asian lineage, and carried 2 unique mutations in its NS5 region. These mutations are characteristic of strains that originated in the Dominican Republic and the USA in 2016.

On the effective depth of viral sequence data

Genome sequence data are of great value in describing evolutionary processes in viral populations. However, in such studies, the extent to which data accurately describes the viral population is a matter of importance. Multiple factors may influence the accuracy of a dataset, including the quantity and nature of the sample collected, and the subsequent steps in viral processing. To investigate this phenomenon, we sequenced replica datasets spanning a range of viruses, and in which the point at which samples were split was different in each case, from a dataset in which independent samples were collected from a single patient to another in which all processing steps up to sequencing were applied to a single sample before splitting the sample and sequencing each replicate. We conclude that neither a high read depth nor a high template number in a sample guarantee the precision of a dataset. Measures of consistency calculated from within a single biological sample may also be insufficient; distortion of the composition of a population by the experimental procedure or genuine within-host diversity between samples may each affect the results. Where it is possible, data from replicate samples should be collected to validate the consistency of short-read sequence data.

Identification of vaccine-derived rotavirus strains in children with acute gastroenteritis in Japan, 2012-2015

Two live attenuated oral rotavirus vaccines, Rotarix and RotaTeq, have been introduced as voluntary vaccination in Japan since 2011 and 2012, respectively. Effectiveness of the vaccines has been confirmed, whereas concerns such as shedding of the vaccine strains and gastroenteritis cases caused by vaccine strains are not well assessed. We aimed to identify the vaccine strains in children with acute gastroenteritis (AGE) to investigate the prevalence of AGE caused by vaccination or horizontal transmission of vaccine strains. A total of 1,824 stool samples were collected from children with AGE at six outpatient clinics in 2012-2015. Among all, 372 group A rotavirus (RVA) positive samples were screened for vaccine components by real-time RT-PCR which were designed to differentiate vaccine strains from rotavirus wild-type strains with high specificity. For samples possessing both vaccine and wild-type strains, analyses by next-generation sequencing (NGS) were conducted to characterize viruses existed in the intestine. As a result, Rotarix-derived strains were identified in 6 of 372 (1.6%) RVA positive samples whereas no RotaTeq strain was detected. Among six samples, four possessed Rotarix-derived strains while two possessed both Rotarix-derived strains and wild-type strains. In addition, other pathogens such as norovirus, enterovirus and E.coli were detected in four samples. The contribution of these vaccine strains to each patient's symptoms was unclear as all of the cases were vaccinated 2-14 days before sample collection. Proportion of average coverage for each segmented gene by NGS strongly suggested the concurrent infection of the vaccine-derived strain and the wild-type strain rather than reassortment of these two strains in one sample. This is the first study to report the prevalence of vaccine-derived strains in patients with RVA AGE in Japan as 1.6% without evidence of horizontal transmission. The results emphasized the importance of continuous monitoring on vaccine strains and their clinical impacts on children.

Whole-genome analysis of human papillomavirus genotypes 52 and 58 isolated from Japanese women with cervical intraepithelial neoplasia and invasive cervical cancer

Background

Human papillomavirus genotypes 52 and 58 (HPV52/58) are frequently detected in patients with cervical intraepithelial neoplasia (CIN) and invasive cervical cancer (ICC) in East Asian countries including Japan. As with other HPV genotypes, HPV52/58 consist of multiple lineages of genetic variants harboring less than 10% differences between complete genome sequences of the same HPV genotype. However, site variations of nucleotide and amino acid sequences across the viral whole-genome have not been fully examined for HPV52/58. The aim of this study was to investigate genetic variations of HPV52/58 prevalent among Japanese women by analyzing the viral whole-genome sequences.

Methods

The entire genomic region of HPV52/58 was amplified by long-range PCR with total cellular DNA extracted from cervical exfoliated cells isolated from Japanese patients with CIN or ICC. The amplified DNA was subjected to next generation sequencing to determine the complete viral genome sequences. Phylogenetic analyses were performed with the whole-genome sequences to assign variant lineages/sublineages to the HPV52/58 isolates. The variability in amino acid sequences of viral proteins was assessed by calculating the Shannon entropy scores at individual amino acid positions of HPV proteins.

Results

Among 52 isolates of HPV52 (CIN1, n = 20; CIN2/3, n = 21; ICC, n = 11), 50 isolates belonged to lineage B (sublineage B2) and two isolates belonged to lineage A (sublineage A1). Among 48 isolates of HPV58 (CIN1, n = 21; CIN2/3, n = 19; ICC, n = 8), 47 isolates belonged to lineage A (sublineages A1/A2/A3) and one isolate belonged to lineage C. Single nucleotide polymorphisms specific for individual variant lineages were determined throughout the viral genome based on multiple sequence alignments of the Japanese HPV52/58 isolates and reference HPV52/58 genomes. Entropy analyses revealed that the E1 protein was relatively variable among the HPV52 isolates, whereas the E7, E4, and L2 proteins showed some variations among the HPV58 isolates.

Conclusions

Among the HPV52/58-positive specimens from Japanese women with CIN/ICC, the variant distributions were strongly biased toward lineage B for HPV52 and lineage A for HPV58 across histological categories. Different patterns of amino acid variations were observed in HPV52 and HPV58 across the viral whole-genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13027-017-0155-4) contains supplementary material, which is available to authorized users.

drVM: a new tool for efficient genome assembly of known eukaryotic viruses from metagenomes

Background: Virus discovery using high-throughput next-generation sequencing has become more commonplace. However, although analysis of deep next-generation sequencing data allows us to identity potential pathogens, the entire analytical procedure requires competency in the bioinformatics domain, which includes implementing proper software packages and preparing prerequisite databases. Simple and user-friendly bioinformatics pipelines are urgently required to obtain complete viral genome sequences from metagenomic data.

Results: This manuscript presents a pipeline, drVM (detect and reconstruct known viral genomes from metagenomes), for rapid viral read identification, genus-level read partition, read normalization, de novo assembly, sequence annotation, and coverage profiling. The first two procedures and sequence annotation rely on known viral genomes as a reference database. drVM was validated via the analysis of over 300 sequencing runs generated by Illumina and Ion Torrent platforms to provide complete viral genome assemblies for a variety of virus types including DNA viruses, RNA viruses, and retroviruses. drVM is available for free download at: https://sourceforge.net/projects/sb2nhri/files/drVM/ and is also assembled as a Docker container, an Amazon machine image, and a virtual machine to facilitate seamless deployment.

Conclusions: drVM was compared with other viral detection tools to demonstrate its merits in terms of viral genome completeness and reduced computation time. This substantiates the platform's potential to produce prompt and accurate viral genome sequences from clinical samples.

Vipie: web pipeline for parallel characterization of viral populations from multiple NGS samples

Background

Next generation sequencing (NGS) technology allows laboratories to investigate virome composition in clinical and environmental samples in a culture-independent way. There is a need for bioinformatic tools capable of parallel processing of virome sequencing data by exactly identical methods: this is especially important in studies of multifactorial diseases, or in parallel comparison of laboratory protocols.

Results

We have developed a web-based application allowing direct upload of sequences from multiple virome samples using custom parameters. The samples are then processed in parallel using an identical protocol, and can be easily reanalyzed. The pipeline performs de-novo assembly, taxonomic classification of viruses as well as sample analyses based on user-defined grouping categories. Tables of virus abundance are produced from cross-validation by remapping the sequencing reads to a union of all observed reference viruses. In addition, read sets and reports are created after processing unmapped reads against known human and bacterial ribosome references. Secured interactive results are dynamically plotted with population and diversity charts, clustered heatmaps and a sortable and searchable abundance table.

Conclusions

The Vipie web application is a unique tool for multi-sample metagenomic analysis of viral data, producing searchable hits tables, interactive population maps, alpha diversity measures and clustered heatmaps that are grouped in applicable custom sample categories. Known references such as human genome and bacterial ribosomal genes are optionally removed from unmapped (‘dark matter’) reads. Secured results are accessible and shareable on modern browsers. Vipie is a freely available web-based tool whose code is open source.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3721-7) contains supplementary material, which is available to authorized users.

Deep-Sequence Identification and Role in Virus Replication of a JC Virus Quasispecies in Patients with Progressive Multifocal Leukoencephalopathy

JC virus (JCV) is a DNA virus causing progressive multifocal leukoencephalopathy (PML) in immunodeficient patients. In the present study, 22 genetic quasispecies with more than 1.5% variant frequency were detected in JCV genomes from six clinical samples of PML by next-generation sequencing. A mutation from A to C at nucleotide (nt) 3495 in JCV Mad1 resulting in a V-to-G amino acid substitution at amino acid (aa) position 392 of the large T antigen (TAg) was identified in all six cases of PML at 3% to 19% variant frequencies. Transfection of JCV Mad1 DNA possessing the V392G substitution in TAg into IMR-32 and human embryonic kidney 293 (HEK293) cells resulted in dramatically decreased production of JCV-encoded proteins. The virus DNA copy number was also reduced in supernatants of the mutant virus-transfected cells. Transfection of the IMR-32 and HEK293 cells with a virus genome containing a revertant mutation recovered viral production and protein expression. Cotransfection with equal amounts of wild-type genome and mutated JCV genome did not reduce the expression of viral proteins or viral replication, suggesting that the mutation did not have any dominant-negative function. Finally, immunohistochemistry demonstrated that TAg was expressed in all six pathological samples in which the quasispecies were detected. In conclusion, the V392G amino acid substitution in TAg identified frequently in PML lesions has a function in suppressing JCV replication, but the frequency of the mutation was restricted and its role in PML lesions was limited.

IMPORTANCE

DNA viruses generally have lower mutation frequency than RNA viruses, and the detection of quasispecies in JCV has rarely been reported. In the present study, a next-generation sequencer identified a JCV quasispecies with an amino acid substitution in the T antigen in patients with PML. In vitro studies showed that the mutation strongly repressed the expression of JC viral proteins and reduced the viral replication. However, because the frequency of the mutation was low in each case, the total expression of virus proteins was sustained in vivo. Thus, JC virus replicates in PML lesions in the presence of a mutant virus which is able to repress virus replication.

Establishing and characterizing a new primary effusion lymphoma cell line harboring Kaposi's sarcoma-associated herpesvirus

BACKGROUND

Primary effusion lymphoma is a rare distinct large B-cell neoplasm that is associated with Kaposi's sarcoma-associated herpesvirus (KSHV) infection. Over recent years, 9 KSHV-positive/Epstein-Barr virus (EBV)-negative PEL cell lines have been established.

METHODS

Tumor cells were collected from the pleural effusion of a 49-year-old male with AIDS. Cells were grown in RPMI1640 culture medium supplemented with 10 % fetus bovine serum. Single cell cloning was performed successfully by a limiting dilution method in a 96-well plate. The cell line obtained was designated SPEL.

RESULTS

SPEL cells showed gourd-shaped morphology with a polarized nucleus, expressing CD38, CD138, and Blimp-1, but not B cell markers such as CD19 and CD20. Polymerase chain reaction analysis revealed that SPEL cells were positive for KSHV but negative for EBV. Tetradecanoylphorbol acetate induced expression of KSHV lytic proteins and the production of KSHV particles in SPEL cells. Subcutaneous inoculation of SPEL cells into severe combined immunodeficiency mice resulted in the formation of solid tumors. Next-generation sequencing revealed the 138 kbp genome sequence of KSHV in SPEL cells. Suberic bishydroxamate, a histone deacetylase inhibitor, induced the expression of KSHV-encoded lytic proteins and cell death in SPEL cells.

CONCLUSIONS

A new KSHV-positive and EBV-negative PEL cell line, SPEL was established. This cell line may contribute to furthering our understanding of the pathogenesis of PEL and KSHV infection.