A modified RNA-Seq approach for whole genome sequencing of RNA viruses from faecal and blood samples

Genetically distinct within-host subpopulations of hepatitis C virus persist after Direct-Acting Antiviral treatment failure

Analysis of viral genetic data has previously revealed distinct within-host population structures in both untreated and interferon-treated chronic hepatitis C virus (HCV) infections. While multiple subpopulations persisted during the infection, each subpopulation was observed only intermittently. However, it was unknown whether similar patterns were also present after Direct-Acting Antiviral (DAA) treatment, where viral populations were often assumed to go through narrow bottlenecks. Here we tested for the maintenance of population structure after DAA treatment failure, and whether there were different evolutionary rates along distinct lineages where they were observed. We analysed whole-genome next-generation sequencing data generated from a randomised study using DAAs (the BOSON study). We focused on samples collected from patients (N=84) who did not achieve sustained virological response (i.e., treatment failure) and had sequenced virus from multiple timepoints. Given the short-read nature of the data, we used a number of methods to identify distinct within-host lineages including tracking concordance in intra-host nucleotide variant (iSNV) frequencies, applying sequenced-based and tree-based clustering algorithms to sliding windows along the genome, and haplotype reconstruction. Distinct viral subpopulations were maintained among a high proportion of individuals post DAA treatment failure. Using maximum likelihood modelling and model comparison, we found an overdispersion of viral evolutionary rates among individuals, and significant differences in evolutionary rates between lineages within individuals. These results suggest the virus is compartmentalised within individuals, with the varying evolutionary rates due to different viral replication rates and/or different selection pressures. We endorse lineage awareness in future analyses of HCV evolution and infections to avoid conflating patterns from distinct lineages, and to recognise the likely existence of unsampled subpopulations.

Metagenomic Surveillance of Viral Gastroenteritis in a Public Health Setting

Norovirus is the primary cause of viral gastroenteritis (GE). To investigate norovirus epidemiology, there is a need for whole-genome sequencing and reference sets consisting of complete genomes. To investigate the potential of shotgun metagenomic sequencing on the Illumina platform for whole-genome sequencing, 71 reverse transcriptase quantitative PCR (RT-qPCR) norovirus positive-feces (threshold cycle [CT], <30) samples from norovirus surveillance within The Netherlands were subjected to metagenomic sequencing. Data were analyzed through an in-house next-generation sequencing (NGS) analysis workflow. Additionally, we assessed the potential of metagenomic sequencing for the surveillance of off-target viruses that are of importance for public health, e.g., sapovirus, rotavirus A, enterovirus, parechovirus, aichivirus, adenovirus, and bocaparvovirus. A total of 60 complete and 10 partial norovirus genomes were generated, representing 7 genogroup I capsid genotypes and 12 genogroup II capsid genotypes. In addition to the norovirus genomes, the metagenomic approach yielded partial or complete genomes of other viruses for 39% of samples from children and 6.7% of samples from adults, including adenovirus 41 (N = 1); aichivirus 1 (N = 1); coxsackievirus A2 (N = 2), A4 (N = 2), A5 (N = 1), and A16 (N = 1); bocaparvovirus 1 (N = 1) and 3 (N = 1); human parechovirus 1 (N = 2) and 3 (N = 1); Rotavirus A (N = 1); and a sapovirus GI.7 (N = 1). The sapovirus GI.7 was initially not detected through RT-qPCR and warranted an update of the primer and probe set. Metagenomic sequencing on the Illumina platform robustly determines complete norovirus genomes and may be used to broaden gastroenteritis surveillance by capturing off-target enteric viruses. IMPORTANCE Viral gastroenteritis results in significant morbidity and mortality in vulnerable individuals and is primarily caused by norovirus. To investigate norovirus epidemiology, there is a need for whole-genome sequencing and reference sets consisting of full genomes. Using surveillance samples sent to the Dutch National Institute for Public Health and the Environment (RIVM), we compared metagenomics against conventional techniques, such as RT-qPCR and Sanger-sequencing, with norovirus as the target pathogen. We determined that metagenomics is a robust method to generate complete norovirus genomes, in parallel to many off-target pathogenic enteric virus genomes, thereby broadening our surveillance efforts. Moreover, we detected a sapovirus that was not detected by our validated gastroenteritis RT-qPCR panel, which exemplifies the strength of metagenomics. Our study shows that metagenomics can be used for public health gastroenteritis surveillance, the generation of reference-sets for molecular epidemiology, and how it compares to current surveillance strategies.

Culture-Independent Workflow for Nanopore MinION-Based Sequencing of Influenza A Virus

Whole-genome sequencing (WGS) of influenza A virus (IAV) is crucial for identifying diverse subtypes and newly evolved variants and for selecting vaccine strains. In developing countries, where facilities are often inadequate, WGS is challenging to perform using conventional next-generation sequencers. In this study, we established a culture-independent, high-throughput native barcode amplicon sequencing workflow that can sequence all influenza subtypes directly from a clinical specimen. All segments of IAV in 19 clinical specimens, irrespective of their subtypes, were amplified simultaneously using a two-step reverse transcriptase PCR (RT-PCR) system. First, the library was prepared using the ligation sequencing kit, barcoded individually using the native barcodes, and sequenced on the MinION MK 1C platform with real-time base-calling. Then, subsequent data analyses were performed with the appropriate tools. WGS of 19 IAV-positive clinical samples was carried out successfully with 100% coverage and 3,975-fold mean coverage for all segments. This easy-to-install and low-cost capacity-building protocol took only 24 h complete from extracting RNA to obtaining finished sequences. Overall, we developed a high-throughput portable sequencing workflow ideal for resource-limited clinical settings, aiding in real-time surveillance, outbreak investigation, and the detection of emerging viruses and genetic reassortment events. However, further evaluation is required to compare its accuracy with other high-throughput sequencing technologies to validate the widespread application of these findings, including WGS from environmental samples. IMPORTANCE The Nanopore MinION-based influenza sequencing approach we are proposing makes it possible to sequence the influenza A virus, irrespective of its diverse serotypes, directly from clinical and environmental swab samples, so that we are not limited to virus culture. This third-generation, portable, multiplexing, and real-time sequencing strategy is highly convenient for local sequencing, particularly in low- and middle-income countries like Bangladesh. Furthermore, the cost-efficient sequencing method could provide new opportunities to respond to the early phase of an influenza pandemic and enable the timely detection of the emerging subtypes in clinical samples. Here, we meticulously described the entire process that might help the researcher who will follow this methodology in the future. Our findings suggest that this proposed method is ideal for clinical and academic settings and will aid in real-time surveillance and in the detection of potential outbreak agents and newly evolved viruses.

Viral genome wide association study identifies novel hepatitis C virus polymorphisms associated with sofosbuvir treatment failure

Persistent hepatitis C virus (HCV) infection is a major cause of chronic liver disease, worldwide. With the development of direct-acting antivirals, treatment of chronically infected patients has become highly effective, although a subset of patients responds less well to therapy. Sofosbuvir is a common component of current de novo or salvage combination therapies, that targets the HCV NS5B polymerase. We use pre-treatment whole-genome sequences of HCV from 507 patients infected with HCV subtype 3a and treated with sofosbuvir containing regimens to detect viral polymorphisms associated with response to treatment. We find three common polymorphisms in non-targeted HCV NS2 and NS3 proteins are associated with reduced treatment response. These polymorphisms are enriched in post-treatment HCV sequences of patients unresponsive to treatment. They are also associated with lower reductions in viral load in the first week of therapy. Using in vitro short-term dose-response assays, these polymorphisms do not cause any reduction in sofosbuvir potency, suggesting an indirect mechanism of action in decreasing sofosbuvir efficacy. The identification of polymorphisms in NS2 and NS3 proteins associated with poor treatment outcomes emphasises the value of systematic genome-wide analyses of viruses in uncovering clinically relevant polymorphisms that impact treatment.

Microbiome analysis, the immune response and transplantation in the era of next generation sequencing

