Quasispecies and its impact on viral hepatitis

The dynamic variation position and predominant quasispecies of hepatitis B virus: Novel predictors of early hepatocarcinoma

To find the predictors of early HCC based on the dynamic changes of HBV quasispecies, this study utilizing the second-generation sequencing (NGS) and high-order multiplex droplet digital PCR (ddPCR) technology to examine the HBV quasispecies in serum of total 247 subjects recruited from high-incidence area of HCC. In the discovery stage, 15 non-synonymous Single Nucleotide Polymorphisms (SNPs) with higher variant proportion in HCC case group were founded (all P<0.05). Furthermore, the variant proportions in some of these SNPs were observed changing regularly within 5 years before the onset of HCC, and 5 of them located in HBX, 2 in HBS and 2 in HBC. The HBV predominant quasispecies and their consensus sequences were identified by genetic evolution analysis, in which the high HBS and HBC quasispecies heterogeneity were found associated with the forming of multifarious quasispecies clones, and the HBX gene had the highest proportion of predominant quasispecies (46.7 % in HBX vs 12.7 % and 13.8 % in HBS and HBC respectively) with the key variations (G1512A, A1630G, T1753C/G/A, A1762T and G1764A) determined. In the validation stage, we confirmed that the combined double mutations of G1512A+A1630G, A1762T+G1764A, and the combined triple mutations of T1753C/G/A + A1762T+G1764A, all expressed higher in early HCC cases when comparing with control group (all P<0.05). We also demonstrated the advantages of ddPCR using in multi-variations detection in large-sample for early HCC surveillance and screening. So we think that the dynamic of key HBV variation positions and their different combinations determined by quasispecies anlysis in this study can act as the novel predictors of early hepatocarcinoma and suitable to popularize and apply in HCC screening.

Low host immune pressure may be associated with the development of hepatocellular carcinoma: a longitudinal analysis of complete genomes of the HBV 1762T, 1764A mutant

Background

It has been reported that hepatitis B virus (HBV) double mutations (A1762T, G1764A) are an aetiological factor of hepatocellular carcinoma (HCC). However, it is unclear who is prone to develop HCC, among those infected with the mutant. Exploring HBV quasispecies, which are strongly influenced by host immune pressure, may provide more information about the association of viral factors and HCC.

Materials and methods

Nine HCC cases and 10 controls were selected from the Long An cohort. Serum samples were collected in 2004 and 2019 from subjects with HBV double mutations and the complete genome of HBV was amplified and sequenced using next-generation sequencing (NGS).

Results

The Shannon entropy values increased from 2004 to 2019 in most cases and controls. There was no significant difference in mean intrahost quasispecies genetic distances between cases and controls. The change in the values of mean intrahost quasispecies genetic distances of the controls between 2004 and 2019 was significantly higher than that of the cases (P<0.05). The viral loads did not differ significantly between cases and controls in 2004(p=0.086) but differed at diagnosed in 2019 (p=0.009). Three mutations occurring with increasing frequency from 2004 to 2019 were identified in the HCC cases, including nt446 C→G, nt514 A→C and nt2857T→C. Their frequency differed significantly between the cases and controls (P<0.05).

Conclusions

The change in the values of mean intrahost quasispecies genetic distances in HCC was smaller, suggesting that HBV in HCC cases may be subject to low host immune pressure. Increasing viral loads during long-term infection are associated with the development of HCC. The novel mutations may increase the risk for HCC.

Compartmentalisation of Hepatitis B virus X gene evolution in hepatocellular carcinoma microenvironment and the genotype-phenotype correlation of tumorigenicity in HBV-related patients with hepatocellular carcinoma

Hepatitis B virus (HBV) exists as quasispecies (QS). However, the evolutionary characteristics of haplotypes of HBV X gene in the hepatocellular carcinoma (HCC) microenvironment remain unclear. Mutations across X gene are essential for the tumorigenicity of HBV X protein (HBx). However, the functional phenotypes of many mutant HBx remain unknown. This study aims to compare the characteristics of X gene evolution between tumour and non-tumour tissues in HCC patients and investigate the tumorigenic phenotype of HBx harbouring mutation T81P/S101P/L123S. This study included 24 HCC patients. Molecular cloning of X gene was performed to analyse characteristics of haplotypes in liver tissues. HCC cell lines stably expressing wild-type or mutant HBx and subcutaneous tumour xenograft mouse model were used to assess HBx-T81P/S101P/L123S tumorigenicity. The mean heterogeneity of HBV QS across X gene in tumour tissues was lower than that in non-tumour tissues. A location bias was observed in X gene clones with genotype C or D in tumour tissues compared to those with genotype B. Mutations in genotype-C or - D clones were mainly clustered in the dimerization region and aa110-aa140 within the transactivation region. A novel mutation combination at residues 81, 101 and 123 was identified in tumour tissues. Further, HBx-T81P/S101P/L123S promotes cell proliferation and increases genomic instability, which was mediated by MYC. This study elucidates the compartmentalized evolution patterns of HBV X gene between intra tumour and non-tumour tissues in HCC patients and provides a new mechanism underlying HBV-driven hepatocarcinogenesis, suggesting a potential viral marker for monitoring HCC.

Targeting non-structural proteins of Hepatitis C virus for predicting repurposed drugs using QSAR and machine learning approaches

Hepatitis C virus (HCV) infection causes viral hepatitis leading to hepatocellular carcinoma. Despite the clinical use of direct-acting antivirals (DAAs) still there is treatment failure in 5–10% cases. Therefore, it is crucial to develop new antivirals against HCV. In this endeavor, we developed the “Anti-HCV” platform using machine learning and quantitative structure–activity relationship (QSAR) approaches to predict repurposed drugs targeting HCV non-structural (NS) proteins. We retrieved experimentally validated small molecules from the ChEMBL database with bioactivity (IC₅₀/EC₅₀) against HCV NS3 (454), NS3/4A (495), NS5A (494) and NS5B (1671) proteins. These unique compounds were divided into training/testing and independent validation datasets. Relevant molecular descriptors and fingerprints were selected using a recursive feature elimination algorithm. Different machine learning techniques viz. support vector machine, k-nearest neighbour, artificial neural network, and random forest were used to develop the predictive models. We achieved Pearson’s correlation coefficients from 0.80 to 0.92 during 10-fold cross validation and similar performance on independent datasets using the best developed models. The robustness and reliability of developed predictive models were also supported by applicability domain, chemical diversity and decoy datasets analyses. The “Anti-HCV” predictive models were used to identify potential repurposing drugs. Representative candidates were further validated by molecular docking which displayed high binding affinities. Hence, this study identified promising repurposed drugs viz. naftifine, butalbital (NS3), vinorelbine, epicriptine (NS3/4A), pipecuronium, trimethaphan (NS5A), olodaterol and vemurafenib (NS5B) etc. targeting HCV NS proteins. These potential repurposed drugs may prove useful in antiviral drug development against HCV.

Sparse logistic regression revealed the associations between HBV PreS quasispecies and hepatocellular carcinoma

BACKGROUND

Chronic infection with hepatitis B virus (HBV) has been proved highly associated with the development of hepatocellular carcinoma (HCC).

AIMS

The purpose of the study is to investigate the association between HBV preS region quasispecies and HCC development, as well as to develop HCC diagnosis model using HBV preS region quasispecies.

METHODS

A total of 104 chronic hepatitis B (CHB) patients and 117 HBV-related HCC patients were enrolled. HBV preS region was sequenced using next generation sequencing (NGS) and the nucleotide entropy was calculated for quasispecies evaluation. Sparse logistic regression (SLR) was used to predict HCC development and prediction performances were evaluated using receiver operating characteristic curves.

RESULTS

Entropy of HBV preS1, preS2 regions and several nucleotide points showed significant divergence between CHB and HCC patients. Using SLR, the classification of HCC/CHB groups achieved a mean area under the receiver operating characteristic curve (AUC) of 0.883 in the training data and 0.795 in the test data. The prediction model was also validated by a completely independent dataset from Hong Kong. The 10 selected nucleotide positions showed significantly different entropy between CHB and HCC patients. The HBV quasispecies also classified three clinical parameters, including HBeAg, HBVDNA, and Alkaline phosphatase (ALP) with the AUC value greater than 0.6 in the test data.

CONCLUSIONS

Using NGS and SLR, the association between HBV preS region nucleotide entropy and HCC development was validated in our study and this could promote the understanding of HCC progression mechanism.

CRISPR/Cas and Hepatitis B Therapy: Technological Advances and Practical Barriers

