Single-molecular real-time deep sequencing reveals the dynamics of multi-drug resistant haplotypes and structural variations in the hepatitis C virus genome

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.

Adaptation of single molecule real time (SMRT) sequence technology for hepatitis C virus genome sequencing and identification of resistance-associated substitutions

BACKGROUND AND OBJECTIVE

Hepatitis C virus (HCV) resistance-associated substitutions (RASs) impact HCV treatment with direct-acting antivirals (DAAs). Therefore, a sensitive sequencing assay for detecting HCV RASs is crucial. PacBio sequencing, a Single molecule real time (SMRT) platform, is capable of long-fragment sequencing. This study aims to assess the prospects of PacBio sequencing by comparing its accuracy with Sanger sequencing at distinct thresholds and its cost-effectiveness with next-generation sequencing (NGS).

METHODS

A total of 28 specimens were selected from the HCV RNA-positive individuals in Linzhou, China. Of these specimens, the NS3 to NS5B regions were amplified and sent for PacBio sequencing. Sequence processing was accomplished by lima, pbmm2 and quasitools software, with threshold setting at 1%, 5%, 10% and 15%.

RESULTS

High-frequency RASs (S122G in NS3, R30Q in NS5A, and C316N in NS5B), rare RASs (V55A, R155P, and S174AT in NS3; Y448H in NS5B), and low-threshold RASs (L31Q and L31QFHY in NS5A; L159F in NS5B) were identified. It was found that the results of HCV RASs at the 10% threshold of PacBio sequencing are comparable to those of Sanger sequencing. In terms of high-throughput sequencing, the price of PacBio sequencing (571 CNY per specimen) is significantly lower than that of NGS (1000 CNY per specimen).

CONCLUSIONS

In summary, these findings suggest that the adoption of the 10% threshold in PacBio sequencing for the analysis of HCV RASs is advisable. Moreover, given its relatively lower cost, PacBio sequencing holds great promise as a valuable tool for large-scale population sequencing.

Easing genomic surveillance: A comprehensive performance evaluation of long-read assemblers across multi-strain mixture data of HIV-1 and Other pathogenic viruses for constructing a user-friendly bioinformatic pipeline

Background

Determining the appropriate computational requirements and software performance is essential for efficient genomic surveillance. The lack of standardized benchmarking complicates software selection, especially with limited resources.

Methods

We developed a containerized benchmarking pipeline to evaluate seven long-read assemblers-Canu, GoldRush, MetaFlye, Strainline, HaploDMF, iGDA, and RVHaplo-for viral haplotype reconstruction, using both simulated and experimental Oxford Nanopore sequencing data of HIV-1 and other viruses. Benchmarking was conducted on three computational systems to assess each assembler's performance, utilizing QUAST and BLASTN for quality assessment.

Results

Our findings show that assembler choice significantly impacts assembly time, with CPU and memory usage having minimal effect. Assembler selection also influences the size of the contigs, with a minimum read length of 2,000 nucleotides required for quality assembly. A 4,000-nucleotide read length improves quality further. Canu was efficient among de novo assemblers but not suitable for multi-strain mixtures, while GoldRush produced only consensus assemblies. Strainline and MetaFlye were suitable for metagenomic sequencing data, with Strainline requiring high memory and MetaFlye operable on low-specification machines. Among reference-based assemblers, iGDA had high error rates, RVHaplo showed the best runtime and accuracy but became ineffective with similar sequences, and HaploDMF, utilizing machine learning, had fewer errors with a slightly longer runtime.

Conclusions

The HIV-64148 pipeline, containerized using Docker, facilitates easy deployment and offers flexibility to select from a range of assemblers to match computational systems or study requirements. This tool aids in genome assembly and provides valuable information on HIV-1 sequences, enhancing viral evolution monitoring and understanding.

Lassa virus NP DEDDh 3'-5' exoribonuclease activity is required for optimal viral RNA replication and mutation control

