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

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.

Whole genome sequencing of hepatitis B virus using tiled amplicon (HEPTILE) and probe based enrichment on Illumina and Nanopore platforms

Hepatitis B virus (HBV) whole genome sequencing (WGS) is currently limited as the DNA viral loads (VL) of many clinical samples are below the threshold required to generate full genomes using current sequencing methods. We developed two pan-genotypic viral enrichment methods, using probe-based capture and tiled amplicon PCR (HEP-TILE) for HBV WGS. We demonstrate using mock samples that both enrichment methods are pan-genotypic (genotypes A-J). Using clinical samples, we demonstrate that HEP-TILE amplification successfully amplifies full genomes at the lowest HBV VL tested (30 IU/ml), and the PCR products can be sequenced using both Nanopore and Illumina platforms. Probe-based capture with Illumina sequencing required VL > 300,000 IU/ml to generate full length HBV genomes. The capture-Illumina and HEP-TILE-Nanopore pipelines had consensus sequencing accuracy of 100% in mock samples with known DNA sequences. Together, these protocols will facilitate the generation of HBV sequence data, enabling a more accurate and representative picture of HBV molecular epidemiology, cast light on persistence and pathogenesis, and enhance understanding of the outcomes of infection and its treatment.

Castanet: a pipeline for rapid analysis of targeted multi-pathogen genomic data

MOTIVATION

Target enrichment strategies generate genomic data from multiple pathogens in a single process, greatly improving sensitivity over metagenomic sequencing and enabling cost-effective, high-throughput surveillance and clinical applications. However, uptake by research and clinical laboratories is constrained by an absence of computational tools that are specifically designed for the analysis of multi-pathogen enrichment sequence data. Here we present an analysis pipeline, Castanet, for use with multi-pathogen enrichment sequencing data. Castanet is designed to work with short-read data produced by existing targeted enrichment strategies, but can be readily deployed on any BAM file generated by another methodology. Also included are an optional graphical interface and installer script.

RESULTS

In addition to genome reconstruction, Castanet reports method-specific metrics that enable quantification of capture efficiency, estimation of pathogen load, differentiation of low-level positives from contamination, and assessment of sequencing quality. Castanet can be used as a traditional end-to-end pipeline for consensus generation, but its strength lies in the ability to process a flexible, pre-defined set of pathogens of interest directly from multi-pathogen enrichment experiments. In our tests, Castanet consensus sequences were accurate reconstructions of reference sequences, including in instances where multiple strains of the same pathogen were present. Castanet performs effectively on standard computers and can process the entire output of a 96-sample enrichment sequencing run (50M reads) using a single batch process command, in $<$2 h.

AVAILABILITY AND IMPLEMENTATION

Source code freely available under GPL-3 license at https://github.com/MultipathogenGenomics/castanet, implemented in Python 3.10 and supported in Ubuntu Linux 22.04. The data underlying this article are available in Europe Nucleotide Archives, at https://www.ebi.ac.uk/ena/browser/view/PRJEB77004.

A proof of concept for a targeted enrichment approach to the simultaneous detection and characterization of rickettsial pathogens from clinical specimens

Infection with either Rickettsia prowazekii or Orientia tsutsugamushi is common, yet diagnostic capabilities are limited due to the short window for positive identification. Until now, although targeted enrichment had been applied to increase sensitivity of sequencing-based detection for various microorganisms, it had not been applied to sequencing of R. prowazekii in clinical samples. Additionally, hybridization-based targeted enrichment strategies had only scarcely been applied to qPCR of any pathogens in clinical samples. Therefore, we tested a targeted enrichment technique as a proof of concept and found that it dramatically reduced the limits of detection of these organisms by both qPCR and high throughput sequencing. The enrichment methodology was first tested in contrived clinical samples with known spiked-in concentrations of R. prowazekii and O. tsutsugamushi DNA. This method was also evaluated using clinical samples, resulting in the simultaneous identification and characterization of O. tsutsugamushi directly from clinical specimens taken from sepsis patients. We demonstrated that the targeted enrichment technique is helpful by lowering the limit of detection, not only when applied to sequencing, but also when applied to qPCR, suggesting the technique could be applied more broadly to include other assays and/or microbes for which there are limited diagnostic or detection modalities.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

Targeted metagenomics reveals association between severity and pathogen co-detection in infants with respiratory syncytial virus

Respiratory syncytial virus (RSV) is the leading cause of hospitalisation for respiratory infection in young children. RSV disease severity is known to be age-dependent and highest in young infants, but other correlates of severity, particularly the presence of additional respiratory pathogens, are less well understood. In this study, nasopharyngeal swabs were collected from two cohorts of RSV-positive infants <12 months in Spain, the UK, and the Netherlands during 2017-20. We show, using targeted metagenomic sequencing of >100 pathogens, including all common respiratory viruses and bacteria, from samples collected from 433 infants, that burden of additional viruses is common (111/433, 26%) but only modestly correlates with RSV disease severity. In contrast, there is strong evidence in both cohorts and across age groups that presence of Haemophilus bacteria (194/433, 45%) is associated with higher severity, including much higher rates of hospitalisation (odds ratio 4.25, 95% CI 2.03-9.31). There is no evidence for association between higher severity and other detected bacteria, and no difference in severity between RSV genotypes. Our findings reveal the genomic diversity of additional pathogens during RSV infection in infants, and provide an evidence base for future causal investigations of the impact of co-infection on RSV disease severity.

A robust, scalable, and cost-efficient approach to whole genome sequencing of RSV directly from clinical samples

Respiratory syncytial virus (RSV) is a leading cause of acute lower respiratory tract infections causing significant morbidity and mortality among children and the elderly; two RSV vaccines and a monoclonal antibody have recently been approved. Thus, there is an increasing need for a detailed and continuous genomic surveillance of RSV circulating in resource-rich and resource-limited settings worldwide. However, robust, cost-effective methods for whole genome sequencing of RSV from clinical samples that are amenable to high-throughput are still scarce. We developed Next-RSV-SEQ, an experimental and computational pipeline to generate whole genome sequences of historic and current RSV genotypes by in-solution hybridization capture-based next generation sequencing. We optimized this workflow by automating library preparation and pooling libraries prior to enrichment in order to reduce hands-on time and cost, thereby augmenting scalability. Next-RSV-SEQ yielded near-complete to complete genome sequences for 98% of specimens with Cp values ≤31, at median on-target reads >93%, and mean coverage depths between ~1,000 and >5,000, depending on viral load. Whole genomes were successfully recovered from samples with viral loads as low as 230 copies per microliter RNA. We demonstrate that the method can be expanded to other respiratory viruses like parainfluenza virus and human metapneumovirus. Next-RSV-SEQ produces high-quality RSV genomes directly from culture isolates and, more importantly, clinical specimens of all genotypes in circulation. It is cost-efficient, scalable, and can be extended to other respiratory viruses, thereby opening new perspectives for a future effective and broad genomic surveillance of respiratory viruses.