The human gastrointestinal tract, skin and mucosal surfaces are inhabited by a complex system of bacteria, viruses, fungi, archaea, protists, and eukaryotic parasites with predominance of bacteria and bacterial viruses (bacteriophages). Collectively these microbes form the microbiota of the microecosystem of humans. Recent advancement in technologies for nucleic acid isolation from various environmental samples, feces and body secretions and advancements in shotgun throughput massive parallel DNA and RNA sequencing along with 16S ribosomal gene sequencing have unraveled the identity of otherwise unknown microbial entities constituting the human microecosystem. The improved transcriptome analysis, technological developments in biochemical analytical methods and availability of complex bioinformatics tools have allowed us to begin to understand the metabolome of the microbiome and the biochemical pathways and potential signal transduction pathways in human cells in response to microbial infections and their products. Also, developments in human whole genome sequencing, targeted gene sequencing of histocompatibility genes and other immune response associated genes by Next Generation Sequencing (NGS) have allowed us to have a better conceptualization of immune responses, and alloimmune responses. These modern technologies have enabled us to dive into the intricate relationship between commensal symbiotic and pathogenic microbiome and immune system. For the most part, the commensal symbiotic microbiota helps to maintain normal immune homeostasis besides providing healthy nutrients, facilitating digestion, and protecting the skin, mucosal and intestinal barriers. However, changes in diets, administration of therapeutic agents like antibiotics, chemotherapeutic agents, immunosuppressants etc. along with certain host factors including human histocompatibility antigens may alter the microbial ecosystem balance by causing changes in microbial constituents, hierarchy of microbial species and even dysbiosis. Such alterations may cause immune dysregulation, breach of barrier protection and lead to immunopathogenesis rather than immune homeostasis. The effects of human microbiome on immunity, health and disease are currently under intense research with cutting edge technologies in molecular biology, biochemistry, and bioinformatics along with tremendous ability to characterize immune response at single cell level. This review will discuss the contemporary status on human microbiome immune system interactions and their potential effects on health, immune homeostasis and allograft transplantation.

Using Whole Genome Sequences to Investigate Adenovirus Outbreaks in a Hematopoietic Stem Cell Transplant Unit

A recent surge in human mastadenovirus (HAdV) cases, including five deaths, amongst a haematopoietic stem cell transplant population led us to use whole genome sequencing (WGS) to investigate. We compared sequences from 37 patients collected over a 20-month period with sequences from GenBank and our own database of HAdVs. Maximum likelihood trees and pairwise differences were used to evaluate genotypic relationships, paired with the epidemiological data from routine infection prevention and control (IPC) records and hospital activity data. During this time period, two formal outbreaks had been declared by IPC, while WGS detected nine monophyletic clusters, seven were corroborated by epidemiological evidence and by comparison of single-nucleotide polymorphisms. One of the formal outbreaks was confirmed, and the other was not. Of the five HAdV-associated deaths, three were unlinked and the remaining two considered the source of transmission. Mixed infection was frequent (10%), providing a sentinel source of recombination and superinfection. Immunosuppressed patients harboring a high rate of HAdV positivity require comprehensive surveillance. As a consequence of these findings, HAdV WGS is being incorporated routinely into clinical practice to influence IPC policy contemporaneously.

Conservation genomics reveals pesticide and pathogen exposure in the declining bumble bee Bombus terricola

In recent years, many pollinators have experienced large population declines, which threaten food security and the stability of natural ecosystems. Bumble bees are particularly important because their ability to “buzz” pollinate and tolerate cooler conditions make them critical pollinators for certain plants and regions. Here, we apply a conservation genomics approach to study the vulnerable Bombus terricola. We sequenced RNA from 30 worker abdomens, 18 of which were collected from agricultural sites and 12 of which were collected from nonagricultural sites. We found transcriptional signatures associated with exposure to insecticides, with gene expression patterns suggesting that bumble bees were exposed to neonicotinoids and/or fipronil—two compounds known to negatively impact bees. We also found transcriptional signatures associated with pathogen infections. In addition to the transcriptomic analysis, we carried out a metatranscriptomic analysis and detected five pathogens in the abdomens of workers, three of which are common in managed honey bee and bumble bee colonies. Our conservation genomics study provides functional support for the role of pesticides and pathogen spillover in the decline of B. terricola. We demonstrate that conservation genomics is an invaluable tool which allows researchers to quantify the effects of multiple stressors that impact pollinator populations in the wild.

Complete genome characterization of human noroviruses allows comparison of minor alleles during acute and chronic infections

Human noroviruses (HuNoVs) circulate globally, affect all age groups and place a substantial burden upon health services. High genetic diversity leading to antigenic variation plays a significant role in HuNoV epidemiology, driving periodic global emergence of epidemic variants. Studies have suggested that immunocompromised individuals may be a reservoir for such epidemic variants, but studies investigating the diversity and emergence of HuNoV variants in immunocompetent individuals are underrepresented. To address this, we sequenced the genomes of HuNoVs present in samples collected longitudinally from one immunocompetent (acute infection) and one immunocompromised (chronic infection) patient. A broadly reactive HuNoV capture-based method was used to concentrate the virus present in these specimens prior to massively parallel sequencing to recover near complete viral genomes. Using a novel bioinformatics pipeline, we demonstrated that persistent minor alleles were present in both acute and chronic infections, and that minor allele frequencies represented a larger proportion of the population during chronic infection. In acute infection, minor alleles were more evenly spread across the genome, although present at much lower frequencies, and therefore difficult to discern from error. By contrast, in the chronic infection, more minor alleles were present in the minor structural protein. No non-synonymous minor alleles were detected in the major structural protein over the short sampling period of the HuNoV chronic infection, suggesting where immune pressure is variable or non-existent, epidemic variants could emerge over longer periods of infection by random chance.

Predicting whole genome sequencing success for archived avian influenza virus (Orthomyxoviridae) samples using real-time and droplet PCRs

Long-term viral archives are valuable sources of research data. Each archive can store hundreds of thousands of diverse sample types. In the current era of whole genome sequencing, archived samples become a rich source of evolutionary and epidemiological data that can span years, and even decades. However, the ability to obtain high quality viral whole genome sequences from samples of various types, age, and quality is inconsistent. A minimum quality threshold that helps predict the best success of obtaining high quality genomic sequences for both recent and archived samples is highly valuable. Real-time reverse transcription PCR (rrt-PCR) and droplet digital PCR (ddPCR) are useful tools to evaluate nucleic acid integrity. We hypothesized that diagnostic rrt-PCR and ddPCR data for avian influenza virus (AIV) can predict viral whole genome sequencing success. To test this hypothesis we used RNA extracted from cloacal and oropharyngeal swabs stored in the USDA-APHIS National Wildlife Disease Program Wildlife Tissue Archive. We determined that a specific rrt-PCR C_q value or ddPCR copies/μL resulted in recovery of complete sequences of all eight AIV gene segments. We used logistic regression to estimate probabilities of whole genome recovery at 0.95 (C_q = 15, copies/μL = 49,350), 0.75 (C_q = 24, copies/μL = 16,800), 0.50 (C_q = 29, copies/μL = <1), and 0.25 (C_q = 235, copies/μL = <1). We also identified values at which we predictably recovered HA and NA segments for diagnosing subtypes (C_q = 27.29; copies/μL = 757.50). This approach will allow researchers to assess the potential success of AIV whole genome recovery from diagnostic samples collected in routine AIV surveillance.

The AC4 Protein of a Cassava Geminivirus Is Required for Virus Infection

Geminiviruses (family Geminiviridae) are among the most devastating plant viruses worldwide, causing severe damage in crops of economic and subsistence importance. These viruses have very compact genomes and many of the encoded proteins are multifunctional. Here, we investigated the role of the East African cassava mosaic Cameroon virus (EACMCV) AC4 on virus infectivity in Nicotiana benthamiana. Results showed that plants inoculated with EACMCV containing a knockout mutation in an AC4 open reading frame displayed symptoms 2 to 3 days later than plants inoculated with wild-type virus, and these plants recovered from infection, whereas plants inoculated with the wild-type virus did not. Curiously, when an additional mutation was made in the knockout mutant, the resulting double mutant virus completely failed to cause any apparent symptoms. Interestingly, the role of AC4 on virus infectivity appeared to be dependent on an encoded N-myristoylation motif that mediates cell membrane binding. We previously showed that EACMCV containing the AC4T38I mutant produced virus progeny characterized by second-site mutations and reversion to wild-type virus. These results were confirmed in this study using additional mutations. Together, these results show involvement of EACMCV AC4 in virus infectivity; they also suggest a role for the combined action of mutation and selection, under prevailing environmental conditions, on begomovirus genetic variation and diversity.