Lassa virus (LASV), a mammarenavirus from Arenaviridae, is the causative agent of Lassa fever (LF) endemic in West Africa. Currently, there are no vaccines or antivirals approved for LF. The RNA-dependent RNA polymerases (RdRp) of RNA viruses are error-prone. As a negative-sense RNA virus, how LASV copes with errors in RNA synthesis and ensures optimal RNA replication are not well elucidated. LASV nucleoprotein (NP) contains a DEDDH 3'-to-5' exoribonuclease motif (ExoN), which is known to be essential for LASV evasion of the interferon response via its ability to degrade virus-derived double-stranded RNA. Herein, we present evidence that LASV NP ExoN has an additional function important for viral RNA replication. We rescued an ExoN-deficient LASV mutant (ExoN- rLASV) by using a reverse genetics system. Our data indicated that abrogation of NP ExoN led to impaired LASV growth and RNA replication in interferon-deficient cells as compared with wild-type rLASV. By utilizing PacBio Single Molecule, Real-Time (SMRT) long-read sequencing technology, we found that rLASV lacking ExoN activity was prone to producing aberrant viral genomic RNA with structural variations. In addition, NP ExoN deficiency enhanced LASV sensitivity to mutagenic nucleoside analogues in virus titration assay. Next-generation deep sequencing analysis showed increased single nucleotide substitution in ExoN- LASV RNA following mutagenic 5-flurouracil treatment. In conclusion, our study revealed that LASV NP ExoN is required for efficient viral RNA replication and mutation control. Among negative-sense RNA viruses, LASV NP is the first example that a viral protein, other than the RdRp, contributes to reduce errors in RNA replication and maintain genomic RNA integrity. These new findings promote our understanding of the basics of LASV infection and inform antiviral and vaccine development.

Evolutional transition of HBV genome during the persistent infection determined by single-molecule real-time sequencing

BACKGROUND

Although HBV infection is a serious health issue worldwide, the landscape of HBV genome dynamics in the host has not yet been clarified. This study aimed to determine the continuous genome sequence of each HBV clone using a single-molecule real-time sequencing platform, and clarify the dynamics of structural abnormalities during persistent HBV infection without antiviral therapy.

PATIENTS AND METHODS

Twenty-five serum specimens were collected from 10 untreated HBV-infected patients. Continuous whole-genome sequencing of each clone was performed using a PacBio Sequel sequencer; the relationship between genomic variations and clinical information was analyzed. The diversity and phylogeny of the viral clones with structural variations were also analyzed.

RESULTS

The whole-genome sequences of 797,352 HBV clones were determined. The deletion was the most common structural abnormality and concentrated in the preS/S and C regions. Hepatitis B e antibody (anti-HBe)-negative samples or samples with high alanine aminotransferase levels have significantly diverse deletions than anti-HBe-positive samples or samples with low alanine aminotransferase levels. Phylogenetic analysis demonstrated that various defective and full-length clones evolve independently and form diverse viral populations.

CONCLUSIONS

Single-molecule real-time long-read sequencing revealed the dynamics of genomic quasispecies during the natural course of chronic HBV infections. Defective viral clones are prone to emerge under the condition of active hepatitis, and several types of defective variants can evolve independently of the viral clones with the full-length genome.

Stem cells as an ideal carrier for gene therapy: A new approach to the treatment of hepatitis C virus

BACKGROUND AND AIM

Various chemical drugs have been approved for the treatment of patients with hepatitis C, but most of these treatments are costly, and also have an inadequate response and many side effects. Also, there is no effective vaccine for hepatitis C due to its high genetic diversity. In recent decades, clinical trials have grown dramatically regarding the benefits of stem cell therapy as a modulator of immune system responses and anti-inflammatory drugs. The most promising point in stem cell therapy and similar therapies is that patients with chronic pain and severe injuries are offered drug-free treatment or surgery. In the present study, we examine the various dimensions of the use of stem cells with the approach of gene therapy carriers as a new treatment method in the treatment of Hepatitis C.

METHODS

Search terms were including gene carrier, stem cell therapy, gene therapy, liver disorders, hepatitis C virus. At first, 1000 article titles related to the mentioned keywords for different diseases were found. After removing duplicate titles and items that did not match the scope of the research, articles that met the criteria for entering the research and had usable information were selected. All abstracts of selected articles were studied by researchers. In the initial review, articles related to the title were identified and categorized based on the type of challenge.

CONCLUSION

Gene therapy, either directly and in vivo or indirectly and in vitro, requires carriers (vectors) to transfer the gene. These carriers are divided into two groups, viral and non-viral. In indirect gene therapy, living cells are isolated from a person's body and genetically modified. Stem cells have the properties to transfer the desired genes to the patient's body, including the ability to proliferate for a long time and differentiate into the tissue cells in which they are located.

Genetic diversity of hepatitis C virus and resistance associated substitutions to direct-acting antiviral treatment in Colombia

