Error baseline rates of five sample preparation methods used to characterize RNA virus populations

Investigating the effect of ribavirin treatment on genetic mutations in Crimean-Congo haemorrhagic fever virus (CCHFV) through next-generation sequencing

Crimean-Congo haemorrhagic fever (CCHF) is the most widespread tick-borne viral haemorrhagic fever affecting humans, and yet a licensed drug against the virus (CCHFV) is still not available. While several studies have suggested the efficacy of ribavirin against CCHFV, current literature remains inconclusive. In this study, we have utilised next-generation sequencing to investigate the mutagenic effect of ribavirin on the CCHFV genome during clinical disease. Samples collected from CCHF patients receiving ribavirin treatment or supportive care only at Sivas Cumhuriyet University Hospital, Turkey, were analysed. By comparing the frequency of mutations in each group, we found little evidence of an overall mutagenic effect. This suggests that ribavirin, administered at the acute stages of CCHFV infection (at the World Health Organization-recommended dose) is unable to induce lethal mutagenesis that would cause an extinction event in the CCHFV population and reduce viremia.

On-Demand Patient-Specific Phenotype-to-Genotype Ebola Virus Characterization

Biosafety, biosecurity, logistical, political, and technical considerations can delay or prevent the wide dissemination of source material containing viable virus from the geographic origin of an outbreak to laboratories involved in developing medical countermeasures (MCMs). However, once virus genome sequence information is available from clinical samples, reverse-genetics systems can be used to generate virus stocks de novo to initiate MCM development. In this study, we developed a reverse-genetics system for natural isolates of Ebola virus (EBOV) variants Makona, Tumba, and Ituri, which have been challenging to obtain. These systems were generated starting solely with in silico genome sequence information and have been used successfully to produce recombinant stocks of each of the viruses for use in MCM testing. The antiviral activity of MCMs targeting viral entry varied depending on the recombinant virus isolate used. Collectively, selecting and synthetically engineering emerging EBOV variants and demonstrating their efficacy against available MCMs will be crucial for answering pressing public health and biosecurity concerns during Ebola disease (EBOD) outbreaks.

Next generation sequencing based in-house HIV genotyping method: validation report

Background

HIV genotyping has had a significant impact on the care and treatment of HIV/AIDS. At a clinical level, the test guides physicians on the choice of treatment regimens. At the surveillance level, it informs policy on consolidated treatment guidelines and microbial resistance control strategies. Until recently, the conventional test has utilized the Sanger sequencing (SS) method. Unlike Next Generation Sequencing (NGS), SS is limited by low data throughput and the inability of detecting low abundant drug-resistant variants. NGS can improve sensitivity and quantitatively identify low-abundance variants; in addition, it has the potential to improve efficiency as well as lowering costs when samples are batched. Despite the NGS benefits, its utilization in clinical drug resistance profiling is faced with mixed reactions. These are largely based on a lack of a consensus regarding the quality control strategy. Nonetheless, transitional views suggest validating the method against the gold-standard SS. Therefore, we present a validation report of an NGS-based in-house HIV genotyping method against the SS method in Uganda.

Results

Since there were no established proficiency test panels for NGS-based HIV genotyping, 15 clinical plasma samples for routine care were utilized. The use of clinical samples allowed for accuracy and precision studies. The workflow involved four main steps; viral RNA extraction, targeted amplicon generation, amplicon sequencing and data analysis. Accuracy of 98% with an average percentage error of 3% was reported for the NGS based assay against the SS platform demonstrating similar performance. The coefficient of variation (CV) findings for both the inter-run and inter-personnel precision showed no variability (CV ≤ 0%) at the relative abundance of ≥ 20%. For both inter-run and inter-personnel, a variation that affected the precision was observed at 1% frequency. Overall, for all the frequencies, CV registered a small range of (0–2%).

Conclusion

The NGS-based in-house HIV genotyping method fulfilled the minimum requirements that support its utilization for drug resistance profiling in a clinical setting of a low-income country. For more inclusive quality control studies, well-characterized wet panels need to be established.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12981-021-00390-8.

Molecular analysis of the 2012 Bundibugyo virus disease outbreak

