A method for near full-length amplification and sequencing for six hepatitis C virus genotypes

The role of the hepatitis B virus genome and its integration in the hepatocellular carcinoma

The integration of Hepatitis B Virus (HBV) is now known to be closely associated with the occurrence of liver cancer and can impact the functionality of liver cells through multiple dimensions. However, despite the detailed understanding of the characteristics of HBV integration and the mechanisms involved, the subsequent effects on cellular function are still poorly understood in current research. This study first systematically discusses the relationship between HBV integration and the occurrence of liver cancer, and then analyzes the status of the viral genome produced by HBV replication, highlighting the close relationship and structure between double-stranded linear (DSL)-HBV DNA and the occurrence of viral integration. The integration of DSL-HBV DNA leads to a certain preference for HBV integration itself. Additionally, exploration of HBV integration hotspots reveals obvious hotspot areas of HBV integration on the human genome. Virus integration in these hotspot areas is often associated with the occurrence and development of liver cancer, and it has been determined that HBV integration can promote the occurrence of cancer by inducing genome instability and other aspects. Furthermore, a comprehensive study of viral integration explored the mechanisms of viral integration and the internal integration mode, discovering that HBV integration may form extrachromosomal DNA (ecDNA), which exists outside the chromosome and can integrate into the chromosome under certain conditions. The prospect of HBV integration as a biomarker was also probed, with the expectation that combining HBV integration research with CRISPR technology will vigorously promote the progress of HBV integration research in the future. In summary, exploring the characteristics and mechanisms in HBV integration holds significant importance for an in-depth comprehension of viral integration.

Single-molecule sequencing of the whole HCV genome revealed envelope deletions in decompensated cirrhosis associated with NS2 and NS5A mutations

BACKGROUND

Defective hepatitis C virus (HCV) genomes with deletion of the envelope region have been occasionally reported by short-read sequencing analyses. However, the clinical and virological details of such deletion HCV have not been fully elucidated.

METHODS

We developed a highly accurate single-molecule sequencing system for full-length HCV genes by combining the third-generation nanopore sequencing with rolling circle amplification (RCA) and investigated the characteristics of deletion HCV through the analysis of 21 patients chronically infected with genotype-1b HCV.

RESULT

In 5 of the 21 patients, a defective HCV genome with approximately 2000 bp deletion from the E1 to NS2 region was detected, with the read frequencies of 34-77%, suggesting the trans-complementation of the co-infecting complete HCV. Deletion HCV was found exclusively in decompensated cirrhosis (5/12 patients), and no deletion HCV was observed in nine compensated patients. Comparing the amino acid substitutions between the deletion and complete HCV (DAS, deletion-associated substitutions), the deletion HCV showed higher amino acid mutations in the ISDR (interferon sensitivity-determining region) in NS5A, and also in the TMS (transmembrane segment) 3 to H (helix) 2 region of NS2.

CONCLUSIONS

Defective HCV genome with deletion of envelope genes is associated with decompensated cirrhosis. The deletion HCV seems susceptible to innate immunity, such as endogenous interferon with NS5A mutations, escaping from acquired immunity with deletion of envelope proteins with potential modulation of replication capabilities with NS2 mutations. The relationship between these mutations and liver damage caused by HCV deletion is worth investigating.

Determination of hepatitis C virus subtype prevalent in Sindh, Pakistan: a phylogenetic analysis

Hepatitis is a major public health issue, affecting 10-17 million people worldwide, with its prevalence continuously increasing. The Hepatitis C virus (HCV) is responsible for liver related diseases, which include liver cirrhosis, hepatocellular carcinoma, and chronic hepatitis. Pakistan is experiencing a serious rise in HCV cases. This study aimed to assess the prevalence and distribution of HCV genotypes in Sindh, Pakistan. Serum samples from HCV-positive patients were collected from various local hospitals in Sindh. These samples were first screened for HCV antibodies using ELISA. Samples that tested positive for HCV RNA underwent further genotyping through sequencing using the standard Sanger method. The genotypes were identified by comparing the sequences with those available in the National Center for Biotechnology Information (NCBI) database, and a phylogenetic tree was constructed. The phylogenetic analysis showed that all isolates in this study were clustered with genotypes 3a and 3b, except for one sequence that was clustered with genotype 1a. No isolates were found to be clustered with reference genomes of genotypes 2, 4, 5, 6, and 7 suggesting that genotype 3a is endemic in this region. The analyzed sequences demonstrated a 98% similarity with reference and isolated sequences. In summary, sequencing of the HCV 5' UTR essential for identifying the predominant genotype of HCV RNA in the Sindh region Further research on the distribution of HCV genotypes in other regions of Pakistan could aid in improving screening processes, identifying more effective treatment options, and developing suitable prevention strategies.

Clinical diagnosis of viral hepatitis: Current status and future strategies

Viral hepatitis (VH) is a significant public health issue with tremendous potential to aggravate into chronic hepatitis, liver cirrhosis, and hepatocellular carcinoma. Recent decade has witnessed remarkable uprising in the drug development and effective treatment of VH. An upsurge is seen in identification of antiviral therapies with low rates of viral resistance, the improvement of Hepatitis B Virus (HBV) vaccination and the development of direct-acting antivirals for Hepatitis C Virus (HCV). But unfortunately, the "2030 worldwide eradication" objective of World Health Organization (WHO) is still unmet. It can be largely attributed to the deficit faced by the healthcare system concerning screening and diagnosis. A timely, accurate and comprehensive screening; encompassing maximum population coverage is essential to combat this disease. However, advancements in VH diagnostics remain inadequate and with a marginal use in routine practice. This paper deliberates upon the lacunae in traditional and prevailing diagnostic methodology of viral hepatitis, especially their inadequacy in meeting the unique situations prevailing low- and middle-income countries (LMIC).

Hepatitis C virus transmission between eight high-income countries among men who have sex with men: a whole-genome analysis

BACKGROUND

Microelimination of the hepatitis C virus (HCV) among men who have sex with men (MSM) could be complicated by continuous external introductions and the emergence of phylogenetic clusters harbouring clinically significant resistance-associated substitutions (RAS). To investigate international clustering and the prevalence and transmission of RAS, we aimed to analyse whole-genome HCV sequences from MSM with a recently acquired infection who participated in a large, international HCV treatment trial.

METHODS

For this whole-genome analysis, we obtained HCV sequences from 128 MSM who had acquired HCV within the past 12 months and were participating in the REACT trial. The participants from whom sequences were obtained were recruited at 24 sites in eight countries. We inferred maximum-likelihood phylogenies and identified transmission clusters for HCV genotypes separately. We constructed time-scaled phylogenies to estimate cluster introduction dates and used a Bayesian Skygrid approach to estimate the effective population size over the past 50 years. We calculated the prevalence of RAS and the extent of RAS transmission in the study population.

FINDINGS

The majority of recent HCV infections were part of international networks that arose in the late 1990s and early 2000s. Sequences obtained in the same country clustered frequently, and in 36% of subclusters since 2015 we found evidence of international transmission. European MSM were more likely than non-European MSM to be in a cluster (odds ratio 11·9 [95% CI 3·6-43·4], p<0·0001). The effective population size decreased rapidly since around 2015 in Europe. RAS associated with substantially diminished cure rates were infrequently detected and transmission of highly resistant viruses was not observed.

INTERPRETATION

Despite antiviral treatment becoming widely available, international transmission of HCV among MSM has still occurred over the past 8 years, which could complicate microelimination of the virus in this population. RAS-enriched clusters and widespread RAS transmission are currently not a threat to elimination goals. These findings support an international approach for HCV microelimination among MSM.

FUNDING

National Institutes of Health and Dr. C.J. Vaillant Fonds.

Multiplex MinION sequencing suggests enteric adenovirus F41 genetic diversity comparable to pre-COVID-19 era

