Mapping the genomic diversity of HCV subtypes 1a and 1b: Implications of structural and immunological constraints for vaccine and drug development

Identification of a novel replication-competent hepatitis C virus variant that confers the sofosbuvir resistance

Despite the excellent antiviral potency of direct-acting antivirals (DAAs) against hepatitis C virus (HCV), emergence of drug-resistant viral mutations remains a potential challenge. Sofobuvir (SOF), a nucleotide analog targeting HCV NS5B - RNA-dependent RNA polymerase (RdRp), constitutes a key component of many anti-HCV cocktail regimens and confers a high barrier for developing drug resistance. The serine to threonine mutation at the amino acid position 282 of NS5B (S282T) is the mostly documented SOF resistance-associated substitution (RAS), but severely hampers the virus fitness. In this study, we first developed new genotype 1b (GT1b) subgenomic replicon cells, denoted PR52D4 and PR52D9, directly from a GT1b clinical isolate. Next, we obtained SOF-resistant and replication-competent PR52D4 replicon by culturing the replicon cells in the presence of SOF. Sequencing analysis showed that the selected replicon harbored two mutations K74R and S282T in NS5B. Reverse genetics analysis showed that while PR52D4 consisting of either single mutation K74R or S282T could not replicate efficiently, the engineering of the both mutations led to a replication-competent and SOF-resistant PR52D4 replicon. Furthermore, we showed that the K74R mutation could also rescue the replication deficiency of the S282T mutation in Con1, another GT1b replicon as well as in JFH1, a GT2a replicon. Structural modeling analysis suggested that K74R might help maintain an active catalytic conformation of S282T by engaging with Y296. In conclusion, we identified the combination of two NS5B mutations S282T and K74R as a novel RAS that confers a substantial resistance to SOF while retains the HCV replication capacity.

Computational strategies to combat COVID-19: useful tools to accelerate SARS-CoV-2 and coronavirus research

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) is a novel virus of the family Coronaviridae. The virus causes the infectious disease COVID-19. The biology of coronaviruses has been studied for many years. However, bioinformatics tools designed explicitly for SARS-CoV-2 have only recently been developed as a rapid reaction to the need for fast detection, understanding and treatment of COVID-19. To control the ongoing COVID-19 pandemic, it is of utmost importance to get insight into the evolution and pathogenesis of the virus. In this review, we cover bioinformatics workflows and tools for the routine detection of SARS-CoV-2 infection, the reliable analysis of sequencing data, the tracking of the COVID-19 pandemic and evaluation of containment measures, the study of coronavirus evolution, the discovery of potential drug targets and development of therapeutic strategies. For each tool, we briefly describe its use case and how it advances research specifically for SARS-CoV-2. All tools are free to use and available online, either through web applications or public code repositories. Contact:evbc@unj-jena.de.

Cytomegalovirus and disseminated histoplasmosis-related hemophagocytic lymphohistiocytosis syndrome in an HIV-patient late presenter with IRIS: a case report

Background

Hemophagocytic lymphohistiocytosis syndrome (HLS) is an immune-mediated life-threatening disease considered as a medical emergency, with a potentially fatal multisystem inflammatory outcome. We present a patient that developed HLS and was able to be diagnosed efficiently with the help of an academic research institute of immunology.

Case presentation

A 21 years old male Mexican with human immunodeficiency virus (HIV), late presenter; who developed cytomegalovirus (CMV) infection and a disseminated histoplasmosis-related HLS, as part of an immune reconstitution inflammatory syndrome (IRIS). The patient required a long course of corticotherapy, intravenous immunoglobulin and massive transfusions (more than 10 units in 24 h, and a total of 83 units), besides amphotericin-B and ganciclovir treatment. An academic research institute of immunology aided in the accurate diagnosis of HLS with the implementation of tests not available within the hospital, thus improving the care provided to the patient. The patient recovered, was discharged, and continue to improve.

Conclusion

The objective of this report is to highlight the importance of having multidisciplinary support, including basic medical sciences groups providing specific tests that are sometimes very difficult to get, which provides a benefit to patients in the well-aimed diagnosis as part of applied translational medicine.

Prevalence of HCV resistance-associated substitutions among treatment-failure patients receiving direct-acting antiviral agents