Hepatitis C virus (HCV) infection is one of the leading risk factors for end-stage liver disease development worldwide. This RNA virus displays high genetic diversity with 8 genotypes and 96 subgenotypes with heterogeneous geographical distribution around the world. In this study, we carried out an active case finding of individuals with a history of transfusion events before 1996 in three cities in Colombia. Then, the characterization of the HCV genotypes, subgenotypes, and resistance associate substitutions (RAS) was performed in samples positives for antibodies anti-HCV + from this study population. In addition, samples from PWID and patients with end-stage liver disease submitted to liver transplantation were included in the phylogenetic and RAS analysis. The 5'UTR, NS5A, and NS5B regions of the HCV genome were amplified in serum or liver explants samples. After the edition, assembly, and alignment of the sequences, genotyping through phylogenetic analysis was performed using IQTREE V2.0.5 based on the maximum likelihood approach. The identification of RAS was carried out by alignments based on the reference sequence (GenBank NC_004102). Two hundred sixty individuals with blood transfusion events before 1996 were recruited. The seroprevalence of antibodies anti-HCV was 2.69% in this population. The HCV genotypes 1, 2, and 4 and subgenotypes 1a, 1b, 2a, 4a and 4d were characterized in samples of the study populations. Three RAS (Q30R, C316N, and Y93H) were identified in samples obtained from 2 individuals who received blood transfusion before 1996 and without previous antiviral treatment and 6 samples obtained from patients with end-stage liver disease. Among the 20 samples analyzed, the HCV genotype 1, subgenotype 1b, was the most frequent (60%). We report the first characterization of HCV subgenotypes 4a and 4d and the first RAS identification in patients in Colombia.

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.

Natural Selection, Intracellular Bottlenecks of Virus Populations, and Viral Superinfection Exclusion

Natural selection acts on cellular organisms by ensuring the genes responsible for an advantageous phenotype consistently reap the phenotypic advantage. This is possible because reproductive cells of these organisms are almost always haploid, separating the beneficial gene from its rival allele at every generation. How natural selection acts on plus-strand RNA viruses is unclear because these viruses frequently load host cells with numerous genome copies and replicate thousands of progeny genomes in each cell. Recent studies suggest that these viruses encode the Bottleneck, Isolate, Amplify, Select (BIAS) mechanism that blocks all but a few viral genome copies from replication, thus creating the environment in which the bottleneck-escaping viral genome copies are isolated from each other, allowing natural selection to reward beneficial mutations and purge lethal errors. This BIAS mechanism also blocks the genomes of highly homologous superinfecting viruses, thus explaining cellular-level superinfection exclusion. Expected final online publication date for the Annual Review of Virology, Volume 9 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Mutational spectrum of hepatitis C virus in patients with chronic hepatitis C determined by single molecule real-time sequencing

The emergence of hepatitis C virus (HCV) with resistance-associated substitution (RAS), produced by mutations in the HCV genome, is a major problem in direct acting antivirals (DAA) treatment. This study aimed to clarify the mutational spectrum in HCV-RNA and the substitution pattern for the emergence of RASs in patients with chronic HCV infection. HCV-RNA from two HCV replicon cell lines and the serum HCV-RNA of four non-liver transplant and four post-liver transplant patients with unsuccessful DAA treatment were analyzed using high-accuracy single-molecule real-time long-read sequencing. Transition substitutions, especially A>G and U>C, occurred prominently under DAAs in both non-transplant and post-transplant patients, with a mutational bias identical to that occurring in HCV replicon cell lines during 10-year culturing. These mutational biases were reproduced in natural courses after DAA treatment. RASs emerged via both transition and transversion substitutions. NS3-D168 and NS5A-L31 RASs resulted from transversion mutations, while NS5A-Y93 RASs was caused by transition substitutions. The fidelity of the RNA-dependent RNA polymerase, HCV-NS5B, produces mutational bias in the HCV genome, characterized by dominant transition mutations, notably A>G and U>C substitutions. However, RASs are acquired by both transition and transversion substitutions, and the RASs-positive HCV clones are selected and proliferated under DAA treatment pressure.

Identifying Different Mutation Sites Leading to Resistance to the Direct-Acting Antiviral (DAA) Sofosbuvir in Hepatitis C Virus Patients from Egypt