Bundibugyo virus (BDBV) is one of four ebolaviruses known to cause disease in humans. Bundibugyo virus disease (BVD) outbreaks occurred in 2007-2008 in Bundibugyo District, Uganda, and in 2012 in Isiro, Province Orientale, Democratic Republic of the Congo. The 2012 BVD outbreak resulted in 38 laboratory-confirmed cases of human infection, 13 of whom died. However, only 4 BDBV specimens from the 2012 outbreak have been sequenced. Here, we provide BDBV sequences from seven additional patients. Analysis of the molecular epidemiology and evolutionary dynamics of the 2012 outbreak with these additional isolates challenges the current hypothesis that the outbreak was the result of a single spillover event. In addition, one patient record indicates that BDBV's initial emergence in Isiro occurred 50 days earlier than previously accepted. Collectively, this work demonstrates how retrospective sequencing can be used to elucidate outbreak origins and provide epidemiological contexts to a medically relevant pathogen.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Stability of SARS-CoV-2 phylogenies

The SARS-CoV-2 pandemic has led to unprecedented, nearly real-time genetic tracing due to the rapid community sequencing response. Researchers immediately leveraged these data to infer the evolutionary relationships among viral samples and to study key biological questions, including whether host viral genome editing and recombination are features of SARS-CoV-2 evolution. This global sequencing effort is inherently decentralized and must rely on data collected by many labs using a wide variety of molecular and bioinformatic techniques. There is thus a strong possibility that systematic errors associated with lab-or protocol-specific practices affect some sequences in the repositories. We find that some recurrent mutations in reported SARS-CoV-2 genome sequences have been observed predominantly or exclusively by single labs, co-localize with commonly used primer binding sites and are more likely to affect the protein-coding sequences than other similarly recurrent mutations. We show that their inclusion can affect phylogenetic inference on scales relevant to local lineage tracing, and make it appear as though there has been an excess of recurrent mutation or recombination among viral lineages. We suggest how samples can be screened and problematic variants removed, and we plan to regularly inform the scientific community with our updated results as more SARS-CoV-2 genome sequences are shared (https://virological.org/t/issues-with-sars-cov-2-sequencing-data/473 and https://virological.org/t/masking-strategies-for-sars-cov-2-alignments/480). We also develop tools for comparing and visualizing differences among very large phylogenies and we show that consistent clade- and tree-based comparisons can be made between phylogenies produced by different groups. These will facilitate evolutionary inferences and comparisons among phylogenies produced for a wide array of purposes. Building on the SARS-CoV-2 Genome Browser at UCSC, we present a toolkit to compare, analyze and combine SARS-CoV-2 phylogenies, find and remove potential sequencing errors and establish a widely shared, stable clade structure for a more accurate scientific inference and discourse.

Amplicon-Based, Next-Generation Sequencing Approaches to Characterize Single Nucleotide Polymorphisms of Orthohantavirus Species

Whole-genome sequencing (WGS) of viruses from patient or environmental samples can provide tremendous insight into the epidemiology, drug resistance or evolution of a virus. However, we face two common hurdles in obtaining robust sequence information; the low copy number of viral genomes in specimens and the error introduced by WGS techniques. To optimize detection and minimize error in WGS of hantaviruses, we tested four amplification approaches and different amplicon pooling methods for library preparation and examined these preparations using two sequencing platforms, Illumina MiSeq and Oxford Nanopore Technologies MinION. First, we tested and optimized primers used for whole segment PCR or one kilobase amplicon amplification for even coverage using RNA isolated from the supernatant of virus-infected cells. Once optimized we assessed two sources of total RNA, virus-infected cells and supernatant from the virus-infected cells, with four variations of primer pooling for amplicons, and six different amplification approaches. We show that 99-100% genome coverage was obtained using a one-step RT-PCR reaction with one forward and reverse primer. Using a two-step RT-PCR with three distinct tiling approaches for the three genomic segments (vRNAs), we optimized primer pooling approaches for PCR amplification to achieve a greater number of aligned reads, average depth of genome, and genome coverage. The single nucleotide polymorphisms identified from MiSeq and MinION sequencing suggested intrinsic mutation frequencies of ~10-5-10-7 per genome and 10-4-10-5 per genome, respectively. We noted no difference in the coverage or accuracy when comparing WGS results with amplicons amplified from RNA extracted from infected cells or supernatant of these infected cells. Our results show that high-throughput diagnostics requiring the identification of hantavirus species or strains can be performed using MiSeq or MinION using a one-step approach. However, the two-step MiSeq approach outperformed the MinION in coverage depth and accuracy, and hence would be superior for assessment of genomes for epidemiology or evolutionary questions using the methods developed herein.

A beginner's guide for FMDV quasispecies analysis: sub-consensus variant detection and haplotype reconstruction using next-generation sequencing

Deep sequencing of viral genomes is a powerful tool to study RNA virus complexity. However, the analysis of next-generation sequencing data might be challenging for researchers who have never approached the study of viral quasispecies by this methodology. In this work we present a suitable and affordable guide to explore the sub-consensus variability and to reconstruct viral quasispecies from Illumina sequencing data. The guide includes a complete analysis pipeline along with user-friendly descriptions of software and file formats. In addition, we assessed the feasibility of the workflow proposed by analyzing a set of foot-and-mouth disease viruses (FMDV) with different degrees of variability. This guide introduces the analysis of quasispecies of FMDV and other viruses through this kind of approach.

A Model for the Production of Regulatory Grade Viral Hemorrhagic Fever Exposure Stocks: From Field Surveillance to Advanced Characterization of SFTSV

Severe fever with thrombocytopenia syndrome virus (SFTSV) is an emerging human pathogen, endemic in areas of China, Japan, and the Korea (KOR). It is primarily transmitted through infected ticks and can cause a severe hemorrhagic fever disease with case fatality rates as high as 30%. Despite its high virulence and increasing prevalence, molecular and functional studies in situ are scarce due to the limited availability of high-titer SFTSV exposure stocks. During the course of field virologic surveillance in 2017, we detected SFTSV in ticks and in a symptomatic soldier in a KOR Army training area. SFTSV was isolated from the ticks producing a high-titer viral exposure stock. Through the use of advanced genomic tools, we present here a complete, in-depth characterization of this viral stock, including a comparison with both the virus in its arthropod source and in the human case, and an in vivo study of its pathogenicity. Thanks to this detailed characterization, this SFTSV viral exposure stock constitutes a quality biological tool for the study of this viral agent and for the development of medical countermeasures, fulfilling the requirements of the main regulatory agencies.

Development and Application of Performance Assessment Criteria for Next-Generation Sequencing-Based HIV Drug Resistance Assays

Next-generation sequencing (NGS)-based HIV drug resistance (HIVDR) assays outperform conventional Sanger sequencing in scalability, sensitivity, and quantitative detection of minority resistance variants. Thus far, HIVDR assays have been applied primarily in research but rarely in clinical settings. One main obstacle is the lack of standardized validation and performance evaluation systems that allow regulatory agencies to benchmark and accredit new assays for clinical use. By revisiting the existing principles for molecular assay validation, here we propose a new validation and performance evaluation system that helps to both qualitatively and quantitatively assess the performance of an NGS-based HIVDR assay. To accomplish this, we constructed a 70-specimen proficiency test panel that includes plasmid mixtures at known ratios, viral RNA from infectious clones, and anonymized clinical specimens. We developed assessment criteria and benchmarks for NGS-based HIVDR assays and used these to assess data from five separate MiSeq runs performed in two experienced HIVDR laboratories. This proposed platform may help to pave the way for the standardization of NGS HIVDR assay validation and performance evaluation strategies for accreditation and quality assurance purposes in both research and clinical settings.

External Quality Assessment for Next-Generation Sequencing-Based HIV Drug Resistance Testing: Unique Requirements and Challenges

Over the past decade, there has been an increase in the adoption of next generation sequencing (NGS) technologies for HIV drug resistance (HIVDR) testing. NGS far outweighs conventional Sanger sequencing as it has much higher throughput, lower cost when samples are batched and, most importantly, significantly higher sensitivities for variants present at low frequencies, which may have significant clinical implications. Despite the advantages of NGS, Sanger sequencing remains the gold standard for HIVDR testing, largely due to the lack of standardization of NGS-based HIVDR testing. One important aspect of standardization includes external quality assessment (EQA) strategies and programs. Current EQA for Sanger-based HIVDR testing includes proficiency testing where samples are sent to labs and the performance of the lab conducting such assays is evaluated. The current methods for Sanger-based EQA may not apply to NGS-based tests because of the fundamental differences in their technologies and outputs. Sanger-based genotyping reports drug resistance mutations (DRMs) data as dichotomous, whereas NGS-based HIVDR genotyping also reports DRMs as numerical data (percent abundance). Here we present an overview of the need to develop EQA for NGS-based HIVDR testing and some unique challenges that may be encountered.

Reverse Genetics of RNA Viruses: ISA-Based Approach to Control Viral Population Diversity without Modifying Virus Phenotype

Reverse genetic systems are essential for the study of RNA viruses. Infectious clones remain the most widely used systems to manipulate viral genomes. Recently, a new PCR-based method called ISA (infectious subgenomic amplicons) has been developed. This approach has resulted in greater genetic diversity of the viral populations than that observed using infectious clone technology. However, for some studies, generation of clonal viral populations is necessary. In this study, we used the tick-borne encephalitis virus as model to demonstrate that utilization of a very high-fidelity, DNA-dependent DNA polymerase during the PCR step of the ISA procedure gives the possibility to reduce the genetic diversity of viral populations. We also concluded that the fidelity of the polymerase is not the only factor influencing this diversity. Studying the impact of genotype modification on virus phenotype is a crucial step for the development of reverse genetic methods. Here, we also demonstrated that the utilization of different PCR polymerases did not affect the phenotype (replicative fitness in cellulo and virulence in vivo) compared to the initial ISA procedure and the use of an infectious clone. In conclusion, we provide here an approach to control the genetic diversity of RNA viruses without modifying their phenotype.

T-705 induces lethal mutagenesis in Ebola and Marburg populations in macaques

Nucleoside analogues (NA) disrupt RNA viral RNA-dependent RNA polymerase (RdRP) function and fidelity for multiple viral families. The mechanism of action (MOA) of T-705 has been attributed alternatively or concurrently to chain termination and lethal mutagenesis depending on the viral species during in vitro studies. In this study, we evaluated the effect of T-705 on the viral population in non-human primates (NHPs) after challenge with Ebola virus (EBOV) or Marburg virus (MARV) to identify the predominant in vivo MOA. We used common virological assays in conjunction with deep sequencing to characterize T-705 effects. T-705 exhibited antiviral activity that was associated with a reduction in specific infectivity and an accumulation of low frequency nucleotide variants in plasma samples collected day 7 post infection. Stranded analysis of deep sequencing data to identify chain termination demonstrated no change in the transcriptional gradient in negative stranded viral reads and minimal changes in positive stranded viral reads in T-705 treated animals, questioning as a MOA in vivo. These findings indicate that lethal mutagenesis is a MOA of T-705 that may serve as an indication of therapeutic activity of NAs for evaluation in clinical settings. This study expands our understanding of MOAs of these compounds for the Filovirus family and provides further evidence that lethal mutagenesis could be a preponderant MOA for this class of therapeutic compounds.

Prospective Cohort Study of Next-Generation Sequencing as a Diagnostic Modality for Unexplained Encephalitis in Children

BACKGROUND

Encephalitis is an inflammatory condition of the brain associated with long-term neurologic sequelae and even death in children. Although viruses are often implicated, an etiology is not identified in the majority of cases. Metagenomics-based next-generation sequencing (mNGS) is a high-throughput sequencing technique that can enhance the detection of novel or low-frequency pathogens.

METHODS

Hospitalized immunocompetent children aged 6 months to 18 years with encephalitis of unidentified etiology were eligible for enrollment. Demographic, historical, and clinical information was obtained, and residual blood and cerebrospinal fluid (CSF) samples were subjected to mNGS. Pathogens were identified by querying the sequence data against the NCBI GenBank database.

RESULTS

Twenty children were enrolled prospectively between 2013 and 2017. mNGS of CSF identified 7 nonhuman nucleic acid sequences of significant frequency in 6 patients, including that of Mycoplasma bovis, parvovirus B19, Neisseria meningitidis, and Balamuthia mandrillaris. mNGS also detected Cladophialophora species, tobacco mosaic virus, and human bocavirus, which were presumed to be contaminants or nonpathogenic organisms. One patient was found to have positive serology results for California encephalitis virus, but mNGS did not detect it. Patients for whom mNGS identified a diagnosis had a significantly higher CSF white blood cell count, a higher CSF protein concentration, and a lower CSF glucose level than patients for whom mNGS did not identify a diagnosis.

CONCLUSION

We describe here the results of a prospective cohort analysis to evaluate mNGS as a diagnostic tool for children with unexplained encephalitis. Although mNGS detected multiple nonpathogenic organisms, it also identified multiple pathogens successfully and was most useful in patients with a CSF abnormality.

Growth-Adaptive Mutations in the Ebola Virus Makona Glycoprotein Alter Different Steps in the Virus Entry Pathway

The Zaire ebolavirus (EBOV) glycoprotein (GP) is cleaved into two subunits (GP1 and GP2) that are both required for virus attachment and entry into cells. Sequence changes in the GP have been proposed to increase pathogenesis and to alter virus growth properties. Mutations in GP acquired during EBOV tissue culture passage have also been reported to change virus growth properties. Here, we report the isolation of six amino acid mutations in EBOV GP that spontaneously appeared during recovery and passage of an EBOV-Makona GP-pseudotyped vesicular stomatitis virus (VSV), two of which also occur during passage of EBOV clinical isolates in tissue culture. Each of the six mutations resulted in increased virus growth in monkey and human cell lines. All mutations are located in the GP2 fusion subunit and increase entry kinetics of EBOV virus-like particles (VLPs). The gain-of-entry function mapped to two mechanistic phenotypes. Mutations in heptad repeat 1 (HR1) decreased the requirement for cathepsin B activity for viral infection. Mutations directly within the fusion loop increased entry kinetics without altering the cathepsin B dependence. Several mutations in the fusion loop were substitutions of residues present in other ebolavirus glycoproteins, illustrating the evolutionary paths for maintaining an optimally functioning fusion loop under selection pressure.IMPORTANCEZaire ebolavirus (EBOV) is the causative agent of the highly lethal Ebola virus disease and poses a significant threat to the global health community. Approved antivirals against EBOV are lacking; however, promising therapies targeting the EBOV glycoprotein are being developed. Efficacy testing of these candidate therapeutics relies on EBOV laboratory stocks, which when grown in tissue culture may acquire mutations in the glycoprotein. These mutations can produce inaccurate results in therapeutic testing. Until recently, distinguishing between tissue culture mutations and naturally occurring polymorphisms in EBOV GP was difficult in the absence of consensus clinical GP sequences. Here, we utilize recombinant VSV (rVSV) pseudotyped with the consensus clinical EBOV Makona GP to identify several mutations that have emerged or have potential to emerge in EBOV GP during tissue culture passage. Identifying these mutations informs the EBOV research community as to which mutations may arise during preparation of laboratory virus stocks.

Stochastic processes constrain the within and between host evolution of influenza virus

The evolutionary dynamics of influenza virus ultimately derive from processes that take place within and between infected individuals. Here we define influenza virus dynamics in human hosts through sequencing of 249 specimens from 200 individuals collected over 6290 person-seasons of observation. Because these viruses were collected from individuals in a prospective community-based cohort, they are broadly representative of natural infections with seasonal viruses. Consistent with a neutral model of evolution, sequence data from 49 serially sampled individuals illustrated the dynamic turnover of synonymous and nonsynonymous single nucleotide variants and provided little evidence for positive selection of antigenic variants. We also identified 43 genetically-validated transmission pairs in this cohort. Maximum likelihood optimization of multiple transmission models estimated an effective transmission bottleneck of 1-2 genomes. Our data suggest that positive selection is inefficient at the level of the individual host and that stochastic processes dominate the host-level evolution of influenza viruses.

Within-Host Evolution of Human Influenza Virus

The rapid global evolution of influenza virus begins with mutations that arise de novo in individual infections, but little is known about how evolution occurs within hosts. We review recent progress in understanding how and why influenza viruses evolve within human hosts. Advances in deep sequencing make it possible to measure within-host genetic diversity in both acute and chronic influenza infections. Factors like antigenic selection, antiviral treatment, tissue specificity, spatial structure, and multiplicity of infection may affect how influenza viruses evolve within human hosts. Studies of within-host evolution can contribute to our understanding of the evolutionary and epidemiological factors that shape influenza virus's global evolution.