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

Synergism between remdesivir and ribavirin leads to SARS-CoV-2 extinction in cell culture

BACKGROUND AND PURPOSE

There is a need for effective anti-COVID-19 treatments, mainly for individuals at risk of severe disease such as the elderly and the immunosuppressed. Drug repositioning has proved effective in identifying drugs that can find a new application for the control of coronavirus disease, in particular COVID-19. The purpose of the present study was to find synergistic antiviral combinations for COVID-19 based on lethal mutagenesis.

EXPERIMENTAL APPROACH

The effect of combinations of remdesivir and ribavirin on the infectivity of SARS-CoV-2 in cell culture has been tested. Viral populations were monitored by ultra-deep sequencing, and the decrease of infectivity as a result of the treatment was measured.

KEY RESULTS

Remdesivir and ribavirin exerted a synergistic inhibitory activity against SARS-CoV-2, quantified both by CompuSyn (Chou-Talalay method) and Synergy Finder (ZIP-score model). In serial passage experiments, virus extinction was readily achieved with remdesivir-ribavirin combinations at concentrations well below their cytotoxic 50 value, but not with the drugs used individually. Deep sequencing of treated viral populations showed that remdesivir, ribavirin, and their combinations evoked significant increases of the number of viral mutations and haplotypes, as well as modification of diversity indices that characterize viral quasi-species.

CONCLUSION AND IMPLICATIONS

SARS-CoV-2 extinction can be achieved by synergistic combination treatments based on lethal mutagenesis. In addition, the results offer prospects of triple drug treatments for effective SARS-CoV-2 suppression.

Genotyping Hepatitis B virus by Next-Generation Sequencing: Detection of Mixed Infections and Analysis of Sequence Conservation

Our aim was to develop an accurate, highly sensitive method for HBV genotype determination and detection of genotype mixtures. We examined the preS and 5' end of the HBV X gene (5X) regions of the HBV genome using next-generation sequencing (NGS). The 1852 haplotypes obtained were subjected to genotyping via the Distance-Based discrimination method (DB Rule) using two sets of 95 reference sequences of genotypes A-H. In clinical samples from 125 patients, the main genotypes were A, D, F and H in Caucasian, B and C in Asian and A and E in Sub-Saharan patients. Genotype mixtures were identified in 28 (22.40%) cases, and potential intergenotypic recombination was observed in 29 (23.20%) cases. Furthermore, we evaluated sequence conservation among haplotypes classified into genotypes A, C, D, and E by computing the information content. The preS haplotypes exhibited limited shared conserved regions, whereas the 5X haplotypes revealed two groups of conserved regions across the genotypes assessed. In conclusion, we developed an NGS-based HBV genotyping method utilizing the DB Rule for genotype classification. We identified two regions conserved across different genotypes at 5X, offering promising targets for RNA interference-based antiviral therapies.

Molecular characterization of the nonstructural 5A (NS5A) region of hepatitis C virus in Thai blood donors

Direct acting antivirals (DAAs) have been developed for hepatitis C virus (HCV) therapy, and they are usually effective, however resistance to DAA regimens has also been reported to have a significant impact. Resistance associated substitutions (RASs) in the NS5A region are known to be correlated with failure of DAA therapy. HCV genotypes 3a and 1 are the most prevalent genotypes in Thailand. This study analyzed the type and frequency of RASs associated with DAA failure, focusing on the NS5A region. Serum samples of HCV genotype 3a, 1a, and 1b infection from Thai blood donors were selected. The NS5A region was amplified using reverse transcription-polymerase chain reaction (RT-PCR). A phylogenetic tree was constructed to identify the genotypes of HCV. Nucleotide sequencing and amino acid sequencing were conducted to determine the prevalence of RASs. Construction of the phylogenetic tree indicated that 29 samples were genotype 3a, 11 samples were genotype 1a, and 9 were genotype 1b. Both HCV genotypes 1a and 3a can be categorized into two subclades. Results showed that the NS5A substitutions A30V/K, A62T/V/I/M/P/S/L, and S98G were present in HCV genotype 3a. In HCV genotype 1a, only NS5A RASs H54Y was detected. NS5A amino acid substitutions Q54H and P58L were found in HCV genotype 1b. In conclusion, NS5A RASs at amino acid positions 30, 62, 54, 58, and 98 are present within HCV genotypes 3a and 1. While keeping in mind that additional information was not available on the anonymous blood donors tested in this study, these findings can contribute to understand the NS5A mutation. Further study with known patients under drug treatment is recommended.