Human adenovirus F41 causes acute gastroenteritis in children, and has recently been associated with an apparent increase in paediatric hepatitis of unknown aetiology in the UK, with further cases reported in multiple countries. Relatively little is known about the genetic diversity of adenovirus F41 in UK children; and it is unclear what, if any, impact the COVID-19 pandemic has had on viral diversity in the UK. Methods that allow F41 to be sequenced from clinical samples without the need for viral culture are required to provide the genomic data to address these questions. Therefore, we evaluated an overlapping-amplicon method of sequencing adenovirus genomes from clinical samples using Oxford Nanopore technology. We applied this method to a small sample of adenovirus-species-F-positive extracts collected as part of standard care in the East of England region in January-May 2022. This method produced genomes with >75 % coverage in 13/22 samples and >50 % coverage in 19/22 samples. We identified two F41 lineages present in paediatric patients in the East of England in 2022. Where F41 genomes from paediatric hepatitis cases were available (n=2), these genomes fell within the diversity of F41 from the UK and continental Europe sequenced before and after the 2020-2021 phase of the COVID-19 pandemic. Our analyses suggest that overlapping amplicon sequencing is an appropriate method for generating F41 genomic data from high-virus-load clinical samples, and currently circulating F41 viral lineages were present in the UK and Europe before the COVID-19 pandemic.

Independent evolution of multi-dominant viral genome species observed in a hepatitis C virus carrier

The viral genome quasispecies composition of hepatitis C virus (HCV) could have important implications to viral pathogenesis and resistance to anti-viral treatment. The purpose of the present study was to profile the HCV RNA quasispecies. We developed a strategy to determine the full-length HCV genome sequences co-existing within a single patient serum by using next-generation sequencing technologies. The isolated viral clones were divided into the groups that can be distinguished by core amino acid 70 substitution. Subsequently, we determined HCV full-length genome sequences of three independent dominant species co-existing in the sequential serum with a 7-year interval. From phylogenetic analysis, these dominant species evolved independently. Our study demonstrated that multiple dominant species co-existed in patient sera and evolved independently.

•

HCV RNA genome forms quasispecies which may contribute viral pathogenesis.

•

A strategy was established to determine the full-length HCV genome sequences co-existing within the sera of a single patient.

•

Multiple dominant viral species co-existed in patient sera and evolved independently.

Quasispecies Fitness Partition to Characterize the Molecular Status of a Viral Population. Negative Effect of Early Ribavirin Discontinuation in a Chronically Infected HEV Patient

The changes occurring in viral quasispecies populations during infection have been monitored using diversity indices, nucleotide diversity, and several other indices to summarize the quasispecies structure in a single value. In this study, we present a method to partition quasispecies haplotypes into four fractions according to their fitness: the master haplotype, rare haplotypes at two levels (those present at <0.1%, and those at 0.1−1%), and a fourth fraction that we term emerging haplotypes, present at frequencies >1%, but less than that of the master haplotype. We propose that by determining the changes occurring in the volume of the four quasispecies fitness fractions together with those of the Hill number profile we will be able to visualize and analyze the molecular changes in the composition of a quasispecies with time. To develop this concept, we used three data sets: a technical clone of the complete SARS-CoV-2 spike gene, a subset of data previously used in a study of rare haplotypes, and data from a clinical follow-up study of a patient chronically infected with HEV and treated with ribavirin. The viral response to ribavirin mutagenic treatment was selection of a rich set of synonymous haplotypes. The mutation spectrum was very complex at the nucleotide level, but at the protein (phenotypic/functional) level the pattern differed, showing a highly prevalent master phenotype. We discuss the putative implications of this observation in relation to mutagenic antiviral treatment.

Analysis of direct-acting antiviral-resistant hepatitis C virus haplotype diversity by single-molecule and long-read sequencing

The method of analyzing individual resistant hepatitis C virus (HCV) by a combination of haplotyping and resistance-associated substitution (RAS) has not been fully elucidated because conventional sequencing has only yielded short and fragmented viral genomes. We performed haplotype analysis of HCV mutations in 12 asunaprevir/daclatasvir treatment-failure cases using the Oxford Nanopore sequencer. This enabled single-molecule long-read sequencing using rolling circle amplification (RCA) for correction of the sequencing error. RCA of the circularized reverse-transcription polymerase chain reaction products successfully produced DNA longer than 30 kilobase pairs (kb) containing multiple tandem repeats of a target 3 kb HCV genome. The long-read sequencing of these RCA products could determine the original sequence of the target single molecule as the consensus nucleotide sequence of the tandem repeats and revealed the presence of multiple viral haplotypes with the combination of various mutations in each host. In addition to already known signature RASs, such as NS3-D168 and NS5A-L31/Y93, there were various RASs specific to a different haplotype after treatment failure. The distribution of viral haplotype changed over time; some haplotypes disappeared without acquiring resistant mutations, and other haplotypes, which were not observed before treatment, appeared after treatment. Conclusion: The combination of various mutations other than the known signature RAS was suggested to influence the kinetics of individual HCV quasispecies in the direct-acting antiviral treatment. HCV haplotype dynamic analysis will provide novel information on the role of HCV diversity within the host, which will be useful for elucidating the pathological mechanism of HCV-related diseases.

Emergence of resistance against direct acting antivirals in chronic HCV patients: A real-world study

Interferon/Ribavirin therapy has been replaced by Direct Acting Antivirals (DAAs) due to emergence of Resistance Associated Variants (RAVs) and decrease Sustain Virologic Response (SVR). Current study investigated treatment response of Sofosbuvir and Ribavirin in chronic HCV patients. Total 256 HCV patients with genotype 1a, 2 and 3a received sofosbuvir/ribavirin according to international standards. HCV RNA presence in serum was used as marker for end treatment response (ETR) and sustain virologic response after 24 weeks of treatment (SVR24) in each case. Response to treatment with SOF + RBV was found statistically significant among different HCV genotypes (GT) as out of 47 HCV GT1 patients 42(89.36%) resulted into good ETR but 4(9.52%) of these relapsed and 5(10.63%) led into virologic failure. 5(100%) HCV GT2 patients resulted into SVR24 whereas, out of 204 HCV GT3 patients 194(95.69%) achieved good ETR however, 8(4.12%) of these relapsed and 10(4.90%) resulted in to virologic failure. Efficacy of therapy was found non-significant in treatment naïve and treatment experienced patients as in this study out of 145 treatment naïve patients 139(95.86%) achieved good ETR where 4(2.87%) relapsed while 6(4.13%) led into virologic break through on the other hand among 111 treatment experienced patients 102(91.89%) resulted into good ETR but 8(7.84%) relapsed whereas 9(8.10%) lead into virologic failure. Current study also propose that various liver and spleen complications/liver cirrhosis are related to response of HCV patients to SOF + RBV therapy whereas, variables like old age, gender is not compromising treatment response to DAAs therapy. Various mild side effects encountered by patients during treatment were fatigue, insomnia, headache, nausea, burning body, diarrhea, cough. Overall, this study reported 89.45% efficacy of SOF + RBV regime in chronic HCV Pakistani patients. Current study suggests hunting for possible reasons of resistance so that SOF + RBV therapy may not share the same fortune as previous therapies in near future.

Evolution of within-host variants of the hepatitis C virus

INTRODUCTION

Comprehensive investigation of the within-host evolution of hepatitis C virus (HCV) variants has been difficult without high coverage deep sequencing data and bioinformatics tools to characterise these variants. With the advent of high throughput, long-read sequencing platforms such as Oxford Nanopore Technology (ONT), capturing within-host evolution of HCV using full genome sequences has become feasible. This study aimed to provide the proof of concept that within-host evolutionary analysis of HCV using near-full-length genomes, is achievable.

METHODS

