Gene Quasispecies Research Articles

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

BACKGROUNDHepatitis 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 mono-infection may provide information on the interactions between these two viruses.AIMTo compare HBV quasispecies complexity in the HBX 5’ region between chronic hepatitis delta (CHD) and chronic HBV mono-infected patients.METHODSTwenty-four untreated patients were included: 7/24 (29.2%) with HBeAg-negative 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 incidence-based indices (number of haplotypes and number of mutations), abundance-based 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.RESULTSCHB 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 (CHB 2.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 APOBEC3G, 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.CONCLUSIONThe 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.

Analysis of HBV X gene quasispecies characteristics by next-generation sequencing and cloning-based sequencing and its association with hepatocellular carcinoma progression.

This study aimed to describe the differences between next-generation sequencing (NGS) and cloning-based sequencing (CBS) in HBX quasispecies research and primitively investigate the relationship between the dominant HBX quasispecies and hepatocellular carcinoma (HCC). A total of 12 serum samples were collected. Serum hepatitis B virus (HBV) DNA was extracted, and the HBV X-region (HBX) was amplified by nested polymerase chain reaction (PCR). The PCR products were simultaneously tested with NGS and CBS to detect quasispecies of the HBX. A total of 9348 eligible quasispecies sequences were obtained by NGS, which were much larger than the 98 of that by CBS. By the phylogenetic tree, the dominant quasispecies sequence of each sample could be found, although they had several nucleotides differences between the dominant quasispecies sequences found by CBS and NGS. By comparing the quasispecies heterogeneity, it was found that the quasispecies complexity value of HBV X-region obtained by NGS was higher than CBS (P < 0.05). The diversity values, including d, dS, dN, an d d N/ dS obtained by NGS were lower than by CBS (all of P < 0.01). The relativity of Spearman(rs) in d, dS, and dN were statistically significant (rs_ d = 0.865, P = 0.001; rs_ dS = 0.722, P = 0.014; and rs_ dN = 0.738, P = 0.011, respectively). There were 21 different bases between the HBX quasispecies of case A and control B. The results of this can be used as guidance when researchers plan to choose a suitable method to study quasispecies, especially the HBV X gene quasispecies. Some high-risk mutations of HBX quasispecies were also found in this study and their relationship with HCC need deeper exploration.

Evolution of H5 highly pathogenic avian influenza: sequence data indicate stepwise changes in the cleavage site.

The genetic composition of an H5 subtype hemagglutinin gene quasispecies, obtained from ostrich tissues that had been infected with H5 subtype influenza virus was analysed using a next generation sequencing approach. The first evidence for the reiterative copying of a poly (U) stretch in the connecting peptide region in the haemagglutinin cleavage site (HACS) by the viral RNA-dependent RNA polymerase (RdRp) is provided. Multiple non-consensus species of RNA were detected in the infected host, corresponding to likely intermediate sequences between the putative low pathogenic precursor nucleotide sequence of the H5 influenza strain and the highly pathogenic avian influenza virus gene sequence. In silico analysis of the identified RNA sequences predicted that the intermediary H5 sequence PQREKRGLF plays an important role in subsequent mutational events that relocate the HACS coding region from stable base-paired RNA regions to a single-stranded bulge, thereby priming the connecting peptide coding region for RdRp slippage.

Impact of HCV core gene quasispecies on hepatocellular carcinoma risk among HALT-C trial patients.

Mutations at positions 70 and/or 91 in the core protein of genotype-1b, hepatitis C virus (HCV) are associated with hepatocellular carcinoma (HCC) risk in Asian patients. To evaluate this in a US population, the relationship between the percentage of 70 and/or 91 mutant HCV quasispecies in baseline serum samples of chronic HCV patients from the HALT-C trial and the incidence of HCC was determined by deep sequencing. Quasispecies percentage cut-points, ≥42% of non-arginine at 70 (non-R70) or ≥98.5% of non-leucine at 91 (non-L91) had optimal sensitivity at discerning higher or lower HCC risk. In baseline samples, 88.5% of chronic HCV patients who later developed HCC and 68.8% of matched HCC-free control patients had ≥42% non-R70 quasispecies (P = 0.06). Furthermore, 30.8% of patients who developed HCC and 54.7% of matched HCC-free patients had quasispecies with ≥98.5% non-L91 (P = 0.06). By Kaplan-Meier analysis, HCC incidence was higher, but not statistically significant, among patients with quasispecies ≥42% non-R70 (P = 0.08), while HCC incidence was significantly reduced among patients with quasispecies ≥98.5% non-L91 (P = 0.01). In a Cox regression model, non-R70 ≥42% was associated with increased HCC risk. This study of US patients indicates the potential utility of HCV quasispecies analysis as a non-invasive biomarker of HCC risk.

Co-evolution of compensatory mutation K43E with mutation M41L in long-term HIV antiretroviral treatment.

Compensatory mutations have been observed to emerge with drug resistance (DR) mutations, but their effects on virological response to treatment have not been fully examined. In this study, we characterized the emergence and depletion dynamics of a compensatory mutation K43E that correlated with primary nucleoside reverse transcriptase inhibitor (NRTI) drug resistance mutations in Chinese HIV patients on antiretroviral treatment. Single Genome Amplification (SGA) was used to obtain the HIV-1 pol gene quasispecies in three patients over 6 years of Antiretroviral Therapy (ART) treatment. SGA sequences were analyzed by Markov Chain Monte Carlo (MCMC) phylogenetic trees with molecular clock to identify and track compensatory mutation K43E correlated with primary DR mutation M41L. We evaluated the relationship between potential compensatory mutation and DR mutations through Ka/Ks ratio, Jaccard similarity coefficient, and compared these with concurrent viral load data. We determined that a known compensatory mutation, K43E, frequently co-occurs with the drug resistance mutation M41L and may offer a significant advantage in the long-term survival of such drug resistant strains.

Detection of haemagglutinin D222 polymorphisms in influenza A(H1N1)pdm09-infected patients by ultra-deep pyrosequencing

This study was aimed at establishing the genetic heterogeneity of influenza virus haemagglutinin (HA) gene quasi-species and the polymorphisms at codon 222, by application of ultra-deep pyrosequencing (UDPS) to respiratory samples from patients hospitalized for influenza A(H1N1)pdm09 infection, presenting with severe or moderate–mild disease. HA diversity was significantly higher in samples collected from patients with severe manifestations than in those from patients with moderate–mild manifestations (p 0.02). D222 polymorphism was detected in 40.7% of patients by UDPS, and in only 7.1% by Sanger sequencing. D222E, D222G, D222N and D222A were observed in 37.0%, 11.1%, 7.4% and 3.7% of patients, respectively; 10.7% of samples harboured more than two variants. The relative frequency of each single variant showed a wide range of intrapatient variation. D222G/N/A were detected, as either minor or predominant variants, only in severe cases, whereas D222E was equally represented in severe and moderate–mild infections. Other amino acid variants were observed at different positions within the analysed HA fragment. Consistent with higher heterogeneity, non-D222 variants were more frequently detected in severe cases than in moderate–mild cases. In addition, seven non-D222 mutations carried by minority variants, not previously described, were observed.

Complexity and Catalytic Efficiency of Hepatitis C Virus (HCV) NS3 and NS4A Protease Quasispecies Influence Responsiveness to Treatment with Pegylated Interferon plus Ribavirin in HCV/HIV-Coinfected Patients

The role of the hepatitis C virus (HCV) NS3/4A protease in ablating the signaling pathway involved in the production of alpha/beta interferon (IFN-α/β) suggests a relationship between NS3/4A proteolytic activity and a patient's response to IFN-based therapy. To identify viral factors associated with the HCV treatment response, we analyzed the pretreatment NS3/4A protease gene quasispecies composition of 56 HCV genotype 1-HIV-1-coinfected patients treated in our clinic with pegylated IFN (pegIFN) plus ribavirin (RBV). The catalytic efficiency of the dominant (i.e., the most abundant) quasispecies was also assayed for Cardif cleavage and correlated with treatment outcome. A total of 1,745 clones were isolated and sequenced. Significantly less nucleotide quasispecies heterogeneity and lower Shannon entropy values were detected within the responder group (P < 0.05). A correlation was also found between the efficiency of NS3/4A protease Cardif cleavage and therapy outcome. Proteases from sustained responder patients were more efficient at processing Cardif (mean ± standard error of the mean [SEM], 0.8960 ± 0.05568; n = 19) than proteases from nonresponders (mean ± SEM, 0.7269 ± 0.05306; n = 37; P < 0.05). Finally, the amino acid p distance (the proportion [p] of nucleotide sites at which two sequences being compared are different) was significantly shorter in patients with an interleukin-28B (IL-28B) risk allele (P < 0.01), suggesting that IL-28B risk allele carriers exert a lower positive selection pressure on the NS3/4A protease. NS3/4A protease efficiency in cleaving Cardif may be associated with the pegIFN-RBV treatment response, as shown in our cohort of HIV-HCV-coinfected patients. Greater NS3/4A nucleotide heterogeneity and higher Shannon entropy values in nonresponders suggest that less HCV quasispecies complexity may favor a better response to pegIFN-RBV.

The 30th anniversary of quasispecies

The meeting on ‘Quasispecies: past, present and future’ took place between 17 and 18 November 2008, in Barcelona, Spain, and was organized by J. Gomez, C. Lopez‐Galindez, M.A. Martinez & A. Mas. ![][1] A meeting was held in Barcelona, Spain, in November 2008 to celebrate the 30th anniversary of the publication of the article that described the extensive genetic heterogeneity of bacteriophage Qβ (Domingo et al , 1978), which is considered to mark the beginning of experimental studies on viral quasispecies. This meeting was held at the impressive fifteenth century building of the ancient Hospital de la Santa Creu in the old town of Barcelona, which is now the headquarters of the Institut d'Estudis Catalans, and was attended by C. Weissmann (Jupiter, FL, USA), M. Billeter (Zurich, Switzerland) and E. Domingo (Madrid, Spain), who were three early protagonists of the phage Qβ work at the University of Zurich in the 1970s. Several speakers presented their results on the theoretical aspects of the population dynamics of cells and viruses, the clinical implications of quasispecies, and extensions of the quasispecies concept to cellular genes and prions. The meeting was introduced by A. Mas (Albacete, Spain), who reflected on the increasing impact that quasispecies have had in the scientific literature over the past three decades, and quoted some of the key references on viral quasispecies (Martell et al , 1992; Meyerhans et al , 1989; Najera et al , 1995; Vignuzzi et al , 2006; for a historical review of the impact of quasispecies in virology, see Holland, 2006). The scientific presentations were opened by Weissmann and Domingo, who were the last and first authors of the 1978 paper, respectively. Their talks conveyed the scientific atmosphere of the 1970s—when molecular biology was carried out with few recombinant‐DNA techniques—to a young audience. Nucleic‐acid sequencing was in … [1]: /embed/graphic-1.gif

Genetic and catalytic efficiency structure of an HCV protease quasispecies†

The HCV nonstructural protein (NS)3/4A serine protease is not only involved in viral polyprotein processing but also efficiently blocks the retinoic-acid-inducible gen I and Toll-like receptor 3 signaling pathways and contributes to virus persistence by enabling HCV to escape the interferon antiviral response. Therefore, the NS3/4A protease has emerged as an ideal target for the control of the disease and the development of new anti-HCV agents. Here, we analyzed, at a high resolution (approximately 100 individual clones), the HCV NS3 protease gene quasispecies from three infected individuals. Nucleotide heterogeneity of 49%, 84%, and 91% were identified, respectively, which created a dense net that linked different parts of the viral population. Minority variants having mutations involved in the acquisition of resistance to current NS3/4A protease inhibitors (PIs) were also found. A vast diversity of different catalytic efficiencies could be distinguished. Importantly, 67% of the analyzed enzymes displayed a detectable protease activity. Moreover, 35% of the minority individual variants showed similar or better catalytic efficiency than the master (most abundant) enzyme. Nevertheless, and in contrast to minority variants, master enzymes always displayed a high catalytic efficiency when different viral polyprotein cleavage sites were tested. Finally, genetic and catalytic efficiency differences were observed when the 3 quasispecies were compared, suggesting that different selective forces were acting in different infected individuals. The rugged HCV protease quasispecies landscape should be able to react to environmental changes that may threaten its survival.

Fitness landscape of human immunodeficiency virus type 1 protease quasispecies.

Here we show, at a high resolution (1%), the human immunodeficiency virus type 1 (HIV-1) protease gene quasispecies landscape from three infected naïve individuals. A huge range of genetic configurations was found (67%, 71%, and 80% of the nucleotide clones from the three individuals, respectively, were different), and these configurations created a dense net that linked different parts of the viral population. Similarly, a vast diversity of different protease activities was also found. Importantly, 65% of the analyzed enzymes had detectable protease activity, and 11% of the minority individual variants showed similar or better fitness than the master (most abundant) enzyme, suggesting that the viral complexity in this genomic region does not exclusively depend on the enzyme's catalytic efficiency. Several high-fitness minority variants had only one substitution compared to the master sequence, supporting the possibility that the rugged HIV-1 protease quasispecies fitness landscape may be formed by a continuous network that can be traversed by single mutational steps without passing through defective or less-adapted proteins.

Evidence for Host-Driven Selection of the HIV Type 1vprGenein Vivoduring HIV Disease Progression in a Transfusion-Acquired Cohort

An epidemiologically linked HIV-1-infected cohort, in which a nonprogressor donor infected two recipients who progressed to AIDS, was examined. Sequence analysis, over time, of HIV-1 vpr gene quasispecies from uncultured peripheral blood cells revealed an insertion of arginine at position 90 altering a highly conserved C-terminal motif, believed to play a role in Vpr nuclear targeting. Full genome analysis from each patient showed no gene defects in other gene regions, implying that the mutational selection was unique to the vpr gene. A detailed analysis of the vpr quasispecies showed very little amino acid diversity in the nonprogressing donor, whereas, following viral transmission, the amino acid diversity increased dramatically over time in tandem with disease progression in the two recipients. Although the R insertion at position 90 was present in all three individuals, the variable degree of additional amino acid changes over time may have influenced HIV disease in the nonprogressor donor and the two progressing recipients. These data provide the first evidence in favor of vpr gene evolution over time, which was host-driven. The status of the nonprogressing donor was consistent with a highly protective B-57 HLA type, which was absent in the two progressing recipients, implying a role for host HLA type and other immunologic selective pressures in vpr gene selection in vivo.

Population genetic analysis of the protease locus of human immunodeficiency virus type 1 quasispecies undergoing drug selection, using a denaturing gradient-heteroduplex tracking assay.

Monitoring the evolution of human immunodeficiency virus type 1 (HIV-1) drug resistance requires measuring the frequency of closely related genetic variants making up the complex viral quasispecies found in vivo. In order to resolve both major and minor (>/=2%) protease gene variants differing by one or more nucleotide substitutions, we analyzed PCR products derived from plasma viral quasispecies by using a combination of denaturing gradient gel electrophoresis and DNA heteroduplex tracking assays. Correct population sampling of the high level of genetic diversity present within viral quasispecies could be documented by parallel analysis of duplicate, independently generated PCR products. The composition of genetically complex protease gene quasispecies remained constant over short periods of time in the absence of treatment and while plasma viremia fell >100-fold following the initiation of protease inhibitor ritonavir monotherapy. Within a month of initiating therapy, a strong reduction in the genetic diversity of plasma viral populations at the selected protease locus was associated with rising plasma viremia and the emergence of drug resistance. The high levels of protease genetic diversity seen before treatment reemerged only months later. In one patient, reduction in genetic diversity at the protease gene was observed concomitantly with an increase in diversity at the envelope gene (E. L. Delwart, P. Heng, A. Neumann, and M. Markowitz, J. Virol. 72:2416-2421, 1998), indicating that opposite population genetic changes can take place in different HIV-1 loci. The rapid emergence of drug-resistant HIV-1 was therefore associated with a strong, although only transient, reduction in genetic diversity at the selected locus. The denaturing gradient-heteroduplex tracking assay is a simple method for the separation and quantitation of very closely related, low-frequency, genetic variants within complex viral populations.

Active HIV-1 redistribution and replication in the brain with HIV encephalitis.

The central nervous system (CNS) is of particular importance in human immunodeficiency virus type 1 (HIV-1) infection. First, the CNS may be difficult to access for anti-retroviral treatment and may become a sanctuary for residual viruses. Second, HIV-1 infection may lead to AIDS dementia complex (ADC) culminating in HIV-1 encephalitis. In order to examine the pattern of drug resistance and the role of encephalitis in enhancing viral redistribution to the CNS, we compared pol gene quasispecies of the spleen and brain in two patients with and two patients without HIV-1 encephalitis, who had been treated with zidovudine (AZT). Although a variable degree of AZT resistance was noted in both the spleen and brain of all patients, phylogenetic analysis indicated that quasispecies developed rather independently in the systemic circulation (spleen) and CNS (brain) of patients without HIV-1 encephalitis, while similar pol gene sequences were obtained from the two compartments of patients with HIV-1 encephalitis. env gene V3 region of patients with HIV-1 encephalitis showed distinct quasispecies in the spleen and brain. Our results suggest that HIV-1 redistribution to CNS is more active in cases with encephalitis and that HIV-1 distributed late to CNS grow actively under certain selective pressure exerted on the V3 region of the env gene.

Coexisting members of HIV-1 p17 gene quasispecies represent proteins with distinct antigenicity and immunogenicity.

Despite the comparatively conserved nature of the HIV-1 gag gene, countless quasispecies of the p17 gene coexist in HIV-1-infected patients. It is not known if the minor genetic differences in quasispecies will affect immune recognition. To characterize the antigenicity and immunogenicity of three different members of HIV-1 p17 quasispecies. Three members of HIV-1 p17 gene quasispecies, one from patient A (clone 9; qsA9) and two from patient E (clones 5 and 8; qsE5 and qsE8), were expressed and purified from Escherichia coli. The antigenicity of the p17 proteins was analysed using sera from HIV-1-infected individuals, and the immunogenicity was evaluated using sera and lymphocytes from primed mice of three different haplotypes. The antigenicity of the qsE5 and qsE8 p17 recombinant proteins were distinct when tested for reactivity with human p17 antibodies. The qsE5 and qsE8 p17 were equally immunogenic in H-2d mice, but not in H-2b and H-2k mice. In H-2b mice the qsE8 protein induced higher levels of anti-p17 IgG2a, IgG2b and IgG3 than the qsE5 protein. Corroborating the IgG subclass pattern, H-2b-restricted qsE5-specific T cells produced higher in vitro levels of interferon-gamma, but not of interleukin (IL)-4, IL-5 and IL-6, than qsE8-specific T cells, suggesting a more pronounced T-helper (TH)1-like response. The p17 gene quasispecies coexisting in the same patient at the same time may represent antigenically and immunogenically distinct proteins despite sequence homologies of above 90%. Subsequently, subtle differences between two p17 protein quasispecies are enough to prime different TH1/TH2 subsets. These findings will have implications for therapeutic HIV-1 immunizations.

Pol gene quasispecies of human immunodeficiency virus: mutations associated with drug resistance in virus from patients undergoing no drug therapy.

The nucleotide sequences of two pol gene regions (codons 41 to 108 and 181 to 219 of reverse transcriptase) of 60 human immunodeficiency virus type 1 genomes obtained directly from primary lymphocytes from infected individuals are reported. In addition, the mutant spectra of several quasispecies have been sampled by repetitive sequencing of molecular clones representing the same pol genomic regions. Average mutation frequencies ranged from 1.6 x 10(-2) to 3.4 x 10(-2) substitutions per nucleotide for independent samples (relative to their consensus nucleotide sequence) and from 3.6 x 10(-3) to 1.1 x 10(-2) substitutions per nucleotide for individual quasispecies distributions. Several mutations leading to amino acid substitutions related to loss of sensitivity to reverse transcriptase inhibitors have been identified in samples from patients not subjected to antiretroviral therapy. Mutation frequencies in the codons previously identified as involved in resistance to reverse transcriptase inhibitors were very similar to the average mutation frequencies in the pol region analyzed. Thus, the finding of mutations related to drug resistance (even in the absence of positive selection by the corresponding drugs) is the expected consequence of the statistical distribution of mutations along the pol gene. The presence of such critical amino acid replacements in human immunodeficiency virus type 1 populations underscores the importance of viral quasispecies as reservoirs of phenotypic virus variants and has a number of implications for AIDS control.