Fitness-Dependent, Mild Mutagenic Activity of Sofosbuvir for Hepatitis C Virus

The concept of a mild mutagen was coined to describe a minor mutagenic activity exhibited by some nucleoside analogues that potentiated their efficacy as antiretroviral agents. In the present study, we report the mild mutagen activity of sofosbuvir (SOF) for hepatitis C virus (HCV). Serial passages of HCV in human hepatoma cells, in the presence of SOF at a concentration well below its cytotoxic concentration 50 (CC50) led to pre-extinction populations whose mutant spectra exhibited a significant increase of C→U transitions, relative to populations passaged in the absence of SOF. This was reflected in an increase in several diversity indices that were used to characterize viral quasispecies. The mild mutagenic activity of SOF was largely absent when it was tested with isogenic HCV populations that displayed high replicative fitness. Thus, SOF can act as a mild mutagen for HCV, depending on HCV fitness. Possible mechanisms by which the SOF mutagenic activity may contribute to its antiviral efficacy are discussed.

Atypical Mutational Spectrum of SARS-CoV-2 Replicating in the Presence of Ribavirin

We report that ribavirin exerts an inhibitory and mutagenic activity on SARS-CoV-2-infecting Vero cells, with a therapeutic index higher than 10. Deep sequencing analysis of the mutant spectrum of SARS-CoV-2 replicating in the absence or presence of ribavirin indicated an increase in the number of mutations, but not in deletions, and modification of diversity indices, expected from a mutagenic activity. Notably, the major mutation types enhanced by replication in the presence of ribavirin were A→G and U→C transitions, a pattern which is opposite to the dominance of G→A and C→U transitions previously described for most RNA viruses. Implications of the inhibitory activity of ribavirin, and the atypical mutational bias produced on SARS-CoV-2, for the search for synergistic anti-COVID-19 lethal mutagen combinations are discussed.

Quantifying In-Host Quasispecies Evolution

What takes decades, centuries or millennia to happen with a natural ecosystem, it takes only days, weeks or months with a replicating viral quasispecies in a host, especially when under treatment. Some methods to quantify the evolution of a quasispecies are introduced and discussed, along with simple simulated examples to help in the interpretation and understanding of the results. The proposed methods treat the molecules in a quasispecies as individuals of competing species in an ecosystem, where the haplotypes are the competing species, and the ecosystem is the quasispecies in a host, and the evolution of the system is quantified by monitoring changes in haplotype frequencies. The correlation between the proposed indices is also discussed, and the R code used to generate the simulations, the data and the plots is provided. The virtues of the proposed indices are finally shown on a clinical case.

Rapid genotyping of targeted viral samples using Illumina short-read sequencing data