Taking gastro-surveillance into the 21st century

Enteric viruses, particularly rotaviruses and noroviruses, are leading causes of gastroenteritis worldwide. Human rotaviruses are ubiquitous and globally almost every child has been infected by 3-5 years of age. Noroviruses affect people of all ages and is the leading cause of foodborne outbreaks. Rota- and noroviruses account for ˜40% and ˜17% of diarrhea-associated hospitalizations, and ˜200,000 deaths annually respectively, with most deaths occurring in developing countries. Two rotavirus vaccines have currently been implemented in ˜95 countries and several norovirus vaccine candidates are currently in development and/or clinical testing. Surveillance of enteric viruses is an important part of outbreak investigations as well as pre- and post-vaccine impact studies but is even in developed countries often limited to investigation of sporadic cases or comprehensive outbreaks. Conventional methods for enteric virus detection and subtyping relies on standard RT-PCR methods, supplemented with Sanger-sequencing. However, for viruses with even moderate mutationrates, PCR-based-typing of only limited parts of the virus genome is challenging and requires regular update of primers. Full-genomecharacterization technologies based on sequence independent methods based on next generation sequencing (NGS), have demonstrated great potential for enteric virus detection and/or typing in both clinical and environmental samples. However, cost-benefits must balance for such methods to be widely accepted for public health purposes. In Europe as also globally, routine use of NGS-methods for surveillance of enteric viruses is currently limited to few national public health laboratories. What important lessons can be learned from these and what is the future of NGS-based surveillance?

Resistance analysis of genotype 3 hepatitis C virus indicates subtypes inherently resistant to nonstructural protein 5A inhibitors

Hepatitis C virus (HCV) genotype (gt) 3 is highly prevalent globally, with non-gt3a subtypes common in Southeast Asia. Resistance-associated substitutions (RASs) have been shown to play a role in treatment failure. However, the role of RASs in gt3 is not well understood. We report the prevalence of RASs in a cohort of direct-acting antiviral treatment-naive, gt3-infected patients, including those with rarer subtypes, and evaluate the effect of these RASs on direct-acting antivirals in vitro. Baseline samples from 496 gt3 patients enrolled in the BOSON clinical trial were analyzed by next-generation sequencing after probe-based enrichment for HCV. Whole viral genomes were analyzed for the presence of RASs to approved direct-acting antivirals. The resistance phenotype of RASs in combination with daclatasvir, velpatasvir, pibrentasvir, elbasvir, and sofosbuvir was measured using the S52 ΔN gt3a replicon model. The nonstructural protein 5A A30K and Y93H substitutions were the most common at 8.9% (n = 44) and 12.3% (n = 61), respectively, and showed a 10-fold and 11-fold increase in 50% effect concentration for daclatasvir compared to the unmodified replicon. Paired RASs (A30K + L31M and A30K + Y93H) were identified in 18 patients (9 of each pair); these combinations were shown to be highly resistant to daclatasvir, velpatasvir, elbasvir, and pibrentasvir. The A30K + L31M combination was found in all gt3b and gt3g samples. Conclusion: Our study reveals high frequencies of RASs to nonstructural protein 5A inhibitors in gt3 HCV; the paired A30K + L31M substitutions occur in all patients with gt3b and gt3g virus, and in vitro analysis suggests that these subtypes may be inherently resistant to all approved nonstructural protein 5A inhibitors for gt3 HCV. (Hepatology 2018).

Interpreting Viral Deep Sequencing Data with GLUE

Using deep sequencing technologies such as Illumina's platform, it is possible to obtain reads from the viral RNA population revealing the viral genome diversity within a single host. A range of software tools and pipelines can transform raw deep sequencing reads into Sequence Alignment Mapping (SAM) files. We propose that interpretation tools should process these SAM files, directly translating individual reads to amino acids in order to extract statistics of interest such as the proportion of different amino acid residues at specific sites. This preserves per-read linkage between nucleotide variants at different positions within a codon location. The samReporter is a subsystem of the GLUE software toolkit which follows this direct read translation approach in its processing of SAM files. We test samReporter on a deep sequencing dataset obtained from a cohort of 241 UK HCV patients for whom prior treatment with direct-acting antivirals has failed; deep sequencing and resistance testing have been suggested to be of clinical use in this context. We compared the polymorphism interpretation results of the samReporter against an approach that does not preserve per-read linkage. We found that the samReporter was able to properly interpret the sequence data at resistance-associated locations in nine patients where the alternative approach was equivocal. In three cases, the samReporter confirmed that resistance or an atypical substitution was present at NS5A position 30. In three further cases, it confirmed that the sofosbuvir-resistant NS5B substitution S282T was absent. This suggests the direct read translation approach implemented is of value for interpreting viral deep sequencing data.

Understanding and overcoming the pitfalls and biases of next-generation sequencing (NGS) methods for use in the routine clinical microbiological diagnostic laboratory

Recent advancements in next-generation sequencing (NGS) have provided the foundation for modern studies into the composition of microbial communities. The use of these NGS methods allows for the detection and identification of (‘difficult-to-culture’) microorganisms using a culture-independent strategy. In the field of routine clinical diagnostics however, the application of NGS is currently limited to microbial strain typing for epidemiological purposes only, even though the implementation of NGS for microbial community analysis may yield clinically important information. This lack of NGS implementation is due to many different factors, including issues relating to NGS method standardization and result reproducibility. In this review article, the authors provide a general introduction to the most widely used NGS methods currently available (i.e., targeted amplicon sequencing and shotgun metagenomics) and the strengths and weaknesses of each method is discussed. The focus of the publication then shifts toward 16S rRNA gene NGS methods, which are currently the most cost-effective and widely used NGS methods for research purposes, and are therefore more likely to be successfully implemented into routine clinical diagnostics in the short term. In this respect, the experimental pitfalls and biases created at each step of the 16S rRNA gene NGS workflow are explained, as well as their potential solutions. Finally, a novel diagnostic microbiota profiling platform (‘MYcrobiota’) is introduced, which was developed by the authors by taking into consideration the pitfalls, biases, and solutions explained in this article. The development of the MYcrobiota, and future NGS methodologies, will help pave the way toward the successful implementation of NGS methodologies into routine clinical diagnostics.

Low-Bias RNA Sequencing of the HIV-2 Genome from Blood Plasma

Accurate determination of the genetic diversity present in the HIV quasispecies is critical for the development of a preventative vaccine: in particular, little is known about viral genetic diversity for the second type of HIV, HIV-2. A better understanding of HIV-2 biology is relevant to the HIV vaccine field because a substantial proportion of infected people experience long-term viral control, and prior HIV-2 infection has been associated with slower HIV-1 disease progression in coinfected subjects. The majority of traditional and next-generation sequencing methods have relied on target amplification prior to sequencing, introducing biases that may obscure the true signals of diversity in the viral population. Additionally, target enrichment through PCR requires a priori sequence knowledge, which is lacking for HIV-2. Therefore, a target enrichment free method of library preparation would be valuable for the field. We applied an RNA shotgun sequencing (RNA-Seq) method without PCR amplification to cultured viral stocks and patient plasma samples from HIV-2-infected individuals. Libraries generated from total plasma RNA were analyzed with a two-step pipeline: (i) de novo genome assembly, followed by (ii) read remapping. By this approach, whole-genome sequences were generated with a 28× to 67× mean depth of coverage. Assembled reads showed a low level of GC bias, and comparison of the genome diversities at the intrahost level showed low diversity in the accessory gene vpx in all patients. Our study demonstrates that RNA-Seq is a feasible full-genome de novo sequencing method for blood plasma samples collected from HIV-2-infected individuals.IMPORTANCE An accurate picture of viral genetic diversity is critical for the development of a globally effective HIV vaccine. However, sequencing strategies are often complicated by target enrichment prior to sequencing, introducing biases that can distort variant frequencies, which are not easily corrected for in downstream analyses. Additionally, detailed a priori sequence knowledge is needed to inform robust primer design when employing PCR amplification, a factor that is often lacking when working with tropical diseases localized in developing countries. Previous work has demonstrated that direct RNA shotgun sequencing (RNA-Seq) can be used to circumvent these issues for hepatitis C virus (HCV) and norovirus. We applied RNA-Seq to total RNA extracted from HIV-2 blood plasma samples, demonstrating the applicability of this technique to HIV-2 and allowing us to generate a dynamic picture of genetic diversity over the whole genome of HIV-2 in the context of low-bias sequencing.