The hepatitis C virus (HCV) is a major global health challenge and a leading cause of morbidity and mortality. Many direct-acting antivirals (DAAs) target essential macromolecules involved in the virus’ life cycle. Although such DAAs achieve great success in reducing the viral load in genotype 1 infections, other genotypes demonstrate different levels of response. This study focused on mutation sites associated with patients with genotype 4a infections that failed to respond to treatment with sofosbuvir. The genotyping of HCV samples from patients with virological failure, and responder patients, was conducted using Geno2Pheno webserver-based full NS5B sequences. We constructed 3D structural models for all the samples and used structural analysis to investigate the effect of amino acid substitution on the observed resistance to SOF-based treatment, and the docking of sofosbuvir into the active sites of the 10 models was performed. Finally, 10 molecular dynamic (MD) simulation experiments were conducted to compare the stability of the 3D models of the resistant samples against the stability of the 3D models of the responder samples. The results highlighted the presence of HCV subtype 4a in all ten samples; in addition, an amino acid (aa) substitution in the palm region may hinder HCV polymerase activity. In this study, we provide evidence that a mutation in the NS5B gene that induces resistance to sofosbuvir in patients with the S282T/C/R mutant virus is present in the Egyptian population. Overall, the docking and MD results support our findings and highlight the significant impact of the identified mutations on the resistance of HCV NS5B RNA-dependent RNA polymerase to direct-acting antivirals (DAAs).

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.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.

SARS-CoV-2 Quasispecies Provides an Advantage Mutation Pool for the Epidemic Variants

The dynamics of quasispecies afford RNA viruses a great fitness on cell tropism and host range. To study the quasispecies features and the intra-host evolution of SARS-CoV-2, we collected nine confirmed patients and sequenced the haplotypes of spike gene using a single-molecule real-time platform. Fourteen samples were extracted from sputum, nasopharyngeal swabs, or stool, which in total produced 283,655 high-quality circular consensus sequences. We observed a stable quasispecies structure that one master mutant (mean abundance ∼0.70), followed by numerous minor mutants (mean abundance ∼1.21 × 10⁻³). Under high selective pressure, minor mutants may obtain a fitness advantage and become the master ones. The later predominant substitution D614G existed in the minor mutants of more than one early patient. An epidemic variant had a possibility to be independently originated from multiple hosts. The mutant spectrums covered ∼85% amino acid variations of public genomes (GISAID; frequency ≥ 0.1) and likely provided an advantage mutation pool for the current/future epidemic variants. Notably, 32 of 35 collected antibody escape substitutions were preexistent in the early quasispecies. Virus populations in different tissues/organs revealed potentially independent replications. The quasispecies complexity of sputum samples was significantly lower than that of nasopharyngeal swabs (P = 0.02). Evolution analysis revealed that three continuous S2 domains (HR1, CH, and CD) had undergone a positive selection. Cell fusion-related domains may play a crucial role in adapting to the intrahost immune system. Our findings suggested that future epidemiologic investigations and clinical interventions should consider the quasispecies information that has missed by routine single consensus genome.

IMPORTANCE RNA virus population in a host does not consist of a consensus single haplotype but rather an ensemble of related sequences termed quasispecies. The dynamics of quasispecies afford SARS-CoV-2 a great ability on genetic fitness during intrahost evolution. The process is likely achieved by changing the genetic characteristics of key functional genes, such as the spike glycoprotein. Previous studies have applied the next-generation sequencing (NGS) technology to evaluate the quasispecies of SARS-CoV-2, and results indicated a low genetic diversity of the spike gene. However, the NGS platform cannot directly obtain the full haplotypes without assembling, and it is also difficult to predict the extremely low-frequency variations. Therefore, we introduced a single-molecule real-time technology to directly obtain the haplotypes of the RNA population and further study the quasispecies features and intrahost evolution of the spike gene.

LUCS: a high-resolution nucleic acid sequencing tool for accurate long-read analysis of individual DNA molecules

Nucleic acid sequence analyses are fundamental to all aspects of biological research, spanning aging, mitochondrial DNA (mtDNA) and cancer, as well as microbial and viral evolution. Over the past several years, significant improvements in DNA sequencing, including consensus sequence analysis, have proven invaluable for high-throughput studies. However, all current DNA sequencing platforms have limited utility for studies of complex mixtures or of individual long molecules, the latter of which is crucial to understanding evolution and consequences of single nucleotide variants and their combinations. Here we report a new technology termed LUCS (Long-molecule UMI-driven Consensus Sequencing), in which reads from third-generation sequencing are aggregated by unique molecular identifiers (UMIs) specific for each individual DNA molecule. This enables in-silico reconstruction of highly accurate consensus reads of each DNA molecule independent of other molecules in the sample. Additionally, use of two UMIs enables detection of artificial recombinants (chimeras). As proof of concept, we show that application of LUCS to assessment of mitochondrial genomes in complex mixtures from single cells was associated with an error rate of 1X10^-4 errors/nucleotide. Thus, LUCS represents a major step forward in DNA sequencing that offers high-throughput capacity and high-accuracy reads in studies of long DNA templates and nucleotide variants in heterogenous samples.