Quasispecies structure,cornerstone of hepatitis B virus infection: Mass sequencing approach

Hepatitis B virus(HBV)is a DNA virus with complex replication,and high replication and mutation rates,leading to a heterogeneous viral population.The population is comprised of genomes that are closely related,but not identical;hence,HBV is considered a viral quasispecies.Quasispecies variability may be somewhat limited by the high degree of overlapping between the HBV coding regions,which is especially important in the P and S gene overlapping regions,but is less significant in the X and preCore/Core genes.Despite this restriction,several clinically and pathologically relevant variants have been characterized along the viral genome.Next-generation sequencing(NGS)approaches enable high-throughput analysis of thousands of clonally amplified regions and are powerful tools for characterizing genetic diversity in viral strains.In the present review,we update the information regarding HBV variability and present a summary of the various NGS approaches available for research in this virus.In addition,we provide an analysis of the clinical implications of HBV variants and their study by NGS.

Characterization of hepatitis B virus X gene quasispecies complexity in mono-infection and hepatitis delta virus superinfection

Hepatitis delta virus(HDV) seems to strongly suppress hepatitis B virus(HBV)replication, although little is known about the mechanism of this interaction. Both these viruses show a dynamic distribution of mutants, resulting in viral quasispecies. Next-generation sequencing is a viable approach for analyzing the composition of these mutant spectra. As the regulatory hepatitis B X protein(HBx) is essential for HBV replication, determination of HBV X gene(HBX)quasispecies complexity in HBV/HDV infection compared to HBV monoinfection may provide information on the interactions between these two viruses.AIM To compare HBV quasispecies complexity in the HBX 5' region between chronic hepatitis delta(CHD) and chronic HBV mono-infected patients.METHODS Twenty-four untreated patients were included: 7/24(29.2%) with HBeAgnegative chronic HBV infection(CI, previously termed inactive carriers), 8/24(33.3%) with HBeAg-negative chronic hepatitis B(CHB) and 9/24(37.5%) with CHD. A serum sample from each patient was first tested for HBV DNA levels.The HBX 5' region [nucleotides(nt) 1255-1611] was then PCR-amplified for subsequent next-generation sequencing(MiSeq, Illumina, United States). HBV quasispecies complexity in the region analyzed was evaluated using incidencebased indices(number of haplotypes and number of mutations), abundancebased indices(Hill numbers of order 1 and 2), and functional indices(mutation frequency and nucleotide diversity). We also evaluated the pattern of nucleotide changes to investigate which of them could be the cause of the quasispecies complexity.RESULTS CHB patients showed higher median HBV-DNA levels [5.4 logIU/mL,interquartile range(IQR) 3.5-7.9] than CHD(3.4 logIU/mL, IQR 3-7.6)(P = n.s.)or CI(3.2 logIU/mL, IQR 2.3-3.5)(P < 0.01) patients. The incidence and abundance indices indicated that HBV quasispecies complexity was significantly greater in CI than CHB. A similar trend was observed in CHD patients, although only Hill numbers of order 2 showed statistically significant differences(CHB2.81, IQR 1.11-4.57 vs CHD 8.87, 6.56-11.18, P = 0.038). There were no significant differences in the functional indices, but CI and CHD patients also showed a trend towards greater complexity than CHB. No differences were found for any HBV quasispecies complexity indices between CHD and CI patients. G-to-A and C-to-T nucleotide changes, characteristic of APOBEC3 G, were higher in CHD and CI than in CHB in genotype A haplotypes, but not in genotype D. The proportion of nt G-to-A vs A-to-G changes and C-to-T vs T-to-C changes in genotype A and D haplotypes in CHD patients showed no significant differences. In CHB and CI the results of these comparisons were dependent on HBV genotype.CONCLUSION The lower-replication CHD and CI groups show a trend to higher quasispecies complexity than the higher-replication CHB group. The mechanisms associated with this greater complexity require elucidation.