Genetic characterization of norovirus GII.4 variants circulating in Canada using a metagenomic technique

BACKGROUND

Human norovirus is the leading cause of viral gastroenteritis globally, and the GII.4 has been the most predominant genotype for decades. This genotype has numerous variants that have caused repeated epidemics worldwide. However, the molecular evolutionary signatures among the GII.4 variants have not been elucidated throughout the viral genome.

METHOD

A metagenomic, next-generation sequencing method, based on Illumina RNA-Seq, was applied to determine norovirus sequences from clinical samples.

RESULTS

Herein, the obtained deep-sequencing data was employed to analyze full-genomic sequences from GII.4 variants prevailing in Canada from 2012 to 2016. Phylogenetic analysis demonstrated that the majority of these sequences belong to New Orleans 2009 and Sydney 2012 strains, and a recombinant sequence was also identified. Genome-wide similarity analyses implied that while the capsid gene is highly diverse among the isolates, the viral protease and polymerase genes remain relatively conserved. Numerous amino acid substitutions were observed at each putative antigenic epitope of the VP1 protein, whereas few amino acid changes were identified in the polymerase protein. Co-infection with other enteric RNA viruses was investigated and the astrovirus genome was identified in one of the samples.

CONCLUSIONS

Overall this study demonstrated the application of whole genome sequencing as an important tool in molecular characterization of noroviruses.

Characterization of hepatitis C virus resistance to grazoprevir reveals complex patterns of mutations following on-treatment breakthrough that are not observed at relapse

PURPOSE

A detailed analysis of hepatitis C virus (HCV) resistance-associated substitutions (RASs) is required to understand why people fail direct-acting antiviral therapies. This study was conducted to assess RASs in an analysis of 2 trials evaluating the second-generation NS3/4A protease inhibitor grazoprevir (GZR) in combination with peginterferon/ribavirin.

PATIENTS AND METHODS

From a total of 113 participants with HCV genotype 1 infection, RASs were evaluated in 25 patients who relapsed and 6 patients with on-treatment virologic breakthrough using consensus Sanger and clonal sequence analysis of NS3/NS4a genes, with in vitro testing of replicon mutants. Next-generation sequencing (NGS) was used in a subset of participants to assess minority variants and the evolution of the whole viral genome.

RESULTS

Baseline RASs did not predict treatment failure. Relapse was most commonly associated with RASs at D168, although additional RASs (Y56, R155 and A156) were also detected, particularly in participants with on-treatment breakthrough. Treatment-emergent RASs usually reverted to wild-type (WT), suggesting these mutations were associated with a negative fitness cost (confirmed using in vitro assays). NGS was the most sensitive assay for the detection of minor variants. Significant viral sequence divergence (up to 5.9% codons) was observed across whole genomes in association with the acquisition and reversion of RASs.

CONCLUSION

Relapse with GZR and peginterferon/ribavirin is commonly associated with single RASs in NS3 that generally revert to WT, while breakthrough follows more complex patterns of viral resistance. NGS suggests that large diverse pools of viral quasispecies that emerge with RASs facilitate rapid viral evolution.

A near full-length open reading frame next generation sequencing assay for genotyping and identification of resistance-associated variants in hepatitis C virus

BACKGROUND

The current treatment options for hepatitis C virus (HCV), based on direct acting antivirals (DAA), are dependent on virus genotype and previous treatment experience. Treatment failures have been associated with detection of resistance-associated substitutions (RASs) in the DAA targets of HCV, the NS3, NS5A and NS5 B proteins.

OBJECTIVE

To develop a next generation sequencing based method that provides genotype and detection of HCV NS3, NS5A, and NS5 B RASs without prior knowledge of sample genotype.

STUDY DESIGN

In total, 101 residual plasma samples from patients with HCV covering 10 different viral subtypes across 4 genotypes with viral loads of 3.84-7.61 Log IU/mL were included. All samples were de-identified and consequently prior treatment status for patients was unknown. Almost full open reading frame amplicons (∼ 9 kb) were generated using RT-PCR with a single primer set. The resulting amplicons were sequenced with high throughput sequencing and analysed using an in-house developed script for detecting RASs.

RESULTS

The method successfully amplified and sequenced 94% (95/101) of samples with an average coverage of 14,035; four of six failed samples were genotype 4a. Samples analysed twice yielded reproducible nucleotide frequencies across all sites. RASs were detected in 21/95 (22%) samples at a 15% threshold. The method identified one patient infected with two genotype 2b variants, and the presence of subgenomic deletion variants in 8 (8.4%) of 95 successfully sequenced samples.

CONCLUSIONS

The presented method may provide identification of HCV genotype, RASs detection, and detect multiple HCV infection without prior knowledge of sample genotype.

Microfluidic PCR Amplification and MiSeq Amplicon Sequencing Techniques for High-Throughput Detection and Genotyping of Human Pathogenic RNA Viruses in Human Feces, Sewage, and Oysters

Detection and genotyping of pathogenic RNA viruses in human and environmental samples are useful for monitoring the circulation and prevalence of these pathogens, whereas a conventional PCR assay followed by Sanger sequencing is time-consuming and laborious. The present study aimed to develop a high-throughput detection-and-genotyping tool for 11 human RNA viruses [Aichi virus; astrovirus; enterovirus; norovirus genogroup I (GI), GII, and GIV; hepatitis A virus; hepatitis E virus; rotavirus; sapovirus; and human parechovirus] using a microfluidic device and next-generation sequencer. Microfluidic nested PCR was carried out on a 48.48 Access Array chip, and the amplicons were recovered and used for MiSeq sequencing (Illumina, Tokyo, Japan); genotyping was conducted by homology searching and phylogenetic analysis of the obtained sequence reads. The detection limit of the 11 tested viruses ranged from 10⁰ to 10³ copies/μL in cDNA sample, corresponding to 10¹–10⁴ copies/mL-sewage, 10⁵–10⁸ copies/g-human feces, and 10²–10⁵ copies/g-digestive tissues of oyster. The developed assay was successfully applied for simultaneous detection and genotyping of RNA viruses to samples of human feces, sewage, and artificially contaminated oysters. Microfluidic nested PCR followed by MiSeq sequencing enables efficient tracking of the fate of multiple RNA viruses in various environments, which is essential for a better understanding of the circulation of human pathogenic RNA viruses in the human population.

HIV-1 full-genome phylogenetics of generalized epidemics in sub-Saharan Africa: impact of missing nucleotide characters in next-generation sequences