After almost a decade of using CRISPR/Cas9 systems to edit target genes, CRISPR/Cas9 and related technologies are rapidly moving to clinical trials. Hepatitis B virus (HBV), which causes severe liver disease, cannot be cleared by modern antivirals, but represents an ideal target for CRISPR/Cas9 systems. Early studies demonstrated very high antiviral potency of CRISPR/Cas9 and supported its use for developing a cure against chronic HBV infection. This review discusses the key issues that must be solved to make CRISPR/Cas9 an anti-HBV therapy.

Viral quasispecies quantitative analysis: a novel approach for appraising the immune tolerant phase of chronic hepatitis B virus infection

Few non-invasive models were established for precisely identifying the immune tolerant (IT) phase from chronic hepatitis B (CHB). This study aimed to develop a novel approach that combined next-generation sequencing (NGS) and machine learning algorithms using our recently published viral quasispecies (QS) analysis package. 290 HBeAg positive patients from whom liver biopsies were taken were enrolled and divided into a training group (n = 148) and a validation group (n = 142). HBV DNA was extracted and QS sequences were obtained by NGS. Hierarchical clustering analysis (HCA) and principal component analysis (PCA) based on viral operational taxonomic units (OTUs) were performed to explore the correlations among QS and clinical phenotypes. Three machine learning algorithms, including K-nearest neighbour, support vector machine, and random forest algorithm, were used to construct diagnostic models for IT phase classification. Based on histopathology, 90 IT patients and 200 CHB patients were diagnosed. HBsAg titres for IT patients were higher than those of CHB patients (p < 0.001). HCA and PCA analysis grouped IT and CHB patients into two distinct clusters. The relative abundance of viral OTUs differed mainly within the BCP/precore/core region and was significantly correlated with liver inflammation and fibrosis. For the IT phase classification, all machine-learning models showed higher AUC values compared to models based on HBsAg, APRI, and FIB-4. The relative abundance of viral OTUs reflects the severity of liver inflammation and fibrosis. The novel QS quantitative analysis approach could be used to diagnose IT patients more precisely and reduce the need for liver biopsy.

Social Networking of Quasi-Species Consortia drive Virolution via Persistence

The emergence of cooperative quasi-species consortia (QS-C) thinking from the more accepted quasispecies equations of Manfred Eigen, provides a conceptual foundation from which concerted action of RNA agents can now be understood. As group membership becomes a basic criteria for the emergence of living systems, we also start to understand why the history and context of social RNA networks become crucial for survival and function. History and context of social RNA networks also lead to the emergence of a natural genetic code. Indeed, this QS-C thinking can also provide us with a transition point between the chemical world of RNA replicators and the living world of RNA agents that actively differentiate self from non-self and generate group identity with membership roles. Importantly the social force of a consortia to solve complex, multilevel problems also depend on using opposing and minority functions. The consortial action of social networks of RNA stem-loops subsequently lead to the evolution of cellular organisms representing a tree of life.

Comparison and Correlation of Genetic Variability of the HBV Pre-S Region in HIV/HBV Co-Infected Patients: Quasispecies Perspective

Background

Human immunodeficiency virus (HIV)/hepatitis B virus (HBV) co-infection can accelerate HBV-induced liver disease. A previous study showed that variation in the HBV pre-S region and quasispecies heterogeneity (Sn, mean genetic distance, dS, dN, and dS/dN) are both related to HBV-induced terminal liver disease in HBV mono-infection. Currently, data are lacking on quasispecies variation of the HBV pre-S region in HIV/HBV co-infection. Investigating the quasispecies variation of the HBV pre-S region and its related factors in HIV/HBV co-infection will help to better explore the pathogenic mechanism of HIV/HBV co-infection.

Methods

According to the HIV antibody results obtained before treatment, chronic HBV-infected patients were divided into HIV/HBV co-infected and HBV mono-infected groups. The clinical characteristics of all patients were collected, and DNA was extracted from the serum. The HBV pre-S region was amplified by nested PCR and was further TA cloned. BioEdit software 7.0 was used for sequence alignment with reference to the standard sequence of the matched HBV genotype. We used 1:1 propensity score matching (PSM) to control for baseline confounding factors between the two groups.

Results

After 1:1 PSM, we identified 100 patients with similar propensities: 50 HIV/HBV co-infected patients and 50 HBV mono-infected patients. HBV quasispecies indices were lower in the HIV/HBV co-infected group than those in the HBV mono-infected group. A significant correlation was observed between all quasispecies indices and soluble cluster of differentiation 163 (sCD163) and interleukin-18 (IL-18) in the HIV/HBV co-infected group; however, this phenomenon was not found in the HBV mono-infected group.

Conclusion

Combined HIV infection reduces quasispecies heterogeneity in the HBV pre-S region, and the quasispecies heterogeneity is related to the sCD163 and IL-18 levels.

Two Sides of a Coin: a Zika Virus Mutation Selected in Pregnant Rhesus Macaques Promotes Fetal Infection in Mice but at a Cost of Reduced Fitness in Nonpregnant Macaques and Diminished Transmissibility by Vectors

Although Zika virus infection of pregnant women can result in congenital Zika syndrome, the factors that cause the syndrome in some but not all infected mothers are still unclear. We identified a mutation that was present in some ZIKV genomes in experimentally inoculated pregnant rhesus macaques and their fetuses. Although we did not find an association between the presence of the mutation and fetal death, we performed additional studies with ZIKV with the mutation in nonpregnant macaques, pregnant mice, and mosquitoes. We observed that the mutation increased the ability of the virus to infect mouse fetuses but decreased its capacity to produce high levels of virus in the blood of nonpregnant macaques and to be transmitted by mosquitoes. This study shows that mutations in mosquito-borne viruses like ZIKV that increase fitness in pregnant vertebrates may not spread in outbreaks when they compromise transmission via mosquitoes and fitness in nonpregnant hosts.

Although fetal death is now understood to be a severe outcome of congenital Zika syndrome, the role of viral genetics is still unclear. We sequenced Zika virus (ZIKV) from a rhesus macaque fetus that died after inoculation and identified a single intrahost substitution, M1404I, in the ZIKV polyprotein, located in nonstructural protein 2B (NS2B). Targeted sequencing flanking position 1404 in 9 additional macaque mothers and their fetuses identified M1404I at a subconsensus frequency in the majority (5 of 9, 56%) of animals and some of their fetuses. Despite its repeated presence in pregnant macaques, M1404I has occurred rarely in humans since 2015. Since the primary ZIKV transmission cycle is human-mosquito-human, mutations in one host must be retained in the alternate host to be perpetuated. We hypothesized that ZIKV I1404 increases viral fitness in nonpregnant macaques and pregnant mice but is less efficiently transmitted by vectors, explaining its low frequency in humans during outbreaks. By examining competitive fitness relative to that of ZIKV M1404, we observed that ZIKV I1404 produced lower viremias in nonpregnant macaques and was a weaker competitor in tissues. In pregnant wild-type mice, ZIKV I1404 increased the magnitude and rate of placental infection and conferred fetal infection, in contrast to ZIKV M1404, which was not detected in fetuses. Although infection and dissemination rates were not different, Aedes aegypti mosquitoes transmitted ZIKV I1404 more poorly than ZIKV M1404. Our data highlight the complexity of arbovirus mutation-fitness dynamics and suggest that intrahost ZIKV mutations capable of augmenting fitness in pregnant vertebrates may not necessarily spread efficiently via mosquitoes during epidemics.

IMPORTANCE Although Zika virus infection of pregnant women can result in congenital Zika syndrome, the factors that cause the syndrome in some but not all infected mothers are still unclear. We identified a mutation that was present in some ZIKV genomes in experimentally inoculated pregnant rhesus macaques and their fetuses. Although we did not find an association between the presence of the mutation and fetal death, we performed additional studies with ZIKV with the mutation in nonpregnant macaques, pregnant mice, and mosquitoes. We observed that the mutation increased the ability of the virus to infect mouse fetuses but decreased its capacity to produce high levels of virus in the blood of nonpregnant macaques and to be transmitted by mosquitoes. This study shows that mutations in mosquito-borne viruses like ZIKV that increase fitness in pregnant vertebrates may not spread in outbreaks when they compromise transmission via mosquitoes and fitness in nonpregnant hosts.

Limited evolution of the yellow fever virus 17d in a mouse infection model