Cross-sectional evaluation of circulating hepatitis B virus RNA and DNA: Different quasispecies?

BACKGROUND Different forms of pregenomic and other hepatitis B virus(HBV)RNA have been detected in patients’sera.These circulating HBV-RNAs may be useful for monitoring covalently closed circular DNA activity,and predicting hepatitis B eantigen seroconversion or viral rebound after nucleos(t)ide analog cessation.Data on serum HBV-RNA quasispecies,however,is scarce.It is therefore important to develop methodologies to thoroughly analyze this quasispecies,ensuring the elimination of any residual HBV-DNA.Studying circulating HBV-RNA quasispecies may facilitate achieving functional cure of HBV infection.AIM To establish a next-generation sequencing(NGS)methodology for analyzing serum HBV-RNA and comparing it with DNA quasispecies.METHODS Thirteen untreated chronic hepatitis B patients,showing different HBV-genotypes and degrees of severity of liver disease were enrolled in the study and a serum sample with HBV-DNA>5 Log10 IU/mL and HBV-RNA>4 Log10 copies/mL was taken from each patient.HBV-RNA was treated with DNAse I to remove any residual DNA,and the region between nucleotides(nt)1255-1611 was amplified using a 3-nested polymerase chain reaction protocol,and analyzed with NGS.Variability/conservation and complexity was compared between HBV-DNA and RNA quasispecies.RESULTS No HBV-DNA contamination was detected in cDNA samples from HBV-RNA quasispecies.HBV quasispecies complexity showed heterogeneous behavior among patients.The Rare Haplotype Load at 1%was greater in DNA than in RNA quasispecies,with no statistically significant differences(P=0.1641).Regarding conservation,information content was equal in RNA and DNA quasispecies in most nt positions[218/357(61.06%)].In 102 of the remaining 139(73.38%),HBV-RNA showed slightly higher variability.Sliding window analysis identified 4 hyper-conserved sequence fragments in each quasispecies,3 of them coincided between the 2 quasispecies:nts 1258-1286,1545-1573 and 1575-1604.The 2 hyper-variable sequence fragments also coincided:nts 1311-1344 and 1461-1485.Sequences between nts 1519-1543 and 1559-1587 were only hyper-conserved in HBV-DNA and RNA,respectively.CONCLUSION Our methodology allowed analyzing HBV-RNA quasispecies complexity and conservation without interference from HBV-DNA.Thanks to this,we have been able to compare both quasispecies in the present study.

Quasispecies groups in the core promoter region of hepatitis B virus

Objectives: To investigate the mutation of the basic core promoter (BCP) of hepatitis B virus (HBV) and clarify the significance of HBV quasispecies groups in patients with chronic HBV infection. Methods: A set of specific primers was synthesized according to the HBV DNA sequence of a Chinese strain. The BCP was amplified by PCR method from the serum of 40 patients with chronic HBV infection, and the PCR products of 2 patients were subcloned into pGEM Teasy vectors. Polyacrylamide gel elec- trophoresis (PAGE) was employed to display the de- letion mutations, and clones with differential length were selected to be sequenced. Sequence comparison was made to find the difference. Results: Two or three bands were displayed by PAGE in 60% patients. The results of sequence anal- ysis showed that there are some kinds of mutations in the BCP region. The substitution always occurs in TATA-like boxes, especially from T to C on 140 site. The deletion mutations were detected in TA1, TA2 and TA3. The 8bp, 20bp deletion mutations fre- quently happened. Conclusions: There is a hot deletion region in the BCP. The deletion and the substitution in the TATA- like box may influence the expression of preC/C pro- tein. The sequencing results indicate that there are HBV quasispecies groups in patients with chronic HBV infection.

Early Evolution of Hepatitis B Virus Quasispecies During IFN-α Treatment