Five treatment naïve subjects with chronic HCV infection were sampled longitudinally from 6 months to 5 years post-infection, with 3-5 sampling timepoints per subject. Near full-length sequences generated using the ONT platform encompassing within-host HCV variants were analysed using an in-house bioinformatic tool. A 200-sequence proxy alignment of the viral variants was made for each subject and timepoint, proportionately representing the observed within-host variants. This alignment was then used in a Bayesian evolutionary analysis using BEAST software suite (v1.8).

RESULTS

The estimated within-host substitution rates ranged between 0.89 and 6.19 × 10^-5 substitutions/site/day. At most timepoints, observed viral lineages were closely related to those from the immediately preceding timepoint, and genetic diversity bottlenecks were observed at intervals in both the acute and chronic phases of infection. The highest within-host mutation rates were observed in the Envelope-P7 and NS5 regions while the Core region was the most conserved.

CONCLUSION

This study demonstrates the feasibility of studying within-host evolution of near-full-length HCV genomes, using long-read sequencing platforms. When considered in conjunction with meta-data such as the host immune response, these methods may offer high resolution insights into immune escape (in vivo or in vitro) to inform vaccine design and to predict spontaneous clearance.

Long-Reads-Based Metagenomics in Clinical Diagnosis With a Special Focus on Fungal Infections

Identification of the causative infectious agent is essential in the management of infectious diseases, with the ideal diagnostic method being rapid, accurate, and informative, while remaining cost-effective. Traditional diagnostic techniques rely on culturing and cell propagation to isolate and identify the causative pathogen. These techniques are limited by the ability and the time required to grow or propagate an agent in vitro and the facts that identification based on morphological traits are non-specific, insensitive, and reliant on technical expertise. The evolution of next-generation sequencing has revolutionized genomic studies to generate more data at a cheaper cost. These are divided into short- and long-read sequencing technologies, depending on the length of reads generated during sequencing runs. Long-read sequencing also called third-generation sequencing emerged commercially through the instruments released by Pacific Biosciences and Oxford Nanopore Technologies, although relying on different sequencing chemistries, with the first one being more accurate both platforms can generate ultra-long sequence reads. Long-read sequencing is capable of entirely spanning previously established genomic identification regions or potentially small whole genomes, drastically improving the accuracy of the identification of pathogens directly from clinical samples. Long-read sequencing may also provide additional important clinical information, such as antimicrobial resistance profiles and epidemiological data from a single sequencing run. While initial applications of long-read sequencing in clinical diagnosis showed that it could be a promising diagnostic technique, it also has highlighted the need for further optimization. In this review, we show the potential long-read sequencing has in clinical diagnosis of fungal infections and discuss the pros and cons of its implementation.

Composition and Dynamics of H1N1 and H7N9 Influenza A Virus Quasispecies in a Co-infected Patient Analyzed by Single Molecule Sequencing Technology

A human co-infected with H1N1 and H7N9 subtypes influenza A virus (IAV) causes a complex infectious disease. The identification of molecular-level variations in composition and dynamics of IAV quasispecies will help to understand the pathogenesis and provide guidance for precision medicine treatment. In this study, using single-molecule real-time sequencing (SMRT) technology, we successfully acquired full-length IAV genomic sequences and quantified their genotypes abundance in serial samples from an 81-year-old male co-infected with H1N1 and H7N9 subtypes IAV. A total of 26 high diversity nucleotide loci was detected, in which the A-G base transversion was the most abundant substitution type (67 and 64%, in H1N1 and H7N9, respectively). Seven significant amino acid variations were detected, such as NA:H275Y and HA: R222K in H1N1 as well as PB2:E627K and NA: K432E in H7N9, which are related to viral drug-resistance or mammalian adaptation. Furtherly, we retrieved 25 H1N1 and 22 H7N9 genomic segment haplotypes from the eight samples based on combining high-diversity nucleotide loci, which provided a more concise overview of viral quasispecies composition and dynamics. Our approach promotes the popularization of viral quasispecies analysis in a complex infectious disease, which will boost the understanding of viral infections, pathogenesis, evolution, and precision medicine.

Assessing hepatitis C virus distribution among vulnerable populations in London using whole genome sequencing: results from the TB-REACH study

Background: Injecting drugs substantially increases the risk of hepatitis C virus (HCV) infection and is common in vulnerable population groups, such as the homeless and prisoners. Capturing accurate data on relative genotype distribution within these groups is essential to inform strategies to reduce HCV transmission. The aim of this study was to utilise a next-generation whole-genome sequencing method recently validated by Public Health England, in order to produce near complete HCV genomes. Methods: In total, 98 HCV positive patients were recruited from homeless hostels and drug treatment services through the National Health Services (NHS) Find and Treat (F&T) Service between May 2011 and June 2013 in London, UK. Samples were sequenced by Next-generation sequencing, with 88 complete HCV genomes constructed by a de novo assembly pipeline. They were analysed phylogenetically for an estimate of their genetic distance. Results: Of the 88 complete HCV genomes, 50/88 (56.8%) were genotype 1; 32/88 (36.4%) genotype 3; 4/88 (4.5%) genotype 2; and 1/88 (1.1%) for genotypes 4 and 6 each. Subtype 1a had the highest number of samples (51.1%), followed by subtype 3a (35.2%), 1b (5.7%), and 2b (3.4%). Samples collected from drug treatment services had the highest number of genotype 1 (69%); genotypes 4 and 6 were only found from samples collected in homeless shelters. Small clusters of highly related genomic sequences were observed both across and within the vulnerable groups sampled. Conclusions: Subsequent phylogenetic analysis provides a first indication that there are related HCV sequences amongst the three vulnerable population groups, reflecting their overlapping social behaviours. This study is the first presentation of whole genome HCV sequences from such vulnerable groups in London and paves the way for similar research in the future.

Validation of a Genotype-Independent Hepatitis C Virus Near-Whole Genome Sequencing Assay

Despite the effectiveness of direct-acting antiviral agents in treating hepatitis C virus (HCV), cases of treatment failure have been associated with the emergence of resistance-associated substitutions. To better guide clinical decision-making, we developed and validated a near-whole-genome HCV genotype-independent next-generation sequencing strategy. HCV genotype 1-6 samples from direct-acting antiviral agent treatment-naïve and -treated HCV-infected individuals were included. Viral RNA was extracted using a NucliSens easyMAG and amplified using nested reverse transcription-polymerase chain reaction. Libraries were prepared using Nextera XT and sequenced on the Illumina MiSeq sequencing platform. Data were processed by an in-house pipeline (MiCall). Nucleotide consensus sequences were aligned to reference strain sequences for resistance-associated substitution identification and compared to NS3, NS5a, and NS5b sequence data obtained from a validated in-house assay optimized for HCV genotype 1. Sequencing success rates (defined as achieving >100-fold read coverage) approaching 90% were observed for most genotypes in samples with a viral load >5 log₁₀ IU/mL. This genotype-independent sequencing method resulted in >99.8% nucleotide concordance with the genotype 1-optimized method, and 100% agreement in genotype assignment with paired line probe assay-based genotypes. The assay demonstrated high intra-run repeatability and inter-run reproducibility at detecting substitutions above 2% prevalence. This study highlights the performance of a freely available laboratory and bioinformatic approach for reliable HCV genotyping and resistance-associated substitution detection regardless of genotype.

Adaptation of Oxford Nanopore technology for hepatitis C whole genome sequencing and identification of within-host viral variants

Background

Hepatitis C (HCV) and many other RNA viruses exist as rapidly mutating quasi-species populations in a single infected host. High throughput characterization of full genome, within-host variants is still not possible despite advances in next generation sequencing. This limitation constrains viral genomic studies that depend on accurate identification of hemi-genome or whole genome, within-host variants, especially those occurring at low frequencies. With the advent of third generation long read sequencing technologies, including Oxford Nanopore Technology (ONT) and PacBio platforms, this problem is potentially surmountable. ONT is particularly attractive in this regard due to the portable nature of the MinION sequencer, which makes real-time sequencing in remote and resource-limited locations possible. However, this technology (termed here ‘nanopore sequencing’) has a comparatively high technical error rate. The present study aimed to assess the utility, accuracy and cost-effectiveness of nanopore sequencing for HCV genomes. We also introduce a new bioinformatics tool (Nano-Q) to differentiate within-host variants from nanopore sequencing.