By infecting mice with the yellow fever virus vaccine strain 17D (YFV-17D; Stamaril®), the dose dependence and evolutionary consequences of neurotropic yellow fever infection was assessed. Highly susceptible AG129 mice were used to allow for a maximal/unlimited expansion of the viral populations. Infected mice uniformly developed neurotropic disease; the virus was isolated from their brains, plaque purified and sequenced. Viral RNA populations were overall rather homogenous [Shannon entropies 0−0.15]. The remaining, yet limited intra-host population diversity (0−11 nucleotide exchanges per genome) appeared to be a consequence of pre-existing clonal heterogeneities (quasispecies) of Stamaril®. In parallel, mice were infected with a molecular clone of YFV-17D which was in vivo launched from a plasmid. Such plasmid-launched YFV-17D had a further reduced and almost clonal evolution. The limited intra-host evolution during unrestricted expansion in a highly susceptible host is relevant for vaccine and drug development against flaviviruses in general. Firstly, a propensity for limited evolution even upon infection with a (very) low inoculum suggests that fractional dosing as implemented in current YF-outbreak control may pose only a limited risk of reversion to pathogenic vaccine-derived virus variants. Secondly, it also largely lowers the chance of antigenic drift and development of resistance to antivirals.

Hepatitis C reinfection with protease inhibitor-resistant hepatitis C virus in an HIV-coinfected MSM

There is an ongoing global hepatitis C virus (HCV) epidemic mainly related to high-risk behaviour in persons who inject drugs (PWID) and in HIV-infected men who have sex with men (MSM) which continues to fuel the HCV epidemic. Treatment of HCV infection with direct antiviral therapy (DAA) has been very successful in the last decade. Main obstacles for HCV elimination are HCV reinfections observed in PWID and HIV-infected MSM. We present here an HIV-infected MSM patient who has been reinfected thrice with HCV. The virus which was investigated from his last reinfection episode reveals transmission of a newly acquired HCV protease inhibitor (PI) resistance, despite not having been exposed to HCV-PIs during his last DAA therapy.

The genetic variability of hepatitis B virus subgenotype F1b precore/core gene is related to the outcome of the acute infection

AIM

To assess the association of viral and host genetic variability with the outcome of acute infection with hepatitis B virus subgenotype F1b (HBV/F1b).

METHODS

The cohort consisted of 26 patients with acute HBV/F1b infection who exhibit different outcomes: spontaneous resolution (n = 10), progression to chronic hepatitis (n = 10) and acute liver failure (n = 6). HLA SNPs (rs3077, rs9277542, rs2856718 and rs7453920) were determined. The S gene and core promoter/precore/core region were direct sequenced, and this latter region was also ultra-deep sequenced. Mean number of mutations, mutation rate, Shannon entropy, positive selection sites and mutational patterns of quasispecies were compared between groups.

RESULTS

HLA SNPs were associated with spontaneous resolution or progression to chronic hepatitis, but not with the development of acute liver failure. The mean number of mutations in the S gene was similar among the three groups. Patients with spontaneous resolution had the lowest number of mutations, mutation rates and Shannon entropy values in the precore/core compared to the other two groups. Ten positive selection sites mapped on HLA-restricted epitopes were related to progression to chronic hepatitis and acute liver failure. Mutations T1753C, A1762T, G1764A, C1766T, T1768A G1896A, G2092T and T2107C were associated with acute liver failure and progression to chronic hepatitis.

CONCLUSION

Highly heterogeneous and complex HBV precore/core carrying specific point mutations, combined with the host HLA background, were associated with a worse clinical outcome of acute HBV/F1b infection.

Attenuation of Live-Attenuated Yellow Fever 17D Vaccine Virus Is Localized to a High-Fidelity Replication Complex

Live-attenuated viral vaccines are highly safe and efficacious but represent complex and often multigenic attenuation mechanisms. Most of these vaccines have been generated empirically by serial passaging of a wild-type (WT) virus in cell culture. One of the safest and most effective live-attenuated vaccines is yellow fever (YF) virus strain 17D, which has been used for over 80 years to control YF disease. The availability of the WT parental strain of 17D, Asibi virus, and large quantities of clinical data showing the effectiveness of the 17D vaccine make this WT parent/vaccine pair an excellent model for investigating RNA virus attenuation. Here, we investigate a mechanism of 17D attenuation and show that the vaccine virus is resistant to the antiviral compound ribavirin. The findings suggest that attenuation is in part due to a low probability of reversion or mutation of the vaccine virus genome to WT, thus maintaining a stable genotype despite external pressures.

The molecular basis of attenuation for live-attenuated vaccines is poorly understood. The yellow fever (YF) 17D vaccine virus was derived from the wild-type, parental strain Asibi virus by serial passage in chicken tissue and has proven to be a very safe and efficacious vaccine. We have previously shown that wild-type Asibi is a typical RNA virus with high genetic diversity, while the 17D vaccine virus has very little genetic diversity. To investigate this further, we treated Asibi and 17D viruses with ribavirin, a GTP analog with strong antiviral activity that increases levels of mutations in the viral genome. As expected, ribavirin treatment introduced mutations into the Asibi virus genome at a very high frequency and decreased viral infectivity while, in contrast, the 17D vaccine virus was resistant to ribavirin, as treatment with the antiviral introduced very few mutations into the genome, and viral infectivity was not lost. The results were confirmed for another YF wild-type parental and vaccine pair, a wild-type French viscerotropic virus and French neurotropic vaccine. Using recombinant Asibi and 17D viruses, ribavirin sensitivity was located to viral nonstructural genes. Thus, two live-attenuated YF vaccine viruses are genetically stable even under intense mutagenic pressure, suggesting that attenuation of live-attenuated YF vaccines is due, at least in part, to fidelity of the replication complex resulting in high genetic stability.

Next generation sequencing identifies baseline viral mutants associated with treatment response to pegylated interferon in HBeAg-positive chronic hepatitis B

Current data of hepatitis B virus (HBV) variants associated with treatment outcome identified by next generation sequencing (NGS) are limited. This study was aimed at determining the role of baseline sequence variations in the enhancer II (EnhII), basal core promotor (BCP) and pre-core (PC) regions of HBV genotype C in patients treated with pegylated interferon (PEG-IFN). Patients with HBeAg-positive chronic hepatitis B (CHB) treated with 48-week PEG-IFN were enrolled. Combined response (CR) at week 96 was defined by HBeAg seroconversion plus HBV DNA < 2000 IU/mL and HBsAg < 1000 IU/mL. Pre-treatment viral mutations were characterized by Sanger sequencing and NGS (Miseq Illumina platform). Among 47 patients (32 male, mean age 32.4 years), CR was achieved in 12 (25.5%) individuals. Overall, NGS was superior to Sanger sequencing in detecting mutations (61.7% vs. 38.3%, P < 0.001). Based on NGS, the prevalence of T1753V (T1753C/A/G) and A1762T/G1764A variants were significantly lower in responders compared to non-responders (8.3% vs. 51.4%, P = 0.009 and 33.3% vs. 68.6%, P = 0.032, respectively). No significant difference between groups was found regarding C1653T and G1896A mutants. The absence of T1753V and A1762T/G1764A mutations were factors associated with CR (OR 11.65, 95%CI 1.36-100.16, P = 0.025, and OR 4.36, 95%CI 1.08-17.63, P = 0.039, respectively). The existence of pre-treatment T1753V, A1762T/G1764A mutations and their combination yielded negative predictive values of 94.7%, 85.7% and 93.8%, respectively. The presence of HBV mutants in the BCP region determined by NGS at baseline was associated with poor treatment outcome in patients with HBeAg-positive CHB receiving PEG-IFN.

Viral quasispecies of hepatitis B virus in patients with YMDD mutation and lamivudine resistance may not predict the efficacy of lamivudine/adefovir rescue therapy

The association between hepatitis B virus (HBV) quasispecies (QS) and the efficacy of nucleos(t)ide analog therapy is currently not well defined, particularly in the case of lamivudine (LAM)/adefovir (ADV) combination rescue therapy for patients with chronic HBV infection (CHB) presenting with LAM resistance. In the present study, 16 CHB patients with the rtM204I/V mutation in the tyrosine-methionine-aspartate-aspartate motif of the C domain of the polymerase gene who switched to LAM/ADV treatment due to LAM resistance were assessed. HBV DNA was isolated from these patients and the reverse transcriptase (RT) region was sequenced. The QS heterogeneity and distribution was analyzed, the mutation sites were recorded and the phylogenetic trees were constructed. The results indicated that QS heterogeneity and distribution in the RT and S regions were not significantly different between responders (RS) and non-RS (NRS) at baseline (P>0.05), except for the higher frequency of a dominant strain in the RT region at the nucleotide level in the RS group (P=0.039). In addition, in NRS, no significant difference in QS heterogeneity or distribution in these regions was identified at six months vs. the baseline. Furthermore, although in the non-responder group the frequency of the LAM resistance-associated mutations (rtM204V/I) decreased at 6 months compared with the baseline, it did not disappear in any of the patients after six months of treatment. Analysis of individual patients did not indicate any consistent selection of specific HBV mutants during LAM/ADV rescue therapy. In conclusion, the baseline HBV QS within the RT and S regions may not be a valid predictor of the response to LAM/ADV rescue treatment in CHB patients with LAM resistance.

