儿童与成人慢性乙型肝炎患者乙型肝炎病毒Core基因区准种特征及正选择压力差异分析

儿童与成人慢性乙型肝炎患者的临床特征差异明显。乙型肝炎病毒(Hepatitis B virus, HBV)病毒准种特征与其致病特性紧密相连,HBV病毒Core 基因区富含免疫表位,该区域的准种特征直接反映病毒变异与病毒应对宿主免疫压力间的动态过程。文章通过扩增170名儿童慢性乙型肝炎患者及121名成人慢性乙型肝炎患者病毒Core基因区,按照病毒基因型以及病毒e抗原(Hepatitis B virus e antigen, HBeAg)状态进行分组,使用序列复杂度、多样性、非同义突变率(Non-synonymous substitution ratio,dN)、同义突变率(Synonymous substitution ratios , dS)等指标衡量不同组别之间的病毒准种特征;使用不同模型计算不同组别中受到正选择压力的位点,进一步结合HBV Core基因区免疫表位信息,进行正选择位点的定位分析。结果发现,儿童乙型肝炎病毒患者体内病毒Core基因区序列复杂性和多样性低于成人患者,且前者Core基因区正选择位点个数显著低于后者,这说明儿童慢性乙型肝炎患者体内病毒受到的选择压力低于成人患者。在儿童及成人慢性感染病人组中,HBeAg阳性病人体内病毒受到的选择压力低于HBeAg阴性病人。儿童及成人慢性感染患者体内病毒存在13个正选择位点,大多数正选择位点位于已知的抗原表位上。本研究从分子进化角度揭示了儿童与成人慢性乙型肝炎病例体内病毒Core基因区序列准种差异,为两类病人显著不同的临床表征提供了群体遗传学的解释。     

PAQ: Partition Analysis of Quasispecies

MOTIVATION: The complexities of genetic data may not be accurately described by any single analytical tool. Phylogenetic analysis is often used to study the genetic relationship among different sequences. Evolutionary models and assumptions are invoked to reconstruct trees that describe the phylogenetic relationship among sequences. Genetic databases are rapidly accumulating large amounts of sequences. Newly acquired sequences, which have not yet been characterized, may require preliminary genetic exploration in order to build models describing the evolutionary relationship among sequences. There are clustering techniques that rely less on models of evolution, and thus may provide nice exploratory tools for identifying genetic similarities. Some of the more commonly used clustering methods perform better when data can be grouped into mutually exclusive groups. Genetic data from viral quasispecies, which consist of closely related variants that differ by small changes, however, may best be partitioned by overlapping groups. RESULTS: We have developed an intuitive exploratory program, Partition Analysis of Quasispecies (PAQ), which utilizes a non-hierarchical technique to partition sequences that are genetically similar. PAQ was used to analyze a data set of human immunodeficiency virus type 1 (HIV-1) envelope sequences isolated from different regions of the brain and another data set consisting of the equine infectious anemia virus (EIAV) regulatory gene rev. Analysis of the HIV-1 data set by PAQ was consistent with phylogenetic analysis of the same data, and the EIAV rev variants were partitioned into two overlapping groups. PAQ provides an additional tool which can be used to glean information from genetic data and can be used in conjunction with other tools to study genetic similarities and genetic evolution of viral quasispecies.     