Objective To investigate the dynamic change of hepatitis B virus quasispecies within complete genome during the early stage of IFN-α treatment and its impact on virological response.Methods Sixteen patients with chronic hepatitis B receiving IFN-α treatment were investigated. HBV DNA was extracted from serum sample at baseline and week 12. The complete genome of HBV was amplified, then cloned and sequenced. The quasispecies heterogeneity of HBV complete genome was depicted at baseline and week 12. Results The quasispecies heterogeneity of the genome except for C-ORF were comparable in three groups at baseline and week 12. The quasispecies diversity at amino acid levels of responders within C-ORF were higher than that of non-responders at baseline. The quasispecies diversity within the C-ORF of partial responders was reduced in the early stage of IFN-α treatment. Furthermore, the mean genetic distance at amino acid levels of partial responders was significantly higher than that of the non-responders at week 12. The evolutionary rate was not different between non-responders and partial responders. Conclusions In the immune clearance phase, the patients who had greater viral quasispecies diversity within C-ORF at amino acid level had more chance to obtain the early virological response during IFN-α treatment.

Rethinking quasispecies theory: From fittest type to cooperative consortia

Recent investigations surprisingly indicate that single RNA "stem-loops" operate solely by chemical laws that act without selective forces, and in contrast, self-ligated consortia of RNA stem-loops operate by biological selection. To understand consortial RNA selection, the concept of single quasi-species and its mutant spectra as drivers of RNA variation and evolution is rethought here. Instead, we evaluate the current RNA world scenario in which consortia of cooperating RNA stem-loops(not individuals) are the basic players. We thus redefine quasispecies as RNA quasispecies consortia(qs-c) and argue that it has essential behavioral motifs that are relevant to the inherent variation, evolution and diversity in biology. We propose that qs-c is an especially innovative force. We apply qs-c thinking to RNA stem-loops and evaluate how it yields altered bulges and loops in the stem-loop regions, not as errors, but as a natural capability to generate diversity. This basic competencenot error-opens a variety of combinatorial possibilities which may alter and create new biological interactions, identities and newly emerged self identity(immunity) functions. Thus RNA stem-loops typically operate as cooperative modules, like members of social groups. Fromsuch qs-c of stem-loop groups we can trace a variety of RNA secondary structures such as ribozymes, viroids, viruses, mobile genetic elements as abundant infection derived agents that provide the stem-loop societies of small and long non-coding RNAs.

Quasispecies of Hepatitis B Virus

Hepatitis B virus(HBV) circulates in blood and replicates in the presence of quasispecies. During HBV replication, HBV DNA polymerase lacks fidelity and proofreading function partly because its exonuclease activity is either absent or deficient. Therefore, HBV genome is mutated with unusually high frequency. And these mutations can affect more than one open reading frame due to overlapping genes. Otherwise, natural substitutions, deletions or insertions involving the Cp/ENⅡ locus in the X gene can significantly alter the extent of viral replication activity. Particular selection pressures such as host immune system and antiviral therapy readily select out escape mutants from this pre-existing quasispecies pool. Antiviral drug resistance in chronic hepatitis B(CHB) can be caused by the viral mutation frequency, the intrinsic mutability of the antiviral target site, the selective pressure exerted by the drug, the magnitude and rate of virus replication, the overall replication fitness of the mutant, the genetic barrier of the compound and the availability of replication space. Potent inhibition of HBV replication could be able to prevent the development of drug resistance because mutagenesis is replication dependent. Viral load may decline to a point where the continued production of quasispecies with the potential to resist new drug treatments no longer occurs, if viral replication can be suppressed for a sufficient length of time.

Variations and mutations in the hepatitis B virus genome and their associations with clinical characteristics