To characterize HIV-1 transmission dynamics in regions where the burden of HIV-1 is greatest, the “Phylogenetics and Networks for Generalised HIV Epidemics in Africa” consortium (PANGEA-HIV) is sequencing full-genome viral isolates from across sub-Saharan Africa. We report the first 3,985 PANGEA-HIV consensus sequences from four cohort sites (Rakai Community Cohort Study, n = 2,833; MRC/UVRI Uganda, n = 701; Mochudi Prevention Project, n = 359; Africa Health Research Institute Resistance Cohort, n = 92). Next-generation sequencing success rates varied: more than 80% of the viral genome from the gag to the nef genes could be determined for all sequences from South Africa, 75% of sequences from Mochudi, 60% of sequences from MRC/UVRI Uganda, and 22% of sequences from Rakai. Partial sequencing failure was primarily associated with low viral load, increased for amplicons closer to the 3′ end of the genome, was not associated with subtype diversity except HIV-1 subtype D, and remained significantly associated with sampling location after controlling for other factors. We assessed the impact of the missing data patterns in PANGEA-HIV sequences on phylogeny reconstruction in simulations. We found a threshold in terms of taxon sampling below which the patchy distribution of missing characters in next-generation sequences (NGS) has an excess negative impact on the accuracy of HIV-1 phylogeny reconstruction, which is attributable to tree reconstruction artifacts that accumulate when branches in viral trees are long. The large number of PANGEA-HIV sequences provides unprecedented opportunities for evaluating HIV-1 transmission dynamics across sub-Saharan Africa and identifying prevention opportunities. Molecular epidemiological analyses of these data must proceed cautiously because sequence sampling remains below the identified threshold and a considerable negative impact of missing characters on phylogeny reconstruction is expected.

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.

Efficient and unbiased metagenomic recovery of RNA virus genomes from human plasma samples

RNA viruses cause significant human pathology and are responsible for the majority of emerging zoonoses. Mainstream diagnostic assays are challenged by their intrinsic diversity, leading to false negatives and incomplete characterisation. New sequencing techniques are expanding our ability to agnostically interrogate nucleic acids within diverse sample types, but in the clinical setting are limited by overwhelming host material and ultra-low target frequency. Through selective host RNA depletion and compensatory protocol adjustments for ultra-low RNA inputs, we are able to detect three major blood-borne RNA viruses – HIV, HCV and HEV. We recovered complete genomes and up to 43% of the genome from samples with viral loads of 10⁴ and 10³ IU/ml respectively. Additionally, we demonstrated the utility of this method in detecting and characterising members of diverse RNA virus families within a human plasma background, some present at very low levels. By applying this method to a patient sample series, we have simultaneously determined the full genome of both a novel subtype of HCV genotype 6, and a co-infecting human pegivirus. This method builds upon earlier RNA metagenomic techniques and can play an important role in the surveillance and diagnostics of blood-borne viruses.

Evolution of selective-sequencing approaches for virus discovery and virome analysis

•

Description of virus enrichment techniques for metagenomics based virome analysis.

•

Usefulness of recently developed virome capture sequencing techniques.

•

Perspective on negative and positive selection approaches for virome analysis.

Recent advances in sequencing technologies have transformed the field of virus discovery and virome analysis. Once mostly confined to the traditional Sanger sequencing based individual virus discovery, is now entirely replaced by high throughput sequencing (HTS) based virus metagenomics that can be used to characterize the nature and composition of entire viromes. To better harness the potential of HTS for the study of viromes, sample preparation methodologies use different approaches to exclude amplification of non-viral components that can overshadow low-titer viruses. These virus-sequence enrichment approaches mostly focus on the sample preparation methods, like enzymatic digestion of non-viral nucleic acids and size exclusion of non-viral constituents by column filtration, ultrafiltration or density gradient centrifugation. However, recently a new approach of virus-sequence enrichment called virome-capture sequencing, focused on the amplification or HTS library preparation stage, was developed to increase the ability of virome characterization. This new approach has the potential to further transform the field of virus discovery and virome analysis, but its technical complexity and sequence-dependence warrants further improvements. In this review we discuss the different methods, their applications and evolution, for selective sequencing based virome analysis and also propose refinements needed to harness the full potential of HTS for virome analysis.

Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples

Genome sequencing has become a powerful tool for studying emerging infectious diseases; however, genome sequencing directly from clinical samples (i.e., without isolation and culture) remains challenging for viruses such as Zika, for which metagenomic sequencing methods may generate insufficient numbers of viral reads. Here we present a protocol for generating coding-sequence-complete genomes, comprising an online primer design tool, a novel multiplex PCR enrichment protocol, optimized library preparation methods for the portable MinION sequencer (Oxford Nanopore Technologies) and the Illumina range of instruments, and a bioinformatics pipeline for generating consensus sequences. The MinION protocol does not require an Internet connection for analysis, making it suitable for field applications with limited connectivity. Our method relies on multiplex PCR for targeted enrichment of viral genomes from samples containing as few as 50 genome copies per reaction. Viral consensus sequences can be achieved in 1-2 d by starting with clinical samples and following a simple laboratory workflow. This method has been successfully used by several groups studying Zika virus evolution and is facilitating an understanding of the spread of the virus in the Americas. The protocol can be used to sequence other viral genomes using the online Primal Scheme primer designer software. It is suitable for sequencing either RNA or DNA viruses in the field during outbreaks or as an inexpensive, convenient method for use in the lab.

A universal primer-independent next-generation sequencing approach for investigations of norovirus outbreaks and novel variants

Norovirus (NoV) is the most common cause of non-bacterial gastroenteritis and is a major agent associated with outbreaks of gastroenteritis. Conventional molecular genotyping analysis of NoV, used for the identification of transmission routes, relies on standard typing methods (STM) by Sanger-sequencing of only a limited part of the NoV genome, which could lead to wrong conclusions. Here, we combined a NoV capture method with next generation sequencing (NGS), which increased the proportion of norovirus reads by ~40 fold compared to NGS without prior capture. Of 15 NoV samples from 6 single-genotype outbreaks, near full-genome coverage (>90%) was obtained from 9 samples. Fourteen polymerase (RdRp) and 15 capsid (cap) genotypes were identified compared to 12 and 13 for the STM, respectively. Analysis of 9 samples from two mixed-genotype outbreaks identified 6 RdRp and 6 cap genotypes (two at >90% NoV genome coverage) compared to 4 and 2 for the STM, respectively. Furthermore, complete or partial sequences from the P2 hypervariable region were obtained from 7 of 8 outbreaks and a new NoV recombinant was identified. This approach could therefore strengthen outbreak investigations and could be applied to other important viruses in stool samples such as hepatitis A and enterovirus.

Characterization of the Genomic Diversity of Norovirus in Linked Patients Using a Metagenomic Deep Sequencing Approach

Norovirus (NoV) is the leading cause of gastroenteritis worldwide. A robust cell culture system does not exist for NoV and therefore detailed characterization of outbreak and sporadic strains relies on molecular techniques. In this study, we employed a metagenomic approach that uses non-specific amplification followed by next-generation sequencing to whole genome sequence NoV genomes directly from clinical samples obtained from 8 linked patients. Enough sequencing depth was obtained for each sample to use a de novo assembly of near-complete genome sequences. The resultant consensus sequences were then used to identify inter-host nucleotide variations that occur after direct transmission, analyze amino acid variations in the major capsid protein, and provide evidence of recombination events. The analysis of intra-host quasispecies diversity was possible due to high coverage-depth. We also observed a linear relationship between NoV viral load in the clinical sample and the number of sequence reads that could be attributed to NoV. The method demonstrated here has the potential for future use in whole genome sequence analyses of other RNA viruses isolated from clinical, environmental, and food specimens.

High-throughput whole genome sequencing of Porcine reproductive and respiratory syndrome virus from cell culture materials and clinical specimens using next-generation sequencing technology

Next-generation sequencing (NGS) technologies have increasingly played crucial roles in biological and medical research, but are not yet in routine use in veterinary diagnostic laboratories. We developed and applied a procedure for high-throughput RNA sequencing of Porcine reproductive and respiratory syndrome virus (PRRSV) from cell culture-derived isolates and clinical specimens. Ten PRRSV isolates with known sequence information, 2 mixtures each with 2 different PRRSV isolates, and 51 clinical specimens (19 sera, 16 lungs, and 16 oral fluids) with various PCR threshold cycle (Ct) values were subjected to nucleic acid extraction, cDNA library preparation (24-plexed), and sequencing. Whole genome sequences were obtained from 10 reference isolates with expected sequences and from sera with a PRRSV real-time reverse transcription PCR Ct ≤ 23.6, lung tissues with Ct ≤ 21, and oral fluids with Ct ≤ 20.6. For mixtures with PRRSV-1 and -2 isolates (57.8% nucleotide identity), NGS was able to distinguish them as well as obtain their respective genome sequences. For mixtures with 2 PRRSV-2 isolates (92.4% nucleotide identity), sequence reads with nucleotide ambiguity at numerous sites were observed, indicating mixed infection; however, individual virus sequences could only be separated when 1 isolate identity and sequence in the mixture is known. The NGS approach described herein offers the prospect of high-throughput sequencing and could be adapted to routine workflows in veterinary diagnostic laboratories, although further improvement of sequencing outcomes from clinical specimens with higher Ct values remains to be investigated.