Results

The Nanopore platform, when the coverage exceeded 300 reads, generated comparable consensus sequences to Illumina sequencing. Using HCV Envelope plasmids (~ 1800 nt) mixed in known proportions, the capacity of nanopore sequencing to reliably identify variants with an abundance as low as 0.1% was demonstrated, provided the autologous reference sequence was available to identify the matching reads. Successful pooling and nanopore sequencing of 52 samples from patients with HCV infection demonstrated its cost effectiveness (AUD$ 43 per sample with nanopore sequencing versus $100 with paired-end short read technology). The Nano-Q tool successfully separated between-host sequences, including those from the same subtype, by bulk sorting and phylogenetic clustering without an autologous reference sequence (using only a subtype-specific generic reference). The pipeline also identified within-host viral variants and their abundance when the parameters were appropriately adjusted.

Conclusion

Cost effective HCV whole genome sequencing and within-host variant identification without haplotype reconstruction are potential advantages of nanopore sequencing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07460-1.

Full-length genome sequencing of RNA viruses-How the approach can enlighten us on hepatitis C and hepatitis E viruses

Among the five main viruses responsible for human hepatitis, hepatitis C virus (HCV) and hepatitis E virus (HEV) are different while sharing similarities. Both viruses can be transmitted by blood or derivatives whereas HEV can also follow environmental or zoonotic routes. These highly variable RNA viruses can cause chronic hepatitis potentially leading to hepatocarcinoma. HCV and HEV can develop new structures and functions under selective pressure to adapt to host immunity, human tissues, treatments or even various animal reservoirs. Elsewhere, with directly acting antiviral treatments, HCV can be eradicated whereas HEV is an emerging pathogen against which specific treatments have to be improved. As a unique molecular tool able to explore viral genomic plasticity, full-length genome (FLG) sequencing has become easier, faster and cheaper. The present review will show how FLG sequencing can explore these RNA viruses with the aim to investigate key genomics data to improve basic knowledge, patients' healthcare and preventive tools.

A machine learning-based treatment prediction model using whole genome variants of hepatitis C virus

In recent years, the development of diagnostics using artificial intelligence (AI) has been remarkable. AI algorithms can go beyond human reasoning and build diagnostic models from a number of complex combinations. Using next-generation sequencing technology, we identified hepatitis C virus (HCV) variants resistant to directing-acting antivirals (DAA) by whole genome sequencing of full-length HCV genomes, and applied these variants to various machine-learning algorithms to evaluate a preliminary predictive model. HCV genomic RNA was extracted from serum from 173 patients (109 with subsequent sustained virological response [SVR] and 64 without) before DAA treatment. HCV genomes from the 109 SVR and 64 non-SVR patients were randomly divided into a training data set (57 SVR and 29 non-SVR) and a validation-data set (52 SVR and 35 non-SVR). The training data set was subject to nine machine-learning algorithms selected to identify the optimized combination of functional variants in relation to SVR status following DAA therapy. Subsequently, the prediction model was tested by the validation-data set. The most accurate learning method was the support vector machine (SVM) algorithm (validation accuracy, 0.95; kappa statistic, 0.90; F-value, 0.94). The second-most accurate learning algorithm was Multi-layer perceptron. Unfortunately, Decision Tree, and Naive Bayes algorithms could not be fitted with our data set due to low accuracy (< 0.8). Conclusively, with an accuracy rate of 95.4% in the generalization performance evaluation, SVM was identified as the best algorithm. Analytical methods based on genomic analysis and the construction of a predictive model by machine-learning may be applicable to the selection of the optimal treatment for other viral infections and cancer.

Single molecule, near full-length genome sequencing of dengue virus

Current methods for dengue virus (DENV) genome amplification, amplify parts of the genome in at least 5 overlapping segments and then combine the output to characterize a full genome. This process is laborious, costly and requires at least 10 primers per serotype, thus increasing the likelihood of PCR bias. We introduce an assay to amplify near full-length dengue virus genomes as intact molecules, sequence these amplicons with third generation “nanopore” technology without fragmenting and use the sequence data to differentiate within-host viral variants with a bioinformatics tool (Nano-Q). The new assay successfully generated near full-length amplicons from DENV serotypes 1, 2 and 3 samples which were sequenced with nanopore technology. Consensus DENV sequences generated by nanopore sequencing had over 99.5% pairwise sequence similarity to Illumina generated counterparts provided the coverage was > 100 with both platforms. Maximum likelihood phylogenetic trees generated from nanopore consensus sequences were able to reproduce the exact trees made from Illumina sequencing with a conservative 99% bootstrapping threshold (after 1000 replicates and 10% burn-in). Pairwise genetic distances of within host variants identified from the Nano-Q tool were less than that of between host variants, thus enabling the phylogenetic segregation of variants from the same host.

Full-genome analysis of hepatitis C virus in HIV-coinfected hemophiliac Japanese patients

AIM

More than 1400 Japanese hemophiliacs acquired HIV infection around 1983 through contaminated blood products imported from the USA, most of whom also acquired hepatitis C virus (HCV) infection. To delineate the HCV genetic relations in HIV-coinfected hemophiliacs, we analyzed stocked plasma samples of the patients seen at the largest referral center for HIV care in Japan.

METHODS

Hepatitis C virus full-genome sequences were amplified and determined using next-generation sequencing, and genotyping and phylogenetic analyses of these sequences were carried out. The results of these hemophiliacs were compared with those of previously studied HIV-coinfected Japanese non-hemophiliacs who had undergone similar analysis of HCV full-genome sequences.

RESULTS

From 1997 to the end of 2017, 72 HIV-infected Japanese hemophiliacs regularly visited our outpatient clinic. Of these, 51 patients had detectable plasma HCV-RNA. The HCV full genome was successfully amplified and sequenced in 50 patients. Not only HCV genotypes 1b (28%) and 2a (6%), which are common in Japan, but also HCV genotypes 1a (32%) and 3a (22%) were identified at high frequency. A single case of intergenotypic recombinant form (2b/1a) and a single case of mixed infection (1a and 3a) were also identified. Each sequence derived from hemophiliacs was more than 0.05 genetic distance away from the other sequences in phylogenetic analysis.

CONCLUSIONS

Various HCV genotypes were identified in Japanese hemophiliacs, a finding that reflects the HCV genotypic distribution in the USA. The genetic distance among them are the results of viral evolution in each patient plus HCV genetic diversity in the USA.

Substitution of the CD81 Binding Site and β-Sandwich Area in E2 of HCV in Cambodia

The high genetic variability of hepatitis C virus (HCV) is the main obstacle to developing a vaccine. E2 has attracted attention for vaccine development because targeting this protein could potentially overcome issues related to the genetic diversity of HCV. In this study, we analyzed HCV genes in the general population of Cambodia and investigated the E2 locus as a candidate for vaccine development. HCV sero-epidemiological surveys were conducted between the period 2010 and 2014, with an HCV RNA–positive rate of 1.3% (11/868). Follow-up blood samples were collected from four anti-HCV– and HCV RNA– positive patients (genotype 1b: 2 cases, 6e: 1 case, 6r: 1 case) after 4.12 years. Analysis of HCV full-length nucleotide sequences in paired specimens revealed that the mutation rates of HCV genotypes 1b and 6e/6r were 1.61–2.03 × 10⁻³ and 2.52–2.74 × 10⁻³ substitutions/site/year, respectively. Non-synonymous substitutions were detected in HVR1, the front layer of the CD81 binding site, and the β-sandwich, but not in the N-terminal region or adjacent to the CD81 binding site. Therefore, we conclude that the CD81 binding site is a promising locus for HCV vaccine development.