IMPORTANCE

Respiratory syncytial virus (RSV) is a leading cause of severe acute respiratory tract infections in children and the elderly, and its prevention has become an increasing priority. Recently, vaccines and a long-acting monoclonal antibody to protect effectively against severe disease have been approved for the first time. Hence, there is an urgent need for genomic surveillance of RSV at the global scale to monitor virus evolution, especially with an eye toward immune evasion. However, robust, cost-effective methods for RSV whole genome sequencing that are suitable for high-throughput of clinical samples are currently scarce. Therefore, we have developed Next-RSV-SEQ, an experimental and computational pipeline that produces reliably high-quality RSV genomes directly from clinical specimens and isolates.

High-performance enrichment-based genome sequencing to support the investigation of hepatitis A virus outbreaks

IMPORTANCE

This proof-of-concept study introduces a hybrid capture oligo panel for whole-genome sequencing of all six human pathogenic hepatitis A virus (HAV) subgenotypes, exhibiting a higher sensitivity than some conventional genotyping assays. The ability of hybrid capture to enrich multiple targets allows for a single, streamlined workflow, thus facilitating the potential harmonization of molecular surveillance of HAV with other enteric viruses. Even challenging sample matrices can be accommodated, making them suitable for broad implementation in clinical and public health laboratories. This innovative approach has significant implications for enhancing multijurisdictional outbreak investigations as well as our understanding of the global diversity and transmission dynamics of HAV.

Longitudinal population-level HIV epidemiologic and genomic surveillance highlights growing gender disparity of HIV transmission in Uganda

HIV incidence in eastern and southern Africa has historically been concentrated among girls and women aged 15-24 years. As new cases decline with HIV interventions, population-level infection dynamics may shift by age and gender. Here, we integrated population-based surveillance of 38,749 participants in the Rakai Community Cohort Study and longitudinal deep-sequence viral phylogenetics to assess how HIV incidence and population groups driving transmission have changed from 2003 to 2018 in Uganda. We observed 1,117 individuals in the incidence cohort and 1,978 individuals in the transmission cohort. HIV viral suppression increased more rapidly in women than men, however incidence declined more slowly in women than men. We found that age-specific transmission flows shifted: whereas HIV transmission to girls and women (aged 15-24 years) from older men declined by about one-third, transmission to women (aged 25-34 years) from men that were 0-6 years older increased by half in 2003 to 2018. Based on changes in transmission flows, we estimated that closing the gender gap in viral suppression could have reduced HIV incidence in women by half in 2018. This study suggests that HIV programmes to increase HIV suppression in men are critical to reduce incidence in women, close gender gaps in infection burden and improve men's health in Africa.

Phylogenetic Analysis of Hepatitis C Virus Infections in a Large Belgian Cohort Using Next-Generation Sequencing of Full-Length Genomes

The hepatitis C virus (HCV) epidemic in Western countries is primarily perpetuated by the sub-populations of men who have sex with men (MSM) and people who inject drugs (PWID). Understanding the dynamics of transmission in these communities is crucial for removing the remaining hurdles towards HCV elimination. We sequenced 269 annotated HCV plasma samples using probe enrichment and next-generation sequencing, obtaining 224 open reading frames of HCV (OR497849-OR498072). Maximum likelihood phylogenies were generated on the four most prevalent subtypes in this study (HCV1a, 1b, 3a, 4d) with a subsequent transmission cluster analysis. The highest rate of clustering was observed for HCV4d samples (13/17 (76.47%)). The second highest rate of clustering was observed in HCV1a samples (42/78 (53.85%)) with significant association with HIV-positive MSM. HCV1b and HCV3a had very low rates of clustering (2/83 (2.41%) and (0/29)). The spread of the prevalent subtype HCV1b appears to have been largely curtailed, and we demonstrate the onwards transmission of HCV1a and HCV4d in the HIV-positive MSM population across municipal borders. More systematic data collection and sequencing is needed to allow a better understanding of the HCV transmission among the community of PWID and overcome the remaining barriers for HCV elimination in Belgium.

Longitudinal population-level HIV epidemiologic and genomic surveillance highlights growing gender disparity of HIV transmission in Uganda

HIV incidence in eastern and southern Africa has historically been concentrated among girls and women aged 15-24 years. As new cases decline with HIV interventions, population-level infection dynamics may shift by age and gender. Here, we integrated population-based surveillance of 38,749 participants in the Rakai Community Cohort Study and longitudinal deep sequence viral phylogenetics to assess how HIV incidence and population groups driving transmission have changed from 2003 to 2018 in Uganda. We observed 1,117 individuals in the incidence cohort and 1,978 individuals in the transmission cohort. HIV viral suppression increased more rapidly in women than men, however incidence declined more slowly in women than men. We found that age-specific transmission flows shifted, while HIV transmission to girls and women (aged 15-24 years) from older men declined by about one third, transmission to women (aged 25-34 years) from men that were 0-6 years older increased by half in 2003 to 2018. Based on changes in transmission flows, we estimated that closing the gender gap in viral suppression could have reduced HIV incidence in women by half in 2018. This study suggests that HIV programs to increase HIV suppression in men are critical to reduce incidence in women, close gender gaps in infection burden and improve men's health in Africa.

An enrichment protocol and analysis pipeline for long read sequencing of the hepatitis B virus transcriptome

Hepatitis B virus (HBV) is one of the smallest human DNA viruses and its 3.2 Kb genome encodes multiple overlapping open reading frames, making its viral transcriptome challenging to dissect. Previous studies have combined quantitative PCR and Next Generation Sequencing to identify viral transcripts and splice junctions, however the fragmentation and selective amplification used in short read sequencing precludes the resolution of full length RNAs. Our study coupled an oligonucleotide enrichment protocol with state-of-the-art long read sequencing (PacBio) to identify the repertoire of HBV RNAs. This methodology provides sequencing libraries where up to 25 % of reads are of viral origin and enable the identification of canonical (unspliced), non-canonical (spliced) and chimeric viral-human transcripts. Sequencing RNA isolated from de novo HBV infected cells or those transfected with 1.3 × overlength HBV genomes allowed us to assess the viral transcriptome and to annotate 5' truncations and polyadenylation profiles. The two HBV model systems showed an excellent agreement in the pattern of major viral RNAs, however differences were noted in the abundance of spliced transcripts. Viral-host chimeric transcripts were identified and more commonly found in the transfected cells. Enrichment capture and PacBio sequencing allows the assignment of canonical and non-canonical HBV RNAs using an open-source analysis pipeline that enables the accurate mapping of the HBV transcriptome.