Norovirus Whole-Genome Sequencing by SureSelect Target Enrichment: a Robust and Sensitive Method

Norovirus full-genome sequencing is challenging due to sequence heterogeneity among genomes. Previous methods have relied on PCR amplification, which is problematic due to primer design, and transcriptome sequencing (RNA-Seq), which nonspecifically sequences all RNA, including host and bacterial RNA, in stool specimens. Target enrichment uses a panel of custom-designed 120-mer RNA baits that are complementary to all publicly available norovirus sequences, with multiple baits targeting each position of the genome, which overcomes the challenge of primer design. Norovirus genomes are enriched from stool RNA extracts to minimize the sequencing of nontarget RNA. SureSelect target enrichment and Illumina sequencing were used to sequence full genomes from 507 norovirus-positive stool samples with reverse transcription-real-time PCR cycle threshold (CT) values of 10 to 43. Sequencing on an Illumina MiSeq system in batches of 48 generated, on average, 81% on-target reads per sample and 100% genome coverage with >12,000-fold read depth. Samples included genotypes GI.1, GI.2, GI.3, GI.6, GI.7, GII.1, GII.2, GII.3, GII.4, GII.5, GII.6, GII.7, GII.13, GII.14, and GII.17. When outliers were accounted for, we generated >80% genome coverage for all positive samples, regardless of CT values. A total of 164 samples were tested in parallel with conventional PCR genotyping of the capsid shell domain; 164/164 samples were successfully sequenced, compared to 158/164 samples that were amplified by PCR. Four of the samples that failed capsid PCR analysis had low titers, which suggests that target enrichment is more sensitive than gel-based PCR. Two samples failed PCR due to primer mismatches; target enrichment uses multiple baits targeting each position, thus accommodating sequence heterogeneity among norovirus genomes.

Highly-Immunogenic Virally-Vectored T-cell Vaccines Cannot Overcome Subversion of the T-cell Response by HCV during Chronic Infection

An effective therapeutic vaccine for the treatment of chronic hepatitis C virus (HCV) infection, as an adjunct to newly developed directly-acting antivirals (DAA), or for the prevention of reinfection, would significantly reduce the global burden of disease associated with chronic HCV infection. A recombinant chimpanzee adenoviral (ChAd3) vector and a modified vaccinia Ankara (MVA), encoding the non-structural proteins of HCV (NSmut), used in a heterologous prime/boost regimen induced multi-specific, high-magnitude, durable HCV-specific CD4+ and CD8+ T-cell responses in healthy volunteers, and was more immunogenic than a heterologous Ad regimen. We now assess the immunogenicity of this vaccine regimen in HCV infected patients (including patients with a low viral load suppressed with interferon/ribavirin therapy), determine T-cell cross-reactivity to endogenous virus, and compare immunogenicity with that observed previously in both healthy volunteers and in HCV infected patients vaccinated with the heterologous Ad regimen. Vaccination of HCV infected patients with ChAd3-NSmut/MVA-NSmut was well tolerated. Vaccine-induced HCV-specific T-cell responses were detected in 8/12 patients; however, CD4+ T-cell responses were rarely detected, and the overall magnitude of HCV-specific T-cell responses was markedly reduced when compared to vaccinated healthy volunteers. Furthermore, HCV-specific cells had a distinct partially-functional phenotype (lower expression of activation markers, granzyme B, and TNFα production, weaker in vitro proliferation, and higher Tim3 expression, with comparable Tbet and Eomes expression) compared to healthy volunteers. Robust anti-vector T-cells and antibodies were induced, showing that there is no global defect in immunity. The level of viremia at the time of vaccination did not correlate with the magnitude of the vaccine-induced T-cell response. Full-length, next-generation sequencing of the circulating virus demonstrated that T-cells were only induced by vaccination when there was a sequence mismatch between the autologous virus and the vaccine immunogen. However, these T-cells were not cross-reactive with the endogenous viral variant epitopes. Conversely, when there was complete homology between the immunogen and circulating virus at a given epitope T-cells were not induced. T-cell induction following vaccination had no significant impact on HCV viral load. In vitro T-cell culture experiments identified the presence of T-cells at baseline that could be expanded by vaccination; thus, HCV-specific T-cells may have been expanded from pre-existing low-level memory T-cell populations that had been exposed to HCV antigens during natural infection, explaining the partial T-cell dysfunction. In conclusion, vaccination with ChAd3-NSmut and MVA-NSmut prime/boost, a potent vaccine regimen previously optimized in healthy volunteers was unable to reconstitute HCV-specific T-cell immunity in HCV infected patients. This highlights the major challenge of overcoming T-cell exhaustion in the context of persistent antigen exposure.

A next generation sequencing-based method to study the intra-host genetic diversity of norovirus in patients with acute and chronic infection

Background

Immunocompromised individuals with chronic norovirus (NoV) infection and elderly patients are hypothesized to be reservoirs where NoV might accumulate mutations and evolve into pandemic strains. Next generation sequencing (NGS) methods can monitor the intra-host diversity of NoV and its evolution but low abundance of viral RNA results in sub-optimal efficiency. In this study, we: 1) established a next generation sequencing-based method for NoV using bacterial rRNA depletion as a viral RNA enrichment strategy, and 2) measured the intra-host genetic diversity of NoV in specimens of patients with acute NoV infection (n = 4) and in longitudinal specimens of an immunocompromised patient with chronic NoV infection (n = 2).

Results

A single Illumina MiSeq dataset resulted in near full-length genome sequences for 5 out of 6 multiplexed samples. Experimental depletion of bacterial rRNA in stool RNA provided up to 1.9 % of NoV reads. The intra-host viral population in patients with acute NoV infection was homogenous and no single nucleotide variants (SNVs) were detected. In contrast, the NoV population from the immunocompromised patient was highly diverse and accumulated SNVs over time (51 SNVs in the first sample and 122 SNVs in the second sample collected 4 months later). The percentages of SNVs causing non-synonymous mutations were 27.5 % and 20.5 % for the first and second samples, respectively. The majority of non-synonymous mutations occurred, in increasing order of frequency, in p22, the major capsid (VP1) and minor capsid (VP2) genes.

Conclusions

The results provide data useful for the selection and improvement of NoV RNA enrichment strategies for NGS. Whole genome analysis using next generation sequencing confirmed that the within-host population of NoV in an immunocompromised individual with chronic NoV infection was more diverse compared to that in individuals with acute infection. We also observed an accumulation of non-synonymous mutations at the minor capsid gene that has not been reported in previous studies and might have a role in NoV adaptation.

Electronic supplementary material

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

Molecular Diagnostic Methods for Detection and Characterization of Human Noroviruses

Human noroviruses are a group of viral agents that afflict people of all age groups. The viruses are now recognized as the most common causative agent of nonbacterial acute gastroenteritis and foodborne viral illness worldwide. However, they have been considered to play insignificant roles in the disease burden of acute gastroenteritis for the past decades until the recent advent of new and more sensitive molecular diagnostic methods. The availability and application of the molecular diagnostic methods have led to enhanced detection of noroviruses in clinical, food and environmental samples, significantly increasing the recognition of noroviruses as an etiologic agent of epidemic and sporadic acute gastroenteritis. This article aims to summarize recent efforts made for the development of molecular methods for the detection and characterization of human noroviruses.

Direct next-generation sequencing of virus-human mixed samples without pretreatment is favorable to recover virus genome

Abstract