An integrated software for virus community sequencing data analysis

BACKGROUND

A virus community is the spectrum of viral strains populating an infected host, which plays a key role in pathogenesis and therapy response in viral infectious diseases. However automatic and dedicated pipeline for interpreting virus community sequencing data has not been developed yet.

RESULTS

We developed Quasispecies Analysis Package (QAP), an integrated software platform to address the problems associated with making biological interpretations from massive viral population sequencing data. QAP provides quantitative insight into virus ecology by first introducing the definition "virus OTU" and supports a wide range of viral community analyses and results visualizations. Various forms of QAP were developed in consideration of broader users, including a command line, a graphical user interface and a web server. Utilities of QAP were thoroughly evaluated with high-throughput sequencing data from hepatitis B virus, hepatitis C virus, influenza virus and human immunodeficiency virus, and the results showed highly accurate viral quasispecies characteristics related to biological phenotypes.

CONCLUSIONS

QAP provides a complete solution for virus community high throughput sequencing data analysis, and it would facilitate the easy analysis of virus quasispecies in clinical applications.

B cell immunodominance in primary hepatitis C virus infection

BACKGROUND & AIMS

Neutralising antibodies (NAbs) play a key role in clearance of HCV. NAbs have been isolated and mapped to several domains on the HCV envelope proteins. However, the immunodominance of these epitopes in HCV infection remains unknown, hindering efforts to elicit optimal epitope-specific responses. Furthermore, it remains unclear which epitope-specific responses are associated with broad NAb (bNAb) activity in primary HCV infection. The aim of this study was to define B cell immunodominance in primary HCV, and its implications on neutralisation breadth and clearance.

METHODS

Using samples from 168 patients with primary HCV infection, the antibody responses targeted 2 immunodominant domains, termed domains B and C. Genotype 1 and 3 infections were associated with responses targeted towards different bNAb domains.

RESULTS

No epitopes were uniquely targeted by clearers compared to those who developed chronic infection. Samples with bNAb activity were enriched for multi-specific responses directed towards the epitopes antigenic region 3, antigenic region 4, and domain D, and did not target non-neutralising domains.

CONCLUSIONS

This study outlines for the first time a clear NAb immunodominance profile in primary HCV infection, and indicates that it is influenced by the infecting virus. It also highlights the need for a vaccination strategy to induce multi-specific responses that do not target non-neutralising domains.

LAY SUMMARY

Neutralising antibodies will likely form a key component of a protective hepatitis C virus vaccine. In this work we characterise the predominant neutralising and non-neutralising antibody (epitope) targets in acute hepatitis C virus infection. We have defined the natural hierarchy of epitope immunodominance, and demonstrated that viral genotype can impact on this hierarchy. Our findings highlight key epitopes that are associated with broadly neutralising antibodies, and the deleterious impact of mounting a response towards some of these domains on neutralising breadth. These findings should guide future efforts to design immunogens aimed at generating neutralising antibodies with a vaccine candidate.

Full-Genome Analysis of Hepatitis C Virus in Japanese and Non-Japanese Patients Coinfected With HIV-1 in Tokyo

BACKGROUND

Acute hepatitis C virus (HCV) infection is increasing among HIV-1-infected individuals in Tokyo. Appropriate clinical management is needed.

SETTING

To delineate the epidemiological status of HCV transmission, we analyzed stocked plasma samples of HCV/HIV-1-coinfected patients seen at the largest referral center for HIV care in Tokyo.

METHODS

HCV full-genome sequences were amplified and determined using next-generation sequencing. HCV genotyping and phylogenetic and phylodynamic analyses of thus obtained sequences were performed and combined with the analysis of HIV-1 reverse transcriptase sequences.

RESULTS

HCV phylogenetic analysis identified 3 dense clusters containing cases of men who have sex with men (MSM) and injection drug users (IDUs). Most of the confirmed acute infection cases were included within these clusters, indicating that the clustered viruses are currently being actively transmitted among HIV-1-infected MSM and IDU. Phylodynamic analysis indicated population expansion of one of these clusters from 2006 to 2008, during which the largest number of HIV-1-infected MSM was diagnosed in Tokyo. HIV-1 reverse transcriptase sequences of HCV-coinfected patients included in the same clusters did not converge together and did not form clusters, but rather diverged in the area of subtype B in the phylogenetic tree, indicating that they acquired HCV infection from individuals different from those from whom they had acquired HIV-1 infection. It is considered that these MSM changed their sexual partners and that IDU changed their drug use groups.

CONCLUSIONS

The results warrant careful monitoring of high-risk groups including MSM and IDU and early introduction of HCV treatment to prevent HCV epidemic.

Molecular characterization of hepatitis C virus in liver disease patients in Botswana: a retrospective cross-sectional study

Background

Hepatitis C virus (HCV) infection is a major cause of chronic liver disease globally. Direct acting antivirals (DAAs) have proven effective in curing HCV. However, the current standard of care (SOC) in Botswana remains PEGylated interferon-α (IFN-α) with ribavirin. Several mutations have been reported to confer resistance to interferon-based treatments. Therefore, there is a need to determine HCV genotypes in Botswana, as these data will guide new treatment guidelines and understanding of HCV epidemiology in Botswana.

Methods

This was a retrospective cross-sectional pilot study utilizing plasma obtained from 55 participants from Princess Marina Hospital in Gaborone, Botswana. The partial core region of HCV was amplified, and genotypes were determined using phylogenetic analysis.

Results

Four genotype 5a and two genotype 4v sequences were identified. Two significant mutations – K10Q and R70Q – were observed in genotype 5a sequences and have been associated with increased risk of hepatocellular carcinoma (HCC), while R70Q confers resistance to interferon-based treatments.

Conclusion

Genotypes 5a and 4v are circulating in Botswana. The presence of mutations in genotype 5 suggests that some patients may not respond to IFN-based regimens. The information obtained in this study, in addition to the World health organization (WHO) recommendations, can be utilized by policy makers to implement DAAs as the new SOC for HCV treatment in Botswana.

Ribavirin induces hepatitis C virus genome mutations in chronic hepatitis patients who failed to respond to prior daclatasvir plus asunaprevir therapy

Ribavirin (RBV) induces nucleotide (nt) substitutions in hepatitis C virus (HCV) genome nonstructural (NS) regions. Although emergence of drug resistance-associated variants is associated with direct-acting antiviral treatment failure, the effect of RBV on genome substitutions in such patients is unknown. Genotype 1b HCV subgenomic replicon cells were treated with RBV for 120 hours. Six patients with chronic genotype 1b with HCV-infected patients who failed to respond to prior daclatasvir plus asunaprevir (DCV/ASV) therapy were treated with 12 weeks of sofosbuvir and ledipasvir plus RBV after 4 weeks of RBV monotherapy. RBV-induced genome mutations in the HCV NS region (nt3493-9301) in replicon cells and in patients during 4 weeks of RBV monotherapy were analyzed by deep sequencing. RBV-associated G-to-A and C-to-U transitions increased in a dose-dependent manner in HCV replicon cells after the RBV treatment. In patients with prior DCV/ASV treatment failures, the median serum HCV RNA level was 6.25 ± 0.31 log IU/mL at the start of RBV therapy and decreased significantly to 5.95 ± 0.4 log IU/mL (P = .03) after 4 weeks of RBV monotherapy. Although predominant HCV genome substitutions rates were similar between nontreatment and RBV-treatment periods (0.042 and 0.031 per base pair, respectively; P  = .248), the frequencies of G-to-A and C-to-U transitions significantly increased after RBV monotherapy. These transitions were enriched, particularly within the HCV NS3 region in all patients. RBV treatment induces G-to-A and C-to-U transitions in the HCV genome even in chronic patients with hepatitis C with prior DCV/ASV treatment failures.