Clinical Manifestations and Laboratory Tests of AECHB and Severe Hepatitis (Liver Failure)

This chapter describes the clinical symptoms and signs of AECHB and HBV ACLF, classification, grading of HBV ACLF and their features, diagnostic principles and standards in liver pathology, biochemistry, and virology of HBV ACLF.

Liver failure is defined as serious damage to the liver cause by a variety of etiologies, leading to liver function disorder or even decompensation, and clinical syndromes with coagulopathy, jaundice, hepatic encephalopathy, and ascites.

Severe hepatitis B can be indicated pathologically by apparent hepatocellular necrosis, including extensive multifocal, confluent, bridging, sub-massive or massive necrosis.

Laboratory tests during the course of severe exacerbation of chronic hepatitis B can reflect pathological changes and liver function in a timely manner, providing objective and informative reference data for evaluation of disease severity and treatment efficacy. Among the most important laboratory tests are those for prothrombin activity, international normalized ratio, and increases in total bilirubin concentration.

Severe hepatitis B is associated with interactions between the virus and host factors. Detection of HBV DNA, HBV genotype, quasispecies and HBV mutation can provide important theoretical bases for the prevention, control or mitigation of the progress of severe hepatitis B.

Noninvasive imaging modalities can be used to visualize the entire liver and parts of it. Measuring liver volume to evaluate liver size and liver reserve capacity is regarded as important in diagnosis, surgical approach and prognostic evaluation of patients with severe exacerbation of chronic hepatitis B and liver failure.

Model for End-Stage Liver Disease (MELD) is the first quantitative method developed to assess whether a patient with liver failure requires a liver transplant. The predictive value of the MELD model has been improved by the MELD-Na, iMELD, and MESO models. Several other valuable prognostic models have been developed. For example, for patients with HBV-ACLF, the established TPPM scoring system was found to be more predictive than MELD score.

Next-generation sequencing analysis of new genotypes appearing during antiviral treatment of chronic hepatitis C reveals that these are selected from pre-existing minor strains

Coinfection with more than one hepatitis C virus (HCV) genotype is common, but its dynamics, particularly during antiviral treatment, remain largely unknown. We employed next-generation sequencing (NGS) to analyse sequential serum and peripheral blood mononuclear cell (PBMC) samples in seven patients with transient presence or permanent genotype change during antiviral treatment with interferon and ribavirin. Specimens were collected right before the therapy initiation and at 2, 4, 6, 8, 12, 20, 24, 36, 44 and 48 weeks during treatment and 6 months after treatment ceased. A mixture of two different genotypes was detected in the pretreatment samples from five patients and the minor genotype constituted 0.02 to 38 %. A transient or permanent change of the predominant genotype was observed in six patients. In three cases genotype 3 was replaced as the predominant genotype by genotype 4, in two cases genotype 3 was replaced by genotype 1, and in one subject genotype 1 was replaced by genotype 4. The PBMC- and serum-derived sequences were frequently discordant with respect to genotype and/or genotype proportions. In conclusion, pre-existing minor HCV genotypes can be selected rapidly during antiviral treatment and become transiently or permanently predominant. In coinfections involving genotype 3, genotype 3 was eliminated first from both the serum and PBMC compartments. The PBMC- and serum-derived HCV sequences were frequently discordant with respect to genotype and/or genotype proportions, suggesting that they constitute separate compartments with their own dynamics.

Characterization of Hepatitis C Virus IRES Quasispecies - From the Individual to the Pool

Hepatitis C virus (HCV) is a single-stranded positive-sense RNA virus from the genus Hepacivirus. The viral genomic +RNA is 9.6 kb long and contains highly structured 5′ and 3′ untranslated regions (UTRs) and codes for a single large polyprotein, which is co- and post-translationally processed by viral and cellular proteases into at least 11 different polypeptides. Most of the 5′ UTR and an initial part of the polyprotein gene are occupied by an internal ribosome entry site (IRES), which mediates cap-independent translation of the viral proteins and allows the virus to overcome cellular antiviral defense based on the overall reduction of the cap-dependent translation initiation. We reconsidered published results concerning a search for possible correlation between patient response to interferon-based antiviral therapy and accumulation of nucleotide changes within the HCV IRES. However, we were unable to identify any such correlation. Rather than searching for individual mutations, we suggest to focus on determination of individual and collective activities of the HCV IRESs found in patient specimens. We developed a combined, fast, and undemanding approach based on high-throughput cloning of the HCV IRES species to a bicistronic plasmid followed by determination of the HCV IRES activity by flow cytometry. This approach can be adjusted for measurement of the individual HCV IRES activity and for estimation of the aggregate ability of the whole HCV population present in the specimen to synthesize viral proteins. To detect nucleotide variations in the individual IRESs, we used denaturing gradient gel electrophoresis (DGGE) analysis that greatly improved identification and classification of HCV IRES variants in the sample. We suggest that determination of the collective activity of the majority of HCV IRES variants present in one patient specimen in a given time represents possible functional relations among variant sequences within the complex population of viral quasispecies better than bare information about their nucleotide sequences. A similar approach might be used for monitoring of sequence variations in quasispecies populations of other RNA viruses in all cases when changes in primary sequence represent changes in measurable and easily quantifiable phenotypes.

A new hepatitis B virus e antigen-negative strain gene used as a reference sequence in an animal model

Infection with hepatitis B virus (HBV) e-antigen (HBeAg)-negative strains is increasingly prevalent. Currently, detailed information of the obtained natural HBV strain is not available except for the B genotype and HBeAg-negative. The aim of the present study was to characterize the natural genetic variation of the HBeAg-negative strain and investigate its function. The genic sequence was determined using Sanger sequencing, and compared to related sequences using alignment and phylogenetic analysis. In vivo, virus-specific serum markers were investigated in CBA/CaJ mice. The sequence had a full genome length of 3215 nucleotides. Sites 122, 125, 127, and 160 in S regions were identified as lysine, threonine, proline, and lysine respectively. The main four point variants including A1762T, G1764A, G1896A, and G1899A were detected in the full-length genome. The genotype of the sequence was B, with sub-genotype B2 and serological subtype adw2. The characterize of the natural genetic variation strain showed no reported drug-resistant variant in P region and no reported immune escape site in S region. The strain will increase viral replication and infection for mutations A1762T and G1764A in the basal core promoter region, and mutations G1896A and G1899A in the pre-core region. The G1896A variant resulted in a premature stop codon and abolished HBeAg expression. HBsAg persisted for 26 weeks and HBeAg was still negative in CBA/CaJ mice. The present sequence is representative of the HBeAg-negative genome and may serve as a valuable reference for studying HBeAg-negative strains. The present findings were successfully verified in CBA/CaJ mice, demonstrating good applicability of the sequence.

Molecular characterization of hepatitis C virus in end-stage renal disease patients under hemodialysis

New direct-acting antiviral (DAA) agents are in development or already approved for the treatment of chronic hepatitis C virus (HCV) infection. The effectiveness of these drugs is related to the previous existence of resistant variants. Certain clinical conditions can allow changes in immunological characteristics of the host and even modify genetic features of viral populations. The aim of this study was to perform HCV molecular characterization from samples of end-stage renal disease patients on hemodialysis (ESRD-HD). Nested PCR and Sanger sequencing were used to obtain genetic information from the NS5B partial region of a cohort composed by 86 treatment-naïve patients. Genomic sequences from the Los Alamos databank were employed for comparative analysis. Bioinformatics methodologies such as phylogenetic reconstructions, informational entropy, and mutation analysis were used to analyze datasets separated by geographical location, HCV genotype, and renal function status. ESRD-HD patients presented HCV genotypes 1a (n = 18), 1b (n = 16), 2a (n = 2), 2b (n = 2), and 3a (n = 4). Control subjects were infected with genotypes 1a (n = 11), 1b (n = 21), 2b (n = 4), and 3a (n = 8). Dataset phylogenetic reconstruction separated HCV subtype 1a into two distinct clades. The entropy analysis from the ESRD-HD group revealed two amino acid positions related to an epitope for cytotoxic T lymphocytes and T helper cells. Genotype 1a was found to be more diverse than subtype 1b. Also, genotype 1a ERSD-HD patients had a higher mean of amino acids changes in comparison to control group patients. The identification of specific mutations on epitopes and high genetic diversity within the NS5B HCV partial protein in hemodialysis patients can relate to host immunological features and geographical distribution patterns. This genetic diversity can affect directly the new DAA's resistance mechanisms.

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.

Deep sequencing analysis of quasispecies in the HBV pre-S region and its association with hepatocellular carcinoma

BACKGROUND

The association between the evolution of hepatitis B virus (HBV) quasispecies and the development of hepatocellular carcinoma (HCC) is unknown.