Targeted genomic sequencing with probe capture for discovery and surveillance of coronaviruses in bats

Public health emergencies like SARS, MERS, and COVID-19 have prioritized surveillance of zoonotic coronaviruses, resulting in extensive genomic characterization of coronavirus diversity in bats. Sequencing viral genomes directly from animal specimens remains a laboratory challenge, however, and most bat coronaviruses have been characterized solely by PCR amplification of small regions from the best-conserved gene. This has resulted in limited phylogenetic resolution and left viral genetic factors relevant to threat assessment undescribed. In this study, we evaluated whether a technique called hybridization probe capture can achieve more extensive genome recovery from surveillance specimens. Using a custom panel of 20,000 probes, we captured and sequenced coronavirus genomic material in 21 swab specimens collected from bats in the Democratic Republic of the Congo. For 15 of these specimens, probe capture recovered more genome sequence than had been previously generated with standard amplicon sequencing protocols, providing a median 6.1-fold improvement (ranging up to 69.1-fold). Probe capture data also identified five novel alpha- and betacoronaviruses in these specimens, and their full genomes were recovered with additional deep sequencing. Based on these experiences, we discuss how probe capture could be effectively operationalized alongside other sequencing technologies for high-throughput, genomics-based discovery and surveillance of bat coronaviruses.

ProbeTools: designing hybridization probes for targeted genomic sequencing of diverse and hypervariable viral taxa

BACKGROUND

Sequencing viruses in many specimens is hindered by excessive background material from hosts, microbiota, and environmental organisms. Consequently, enrichment of target genomic material is necessary for practical high-throughput viral genome sequencing. Hybridization probes are widely used for enrichment in many fields, but their application to viral sequencing faces a major obstacle: it is difficult to design panels of probe oligo sequences that broadly target many viral taxa due to their rapid evolution, extensive diversity, and genetic hypervariability. To address this challenge, we created ProbeTools, a package of bioinformatic tools for generating effective viral capture panels, and for assessing coverage of target sequences by probe panel designs in silico. In this study, we validated ProbeTools by designing a panel of 3600 probes for subtyping the hypervariable haemagglutinin (HA) and neuraminidase (NA) genome segments of avian-origin influenza A viruses (AIVs). Using in silico assessment of AIV reference sequences and in vitro capture on egg-cultured viral isolates, we demonstrated effective performance by our custom AIV panel and ProbeTools' suitability for challenging viral probe design applications.

RESULTS

Based on ProbeTool's in silico analysis, our panel provided broadly inclusive coverage of 14,772 HA and 11,967 NA reference sequences. For each reference sequence, we calculated the percentage of nucleotide positions covered by our panel in silico; 90% of HA and NA references sequences had at least 90.8 and 95.1% of their nucleotide positions covered respectively. We also observed effective in vitro capture on a representative collection of 23 egg-cultured AIVs that included isolates from wild birds, poultry, and humans and representatives from all HA and NA subtypes. Forty-two of forty-six HA and NA segments had over 98.3% of their nucleotide positions significantly enriched by our custom panel. These in vitro results were further used to validate ProbeTools' in silico coverage assessment algorithm; 89.2% of in silico predictions were concordant with in vitro results.

CONCLUSIONS

ProbeTools generated an effective panel for subtyping AIVs that can be deployed for genomic surveillance, outbreak prevention, and pandemic preparedness. Effective probe design against hypervariable AIV targets also validated ProbeTools' design and coverage assessment algorithms, demonstrating their suitability for other challenging viral capture applications.

Implementation of a controlled human infection model for evaluation of HCV vaccine candidates

Hepatitis C virus (HCV) remains a major global health concern. Directly acting antiviral (DAA) drugs have transformed the treatment of HCV. However, it has become clear that, without an effective HCV vaccine, it will not be possible to meet the World Health Organization targets of HCV viral elimination. Promising new vaccine technologies that generate high magnitude antiviral T and B cell immune responses and significant new funding have recently become available, stimulating the HCV vaccine pipeline. In the absence of an immune competent animal model for HCV, the major block in evaluating new HCV vaccine candidates will be the assessment of vaccine efficacy in humans. The development of a controlled human infection model (CHIM) for HCV could overcome this block, enabling the head-to-head assessment of vaccine candidates. The availability of highly effective DAA means that a CHIM for HCV is possible for the first time. In this review, we highlight the challenges and issues with currently available strategies to assess HCV vaccine efficacy including HCV "at-risk" cohorts and animal models. We describe the development of CHIM in other infections that are increasingly utilized by trialists and explore the ethical and safety concerns specific for an HCV CHIM. Finally, we propose an HCV CHIM study design including the selection of volunteers, the development of an infectious inoculum, the evaluation of host immune and viral parameters, and the definition of study end points for use in an HCV CHIM. Importantly, the study design (including number of volunteers required, cost, duration of study, and risk to volunteers) varies significantly depending on the proposed mechanism of action (sterilizing/rapid viral clearance vs. delayed viral clearance) of the vaccine under evaluation. We conclude that an HCV CHIM is now realistic, that safety and ethical concerns can be addressed with the right study design, and that, without an HCV CHIM, it is difficult to envisage how the development of an HCV vaccine will be possible.

Probe Capture Enrichment Methods for HIV and HCV Genome Sequencing and Drug Resistance Genotyping

Human immunodeficiency virus (HIV) infections remain a significant public health concern worldwide. Over the years, sophisticated sequencing technologies such as next-generation sequencing (NGS) have emerged and been utilized to monitor the spread of HIV drug resistance (HIVDR), identify HIV drug resistance mutations, and characterize transmission dynamics. Similar applications also apply to the Hepatitis C virus (HCV), another bloodborne viral pathogen with significant intra-host genetic diversity. Several advantages to using NGS over conventional Sanger sequencing include increased data throughput, scalability, cost-effectiveness when batched sample testing is performed, and sensitivity for quantitative detection of minority resistant variants. However, NGS alone may fail to detect genomes from pathogens present in low copy numbers. As with all sequencing platforms, the primary determinant in achieving quality sequencing data is the quality and quantity of the initial template input. Samples containing degraded RNA/DNA and/or low copy number have been a consistent sequencing challenge. To overcome this limitation probe capture enrichment is a method that has recently been employed to target, enrich, and sequence the genome of a pathogen present in low copies, and for compromised specimens that contain poor quality nucleic acids. It involves the hybridization of sequence-specific DNA or RNA probes to a target sequence, which is followed by an enrichment step via PCR to increase the number of copies of the targeted sequences after which the samples are subjected to NGS procedures. This method has been performed on pathogens such as bacteria, fungus, and viruses and allows for the sequencing of complete genomes, with high coverage. Post NGS, data analysis can be performed through various bioinformatics pipelines which can provide information on genetic diversity, genotype, virulence, and drug resistance. This article reviews how probe capture enrichment helps to increase the likelihood of sequencing HIV and HCV samples that contain low viral loads and/or are compromised.