Clearance of hepatitis C virus is associated with early and potent but narrowly-directed, Envelope-specific antibodies

Hepatitis C virus (HCV) is one of very few viruses that are either naturally cleared, or alternatively persist to cause chronic disease. Viral diversity and escape, as well as host adaptive immune factors, are believed to control the outcome. To date, there is limited understanding of the critical, early host-pathogen interactions. The asymptomatic nature of early HCV infection generally prevents identification of the transmitted/founder (T/F) virus, and thus the study of host responses directed against the autologous T/F strain. In this study, 14 rare subjects identified from very early in infection (4–45 days) with varied disease outcomes (n = 7 clearers) were examined in regard to the timing, breadth, and magnitude of the neutralizing antibody (nAb) response, as well as evolution of the T/F strain. Clearance was associated with earlier onset and more potent nAb responses appearing at a mean of 71 days post-infection (DPI), but these responses were narrowly directed against the autologous T/F virus or closely related variants. In contrast, a delayed onset of nAbs (mean 425 DPI) was observed in chronic progressors that appear to have targeted longitudinal variants rather than the T/F strain. The nAb responses in the chronic progressors mapped to known CD81 binding epitopes, and were associated with rapid emergence of new viral variants with reduced CD81 binding. We propose that the prolonged period of viremia in the absence of nAbs in these subjects was associated with an increase in viral diversity, affording the virus greater options to escape nAb pressure once it emerged. These findings indicate that timing of the nAb response is essential for clearance. Further investigation of the specificities of the early nAbs and the factors regulating early induction of protective nAbs is needed.

Exploring the hepatitis C virus genome using single molecule real-time sequencing

Single molecular real-time (SMRT) sequencing, also called third-generation sequencing, is a novel sequencing technique capable of generating extremely long contiguous sequence reads. While conventional short-read sequencing cannot evaluate the linkage of nucleotide substitutions distant from one another, SMRT sequencing can directly demonstrate linkage of nucleotide changes over a span of more than 20 kbp, and thus can be applied to directly examine the haplotypes of viruses or bacteria whose genome structures are changing in real time. In addition, an error correction method (circular consensus sequencing) has been established and repeated sequencing of a single-molecule DNA template can result in extremely high accuracy. The advantages of long read sequencing enable accurate determination of the haplotypes of individual viral clones. SMRT sequencing has been applied in various studies of viral genomes including determination of the full-length contiguous genome sequence of hepatitis C virus (HCV), targeted deep sequencing of the HCV NS5A gene, and assessment of heterogeneity among viral populations. Recently, the emergence of multi-drug resistant HCV viruses has become a significant clinical issue and has been also demonstrated using SMRT sequencing. In this review, we introduce the novel third-generation PacBio RSII/Sequel systems, compare them with conventional next-generation sequencers, and summarize previous studies in which SMRT sequencing technology has been applied for HCV genome analysis. We also refer to another long-read sequencing platform, nanopore sequencing technology, and discuss the advantages, limitations and future perspectives in using these third-generation sequencers for HCV genome analysis.

Single molecule real-time (SMRT) sequencing comes of age: applications and utilities for medical diagnostics

Short read massive parallel sequencing has emerged as a standard diagnostic tool in the medical setting. However, short read technologies have inherent limitations such as GC bias, difficulties mapping to repetitive elements, trouble discriminating paralogous sequences, and difficulties in phasing alleles. Long read single molecule sequencers resolve these obstacles. Moreover, they offer higher consensus accuracies and can detect epigenetic modifications from native DNA. The first commercially available long read single molecule platform was the RS system based on PacBio's single molecule real-time (SMRT) sequencing technology, which has since evolved into their RSII and Sequel systems. Here we capsulize how SMRT sequencing is revolutionizing constitutional, reproductive, cancer, microbial and viral genetic testing.

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.

Efficacy and safety of glecaprevir/pibrentasvir for chronic hepatitis C virus genotypes 1-6 infection: A systematic review and meta-analysis

This systematic review and meta-analysis investigated the efficacy and safety of glecaprevir and pibrentasvir (G/P) for chronic hepatitis C virus (HCV) infection. Pubmed, Embase, Cochrane Library and Scopus were searched to identify relevant studies through August 2018. Data from eligible studies were pooled and sustained virological response rates at 12 weeks' post-treatment (SVR12) were calculated. Thirteen studies with 3082 patients were included and the overall SVR12 rate was 97.8%. The SVR12 rates of subgroups were: G/P 300 mg/120 mg and 200 mg/120 mg: 97.9% and 98.3%; HCV genotype (GT)1, GT2, GT3 and GT4-6: 99.8%, 99.2%, 96.1% and 100%; G/P and G/P plus ribavirin (RBV): 97.9% and 98.2%; G/P (300 mg/120 mg) for 8 weeks, 12 weeks and 16 weeks: 98.8%, 98.5% and 95.6%; treatment-naïve and treatment-experienced patients: 96.7% and 98.3%; patients without and with compensated cirrhosis: 99.4% and 98.8%; patients without and with human immunodeficiency virus (HIV) co-infection: 97.8% and 99.4%; and patients without and with severe renal impairment (SRI): 97.8% and 99.4%. Virological failure and relapse and serious drug-related adverse events were rare. These results indicate that 8- or 12-week G/P treatment achieved high SVR12 rates in HCV GTs 1-6 patients without or with compensated cirrhosis, with good safety profiles, irrespective of dose, RBV use, treatment-experience, HIV co-infection and renal impairment. Due to the limited number of evaluated patients with GT3 infection, further studies are needed to define optimal treatment duration for GT3 cirrhosis patients and patients with prior treatment experience of direct-acting antivirals.

Genomic variability of within-host hepatitis C variants in acute infection

Interactions between the host immune system and the viral variants determine persistence of hepatitis C virus (HCV) infection after the acute phase of infection. This study describes the genetic variability of within-host HCV viral variants in acute infection and correlates it with host- and virus-related traits and infection outcome. Next generation sequence data (Illumina, MiSeq platform) of viral genomes from 116 incident acute infections (within 180 days of infection) were analysed to determine all the single nucleotide polymorphism (SNP) frequencies above a threshold of 0.1%. The variability of the SNPs for the full open reading frame of the genome as well as for each protein coding region were compared using mean standardized Shannon entropy (SE) values calculated separately for synonymous and nonsynonymous mutations. The envelope glycoproteins regions (E1 and E2) had the highest SE values (indicating greater variability) followed by the NS5B region. Nonsynonymous mutations rather than synonymous mutations were the main contributors to genomic variability in acute infection. The mean difference of Shannon entropy was also compared between subjects after categorizing the samples according to host and virus-related traits. Host IFNL3 allele CC polymorphism at rs12979860 (vs others) and viral genotype 1a (vs 3a) were associated with higher genomic variability across the viral open reading frame. Time since infection, host gender or continent of origin was not associated with the viral genomic variability. Viral genomic variability did not predict spontaneous clearance.

Application of deep sequencing methods for inferring viral population diversity

The first deep sequencing method was announced in 2005. Due to an increasing number of sequencing data and a reduction in the costs of each sequencing dataset, this innovative technique was soon applied to genetic investigations of viral genome diversity in various viruses, particularly RNA viruses. These deep sequencing findings documented viral epidemiology and evolution and provided high-resolution data on the genetic changes in viral populations. Here, we review deep sequencing platforms that have been applied in viral quasispecies studies. Further, we discuss recent deep sequencing studies on viral inter- and intrahost evolution, drug resistance, and humoral immune selection, especially in emerging and re-emerging viruses. Deep sequencing methods are becoming the standard for providing comprehensive results of viral population diversity, and their applications are discussed.

Dynamic interactions between RNA viruses and human hosts unravelled by a decade of next generation sequencing

BACKGROUND

Next generation sequencing (NGS) methods have significantly contributed to a paradigm shift in genomic research for nearly a decade now. These methods have been useful in studying the dynamic interactions between RNA viruses and human hosts.