Evidence for the non-quasispecies evolution of RNA viruses [corrected

The quasispecies model of RNA virus evolution differs from those formulated in conventional population genetics in that neutral mutations do not lead to genetic drift of the population, and natural selection acts on the mutant distribution as a whole rather than on individual variants. By computer simulation, we show that this model could be inappropriate for many RNA viruses because the neutral sequence space may be too large to allow the formation of a quasispecies distribution. This view is supported by our analysis of gene sequences from vesicular stomatitis virus, which is considered a prototype RNA virus quasispecies. Our results are relevant to the evolution of RNA systems in general.     

Compartmentalization of hepatitis C virus quasispecies in blood mononuclear cells of patients with mixed cryoglobulinemic syndrome

The aim of this study was to investigate the quasispecies heterogeneity of hepatitis C virus (HCV) in the plasma, cryoprecipitate, and peripheral lymphocytes of chronically infected HCV patients with mixed cryoglobulinemia (MC). We studied 360 clones from 10 HCV-positive patients with MC and 8 age-, gender- and HCV genotype-matched subjects with chronic HCV infection but without MC. A partial nucleotide sequence encompassing the E1/E2 region, including hypervariable region 1 (HVR1), was amplified and cloned from plasma, cryoprecipitates, and peripheral blood mononuclear cells (PBMC), and the genetic diversity and complexity and synonymous and nonsynonymous substitution rates were determined. Heterogeneous selection pressure at codon sites was evaluated. Compartmentalization was estimated by phylogenetic and phenetic (Mantel's test) approaches. The patients with MC had 3.3 times lower nonsynonymous substitution rates (1.7 versus 5.7 substitutions/100 sites). Among the subjects with HCV genotype 1, the MC patients had significantly less complexity than the controls, whereas the diversity and complexity were similar in the genotype 2 patients and controls. Site-specific selection analysis confirmed the low frequency of MC patients showing positive selection. There was a significant correlation between positive selection and the infecting HCV genotype. The quasispecies were less heterogeneous in PBMC than in plasma. Significant compartmentalization of HCV quasispecies was observed in the PBMC of four of nine subjects (three with MC) and seven of nine cryoprecipitates. In one subject with MC, we detected a 5-amino-acid insertion at codons 385 to 389 of HVR1. Our results suggest reduced quasispecies heterogeneity in MC patients that is related to a low selection pressure which is probably due to an impaired immune response, the HCV genotype, and/or the duration of the infection. The frequent HCV quasispecies compartmentalization in patients' PBMC suggests a possible pathogenetic significance.     

Human immunodeficiency virus type 1 evolution in vivo tracked by DNA heteroduplex mobility assays.

High mutation rates and strong selective pressures imposed on human immunodeficiency viruses in vivo result in the formation of pools of genetic variants known as quasispecies. DNA heteroduplex mobility and tracking analyses were used to monitor the generation of HIV sequence diversity, to estimate quasispecies complexity, and to assess the turnover of genetic variants to approach an understanding of the relationship between viral quasispecies evolution in vivo and disease progression. Proviral DNA pools were nearly homogeneous soon after sexual transmission. The emergence and clearance of individual variants then occurred at different rates in different individuals. High quasispecies complexity was found in long-term-infected, asymptomatic individuals, while rapid CD4+ cell decline and AIDS were often, but not always, associated with lower quasispecies complexity. Proviral genetic variation was often low following in vitro culture, because of the outgrowth of one or a few variants that often became more abundant only later as proviruses in peripheral blood mononuclear cells. These studies provide insight into the dynamics of human immunodeficiency virus sequence changes in vivo and illustrate the utility of heteroduplex analysis for the study of phenomena associated with rapid genetic changes.     

Ecological connectivity shapes quasispecies structure of RNA viruses in an Antarctic lake

RNA viruses exist as complex mixtures of genotypes, known as quasispecies, where the evolution potential resides in the whole community of related genotypes. Quasispecies structure and dynamics have been studied in detail for virus infecting animals and plants but remain unexplored for those infecting micro‐organisms in environmental samples. We report the first metagenomic study of RNA viruses in an Antarctic lake (Lake Limnopolar, Livingston Island). Similar to low‐latitude aquatic environments, this lake harbours an RNA virome dominated by positive single‐strand RNA viruses from the order Picornavirales probably infecting micro‐organisms. Antarctic picorna‐like virus 1 (APLV1), one of the most abundant viruses in the lake, does not incorporate any mutation in the consensus sequence from 2006 to 2010 and shows stable quasispecies with low‐complexity indexes. By contrast, APLV2‐APLV3 are detected in the lake water exclusively in summer samples and are major constituents of surrounding cyanobacterial mats. Their quasispecies exhibit low complexity in cyanobacterial mat, but their run‐off‐mediated transfer to the lake results in a remarkable increase of complexity that may reflect the convergence of different viral quasispecies from the catchment area or replication in a more diverse host community. This is the first example of viral quasispecies from natural aquatic ecosystems and points to ecological connectivity as a modulating factor of quasispecies complexity.     

What is a quasispecies?

A quasispecies is a well-defined distribution of mutants that is generated by a mutation-selection process. Selection does not act on a single mutant but on the quasispecies as a whole. Experimental systems have been designed to study quasispecies evolution under laboratory conditions. More recently, virus populations have been called quasispecies to indicate their extensive genetic heterogeneity. The most prominent examples are probably the human immunodeficiency viruses HIV-1 and HIV-2. The quasispecies nature of HIV has formed the basis of a model that provides a mechanism for the pathogenesis of acquired immunodeficiency syndrome (AIDS) in humans. This article focuses on the nature of the quasispecies concept and its implications for evolutionary biology and virology.     

Tracking hepatitis C virus quasispecies major and minor variants in symptomatic and asymptomatic liver transplant recipients.

To evaluate the possibility that distinct viral quasispecies play a role in the pathogenesis of progressive hepatitis C virus (HCV) infection, we performed a detailed evaluation of HCV quasispecies before and after liver transplantation in five patients infected with HCV genotype 1, three of whom developed severe recurrent hepatitis C and two of whom developed asymptomatic posttransplant infections with high-titered viremia. HCV quasispecies were characterized by using a combination of nucleotide sequencing plus heteroduplex tracking assay of the second envelope gene hypervariable region (HVR). An average of 30 HVR clones were analyzed per specimen; an average of five specimens were analyzed per patient over a 6- to 24-month study period. The complexity of HCV quasispecies in pretransplant serum varied, ranging from one to nine genetically distinct variants for the five patients. However, in all five cases, relatively homogenous quasispecies variants emerged after liver transplantation. In the three patients who developed recurrent hepatitis, quasispecies major variants present in pretransplant serum were efficiently propagated immediately after liver transplantation and were propagated throughout the course of acute and chronic hepatitis. In contrast, in the two asymptomatic cases, we observed rapid depletion of pretransplant quasispecies major variants from posttransplant serum, followed by emergence of new quasispecies variants by posttransplant day 30. Genetic analysis suggested that in these cases, the new quasispecies variants were derived from minor variants present at relatively low clonal frequency (less than 5% of HVR clones) within the pretransplant quasispecies populations. These data demonstrate that quasispecies tracking patterns are associated with the rapidity and severity of HCV-associated liver disease after liver transplantation. Further characterization of HCV quasispecies in animal model systems is warranted.     

Impact of highly active antiretroviral therapy and immunologic status on hepatitis C virus quasispecies diversity in human immunodeficiency virus/hepatitis C virus-coinfected patients

This study analyzes the effect of highly active antiretroviral therapy (HAART), and thus immunologic status, on hepatitis C virus (HCV) load and quasispecies diversity in patients coinfected with the human immunodeficiency virus (HIV) and HCV. Three cohorts of coinfected patients were analyzed retrospectively over a period of 7 to 10 months: group A was antiretroviral drug na?ve at baseline and then on HAART for the remainder of the study, group B did not receive antiretroviral therapy at any point, and group C was on HAART for the entire study. HCV quasispecies diversity was analyzed by sequencing hypervariable region 1. In a longitudinal analysis, there was no significant change from baseline in any immunologic, virologic, or quasispecies parameter in any of the three groups. However, in comparison to groups A and B, group C had significantly higher CD4+- and CD8+-cell counts, a trend toward a higher HCV load, and significantly increased number of HCV clones, entropy, genetic distance, and ratio of nonsynonymous substitutions per nonsynonymous site to synonymous substitutions per synonymous site (Ka/Ks). In addition, CD4+-cell count was positively correlated with HCV load, genetic distance, and Ka. Interestingly, patients infected with HCV genotype 2 or 3 had a significantly higher CD4+-cell count, HCV load, genetic distance, and Ka/Ks than those infected with genotype 1. These results suggest that there is no immediate effect of HAART on HCV but that, with prolonged HAART, immune restoration results in an increase in HCV load and quasispecies diversity.     

Effect of Retreatment with Interferon Alone or Interferon plus Ribavirin on Hepatitis C Virus Quasispecies Diversification in Nonresponder Patients with Chronic Hepatitis C

Alpha interferon (IFN-alpha) treatment is effective on a long-term basis in only 15 to 25% of patients with chronic hepatitis C. The results of recent trials indicate that response rates can be significantly increased when IFN-alpha is given in combination with ribavirin. However, a large number of patients do not respond even to combination therapy. Nonresponsiveness to IFN is characterized by evolution of the hepatitis C virus (HCV) quasispecies. Little is known about the changes occurring within the HCV genomes when nonresponder patients are retreated with IFN or with IFN plus ribavirin. In the present study we have examined the genetic divergence of HCV quasispecies during unsuccessful retreatment with IFN or IFN plus ribavirin. Fifteen nonresponder patients with HCV-1 (4 patients with HCV-1a and 11 patients with HCV-1b) infection were studied while being retreated for 2 months (phase 1) with IFN-alpha (6 MU given three times a week), followed by IFN plus ribavirin or IFN alone for an additional 6 months (phase 2). HCV quasispecies diversification in the E2 hypervariable region-1 (HVR1) and in the putative NS5A IFN sensitivity determining region (ISDR) were analyzed for phase 1 and phase 2 by using the heteroduplex tracking assay and clonal frequency analysis techniques. A major finding of this study was the relatively rapid evolution of the HCV quasispecies observed in both treatment groups during the early phase 1 compared to the late phase 2 of treatment. The rate of quasispecies diversification in HVR1 was significantly higher during phase 1 versus phase 2 both in patients who received IFN plus ribavirin (P = 0.017) and in patients who received IFN alone (P = 0. 05). A trend toward higher rates of quasispecies evolution in the ISDR was also observed during phase 1 in both groups, although the results did not reach statistical significance. However, the NS5A quasispecies appeared to be rather homogeneous and stable in most nonresponder patients, suggesting the presence of a single well-fit major variant, resistant to antiviral treatment, in agreement with published data which have identified an IFN sensitivity determinant region within the NS5A. During the entire 8 months of retreatment, there was no difference in the rate of fixation of mutation between patients who received combination therapy and patients who were treated with IFN alone, suggesting that ribavirin had no major effects on the evolution of the HCV quasispecies after the initial 2 months of IFN therapy.     

The effect of transcutaneous application of carbon dioxide (CO₂) on skeletal muscle

Quasispecies is a remarkable characteristic of hepatitis C virus (HCV) and has profound roles in HCV biology and clinical practice. The understanding of HCV quasispecies behavior, in particular in acute HCV infection, is valuable for vaccine development and therapeutic interference. However, acute HCV infection is seldom encountered in clinic practice due to its silent onset. In the present study, we reported a unique case of de novo HCV infection associated with the transplantation of bone marrow from a HCV-positive donor. HCV quasispecies diversity was determined in both the donor and the recipient over a 4-year follow-up, accompanied with simultaneous measurement of HCV neutralizing antibody. Detailed genetic and phylogenetic analyses revealed a divergent quasispecies evolution, which was not related to dynamic changes of HCV neutralizing antibody. Instead, our data suggested an essential role of the fitness adaptation of founder viral population in driving such an evolutionary pattern.     

Origin of human immunodeficiency virus type 1 quasispecies emerging after antiretroviral treatment interruption in patients with therapeutic failure

The emergence of antiretroviral (ARV) drug-resistant human immunodeficiency virus type 1 (HIV-1) quasispecies is a major cause of treatment failure. These variants are usually replaced by drug-sensitive ones when the selective pressure of the drugs is removed, as the former have reduced fitness in a drug-free environment. This was the rationale for the design of structured ARV treatment interruption (STI) studies for the management of HIV-1 patients with treatment failure. We have studied the origin of drug-sensitive HIV-1 quasispecies emerging after STI in patients with treatment failure due to ARV drug resistance. Plasma and peripheral blood mononuclear cell samples were obtained the day of treatment interruption (day 0) and 30 and 60 days afterwards. HIV-1 pol and env were partially amplified, cloned, and sequenced. At day 60 drug-resistant variants were replaced by completely or partially sensitive quasispecies. Phylogenetic analyses of pol revealed that drug-sensitive variants emerging after STI were not related to their immediate temporal ancestors but formed a separate cluster, demonstrating that STI leads to the recrudescence and reemergence of a sequestrated viral population rather than leading to the back mutation of drug-resistant forms. No evidence for concomitant changes in viral tropism was seen, as deduced from env sequences. This study demonstrates the important role that the reemergence of quasispecies plays in HIV-1 population dynamics and points out the difficulties that may be found when recycling ARV therapies with patients with treatment failure.     

In silico identification of functional divergence between the multiple groEL gene paralogs in Chlamydiae

Background

RNA molecules, through their dual appearance as sequence and structure, represent a suitable model to study evolutionary properties of quasispecies. The essential ingredient in this model is the differentiation between genotype (molecular sequences which are affected by mutation) and phenotype (molecular structure, affected by selection). This framework allows a quantitative analysis of organizational properties of quasispecies as they adapt to different environments, such as their robustness, the effect of the degeneration of the sequence space, or the adaptation under different mutation rates and the error threshold associated.

Results

We describe and analyze the structural properties of molecular quasispecies adapting to different environments both during the transient time before adaptation takes place and in the asymptotic state, once optimization has occurred. We observe a minimum in the adaptation time at values of the mutation rate relatively far from the phenotypic error threshold. Through the definition of a consensus structure, it is shown that the quasispecies retains relevant structural information in a distributed fashion even above the error threshold. This structural robustness depends on the precise shape of the secondary structure used as target of selection. Experimental results available for natural RNA populations are in qualitative agreement with our observations.

Conclusion

Adaptation time of molecular quasispecies to a given environment is optimized at values of the mutation rate well below the phenotypic error threshold. The optimal value results from a trade-off between diversity generation and fixation of advantageous mutants. The critical value of the mutation rate is a function not only of the sequence length, but also of the specific properties of the environment, in this case the selection pressure and the shape of the secondary structure used as target phenotype. Certain functional motifs of RNA secondary structure that withstand high mutation rates (as the ubiquitous hairpin motif) might appear early in evolution and be actually frozen evolutionary accidents.     

Ultra-Deep Pyrosequencing (UDPS) Data Treatment to Study Amplicon HCV Minor Variants

We have investigated the reliability and reproducibility of HCV viral quasispecies quantification by ultra-deep pyrosequencing (UDPS) methods. Our study has been divided in two parts. First of all, by UDPS sequencing of clone mixes samples we have established the global noise level of UDPS and fine tuned a data treatment workflow previously optimized for HBV sequence analysis. Secondly, we have studied the reproducibility of the methodology by comparing 5 amplicons from two patient samples on three massive sequencing platforms (FLX+, FLX and Junior) after applying the error filters developed from the clonal/control study. After noise filtering the UDPS results, the three replicates showed the same 12 polymorphic sites above 0.7%, with a mean CV of 4.86%. Two polymorphic sites below 0.6% were identified by two replicates and one replicate respectively. A total of 25, 23 and 26 haplotypes were detected by GS-Junior, GS-FLX and GS-FLX+. The observed CVs for the normalized Shannon entropy (Sn), the mutation frequency (Mf), and the nucleotidic diversity (Pi) were 1.46%, 3.96% and 3.78%. The mean absolute difference in the two patients (5 amplicons each), in the GS-FLX and GS-FLX+, were 1.46%, 3.96% and 3.78% for Sn, Mf and Pi. No false polymorphic site was observed above 0.5%.Our results indicate that UDPS is an optimal alternative to molecular cloning for quantitative study of HCV viral quasispecies populations, both in complexity and composition. We propose an UDPS data treatment workflow for amplicons from the RNA viral quasispecies which, at a sequencing depth of at least 10,000 reads per strand, enables to obtain sequences and frequencies of consensus haplotypes above 0.5% abundance with no erroneous mutations, with high confidence, resistant mutants as minor variants at the level of 1%, with high confidence that variants are not missed, and highly confident measures of quasispecies complexity.     

The relationship between the error catastrophe, survival of the flattest, and natural selection

Background  

The quasispecies model is a general model of evolution that is generally applicable to replication up to high mutation rates. It predicts that at a sufficiently high mutation rate, quasispecies with higher mutational robustness can displace quasispecies with higher replicative capacity, a phenomenon called "survival of the flattest". In some fitness landscapes it also predicts the existence of a maximum mutation rate, called the error threshold, beyond which the quasispecies enters into error catastrophe, losing its genetic information. The aim of this paper is to study the relationship between survival of the flattest and the transition to error catastrophe, as well as the connection between these concepts and natural selection.     

Combinatorial analysis and algorithms for quasispecies reconstruction using next-generation sequencing

Background  

Next-generation sequencing (NGS) offers a unique opportunity for high-throughput genomics and has potential to replace Sanger sequencing in many fields, including de-novo sequencing, re-sequencing, meta-genomics, and characterisation of infectious pathogens, such as viral quasispecies. Although methodologies and software for whole genome assembly and genome variation analysis have been developed and refined for NGS data, reconstructing a viral quasispecies using NGS data remains a challenge. This application would be useful for analysing intra-host evolutionary pathways in relation to immune responses and antiretroviral therapy exposures. Here we introduce a set of formulae for the combinatorial analysis of a quasispecies, given a NGS re-sequencing experiment and an algorithm for quasispecies reconstruction. We require that sequenced fragments are aligned against a reference genome, and that the reference genome is partitioned into a set of sliding windows (amplicons). The reconstruction algorithm is based on combinations of multinomial distributions and is designed to minimise the reconstruction of false variants, called in-silico recombinants.     

Monitoring sequence space as a test for the target of selection in viruses

An essential feature of viral quasispecies, predicted from quasispecies theory, is that the target of selection is the mutant distribution as a whole. To test molecularly the mutant composition selected from a viral quasispecies we reconstructed a mutant distribution using 19 antigenic variants of foot-and-mouth disease virus (FMDV). Each variant was marked by a specific amino acid replacement at a major antigenic site of the virus that conferred resistance to a monoclonal antibody (mAb). The variants were introduced in the mutant spectrum of a biological FMDV clone, at a frequency commonly found in FMDV quasispecies. The reconstructed quasispecies (and a number of control populations) were allowed to replicate in the presence or absence of the mAb. The mutant distribution that became dominant as a result of antibody selection included at least ten of the 19 mutants initially used to reconstruct the quasispecies. No such biased mutant repertoire was found in control populations. The results show that a mutant distribution was selected, and are incompatible with selection of an individual genome, which then generated multiple mutants upon further replication. An ample representation of variants immediately following a selection event should contribute to subsequent adaptability of the virus.     

Memory in retroviral quasispecies: experimental evidence and theoretical model for human immunodeficiency virus

Viral quasispecies may possess a molecular memory of their past evolutionary history, imprinted on minority components of the mutant spectrum. Here we report experimental evidence and a theoretical model for memory in retroviral quasispecies in vivo. Apart from replicative memory associated with quasispecies dynamics, retroviruses may harbour a "cellular" or "anatomical" memory derived from their integrative cycle and the presence of viral reservoirs in body compartments. Three independent sets of data exemplify the two kinds of memory in human immunodeficiency virus type 1 (HIV-1). The data provide evidence of re-emergence of sequences that were hidden in cellular or anatomical compartments for extended periods of infection, and recovery of a quasispecies from pre-existing genomes. We develop a three-component model that incorporates the essential features of the quasispecies dynamics of retroviruses exposed to selective pressures. Significantly, a numerical study based on this model is in agreement with the experimental data, further supporting the existence of both replicative and reservoir memory in retroviral quasispecies.     

利用单基因组扩增法对马传染性贫血病毒疫苗株异质性的分析

前期研究发现,马传染性贫血病毒（Equine infectious anemia virus EIAV）中国弱毒疫苗株并非单一病毒,而是由多种准种（quasispecies）组成的种群。阐明该疫苗株的具体构成,对于确定优势疫苗株和分析其在体内的进化具有重要意义。本研究比较了传统RNA病毒测序法（即bulk PCR）和单基因组扩增法（Single-genome Amplification, SGA）在扩增EIAV疫苗株囊膜表面蛋白gp90基因 V3~V5区序列上的差异。结果发现,利用SGA法和bulk PCR法获得的序列在组内差异率分别为1.84%和1.88%。进一步序列比较发现,SGA法扩增的序列中除了含有与bulk PCR法中同源性较高的序列外,还存在bulk PCR法未检出的含强毒株LN40特异性位点,以及单个氨基酸缺失的序列。上述序列的存在为该疫苗株“多克隆构成”假说提供了佐证。此外,在对抽样偏差分析中发现,由于疫苗株中各种病毒准种在量上的差异,使得传统bulk PCR法不能有效的扩增组成比例较低的病毒准种,而导致测得的序列组成不能完全代表实际情况。SGA法通过对单基因组分的扩增和测序,可避免bulk PCR法的以上缺陷,在分析以准种形式存在的RNA病毒序列方面具有独特的优势。