Hepatitis B virus(HBV) infection is major global issue, because chronic HBV infection is strongly associated with liver cancer. HBV spread worldwide with variousmutations and variations. This variability, called quasispecies, is derived from no proof-reading capacity of viral reverse transcriptase. So far, thousands of studies reported that the variety of genome is closely related to the geographic distribution and clinical characteristics. Recent technological advances including capillary sequencer and next generation sequencer have made in easier to analyze mutations. The variety of HBV genome is related to not only antigenicity of HBs-antigen but also resistance to antiviral therapies. Understanding of these variations is important for the development of diagnostic tools and the appropriate therapy for chronic hepatitis B. In this review, recent publications in relation to HBV mutations and variations are updated and summarized.

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 singlestranded, 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.

Clinical relevance of hepatitis B virus variants

The hepatitis B virus(HBV) is a global public health problem with more than 240 million people chronically infected worldwide, who are at risk for end-stage liver disease and hepatocellular carcinoma. There are an estimated 600000 deaths annually from complications of HBV-related liver disease. Antiviral therapy with nucleos/tide analogs(NA) targeting the HBV polymerase(P) can inhibit disease progression by long-term suppression of HBV replication. However, treatment may fail with first generation NA therapy due to the emergence of drugresistant mutants, as well as incomplete medication adherence. The HBV replicates via an error-prone reverse transcriptase leading to quasispecies. Due to overlapping open reading frames mutations within the HBV P can cause concomitant changes in the HBV surface gene(S) and vice versa. HBV quasispecies diversity is associated with response to antiviral therapy, disease severity and long-term clinical outcomes. Specific mutants have been associated with antiviral drug resistance, immune escape, liver fibrosis development and tumorgenesis. An understanding of HBV variants and their clinical relevance may be important for monitoring chronic hepatitis B disease progression and treatment response. In this review, we will discuss HBV molecular virology, mechanism of variant development, and their potential clinical impact.

Polymorphisms and resistance mutations of hepatitis C virus on sequences in the European hepatitis C virus database