SCOPE OF THE REVIEW

In this review, we summarise and discuss key applications of NGS in studying the host - pathogen interactions in RNA viral infections of humans with examples.

MAJOR CONCLUSIONS

Use of NGS to study globally relevant RNA viral infections have revolutionized our understanding of the within host and between host evolution of these viruses. These methods have also been useful in clinical decision-making and in guiding biomedical research on vaccine design.

GENERAL SIGNIFICANCE

NGS has been instrumental in viral genomic studies in resolving within-host viral genomic variants and the distribution of nucleotide polymorphisms along the full-length of viral genomes in a high throughput, cost effective manner. In the future, novel advances such as long read, single molecule sequencing of viral genomes and simultaneous sequencing of host and pathogens may become the standard of practice in research and clinical settings. This will also bring on new challenges in big data analysis.

Quasispecies dynamics and clinical significance of hepatitis C virus (HCV) antiviral resistance

Hepatitis C virus (HCV) follows quasispecies dynamics in infected hosts and this influences its biology, how the virus diversifies into several genotypes and many subtypes, and how viral populations respond to antiviral therapies. Despite current antiviral combinations being able to cure a great percentage of HCV-infected patients, the presence of resistance-associated substitutions (RASs) diminishes the success of antiviral therapies, which is a main concern in the re-treatment of patients treated with direct-acting antiviral agents. Current methodologies such as ultra deep sequencing are ideal tools to obtain a detailed representation of the mutant spectrum composition circulating in infected patients. Such knowledge should allow optimisation of rescue treatments. A new mechanism of antiviral resistance not based on the selection of RASs but on high viral fitness is discussed.

Pipeline for specific subtype amplification and drug resistance detection in hepatitis C virus

BACKGROUND

Despite the high sustained virological response rates achieved with current directly-acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 5-10% of treated patients do not respond to current antiviral therapies, and basal resistance to DAAs is increasingly detected among treatment-naïve infected individuals. Identification of amino acid substitutions (including those in minority variants) associated with treatment failure requires analytical designs that take into account the high diversification of HCV in more than 86 subtypes according to the ICTV website (June 2017).

METHODS

The methodology has involved five sequential steps: (i) to design 280 oligonucleotide primers (some including a maximum of three degenerate positions), and of which 120 were tested to amplify NS3, NS5A-, and NS5B-coding regions in a subtype-specific manner, (ii) to define a reference sequence for each subtype, (iii) to perform experimental controls to define a cut-off value for detection of minority amino acids, (iv) to establish bioinformatics' tools to quantify amino acid replacements, and (v) to validate the procedure with patient samples.

RESULTS

A robust ultra-deep sequencing procedure to analyze HCV circulating in serum samples from patients infected with virus that belongs to the ten most prevalent subtypes worldwide: 1a, 1b, 2a, 2b, 2c, 2j, 3a, 4d, 4e, 4f has been developed. Oligonucleotide primers are subtype-specific. A cut-off value of 1% mutant frequency has been established for individual mutations and haplotypes.

CONCLUSION

The methodological pipeline described here is adequate to characterize in-depth mutant spectra of HCV populations, and it provides a tool to understand HCV diversification and treatment failures. The pipeline can be periodically extended in the event of HCV diversification into new genotypes or subtypes, and provides a framework applicable to other RNA viral pathogens, with potential to couple detection of drug-resistant mutations with treatment planning.

Amplification and next generation sequencing of near full-length human enteroviruses for identification and characterisation from clinical samples

More than 100 different enterovirus (EV) genotypes infect humans and contribute to substantial morbidity. However, current methods for characterisation of full-length genomes are based on Sanger sequencing of short genomic regions, which are labour-intensive and do not enable comprehensive characterisation of viral populations. Here, we describe a simple and sensitive protocol for the amplification and sequencing of near full-length genomes of human EV species using next generation sequencing. EV genomes were amplified from 89% of samples tested, with C_t values ranging between 15.7 and 39.3. These samples included 7 EV-A genotypes (CVA2, 5–7, 10, 16 and EV71), 19 EV-B genotypes (CVA9, CVB1-6, ECHO3, 4, 6, 7, 9, 11, 16, 18, 25, 29, 30, and EV69), 3 EV-C genotypes (CVA19 and PV2, 3) and 1 EV-D genotype (EV70). We characterised 70 EVs from 58 clinical stool samples and eight reference strains, with a minimum of 100X depth. We found evidence of co-infection in four clinical specimens, each containing two distinct EV genotypes (CVB3/ECHO7, CVB3/ECHO18 and ECHO9/30). Characterisation of the complete genome provided conclusive genotyping of EVs, which can be applied to investigate the intra-host virus evolution of EVs, and allows further identification and investigation of EV outbreaks.

Long reads: their purpose and place

In recent years long-read technologies have moved from being a niche and specialist field to a point of relative maturity likely to feature frequently in the genomic landscape. Analogous to next generation sequencing, the cost of sequencing using long-read technologies has materially dropped whilst the instrument throughput continues to increase. Together these changes present the prospect of sequencing large numbers of individuals with the aim of fully characterizing genomes at high resolution. In this article, we will endeavour to present an introduction to long-read technologies showing: what long reads are; how they are distinct from short reads; why long reads are useful and how they are being used. We will highlight the recent developments in this field, and the applications and potential of these technologies in medical research, and clinical diagnostics and therapeutics.

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.

Improving eukaryotic genome annotation using single molecule mRNA sequencing

BACKGROUND

The advantages of Pacific Biosciences (PacBio) single-molecule real-time (SMRT) technology include long reads, low systematic bias, and high consensus read accuracy. Here we use these attributes to improve on the genome annotation of the parasitic hookworm Ancylostoma ceylanicum using PacBio RNA-Seq.

RESULTS

We sequenced 192,888 circular consensus sequences (CCS) derived from cDNAs generated using the CloneTech SMARTer system. These SMARTer-SMRT libraries were normalized and size-selected providing a robust population of expressed structural genes for subsequent genome annotation. We demonstrate PacBio mRNA sequences based genome annotation improvement, compared to genome annotation using conventional sequencing-by-synthesis alone, by identifying 1609 (9.2%) new genes, extended the length of 3965 (26.7%) genes and increased the total genomic exon length by 1.9 Mb (12.4%). Non-coding sequence representation (primarily from UTRs based on dT reverse transcription priming) was particularly improved, increasing in total length by fifteen-fold, by increasing both the length and number of UTR exons. In addition, the UTR data provided by these CCS allowed for the identification of a novel SL2 splice leader sequence for A. ceylanicum and an increase in the number and proportion of functionally annotated genes. RNA-seq data also confirmed some of the newly annotated genes and gene features.

CONCLUSION

Overall, PacBio data has supported a significant improvement in gene annotation in this genome, and is an appealing alternative or complementary technique for genome annotation to the other transcript sequencing technologies.

A comparative study on the characterization of hepatitis B virus quasispecies by clone-based sequencing and third-generation sequencing

Hepatitis B virus (HBV) has a high mutation rate due to the extremely high replication rate and the proofreading deficiency during reverse transcription. The generated variants with genetic heterogeneity are described as viral quasispecies (QS). Clone-based sequencing (CBS) is thought to be the ‘gold standard’ for assessing QS complexity and diversity of HBV, but an important issue about CBS is cost-effectiveness and laborious. In this study, we investigated the utility of the third-generation sequencing (TGS) DNA sequencing to characterize genetic heterogeneity of HBV QS and assessed the possible contribution of TGS technology in HBV QS studies. Parallel experiments including 3 control samples, which consisted of HBV full gene genotype B and genotype C plasmids, and 10 patients samples were performed by using CBS and TGS to analyze HBV whole-genome QS. Characterization of QS heterogeneity was conducted by using comprehensive statistical analysis. The results showed that TGS had a high consistency with CBS when measuring the complexity and diversity of QS. In addition, to detect rare variants, there were strong advantages conferred by TGS. In summary, TGS was considered to be practicable in HBV QS studies and it might have a relevant role in the clinical management of HBV infection in the future.