Next-generation sequencing (NGS) enables the recovery of pathogen genomes from clinical samples without the need for culturing. Depletion of host/microbiota components (e.g., ribosomal RNA and poly-A RNA) and whole DNA/cDNA amplification are routine methods to improve recovery results. Using mixtures of human and influenza A virus (H1N1) RNA as a model, we found that background depletion and whole transcriptome amplification introduced biased distributions of read coverage over the H1N1 genome, thereby hampering genome assembly. Influenza serotyping was also affected by pretreatments. We propose that direct sequencing of noncultured samples without pretreatment is a favorable option for pathogen genome recovery applications.

Reviewer

This article was reviewed by Sebastian Maurer-Stroh.

Electronic supplementary material

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

Quantifying Next Generation Sequencing Sample Pre-Processing Bias in HIV-1 Complete Genome Sequencing

Genetic analyses play a central role in infectious disease research. Massively parallelized “mechanical cloning” and sequencing technologies were quickly adopted by HIV researchers in order to broaden the understanding of the clinical importance of minor drug-resistant variants. These efforts have, however, remained largely limited to small genomic regions. The growing need to monitor multiple genome regions for drug resistance testing, as well as the obvious benefit for studying evolutionary and epidemic processes makes complete genome sequencing an important goal in viral research. In addition, a major drawback for NGS applications to RNA viruses is the need for large quantities of input DNA. Here, we use a generic overlapping amplicon-based near full-genome amplification protocol to compare low-input enzymatic fragmentation (Nextera™) with conventional mechanical shearing for Roche 454 sequencing. We find that the fragmentation method has only a modest impact on the characterization of the population composition and that for reliable results, the variation introduced at all steps of the procedure—from nucleic acid extraction to sequencing—should be taken into account, a finding that is also relevant for NGS technologies that are now more commonly used. Furthermore, by applying our protocol to deep sequence a number of pre-therapy plasma and PBMC samples, we illustrate the potential benefits of a near complete genome sequencing approach in routine genotyping.

ve-SEQ: Robust, unbiased enrichment for streamlined detection and whole-genome sequencing of HCV and other highly diverse pathogens

The routine availability of high-depth virus sequence data would allow the sensitive detection of resistance-associated variants that can jeopardize HIV or hepatitis C virus (HCV) treatment. We introduce ve-SEQ, a high-throughput method for sequence-specific enrichment and characterization of whole-virus genomes at up to 20% divergence from a reference sequence and 1,000-fold greater sensitivity than direct sequencing. The extreme genetic diversity of HCV led us to implement an algorithm for the efficient design of panels of oligonucleotide probes to capture any sequence among a defined set of targets without detectable bias. ve-SEQ enables efficient detection and sequencing of any HCV genome, including mixtures and intra-host variants, in a single experiment, with greater tolerance of sequence diversity than standard amplification methods and greater sensitivity than metagenomic sequencing, features that are directly applicable to other pathogens or arbitrary groups of target organisms, allowing the combination of sensitive detection with sequencing in many settings.

Full genome analysis of bovine astrovirus from fecal samples of cattle in Japan: identification of possible interspecies transmission of bovine astrovirus

A viral metagenomics approach was used to investigate fecal samples of Japanese calves with and without diarrhea. Of the different viral pathogens detected, read counts gave nearly complete astrovirus-related RNA sequences in 15 of the 146 fecal samples collected in three distinct areas (Hokkaido, Ishikawa, and Kagoshima Prefectures) between 2009 and 2015. Due to the lack of genetic information about bovine astroviruses (BoAstVs) in Japan, these sequences were analyzed in this study. Nine of the 15 Japanese BoAstVs were closely related to Chinese BoAstVs and clustered into a lineage (tentatively named lineage 1) in all phylogenetic trees. Three of 15 strains were phylogenetically separate from lineage 1, showing low sequence identities, and clustered instead with an American strain isolated from cattle with respiratory disease (tentatively named lineage 2). Interestingly, two of 15 strains clustered with lineage 1 in the open reading frame (ORF)1a and ORF1b regions, while they clustered with lineage 2 in the ORF2 region. Remarkably, one of 15 strains exhibited low amino acid sequence similarity to other BoAstVs and was clustered separately with porcine astrovirus type 5 in all trees, and ovine astrovirus in the ORF2 region, suggesting past interspecies transmission.

Characterization of Hepatitis C Virus Recombination in Cameroon by Use of Nonspecific Next-Generation Sequencing

The importance of recombination in the evolution and genetic diversity of the hepatitis C virus (HCV) is currently uncertain. Only a small number of intergenotypic recombinants have been identified so far, and each has core and envelope genes classified as belonging to genotype 2. Here, we investigated two putative genotype 4/1 recombinants from southern Cameroon using a number of approaches, including standard Sanger sequencing, genotype-specific PCR amplification, and non-HCV-specific Illumina RNA sequencing (RNA-seq). Recombination between genotypes 1 and 4 was confirmed in both samples, and the parental lineages of each recombinant belong to HCV subtypes that are cocirculating at a high prevalence in Cameroon. Using the RNA-seq approach, we obtained a complete genome for one sample, which contained a recombination breakpoint at the E2/P7 gene junction. We developed and applied a new method, called Deep SimPlot, which can be used to visualize and identify viral recombination directly from the short sequence reads created by next-generation sequencing in conjunction with a consensus sequence.

Laboratory diagnosis of noroviruses: present and future

Norovirus is an important cause of gastroenteritis outbreaks globally and the most prevalent cause of sporadic gastroenteritis in many regions. Rapid and accurate identification of causative viral agents is critical for outbreak investigation, disease surveillance, and management. Because norovirus is not cultivable and has a highly diversified and variable genome, it is difficult to develop diagnostic assays. Detection methods have evolved from electron microscopy to conventional end-point reverse transcription polymerase chain reaction (RT-PCR), immunoassay, real-time RT-PCR, other molecular technologies, and nanotechnology array-based assays. The status and features of various testing methods are summarized in this review.

Hepatitis C

Hepatitis C virus (HCV) infection is a major health problem worldwide. The effects of chronic infection include cirrhosis, end-stage liver disease, and hepatocellular carcinoma. As a result of shared routes of transmission, co-infection with HIV is a substantial problem, and individuals infected with both viruses have poorer outcomes than do peers infected with one virus. No effective vaccine exists, although persistent HCV infection is potentially curable. The standard of care has been subcutaneous interferon alfa and oral ribavirin for 24-72 weeks. This treatment results in a sustained virological response in around 50% of individuals, and is complicated by clinically significant adverse events. In the past 10 years, advances in HCV cell culture have enabled an improved understanding of HCV virology, which has led to development of many new direct-acting antiviral drugs that target key components of virus replication. These direct-acting drugs allow for simplified and shortened treatments for HCV that can be given as oral regimens with increased tolerability and efficacy than interferon and ribavirin. Remaining obstacles include access to appropriate care and treatment, and development of a vaccine.

A universal protocol to generate consensus level genome sequences for foot-and-mouth disease virus and other positive-sense polyadenylated RNA viruses using the Illumina MiSeq

BACKGROUND

Next-Generation Sequencing (NGS) is revolutionizing molecular epidemiology by providing new approaches to undertake whole genome sequencing (WGS) in diagnostic settings for a variety of human and veterinary pathogens. Previous sequencing protocols have been subject to biases such as those encountered during PCR amplification and cell culture, or are restricted by the need for large quantities of starting material. We describe here a simple and robust methodology for the generation of whole genome sequences on the Illumina MiSeq. This protocol is specific for foot-and-mouth disease virus (FMDV) or other polyadenylated RNA viruses and circumvents both the use of PCR and the requirement for large amounts of initial template.

RESULTS

The protocol was successfully validated using five FMDV positive clinical samples from the 2001 epidemic in the United Kingdom, as well as a panel of representative viruses from all seven serotypes. In addition, this protocol was successfully used to recover 94% of an FMDV genome that had previously been identified as cell culture negative. Genome sequences from three other non-FMDV polyadenylated RNA viruses (EMCV, ERAV, VESV) were also obtained with minor protocol amendments. We calculated that a minimum coverage depth of 22 reads was required to produce an accurate consensus sequence for FMDV O. This was achieved in 5 FMDV/O/UKG isolates and the type O FMDV from the serotype panel with the exception of the 5' genomic termini and area immediately flanking the poly(C) region.