SNP and Phylogenetic Characterization of Low Viral Load SARS-CoV-2 Specimens by Target Enrichment

Background: Surveillance of SARS-CoV-2 across the globe has enabled detection of new variants and informed the public health response. With highly sensitive methods like qPCR widely adopted for diagnosis, the ability to sequence and characterize specimens with low titers needs to keep pace.

Methods: Nucleic acids extracted from nasopharyngeal swabs collected from four sites in the United States in early 2020 were converted to NGS libraries to sequence SARS-CoV-2 genomes using metagenomic and xGen target enrichment approaches. Single nucleotide polymorphism (SNP) analysis and phylogeny were used to determine clade assignments and geographic origins of strains.

Results: SARS-CoV-2-specific xGen enrichment enabled full genome coverage for 87 specimens with Ct values &lt;29, corresponding to viral loads of &gt;10,000 cp/ml. For samples with viral loads between 103 and 106 cp/ml, the median genome coverage for xGen was 99.1%, sequence depth was 605X, and the “on-target” rate was 57 ± 21%, compared to 13%, 2X and 0.001 ± 0.016%, respectively, for metagenomic sequencing alone. Phylogenetic analysis revealed the presence of most clades that existed at the time of the study, though clade GH dominated in the Midwest.

Conclusions: Even as vaccines are being widely distributed, a high case load of SARS-CoV-2 infection persists around the world. Viral genetic surveillance has succeeded in warning the public of new variants in circulation and ensured that diagnostic tools remain resilient to a steadily increasing number of mutations. Target capture offers a means of characterizing low viral load samples which would normally pose a challenge for metagenomic sequencing.

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.

High Cure Rates for Hepatitis C Virus Genotype 6 in Advanced Liver Fibrosis With 12 Weeks Sofosbuvir and Daclatasvir: The Vietnam SEARCH Study

BACKGROUND

Genotype 6 is the most genetically diverse lineage of hepatitis C virus, and it predominates in Vietnam. It can be treated with sofosbuvir with daclatasvir (SOF/DCV), the least expensive treatment combination globally. In regional guidelines, longer treatment durations of SOF/DCV (24 weeks) are recommended for cirrhotic individuals, compared with other pangenotypic regimens (12 weeks), based on sparse data. Early on-treatment virological response may offer means of reducing length and cost of therapy in patients with liver fibrosis.

METHODS

In this prospective trial in Vietnam, genotype 6-infected adults with advanced liver fibrosis or compensated cirrhosis were treated with SOF/DCV. Day 14 viral load was used to guide duration of therapy: participants with viral load <500 IU/mL at day 14 were treated with 12 weeks of SOF/DCV and those ≥500 IU/mL received 24 weeks. Primary endpoint was sustained virological response (SVR).

RESULTS

Of 41 individuals with advanced fibrosis or compensated cirrhosis who commenced treatment, 51% had genotype 6a and 34% had 6e. The remainder had 6h, 6k, 6l, or 6o. One hundred percent had viral load <500 IU/mL by day 14, meaning that all received 12 weeks of SOF/DCV. One hundred percent achieved SVR12 despite a high frequency of putative NS5A inhibitor resistance-associated substitutions at baseline.

CONCLUSIONS

Prescribing 12 weeks of SOF/DCV results in excellent cure rates in this population. These data support the removal of costly genotyping in countries where genotype 3 prevalence is <5%, in keeping with World Health Organization guidelines. NS5A resistance-associated mutations in isolation do not affect efficacy of SOF/DCV therapy. Wider evaluation of response-guided therapy is warranted.

Enrichment of low abundance DNA/RNA by oligonucleotide-clicked iron oxide nanoparticles

Detection of low abundance target DNA/RNA for clinical or research purposes is challenging because the target sequences can be hidden under a large background of human genomic or non-human metagenomic sequences. We describe a probe-based capture method to enrich for target sequences with DNA-clicked iron oxide nanoparticles. Our method was tested against commercial capture assays using streptavidin beads, on a set of probes derived from a common genotype of the hepatitis C virus. We showed that our method is more specific and sensitive, most likely due to the combination of an inert silica coating and a high density of DNA probes clicked to the nanoparticles. This facilitates target capture below the limits of detection for TaqMan qPCR, and we believe that this method has the potential to transform management of infectious diseases.

An enrichment method for capturing mumps virus whole genome sequences directly from clinical specimens

Recently there has been a significant increase in the number of mumps outbreaks occurring in highly vaccinated populations. These outbreaks are often due to a mumps genotype G virus, where sequencing of the SH gene does not reveal enough genetic diversity to sufficient to resolve outbreaks. This has elevated the need to be able to sequence complete mumps viruses from clinical samples without laborious methods. Here we describe a probe enrichment method that allows for whole genome sequencing of the mumps virus directly from clinical specimens. Using 136 clinical samples, we show this method allows for a significant increase in the percentage of viral sequencing reads, resulting in the capture of mumps genomes. This method will be an asset in investigating future mumps outbreaks.

SARS-CoV-2 within-host diversity and transmission

Extensive global sampling and sequencing of the pandemic virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have enabled researchers to monitor its spread and to identify concerning new variants. Two important determinants of variant spread are how frequently they arise within individuals and how likely they are to be transmitted. To characterize within-host diversity and transmission, we deep-sequenced 1313 clinical samples from the United Kingdom. SARS-CoV-2 infections are characterized by low levels of within-host diversity when viral loads are high and by a narrow bottleneck at transmission. Most variants are either lost or occasionally fixed at the point of transmission, with minimal persistence of shared diversity, patterns that are readily observable on the phylogenetic tree. Our results suggest that transmission-enhancing and/or immune-escape SARS-CoV-2 variants are likely to arise infrequently but could spread rapidly if successfully transmitted.

Targeted capture and sequencing of Orientia tsutsugamushi genomes from chiggers and humans