Direct-acting antiviral (DAA) failure, which is mainly associated with the selection of resistance-associated substitutions (RASs), is not rare in HCV treatment. RAS data collected from published literature and RAS prevalence were integrated using meta-analysis. DAA-failure-associated RASs were identified by comparing their prevalence between DAA-failure and DAA-naïve patients. Prevalences of emerging RASs that occurred during treatment were also estimated. A total of 2932 DAA-naïve patients and 1466 DAA-failure patients were included. Significant differences in the prevalence of RASs were found in 76 scenarios that involved 34 RASs (11 in NS3, 18 in NS5A and 5 in NS5B), 4 genotypes (GTs) (GT1a, GT1b and GT3-4) and 14 DAAs (6 NS3 protease inhibitors [PIs], 6 NS5A inhibitors and 2 NS5B inhibitors). For NS3, the DAA-failure-associated RASs included V36L, Y56H, Q80K/R, R155K, A156T and D168A/E/L/T/V/Y. Substitutions at R155 and D168 were dominant for most NS3 PIs. For NS5A, DAA-failure-associated RASs included K24R, Q30R, L31M, and P32L in GT1a; R30Q/H, L31F/I/M/V, P58S, and Y93H in GT1b; A30K, L31M and Y93H in GT3; and M31V and Y93H in GT4. Y93H was the most prevalent RAS for NS5A inhibitors. DAA-failure-associated RASs were found at only five positions in NS5B. The majority of DAA-failure patients relapsed. A significant difference was detected for only four RAS sites between relapse patients and nonresponse/breakthrough patients. The RAS prevalence in DAA-failure patients varied among the HCV GTs and DAA regimens. The identified treatment-selected resistance patterns for broadly used DAA regimens will enable the selection of optimized retreatment options.

Factors Influencing the Prevalence of Resistance-Associated Substitutions in NS5A Protein in Treatment-Naive Patients with Chronic Hepatitis C

Direct-acting antivirals (DAAs) revolutionized treatment of hepatitis C virus (HCV) infection. Resistance-associated substitutions (RASs) present at the baseline impair response to DAA due to rapid selection of resistant HCV strains. NS5A is indispensable target of the current DAA treatment regimens. We evaluated prevalence of RASs in NS5A in DAA-naïve patients infected with HCV 1a (n = 19), 1b (n = 93), and 3a (n = 90) before systematic DAA application in the territory of the Russian Federation. Total proportion of strains carrying at least one RAS constituted 35.1% (71/202). In HCV 1a we detected only M28V (57.9%) attributed to a founder effect. Common RASs in HCV 1b were R30Q (7.5%), L31M (5.4%), P58S (4.4%), and Y93H (5.4%); in HCV 3a, A30S (31.0%), A30K (5.7%), S62L (8.9%), and Y93H (2.2%). Prevalence of RASs in NS5A of HCV 1b and 3a was similar to that worldwide, including countries practicing massive DAA application, i.e., it was not related to treatment. NS5A with and without RASs exhibited different co-variance networks, which could be attributed to the necessity to preserve viral fitness. Majority of RASs were localized in polymorphic regions subjected to immune pressure, with selected substitutions allowing immune escape. Altogether, this explains high prevalence of RAS in NS5A and low barrier for their appearance in DAA-inexperienced population.

Pooled analysis of HCV genotype 1 resistance-associated substitutions in NS5A, NS3 and NS5B pre-and post-treatment with 12 weeks of daclatasvir, asunaprevir and beclabuvir

BACKGROUND

Daclatasvir (DCV; non-structural [NS]5A inhibitor) plus asunaprevir (ASV; NS3 inhibitor) plus beclabuvir (BCV; non-nucleoside NS5B inhibitor) is an approved regimen for hepatitis C virus (HCV) genotype (GT)-1 treatment in Japan. A comprehensive analysis of pre-treatment and treatment-emergent HCV resistance to this regimen ± ribavirin (RBV) was performed.

METHODS

Data were pooled from five Phase 2/3 studies of DCV+ASV+BCV±RBV given for 12 weeks to GT-1a- or GT-1b-infected patients. The prevalence and impact of pre-treatment resistance-associated substitutions (RAS) in NS5A, NS3, and NS5B on sustained virological response (SVR) was assessed, as were emergent RAS and their post-treatment persistence.

RESULTS