METHODS

We used deep sequencing to examine the dynamics of HBV quasispecies and their relationship to HCC development. Thirty-two chronic hepatitis B (CHB) patients with HCC (HCC group) and 32 matched CHB patients without HCC (controls) were recruited. Fourteen patients from each group had serial sera available up to 9 years before the time of the present study. Deep sequencing of the HBV pre-S regions was performed. HBV quasispecies complexity, diversity, and intrapatient prevalence of pre-S deletions/mutations were analyzed.

RESULTS

Compared with control patients, HCC patients had a significant greater quasispecies complexity (p = 0.04 at the nucleotide level), greater diversity (p = 0.004 and 0.009 at the nucleotide level and the amino acid level respectively), and a trend of greater complexity at the amino acid level (p = 0.065). HCC patients had a higher intrapatient prevalence of pre-S deletions and point mutations (at codons 4, 27, and 167) compared with the control patients (all p < 0.05). Longitudinal observation in the sera of 14 HCC patients showed that quasispecies complexity (p = 0.027 and 0.024 at the nucleotide level and the amino acid level respectively) and diversity (p = 0.035 and 0.031 at the nucleotide level and the amino acid level respectively) increased as the disease progressed to HCC.

CONCLUSIONS

Increased HBV quasispecies complexity and diversity in the pre-S region, probably reflecting enhanced virus-host interplay, was associated with disease progression from CHB to HCC.

Comparison of error correction algorithms for Ion Torrent PGM data: application to hepatitis B virus

Ion Torrent Personal Genome Machine (PGM) technology is a mid-length read, low-cost and high-speed next-generation sequencing platform with a relatively high insertion and deletion (indel) error rate. A full systematic assessment of the effectiveness of various error correction algorithms in PGM viral datasets (e.g., hepatitis B virus (HBV)) has not been performed. We examined 19 quality-trimmed PGM datasets for the HBV reverse transcriptase (RT) region and found a total error rate of 0.48% ± 0.12%. Deletion errors were clearly present at the ends of homopolymer runs. Tests using both real and simulated data showed that the algorithms differed in their abilities to detect and correct errors and that the error rate and sequencing depth significantly affected the performance. Of the algorithms tested, Pollux showed a better overall performance but tended to over-correct 'genuine' substitution variants, whereas Fiona proved to be better at distinguishing these variants from sequencing errors. We found that the combined use of Pollux and Fiona gave the best results when error-correcting Ion Torrent PGM viral data.

HBV DNA genome co-transfection procedure for the evaluation of relative fitness

Hepatitis B virus (HBV) has a high mutation rate and exists as a mixture of genetically different but closely related variants. We present a HBV DNA co-transfection fitness assay and use it to evaluate the relative fitness of different HBV variants in two scenarios: seroconversion process and occupation of an ecological niche. In the seroconversion experiment, subgenotype D1 (sgtD1) deletion (1763-1770) had significantly lower fitness comparing with both sgtD1 wild type and sgtD1mut G1896A, while, in the case of occupation of ecological niche experiment, the results showed the same relative fitness between all of the genotype combinations, except F1b-F4. In this case sgtF1b clearly overgrow sgtF4, which is in accordance with the observation that F1b is the most prevalent in the new infections in Argentina. In summary, we present a method aimed to evaluate HBV viral fitness which improve the analysis of the relative frequency of viral variants during the HBV infection process.

Evolution of full-length genomes of HBV quasispecies in sera of patients with a coexistence of HBsAg and anti-HBs antibodies

Although the evolutionary changes of viral quasispecies are correlated to the pathological status of a disease, little is known in the coexistence of hepatitis B surface antigen (HBsAg) and antibodies to these antigens (anti-HBs). To examine evolutionary changes in hepatitis B virus (HBV) and their relationship to the coexistence of HBsAg and anti-HBs antibodies, HBV genomes in patients with a coexistence of HBsAg and anti-HBs antibodies (experimental group) and HBsAg positive without anti-HBs (control group) were assessed. Our results showed that quasispecies diversity was significantly higher in the experimental group for large HBsAg (LHBsAg), middle HBsAg (MHBsAg), and HBsAg genes. LHBsAg harbored dN/dS values eight times higher in the experimental group; however, the mean dN/dS ratios in genes HbxAg, Pol and PreC/C of the experimental patients had an opposite trend. Phylogenetic trees in the experimental group were more complex than the control group. More positive selection sites, mutations and deletions were observed in the experimental group in specific regions. Furthermore, several amino acid variants in epitopes were potentially associated with the immune evasion. In conclusion, cumulative evolutionary changes in HBV genome that facilitate immune evasion provide insights into the genetic mechanism of a coexistence of HBsAg and anti-HBs antibodies.

Evidence of bottleneck effect on hepatitis C virus transmission between a couple under interferon based therapy

Issues on the correlation of viral genetic diversity and treatment response to the hepatitis C infection remain uncertain. The bottleneck effect dictates the characteristics of the viral population that will establish the infection in a new host and is related to how the immune system and treatment will be effective against the virus. Here we evaluated the phylogenetic characteristics of quasispecies population and the treatment response pattern of a HCV infected couple. We also analyzed whether the viral population of these patients indicated that they were exposed to the same source for primer infection. This study included two patients (P10 and P11) HCV genotype 1b infected. The couple presented horizontal transmission. Viral RNA was isolated from serum samples collected before, during and after treatment, at specific time points. The HCV NS5A gene sequence was amplified, cloned and sequenced. Genetic and evolutionary analyses were performed to compare the quasispecies population of these two patients and local control patients. Genetic distance and diversity were calculated. Phylogenetic analyses were performed by using maximum likelihood and Bayesian methodologies. The analysis of the baseline samples showed that the genetic distance of the viral populations of patients P10 and P11 was significantly lower than when these patients and the control group based on sequences from local patients were analyzed, supporting the horizontal transmission hypothesis. Phylogenetic analysis with sequences from all the time point samples also demonstrated two patterns of evolution depending on the treatment response. The Bayesian analysis showed that one isolate corresponding to the baseline sample of P10 was grouped into the P11 clade, suggesting a way of infection and a bottleneck effect. Our data suggests that the patient P11 viral population may be originated from variants from P10 patient and consequently showing that clinical differences between treatment responses can emerge from the bottleneck effect on viral populations.

Hepatitis B virus full-length genomic mutations and quasispecies in hepatocellular carcinoma

BACKGROUND AND AIM

Hepatitis B virus (HBV) full-length genomic mutations and quasispecies characteristics in hepatocellular carcinoma (HCC) were investigated.

METHODS

Hepatitis B virus DNA was extracted from the tumor and non-tumor tissues of 16 HCC patients. Overlapping DNA fragments covering the entire HBV genome were amplified and sequenced. To study HBV sequence at the quasispecies level, the preS region was amplified and clonally sequenced. HBV mutation profiles, quasispecies complexity and diversity, and phylogenetic characteristics were assessed.

RESULTS

Fourteen patients had full-length HBV amplification. Hot-spot mutations at HBx aa130-131 and pre-S deletions were detected in 13 (93%) and 6 (43%) patients, respectively. Deletions in the X/preC/C regions were more frequently detected in the tumor than in the non-tumor tissues (P = 0.031). Compared with the non-tumor tissues, the tumor tissues had a lower quasispecies complexity (P = 0.014 and 0.043, at the nucleotide and amino acid levels, respectively) and diversity (P = 0.048 and 0.022, at the nucleotide and amino acid levels, respectively). Phylogenetic analysis showed that HBV sequences derived from tumor and non-tumor tissues were separately clustered, suggesting the occurrence of compartmentalization, which was confirmed by the correlation coefficient testing on both the number and length of branches of viral populations (all P < 0.02).

CONCLUSIONS

Hepatitis B virus mutation patterns in HCC tumor tissues and non-tumor tissues were different. HBV quasispecies within the preS region were compartmentalized, and tumor tissues had a lower genome complexity and diversity. Our study suggests HBV evolution is conditioned by the differential host cellular environment in HCC tumors.

Current molecular methods for the detection of hepatitis B virus quasispecies

Chronic HBV infection affects more than 240 million people worldwide and is associated with a broad range of clinical manifestations including liver cirrhosis, liver failure and hepatocellular carcinoma. Because of the lack of an efficient cure for chronic hepatitis B, the main goal of antiviral therapy is the prevention of liver disease progression coupled with prolonged survival of patients. Because HBV viral load has been shown to be a crucial determinant of the progression of liver damage, these goals can be achieved as long as HBV replication can be suppressed. Unfortunately, long-term therapy with the low-to-moderate genetic barrier drugs, which are still recommended in a majority of developing countries, are strongly associated with HBV resistance development and treatment failure. In such cases, the precise and accurate determination of drug-resistant variants in an individual patient before treatment is important for a proper choice of first-line potent therapy. Nowadays, a number of techniques are available to study HBV quasispecies evolution. This review describes the advantages and limitations of various assays detecting drug-resistant HBV variants. Copyright © 2016 John Wiley & Sons, Ltd.