The most important information about microorganisms might be their accurate genome sequence. Using current Next Generation Sequencing methods, sequencing data can be generated at an unprecedented pace. However, we still lack tools for the automated and accurate reference-based genotyping of viral sequencing reads. This paper presents our pipeline designed to reconstruct the dominant consensus genome of viral samples and analyze their within-host variability. We benchmarked our approach on numerous datasets and showed that the consensus genome of samples could be obtained reliably without further manual data curation. Our pipeline can be a valuable tool for fast identifying viral samples. The pipeline is publicly available on the project’s GitHub page (https://github.com/laczkol/QVG).

Ultra Deep Sequencing of Circulating Cell-Free DNA as a Potential Tool for Hepatocellular Carcinoma Management

Simple Summary

In this unicentric prospective study, we analyzed the most prevalent mutations in HCC (TERT promoter, TP53, CTNNB1, AXIN1 and ARID1A) in plasma cfDNA by next-generation sequencing, aiming to elucidate their value as prognostic noninvasive biomarkers. Total cfDNA (cut-off value 2 ng/µL), number of mutated genes and number of detectable mutations on cfDNA were significantly related to mortality. Number of mutated genes and number of detected mutations in cfDNA and the ratio between number of mutations and total amount of cfDNA were also significantly associated with recurrence. Detection of more than four mutations in cfDNA correlated with a higher risk of death. Dynamic changes in cfDNA mutations were detected prior to radiological detection of HCC recurrence. We believe that these results support the proof of principle and launching of validation studies to confirm that total cfDNA and detection of prevalent HCC mutations could have prognostic implications in early-stage HCC patients.

Abstract

Background: Cell-free DNA (cfDNA) concentrations have been described to be inversely correlated with prognosis in cancer. Mutations in HCC-associated driver genes in cfDNA have been reported, but their relation with patient’s outcome has not been described. Our aim was to elucidate whether mutations found in cfDNA could be representative from those present in HCC tissue, providing the rationale to use the cfDNA to monitor HCC. Methods: Tumoral tissue, paired nontumor adjacent tissue and blood samples were collected from 30 HCC patients undergoing curative therapies. Deep sequencing targeting HCC driver genes was performed. Results: Patients with more than 2 ng/µL of cfDNA at diagnosis had higher mortality (mean OS 24.6 vs. 31.87 months, p = 0.01) (AUC = 0.782). Subjects who died during follow-up, had a significantly higher number of mutated genes (p = 0.015) and number of mutations (p = 0.015) on cfDNA. Number of mutated genes (p = 0.001), detected mutations (p = 0.001) in cfDNA and ratio (number of mutations/cfDNA) (p = 0.003) were significantly associated with recurrence. However, patients with a ratio (number of mutations/cfDNA) above 6 (long-rank p = 0.0003) presented a higher risk of recurrence than those with a ratio under 6. Detection of more than four mutations in cfDNA correlated with higher risk of death (long-rank p = 0.042). Conclusions: In summary, cfDNA and detection of prevalent HCC mutations could have prognostic implications in early-stage HCC patients

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Discovery of Novel 3-Hydroxyquinazoline-2,4(1H,3H)-Dione Derivatives: A Series of Metal Ion Chelators with Potent Anti-HCV Activities

Millions of people worldwide suffer from acute or chronic liver inflammation caused by the hepatitis C virus (HCV). Metal ion chelators have achieved widespread success in the development of antiviral drugs. Some inhibitors with metal ion chelating structures have been proven to have good inhibitory activities on non-structural protein 5B (NS5B) polymerase. However, most of the reported metal ion chelators showed poor anti-HCV potency at the cellular level. Hence, we designed and synthesized a series of 3-hydroxyquinazoline-2,4(1H,3H)-dione derivatives with novel metal ion chelating structures. Typical compounds such as 21h, 21k, and 21t showed better anti-HCV activities than ribavirin with EC50 values less than 10 μM. 21t is currently known as one of the metal ion chelators with the best anti-HCV potency (EC50 = 2.0 μM) at the cellular level and has a better therapeutic index (TI > 25) as compared to ribavirin and the reported compound 6. In the thermal shift assay, the representative compounds 21e and 21k increased the melting temperature (Tm) of NS5B protein solution by 1.6 °C and 2.1 °C, respectively, at the test concentration, indicating that these compounds may exert an anti-HCV effect by targeting NS5B. This speculation was also supported by our molecular docking studies and ultraviolet-visible (UV-Vis) spectrophotometry assay, in which the possibility of binding of 3-hydroxyquinazoline-2,4(1H,3H)-diones with Mg2+ in the NS5B catalytic center was observed.

Vaccine breakthrough infections with SARS-CoV-2 Alpha mirror mutations in Delta Plus, Iota, and Omicron

Replication of SARS-CoV-2 in the human population is defined by distributions of mutants that are present at different frequencies within the infected host and can be detected by ultra-deep sequencing techniques. In this study, we examined the SARS-CoV-2 mutant spectra of amplicons from the spike-coding (S-coding) region of 5 nasopharyngeal isolates derived from patients with vaccine breakthrough. Interestingly, all patients became infected with the Alpha variant, but amino acid substitutions that correspond to the Delta Plus, Iota, and Omicron variants were present in the mutant spectra of the resident virus. Deep sequencing analysis of SARS-CoV-2 from patients with vaccine breakthrough revealed a rich reservoir of mutant types and may also identify tolerated substitutions that can be represented in epidemiologically dominant variants.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

Population Disequilibrium as Promoter of Adaptive Explorations in Hepatitis C Virus

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium—revealed by the changing composition of the mutant spectrum—may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed.

Sophisticated viral quasispecies with a genotype-related pattern of mutations in the hepatitis B X gene of HBeAg-ve chronically infected patients

Patients with HBeAg-negative chronic infection (CI) have not been extensively studied because of low viremia. The HBx protein, encoded by HBX, has a key role in viral replication. Here, we analyzed the viral quasispecies at the 5' end of HBX in CI patients and compared it with that of patients in other clinical stages. Fifty-eight HBeAg-negative patients were included: 16 CI, 19 chronic hepatitis B, 16 hepatocellular carcinoma and 6 liver cirrhosis. Quasispecies complexity and conservation were determined in the region between nucleotides 1255 and 1611. Amino acid changes detected were tested in vitro. CI patients showed higher complexity in terms of mutation frequency and nucleotide diversity and higher quasispecies conservation (p < 0.05). A genotype D-specific pattern of mutations (A12S/P33S/P46S/T36D-G) was identified in CI (median frequency, 81.7%), which determined a reduction in HBV DNA release of up to 1.5 log in vitro. CI patients showed a more complex and conserved viral quasispecies than the other groups. The genotype-specific pattern of mutations could partially explain the low viremia observed in these patients.

Amino Acid Substitutions Associated with Treatment Failure for Hepatitis C Virus Infection

Despite the high virological response rates achieved with current directly acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. The identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultradeep sequencing (UDS) methods, for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV, which are not bona fide resistance-associated substitutions (RAS), from 220 patients who failed therapy. They were present frequently in basal and posttreatment virus of patients who failed different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. They also have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value for treatment response, are discussed.

Dissimilar Conservation Pattern in Hepatitis C Virus Mutant Spectra, Consensus Sequences, and Data Banks

The influence of quasispecies dynamics on long-term virus diversification in nature is a largely unexplored question. Specifically, whether intra-host nucleotide and amino acid variation in quasispecies fit the variation observed in consensus sequences or data bank alignments is unknown. Genome conservation and dynamics simulations are used for the computational design of universal vaccines, therapeutic antibodies and pan-genomic antiviral agents. The expectation is that selection of escape mutants will be limited when mutations at conserved residues are required. This strategy assumes long-term (epidemiologically relevant) conservation but, critically, does not consider short-term (quasispecies-dictated) residue conservation. We calculated mutant frequencies of individual loci from mutant spectra of hepatitis C virus (HCV) populations passaged in cell culture and from infected patients. Nucleotide or amino acid conservation in consensus sequences of the same populations, or in the Los Alamos HCV data bank did not match residue conservation in mutant spectra. The results relativize the concept of sequence conservation in viral genetics and suggest that residue invariance in data banks is an insufficient basis for the design of universal viral ligands for clinical purposes. Our calculations suggest relaxed mutational restrictions during quasispecies dynamics, which may contribute to higher calculated short-term than long-term viral evolutionary rates.

Development and validation of a simple and rapid method for hepatitis C virus genotyping based on one-step RT-qPCR

Hepatitis C virus (HCV) infections caused by different subtypes require different treatments; therefore, rapid and cost-effective genotyping methods for the diagnosis of HCV are greatly needed. In the present study, a new method to diagnose HCV subtypes that depends on a one-step quantitative reverse transcription PCR (RT-qPCR) and TaqMan fluorescence probe technique is described. Five pairs of primers and five probes were designed, which were able to detect five genotypes in three reaction tubes. One reaction was used to detect the 1b subtype, one was used to detect the 2a and 6a subtypes, and the other was used to detect the 3a and 3b subtypes. Rigorous performance validation was implemented for five aspects: Precision, sensitivity, accuracy, specificity and anti-interference. The HCV subtype that infected 289 patients was evaluated in the present study via RT-qPCR and verified by sequencing. The results revealed that the 1b subtype accounted for 45% of infections, the 2a subtype accounted for 9% of infections, the 3a subtype accounted for 13% of infections, the 3b subtype accounted for 18% of infections, and the 6a subtype accounted for 15% of infections. The analytical sensitivity for the detection of each of the five HCV subtypes was 1,000 IU/ml. The new method performed well in the performance validation mentioned above, indicating its effectiveness as a HCV genotyping method. RT-qPCR has mitigated some of the former challenges of existing HCV genotyping methods, including the time commitment, expense, and inaccuracy of such methods. The performance validation of this new method showed that RT-qPCR is reliable enough to be widely applied in China for HCV genotyping.

Discovery of novel Hepatitis C virus inhibitor targeting multiple allosteric sites of NS5B polymerase

HCV is a viral infection posing a severe global threat when left untreated progress to end-stage liver disease, including cirrhosis and HCC. The NS5B polymerase of HCV is the most potent target that harbors four allosteric binding sites that could interfere with the HCV infection. We present the discovery of a novel synthetic compound that harbors the potential of NS5B polymerase inhibition. All eight compounds belonging to the benzothiazine family of heterocycles displayed no cellular cytotoxicity in HepG2 cells at nontoxic dose concentration (200 μM). Subsequently, among eight compounds of the series, merely compound 5b exhibited significant inhibition of the expression of the HCV NS5B gene as compared to DMSO control in semi-quantitative PCR. Based on our western blot result, 5b at the range of 50, 100 and 200 μM induced 20, 40, and 70% inhibition of NS5B protein respectively. To estimate the binding potential, 5b was docked at respective allosteric sites followed by molecular dynamics (MD) simulations for a period of 20 ns. In addition, binding free energy calculation by MM-GB/PBSA method revealed a conserved interaction profile of residues lining the allosteric sites in agreement with the reported NS5B co-crystallized inhibitors. The presented results provide important information about a novel compound 5b which may facilitate the the discovery of novel inhibitors that tends to target multiple sites on NS5B polymerase.

A new implication of quasispecies dynamics: Broad virus diversification in absence of external perturbations

RNA genetic elements include many important animal and plant pathogens. They share high mutability, a trait that has multiple implications for the interactions with their host organisms. Here we review evidence of a new adaptive feature of RNA viruses that we term "broadly diversifying selection". It constitutes a new type of positive selection without participation of any external selective agent, and which is built upon a progressive increase of the number of different genomes that dominate the population. The evidence was provided by analyses of mutant spectrum composition of two important viral pathogens, foot-and-mouth disease virus (FMDV) and hepatitis C virus (HCV) after prolonged replication in their respective cell culture environment. Despite being fueled by mutations that arise randomly and in absence of an external guiding selective force, this type of selection prepares the viral population for a response to selective forces still to occur. Since current evidence suggests that broadly diversifying selection is favored by elevated mutation rates and population sizes, it may constitute a more general behavior, relevant also to the adaptive dynamics of microbial populations and cancer cells.

Broad and Dynamic Diversification of Infectious Hepatitis C Virus in a Cell Culture Environment

Previous studies documented that long-term hepatitis C virus (HCV) replication in human hepatoma Huh-7.5 cells resulted in viral fitness gain, expansion of the mutant spectrum, and several phenotypic alterations. In the present work, we show that mutational waves (changes in frequency of individual mutations) occurred continuously and became more prominent as the virus gained fitness. They were accompanied by an increasing proportion of heterogeneous genomic sites that affected 1 position in the initial HCV population and 19 and 69 positions at passages 100 and 200, respectively. Analysis of biological clones of HCV showed that these dynamic events affected infectious genomes, since part of the fluctuating mutations became incorporated into viable genomes. While 17 mutations were scored in 3 biological clones isolated from the initial population, the number reached 72 in 3 biological clones from the population at passage 200. Biological clones differed in their responses to antiviral inhibitors, indicating a phenotypic impact of viral dynamics. Thus, HCV adaptation to a specific constant environment (cell culture without external influences) broadens the mutant repertoire and does not focus the population toward a limited number of dominant genomes. A retrospective examination of mutant spectra of foot-and-mouth disease virus passaged in cell cultures suggests a parallel behavior here described for HCV. We propose that virus diversification in a constant environment has its basis in the availability of multiple alternative mutational pathways for fitness gain. This mechanism of broad diversification should also apply to other replicative systems characterized by high mutation rates and large population sizes.IMPORTANCE The study shows that extensive replication of an RNA virus in a constant biological environment does not limit exploration of sequence space and adaptive options. There was no convergence toward a restricted set of adapted genomes. Mutational waves and mutant spectrum broadening affected infectious genomes. Therefore, profound modifications of mutant spectrum composition and consensus sequence diversification are not exclusively dependent on environmental alterations or the intervention of population bottlenecks.

Deep-sequencing reveals broad subtype-specific HCV resistance mutations associated with treatment failure

A percentage of hepatitis C virus (HCV)-infected patients fail direct acting antiviral (DAA)-based treatment regimens, often because of drug resistance-associated substitutions (RAS). The aim of this study was to characterize the resistance profile of a large cohort of patients failing DAA-based treatments, and investigate the relationship between HCV subtype and failure, as an aid to optimizing management of these patients. A new, standardized HCV-RAS testing protocol based on deep sequencing was designed and applied to 220 previously subtyped samples from patients failing DAA treatment, collected in 39 Spanish hospitals. The majority had received DAA-based interferon (IFN) α-free regimens; 79% had failed sofosbuvir-containing therapy. Genomic regions encoding the nonstructural protein (NS) 3, NS5A, and NS5B (DAA target regions) were analyzed using subtype-specific primers. Viral subtype distribution was as follows: genotype (G) 1, 62.7%; G3a, 21.4%; G4d, 12.3%; G2, 1.8%; and mixed infections 1.8%. Overall, 88.6% of patients carried at least 1 RAS, and 19% carried RAS at frequencies below 20% in the mutant spectrum. There were no differences in RAS selection between treatments with and without ribavirin. Regardless of the treatment received, each HCV subtype showed specific types of RAS. Of note, no RAS were detected in the target proteins of 18.6% of patients failing treatment, and 30.4% of patients had RAS in proteins that were not targets of the inhibitors they received. HCV patients failing DAA therapy showed a high diversity of RAS. Ribavirin use did not influence the type or number of RAS at failure. The subtype-specific pattern of RAS emergence underscores the importance of accurate HCV subtyping. The frequency of "extra-target" RAS suggests the need for RAS screening in all three DAA target regions.

Deep-sequencing study of HCV G4a resistance-associated substitutions in Egyptian patients failing DAA treatment

PURPOSE

 To study resistance-associated substitutions using next-generation sequencing in Egyptian hepatitis C virus-infected patients failing direct-acting antiviral treatment.

METHODS

 The current study describes three cases of treatment failure in patients referred to Zagazig Viral Hepatitis Treatment Center (ZVHTC), Sharkia Governorate, Egypt. RAS were identified and characterized using deep sequencing. The first patient had breakthrough while receiving a daclatasvir (DCV)+sofosbuvir (SOF) regimen, patient 2 relapsed after treatment with DCV+SOF+ribavirin (RBV), and patient 3 relapsed after DCV+SOF therapy. A serum sample was collected from each patient at failure and sent to Vall d'Hebron Research Institute at Hospital Universitari Vall d'Hebron in Barcelona (Spain) for deep-sequencing study to identify and characterize the RAS present in the samples.

RESULTS

 The following were identified: L28M, L30S and L28M+L30S in patient 1, L30R in patient 2, and R155C, D168E, L28M, L30H, L30S, L28M+L30H, and L28M+L30S in patient 3.

CONCLUSION

To the best of our knowledge, this is the first report from Egypt of patients failing DAA-based therapy, describing the associated RAS. This information will be of help to understand the natural history of HCV in Egyptian patients and guide the proper choice of retreatment protocols.

Hepatitis C virus early kinetics and resistance-associated substitution dynamics during antiviral therapy with direct-acting antivirals

The emergence of resistance-associated substitutions (RASs) can compromise the high efficacy of direct-acting antivirals (DAAs). Little is known about RASs selection at very early time points during DAA treatment. Therefore, we analyzed the potential emergence of RASs immediately after therapy initiation. Samples of 71 patients treated with different DAAs were collected at baseline, during therapy (hours 4 and 8; days 1-7; weeks 2-4) or until target not detected. HCV-RNA levels were determined by qPCR, and RASs were detected by deep sequencing. Sixty-three (89%) patients achieved a sustained virological response (SVR), 7 (10%) relapsed, and 1 (1%) experienced a breakthrough. Almost all non-SVR (7/8, 88%) showed RASs either at baseline or relapse. High-frequency RASs detected at baseline (Y93H and L159F+C316N) remained detectable at early time points during therapy and reappeared as most prevalent substitutions at relapse. Conversely, emergent RASs at relapse (Q80R, D168E/V, R155K and L31V) were not observed during the first hours-days, before HCV-RNA became undetectable. HCV-RNA decay and genetic evolution of the quasispecies followed a similar pattern during the first hours of therapy in SVR and non-SVR patients. In conclusion, the absence of early RASs selection and the similar dynamics of HCV kinetics and quasispecies in SVR and non-SVR patients after therapy initiation suggest that RASs selection may occur at later stages in the remaining reservoir, where viral populations persist hidden at very low replication levels. Nevertheless, we cannot completely exclude very early selection, when RASs are present below the sensitivity limit of deep sequencing.

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.

The EIMB Hydrogel Microarray Technology: Thirty Years Later

 Biological microarrays (biochips) are analytical tools that can be used to implement complex integrative genomic and proteomic approaches to the solution of problems of personalized medicine (e.g., patient examination in order to reveal the disease long before the manifestation of clinical symptoms, assess the severity of pathological or infectious processes, and choose a rational treatment). The efficiency of biochips is predicated on their ability to perform multiple parallel specific reactions and to allow one to study the interactions of biopolymer molecules, such as DNA, proteins, glycans, etc. One of the pioneers of microarray technology was the Engelhardt Institute of Molecular Biology of the Russian Academy of Sciences (EIMB), with its suggestion to immobilize molecular probes in the three-dimensional structure of a hydrophilic gel. Since the first experiments on sequencing by hybridization on oligonucleotide microarrays conducted some 30 years ago, the hydrogel microarrays designed at the EIMB have come a long and successful way from basic research to clinical laboratory diagnostics. This review discusses the key aspects of hydrogel microarray technology and a number of state-ofthe-art approaches for a multiplex analysis of DNA and the protein biomarkers of socially significant diseases, including the molecular genetic, immunological, and epidemiological aspects of pathogenesis.