Baseline NS5A RAS (GT-1a: M28T, Q30H/L/R/S, L31M, Y93C/H; GT-1b: L31I/M, Y93C/H) were present in 5% (26/561) of GT-1a and 16% (85/537) of GT-1b sequences. SVR12 for GT-1b without RBV was 100% (82/82) with RAS and >99% (427/428) without RAS. For GT-1a, SVR12 without RAS was 97% (85/88) with RBV and 92% (410/447) without RBV; SVR12 with RAS was 100% (2/2) with RBV and 54% (13/24) without RBV. Baseline NS3 (at R155 or D168) and NS5B (at P495) RAS were rare (≤1%). Treatment-emergent NS5A RAS (mostly Q30E/H/K/R±Y93H/N) in GT-1a persisted 60 weeks post-treatment, while NS3 RAS (mostly R155K) and NS5B-P495L/S were no longer detected after 48 or 24 weeks, respectively.

CONCLUSIONS

DCV+ASV+BCV±RBV was highly efficacious in HCV GT-1 infection, including HCV GT-1b with NS5A RAS. The fitness of treatment-emergent RAS post-treatment was NS5A > NS3 > NS5B; NS3 and NS5B RAS were generally replaced by wild-type sequence within 48 weeks.

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.

Evolutionary dynamics of Hepatitis C virus in a chronic HIV co-infected patient and its correlation with the immune status

The HCV evolutionary dynamics play a key role in the infection onset, maintenance of chronicity, pathogenicity, and drug resistance variants fixation, and are thought to be one of the main caveats in the development of an effective vaccine. Previous studies in HCV/HIV co-infected patients suggest that a decline in the immune status is related with increases in the HCV intra-host genetic diversity. However, these findings are based on single point sequence diversity measures or coalescence analyses in several virus-host interactions. In this work, we describe the molecular evolution of HCV-E2 region in a single HIV-co-infected patient with two clearly defined immune conditions. The phylogenetic analysis of the HCV-1a sequences from the studied patient showed that he was co-infected with three different viral lineages. These lineages were not evenly detected throughout time. The sequence diversity and coalescence analyses of these lineages suggested the action of different evolutionary patterns in different immune conditions: a slow rate, drift-like process in an immunocompromised condition (low levels of CD4+ T lymphocytes); and a fast rate, variant-switch process in an immunocompetent condition (high levels of CD4+ T lymphocytes).

Intrapatient viral diversity and treatment outcome in patients with genotype 3a hepatitis C virus infection on sofosbuvir-containing regimens

Treatment with the direct-acting antiviral agent (DAA) sofosbuvir (SOF), an NS5B inhibitor, and velpatasvir (VEL), an NS5A inhibitor, demonstrates viral cure rates of ≥95% in hepatitis C virus (HCV) genotypes (GT) 1-6. Here, we investigated intrapatient HCV diversity in NS5A and NS5B using Shannon entropy to examine the relationship between viral diversity and treatment outcome. At baseline, HCV diversity was lowest in patients infected with HCV GT3 as compared to the other GTs, and viral diversity was greater in NS5A than NS5B (P < .0001). Treatment outcome with SOF/VEL or the comparator regimen of SOF with ribavirin (RBV) was not correlated with baseline diversity. However, among persons treated with SOF/VEL, a decrease in diversity from baseline was observed at relapse in the majority virologic failures, consistent with a viral bottleneck event at relapse. In contrast, an increase in diversity was observed in 27% of SOF+RBV virologic failures. We investigated whether the increase in diversity was due to an increase in the transition rate, one mode of potential RBV-mediated mutagenesis; however, we found no evidence of this mechanism. Overall, we did not observe that viral diversity at baseline influenced treatment outcome, but the diversity changes observed at relapse can improve our understanding of RBV viral suppression in vivo.

Exploring resistance pathways for first-generation NS3/4A protease inhibitors boceprevir and telaprevir using Bayesian network learning

Resistance-associated variants (RAVs) have been shown to influence treatment response to direct-acting antivirals (DAAs) and first generation NS3/4A protease inhibitors (PIs) in particular. Interpretation of hepatitis C virus (HCV) genotypic drug resistance remains a challenge, especially in patients who previously failed DAA therapy and need to be retreated with a second DAA based regimen. Bayesian network (BN) learning on HCV sequence data from patients treated with DAAs could provide insight in resistance pathways against PIs for HCV subtypes 1a and 1b, in a similar way as applied before for HIV. The publicly available 'Rega-BN' tool chain was developed to study associative analyses for various pathogens. Our first analysis, comparing sequences from PI-naïve and PI-experienced patients, determined that NS3 substitutions R155K and V36M arise with PI-exposure in HCV1a infected patients, and were defined as major and minor resistance-associated variants respectively. NS3 variant 174H was newly identified as potentially related to PI resistance. In a second analysis, NS3 sequences from PI-naïve patients who cleared the virus during PI therapy and from PI-naïve patients who failed PI therapy were compared, showing that NS3 baseline variant 67S predisposes to treatment-failure and variant 72I to treatment success. This approach has the potential to better characterize the role of more RAVs, if sufficient therapy annotated sequence data becomes available in curated public databases. In addition, polymorphisms present in baseline sequences that predispose patients to therapy failure can be identified using this approach.