Hepatitis C virus: Promising discoveries and new treatments

Despite advances in therapy, hepatitis C virus (HCV) infection remains an important global health issue. It is estimated that a significant part of the world population is chronically infected with the virus, and many of those affected may develop cirrhosis or liver cancer. The virus shows considerable variability, a characteristic that directly interferes with disease treatment. The response to treatment varies according to HCV genotype and subtype. The continuous generation of variants (quasispecies) allows the virus to escape control by antivirals. Historically, the combination of ribavirin and interferon therapy has represented the only treatment option for the disease. Currently, several new treatment options are emerging and are available to a large part of the affected population. In addition, the search for new substances with antiviral activity against HCV continues, promising future improvements in treatment. Researchers should consider the mutation capacity of the virus and the other variables that affect treatment success.

Natural Polymorphisms Conferring Resistance to HCV Protease and Polymerase Inhibitors in Treatment-Naïve HIV/HCV Co-Infected Patients in China

BACKGROUND

The advent of direct-acting agents (DAAs) has improved treatment of HCV in HIV co-infection, but may be limited by primary drug resistance. This study reports the prevalence of natural polymorphisms conferring resistance to NS3/4A protease inhibitors and NS5B polymerase inhibitors in treatment-naïve HIV/HCV co-infected individuals in China.

METHODS

Population based NS3/4A sequencing was completed for 778 treatment-naïve HIV/HCV co-infected patients from twelve provinces. NS3 sequences were amplified by nested PCR using in-house primers for genotypes 1-6. NS5B sequencing was completed for genotyping in 350 sequences. Resistance-associated variants (RAVs) were identified in positions associated with HCV resistance.

RESULTS

Overall, 72.8% (566/778) of all HCV sequences had at least one RAV associated with HCV NS3/4A protease inhibitor resistance. Variants were found in 3.6% (7/193) of genotype 1, 100% (23/23) of genotype 2, 100% (237/237) of genotype 3 and 92% (299/325) of genotype 6 sequences. The Q80K variant was present in 98.4% of genotype 6a sequences. High-level RAVs were rare, occurring in only 0.8% of patients. 93% (64/69) patients with genotype 1b also carried the C316N variant associated with NS5B low-level resistance.

CONCLUSIONS

The low frequency of high-level RAVs associated with primary HCV DAA resistance among all genotypes in HIV/HCV co-infected patients is encouraging. Further phenotypic studies and clinical research are needed.

The Amino Acid Substitution Q65H in the 2C Protein of Swine Vesicular Disease Virus Confers Resistance to Golgi Disrupting Drugs

Swine vesicular disease virus (SVDV) is a porcine pathogen and a member of the species Enterovirus B within the Picornaviridae family. Brefeldin A (BFA) is an inhibitor of guanine nucleotide exchange factors of Arf proteins that induces Golgi complex disassembly and alters the cellular secretory pathway. Since BFA has been shown to inhibit the RNA replication of different enteroviruses, including SVDV, we have analyzed the effect of BFA and of golgicide A (GCA), another Golgi disrupting drug, on SVDV multiplication. BFA and GCA similarly inhibited SVDV production. To investigate the molecular basis of the antiviral effect of BFA, SVDV mutants with increased resistance to BFA were isolated. A single amino acid substitution, Q65H, in the non-structural protein 2C was found to be responsible for increased resistance to BFA. These results provide new insight into the relationship of enteroviruses with the components of the secretory pathway and on the role of SVDV 2C protein in this process.

Geneticin Stabilizes the Open Conformation of the 5' Region of Hepatitis C Virus RNA and Inhibits Viral Replication

The aminoglycoside Geneticin (G418) is known to inhibit cell culture proliferation, via virus-specific mechanisms, of two different virus genera from the family Flaviviridae. Here, we tried to determine whether Geneticin can selectively alter the switching of the nucleotide 1 to 570 RNA region of hepatitis C virus (HCV) and, if so, whether this inhibits viral growth. Two structure-dependent RNases known to specifically cleave HCV RNA were tested in the presence or absence of the drug. One was the Synechocystis sp. RNase P ribozyme, which cleaves the tRNA-like domain around the AUG start codon under high-salt buffer conditions; the second was Escherichia coli RNase III, which recognizes a double-helical RNA switch element that changes the internal ribosome entry site (IRES) from a closed (C) conformation to an open (O) one. While the drug did not affect RNase P activity, it did inhibit RNase III in the micromolar range. Kinetic studies indicated that the drug favors the switch from the C to the O conformation of the IRES by stabilizing the distal double-stranded element and inhibiting further processing of the O form. We demonstrate that, because the RNA in this region is highly conserved and essential for virus survival, Geneticin inhibits HCV Jc1 NS3 expression, the release of the viral genomic RNA, and the propagation of HCV in Huh 7.5 cells. Our study highlights the crucial role of riboswitches in HCV replication and suggests the therapeutic potential of viral-RNA-targeted antivirals.

Dynamics of Genotypic Mutations of the Hepatitis B Virus Associated With Long-Term Entecavir Treatment Determined With Ultradeep Pyrosequencing: A Retrospective Observational Study

The aim of the study is to explore the evolution of genotypic mutations within the reverse transcriptase region in partial virological responders (PVRs) receiving long-term entecavir (ETV) treatment. A total of 32 patients were classified as completely virological responders (CVRs) (n = 12) or PVRs (n = 20). Five partial responders were hepatitis B virus (HBV)-DNA positive after long-term therapy, which lasted for >3 years. A total of 71 serum samples from these 32 patients were assayed by ultra-deep pyrosequencing (UDPS): 32 samples were from all patients at baseline, and 39 were from PVRs with sequential inter-treatment. Approximately 84,708 sequences were generated per sample. At baseline, the quasispecies heterogeneity did not significantly differ between the 2 groups. The frequencies of substitutions indicating pre-existence of nucleos(t)ide analog resistant (NAr) mutants ranged from 0.10% to 6.70%, which did not statistically differ between groups either. However, the substitutions associated with the NAr mutants were significantly different from those associated with the non-NAr mutants in 13 patients; 6 of these patients were PVRs and the others were CVRs. Five patients were HBV DNA positive after regular ETV monotherapy for >3 years, and 4 of these patients underwent mild NAr substitution fluctuations (<20%). One patient developed virological breakthrough while bearing single, double, and triple (rtL180 M, rtM204 V, rtS202G) substitutions. In addition to the common substitutions, unknown amino acid substitutions, such as rtL145 M/S, rtF151Y/L, rtR153Q, rtI224 V, rtN248H, rtS223A, rtS256C, need to be further verified. NAr substitutions are observed at frequencies of 0.10% to 6.7% before therapy. Long-term ETV therapy generally results in virological responses, as long as the proportion of resistance mutations remains at a relatively low level. Genotypic resistance to ETV is detected in all PVRs receiving long-term ETV therapy.

Evolution of Newcastle Disease Virus Quasispecies Diversity and Enhanced Virulence after Passage through Chicken Air Sacs

It has been reported that lentogenic Newcastle disease virus (NDV) isolates have the potential to become velogenic after their transmission and circulation in chickens, but the underlying mechanism is unclear. In this study, a highly velogenic NDV variant, JS10-A10, was generated from the duck-origin lentogenic isolate JS10 through 10 consecutive passages in chicken air sacs. The velogenic properties of this selected variant were determined using mean death time (MDT) assays, intracerebral pathogenicity index (ICPI), the intravenous pathogenicity index (IVPI), histopathology, and the analysis of host tissue tropism. In contrast, JS10 remained lentogenic after 20 serial passages in chicken eggs (JS10-E20). The JS10, JS10-A10, and JS10-E20 genomes were sequenced and found to be nearly identical, suggesting that both JS10-A10 and JS10-E20 were directly generated from JS10. To investigate the mechanism for virulence enhancement, the partial genome covering the F0 cleavage site of JS10 and its variants were analyzed using ultradeep pyrosequencing (UDPS) and the proportions of virulence-related genomes in the quasispecies were calculated. Velogenic NDV genomes accumulated as a function of JS10 passaging through chicken air sacs. Our data suggest that lentogenic NDV strains circulating among poultry might be a risk factor to future potential velogenic NDV outbreaks in chickens.

IMPORTANCE