Scrub typhus is a febrile disease caused by Orientia tsutsugamushi, transmitted by larval stage Trombiculid mites (chiggers), whose primary hosts are small mammals. The phylogenomics of O. tsutsugamushi in chiggers, small mammals and humans remains poorly understood. To combat the limitations imposed by the low relative quantities of pathogen DNA in typical O. tsutsugamushi clinical and ecological samples, along with the technical, safety and cost limitations of cell culture, a novel probe-based target enrichment sequencing protocol was developed. The method was designed to capture variation among conserved genes and facilitate phylogenomic analysis at the scale of population samples. A whole-genome amplification step was incorporated to enhance the efficiency of sequencing by reducing duplication rates. This resulted in on-target capture rates of up to 93% for a diverse set of human, chigger, and rodent samples, with the greatest success rate in samples with real-time PCR C_t values below 35. Analysis of the best-performing samples revealed phylogeographic clustering at local, provincial and international scales. Applying the methodology to a comprehensive set of samples could yield a more complete understanding of the ecology, genomic evolution and population structure of O. tsutsugamushi and other similarly challenging organisms, with potential benefits in the development of diagnostic tests and vaccines.

Simultaneous Viral Whole-Genome Sequencing and Differential Expression Profiling in Respiratory Syncytial Virus Infection of Infants

Targeted metagenomics using strand-specific libraries with target enrichment is a sensitive, generalized approach to pathogen sequencing and transcriptome profiling. Using this method, we recovered 13 (76%) complete human respiratory syncytial virus (RSV) genomes from 17 clinical respiratory samples, reconstructed the phylogeny of the infecting viruses, and detected differential gene expression between 2 RSV subgroups, specifically, a lower expression of the P gene and a higher expression of the M2 gene in RSV-A than in RSV-B. This methodology can help to relate viral genetics to clinical phenotype and facilitate ongoing population-level RSV surveillance and vaccine development. Clinical Trials Registration. NCT03627572 and NCT03756766.

Performance of a high-throughput next-generation sequencing method for analysis of HIV drug resistance and viral load

OBJECTIVES

To evaluate the performance of a high-throughput research assay for HIV drug resistance testing based on whole genome next-generation sequencing (NGS) that also quantifies HIV viral load.

METHODS

Plasma samples (n = 145) were obtained from HIV-positive MSM (HPTN 078). Samples were analysed using clinical assays (the ViroSeq HIV-1 Genotyping System and the Abbott RealTime HIV-1 Viral Load assay) and a research assay based on whole-genome NGS (veSEQ-HIV).

RESULTS

HIV protease and reverse transcriptase sequences (n = 142) and integrase sequences (n = 138) were obtained using ViroSeq. Sequences from all three regions were obtained for 100 (70.4%) of the 142 samples using veSEQ-HIV; results were obtained more frequently for samples with higher viral loads (93.5% for 93 samples with >5000 copies/mL; 50.0% for 26 samples with 1000-5000 copies/mL; 0% for 23 samples with <1000 copies/mL). For samples with results from both methods, drug resistance mutations (DRMs) were detected in 33 samples using ViroSeq and 42 samples using veSEQ-HIV (detection threshold: 5.0%). Overall, 146 major DRMs were detected; 107 were detected by both methods, 37 were detected by veSEQ-HIV only (frequency range: 5.0%-30.6%) and two were detected by ViroSeq only. HIV viral loads estimated by veSEQ-HIV strongly correlated with results from the Abbott RealTime Viral Load assay (R2 = 0.85; n = 142).

CONCLUSIONS

The NGS-based veSEQ-HIV method provided results for most samples with higher viral loads, was accurate for detecting major DRMs, and detected mutations at lower levels compared with a method based on population sequencing. The veSEQ-HIV method also provided HIV viral load data.

Technical Validation of a Hepatitis C Virus Whole Genome Sequencing Assay for Detection of Genotype and Antiviral Resistance in the Clinical Pathway

Choice of direct acting antiviral (DAA) therapy for Hepatitis C Virus (HCV) in the United Kingdom and similar settings usually requires knowledge of the genotype and, in some cases, antiviral resistance (AVR) profile of the infecting virus. To determine these, most laboratories currently use Sanger technology, but next-generation sequencing (NGS) offers potential advantages in throughput and accuracy. However, NGS poses unique technical challenges, which require idiosyncratic development and technical validation approaches. This applies particularly to virology, where sequence diversity is high and the amount of starting genetic material is low, making it difficult to distinguish real data from artifacts. We describe the development and technical validation of a sequence capture-based HCV whole genome sequencing (WGS) assay to determine viral genotype and AVR profile. We use clinical samples of known subtypes and viral loads, and simulated FASTQ datasets to validate the analytical performances of both the wet laboratory and bioinformatic pipeline procedures. We show high concordance of the WGS assay compared to current "gold standard" Sanger assays. Specificity was 92.3 and 96.1% for AVR and genotyping, respectively. Discordances were due to the inability of Sanger assays to assign the correct subtype or accurately call mixed drug-resistant variants. We show high repeatability and reproducibility with >99.8% sequence similarity between sequence runs as well as high precision for variant frequency detection at >98.8% in the 95th percentile. Post-sequencing bioinformatics quality control workflows allow the accurate distinction between mixed infections, cross-contaminants and recombinant viruses at a threshold of >5% for the minority population. The sequence capture-based HCV WGS assay is more accurate than legacy AVR and genotyping assays. The assay has now been implemented in the clinical pathway of England's National Health Service HCV treatment programs, representing the first validated HCV WGS pipeline in clinical service. The data generated will additionally provide granular national-level genomic information for public health policy making and support the WHO HCV elimination strategy.

A Comprehensive Genomics Solution for HIV Surveillance and Clinical Monitoring in Low-Income Settings

Viral genetic sequencing can be used to monitor the spread of HIV drug resistance, identify appropriate antiretroviral regimes, and characterize transmission dynamics. Despite decreasing costs, next-generation sequencing (NGS) is still prohibitively costly for routine use in generalized HIV epidemics in low- and middle-income countries. Here, we present veSEQ-HIV, a high-throughput, cost-effective NGS sequencing method and computational pipeline tailored specifically to HIV, which can be performed using leftover blood drawn for routine CD4 cell count testing. This method overcomes several major technical challenges that have prevented HIV sequencing from being used routinely in public health efforts; it is fast, robust, and cost-efficient, and generates full genomic sequences of diverse strains of HIV without bias. The complete veSEQ-HIV pipeline provides viral load estimates and quantitative summaries of drug resistance mutations; it also exploits information on within-host viral diversity to construct directed transmission networks. We evaluated the method's performance using 1,620 plasma samples collected from individuals attending 10 large urban clinics in Zambia as part of the HPTN 071-2 study (PopART Phylogenetics). Whole HIV genomes were recovered from 91% of samples with a viral load of >1,000 copies/ml. The cost of the assay (30 GBP per sample) compares favorably with existing VL and HIV genotyping tests, proving an affordable option for combining HIV clinical monitoring with molecular epidemiology and drug resistance surveillance in low-income settings.

Genome-wide capture sequencing to detect hepatitis C virus at the end of antiviral therapy

Background