The epidemic dynamics of hepatitis C virus subtypes 4a and 4d in Saudi Arabia

The relatedness between viral variants sampled at different locations through time can provide information pertinent to public health that cannot readily be obtained through standard surveillance methods. Here, we use virus genetic data to identify the transmission dynamics that drive the hepatitis C virus subtypes 4a (HCV4a) and 4d (HCV4d) epidemics in Saudi Arabia. We use a comprehensive dataset of newly generated and publicly available sequence data to infer the HCV4a and HCV4d evolutionary histories in a Bayesian statistical framework. We also introduce a novel analytical method for an objective assessment of the migration intensity between locations. We find that international host mobility patterns dominate over within country spread in shaping the Saudi Arabia HCV4a epidemic, while this may be different for the HCV4d epidemic. This indicates that the subtypes 4a and 4d burden can be most effectively reduced by combining the prioritized screening and treatment of Egyptian immigrants with domestic prevention campaigns. Our results highlight that the joint investigation of evolutionary and epidemiological processes can provide valuable public health information, even in the absence of extensive metadata information.

Implications of hepatitis C virus subtype 1a migration patterns for virus genetic sequencing policies in Italy

Background

In-depth phylogeographic analysis can reveal migration patterns relevant for public health planning. Here, as a model, we focused on the provenance, in the current Italian HCV subtype 1a epidemic, of the NS3 resistance-associated variant (RAV) Q80K, known to interfere with the action of NS3/4A protease inhibitor simeprevir. HCV1a migration patterns were analysed using Bayesian phylodynamic tools, capitalising on newly generated and publicly available time and geo-referenced NS3 encoding virus genetic sequence data.

Results

Our results showed that both immigration and local circulation fuel the current Italian HCV1a epidemic. The United States and European continental lineages dominate import into Italy, with the latter taking the lead from the 1970s onwards. Since similar migration patterns were found for Q80K and other lineages, no clear differentiation of the risk for failing simeprevir can be made between patients based on their migration and travel history. Importantly, since HCV only occasionally recombines, these results are readily transferable to the genetic sequencing policy concerning NS5A RAVs.

Conclusions

The patient migration and travel history cannot be used to target only part of the HCV1a infected population for drug resistance testing before start of antiviral therapy. Consequently, it may be cost-effective to expand genotyping efforts to all HCV1a infected patients eligible for simeprevir-based therapies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0913-3) contains supplementary material, which is available to authorized users.

Comparative analysis of variation and selection in the HCV genome

Genotype 1 of the hepatitis C virus (HCV) is the most prevalent of the variants of this virus. Its two main subtypes, HCV-1a and HCV-1b, are associated to differences in epidemic features and risk groups, despite sharing similar features in most biological properties. We have analyzed the impact of positive selection on the evolution of these variants using complete genome coding regions, and compared the levels of genetic variability and the distribution of positively selected sites. We have also compared the distributions of positively selected and conserved sites considering different factors such as RNA secondary structure, the presence of different epitopes (antibody, CD4 and CD8), and secondary protein structure. <10% of the genome was found to be under positive selection, and purifying selection was the main evolutionary process acting in both subtypes. We found differences in the number of positively selected sites between subtypes in several genes (Core, HVR2 in E2, P7, helicase in NS3 and NS4a). Heterozygosity values in positively selected sites and the rate of non-synonymous substitutions were significantly higher in subtype HCV-1b. Logistic regression analyses revealed that similar selective forces act at the genome level in both subtypes: RNA secondary structure and CD4 T-cell epitopes are associated with conserved sites, while CD8 T-cell epitopes are associated with positive selection in both subtypes. These results indicate that similar selective constraints are acting along HCV-1a and HCV-1 b genomes, despite some differences in the distribution of positively selected sites at independent genes.