CONCLUSIONS

We have developed a universal WGS method for FMDV and other polyadenylated RNA viruses. This method works successfully from a limited quantity of starting material and eliminates the requirement for genome-specific PCR amplification. This protocol has the potential to generate consensus-level sequences within a routine high-throughput diagnostic environment.

Occurrence of norovirus and hepatitis A virus in wild mussels collected from the Baltic Sea

The aim of the study was to define the occurrence of human noroviruses of genogroup I and II (NoV GI and NoV GII) and hepatitis A virus (HAV) in the Baltic Sea mussels. The shellfish samples were taken at the sampling sites located on the Polish coast. In total, 120 shellfish were tested as pooled samples using RT-PCR and hybridisation with virus specific probes. NoV GI was detected in 22 (18.3%), NoV GII in 28 (23.3%), and HAV in 9 (7.5%) of the shellfish. The nucleotide sequence analysis of the detected NoV GII strains showed a 97.3-99.3% similarity to GII.4 virus strain. This is the first report describing the NoV and HAV occurrence in wild Baltic mussels and their possible role as bioindicators of seawater contamination with human enteric viruses.

Detection of circulating norovirus genotypes: hitting a moving target

BACKGROUND

Although national surveillance programs are in place to monitor norovirus epidemiology, the emergence of new strains and the genetic diversity among genotypes can be challenging for clinical laboratories. This study evaluated the analytical and clinical performance characteristics of one real-time RT-PCR and two end-point RT-PCRs commonly used in microbiology laboratories.

METHODS

Lower limit of detection (LoD) was determined using 10-fold dilutions of noroviruses belonging to different genotypes. The clinical performance of the real-time and end-point RT-PCRs was assessed in parallel using nucleic acids extracted from 186 stool specimens.

RESULTS

The real-time RT-PCR was highly sensitive and specific for the detection of norovirus genotypes that are currently circulating in Canada. In contrast, the two end-point RT-PCRs displayed poor analytical sensitivity or complete failure to detect certain norovirus genotypes, which was correlated to sequence mismatches in the primer-binding sites. In an attempt to improve norovirus detection with the end-point RT-PCRs, both assays were processed concurrently and detection from either assay was considered a positive result. Concurrent testing resulted in only a modest increase in clinical sensitivity (75.0%) compared to each assay alone (62.5% and 71.9%). However, the false positivity rate increased from 1.98% and 3.36% for the assays alone to 5.47% with concurrent testing.

CONCLUSIONS

This study emphasizes the benefits of a real-time method and provides support for routine surveillance to monitor norovirus epidemiology and ongoing proficiency testing to ensure detection of circulating norovirus genotypes.

Comprehensive human virus screening using high-throughput sequencing with a user-friendly representation of bioinformatics analysis: a pilot study

High-throughput sequencing (HTS) provides the means to analyze clinical specimens in unprecedented molecular detail. While this technology has been successfully applied to virus discovery and other related areas of research, HTS methodology has yet to be exploited for use in a clinical setting for routine diagnostics. Here, a bioinformatics pipeline (ezVIR) was designed to process HTS data from any of the standard platforms and to evaluate the entire spectrum of known human viruses at once, providing results that are easy to interpret and customizable. The pipeline works by identifying the most likely viruses present in the specimen given the sequencing data. Additionally, ezVIR can generate optional reports for strain typing, can create genome coverage histograms, and can perform cross-contamination analysis for specimens prepared in series. In this pilot study, the pipeline was challenged using HTS data from 20 clinical specimens representative of those most often collected and analyzed in daily practice. The specimens (5 cerebrospinal fluid, 7 bronchoalveolar lavage fluid, 5 plasma, 2 serum, and 1 nasopharyngeal aspirate) were originally found to be positive for a diverse range of DNA or RNA viruses by routine molecular diagnostics. The ezVIR pipeline correctly identified 14 of 14 specimens containing viruses with genomes of <40,000 bp, and 4 of 6 specimens positive for large-genome viruses. Although further validation is needed to evaluate sensitivity and to define detection cutoffs, results obtained in this pilot study indicate that the overall detection success rate, coupled with the ease of interpreting the analysis reports, makes it worth considering using HTS for clinical diagnostics.

Identification of novel bovine group A rotavirus G15P[14] strain from epizootic diarrhea of adult cows by de novo sequencing using a next-generation sequencer

There are few reports describing diarrhea of adult cattle caused by group A rotaviruses. Here, we report the identification of a novel bovine group A rotavirus from diarrhea of adult cows. A group A rotavirus was detected from an epizootic outbreak of diarrhea in adult cows with a decrease in milk production in Japan in 2013. The comprehensive genomic analyses from fecal samples by viral metagenomics using a next-generation sequencer revealed that it had an unreported genotype combination G15P[14]. The genome constellation of this strain, namely, RVA/Cow-wt/JPN/Tottori-SG/2013/G15P[14] was G15-P[14]-I2-R2-C2-M2-A3-N2-T6-E2-H3 representing VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5, respectively. Each gene segment of Tottori-SG was most closely related to Japanese bovine group A rotaviruses suggesting that Tottori-SG might have derived from multiple reassortment events from group A rotavirus strains circulating among Japanese cattle. No other diarrhea pathogen of adult cattle was detected by routine diagnosis and metagenomics. Viral metagenomics, using a next-generation sequencer, is useful to characterize group A rotaviruses from fecal samples and offers unbiased comprehensive investigations of pathogen.

Diagnostic metagenomics: potential applications to bacterial, viral and parasitic infections

The term ‘shotgun metagenomics’ is applied to the direct sequencing of DNA extracted from a sample without culture or target-specific amplification or capture. In diagnostic metagenomics, this approach is applied to clinical samples in the hope of detecting and characterizing pathogens. Here, I provide a conceptual overview, before reviewing several recent promising proof-of-principle applications of metagenomics in virus discovery, analysis of outbreaks and detection of pathogens in contemporary and historical samples. I also evaluate future prospects for diagnostic metagenomics in the light of relentless improvements in sequencing technologies.

Whole genome sequencing and de novo assembly identifies Sydney-like variant noroviruses and recombinants during the winter 2012/2013 outbreak in England

BACKGROUND

Norovirus is the commonest cause of epidemic gastroenteritis among people of all ages. Outbreaks frequently occur in hospitals and the community, costing the UK an estimated £110 m per annum. An evolutionary explanation for periodic increases in norovirus cases, despite some host-specific post immunity is currently limited to the identification of obvious recombinants. Our understanding could be significantly enhanced by full length genome sequences for large numbers of intensively sampled viruses, which would also assist control and vaccine design. Our objective is to develop rapid, high-throughput, end-to-end methods yielding complete norovirus genome sequences. We apply these methods to recent English outbreaks, placing them in the wider context of the international norovirus epidemic of winter 2012.

METHOD

Norovirus sequences were generated from 28 unique clinical samples by Illumina RNA sequencing (RNA-Seq) of total faecal RNA. A range of de novo sequence assemblers were attempted. The best assembler was identified by validation against three replicate samples and two norovirus qPCR negative samples, together with an additional 20 sequences determined by PCR and fractional capillary sequencing. Phylogenetic methods were used to reconstruct evolutionary relationships from the whole genome sequences.

RESULTS

Full length norovirus genomes were generated from 23/28 samples. 5/28 partial norovirus genomes were associated with low viral copy numbers. The de novo assembled sequences differed from sequences determined by capillary sequencing by <0.003%. Intra-host nucleotide sequence diversity was rare, but detectable by mapping short sequence reads onto its de novo assembled consensus. Genomes similar to the Sydney 2012 strain caused 78% (18/23) of cases, consistent with its previously documented association with the winter 2012 global outbreak. Interestingly, phylogenetic analysis and recombination detection analysis of the consensus sequences identified two related viruses as recombinants, containing sequences in prior circulation to Sydney 2012 in open reading frame (ORF) 2.

CONCLUSION

Our approach facilitates the rapid determination of complete norovirus genomes. This method provides high resolution of full norovirus genomes which, when coupled with detailed epidemiology, may improve the understanding of evolution and control of this important healthcare-associated pathogen.