Viral relapse is a major concern in hepatitis C virus (HCV) antiviral therapy. Currently, there are no satisfactory methods to predict viral relapse, especially in the era of direct acting antivirals in which the virus often quickly becomes undetectable using PCR-based approaches that focus on a small viral region. Next-generation sequencing (NGS) provides an alternative option for viral detection in a genome-wide manner. However, owing to the overwhelming dominance of human genetic content in clinical specimens, direct detection of HCV by NGS has a low sensitivity and hence viral enrichment is required.

Methods

Based on template-dependent multiple displacement amplification (tdMDA), an improved method for whole genome amplification (Wang et al., 2017. Biotechniques 63, 21–27), we evaluated two strategies to enhance the sensitivity of NGS-based HCV detection: duplex-specific nuclease (DSN)-mediated depletion of human sequences and HCV probe-based capture sequencing.

Results

In DSN-mediated depletion, human sequences were significantly reduced in the two HCV serum samples tested, 65.3% → 55.6% → 33.7% (#4727) and 68.6% → 56% → 21% (#4970), respectively for no normalization, self- and driver-applied normalization. However, this approach was associated with a loss of HCV sequences perhaps due to its micro-homology with the human genome. In capture sequencing, HCV-mapped sequencing reads occupied 96.8% (#4727) and 22.14% (#4970) in NGS data, equivalent to 1936x and 7380x enrichment, respectively. Capture sequencing was then applied to ten serum samples collected at the end of HCV antiviral therapy. Interestingly, the number of HCV-mapped reads was significantly higher in relapsed patients (n = 5) than those from patients with sustained virological response (SVR) (n = 5), 102.4 ± 72.3 vs. 2.6 ± 0.55, p = 0.014.

Conclusions

Our data provides concept evidence for a highly sensitive HCV detection by capture sequencing. The abundance difference of HCV sequencing reads at the end of HCV antiviral therapy could be applied to predict treatment outcomes.

Advanced molecular surveillance approaches for characterization of blood borne hepatitis viruses

Defining genetic diversity of viral infections directly from patient specimens is the ultimate goal of surveillance. Simple tools that can provide full-length sequence information on blood borne viral hepatitis viruses: hepatitis C, hepatitis B and hepatitis D viruses (HCV, HBV and HDV) remain elusive. Here, an unbiased metagenomic next generation sequencing approach (mNGS) was used for molecular characterization of HCV infections (n = 99) from Israel which yielded full-length HCV sequences in 89% of samples, with 7 partial sequences sufficient for classification. HCV genotypes were primarily 1b (68%) and 1a (19%), with minor representation of genotypes 2c (1%) and 3a (8%). HBV/HDV coinfections were characterized by suppressed HBV viral loads, resulting in sparse mNGS coverage. A probe-based enrichment approach (xGen) aiming to increase HBV and HDV coverage was validated on a panel of diverse genotypes, geography and titers. The method extended HBV genome coverage a median 61% (range 8–84%) and provided orders of magnitude boosts in reads and sequence depth for both viruses. When HBV-xGen was applied to Israeli samples, coverage was improved by 28–73% in 4 samples and identified HBV genotype A1, A2, D1 specimens and a dual B/D infection. Abundant HDV reads in mNGS libraries yielded 18/26 (69%) full genomes and 8 partial sequences, with HDV-xGen only providing minimal extension (3–11%) of what were all genotype 1 genomes. Advanced molecular approaches coupled to virus-specific capture probes promise to enhance surveillance of viral infections and aid in monitoring the spread of local subtypes.

Transfusion-transmitted hepatitis C: A cluster of cases in transfusion-dependent thalassaemia patients in Sri Lanka

OBJECTIVES

To report the clinical and virologic epidemiology of a recent epidemic of hepatitis C in thalassaemia patients in Sri Lanka.

BACKGROUND

Transfusion-dependent thalassaemia patients remain at risk for hepatitis C virus (HCV). Here, we report a cluster of recent HCV infections in Sri Lankan thalassaemia patients and examine the phylogenetic relationship of viral sequences.

METHODS

We conducted two prospective cross-sectional surveys of 513 patients in four Sri Lankan thalassaemia centres in 2014/2015 and re-surveyed one centre in 2016. We screened for anti-HCV antibodies using the CTK Biotech enzyme-linked immunosorbent assay (ELISA) kits and confirmed active infection by reverse transcription-polymerase chain reaction (RT-PCR) for HCV-RNA. HCV genomes were sequenced by unbiased target enrichment.

RESULTS

Anti-HCV antibodies were found in 116/513 (22.6%) of patients initially tested. Active hepatitis C infection was found in 26 patients with no cases of active hepatitis B infection. Of 26 patients with HCV, two were infected with genotype 1(a), and the rest had 3(a). In a single centre (Ragama), 122 patients (120 new cases and two previously tested, but negative) were retested for anti-HCV antibodies. 32/122 (26.2%) patients were seropositive. Twenty-three (23/122; 18.8%) of these new cases were confirmed by HCV PCR (all genotype 3[a]).

CONCLUSION

There is a significant cluster of recent HCV cases in multiply transfused thalassaemia patients in several centres in Sri Lanka. Most of the viruses shared a close phylogenetic relationship. The results are consistent with recent continuing transfusion-transmitted HCV infection. Routine surveillance for HCV of chronically transfused patients is required irrespective of screening of blood products.

Human pegivirus 2 exhibits minimal geographic and temporal genetic diversity

We applied an NGS based target capture approach to amplify HPgV-2 sequences from metagenomic libraries and enable full genome characterization. Despite expanded geographical sampling, sequence variability remains low, with diversity concentrated in approximately 3.3% of all amino acids. Serial samples from one HPgV-2 positive individual co-infected with comparable titers of HIV, HCV, and GBV-C showed that HPgV-2 remains highly stable over several weeks compared to other RNA viruses, despite a similarly error-prone polymerase. The consistent epidemiological association with and structural similarities to HCV, and the weak positive correlation of HCV and HPgV-2 titers shown here, suggests it may benefit from co-infection. While minimal selective pressure on HPgV-2 to evolve could suggest fitness, the rarity of HPgV-2 and the tight phylogenetic clustering of global strains likely indicates origination from a common source and a virus that is ill-suited to its host. Sporadic infections may explain the limited genetic diversity observed worldwide.

Probe capture enrichment next-generation sequencing of complete foot-and-mouth disease virus genomes in clinical samples

Next-generation sequencing (NGS) techniques offer an unprecedented "step-change" increase in the quantity and quality of sequence data rapidly generated from a sample and can be applied to obtain ultra-deep coverage of viral genomes. This is not possible with the routinely used Sanger sequencing method that gives the consensus reads, or by cloning approaches. In this study, a targeted-enrichment methodology for the simultaneous acquisition of complete foot-and-mouth disease virus (FMDV) genomes directly from clinical samples is presented. Biotinylated oligonucleotide probes (120 nt) were used to capture and enrich viral RNA following library preparation. To create a virus capture panel targeting serotype O and A simultaneously, 18 baits targeting the highly conserved regions of the 8.3 kb FMDV genome were synthesised, with 14 common to both serotypes, 2 specific to serotype O and 2 specific to serotype A. These baits were used to capture and enrich FMDV RNA (as cDNA) from samples collected during one pathogenesis and two vaccine efficacy trials, where pigs were infected with serotype O or A viruses. After enrichment, FMDV-specific sequencing reads increased by almost 3000-fold. The sequence data were used in variant call analysis to identify single nucleotide polymorphisms (SNPs). This methodology was robust in its ability to capture diverse sequences, was shown to be highly sensitive, and can be easily scaled for large-scale epidemiological studies.

Next-generation sequencing for the clinical management of hepatitis C virus infections: does one test fits all purposes?

While the prospect of viral cure is higher than ever for individuals infected with the hepatitis C virus (HCV) due to ground-breaking progress in antiviral treatment, success rates are still negatively influenced by HCV's high genetic variability. This genetic diversity is represented in the circulation of various genotypes and subtypes, mixed infections, recombinant forms and the presence of numerous drug resistant variants among infected individuals. Common misclassifications by commercial genotyping assays in combination with the limitations of currently used targeted population sequencing approaches have encouraged researchers to exploit alternative methods for the clinical management of HCV infections. Next-generation sequencing (NGS), a revolutionary and powerful tool with a variety of applications in clinical virology, can characterize viral diversity and depict viral dynamics in an ultra-wide and ultra-deep manner. The level of detail it provides makes it the method of choice for the diagnosis and clinical assessment of HCV infections. The sequence library provided by NGS is of a higher magnitude and sensitivity than data generated by conventional methods. Therefore, these technologies are helpful to guide clinical practice and at the same time highly valuable for epidemiological studies. The decreasing costs of NGS to determine genotypes, mixed infections, recombinant strains and drug resistant variants will soon make it feasible to employ NGS in clinical laboratories, to assist in the daily care of patients with HCV.

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).