An avirulent isolate, JS10, was passaged through chicken air sacs and embryos, and the pathogenicity of the variants was assessed. A virulent variant, JS10-A10, was generated from consecutive passage in air sacs. We developed a deep-sequencing approach to detect low-frequency viral variants across the NDV genome. We observed that virulence enhancement of JS10 was due to the selective accumulation of velogenic quasispecies and the concomitant disappearance of lentogenic quasispecies. Our results suggest that because it is difficult to avoid contact between natural waterfowl reservoirs and sensitive poultry operations, circulating lentogenic NDV strains may represent a potential reservoir for emergent velogenic NDV strains that could cause outbreaks in chickens.

Increased intrahepatic quasispecies heterogeneity correlates with off-treatment sustained response to nucleos(t)ide analogues in e antigen-positive chronic hepatitis B patients

Finite treatment with nucleos(t)ide analogues (NAs) remains a great challenge for chronic hepatitis B in the clinic. This study aimed to investigate the relationship between intrahepatic quasispecies heterogeneity and the NAs off-treatment outcomes in a prospective cohort. Eighteen HBeAg-positive patients with chronic hepatitis B who achieved the cessation criteria underwent liver biopsy, and stopped treatment thereafter. Patients were followed up prospectively for 1 year. The reverse transcriptase (RT) gene of intrahepatic hepatitis B virus (HBV) was cloned and sequenced. Intrahepatic quasispecies heterogeneity and specific gene mutations were analysed using bioinformatic methods. Ten patients achieved sustained response, and eight patients developed viral relapse. The intrahepatic quasispecies Shannon entropy and nucleotide diversity within either RT or the surface (S) region of patients with sustained response were significantly higher (p < 0.05) than those of patients who had a viral relapse. Intrahepatic quasispecies Shannon entropy at the nucleotide level predicted the sustained off-treatment response (area under receiver operating characteristics curve 0.925; 95% CI 0.807-1.000; p 0.003). More positive selection sites and N-glycosylation mutations within the S region were found in patients with sustained response than in the patients with viral relapse (p < 0.01). Most of the positive selection sites in patients with sustained response were located in reported HLA-I-restricted or HLA-II-restricted epitopes. Intrahepatic quasispecies heterogeneity at the end of treatment was correlated with off-treatment outcomes in HBeAg-positive patients with chronic hepatitis B. More immune escape mutations were found within the S region in patients with sustained response. The higher intrahepatic quasispecies heterogeneity indicated a more robust immune control over HBV, which in turn maintained a sustained response after withdrawal of NAs.

Statistical properties and error threshold of quasispecies on single-peak Gaussian-distributed fitness landscapes

The stochastic Eigen model proposed by Feng et al. (2007) (Journal of Theoretical Biology, 246, 28) showed that error threshold is no longer a phase transition point but a crossover region whose width depends on the strength of the random fluctuation in an environment. The underlying cause of this phenomenon has not yet been well examined. In this article, we adopt a single peak Gaussian distributed fitness landscape instead of a constant one to investigate and analyze the change of the error threshold and the statistical property of the quasi-species population. We find a roughly linear relation between the width of the error threshold and the fitness fluctuation strength. For a given quasi-species, the fluctuation of the relative concentration has a minimum with a normal distribution of the relative concentration at the maximum of the averaged relative concentration, it has however a largest value with a bimodal distribution of the relative concentration near the error threshold. The above results deepen our understanding of the quasispecies and error threshold and are heuristic for exploring practicable antiviral strategies.

Characterization of Full-Length Genomes of Hepatitis B Virus Quasispecies in Sera of Patients at Different Phases of Infection

Hepatitis B virus (HBV) infection results in different clinical presentation due to different levels of immune response. Our study aimed to characterize HBV full-length genome quasispecies (QS) in patients with different phases of infection to better understand its pathogenesis. Forty treatment-naive HBV-infected patients were enrolled, including 10 cases of acute hepatitis B (AHB), 9 cases of immunotolerant (IT) HBV carriers, 11 cases of chronic hepatitis B (CHB), and 10 cases of acute-on-chronic liver failure (ACLF). The present study was conducted by clone-based sequencing. QS heterogeneity within each open reading frame was calculated. The mutation frequency index (MFI) and amino acid variations within the large HBsAg, HBcAg, and HBxAg regions were analyzed based on the different infection phases. In total, 606 HBV full-length sequences were obtained. HBV QS had higher heterogeneity in ACLF and CHB than that in IT among chronically infected individuals. AHB patients had the lower QS heterogeneity at onset than those with chronic infection. ACLF patients had the highest frequency of mutations in the core promoter and precore region. A triple mutation (A1762T/G1764A/G1896A) was observed more frequently in genotype C than in genotype B. The MFI indicated that specific peptides of the studied regions had more frequent mutations in ACLF. Furthermore, several amino acid variations, known as T- and B-cell epitopes, were potentially associated with the immunoactive phase of infection. More HBV genome mutations and deletions were observed in patients with more severe diseases, particularly in specific regions of the core and preS regions, the clinical significance and mechanism of which need to be further investigated.

Hepatitis C virus genetic variability and evolution

Hepatitis C virus (HCV) has infected over 170 million people worldwide and creates a huge disease burden due to chronic, progressive liver disease. HCV is a single-stranded, positive sense, RNA virus, member of the Flaviviridae family. The high error rate of RNA-dependent RNA polymerase and the pressure exerted by the host immune system, has driven the evolution of HCV into 7 different genotypes and more than 67 subtypes. HCV evolves by means of different mechanisms of genetic variation. On the one hand, its high mutation rates generate the production of a large number of different but closely related viral variants during infection, usually referred to as a quasispecies. The great quasispecies variability of HCV has also therapeutic implications since the continuous generation and selection of resistant or fitter variants within the quasispecies spectrum might allow viruses to escape control by antiviral drugs. On the other hand HCV exploits recombination to ensure its survival. This enormous viral diversity together with some host factors has made it difficult to control viral dispersal. Current treatment options involve pegylated interferon-α and ribavirin as dual therapy or in combination with a direct-acting antiviral drug, depending on the country. Despite all the efforts put into antiviral therapy studies, eradication of the virus or the development of a preventive vaccine has been unsuccessful so far. This review focuses on current available data reported to date on the genetic mechanisms driving the molecular evolution of HCV populations and its relation with the antiviral therapies designed to control HCV infection.

Trans-splicing group I intron targeting hepatitis C virus IRES mediates cell death upon viral infection in Huh7.5 cells

The HCV-IRES sequence is vital for both protein translation and genome replication and serves as a potential target for anti-HCV therapy. We constructed a series of anti-HCV group I introns (αHCV-GrpIs) to attack conserved target sites within the HCV IRES. These αHCV-GrpIs were designed to mediate a trans-splicing reaction that replaces the viral RNA genome downstream of the 5' splice site with a 3' exon that encodes an apoptosis-inducing gene. Pro-active forms of the apoptosis inducing genes BID, Caspase 3, Caspase 8, or tBax were modified by incorporation of the HCV NS5A/5B cleavage sequence in place of their respective endogenous cleavage sites to ensure that only HCV infected cells would undergo apoptosis following splicing and expression. Huh7.5 cells transfected with each intron were challenged at MOI 0.1 with HCV-Jc1FLAG2 which expresses a Gaussia Luciferase (GLuc) marker. Virus-containing supernatants were then assayed for GLuc expression as a measure of viral replication inhibition. Cellular extracts were analyzed for the presence of correct splice products by RT-PCR and DNA sequencing. We also measured levels of Caspase 3 activity as a means of quantifying apoptotic cell death. Each of these αHCV-GrpI introns was able to correctly splice their 3' apoptotic exons onto the virus RNA genome at the targeted Uracil, and resulted in greater than 80% suppression of the GLuc marker. A more pronounced suppression effect was observed with TCID₅₀ virus titrations, which demonstrated that these αHCV-GrpIs were able to suppress viral replication by more than 2 logs, or greater than 99%. Robust activation of the apoptotic factor within the challenged cells was evidenced by a significant increase of Caspase 3 activity upon viral infection compared to non-challenged cells. This novel genetic intervention tool may prove beneficial in certain HCV subjects.

Inhibition of hepatitis C virus by an M1GS ribozyme derived from the catalytic RNA subunit of Escherichia coli RNase P

BACKGROUND

Hepatitis C virus (HCV) is a human pathogen causing chronic liver disease in about 200 million people worldwide. However, HCV resistance to interferon treatment is one of the important clinical implications, suggesting the necessity to seek new therapies. It has already been shown that some forms of the catalytic RNA moiety from E. coli RNase P, M1 RNA, can be introduced into the cytoplasm of mammalian cells for the purpose of carrying out targeted cleavage of mRNA molecules. Our study is to use an engineering M1 RNA (i.e. M1GS) for inhibiting HCV replication and demonstrates the utility of this ribozyme for antiviral applications.

RESULTS