AIM To evaluate the occurrence of resistant mutations in treatment-na?ve hepatitis C virus(HCV) sequences deposited in the European hepatitis C virus database(euH CVdb). METHODS The sequences were downloaded from the eu HCVdb(https://euhcvdb.ibcp.fr/eu HCVdb/). The search was performed for full-length NS3 protease, NS5 A and NS5 B polymerase sequences of HCV, separated by genotypes 1a, 1b, 2a, 2b and 3a, and resulted in 798 NS3, 708 NS5 A and 535 NS5 B sequences from HCV genotypes1a, 1b, 2a, 2b and 3a, after the exclusion of sequences containing errors and/or gaps or incomplete sequences, and sequences from patients previously treated with direct antiviral agents(DAA). The sequence alignment was performed with MEGA 6.06 MAC and the resulting protein sequences were then analyzed using the BioE dit 7.2.5. for mutations associated with resistance. Only positions that have been described as being associated with failure in treatment in in vivo studies, and/or as conferring a more than 2-fold change in replication in comparison to the wildtype reference strain in in vitro phenotypic assays were included in the analysis.RESULTS The Q80 K variant in the NS3 gene was the most prevalent mutation, being found in 44.66% of subtype 1a and 0.25% of subtype 1b. Other frequent mutations observed in more than 2% of the NS3 sequences were: I170V(3.21%) in genotype 1a, and Y56F(15.93%), V132I(23.28%) and I170V(65.20%) in genotype 1b. For the NS5 A, 2.21% of the genotype 1a sequences have the P58 S mutation, 5.95% of genotype 1b sequences have the R30 Q mutation, 15.79% of subtypes 2a sequences have the Q30 R mutation, 23.08% of subtype 2b sequences have a L31 M mutation, and in subtype 3a sequences, 23.08% have the M31 L resistant variants. For the NS5 B, the V321 L RAV was identified in 0.60% of genotype 1a and in 0.32% of genotype 1b sequences, and the N142 T variant was observed in 0.32% of subtype 1b sequences. The C316 Y, S556 G, D559 N RAV were identified in 0.33%, 7.82% and 0.32% of genotype 1b sequences, respectively, and were not observed in other genotypes.CONCLUSION HCV mutants resistant to DAAs are found in low frequency, nevertheless they could be selected and therapy could fail due resistance substitutions in HCV genome.

Bayesian Coalescent Analysis of the Intra-Host Evolution of Hepatitis C Virus: Memory Genomes and Clinical Implications

Genetic variability plays a key role in the biology and medical treatment of RNA viruses. As an RNA virus, Hepatitis C virus (HCV) replicate as complex distributions of closely related genomes termed viral quasispecies. The behavior of the evolving HCV quasispecies population is influenced by the ensemble of mutants that compose the viral population. One such influence is the presence of minority subpopulations, termed memory genomes, in the mutant spectra. Biologically relevant mutants have been previously observed to be present as memory genomes in RNA viral populations. For that reason, an in-depth analysis of HCV quasispecies populations is crucial for our understanding viral evolution, drug resistance and therapy outcome. Recently developed next-generation sequencing (NGS) platforms make it possible to investigate viral quasispecies at much greater detail. In order to gain insight into these matters, we have performed a Bayesian coalescent analysis of hypervariable region 1 (HVR1) sequences of a HCV quasispecies population circulating in a chronic patient, recently obtained by ultra-deep sequencing. The results of these studies revealed a mean rate of evolution of HCV HVR1 of the intra-host quasispecies population of 4.80 × 10-2 amino acid substitutions/site/year. A sharp and rapid diversification of the HCV quasispecies isolated from the patient in three different sub-populations was observed. The most abundant sequence in the quasispecies population was not found to be the center of a tight and complex network around this sequence, suggesting that the quasispecies population as a whole efficiently explore a wide sequence space. Co-evolution of relevant amino acid sites had been identified in the HVR1. This speaks of the possible roll of these residues in HVR1 to allow the virus to shift between combinations of residues to escape the immune system while retaining its structure and functions. The results of these studies highlight the importance of minority genomes in HCV population history and evolution, the mutant clouds as reservoirs of phenotypic and genetic variants for virus adaptability, as well as the roll of the mutant spectra to overcome selective constraints.

Molecular Mechanisms of Raccoon Rabies Virus Progression in Its Natural Host

Rabies virus presents a global public health problem. Our current understanding of the molecular determinants of rabies virulence stems from rodent models and laboratory strains of the virus, however, it is unclear how well rodent models represent viral response in natural reservoirs. Here, we examined interactions between the raccoon variant of rabies virus (RRV) and its natural host, raccoons, to gain a better understanding of molecular determinants of virulence in this system. We found expression patterns of RRV genes under tight control until the virus reached the central nervous system where replication increased significantly. Further, our examination of viral variants within an individual revealed that variant diversity may have an effect on virulence. We found that a mutation at a region of a T helper cell epitope on the nucleoprotein was associated with viral challenge outcomes and could be associated with RRV pathogenicity.

Estimate of parameters changing of the Lymantria dispar artificial population in different environmental conditions using the quasispecies model

This paper gives an interpretation of the quasispecies concept for an experiment of growing the Lymantria dispar artificial population.The changes estimate in the population parameters for a significant change in environmental conditions at the stage of hatching caterpillars from the eggs was given.Three numerical and two random parameters were used for quantitative analysis.They characterize the development of individuals of the population at this stage.The parameter changes significance was verified with using statistical tests.The results obtained do not contradict the quasispecies model.The use of fuzzy sets describes the transition from one quasispecies to another at a significant change in environmental conditions.

Mathematical formalization of the process of crossing individuals of different populations within the framework of the quasispecies model

This paper presents a mathematically formalized description of the process of crossing individuals of different populations as a break of the first kind of functional dependencies of the population state vector components.In this case,the validity of the quasispecies model as a model of population development was assumed.The practical interpretation of a mathematically formalized description is considered on the example of an artificial population of Lymantria dispar.The results of the work can be used to assess the development parameters of both natural and artificial populations.

Criticality, adaptability and early-warning signals in time series in a discrete quasispecies model

Complex systems from different fields of knowledge often do not allow a mathematical description or modeling, because of their intricate structure composed of numerous interacting components. As an alternative approach, it is possible to study the way in which observables associated with the system fluctuate in time. These time series may provide valuable information about the underlying dynamics. It has been suggested that complex dynamic systems, ranging from ecosystems to financial markets and the climate, produce generic early-warning signals at the ＂tipping points,＂ where they announce a sudden shift toward a different dynamical regime, such as a population extinction, a systemic market crash, or abrupt shifts in the weather. On the other hand, the framework of Self- Organized Criticality （SOC）, suggests that some complex systems, such as life itself, may spontaneously converge toward a critical point. As a particular example, the quasispecies model suggests that RNA viruses self-organize their mutation rate near the error-catastrophe threshold, where robustness and evolvability are balanced in such a way that survival is optimized. In this paper, we study the time series associated to a classical discrete quasispecies model for different mutation rates, and identify early-warning signals for critical mutation rates near the error-catastrophe threshold, such as irregularities in the kurtosis and a significant increase in the autocorrelation range, reminiscent of 1/f noise. In the present context, we find that the early-warning signals, rather than broadcasting the collapse of the system, are the fingerprint of survival optimization.

Evolutionary study of hepatitis C virus envelope genes during primary infection

Background Hepatitis C virus （HCV） envelope genes encoding glycoproteins E1 and E2 exhibits a high degree of variability that gives rise to differing phenotypic traits; including alterations in receptor-binding affinity and immune recognition and escape. This study aims to elucidate the relationship of the evolutionary patterns for HCV envelope glycoproteins to viral persistence. Methods HCV quasispecies were characterized in specimens collected every two to six months from a cohort of acutely HCV-infected subjects. We evaluated two individuals who spontaneously cleared viremia and three individuals with persistent viremia by cloning 33 1-kb amplicons that spanned E1 and the 5＇ half of E2; including hypervariable region 1 （HVR1）. To detect representative variants for sequencing thirty-three cloned cDNAs representing each specimen were assessed by a method that combined analysis of a single-stranded conformational polymorphism （SSCP） method and heteroduplex analysis （HDA）. For each patient, the rates of both synonymous and nonsynonymous substitutions for the El, HVR1 and E2 regions outside HVR1 were evaluated. The amino acid sequences and predicted antigenic profiles were analyzed. Results The genetic diversity within HVR1 was consistently higher than that in the E1 and E2 regions outside HVR1 in individuals with persistent viremia, but did not change markedly over time in those with clearance of viremia. For individuals with persistent viremia, the rate of nonsynonymous substitutions within the HVR1 region predominated and gradually increased, compared to that in the E1 and E2 regions outside HVR1. By contrast, the rates of both nonsynonymous and synonymous substitutions for the E1 and E2 regions, including HVR1, were consistently lower in individuals with clearance of viremia. HVR1 had a higher antigenic variable and lower positive charge in subjects with persistent viremia. All cysteine residues and N-linked glycosylation sites, some of which were known to play a major role in protein folding and others play a role in HCV entry, were 100% conserved among the sequenced cloned cDNAs from the two outcome groups. Conclusion HCV persistence may be associated with positive selection pressures on HVR1, rather than functional constraints in the envelope reaion.

The effect analysis of the diapause duration on the hatching parameters of the artificial semi-annual population of Lymantria dispar caterpillars

The paper investigates the existence of a relationship between the diapause duration and the temperature threshold of hatching,the variances in temperature and hatching time of the artificial six-month population of Lymantria dispar caterpillars.A half-year population gives two generations a year under artificial conditions.A comparative analysis of these parameters with similar parameters of the hybrid population of L.dispar,as well as studies by other authors,has been carried out.The Eigen quasispecies model was used for the mathematical formalization of the evolution process.The results of this work can be used to study the evolution of L.dispar.