Application of a targeted-enrichment methodology for full-genome sequencing of Dengue 1-4, Chikungunya and Zika viruses directly from patient samples

The frequency of epidemics caused by Dengue viruses 1-4, Zika virus and Chikungunya viruses have been on an upward trend in recent years driven primarily by uncontrolled urbanization, mobility of human populations and geographical spread of their shared vectors, Aedes aegypti and Aedes albopictus. Infections by these viruses present with similar clinical manifestations making them challenging to diagnose; this is especially difficult in regions of the world hyperendemic for these viruses. In this study, we present a targeted-enrichment methodology to simultaneously sequence the complete viral genomes for each of these viruses directly from clinical samples. Additionally, we have also developed a customized computational tool (BaitMaker) to design these enrichment baits. This methodology is robust in its ability to capture diverse sequences and is amenable to large-scale epidemiological studies. We have applied this methodology to two large cohorts: a febrile study based in Colombo, Sri Lanka taken during the 2009-2015 dengue epidemic (n = 170) and another taken during the 2016 outbreak of Zika virus in Singapore (n = 162). Results from these studies indicate that we were able to cover an average of 97.04% ± 0.67% of the full viral genome from samples in these cohorts. We also show detection of one DENV3/ZIKV co-infected patient where we recovered full genomes for both viruses.

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.

Capturing sequence diversity in metagenomes with comprehensive and scalable probe design

Metagenomic sequencing has the potential to transform microbial detection and characterization, but new tools are needed to improve its sensitivity. Here we present CATCH, a computational method to enhance nucleic acid capture for enrichment of diverse microbial taxa. CATCH designs optimal probe sets, with a specified number of oligonucleotides, that achieve full coverage of, and scale well with, known sequence diversity. We focus on applying CATCH to capture viral genomes in complex metagenomic samples. We design, synthesize, and validate multiple probe sets, including one that targets the whole genomes of the 356 viral species known to infect humans. Capture with these probe sets enriches unique viral content on average 18-fold, allowing us to assemble genomes that could not be recovered without enrichment, and accurately preserves within-sample diversity. We also use these probe sets to recover genomes from the 2018 Lassa fever outbreak in Nigeria and to improve detection of uncharacterized viral infections in human and mosquito samples. The results demonstrate that CATCH enables more sensitive and cost-effective metagenomic sequencing.

Universal Target Capture of HIV Sequences From NGS Libraries

Background: Global surveillance of viral sequence diversity is needed to keep pace with the constant evolution of HIV. Recent next generation sequencing (NGS) methods have realized the goal of sequencing circulating virus directly from patient specimens. Yet, a simple, universal approach that maximizes sensitivity and sequencing capacity remains elusive. Here we present a novel HIV enrichment strategy to yield near complete genomes from low viral load specimens. Methodology: A non-redundant biotin-labeled probe set (HIV-xGen; n = 652) was synthesized to tile all HIV-1 (groups M, N, O, and P) and HIV-2 (A and B) strains. Illumina Nextera barcoded libraries of either gene-specific or randomly primed cDNA derived from infected plasma were hybridized to probes in a single pool and unbound sequences were washed away. Captured viral cDNA was amplified by Illumina adaptor primers, sequenced on a MiSeq, and NGS reads were demultiplexed for alignment with CLC Bio software. Results: HIV-xGen probes selectively captured and amplified reads spanning the entirety of the HIV phylogenetic tree. HIV sequences clearly present in unenriched libraries of specimens but previously not observed due to high host background levels, insufficient sequencing depth or the extent of multiplexing, were now enriched by >1,000-fold. Thus, xGen selection not only substantially increased the depth of existing sequence, but also extended overall genome coverage by an average of 40%. We characterized 50 new, diverse HIV strains from clinical specimens and demonstrated a viral load cutoff of approximately log 3.5 copies/ml for full length coverage. Genome coverage was <20% for 5/10 samples with viral loads <log 3.5 copies/ml and >90% for 35/40 samples with higher viral loads. Conclusions: Characterization of >20 complete genomes at a time is now possible from a single probe hybridization and MiSeq run. With the versatility to capture all HIV strains and the sensitivity to detect low titer specimens, HIV-xGen will serve as an important tool for monitoring HIV sequence diversity.

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.

Genomic Epidemiology Reconstructs the Introduction and Spread of Zika Virus in Central America and Mexico