By analyzing the sequence and structure of the 5' untranslated region of HCV RNA, a putative cleavage site (C67-G68) was selected for ribozyme designing. Based on the flanking sequence of this site, a targeting M1GS ribozyme (M1GS-HCV/C67) was constructed by linking a custom guide sequence (GS) to the 3' termini of catalytic RNA subunit (M1 RNA) of RNase P from Escherichia coli through an 88 nt-long bridge sequence. In vitro cleavage assays confirmed that the engineered M1GS ribozyme cleaved the targeted RNA specifically. Moreover, ~85% reduction in the expression levels of HCV proteins and >1000-fold reduction in viral growth were observed in supernatant of cultured cells that transfected the functional ribozyme. In contrast, the HCV core expression and viral growth were not significantly affected by a "disabled" ribozyme (i.e. M1GS-HCV/C67*). Moreover, cholesterol-conjugated M1GS ribozyme (i.e. Chol-M1GS-HCV/C67) showed almost the same bioactivities with M1GS-HCV/C67, demonstrating the potential to improve in vivo pharmacokinetic properties of M1GS-based RNA therapeutics.

CONCLUSION

Our results provide direct evidence that the M1GS ribozyme can function as an antiviral agent and effectively inhibit gene expression and multiplication of HCV.

Clinical implications of evolutionary patterns of homologous, full-length hepatitis B virus quasispecies in different hosts after perinatal infection

Chronic hepatitis B virus (HBV) infection via perinatal transmission is common in the Asia-Pacific region, but related quasispecies (QS) characteristics are not yet defined. To investigate the homologous, full-length HBV QS after perinatal infection and their clinical implications, five pairs of mother-daughter patients with chronic HBV infection (one patient with liver cirrhosis, one with immune tolerance, and eight with chronic hepatitis) were included. Full-length HBV were amplified by PCR from serum samples before antiviral treatment and cloned; an average of 17 clones per sample were sequenced, and the QS complexities, diversities, and evolution patterns were analyzed. Full-length HBV sequence similarities within mother-daughter pairs were 91.3 to 98.3%. The QS complexities of full-length HBV were similar between mothers and daughters (median of 0.9734 compared to 0.9688, P>0.05), as were the diversities (median of 3.396×10(-3) compared to 4.617×10(-3) substitutions/site, P>0.05). However, the distribution patterns of QS complexities and diversities within specific genes were different, and QS genetic distances of the mothers were higher than those of daughters, both more evident in pairs with different antiviral responses and different immune phases or stages. The nucleotide substitution rate of full-length HBV was 14.388×10(-5) substitutions/site/year, whereas the preC/C gene rate was the highest. Mutations and indels were more common in clones from the mothers, which decreased the affinity of epitopes by 6- to 89-fold. The various genes from homologous HBV genomes evolved in different patterns despite numerically comparable full-length QS complexities and diversities. The different patterns may correlate with the immune stages of chronic HBV infection, which warrants further study.

Comparison of next-generation sequencing and clone-based sequencing in analysis of hepatitis B virus reverse transcriptase quasispecies heterogeneity

We previously reported that, based on clone-based sequencing (CBS), hepatitis B virus (HBV) heterogeneity within the reverse transcriptase (RT) region was a predictor of antiviral efficacy. Here, by comparing ultradeep pyrosequencing (UDPS), i.e., next-generation sequencing (NGS), with CBS in characterizing the genetic heterogeneity of HBV quasispecies within the RT region, we evaluated the performance of UDPS in the analysis of HBV viral populations. HBV genomic DNA was extracted from serum samples from 31 antiviral treatment-naive patients with chronic hepatitis B. The RT region quasispecies were analyzed in parallel using CBS and UDPS. Characterization of quasispecies heterogeneity was conducted using bioinformatics analysis. Quasispecies complexity values were calculated with the formula Sn = -Σi(pilnpi)/lnN. The number of qualified strains obtained by UDPS was much larger than that obtained by CBS (P < 0.001). Pearson analysis showed that there was a positive correlation of quasispecies complexity values at the nucleotide level for the two methods (P < 0.05), while the complexity value derived from UDPS data was higher than that derived from CBS data (P < 0.001). Study of the prevalences of variations within the RT region showed that CBS detected an average of 9.7 ± 1.1 amino acid substitutions/sample and UDPS detected an average of 16.2 ± 1.4 amino acid substitutions/sample. The phylogenetic analysis based on UDPS data showed more genetic entities than did that based on CBS data. Viral heterogeneity determination by the UDPS technique is more sensitive and efficient in terms of low-abundance variation detection and quasispecies simulation than that by the CBS method, although imperfect, and thus sheds light on the future clinical application of NGS in HBV quasispecies studies.

Multi-agent model of hepatitis C virus infection

OBJECTIVES

The objective of this study is to design a method for modeling hepatitis C virus (HCV) infection using multi-agent simulation and to verify it in practice.

METHODS AND MATERIALS

In this paper, first, the modeling of HCV infection using a multi-agent system is compared with the most commonly used model type, which is based on differential equations. Then, the implementation and results of the model using a multi-agent simulation is presented. To find the values of the parameters used in the model, a method using inverted simulation flow and genetic algorithm is proposed. All of the data regarding HCV infection are taken from the paper describing the model based on the differential equation to which the proposed method is compared.

RESULTS

Important advantages of the proposed method are noted and demonstrated: these include flexibility, clarity, re-usability and the possibility to model more complex dependencies. Then, the simulation framework that uses the proposed approach is successfully implemented in C++ and is verified by comparing it to the approach based on differential equations. The verification proves that an objective function that performs the best is the function that minimizes the maximal differences in the data. Finally, an analysis of one of the already known models is performed, and it is proved that it incorrectly models a decay in the hepatocytes number by 40%.

CONCLUSIONS

The proposed method has many advantages in comparison to the currently used model types and can be used successfully for analyzing HCV infection. With almost no modifications, it can also be used for other types of viral infections.

Phylogeny and molecular evolution of the hepatitis C virus

The hepatitis C virus (HCV) is a globally prevalent human pathogen that causes persistent liver infections in most infected individuals. HCV is classified into seven phylogenetically distinct genotypes, which have different geographical distributions and levels of genetic diversity. Some of these genotypes are endemic and highly divergent, whereas others disseminate rapidly on an epidemic scale but display lower variability. HCV phylogeny has an important impact on disease epidemiology and clinical practice because the viral genotype may determine the pathogenesis and severity of the resultant chronic liver disease. In addition, there is a clear association between the HCV genotype and its susceptibility to antiviral treatment. Similarly to other RNA viruses, in a single host, HCV exists as a combination of related but genetically different variants. The whole formation is the actual target of selection exerted by a host organism and antiviral therapeutics. The genetic structure of the viral population is largely shaped by mutations that are constantly introduced during an error-prone replication. However, it appears that genetic recombination may also contribute to this process. This heterogeneous collection of variants has a significant ability to evolve towards the fitness optimum. Interestingly, negative selection, which restricts diversity, emerges as an essential force that drives HCV evolution. It is becoming clear that HCV evolves to become stably adapted to the host environment. In this article we review the HCV phylogeny and molecular evolution in the context of host-virus interactions.

Distinct quasispecies characteristics and positive selection within precore/core gene in hepatitis B virus HBV associated acute-on-chronic liver failure

BACKGROUND AND AIM

The cause of hepatitis B virus associated acute-on-chronic liver failure (ACLF) remains unclear. Quasispecies can contribute to virus persistence and pathogenesis. We used a bioinformatics-based molecular evolution approach to compare quasispecies characteristics and positive selection sites within HBV precore/core gene between ACLF and chronic hepatitis B (CHB) patients.

METHODS

HBV precore/core gene were amplified from 11 ACLF and 10 CHB patients harboring HBV genotype B; following DNA cloning and sequencing quasispecies complexity, diversity, and positive selection sites within the precore/core gene were determined by bioinformatics analysis, and compared between the patient groups.

RESULTS

Both quasispecies complexity (P=0.022 at nucleotide level and 0.008 at amino acid level) and diversity (P<0.05) were found to be significantly greater in ACLF than in CHB. The frequency of G1896/A mutation in ACLF (175/298 clones, 58.7%) was also significantly higher than in CHB (100/230 clones, 43.5%) (P=0.0005). Moreover, analysis of positive selection revealed that significantly more patients with such sites were present in ACLF than in CHB (8/11 VS 2/10, P=0.03); the majority of these positive selection sites lay within HLA-restricted epitopes.

CONCLUSIONS

The ACLF patients showed distinct quasispecies characteristics with higher complexity and diversity within the HBV precore/core gene. The increased HBV quasispecies complexity and diversity, together with a distinct set of positive selection sites, is likely associated with the development of ACLF.