Exploiting Genetic Information to Trace Plant Virus Dispersal in Landscapes

During the past decade, knowledge of pathogen life history has greatly benefited from the advent and development of molecular epidemiology. This branch of epidemiology uses information on pathogen variation at the molecular level to gain insights into a pathogen's niche and evolution and to characterize pathogen dispersal within and between host populations. Here, we review molecular epidemiology approaches that have been developed to trace plant virus dispersal in landscapes. In particular, we highlight how virus molecular epidemiology, nourished with powerful sequencing technologies, can provide novel insights at the crossroads between the blooming fields of landscape genetics, phylogeography, and evolutionary epidemiology. We present existing approaches and their limitations and contributions to the understanding of plant virus epidemiology.

Genetic diversity of hepatitis C virus in Ethiopia

Hepatitis C virus (HCV) is genetically highly divergent and classified in seven major genotypes and approximately hundred subtypes. These genotypes/subtypes have different geographic distribution and response to antiviral therapy. In Ethiopia, however, little is known about their molecular epidemiology and genetic diversity. The aim of this study was to investigate the distribution and genetic diversity of HCV genotypes/subtypes in Ethiopia, using 49 HCV RNA positive samples. HCV genotypes and subtypes were determined based on the sequences of the core and the nonstructural protein 5B (NS5B) genomic regions. Phylogenetic analysis revealed that the predominant was genotype 4 (77.6%) followed by 2 (12.2%), 1 (8.2%), and 5 (2.0%). Seven subtypes were identified (1b, 1c, 2c, 4d, 4l, 4r and 4v), with 4d (34.7%), 4r (34.7%) and 2c (12.2%) as the most frequent subtypes. Consistent with the presence of these subtypes was the identification of a potential recombinant virus. One strain was typed as genotype 2c in the NS5B region sequence and genotype 4d in the core region. In conclusion, genotype 4 HCV viruses, subtypes 4d and 4r, are most prevalent in Ethiopia. This genotype is considered to be difficult to treat, thus, our finding has an important impact on the development of treatment strategies and patient management in Ethiopia.

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.

Dynamic evolution of hepatitis C virus resistance-associated substitutions in the absence of antiviral treatment

Resistance against new hepatitis C virus (HCV) antivirals is an area of increasing interest. Resistance-associated substitutions (RASs) have been identified in treatment-naïve individuals, but pressures driving treatment-independent RAS emergence are poorly understood. We analysed the longitudinal evolution of RASs in twelve participants with early acute HCV infections. Full-genome deep sequences were analysed for changes in RAS frequency within NS3, NS5A and NS5B-coding regions over the course of the infection. Emergence of RASs relevant only to the polymerase non-nucleoside inhibitors (NNI) was detected, and these lay within CD8+ T-cell epitopes. Conversely, the loss of NNI RASs over time appeared likely to be driven by viral fitness constraints. These results highlight the importance of monitoring CD8+ T cell epitope-associated RASs in populations with dominant HLA types.

Use of genome sequencing to identify hepatitis C virus transmission in a renal healthcare setting

BACKGROUND

Hepatitis C virus (HCV) infection is a major health burden worldwide. A patient with no history of HCV infection while on a renal unit was found to seroconvert to HCV.

AIM

To report the use of sequencing to postulate how transmission of HCV occurred in a healthcare setting, and how this guided our outbreak investigation.

FINDINGS

Based on infection control inspections the transmission event was surmised to be due to ward environmental contamination with blood and subsequent inoculation from intravenous interventions on the patient acquiring HCV. We discuss the interventions put in place in response to the outbreak investigation findings.

CONCLUSION

Sequencing of healthcare-acquired HCV infections should be undertaken as routine practice in outbreak investigations.

Phylogenetic analysis of full-length, early infection, hepatitis C virus genomes among people with intravenous drug use: the InC3 Study

Cross-continental phylogenetic analysis is important to understand subtle molecular differences of currently circulating hepatitis C virus (HCV) subtypes. Existence of such differences can be crucial in pursuing a universal hepatitis C vaccine. We characterized molecular epidemiology of early HCV infections identified across nine cohorts [North America (n=4), Australia (n=4) and Europe (n=1)] in the International Collaborative of Incident HIV and Hepatitis C in Injecting Cohorts (InC³ ). One hundred and ninety-two full-length HCV genomes were amplified from plasma of incident infections and subjected to next generation sequencing to establish the largest cross-continental, full-length acute HCV genomic data set available to date. Genomes from the most common subtypes (1a: n=94, 2b: n=15 and 3a: n=68) were used in phylogenetic analysis. Using full genome trees, 78 sequences (44%) were found to lie within 29 phylogenetic clusters/pairs defined on the basis of molecular similarity of consensus sequences. Of these, 26 each had exclusively Australian or North American sequences indicating a strong geographical bias for molecular similarity. On further analysis of behavioural and demographic associations, binary logistic regression analysis showed that older age and non-Caucasian ethnicity were significantly associated with clustering. HCV probably evolves in micro-epidemics within geographically isolated communities.

Analysis of resistance-associated substitutions in acute hepatitis C virus infection by deep sequencing across six genotypes and three continents

Several direct-acting antivirals (DAAs) have been approved for the treatment of chronic hepatitis C virus (HCV) infections, opening the door to highly effective interferon-free treatment regimens. Resistance-associated substitutions (RASs) have been reported both in treatment-naïve patients and following treatment with protease (NS3), phosphoprotein (NS5A) and polymerase (NS5B) inhibitors. The prevalence of naturally occurring RASs in untreated HCV-infected individuals has mostly been analysed in those infected with genotype 1 (GT1), in the late phase of infection, and only within limited regions of the genome. Furthermore, the geographic distribution of RASs remains poorly characterized. In this study, we used next-generation sequencing to analyse full-length HCV genomes for the prevalence of RASs in acute HCV infections identified in nine international prospective cohorts. RASs were analysed in 179 participants infected with all six major HCV genotypes (GT1-GT6), and the geographic distribution of RASs was assessed in 107 GT1a and GT3a samples. While RASs were detected at varied frequencies across the three genomic regions, and between genotypes, RASs relevant to multiple DAAs in the leading IFN-free regimens were rarely detected in combination. Low-frequency RASs (<10% of the viral population) were also shown to have a GT-specific distribution. The main RASs with geographic associations were NS3 Q80K in GT1a samples and NS5B N142T in GT3a. These data provide the backdrop for prospective surveillance of RASs during DAA treatment scale-up.

Understanding the complex evolution of rapidly mutating viruses with deep sequencing: Beyond the analysis of viral diversity

With the advent of affordable deep sequencing technologies, detection of low frequency variants within genetically diverse viral populations can now be achieved with unprecedented depth and efficiency. The high-resolution data provided by next generation sequencing technologies is currently recognised as the gold standard in estimation of viral diversity. In the analysis of rapidly mutating viruses, longitudinal deep sequencing datasets from viral genomes during individual infection episodes, as well as at the epidemiological level during outbreaks, now allow for more sophisticated analyses such as statistical estimates of the impact of complex mutation patterns on the evolution of the viral populations both within and between hosts. These analyses are revealing more accurate descriptions of the evolutionary dynamics that underpin the rapid adaptation of these viruses to the host response, and to drug therapies. This review assesses recent developments in methods and provide informative research examples using deep sequencing data generated from rapidly mutating viruses infecting humans, particularly hepatitis C virus (HCV), human immunodeficiency virus (HIV), Ebola virus and influenza virus, to understand the evolution of viral genomes and to explore the relationship between viral mutations and the host adaptive immune response. Finally, we discuss limitations in current technologies, and future directions that take advantage of publically available large deep sequencing datasets.