The Zika virus (ZIKV) epidemic in the Americas established ZIKV as a major public health threat and uncovered its association with severe diseases, including microcephaly. However, genetic epidemiology in some at-risk regions, particularly Central America and Mexico, remains limited. We report 61 ZIKV genomes from this region, generated using metagenomic sequencing with ZIKV-specific enrichment, and combine phylogenetic, epidemiological, and environmental data to reconstruct ZIKV transmission. These analyses revealed multiple independent ZIKV introductions to Central America and Mexico. One introduction, likely from Brazil via Honduras, led to most infections and the undetected spread of ZIKV through the region from late 2014. Multiple lines of evidence indicate biannual peaks of ZIKV transmission in the region, likely driven by varying local environmental conditions for mosquito vectors and herd immunity. The spatial and temporal heterogeneity of ZIKV transmission in Central America and Mexico challenges arbovirus surveillance and disease control measures.

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.

PHYLOSCANNER: Inferring Transmission from Within- and Between-Host Pathogen Genetic Diversity

A central feature of pathogen genomics is that different infectious particles (virions and bacterial cells) within an infected individual may be genetically distinct, with patterns of relatedness among infectious particles being the result of both within-host evolution and transmission from one host to the next. Here, we present a new software tool, phyloscanner, which analyses pathogen diversity from multiple infected hosts. phyloscanner provides unprecedented resolution into the transmission process, allowing inference of the direction of transmission from sequence data alone. Multiply infected individuals are also identified, as they harbor subpopulations of infectious particles that are not connected by within-host evolution, except where recombinant types emerge. Low-level contamination is flagged and removed. We illustrate phyloscanner on both viral and bacterial pathogens, namely HIV-1 sequenced on Illumina and Roche 454 platforms, HCV sequenced with the Oxford Nanopore MinION platform, and Streptococcus pneumoniae with sequences from multiple colonies per individual. phyloscanner is available from https://github.com/BDI-pathogens/phyloscanner.

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.

Sequencing of hepatitis C virus for detection of resistance to direct-acting antiviral therapy: A systematic review

The significance of the clinical impact of direct‐acting antiviral (DAA) resistance‐associated substitutions (RASs) in hepatitis C virus (HCV) on treatment failure is unclear. No standardized methods or guidelines for detection of DAA RASs in HCV exist. To facilitate further evaluations of the impact of DAA RASs in HCV, we conducted a systematic review of RAS sequencing protocols, compiled a comprehensive public library of sequencing primers, and provided expert guidance on the most appropriate methods to screen and identify RASs. The development of standardized RAS sequencing protocols is complicated due to a high genetic variability and the need for genotype‐ and subtype‐specific protocols for multiple regions. We have identified several limitations of the available methods and have highlighted areas requiring further research and development. The development, validation, and sharing of standardized methods for all genotypes and subtypes should be a priority. (Hepatology Communications 2017;1:379–390)

Genome-to-genome analysis highlights the effect of the human innate and adaptive immune systems on the hepatitis C virus

Outcomes of hepatitis C virus (HCV) infection and treatment depend on viral and host genetic factors. We use human genome-wide genotyping arrays and new whole-genome HCV viral sequencing technologies to perform a systematic genome-to-genome study of 542 individuals chronically infected with HCV, predominately genotype 3. We show that both HLA alleles and interferon lambda innate immune system genes drive viral genome polymorphism, and that IFNL4 genotypes determine HCV viral load through a mechanism that is dependent on a specific polymorphism in the HCV polyprotein. We highlight the interplay between innate immune responses and the viral genome in HCV control.

Clinical and biological insights from viral genome sequencing

Whole-genome sequencing (WGS) of pathogens is becoming increasingly important not only for basic research but also for clinical science and practice. In virology, WGS is important for the development of novel treatments and vaccines, and for increasing the power of molecular epidemiology and evolutionary genomics. In this Opinion article, we suggest that WGS of viruses in a clinical setting will become increasingly important for patient care. We give an overview of different WGS methods that are used in virology and summarize their advantages and disadvantages. Although there are only partially addressed technical, financial and ethical issues in regard to the clinical application of viral WGS, this technique provides important insights into virus transmission, evolution and pathogenesis.

The Promise of Whole Genome Pathogen Sequencing for the Molecular Epidemiology of Emerging Aquaculture Pathogens

Aquaculture is the fastest growing food-producing sector, and the sustainability of this industry is critical both for global food security and economic welfare. The management of infectious disease represents a key challenge. Here, we discuss the opportunities afforded by whole genome sequencing of bacterial and viral pathogens of aquaculture to mitigate disease emergence and spread. We outline, by way of comparison, how sequencing technology is transforming the molecular epidemiology of pathogens of public health importance, emphasizing the importance of community-oriented databases and analysis tools.

Characterization of the Specificity, Functionality, and Durability of Host T-Cell Responses Against the Full-Length Hepatitis E Virus

The interplay between host antiviral immunity and immunopathology during hepatitis E virus (HEV) infection determines important clinical outcomes. We characterized the specificity, functionality, and durability of host T-cell responses against the full-length HEV virus and assessed a novel "Quantiferon" assay for the rapid diagnosis of HEV infection. Eighty-nine volunteers were recruited from Oxford, Truro (UK), and Toulouse (France), including 44 immune-competent patients with acute HEV infection, 18 HEV-exposed immunosuppressed organ-transplant recipients (8 with chronic HEV), and 27 healthy volunteers. A genotype 3a peptide library (616 overlapping peptides spanning open reading frames [ORFs] 1-3) was used in interferon-gamma (IFN-γ) T-cell ELISpot assays. CD4+ /CD8+ T-cell subsets and polyfunctionality were defined using ICCS and SPICE analysis. Quantification of IFN-γ used whole-blood stimulation with recombinant HEV-capsid protein in the QuantiFERON kit. HEV-specific T-cell responses were detected in 41/44 immune-competent HEV exposed volunteers (median magnitude: 397 spot-forming units/106 peripheral blood mononuclear cells), most frequently targeting ORF2. High-magnitude, polyfunctional CD4 and CD8+ T cells were detected during acute disease and maintained to 12 years, but these declined over time, with CD8+ responses becoming more monofunctional. Low-level responses were detectable in immunosuppressed patients. Twenty-three novel HEV CD4+ and CD8+ T-cell targets were mapped predominantly to conserved genomic regions. QuantiFERON testing demonstrated an inverse correlation between IFN-γ production and the time from clinical presentation, providing 100% specificity, and 71% sensitivity (area under the receiver operator characteristic curve of 0.86) for HEV exposure at 0.3 IU/mL.

CONCLUSION

Robust HEV-specific T-cell responses generated during acute disease predominantly target ORF2, but decline in magnitude and polyfunctionality over time. Defining HEV T-cell targets will be important for the investigation of HEV-associated autoimmune disease. (Hepatology 2016;64:1934-1950).